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

Inheritance diagram for RigidBodyData:

Collaboration diagram for RigidBodyData:

Public Types
using	TPack3 = PackPtr< double *, 3 >

using	TPack1 = PackPtr< double *, 1 >

Public Member Functions
	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
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

Detailed Description

Definition at line 29 of file RigidBodies.hpp.

Member Typedef Documentation

◆ TPack1

using RigidBodyData::TPack1 = PackPtr<double *, 1>

Definition at line 72 of file RigidBodies.hpp.

◆ TPack3

using RigidBodyData::TPack3 = PackPtr<double *, 3>

Definition at line 71 of file RigidBodies.hpp.

Constructor & Destructor Documentation

◆ RigidBodyData() [1/2]

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

inline

Definition at line 57 of file RigidBodies.hpp.

      : originCoords(c_coords), rollerDisp(roller_disp), iD(id),
        defaultOrientation(0., 0., 1.), oRientation(0., 0., 1.) {
    BodyDispScaled = rollerDisp;
    BodyDispScaled.clear();
    BodyDirectionScaled = VectorDouble({0, 0, 0});
    rotationMat(i, j) = kronecker_delta(i, j);
    std::fill(dataForTags.begin(), dataForTags.end(), 0);
    bodyType = PLANE;
  }

◆ RigidBodyData() [2/2]

RigidBodyData::RigidBodyData ( )

delete

◆ ~RigidBodyData()

virtual RigidBodyData::~RigidBodyData ( )

inlinevirtual

Definition at line 69 of file RigidBodies.hpp.

69{}

Member Function Documentation

◆ computeRotationMatrix()

MoFEMErrorCode RigidBodyData::computeRotationMatrix ( )

inline

Definition at line 91 of file RigidBodies.hpp.

                                                {
    // https://math.stackexchange.com/questions/180418
    MoFEMFunctionBeginHot;
 
    auto compute_rot = [&]() {
      MoFEMFunctionBeginHot;
      oRientation.normalize();
      Tensor1<double, 3> v;
      v(i) = levi_civita(i, j, k) * defaultOrientation(j) * oRientation(k);
      const double s = v.l2();
      if (abs(s) < 1e-16)
        MoFEMFunctionReturnHot(0);
      const double c = defaultOrientation(i) * oRientation(i);
      const double coef = 1. / (1. + c);
      Tensor2<double, 3, 3> t_v;
      t_v(i, j) = levi_civita(i, j, k) * v(k);
      rotationMat(i, j) = kronecker_delta(i, j) + t_v(i, j) +
                          (coef * t_v(i, k)) * t_v(k, j);
      // Tensor1<double, 3> test;
      // test(i) = rotationMat(j, i) * defaultOrientation(j);
      MoFEMFunctionReturnHot(0);
    };
 
    // if we have not changed the direction, compute rotation only once
 
    if (BodyDirectionScaled.empty())
      MoFEMFunctionReturnHot(0);
 
    if (methodOpForRollerDirection) {
      for (int i = 0; i != 3; ++i)
        oRientation(i) = BodyDirectionScaled(i);
      CHKERR compute_rot();
    } else {
      BodyDirectionScaled.resize(0);
      CHKERR compute_rot();
    }
 
    MoFEMFunctionReturnHot(0);
  };

◆ getBodyOffset()

Tensor1< double, 3 > RigidBodyData::getBodyOffset ( )

inline

Definition at line 131 of file RigidBodies.hpp.

                                            {
    VectorDouble pos = originCoords + BodyDispScaled;
    return Tensor1<double, 3>(pos(0), pos(1), pos(2));
  };

◆ getBodyOptions()

virtual MoFEMErrorCode RigidBodyData::getBodyOptions ( )

pure virtual

Implemented in PlaneRigidBody, SphereRigidBody, CylinderRigidBody, TorusRigidBody, RollerRigidBody, ConeRigidBody, and STLRigidBody.

◆ getdGap()

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

pure virtual

Implemented in PlaneRigidBody, SphereRigidBody, CylinderRigidBody, TorusRigidBody, RollerRigidBody, ConeRigidBody, and STLRigidBody.

