SmallStrainTranverslyIsotropic< TYPE > Struct Template Reference

Hook equation. More...

#include <users_modules/basic_finite_elements/nonlinear_elastic_materials/src/SmallTransverselyIsotropic.hpp>

Inheritance diagram for SmallStrainTranverslyIsotropic< TYPE >:

Collaboration diagram for SmallStrainTranverslyIsotropic< TYPE >:

Public Member Functions
	SmallStrainTranverslyIsotropic ()
MoFEMErrorCode	calculateStrain ()
MoFEMErrorCode	calculateLocalStiffnesMatrix ()
MoFEMErrorCode	calculateAxisAngleRotationalMatrix ()
	Function to Calculate the Rotation Matrix at a given axis and angle of rotation. More...
MoFEMErrorCode	stressTransformation ()
	Function to Calculate Stress Transformation Matrix This function computes the stress transformation Matrix at a give axis and angle of rotation One can also output the axis/angle rotational Matrix. More...
MoFEMErrorCode	strainTransformation ()
	Function to Calculate Strain Transformation Matrix This function computes the strain transformation Matrix at a give axis and angle of rotation One can also output the axis/angle rotational Matrix. More...
MoFEMErrorCode	calculateGlobalStiffnesMatrix ()
virtual MoFEMErrorCode	calculateAngles ()
virtual MoFEMErrorCode	calculateP_PiolaKirchhoffI (const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Calculate global stress. More...
virtual MoFEMErrorCode	calculateElasticEnergy (const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	calculate density of strain energy More...
MoFEMErrorCode	calculateFibreAngles ()
Public Member Functions inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
	FunctionsToCalculatePiolaKirchhoffI ()
virtual	~FunctionsToCalculatePiolaKirchhoffI ()=default
MoFEMErrorCode	calculateC_CauchyDeformationTensor ()
MoFEMErrorCode	calculateE_GreenStrain ()
MoFEMErrorCode	calculateS_PiolaKirchhoffII ()
virtual MoFEMErrorCode	calculateP_PiolaKirchhoffI (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Function overload to implement user material. More...
virtual MoFEMErrorCode	calculateCauchyStress (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Function overload to implement user material. More...
virtual MoFEMErrorCode	setUserActiveVariables (int &nb_active_variables)
	add additional active variables More...
virtual MoFEMErrorCode	setUserActiveVariables (VectorDouble &activeVariables)
	Add additional independent variables More complex physical models depend on gradient of defamation and some additional variables. For example can depend on temperature. This function adds additional independent variables to the model. More...
virtual MoFEMErrorCode	calculateElasticEnergy (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Calculate elastic energy density. More...
virtual MoFEMErrorCode	calculatesIGma_EshelbyStress (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Calculate Eshelby stress. More...
virtual MoFEMErrorCode	getDataOnPostProcessor (std::map< std::string, std::vector< VectorDouble > > &field_map, std::map< std::string, std::vector< MatrixDouble > > &grad_map)
	Do operations when pre-process. More...

Public Attributes
ublas::matrix< TYPE >	sTrain
ublas::vector< TYPE >	voightStrain
TYPE	tR
double	nu_p
double	nu_pz
double	E_p
double	E_z
double	G_zp
ublas::symmetric_matrix< TYPE, ublas::upper >	localStiffnessMatrix
ublas::matrix< TYPE >	aARotMat
ublas::vector< TYPE >	axVector
TYPE	axAngle
ublas::matrix< TYPE >	stressRotMat
ublas::matrix< TYPE >	strainRotMat
ublas::matrix< TYPE >	dR
ublas::matrix< TYPE >	globalStiffnessMatrix
ublas::vector< TYPE >	voigtStress
VectorDouble	normalizedPhi
VectorDouble	axVectorDouble
double	axAngleDouble
Public Attributes inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
FTensor::Index< 'i', 3 >	i
FTensor::Index< 'j', 3 >	j
FTensor::Index< 'k', 3 >	k
double	lambda
double	mu
MatrixBoundedArray< TYPE, 9 >	F
MatrixBoundedArray< TYPE, 9 >	C
MatrixBoundedArray< TYPE, 9 >	E
MatrixBoundedArray< TYPE, 9 >	S
MatrixBoundedArray< TYPE, 9 >	invF
MatrixBoundedArray< TYPE, 9 >	P
MatrixBoundedArray< TYPE, 9 >	sIGma
MatrixBoundedArray< TYPE, 9 >	h
MatrixBoundedArray< TYPE, 9 >	H
MatrixBoundedArray< TYPE, 9 >	invH
MatrixBoundedArray< TYPE, 9 >	sigmaCauchy
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_F
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_C
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_E
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_S
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_invF
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_P
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_sIGma
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_h
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_H
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_invH
FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_sigmaCauchy
TYPE	J
TYPE	eNergy
TYPE	detH
TYPE	detF
int	gG
	Gauss point number. More...
CommonData *	commonDataPtr
MoFEM::VolumeElementForcesAndSourcesCore::UserDataOperator *	opPtr
	pointer to finite element tetrahedral operator More...

Additional Inherited Members
Static Protected Member Functions inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
static auto	resizeAndSet (MatrixBoundedArray< TYPE, 9 > &m)

Detailed Description

template<typename TYPE>
struct SmallStrainTranverslyIsotropic< TYPE >

Hook equation.

Definition at line 19 of file SmallTransverselyIsotropic.hpp.

Constructor & Destructor Documentation

SmallStrainTranverslyIsotropic()

template<typename TYPE >

SmallStrainTranverslyIsotropic< TYPE >::SmallStrainTranverslyIsotropic ( )

inline

Definition at line 21 of file SmallTransverselyIsotropic.hpp.

: NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI<TYPE>() {}

Member Function Documentation

calculateAngles()

template<typename TYPE >

virtual MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::calculateAngles ( )

inlinevirtual

Reimplemented in SmallStrainTranverslyIsotropicDouble.

Definition at line 236 of file SmallTransverselyIsotropic.hpp.

                                           {
    MoFEMFunctionBeginHot;
    MoFEMFunctionReturnHot(0);
  }

calculateAxisAngleRotationalMatrix()

template<typename TYPE >

MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::calculateAxisAngleRotationalMatrix ( )

inline

Function to Calculate the Rotation Matrix at a given axis and angle of rotation.

This function computes the rotational matrix for a given axis of rotation and angle of rotation about that angle

Parameters

axVector	A vector representing the axis of rotation
axAngle	Angle of rotation along the axis (in radians)

Definition at line 82 of file SmallTransverselyIsotropic.hpp.

                                                      {
    MoFEMFunctionBeginHot;
 
    aARotMat.resize(3,3,false);
    aARotMat.clear();
 
    TYPE norm = sqrt(pow(axVector[0],2) + pow(axVector[1],2) + pow(axVector[2],2));
 
    aARotMat(0,0) = 1-((1-cos(axAngle))*(pow(axVector[1],2)+pow(axVector[2],2))/pow(norm,2));
    aARotMat(1,1) = 1-((1-cos(axAngle))*(pow(axVector[0],2)+pow(axVector[2],2))/pow(norm,2));
    aARotMat(2,2) = 1-((1-cos(axAngle))*(pow(axVector[0],2)+pow(axVector[1],2))/pow(norm,2));
 
    aARotMat(0,1) = ((1-cos(axAngle))*axVector[0]*axVector[1]-norm*axVector[2]*sin(axAngle))/pow(norm,2);
    aARotMat(1,0) = ((1-cos(axAngle))*axVector[0]*axVector[1]+norm*axVector[2]*sin(axAngle))/pow(norm,2);
 
    aARotMat(0,2) = ((1-cos(axAngle))*axVector[0]*axVector[2]+norm*axVector[1]*sin(axAngle))/pow(norm,2);
    aARotMat(2,0) = ((1-cos(axAngle))*axVector[0]*axVector[2]-norm*axVector[1]*sin(axAngle))/pow(norm,2);
 
    aARotMat(1,2) = ((1-cos(axAngle))*axVector[1]*axVector[2]-norm*axVector[0]*sin(axAngle))/pow(norm,2);
    aARotMat(2,1) = ((1-cos(axAngle))*axVector[1]*axVector[2]+norm*axVector[0]*sin(axAngle))/pow(norm,2);
 
    MoFEMFunctionReturnHot(0);
  }

calculateElasticEnergy()

template<typename TYPE >

virtual MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::calculateElasticEnergy ( const NonlinearElasticElement::BlockData  block_data, boost::shared_ptr< const NumeredEntFiniteElement >  fe_ptr  )

inlinevirtual

calculate density of strain energy

Reimplemented from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >.

Definition at line 284 of file SmallTransverselyIsotropic.hpp.

  {
    MoFEMFunctionBeginHot;
 
    ierr = calculateAngles(); CHKERRG(ierr);
    ierr = calculateStrain(); CHKERRG(ierr);
    ierr = calculateLocalStiffnesMatrix(); CHKERRG(ierr);
    ierr = calculateAxisAngleRotationalMatrix(); CHKERRG(ierr);
    ierr = stressTransformation(); CHKERRG(ierr);
    axAngle = -axAngle;
    ierr = calculateAxisAngleRotationalMatrix(); CHKERRG(ierr);
    ierr = strainTransformation(); CHKERRG(ierr);
    ierr = calculateGlobalStiffnesMatrix(); CHKERRG(ierr);
 
    voigtStress.resize(6,false);
    noalias(voigtStress) = prod(globalStiffnessMatrix,voightStrain);
    this->eNergy = 0.5*inner_prod(voigtStress,voightStrain);
    MoFEMFunctionReturnHot(0);
  }

calculateFibreAngles()

template<typename TYPE >

MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::calculateFibreAngles ( )

inline

Definition at line 310 of file SmallTransverselyIsotropic.hpp.

                                        {
    MoFEMFunctionBeginHot;
 
    try {
 
      int gg = this->gG; // number of integration point
      MatrixDouble &phi = (this->commonDataPtr->gradAtGaussPts["POTENTIAL_FIELD"][gg]);
      normalizedPhi.resize(3,false);
      double nrm2_phi = sqrt(pow(phi(0,0),2)+pow(phi(0,1),2)+pow(phi(0,2),2));
      for(int ii = 0;ii<3;ii++) {
        normalizedPhi[ii] = -phi(0,ii)/nrm2_phi;
      }
 
      axVectorDouble.resize(3,false);
      const double zVec[3]={ 0.0,0.0,1.0 };
      axVectorDouble[0] = normalizedPhi[1]*zVec[2]-normalizedPhi[2]*zVec[1];
      axVectorDouble[1] = normalizedPhi[2]*zVec[0]-normalizedPhi[0]*zVec[2];
      axVectorDouble[2] = normalizedPhi[0]*zVec[1]-normalizedPhi[1]*zVec[0];
      double nrm2_ax_vector = norm_2(axVectorDouble);
      const double eps = 1e-12;
      if(nrm2_ax_vector<eps) {
        axVectorDouble[0] = 0;
        axVectorDouble[1] = 0;
        axVectorDouble[2] = 1;
        nrm2_ax_vector = 1;
      }
      axAngleDouble = asin(nrm2_ax_vector);
 
    } catch (const std::exception& ex) {
      std::ostringstream ss;
      ss << "throw in method: " << ex.what() << std::endl;
      SETERRQ(PETSC_COMM_SELF,1,ss.str().c_str());
    }
 
    MoFEMFunctionReturnHot(0);
  }

calculateGlobalStiffnesMatrix()

template<typename TYPE >

MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::calculateGlobalStiffnesMatrix ( )

inline

Definition at line 225 of file SmallTransverselyIsotropic.hpp.

                                                 {
    MoFEMFunctionBeginHot;
 
    dR.resize(6,6,false);
    noalias(dR) = prod(localStiffnessMatrix,strainRotMat);
    globalStiffnessMatrix.resize(6,6,false);
    noalias(globalStiffnessMatrix) = prod(stressRotMat,dR);
 
    MoFEMFunctionReturnHot(0);
  }

calculateLocalStiffnesMatrix()

template<typename TYPE >

MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::calculateLocalStiffnesMatrix ( )

inline

Definition at line 50 of file SmallTransverselyIsotropic.hpp.

                                                {
    MoFEMFunctionBeginHot;
    double nu_zp=(nu_pz*E_z)/E_p;
    double delta=((1+nu_p)*(1-nu_p-(2*nu_pz*nu_zp)))/(E_p*E_p*E_z);
 
    localStiffnessMatrix.resize(6);
    localStiffnessMatrix.clear();
    localStiffnessMatrix(0,0)=localStiffnessMatrix(1,1)=(1-nu_pz*nu_zp)/(E_p*E_z*delta);
    localStiffnessMatrix(2,2)=(1-nu_p*nu_p)/(E_p*E_p*delta);
 
    localStiffnessMatrix(0,1)=localStiffnessMatrix(1,0)=(nu_p+nu_zp*nu_pz)/(E_p*E_z*delta);
    localStiffnessMatrix(0,2)=localStiffnessMatrix(2,0)=localStiffnessMatrix(1,2)=localStiffnessMatrix(2,1)=(nu_zp+nu_p*nu_zp)/(E_p*E_z*delta);
 
    localStiffnessMatrix(3,3)=E_p/(2*(1+nu_p));
    localStiffnessMatrix(4,4)=localStiffnessMatrix(5,5)=G_zp;
 
    MoFEMFunctionReturnHot(0);
  }

calculateP_PiolaKirchhoffI()

template<typename TYPE >

virtual MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::calculateP_PiolaKirchhoffI ( const NonlinearElasticElement::BlockData  block_data, boost::shared_ptr< const NumeredEntFiniteElement >  fe_ptr  )

inlinevirtual

Calculate global stress.

This is small strain approach, i.e. Piola stress is like a Cauchy stress, since configurations are notation distinguished.

Reimplemented from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >.

Definition at line 250 of file SmallTransverselyIsotropic.hpp.

    {
    MoFEMFunctionBeginHot;
    ierr = calculateAngles(); CHKERRG(ierr);
    ierr = calculateStrain(); CHKERRG(ierr);
    ierr = calculateLocalStiffnesMatrix(); CHKERRG(ierr);
    ierr = calculateAxisAngleRotationalMatrix(); CHKERRG(ierr);
    ierr = stressTransformation(); CHKERRG(ierr);
    axAngle = -axAngle;
    ierr = calculateAxisAngleRotationalMatrix(); CHKERRG(ierr);
    ierr = strainTransformation(); CHKERRG(ierr);
    ierr = calculateGlobalStiffnesMatrix(); CHKERRG(ierr);
 
    voigtStress.resize(6,false);
    noalias(voigtStress) = prod(globalStiffnessMatrix,voightStrain);
    this->P.resize(3,3,false);
    this->P(0,0) = voigtStress[0];
    this->P(1,1) = voigtStress[1];
    this->P(2,2) = voigtStress[2];
    this->P(0,1) = voigtStress[3];
    this->P(1,2) = voigtStress[4];
    this->P(0,2) = voigtStress[5];
    this->P(1,0) = this->P(0,1);
    this->P(2,1) = this->P(1,2);
    this->P(2,0) = this->P(0,2);
 
    MoFEMFunctionReturnHot(0);
  }

calculateStrain()

template<typename TYPE >

MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::calculateStrain ( )

inline

Definition at line 29 of file SmallTransverselyIsotropic.hpp.

                                   {
    MoFEMFunctionBeginHot;
    sTrain.resize(3,3,false);
    noalias(sTrain) = this->F;
    for(int dd = 0;dd<3;dd++) {
      sTrain(dd,dd) -= 1;
    }
    sTrain += trans(sTrain);
    sTrain *= 0.5;
    voightStrain.resize(6,false);
    voightStrain[0] = sTrain(0,0);
    voightStrain[1] = sTrain(1,1);
    voightStrain[2] = sTrain(2,2);
    voightStrain[3] = 2*sTrain(0,1);
    voightStrain[4] = 2*sTrain(1,2);
    voightStrain[5] = 2*sTrain(2,0);
    MoFEMFunctionReturnHot(0);
  }

strainTransformation()

template<typename TYPE >

MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::strainTransformation ( )

inline

Function to Calculate Strain Transformation Matrix
This function computes the strain transformation Matrix at a give axis and angle of rotation
One can also output the axis/angle rotational Matrix.

Definition at line 171 of file SmallTransverselyIsotropic.hpp.

                                        {
    MoFEMFunctionBeginHot;
 
    strainRotMat.resize(6,6,false);
    strainRotMat.clear();
 
    strainRotMat(0, 0) = aARotMat(0,0) * aARotMat(0,0);
    strainRotMat(0, 1) = aARotMat(1,0) * aARotMat(1,0);
    strainRotMat(0, 2) = aARotMat(2,0) * aARotMat(2,0);
    strainRotMat(0, 3) = aARotMat(1,0) * aARotMat(0,0);
    strainRotMat(0, 4) = aARotMat(2,0) * aARotMat(1,0);
    strainRotMat(0, 5) = aARotMat(0,0) * aARotMat(2,0);
 
    strainRotMat(1, 0) = aARotMat(0,1) * aARotMat(0,1);
    strainRotMat(1, 1) = aARotMat(1,1) * aARotMat(1,1);
    strainRotMat(1, 2) = aARotMat(2,1) * aARotMat(2,1);
    strainRotMat(1, 3) = aARotMat(1,1) * aARotMat(0,1);
    strainRotMat(1, 4) = aARotMat(2,1) * aARotMat(1,1);
    strainRotMat(1, 5) = aARotMat(0,1) * aARotMat(2,1);
 
    strainRotMat(2, 0) = aARotMat(0,2) * aARotMat(0,2);
    strainRotMat(2, 1) = aARotMat(1,2) * aARotMat(1,2);
    strainRotMat(2, 2) = aARotMat(2,2) * aARotMat(2,2);
    strainRotMat(2, 3) = aARotMat(1,2) * aARotMat(0,2);
    strainRotMat(2, 4) = aARotMat(2,2) * aARotMat(1,2);
    strainRotMat(2, 5) = aARotMat(0,2) * aARotMat(2,2);
 
    strainRotMat(3, 0) = 2.0 * aARotMat(0,1) * aARotMat(0,0);
    strainRotMat(3, 1) = 2.0 * aARotMat(1,1) * aARotMat(1,0);
    strainRotMat(3, 2) = 2.0 * aARotMat(2,1) * aARotMat(2,0);
    strainRotMat(3, 3) =     ( aARotMat(1,1) * aARotMat(0,0) + aARotMat(0,1) * aARotMat(1,0) );
    strainRotMat(3, 4) =     ( aARotMat(2,1) * aARotMat(1,0) + aARotMat(1,1) * aARotMat(2,0) );
    strainRotMat(3, 5) =     ( aARotMat(0,1) * aARotMat(2,0) + aARotMat(2,1) * aARotMat(0,0) );
 
    strainRotMat(4, 0) = 2.0 * aARotMat(0,2) * aARotMat(0,1);
    strainRotMat(4, 1) = 2.0 * aARotMat(1,2) * aARotMat(1,1);
    strainRotMat(4, 2) = 2.0 * aARotMat(2,2) * aARotMat(2,1);
    strainRotMat(4, 3) =     ( aARotMat(1,2) * aARotMat(0,1) + aARotMat(0,2) * aARotMat(1,1) );
    strainRotMat(4, 4) =     ( aARotMat(2,2) * aARotMat(1,1) + aARotMat(1,2) * aARotMat(2,1) );
    strainRotMat(4, 5) =     ( aARotMat(0,2) * aARotMat(2,1) + aARotMat(2,2) * aARotMat(0,1) );
 
    strainRotMat(5, 0) = 2.0 * aARotMat(0,0) * aARotMat(0,2);
    strainRotMat(5, 1) = 2.0 * aARotMat(1,0) * aARotMat(1,2);
    strainRotMat(5, 2) = 2.0 * aARotMat(2,0) * aARotMat(2,2);
    strainRotMat(5, 3) =     ( aARotMat(1,0) * aARotMat(0,2) + aARotMat(0,0) * aARotMat(1,2) );
    strainRotMat(5, 4) =     ( aARotMat(2,0) * aARotMat(1,2) + aARotMat(1,0) * aARotMat(2,2) );
    strainRotMat(5, 5) =     ( aARotMat(0,0) * aARotMat(2,2) + aARotMat(2,0) * aARotMat(0,2) );
 
    MoFEMFunctionReturnHot(0);
  }

stressTransformation()

template<typename TYPE >

MoFEMErrorCode SmallStrainTranverslyIsotropic< TYPE >::stressTransformation ( )

inline

Function to Calculate Stress Transformation Matrix This function computes the stress transformation Matrix at a give axis and angle of rotation
One can also output the axis/angle rotational Matrix.

Definition at line 113 of file SmallTransverselyIsotropic.hpp.

                                        {
    MoFEMFunctionBeginHot;
 
    stressRotMat.resize(6,6,false);
    stressRotMat.clear();
 
    stressRotMat(0, 0) =       aARotMat(0,0) * aARotMat(0,0);
    stressRotMat(0, 1) =       aARotMat(1,0) * aARotMat(1,0);
    stressRotMat(0, 2) =       aARotMat(2,0) * aARotMat(2,0);
    stressRotMat(0, 3) = 2.0 * aARotMat(1,0) * aARotMat(0,0);
    stressRotMat(0, 4) = 2.0 * aARotMat(2,0) * aARotMat(1,0);
    stressRotMat(0, 5) = 2.0 * aARotMat(0,0) * aARotMat(2,0);
 
    stressRotMat(1, 0) =       aARotMat(0,1) * aARotMat(0,1);
    stressRotMat(1, 1) =       aARotMat(1,1) * aARotMat(1,1);
    stressRotMat(1, 2) =       aARotMat(2,1) * aARotMat(2,1);
    stressRotMat(1, 3) = 2.0 * aARotMat(1,1) * aARotMat(0,1);
    stressRotMat(1, 4) = 2.0 * aARotMat(2,1) * aARotMat(1,1);
    stressRotMat(1, 5) = 2.0 * aARotMat(0,1) * aARotMat(2,1);
 
    stressRotMat(2, 0) =       aARotMat(0,2) * aARotMat(0,2);
    stressRotMat(2, 1) =       aARotMat(1,2) * aARotMat(1,2);
    stressRotMat(2, 2) =       aARotMat(2,2) * aARotMat(2,2);
    stressRotMat(2, 3) = 2.0 * aARotMat(1,2) * aARotMat(0,2);
    stressRotMat(2, 4) = 2.0 * aARotMat(2,2) * aARotMat(1,2);
    stressRotMat(2, 5) = 2.0 * aARotMat(0,2) * aARotMat(2,2);
 
    stressRotMat(3, 0) =   aARotMat(0,1) * aARotMat(0,0);
    stressRotMat(3, 1) =   aARotMat(1,1) * aARotMat(1,0);
    stressRotMat(3, 2) =   aARotMat(2,1) * aARotMat(2,0);
    stressRotMat(3, 3) = ( aARotMat(1,1) * aARotMat(0,0) + aARotMat(0,1) * aARotMat(1,0) );
    stressRotMat(3, 4) = ( aARotMat(2,1) * aARotMat(1,0) + aARotMat(1,1) * aARotMat(2,0) );
    stressRotMat(3, 5) = ( aARotMat(0,1) * aARotMat(2,0) + aARotMat(2,1) * aARotMat(0,0) );
 
    stressRotMat(4, 0) =   aARotMat(0,2) * aARotMat(0,1);
    stressRotMat(4, 1) =   aARotMat(1,2) * aARotMat(1,1);
    stressRotMat(4, 2) =   aARotMat(2,2) * aARotMat(2,1);
    stressRotMat(4, 3) = ( aARotMat(1,2) * aARotMat(0,1) + aARotMat(0,2) * aARotMat(1,1) );
    stressRotMat(4, 4) = ( aARotMat(2,2) * aARotMat(1,1) + aARotMat(1,2) * aARotMat(2,1) );
    stressRotMat(4, 5) = ( aARotMat(0,2) * aARotMat(2,1) + aARotMat(2,2) * aARotMat(0,1) );
 
    stressRotMat(5, 0) =   aARotMat(0,0) * aARotMat(0,2);
    stressRotMat(5, 1) =   aARotMat(1,0) * aARotMat(1,2);
    stressRotMat(5, 2) =   aARotMat(2,0) * aARotMat(2,2);
    stressRotMat(5, 3) = ( aARotMat(1,0) * aARotMat(0,2) + aARotMat(0,0) * aARotMat(1,2) );
    stressRotMat(5, 4) = ( aARotMat(2,0) * aARotMat(1,2) + aARotMat(1,0) * aARotMat(2,2) );
    stressRotMat(5, 5) = ( aARotMat(0,0) * aARotMat(2,2) + aARotMat(2,0) * aARotMat(0,2) );
 
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

aARotMat

template<typename TYPE >

ublas::matrix<TYPE> SmallStrainTranverslyIsotropic< TYPE >::aARotMat

Definition at line 69 of file SmallTransverselyIsotropic.hpp.

axAngle

template<typename TYPE >

TYPE SmallStrainTranverslyIsotropic< TYPE >::axAngle

Definition at line 71 of file SmallTransverselyIsotropic.hpp.

axAngleDouble

template<typename TYPE >

double SmallStrainTranverslyIsotropic< TYPE >::axAngleDouble

Definition at line 308 of file SmallTransverselyIsotropic.hpp.

axVector

template<typename TYPE >

ublas::vector<TYPE> SmallStrainTranverslyIsotropic< TYPE >::axVector

Definition at line 70 of file SmallTransverselyIsotropic.hpp.

axVectorDouble

template<typename TYPE >

VectorDouble SmallStrainTranverslyIsotropic< TYPE >::axVectorDouble

Definition at line 307 of file SmallTransverselyIsotropic.hpp.

dR

template<typename TYPE >

ublas::matrix<TYPE> SmallStrainTranverslyIsotropic< TYPE >::dR

Definition at line 222 of file SmallTransverselyIsotropic.hpp.

E_p

template<typename TYPE >

double SmallStrainTranverslyIsotropic< TYPE >::E_p

Definition at line 48 of file SmallTransverselyIsotropic.hpp.

E_z

template<typename TYPE >

double SmallStrainTranverslyIsotropic< TYPE >::E_z

Definition at line 48 of file SmallTransverselyIsotropic.hpp.

G_zp

template<typename TYPE >

double SmallStrainTranverslyIsotropic< TYPE >::G_zp

Definition at line 48 of file SmallTransverselyIsotropic.hpp.

globalStiffnessMatrix

template<typename TYPE >

ublas::matrix<TYPE> SmallStrainTranverslyIsotropic< TYPE >::globalStiffnessMatrix

Definition at line 223 of file SmallTransverselyIsotropic.hpp.

localStiffnessMatrix

template<typename TYPE >

ublas::symmetric_matrix<TYPE,ublas::upper> SmallStrainTranverslyIsotropic< TYPE >::localStiffnessMatrix

Definition at line 49 of file SmallTransverselyIsotropic.hpp.

normalizedPhi

template<typename TYPE >

VectorDouble SmallStrainTranverslyIsotropic< TYPE >::normalizedPhi

Definition at line 306 of file SmallTransverselyIsotropic.hpp.

nu_p

template<typename TYPE >

double SmallStrainTranverslyIsotropic< TYPE >::nu_p

Definition at line 48 of file SmallTransverselyIsotropic.hpp.

nu_pz

template<typename TYPE >

double SmallStrainTranverslyIsotropic< TYPE >::nu_pz

Definition at line 48 of file SmallTransverselyIsotropic.hpp.

sTrain

template<typename TYPE >

ublas::matrix<TYPE> SmallStrainTranverslyIsotropic< TYPE >::sTrain

Definition at line 25 of file SmallTransverselyIsotropic.hpp.

strainRotMat

template<typename TYPE >

ublas::matrix<TYPE> SmallStrainTranverslyIsotropic< TYPE >::strainRotMat

Definition at line 164 of file SmallTransverselyIsotropic.hpp.

stressRotMat

template<typename TYPE >

ublas::matrix<TYPE> SmallStrainTranverslyIsotropic< TYPE >::stressRotMat

Definition at line 106 of file SmallTransverselyIsotropic.hpp.

tR

template<typename TYPE >

TYPE SmallStrainTranverslyIsotropic< TYPE >::tR

Definition at line 27 of file SmallTransverselyIsotropic.hpp.

voightStrain

template<typename TYPE >

ublas::vector<TYPE> SmallStrainTranverslyIsotropic< TYPE >::voightStrain

Definition at line 26 of file SmallTransverselyIsotropic.hpp.

voigtStress

template<typename TYPE >

ublas::vector<TYPE> SmallStrainTranverslyIsotropic< TYPE >::voigtStress

Definition at line 241 of file SmallTransverselyIsotropic.hpp.

---

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

• users_modules/basic_finite_elements/nonlinear_elastic_materials/src/SmallTransverselyIsotropic.hpp