#include <users_modules/multifield_plasticity/src/RigidBodies.hpp>

Inheritance diagram for STLRigidBody:

Collaboration diagram for STLRigidBody:

Public Member Functions
	STLRigidBody (VectorDouble c_coords, VectorDouble roller_disp, int id)

MoFEMErrorCode	getRollerDataForTag ()

Tensor1< double, 3 >	getNormal (Tensor1< TPack3, 3 > &t_coords, Tensor1< TPack1, 3 > &t_disp)

Tensor1< double, 3 >	getNormal (Tensor1< TPack3, 3 > &t_coords, Tensor1< double, 3 > &t_disp)

Tensor1< double, 3 >	getNormal (Tensor1< double, 3 > &t_coords, Tensor1< double, 3 > &t_disp)

Tensor2< double, 3, 3 >	getDiffNormal (Tensor1< TPack3, 3 > &t_coords, Tensor1< TPack1, 3 > &t_disp, Tensor1< TPack1, 3 > &t_normal)

double	getGap (Tensor1< TPack3, 3 > &t_coords)

Tensor1< double, 3 >	getdGap (Tensor1< TPack3, 3 > &t_coords, Tensor1< TPack1, 3 > &t_normal)

template<typename T1 , typename T2 >
Tensor1< double, 3 > &	getNormalImpl (Tensor1< T1, 3 > &t_coords, Tensor1< T2, 3 > &t_disp)

template<typename T1 >
double	getGapImpl (Tensor1< T1, 3 > &t_coords)

template<typename T1 , typename T2 >
Tensor1< double, 3 >	getdGapImpl (Tensor1< T1, 3 > &t_coords, Tensor1< T2, 3 > &t_normal)

template<typename T1 , typename T2 , typename T3 >
Tensor2< double, 3, 3 >	getDiffNormalImpl (Tensor1< T1, 3 > &t_coords, Tensor1< T2, 3 > &t_disp, Tensor1< T3, 3 > &t_normal)

MoFEMErrorCode	getBodyOptions ()

Public Member Functions inherited from RigidBodyData
	RigidBodyData (VectorDouble c_coords, VectorDouble roller_disp, int id)

	RigidBodyData ()=delete

virtual	~RigidBodyData ()

virtual Tensor1< double, 3 >	getNormal (Tensor1< TPack3, 3 > &t_coords, Tensor1< TPack1, 3 > &t_disp)=0

virtual Tensor1< double, 3 >	getNormal (Tensor1< TPack3, 3 > &t_coords, Tensor1< double, 3 > &t_disp)=0

virtual Tensor1< double, 3 >	getNormal (Tensor1< double, 3 > &t_coords, Tensor1< double, 3 > &t_disp)=0

virtual Tensor2< double, 3, 3 >	getDiffNormal (Tensor1< TPack3, 3 > &t_coords, Tensor1< TPack1, 3 > &t_disp, Tensor1< TPack1, 3 > &t_normal)=0

virtual double	getGap (Tensor1< TPack3, 3 > &t_coords)=0

virtual Tensor1< double, 3 >	getdGap (Tensor1< TPack3, 3 > &t_coords, Tensor1< TPack1, 3 > &t_normal)=0

virtual MoFEMErrorCode	getBodyOptions ()=0

MoFEMErrorCode	computeRotationMatrix ()

Tensor1< double, 3 >	getBodyOffset ()

template<typename T1 , typename T2 >
Tensor1< double, 3 > &	getPointCoords (Tensor1< T1, 3 > &t_coords, Tensor1< T2, 3 > &t_disp)

MoFEMErrorCode	saveBasicDataOnTag (moab::Interface &moab_debug, EntityHandle &vertex)

virtual MoFEMErrorCode	getRollerDataForTag ()=0

Public Attributes
moab::Core	mbInstance

moab::Interface *	bodyMesh

shared_ptr< OrientedBoxTreeTool >	orientedBox

EntityHandle	treeMeshset

Tag	tagNormal

EntityHandle	triangle

Public Attributes inherited from RigidBodyData
const int	iD

Index< 'i', 3 >	i

Index< 'j', 3 >	j

Index< 'k', 3 >	k

VectorDouble	originCoords

VectorDouble	rollerDisp

VectorDouble	BodyDispScaled

VectorDouble	BodyDirectionScaled

double	gAp

Tensor1< double, 3 >	tNormal

Tensor1< double, 3 >	dGap

Tensor1< double, 3 >	pointCoords

Tensor1< double, 3 >	closestPoint

Tensor1< double, 3 >	defaultOrientation

Tensor1< double, 3 >	oRientation