◆ getDiffNormal()

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

pure virtual

Implemented in PlaneRigidBody, SphereRigidBody, CylinderRigidBody, TorusRigidBody, RollerRigidBody, ConeRigidBody, and STLRigidBody.

◆ getGap()

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

pure virtual

Implemented in PlaneRigidBody, SphereRigidBody, CylinderRigidBody, TorusRigidBody, RollerRigidBody, ConeRigidBody, and STLRigidBody.

◆ getNormal() [1/3]

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

pure virtual

Implemented in PlaneRigidBody, SphereRigidBody, CylinderRigidBody, TorusRigidBody, RollerRigidBody, ConeRigidBody, and STLRigidBody.

◆ getNormal() [2/3]

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

pure virtual

Implemented in PlaneRigidBody, SphereRigidBody, CylinderRigidBody, TorusRigidBody, RollerRigidBody, ConeRigidBody, and STLRigidBody.

◆ getNormal() [3/3]

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

pure virtual

Implemented in PlaneRigidBody, SphereRigidBody, CylinderRigidBody, TorusRigidBody, RollerRigidBody, ConeRigidBody, and STLRigidBody.

◆ getPointCoords()

template<typename T1 , typename T2 >

Tensor1< double, 3 > & RigidBodyData::getPointCoords	(	Tensor1< T1, 3 > &	t_coords,
		Tensor1< T2, 3 > &	t_disp
	)

inline

Definition at line 137 of file RigidBodies.hpp.

                                                                    {
    auto t_off = getBodyOffset();
    pointCoords(i) = t_coords(i) - t_off(i) + t_disp(i);
    return pointCoords;
  };

◆ getRollerDataForTag()

virtual MoFEMErrorCode RigidBodyData::getRollerDataForTag ( )

pure virtual

Implemented in PlaneRigidBody, SphereRigidBody, CylinderRigidBody, TorusRigidBody, RollerRigidBody, ConeRigidBody, and STLRigidBody.

◆ saveBasicDataOnTag()

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

inline

