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

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 314 of file EshelbianADOL-C.cpp.

Constructor & Destructor Documentation

◆ HMHPMooneyRivlinWriggersEq63()

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

inline

Definition at line 319 of file EshelbianADOL-C.cpp.

321 : PhysicalEquations(numberOfActiveVariables, numberOfDependentVariables),

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

Member Function Documentation

◆ getOptions()

MoFEMErrorCode EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::getOptions ( )

inline

Definition at line 331 of file EshelbianADOL-C.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 384 of file EshelbianADOL-C.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();
     // 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 325 of file EshelbianADOL-C.cpp.

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

Member Data Documentation

◆ A

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::A

Definition at line 379 of file EshelbianADOL-C.cpp.

◆ alpha

double EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::alpha

Definition at line 353 of file EshelbianADOL-C.cpp.

◆ B

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::B

Definition at line 380 of file EshelbianADOL-C.cpp.

◆ beta

double EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::beta

Definition at line 354 of file EshelbianADOL-C.cpp.

◆ detF

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::detF

Definition at line 363 of file EshelbianADOL-C.cpp.

◆ detH

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::detH

Definition at line 362 of file EshelbianADOL-C.cpp.

◆ energy

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::energy

Definition at line 377 of file EshelbianADOL-C.cpp.

◆ epsilon

double EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::epsilon

Definition at line 356 of file EshelbianADOL-C.cpp.

◆ lambda

double EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::lambda

Definition at line 355 of file EshelbianADOL-C.cpp.

◆ numberOfActiveVariables

constexpr int EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::numberOfActiveVariables = 9

staticconstexpr

Definition at line 316 of file EshelbianADOL-C.cpp.

◆ numberOfDependentVariables

constexpr int EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::numberOfDependentVariables = 9

staticconstexpr

Definition at line 317 of file EshelbianADOL-C.cpp.

◆ phi

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::phi

Definition at line 378 of file EshelbianADOL-C.cpp.

◆ tBCof

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tBCof

Definition at line 374 of file EshelbianADOL-C.cpp.

◆ tBF

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tBF

Definition at line 373 of file EshelbianADOL-C.cpp.

◆ tBj

adouble EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tBj

Definition at line 375 of file EshelbianADOL-C.cpp.

◆ tCof

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tCof

Definition at line 370 of file EshelbianADOL-C.cpp.

◆ tF

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tF

Definition at line 360 of file EshelbianADOL-C.cpp.

◆ th

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::th

Definition at line 358 of file EshelbianADOL-C.cpp.

◆ tH

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tH

Definition at line 359 of file EshelbianADOL-C.cpp.

◆ tInvF

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tInvF

Definition at line 365 of file EshelbianADOL-C.cpp.

◆ tInvH

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tInvH

Definition at line 364 of file EshelbianADOL-C.cpp.

◆ tP

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tP

Definition at line 367 of file EshelbianADOL-C.cpp.

◆ tPulledP

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tPulledP

Definition at line 382 of file EshelbianADOL-C.cpp.

◆ tSigma

ATensor2 EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::tSigma

Definition at line 368 of file EshelbianADOL-C.cpp.

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

users_modules/eshelbian_plasticit/src/impl/EshelbianADOL-C.cpp

Public Member Functions

Public Attributes

Static Public Attributes

Additional Inherited Members