Inheritance diagram for EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63:

Collaboration diagram for EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63:

Public Member Functions
	HMHPMooneyRivlinWriggersEq63 (const double alpha, const double beta, const double lambda)

virtual OpJacobian *	returnOpJacobian (const int tag, const bool eval_rhs, const bool eval_lhs, boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)

MoFEMErrorCode	getOptions ()

MoFEMErrorCode	recordTape (const int tape, DTensor2Ptr *t_h_ptr)

Public Member Functions inherited from EshelbianPlasticity::PhysicalEquations
	PhysicalEquations ()=delete

	PhysicalEquations (const int size_active, const int size_dependent)

virtual	~PhysicalEquations ()=default

virtual VolUserDataOperator *	returnOpSpatialPhysical (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const double alpha_u)

virtual VolUserDataOperator *	returnOpSpatialPhysical_du_du (std::string row_field, std::string col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const double alpha)

virtual VolUserDataOperator *	returnOpCalculateEnergy (boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< double > total_energy_ptr)

virtual VolUserDataOperator *	returnOpCalculateStretchFromStress (boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)

DTensor2Ptr	get_P ()

DTensor3Ptr	get_P_dh0 ()

DTensor3Ptr	get_P_dh1 ()

DTensor3Ptr	get_P_dh2 ()

DTensor2Ptr	get_h ()

Public Attributes
double	alpha

double	beta

double	lambda

double	epsilon

ATensor2	th

ATensor2	tH

ATensor2	tF

adouble	detH

adouble	detF

ATensor2	tInvH

ATensor2	tInvF

ATensor2	tP

ATensor2	tSigma

ATensor2	tCof

ATensor2	tBF

ATensor2	tBCof

adouble	tBj

adouble	energy

adouble	phi

adouble	A

adouble	B

ATensor2	tPulledP

Public Attributes inherited from EshelbianPlasticity::PhysicalEquations
std::vector< double >	activeVariables

std::vector< double >	dependentVariablesPiola

std::vector< double >	dependentVariablesPiolaDirevatives

Static Public Attributes
static constexpr int	numberOfActiveVariables = 9

static constexpr int	numberOfDependentVariables = 9

Additional Inherited Members
Public Types inherited from EshelbianPlasticity::PhysicalEquations
typedef FTensor::Tensor1< adouble, 3 >	ATensor1

typedef FTensor::Tensor2< adouble, 3, 3 >	ATensor2

typedef FTensor::Tensor3< adouble, 3, 3, 3 >	ATensor3

typedef FTensor::Tensor1< double, 3 >	DTensor1

typedef FTensor::Tensor2< double, 3, 3 >	DTensor2

typedef FTensor::Tensor3< double, 3, 3, 3 >	DTensor3

typedef FTensor::Tensor0< FTensor::PackPtr< double *, 1 > >	DTensor0Ptr

typedef FTensor::Tensor2< FTensor::PackPtr< double *, 1 >, 3, 3 >	DTensor2Ptr

typedef FTensor::Tensor3< FTensor::PackPtr< double *, 1 >, 3, 3, 3 >	DTensor3Ptr

Static Public Member Functions inherited from EshelbianPlasticity::PhysicalEquations
template<int S>
static DTensor2Ptr	get_VecTensor2 (std::vector< double > &v)

template<int S>
static DTensor0Ptr	get_VecTensor0 (std::vector< double > &v)

template<int S0>
static DTensor3Ptr	get_vecTensor3 (std::vector< double > &v, const int nba)

Detailed Description

Definition at line 11 of file HMHPMooneyRivlinWriggersEq63.cpp.

Constructor & Destructor Documentation

◆ HMHPMooneyRivlinWriggersEq63()

EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::HMHPMooneyRivlinWriggersEq63	(	const double	alpha,
		const double	beta,
		const double	lambda
	)

inline

Definition at line 16 of file HMHPMooneyRivlinWriggersEq63.cpp.

18 : PhysicalEquations(numberOfActiveVariables, numberOfDependentVariables),

19 alpha(alpha), beta(beta), lambda(lambda), epsilon(0) {}

Member Function Documentation

◆ getOptions()

MoFEMErrorCode EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::getOptions ( )

inline

Definition at line 28 of file HMHPMooneyRivlinWriggersEq63.cpp.

                               {
     MoFEMFunctionBegin;
     CHKERR PetscOptionsBegin(PETSC_COMM_WORLD, "mooneyrivlin_", "", "none");
  
     alpha = 1;
     CHKERR PetscOptionsScalar("-alpha", "Alpha", "", alpha, &alpha, PETSC_NULL);
     beta = 1;
     CHKERR PetscOptionsScalar("-beta", "Beta", "", beta, &beta, PETSC_NULL);
  
     lambda = 1;
     CHKERR PetscOptionsScalar("-lambda", "Lambda", "", lambda, &lambda,
                               PETSC_NULL);
  
     epsilon = 0;
     CHKERR PetscOptionsScalar("-epsilon", "Epsilon", "", epsilon, &epsilon,
                               PETSC_NULL);
  
     ierr = PetscOptionsEnd();
     CHKERRG(ierr);
     MoFEMFunctionReturn(0);
   }

