MatOps::MatNeohookean< DIM > Struct Template Reference

Inheritance diagram for MatOps::MatNeohookean< DIM >:

Collaboration diagram for MatOps::MatNeohookean< DIM >:

Public Types
using	A = MatNeohookeanGeneric< MatElasticImpl< DIM > >
Public Types inherited from MatOps::MatNeohookeanGeneric< MatElasticImpl< DIM > >
using	A = MatElasticImpl< DIM >
Public Types inherited from MatOps::PhysicalEquations
using	HookFunction = std::function< MoFEMErrorCode(boost::shared_ptr< MatOpsData >, int, EntityHandle, int)>

Public Member Functions
MoFEMErrorCode	setParams (FEMethod *fe_ptr, int gg) override
MoFEMErrorCode	recordTape () override
Public Member Functions inherited from MatOps::MatNeohookeanGeneric< MatElasticImpl< DIM > >
MoFEMErrorCode	getOptions (MoFEM::Interface *m_field_ptr=nullptr) override
MoFEMErrorCode	addBlockParameters (MoFEM::Interface &m_field, const std::string &block_name, int block_id, const Range &block_entities, const std::vector< double > &block_data) override
Public Member Functions inherited from MatOps::MatElasticImpl< DIM >
ForcesAndSourcesCore::UserDataOperator *	createOp (boost::shared_ptr< PhysicalEquations > physical_ptr, bool eval_stress, bool eval_tangent, bool update) override
ForcesAndSourcesCore::UserDataOperator *	createOp (boost::shared_ptr< PhysicalEquations > physical_ptr, bool eval_stress, bool eval_tangent, bool update)
ForcesAndSourcesCore::UserDataOperator *	createOp (boost::shared_ptr< PhysicalEquations > physical_ptr, bool eval_stress, bool eval_tangent, bool update)
Public Member Functions inherited from MatOps::MatElastic
	PhysicalEquations ()=delete
	PhysicalEquations (boost::shared_ptr< MatOpsData > mat_ops_data_ptr, int tag, boost::shared_ptr< std::map< int, Range > > tag_vs_range_ptr=nullptr)
Public Member Functions inherited from MatOps::PhysicalEquations
	PhysicalEquations ()=delete
	PhysicalEquations (boost::shared_ptr< MatOpsData > mat_ops_data_ptr, int tag, boost::shared_ptr< std::map< int, Range > > tag_vs_range_ptr=nullptr)
virtual	~PhysicalEquations ()=default
virtual MoFEMErrorCode	evaluateVariable (int tag, EntityHandle entity, int gg)
virtual MoFEMErrorCode	evaluateDerivatives (int tag, EntityHandle entity, int gg)
virtual MoFEMErrorCode	updateState (int tag, EntityHandle entity, int gg)

Additional Inherited Members
Static Public Member Functions inherited from MatOps::MatNeohookeanGeneric< MatElasticImpl< DIM > >
static double	getShearModulus (const double c10)
static double	getAbaqusD1 (const double K)
static MoFEMErrorCode	validateMaterialParameters (const double c10, const double K, const char *source)
Public Attributes inherited from MatOps::PhysicalEquations
HookFunction	hookEvaluateVariable = AdolCEvaluation::evaluateVariable
HookFunction	hookEvaluateDerivatives = AdolCEvaluation::evaluateDerivatives
HookFunction	hookUpdateState = AdolCEvaluation::evaluateVariable
int	tAg
boost::shared_ptr< std::map< int, Range > >	tagVsRangePtr
std::vector< std::pair< Range, std::vector< double > > >	paramVecByRange
boost::shared_ptr< MatOpsData >	matOpsDataPtr
Static Public Attributes inherited from MatOps::MatElastic
static bool	useDeformationGradient
Protected Attributes inherited from MatOps::MatNeohookeanGeneric< MatElasticImpl< DIM > >
double	C10
double	K
std::vector< double >	defaultMaterialParameters
std::vector< double >	commandLineParameters

Detailed Description

template<int DIM>
struct MatOps::MatNeohookean< DIM >

Definition at line 167 of file MatNeohookean.cpp.

Member Typedef Documentation

A

template<int DIM>

using MatOps::MatNeohookean< DIM >::A = MatNeohookeanGeneric<MatElasticImpl<DIM> >

Definition at line 168 of file MatNeohookean.cpp.

Member Function Documentation

recordTape()

template<int DIM>