Definition at line 144 of file RigidBodies.hpp.

                                                          {
    MoFEMFunctionBeginHot;
 
    Tag tag0;
    int def_int = 0;
    std::array<double, 9> def_real;
    std::fill(def_real.begin(), def_real.end(), 0);
 
    CHKERR getRollerDataForTag();
 
    CHKERR moab_debug.tag_get_handle("_ROLLER_ID", 1, MB_TYPE_INTEGER, tag0,
                                     MB_TAG_CREAT | MB_TAG_SPARSE, &def_int);
    CHKERR moab_debug.tag_set_data(tag0, &vertex, 1, &this->iD);
 
    CHKERR moab_debug.tag_get_handle("_BODY_TYPE", 1, MB_TYPE_INTEGER, tag0,
                                     MB_TAG_CREAT | MB_TAG_SPARSE, &def_int);
    CHKERR moab_debug.tag_set_data(tag0, &vertex, 1, &this->bodyType);
 
    CHKERR computeRotationMatrix();
    CHKERR moab_debug.tag_get_handle("_ORIENTATION", 3, MB_TYPE_DOUBLE, tag0,
                                     MB_TAG_CREAT | MB_TAG_SPARSE, &def_real);
    CHKERR moab_debug.tag_set_data(tag0, &vertex, 1, &this->oRientation(0));
 
    CHKERR moab_debug.tag_get_handle("_ROLLER_DATA", 9, MB_TYPE_DOUBLE, tag0,
                                     MB_TAG_CREAT | MB_TAG_SPARSE, &def_real);
    CHKERR moab_debug.tag_set_data(tag0, &vertex, 1, this->dataForTags.data());
 
 
    // tags
    // int roller_id
    // int body_type
    // vec3 orientation
    // tensor9 roller_data
    //  0 radius1 = 5
    //  1 inner_radius = 1
    //  2 height = 1
    //  3 angle = 45
    //  4 angle_a = 45
    //  5 angle_b = 45
    //
 
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

◆ BodyDirectionScaled

VectorDouble RigidBodyData::BodyDirectionScaled

Definition at line 37 of file RigidBodies.hpp.

◆ BodyDispScaled

VectorDouble RigidBodyData::BodyDispScaled

Definition at line 36 of file RigidBodies.hpp.

◆ bodyType

int RigidBodyData::bodyType

Definition at line 54 of file RigidBodies.hpp.

◆ closestPoint

Tensor1<double, 3> RigidBodyData::closestPoint

Definition at line 42 of file RigidBodies.hpp.

◆ dataForTags

array<double, 9> RigidBodyData::dataForTags

Definition at line 55 of file RigidBodies.hpp.

◆ defaultOrientation

Tensor1<double, 3> RigidBodyData::defaultOrientation

Definition at line 43 of file RigidBodies.hpp.

◆ dGap

Tensor1<double, 3> RigidBodyData::dGap

Definition at line 40 of file RigidBodies.hpp.

◆ diffNormal

Tensor2<double, 3, 3> RigidBodyData::diffNormal

Definition at line 46 of file RigidBodies.hpp.

◆ gAp

double RigidBodyData::gAp

Definition at line 38 of file RigidBodies.hpp.

◆ i

Index<'i', 3> RigidBodyData::i

Definition at line 31 of file RigidBodies.hpp.

◆ iD

const int RigidBodyData::iD

Definition at line 30 of file RigidBodies.hpp.

◆ j

Index<'j', 3> RigidBodyData::j

Definition at line 32 of file RigidBodies.hpp.

◆ k

Index<'k', 3> RigidBodyData::k

Definition at line 33 of file RigidBodies.hpp.

◆ methodOpForRollerDirection

boost::shared_ptr<TimeAccelerogram> RigidBodyData::methodOpForRollerDirection

Definition at line 52 of file RigidBodies.hpp.

◆ methodOpForRollerPosition

boost::shared_ptr<TimeAccelerogram> RigidBodyData::methodOpForRollerPosition

Definition at line 51 of file RigidBodies.hpp.

◆ oRientation

Tensor1<double, 3> RigidBodyData::oRientation

Definition at line 44 of file RigidBodies.hpp.

◆ originCoords

VectorDouble RigidBodyData::originCoords

Definition at line 34 of file RigidBodies.hpp.

◆ pointCoords

Tensor1<double, 3> RigidBodyData::pointCoords

Definition at line 41 of file RigidBodies.hpp.

◆ rollerDisp

VectorDouble RigidBodyData::rollerDisp

Definition at line 35 of file RigidBodies.hpp.

◆ rotationMat

Tensor2<double, 3, 3> RigidBodyData::rotationMat

Definition at line 45 of file RigidBodies.hpp.

◆ tNormal

Tensor1<double, 3> RigidBodyData::tNormal

Definition at line 39 of file RigidBodies.hpp.

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

users_modules/multifield_plasticity/src/RigidBodies.hpp

Public Types

Public Member Functions

Public Attributes

Detailed Description

Member Typedef Documentation

◆ TPack1

◆ TPack3

Constructor & Destructor Documentation

◆ RigidBodyData() [1/2]

◆ RigidBodyData() [2/2]

◆ ~RigidBodyData()

Member Function Documentation

◆ computeRotationMatrix()

◆ getBodyOffset()

◆ getBodyOptions()

◆ getdGap()

◆ getDiffNormal()

◆ getGap()

◆ getNormal() [1/3]

◆ getNormal() [2/3]

◆ getNormal() [3/3]

◆ getPointCoords()

◆ getRollerDataForTag()

◆ saveBasicDataOnTag()

Member Data Documentation

◆ BodyDirectionScaled

◆ BodyDispScaled

◆ bodyType

◆ closestPoint

◆ dataForTags

◆ defaultOrientation

◆ dGap

◆ diffNormal

◆ gAp

◆ i

◆ iD

◆ j

◆ k

◆ methodOpForRollerDirection

◆ methodOpForRollerPosition

◆ oRientation

◆ originCoords

◆ pointCoords

◆ rollerDisp

◆ rotationMat

◆ tNormal