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
MatOpsData::MatrixPtr	getEnergyDensityPtr () const
Public Member Functions inherited from MatOps::MatPiolaResponse
	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 constexpr auto	energyDensityDataName = "energy"
Static Public Attributes inherited from MatOps::MatPiolaResponse
static bool	useDeformationGradient
	Otherwise the displacement gradient is supplied and I is added.
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 159 of file MatNeohookean.cpp.

Member Typedef Documentation

A

template<int DIM>

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

Definition at line 160 of file MatNeohookean.cpp.

Member Function Documentation

recordTape()

template<int DIM>

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

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 183 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->insertDependentData("energy", MatrixDouble());
    A::matOpsDataPtr->insertDependentDerivativesData("P_dF", MatrixDouble());
    A::matOpsDataPtr->insertDependentDerivativesData("energy_dF",
                                                     MatrixDouble());
 
    A::matOpsDataPtr->getActiveDataPtr("F")->resize(DIM, DIM, false);
    A::matOpsDataPtr->getDependentDataPtr("P")->resize(DIM, DIM, false);
    A::matOpsDataPtr->getDependentDataPtr("energy")->resize(1, 1, false);
    A::matOpsDataPtr->getDependentDerivativesDataPtr("energy_dF")
        ->resize(1, 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 energy_ptr = A::matOpsDataPtr->getDependentDataPtr("energy");
 
    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, ta_energy;
    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_energy = p_c10 * (jacobian_to_minus_two_thirds * I1 - 3.) +
                0.5 * p_K * (det_aF - 1.) * (det_aF - 1.);
 
    ta_P(i, I) >>= t_P(i, I);
    ta_energy >>= (*energy_ptr)(0, 0);
 
    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 163 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