MoFEMErrorCode MatOps::MatNeohookean< DIM >::recordTape ( )

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 191 of file MatNeohookean.cpp.

                                       {
    MoFEMFunctionBegin;
 
    A::matOpsDataPtr->insertCommonData("grad", MatrixDouble());
    A::matOpsDataPtr->insertCommonData("P", MatrixDouble());
    A::matOpsDataPtr->insertCommonData("P_dF", MatrixDouble());
 
    A::matOpsDataPtr->insertActiveData("F", MatrixDouble());
    A::matOpsDataPtr->insertDependentData("P", MatrixDouble());
    A::matOpsDataPtr->insertDependentDerivativesData("P_dF", MatrixDouble());
 
    A::matOpsDataPtr->getActiveDataPtr("F")->resize(DIM, DIM, false);
    A::matOpsDataPtr->getDependentDataPtr("P")->resize(DIM, DIM, false);
 
    auto t_F = getFTensor2FromPtr<DIM, DIM>(
        A::matOpsDataPtr->getActiveDataPtr("F")->data().data());
    auto t_P = getFTensor2FromPtr<DIM, DIM>(
        A::matOpsDataPtr->getDependentDataPtr("P")->data().data());
 
    auto t_kd = FTensor::Kronecker_Delta<int>();
 
    FTENSOR_INDEX(DIM, i);
    FTENSOR_INDEX(DIM, j);
    FTENSOR_INDEX(DIM, I);
    FTENSOR_INDEX(DIM, J);
 
    t_F(i, J) = 0;
 
    FTensor::Tensor2<adouble, DIM, DIM> ta_F;
    FTensor::Tensor2<adouble, DIM, DIM> ta_P;
 
    adouble det_aF, I1;
    FTensor::Tensor2<adouble, DIM, DIM> ta_invF;
 
    trace_on(A::tAg);
    auto p_c10 = mkparam(this->C10);
    auto p_K = mkparam(this->K);
 
    ta_F(i, J) <<= t_F(i, J);
    // assume that gradient from approximated displacement, that why we add diagonal
    if (!MatElastic::useDeformationGradient) {
      ta_F(i, J) += t_kd(i, J);
    }
 
    det_aF = determinantTensor(ta_F);
    I1 = ta_F(i, I) * ta_F(i, I);
    if constexpr (DIM == 2) {
      // Abaqus plane strain uses the three-dimensional potential with F33=1.
      I1 += 1.;
    }
    CHKERR invertTensor(ta_F, det_aF, ta_invF);
 
    // Abaqus compressible Neo-Hookean potential (reduced polynomial N=1):
    //
    //   W = C10 * (I1_bar - 3) + (1 / D1) * (J - 1)^2,
    //   I1_bar = J^(-2/3) * I1,  mu0 = 2*C10,  D1 = 2/K0.
    //
    // Its first Piola stress is
    //
    //   P = 2*C10*J^(-2/3)*(F - I1/3*F^(-T))
    //       + K0*J*(J - 1)*F^(-T).
    const auto jacobian_to_minus_two_thirds = pow(det_aF, -2. / 3.);
    ta_P(i, I) =
        2. * p_c10 * jacobian_to_minus_two_thirds *
            (ta_F(i, I) - (I1 / 3.) * ta_invF(I, i)) +
        p_K * det_aF * (det_aF - 1.) * ta_invF(I, i);
 
    ta_P(i, I) >>= t_P(i, I);
 
    trace_off();
 
    MoFEMFunctionReturn(0);
  }

setParams()

template<int DIM>

MoFEMErrorCode MatOps::MatNeohookean< DIM >::setParams ( FEMethod *  fe_ptr, int  gg  )

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 171 of file MatNeohookean.cpp.

                                                              {
    (void)gg;
    const auto ent = fe_ptr->getFEEntityHandle();
    for (auto &[range, param_vec] : A::paramVecByRange) {
      if (std::find(range.begin(), range.end(), ent) != range.end()) {
        set_param_vec(A::tAg, param_vec.size(), param_vec.data());
        return 0;
      }
    }
    // command line overrides block parameters
    if (!this->commandLineParameters.empty()) {
      set_param_vec(A::tAg, this->commandLineParameters.size(),
                    this->commandLineParameters.data());
      return 0;
    }
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Neohookean material parameters not found for entity %lu",
            static_cast<unsigned long>(ent));
  }

---

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

• src/materials/impl/MatNeohookean.cpp