EshelbianPlasticity::HMHStVenantKirchhoff Struct Reference

Inheritance diagram for EshelbianPlasticity::HMHStVenantKirchhoff:

Collaboration diagram for EshelbianPlasticity::HMHStVenantKirchhoff:

Public Member Functions
	HMHStVenantKirchhoff (const double lambda, const double mu)
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	lambda
double	mu
double	sigmaY
ATensor2	th
ATensor2	tH
ATensor2	tP
ATensor2	tSigma
adouble	phi
adouble	s2
adouble	f
adouble	energy
adouble	detH
adouble	detF
adouble	trE
ATensor2	tInvH
ATensor2	tF
ATensor2	tInvF
ATensor2	tC
ATensor2	tE
ATensor2	tS
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 13 of file HMHHStVenantKirchhoff.cpp.

Constructor & Destructor Documentation

HMHStVenantKirchhoff()

EshelbianPlasticity::HMHStVenantKirchhoff::HMHStVenantKirchhoff ( const double lambda, const double mu )

inline

Definition at line 18 of file HMHHStVenantKirchhoff.cpp.

      : PhysicalEquations(numberOfActiveVariables, numberOfDependentVariables),
        lambda(lambda), mu(mu) {}

Member Function Documentation

getOptions()

MoFEMErrorCode EshelbianPlasticity::HMHStVenantKirchhoff::getOptions ( )

inline

Definition at line 22 of file HMHHStVenantKirchhoff.cpp.

                              {
    MoFEMFunctionBegin;
 
    double E = 1;
    double nu = 0;
 
    PetscOptionsBegin(PETSC_COMM_WORLD, "stvenant_", "", "none");
 
    CHKERR PetscOptionsScalar("-young_modulus", "Young modulus", "", E, &E,
                              PETSC_NULLPTR);
    CHKERR PetscOptionsScalar("-poisson_ratio", "poisson ratio", "", nu, &nu,
                              PETSC_NULLPTR);
 
    PetscOptionsEnd();
 
    MOFEM_LOG("EP", Sev::inform) << "St Venant Kirchhoff model parameters: "
                                 << "E = " << E << ", nu = " << nu;
 
    lambda = LAMBDA(E, nu);
    mu = MU(E, nu);
 
    MoFEMFunctionReturn(0);
  }

recordTape()

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

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 71 of file HMHHStVenantKirchhoff.cpp.

                                                                  {
    MoFEMFunctionBegin;
    FTensor::Index<'i', 3> i;
    FTensor::Index<'j', 3> j;
    FTensor::Index<'I', 3> I;
    FTensor::Index<'J', 3> J;
    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);
 
    // Deformation and strain
    tC(I, J) = tF(i, I) * tF(i, J);
    tE(I, J) = tC(I, J);
    tE(N0, N0) -= 1;
    tE(N1, N1) -= 1;
    tE(N2, N2) -= 1;
    tE(I, J) *= 0.5;
 
    // Energy
    trE = tE(I, I);
    energy = 0.5 * lambda * trE * trE + mu * (tE(I, J) * tE(I, J));
 
    // Stress Piola II
    tS(I, J) = tE(I, J);
    tS(I, J) *= 2 * mu;
    tS(N0, N0) += lambda * trE;
    tS(N1, N1) += lambda * trE;
    tS(N2, N2) += lambda * trE;
    // Stress Piola I
    tP(i, J) = tF(i, I) * tS(I, J);
 
    // Set dependent variables to ADOL-C
    tP(i, j) >>= r_P(i, j);
 
    trace_off();
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

detF

adouble EshelbianPlasticity::HMHStVenantKirchhoff::detF

Definition at line 60 of file HMHHStVenantKirchhoff.cpp.

detH

adouble EshelbianPlasticity::HMHStVenantKirchhoff::detH

Definition at line 59 of file HMHHStVenantKirchhoff.cpp.

energy

adouble EshelbianPlasticity::HMHStVenantKirchhoff::energy

Definition at line 57 of file HMHHStVenantKirchhoff.cpp.

f

adouble EshelbianPlasticity::HMHStVenantKirchhoff::f

Definition at line 56 of file HMHHStVenantKirchhoff.cpp.

lambda

double EshelbianPlasticity::HMHStVenantKirchhoff::lambda

Definition at line 46 of file HMHHStVenantKirchhoff.cpp.

mu

double EshelbianPlasticity::HMHStVenantKirchhoff::mu

Definition at line 47 of file HMHHStVenantKirchhoff.cpp.

numberOfActiveVariables

int EshelbianPlasticity::HMHStVenantKirchhoff::numberOfActiveVariables = 9

staticconstexpr

Definition at line 15 of file HMHHStVenantKirchhoff.cpp.

numberOfDependentVariables

int EshelbianPlasticity::HMHStVenantKirchhoff::numberOfDependentVariables = 9

staticconstexpr

Definition at line 16 of file HMHHStVenantKirchhoff.cpp.

phi

adouble EshelbianPlasticity::HMHStVenantKirchhoff::phi

Definition at line 54 of file HMHHStVenantKirchhoff.cpp.

s2

adouble EshelbianPlasticity::HMHStVenantKirchhoff::s2

Definition at line 55 of file HMHHStVenantKirchhoff.cpp.

sigmaY

double EshelbianPlasticity::HMHStVenantKirchhoff::sigmaY

Definition at line 48 of file HMHHStVenantKirchhoff.cpp.

tC

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tC

Definition at line 65 of file HMHHStVenantKirchhoff.cpp.

tE

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tE

Definition at line 66 of file HMHHStVenantKirchhoff.cpp.

tF

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tF

Definition at line 63 of file HMHHStVenantKirchhoff.cpp.

tH

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tH

Definition at line 51 of file HMHHStVenantKirchhoff.cpp.

th

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::th

Definition at line 50 of file HMHHStVenantKirchhoff.cpp.

tInvF

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tInvF

Definition at line 64 of file HMHHStVenantKirchhoff.cpp.

tInvH

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tInvH

Definition at line 62 of file HMHHStVenantKirchhoff.cpp.

tP

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tP

Definition at line 52 of file HMHHStVenantKirchhoff.cpp.

tPulledP

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tPulledP

Definition at line 69 of file HMHHStVenantKirchhoff.cpp.

trE

adouble EshelbianPlasticity::HMHStVenantKirchhoff::trE

Definition at line 61 of file HMHHStVenantKirchhoff.cpp.

tS

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tS

Definition at line 67 of file HMHHStVenantKirchhoff.cpp.

tSigma

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tSigma

Definition at line 53 of file HMHHStVenantKirchhoff.cpp.

---

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

• users_modules/eshelbian_plasticity/src/impl/HMHHStVenantKirchhoff.cpp