◆ recordTape()

MoFEMErrorCode EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::recordTape	(	const int	tape,
		DTensor2Ptr *	t_h_ptr
	)

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 79 of file HMHPMooneyRivlinWriggersEq63.cpp.

                                                                   {
     MoFEMFunctionBegin;
  
     FTensor::Index<'i', 3> i;
     FTensor::Index<'j', 3> j;
     FTensor::Index<'k', 3> k;
     FTensor::Index<'I', 3> I;
     FTensor::Index<'J', 3> J;
     FTensor::Index<'K', 3> K;
     FTensor::Number<0> N0;
     FTensor::Number<1> N1;
     FTensor::Number<2> N2;
  
     CHKERR getOptions();
  
     auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
  
     auto ih = get_h();
     if (t_h_ptr)
       ih(i, j) = (*t_h_ptr)(i, j);
     else {
       ih(i, j) = t_kd(i, j);
     }
  
     auto r_P = get_P();
  
     enableMinMaxUsingAbs();
     trace_on(tape);
  
     // Set active variables to ADOL-C
     th(i, j) <<= get_h()(i, j);
  
     tF(i, I) = th(i, I);
     CHKERR determinantTensor3by3(tF, detF);
     CHKERR invertTensor3by3(tF, detF, tInvF);
  
     tCof(i, I) = detF * tInvF(I, i);
  
     A = tF(k, K) * tF(k, K);
     B = tCof(k, K) * tCof(k, K);
  
     tBF(i, I) = 4 * alpha * (A * tF(i, I));
     tBCof(i, I) = 4 * beta * (B * tCof(i, I));
     tBj = (-12 * alpha - 24 * beta) / detF +
           0.5 * (lambda / epsilon) *
               (pow(detF, epsilon - 1) - pow(detF, -epsilon - 1));
  
     tP(i, I) = tBF(i, I);
     tP(i, I) += (levi_civita(i, j, k) * tBCof(j, J)) *
                 (levi_civita(I, J, K) * tF(k, K));
     tP(i, I) += tCof(i, I) * tBj;
  
     // Set dependent variables to ADOL-C
     tP(i, j) >>= r_P(i, j);
  
     trace_off();
  
     MoFEMFunctionReturn(0);
   }

◆ returnOpJacobian()

virtual OpJacobian* EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::returnOpJacobian	(	const int	tag,
		const bool	eval_rhs,
		const bool	eval_lhs,
		boost::shared_ptr< DataAtIntegrationPts >	data_ptr,
		boost::shared_ptr< PhysicalEquations >	physics_ptr
	)

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 22 of file HMHPMooneyRivlinWriggersEq63.cpp.

                                                                    {
     return (new OpHMHH(tag, eval_rhs, eval_lhs, data_ptr, physics_ptr));
   }

Member Data Documentation

◆ A

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::A

Definition at line 74 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ alpha

double EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::alpha

Definition at line 50 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ B

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::B

Definition at line 75 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ beta

double EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::beta

Definition at line 51 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ detF

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::detF

Definition at line 60 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ detH

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::detH

Definition at line 59 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ energy

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::energy

Definition at line 72 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ epsilon

double EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::epsilon

Definition at line 53 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ lambda

double EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::lambda

Definition at line 52 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ numberOfActiveVariables

constexpr int EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::numberOfActiveVariables = 9

staticconstexpr

Definition at line 13 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ numberOfDependentVariables

constexpr int EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::numberOfDependentVariables = 9

staticconstexpr

Definition at line 14 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ phi

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::phi

Definition at line 73 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tBCof

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tBCof

Definition at line 69 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tBF

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tBF

Definition at line 68 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tBj

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tBj

Definition at line 70 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tCof

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tCof

Definition at line 67 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tF

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tF

Definition at line 57 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ th

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::th

Definition at line 55 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tH

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tH

Definition at line 56 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tInvF

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tInvF

Definition at line 62 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tInvH

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tInvH

Definition at line 61 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tP

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tP

Definition at line 64 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tPulledP

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tPulledP

Definition at line 77 of file HMHPMooneyRivlinWriggersEq63.cpp.

◆ tSigma

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tSigma

Definition at line 65 of file HMHPMooneyRivlinWriggersEq63.cpp.

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

users_modules/eshelbian_plasticit/src/impl/HMHPMooneyRivlinWriggersEq63.cpp

Public Member Functions

Public Attributes

Static Public Attributes

Additional Inherited Members