Tensor2< double, 3, 3 >	rotationMat

Tensor2< double, 3, 3 >	diffNormal

boost::shared_ptr< TimeAccelerogram >	methodOpForRollerPosition

boost::shared_ptr< TimeAccelerogram >	methodOpForRollerDirection

int	bodyType

array< double, 9 >	dataForTags

Additional Inherited Members
Public Types inherited from RigidBodyData
using	TPack3 = PackPtr< double *, 3 >

using	TPack1 = PackPtr< double *, 1 >

Detailed Description

Definition at line 1243 of file RigidBodies.hpp.

Constructor & Destructor Documentation

◆ STLRigidBody()

STLRigidBody::STLRigidBody	(	VectorDouble	c_coords,
		VectorDouble	roller_disp,
		int	id
	)

inline

Definition at line 1251 of file RigidBodies.hpp.

      : RigidBodyData(c_coords, roller_disp, id) {
    bodyMesh = &mbInstance;
  }

Member Function Documentation

◆ getBodyOptions()

MoFEMErrorCode STLRigidBody::getBodyOptions ( )

inlinevirtual

Implements RigidBodyData.

Definition at line 1334 of file RigidBodies.hpp.

                                  {
    MoFEMFunctionBegin;
    CHKERR PetscOptionsBegin(PETSC_COMM_WORLD, "", "", "");
    PetscBool rflg;
    string param = "-body_file_name" + to_string(iD);
    char mesh_file_name[255];
    CHKERR PetscOptionsString(param.c_str(), "file name", "", "mesh.vtk",
                              mesh_file_name, 255, &rflg);
    if (string(mesh_file_name).substr(string(mesh_file_name).length() - 3) !=
        "vtk")
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "body surface should be vtk surface");
    bodyMesh->load_file(mesh_file_name, 0, "");
    Range tris;
    CHKERR bodyMesh->get_entities_by_dimension(0, 2, tris, false);
    orientedBox =
        make_shared<OrientedBoxTreeTool>(bodyMesh, "ROOTSETSURF", true);
    CHKERR orientedBox->build(tris, treeMeshset);
    char normal_tag_name[255] = "Normals";
    if (bodyMesh->tag_get_handle(normal_tag_name, tagNormal) != MB_SUCCESS) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Normal tag name %s cannot be found. Please "
              "provide the correct name. or run ./calculate_normals");
    }
 
    ierr = PetscOptionsEnd();
    CHKERRG(ierr);
    MoFEMFunctionReturn(0);
  }

◆ getdGap()

Tensor1< double, 3 > STLRigidBody::getdGap	(	Tensor1< TPack3, 3 > &	t_coords,
		Tensor1< TPack1, 3 > &	t_normal
	)

inlinevirtual

Implements RigidBodyData.

Definition at line 1285 of file RigidBodies.hpp.

                                                           {
    return getdGapImpl(t_coords, t_normal);
  }

◆ getdGapImpl()

template<typename T1 , typename T2 >

Tensor1< double, 3 > STLRigidBody::getdGapImpl	(	Tensor1< T1, 3 > &	t_coords,
		Tensor1< T2, 3 > &	t_normal
	)

inline

Definition at line 1307 of file RigidBodies.hpp.

                                                                  {
    return dGap;
  };

◆ getDiffNormal()

Tensor2< double, 3, 3 > STLRigidBody::getDiffNormal	(	Tensor1< TPack3, 3 > &	t_coords,
		Tensor1< TPack1, 3 > &	t_disp,
		Tensor1< TPack1, 3 > &	t_normal
	)

inlinevirtual

Implements RigidBodyData.

Definition at line 1275 of file RigidBodies.hpp.

                                                                    {
    return getDiffNormalImpl(t_coords, t_disp, t_normal);
  }

◆ getDiffNormalImpl()

template<typename T1 , typename T2 , typename T3 >

Tensor2< double, 3, 3 > STLRigidBody::getDiffNormalImpl	(	Tensor1< T1, 3 > &	t_coords,
		Tensor1< T2, 3 > &	t_disp,
		Tensor1< T3, 3 > &	t_normal
	)

inline

Definition at line 1313 of file RigidBodies.hpp.

                                                                           {
    Tensor2<double, 3, 3> diff_normal;
    Tensor1<double, 3> norm_diff;
    const double clean_gap = gAp;
    const double eps = 1e-8;
    for (int ii = 0; ii != 3; ++ii) {
      Tensor1<double, 3> pert_disp{t_disp(0), t_disp(1), t_disp(2)};
      pert_disp(ii) += eps;
      auto pert_normal = getNormal(t_coords, pert_disp);
      norm_diff(i) = (pert_normal(i) - t_normal(i)) / eps;
      dGap(ii) = (gAp - clean_gap) / eps;
 
      for (int jj = 0; jj != 3; ++jj) {
        diff_normal(jj, ii) = norm_diff(jj);
      }
    }
    return diff_normal;
  };

