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

Inheritance diagram for TorusRigidBody:

Collaboration diagram for TorusRigidBody:

Public Member Functions
	TorusRigidBody (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 , typename T3 >
Tensor2< double, 3, 3 >	getDiffNormalImpl (Tensor1< T1, 3 > &t_coords, Tensor1< T2, 3 > &t_disp, Tensor1< T3, 3 > &t_normal)

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

virtual 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
double	rAdius

double	innerRadius

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

Constructor & Destructor Documentation

◆ TorusRigidBody()

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

inline

Definition at line 598 of file RigidBodies.hpp.

      : RigidBodyData(c_coords, roller_disp, id), rAdius(1), innerRadius(0.5) {
    // , cP(R, r, 1)
    // compute rotation matrix here
  }

Member Function Documentation

◆ getBodyOptions()

virtual MoFEMErrorCode TorusRigidBody::getBodyOptions ( )

inlinevirtual

Implements RigidBodyData.

Definition at line 715 of file RigidBodies.hpp.

                                          {
    MoFEMFunctionBegin;
    PetscBool rflg;
    int nb_dirs = 3;
    CHKERR PetscOptionsBegin(PETSC_COMM_WORLD, "", "", "");
    string param = "-inner_radius" + to_string(iD);
    CHKERR PetscOptionsScalar(param.c_str(), "set torus small radius", "", innerRadius,
                              &innerRadius, PETSC_NULL);
    param = "-radius" + to_string(iD);
    CHKERR PetscOptionsScalar(param.c_str(), "set torus large Radius", "", rAdius,
                              &rAdius, PETSC_NULL);
    param = "-direction" + to_string(iD);
    CHKERR PetscOptionsGetRealArray(NULL, NULL, param.c_str(), &oRientation(0),
                                    &nb_dirs, &rflg);
    ierr = PetscOptionsEnd();
    CHKERRG(ierr);
 
    if (rflg)
      CHKERR computeRotationMatrix();
 
    // cP = make_shared<ArbitrarySurfaceData>(R, r);
 
    MoFEMFunctionReturn(0);
  }

◆ getdGap()

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

inlinevirtual

Implements RigidBodyData.

Definition at line 634 of file RigidBodies.hpp.

                                                           {
    return getdGapImpl(t_coords, t_normal);
  }

◆ getdGapImpl()

template<typename T1 , typename T2 >

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

inline

Definition at line 710 of file RigidBodies.hpp.

                                                                  {
    return dGap;
  };

◆ getDiffNormal()

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

inlinevirtual

Implements RigidBodyData.

Definition at line 625 of file RigidBodies.hpp.

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

◆ getDiffNormalImpl()

template<typename T1 , typename T2 , typename T3 >

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

inline

Definition at line 688 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 TorusRigidBody::getGap ( Tensor1< TPack3, 3 > & t_coords )

inlinevirtual

Implements RigidBodyData.

Definition at line 631 of file RigidBodies.hpp.

                                                     {
    return getGapImpl(t_coords);
  }

◆ getGapImpl()

template<typename T1 >

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

inline

Definition at line 683 of file RigidBodies.hpp.

                                                                            {
    return gAp;
  };

◆ getNormal() [1/3]

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

inlinevirtual

Implements RigidBodyData.

Definition at line 620 of file RigidBodies.hpp.

                                                           {
    return getNormalImpl(t_coords, t_disp);
  }

◆ getNormal() [2/3]

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

inlinevirtual

Implements RigidBodyData.

Definition at line 616 of file RigidBodies.hpp.

                                                           {
    return getNormalImpl(t_coords, t_disp);
  }

◆ getNormal() [3/3]

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

inlinevirtual

Implements RigidBodyData.

Definition at line 612 of file RigidBodies.hpp.

                                                           {
    return getNormalImpl(t_coords, t_disp);
  }

◆ getNormalImpl()

template<typename T1 , typename T2 >

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

inline

Definition at line 640 of file RigidBodies.hpp.

                                                                   {
 
    auto t_d = getPointCoords(t_coords, t_disp);
    CHKERR computeRotationMatrix();
    pointCoords(i) = rotationMat(i, j) * t_d(j);
    t_d(i) = pointCoords(i);
 
    // cP->tt = 1;
    // cP->rr = 1.;
    // cP->R1 = 0;
    // cP->PrintTorusDebug();
    // cP->PrintTorusAnimDebug();
 
    // auto closest_pt_toroid = cP->getClosestPointToToroid(t_d);
    // closestPoint(i) = closest_pt_toroid(i);
 
    // tNormal(j) = t_d(j) - closestPoint(j);
    // if (cP->isPointInsideTorus(t_d))
    //   tNormal(i) *= -1;
    // tNormal.normalize();
 
    // gAp = tNormal(i) * (t_d(i) - closestPoint(i));
 
    // https: // iquilezles.org/www/articles/distfunctions/distfunctions.htm
    const double x = t_d(0);
    const double y = t_d(1);
    const double z = t_d(2);
    const double RR = rAdius - innerRadius;
    const double prod = sqrt(pow(x, 2) + pow(y, 2));
    const double prod2 = sqrt(pow(-RR + prod, 2) + pow(z, 2));
    gAp = -innerRadius + prod2;
    tNormal(0) = (x * (-RR + prod)) / (prod * prod2);
    tNormal(1) = (y * (-RR + prod)) / (prod * prod2);
    tNormal(2) = z / prod2;
    tNormal.normalize();
 
    auto normal_copy = tNormal;
    tNormal(i) = rotationMat(j, i) * normal_copy(j);
 
    return tNormal;
  };

◆ getRollerDataForTag()

MoFEMErrorCode TorusRigidBody::getRollerDataForTag ( )

inlinevirtual

Implements RigidBodyData.

Definition at line 604 of file RigidBodies.hpp.

                                       {
    MoFEMFunctionBeginHot;
    bodyType = TORUS;
    dataForTags[0] = rAdius;
    dataForTags[1] = innerRadius;
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

◆ innerRadius

double TorusRigidBody::innerRadius

Definition at line 597 of file RigidBodies.hpp.

◆ rAdius

double TorusRigidBody::rAdius

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

◆ TorusRigidBody()

Member Function Documentation

◆ getBodyOptions()

◆ getdGap()

◆ getdGapImpl()

◆ getDiffNormal()

◆ getDiffNormalImpl()

◆ getGap()

◆ getGapImpl()

◆ getNormal() [1/3]

◆ getNormal() [2/3]

◆ getNormal() [3/3]

◆ getNormalImpl()

◆ getRollerDataForTag()

Member Data Documentation

◆ innerRadius

◆ rAdius