#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 ()

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)

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

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

string	positionDataParamName

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 1208 of file RigidBodies.hpp.

Constructor & Destructor Documentation

◆ STLRigidBody()

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

Definition at line 1216 of file RigidBodies.hpp.

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

Member Function Documentation

◆ getBodyOptions()

MoFEMErrorCode STLRigidBody::getBodyOptions ( )

virtual

Implements RigidBodyData.

Definition at line 1299 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
	)

virtual

Implements RigidBodyData.

Definition at line 1250 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
	)

Definition at line 1272 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
	)

virtual

Implements RigidBodyData.

Definition at line 1240 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
	)

Definition at line 1278 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 )

virtual

Implements RigidBodyData.

Definition at line 1246 of file RigidBodies.hpp.

                                                      {
     return getGapImpl(t_coords);
   }

◆ getGapImpl()

template<typename T1 >

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

Definition at line 1267 of file RigidBodies.hpp.

                                                                             {
     return gAp;
   };

◆ getNormal() [1/3]

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

virtual

Implements RigidBodyData.

Definition at line 1235 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
	)

virtual

Implements RigidBodyData.

Definition at line 1231 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
	)

virtual

Implements RigidBodyData.

Definition at line 1227 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
	)

Definition at line 1256 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 ( )

virtual

Implements RigidBodyData.

Definition at line 1221 of file RigidBodies.hpp.

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

Member Data Documentation

◆ bodyMesh

moab::Interface* STLRigidBody::bodyMesh

Definition at line 1210 of file RigidBodies.hpp.

◆ mbInstance

moab::Core STLRigidBody::mbInstance

Definition at line 1209 of file RigidBodies.hpp.

◆ orientedBox

shared_ptr<OrientedBoxTreeTool> STLRigidBody::orientedBox

Definition at line 1211 of file RigidBodies.hpp.

◆ tagNormal

Tag STLRigidBody::tagNormal

Definition at line 1213 of file RigidBodies.hpp.

◆ treeMeshset

EntityHandle STLRigidBody::treeMeshset

Definition at line 1212 of file RigidBodies.hpp.

◆ triangle

EntityHandle STLRigidBody::triangle

Definition at line 1214 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