Inheritance diagram for EshelbianPlasticity::HMHStVenantKirchhoff:

Collaboration diagram for EshelbianPlasticity::HMHStVenantKirchhoff:

Public Member Functions
	HMHStVenantKirchhoff (const double lambda, const double mu)

MoFEMErrorCode	getOptions ()

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	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	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.

19 : PhysicalEquations(numberOfActiveVariables, numberOfDependentVariables),

20 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;
  
     CHKERR PetscOptionsBegin(PETSC_COMM_WORLD, "stvenant_", "", "none");
  
     CHKERR PetscOptionsScalar("-young_modulus", "Young modulus", "", E, &E,
                               PETSC_NULL);
     CHKERR PetscOptionsScalar("-poisson_ratio", "poisson ratio", "", nu, &nu,
                               PETSC_NULL);
  
     ierr = PetscOptionsEnd();
     CHKERRG(ierr);
  
     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 79 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);
   }

◆ returnOpJacobian()

virtual OpJacobian* EshelbianPlasticity::HMHStVenantKirchhoff::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 48 of file HMHHStVenantKirchhoff.cpp.

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

Member Data Documentation

◆ detF

adouble EshelbianPlasticity::HMHStVenantKirchhoff::detF

Definition at line 68 of file HMHHStVenantKirchhoff.cpp.

◆ detH

adouble EshelbianPlasticity::HMHStVenantKirchhoff::detH

Definition at line 67 of file HMHHStVenantKirchhoff.cpp.

◆ energy

adouble EshelbianPlasticity::HMHStVenantKirchhoff::energy

Definition at line 65 of file HMHHStVenantKirchhoff.cpp.

◆ f

adouble EshelbianPlasticity::HMHStVenantKirchhoff::f

Definition at line 64 of file HMHHStVenantKirchhoff.cpp.

◆ lambda

double EshelbianPlasticity::HMHStVenantKirchhoff::lambda

Definition at line 54 of file HMHHStVenantKirchhoff.cpp.

◆ mu

double EshelbianPlasticity::HMHStVenantKirchhoff::mu

Definition at line 55 of file HMHHStVenantKirchhoff.cpp.

◆ numberOfActiveVariables

constexpr int EshelbianPlasticity::HMHStVenantKirchhoff::numberOfActiveVariables = 9

staticconstexpr

Definition at line 15 of file HMHHStVenantKirchhoff.cpp.

◆ numberOfDependentVariables

constexpr int EshelbianPlasticity::HMHStVenantKirchhoff::numberOfDependentVariables = 9

staticconstexpr

Definition at line 16 of file HMHHStVenantKirchhoff.cpp.

◆ phi

adouble EshelbianPlasticity::HMHStVenantKirchhoff::phi

Definition at line 62 of file HMHHStVenantKirchhoff.cpp.

◆ s2

adouble EshelbianPlasticity::HMHStVenantKirchhoff::s2

Definition at line 63 of file HMHHStVenantKirchhoff.cpp.

◆ sigmaY

double EshelbianPlasticity::HMHStVenantKirchhoff::sigmaY

Definition at line 56 of file HMHHStVenantKirchhoff.cpp.

◆ tC

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tC

Definition at line 73 of file HMHHStVenantKirchhoff.cpp.

◆ tE

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tE

Definition at line 74 of file HMHHStVenantKirchhoff.cpp.

◆ tF

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tF

Definition at line 71 of file HMHHStVenantKirchhoff.cpp.

◆ th

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::th

Definition at line 58 of file HMHHStVenantKirchhoff.cpp.

◆ tH

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tH

Definition at line 59 of file HMHHStVenantKirchhoff.cpp.

◆ tInvF

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tInvF

Definition at line 72 of file HMHHStVenantKirchhoff.cpp.

◆ tInvH

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tInvH

Definition at line 70 of file HMHHStVenantKirchhoff.cpp.

◆ tP

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tP

Definition at line 60 of file HMHHStVenantKirchhoff.cpp.

◆ tPulledP

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tPulledP

Definition at line 77 of file HMHHStVenantKirchhoff.cpp.

◆ trE

adouble EshelbianPlasticity::HMHStVenantKirchhoff::trE

Definition at line 69 of file HMHHStVenantKirchhoff.cpp.

◆ tS

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tS

Definition at line 75 of file HMHHStVenantKirchhoff.cpp.

◆ tSigma

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tSigma

Definition at line 61 of file HMHHStVenantKirchhoff.cpp.

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

users_modules/eshelbian_plasticit/src/impl/HMHHStVenantKirchhoff.cpp

Public Member Functions

Public Attributes

Static Public Attributes

Additional Inherited Members