◆ getGap()

double STLRigidBody::getGap ( Tensor1< TPack3, 3 > & t_coords )

inlinevirtual

Implements RigidBodyData.

Definition at line 1281 of file RigidBodies.hpp.

                                                     {
    return getGapImpl(t_coords);
  }

◆ getGapImpl()

template<typename T1 >

double STLRigidBody::getGapImpl ( Tensor1< T1, 3 > & t_coords )

inline

Definition at line 1302 of file RigidBodies.hpp.

                                                                            {
    return gAp;
  };

◆ getNormal() [1/3]

Tensor1< double, 3 > STLRigidBody::getNormal	(	Tensor1< double, 3 > &	t_coords,
		Tensor1< double, 3 > &	t_disp
	)

inlinevirtual

Implements RigidBodyData.

Definition at line 1270 of file RigidBodies.hpp.

                                                           {
    return getNormalImpl(t_coords, t_disp);
  }

◆ getNormal() [2/3]

Tensor1< double, 3 > STLRigidBody::getNormal	(	Tensor1< TPack3, 3 > &	t_coords,
		Tensor1< double, 3 > &	t_disp
	)

inlinevirtual

Implements RigidBodyData.

Definition at line 1266 of file RigidBodies.hpp.

                                                           {
    return getNormalImpl(t_coords, t_disp);
  }

◆ getNormal() [3/3]

Tensor1< double, 3 > STLRigidBody::getNormal	(	Tensor1< TPack3, 3 > &	t_coords,
		Tensor1< TPack1, 3 > &	t_disp
	)

inlinevirtual

Implements RigidBodyData.

Definition at line 1262 of file RigidBodies.hpp.

                                                           {
    return getNormalImpl(t_coords, t_disp);
  }

◆ getNormalImpl()

template<typename T1 , typename T2 >

Tensor1< double, 3 > & STLRigidBody::getNormalImpl	(	Tensor1< T1, 3 > &	t_coords,
		Tensor1< T2, 3 > &	t_disp
	)

inline

Definition at line 1291 of file RigidBodies.hpp.

                                                                   {
    auto t_d = getPointCoords(t_coords, t_disp);
 
    CHKERR orientedBox->closest_to_location(&t_d(0), treeMeshset,
                                            &closestPoint(0), triangle);
    CHKERR bodyMesh->tag_get_data(tagNormal, &triangle, 1, &tNormal(0));
    gAp = tNormal(i) * (pointCoords(i) - closestPoint(i));
    return tNormal;
  };

◆ getRollerDataForTag()

MoFEMErrorCode STLRigidBody::getRollerDataForTag ( )

inlinevirtual

Implements RigidBodyData.

Definition at line 1256 of file RigidBodies.hpp.

                                       {
    MoFEMFunctionBeginHot;
    bodyType = STL;
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

◆ bodyMesh

moab::Interface* STLRigidBody::bodyMesh

Definition at line 1245 of file RigidBodies.hpp.

◆ mbInstance

moab::Core STLRigidBody::mbInstance

Definition at line 1244 of file RigidBodies.hpp.

◆ orientedBox

shared_ptr<OrientedBoxTreeTool> STLRigidBody::orientedBox

Definition at line 1246 of file RigidBodies.hpp.

◆ tagNormal

Tag STLRigidBody::tagNormal

Definition at line 1248 of file RigidBodies.hpp.

◆ treeMeshset

EntityHandle STLRigidBody::treeMeshset

Definition at line 1247 of file RigidBodies.hpp.

◆ triangle

EntityHandle STLRigidBody::triangle

Definition at line 1249 of file RigidBodies.hpp.

The documentation for this struct was generated from the following file:

users_modules/multifield_plasticity/src/RigidBodies.hpp

Public Member Functions

Public Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ STLRigidBody()

Member Function Documentation

◆ getBodyOptions()

◆ getdGap()

◆ getdGapImpl()

◆ getDiffNormal()

◆ getDiffNormalImpl()

◆ getGap()

◆ getGapImpl()

◆ getNormal() [1/3]

◆ getNormal() [2/3]

◆ getNormal() [3/3]

◆ getNormalImpl()

◆ getRollerDataForTag()

Member Data Documentation

◆ bodyMesh

◆ mbInstance

◆ orientedBox

◆ tagNormal

◆ treeMeshset

◆ triangle