Inheritance diagram for EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63:

Collaboration diagram for EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63:

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

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 UserDataOperator *	returnOpJacobian (const int tag, const bool eval_rhs, const bool eval_lhs, boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)

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

virtual VolUserDataOperator *	returnOpSpatialPhysicalExternalStrain (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< ExternalStrainVec > external_strain_vec_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv)

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)

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

virtual VolUserDataOperator *	returnOpSetScale (boost::shared_ptr< double > scale_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) {}

EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::epsilon

double epsilon

Definition HMHPMooneyRivlinWriggersEq63.cpp:44

EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::alpha

double alpha

Definition HMHPMooneyRivlinWriggersEq63.cpp:41

EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::numberOfActiveVariables

static constexpr int numberOfActiveVariables

Definition HMHPMooneyRivlinWriggersEq63.cpp:13

EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::numberOfDependentVariables

static constexpr int numberOfDependentVariables

Definition HMHPMooneyRivlinWriggersEq63.cpp:14

EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::lambda

double lambda

Definition HMHPMooneyRivlinWriggersEq63.cpp:43

EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::beta

double beta

Definition HMHPMooneyRivlinWriggersEq63.cpp:42

EshelbianPlasticity::PhysicalEquations::PhysicalEquations

PhysicalEquations()=delete

Member Function Documentation

◆ getOptions()

MoFEMErrorCode EshelbianPlasticity::HMHPMooneyRivlinWriggersEq63::getOptions ( )

inline

Definition at line 21 of file HMHPMooneyRivlinWriggersEq63.cpp.

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

◆ recordTape()

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

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 70 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);
  }