Inheritance diagram for MatOps::MatMetaElastic< DIM >:

Collaboration diagram for MatOps::MatMetaElastic< DIM >:

Public Types
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	getOptions (MoFEM::Interface *m_field_ptr=nullptr) override

MoFEMErrorCode	setParams (FEMethod *fe_ptr, int gg) override

MoFEMErrorCode	recordTape () override

MoFEMErrorCode	evaluateVariable (int tag, EntityHandle entity, int gg) override

MoFEMErrorCode	evaluateDerivatives (int tag, EntityHandle entity, int gg) override

MoFEMErrorCode	updateState (int tag, EntityHandle entity, int gg) override

std::map< int, boost::shared_ptr< PhysicalEquations > > &	getMetaElasticMap ()

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	addBlockParameters (MoFEM::Interface &m_field, const std::string &block_name, int block_id, const Range &block_entities, const std::vector< double > &block_data)

Protected Attributes
std::map< int, boost::shared_ptr< PhysicalEquations > >	metaElasticAtPts

Static Protected Attributes
static constexpr int	modelType

Additional Inherited Members
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

Detailed Description

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

Definition at line 13 of file MatMetaElastic.cpp.

Member Typedef Documentation

◆ A

template<int DIM>

using MatOps::MatMetaElastic< DIM >::A = MatElasticImpl<DIM>

Definition at line 16 of file MatMetaElastic.cpp.

Member Function Documentation

◆ evaluateDerivatives()

template<int DIM>

MoFEMErrorCode MatOps::MatMetaElastic< DIM >::evaluateDerivatives	(	int	tag,
		EntityHandle	entity,
		int	gg
	)

inlineoverridevirtual

Reimplemented from MatOps::PhysicalEquations.

Definition at line 94 of file MatMetaElastic.cpp.

                                                      {
    return metaElasticAtPts.at(tag)->evaluateDerivatives(tag, entity, gg);
  }

◆ evaluateVariable()

template<int DIM>

MoFEMErrorCode MatOps::MatMetaElastic< DIM >::evaluateVariable	(	int	tag,
		EntityHandle	entity,
		int	gg
	)

inlineoverridevirtual

Reimplemented from MatOps::PhysicalEquations.

Definition at line 89 of file MatMetaElastic.cpp.

                                                   {
    return metaElasticAtPts.at(tag)->evaluateVariable(tag, entity, gg);
  }

◆ getMetaElasticMap()

template<int DIM>

std::map< int, boost::shared_ptr< PhysicalEquations > > & MatOps::MatMetaElastic< DIM >::getMetaElasticMap ( )

inline

Definition at line 103 of file MatMetaElastic.cpp.

                                                                       {
    return metaElasticAtPts;
  }

◆ getOptions()

template<int DIM>

MoFEMErrorCode MatOps::MatMetaElastic< DIM >::getOptions ( MoFEM::Interface * m_field_ptr = nullptr )

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 18 of file MatMetaElastic.cpp.

                                                                            {
    MoFEMFunctionBegin;
 
    char umat_library[PETSC_MAX_PATH_LEN] = "./umat.so";
    CHKERR PetscOptionsGetString(PETSC_NULLPTR, PETSC_NULLPTR,
                                 "-umat_library", umat_library,
                                 sizeof(umat_library), PETSC_NULLPTR);
 
    if (std::filesystem::exists(std::filesystem::path(umat_library))) {
      metaElasticAtPts.insert(
          {MatOps::MatOpsTagsRegistry::setTagName("Umat"),
           createMatOpsPhysicalEquationsPtr<UMAT, modelType>(
               this->matOpsDataPtr,
               MatOps::MatOpsTagsRegistry::setTagName("Umat"))});
    }
 
    for (auto &m : metaElasticAtPts) {
      CHKERR m.second->getOptions(m_field_ptr);
    }
    A::tagVsRangePtr = boost::make_shared<std::map<int, Range>>();
    for (const auto &m : metaElasticAtPts) {
      A::tagVsRangePtr->insert(
          {m.first, Range()}); // Initialize with empty range
      for (auto &rp : m.second->paramVecByRange) {
        A::tagVsRangePtr->at(m.first).merge(
            rp.first); // Merge ranges for the same tag
      }
    }
    MoFEMFunctionReturn(0);
  }

◆ recordTape()

template<int DIM>

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

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 81 of file MatMetaElastic.cpp.

                                       {
    MoFEMFunctionBegin;
    for (auto &m : metaElasticAtPts) {
      CHKERR m.second->recordTape();
    }
    MoFEMFunctionReturn(0);
  }

◆ setParams()

template<int DIM>

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

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 49 of file MatMetaElastic.cpp.

                                                              {
    MoFEMFunctionBegin;
    const auto ent = fe_ptr->getFEEntityHandle();
 
    bool has_ranges = false;
    if (A::tagVsRangePtr) {
      for (const auto &tag_range : *A::tagVsRangePtr) {
        has_ranges = has_ranges || !tag_range.second.empty();
      }
    }
 
    if (!has_ranges) {
      for (auto &m : metaElasticAtPts) {
        CHKERR m.second->setParams(fe_ptr, gg);
      }
      MoFEMFunctionReturnHot(0);
    }
 
    for (auto &m : metaElasticAtPts) {
      if (A::tagVsRangePtr->at(m.first).find(ent) !=
          A::tagVsRangePtr->at(m.first).end()) {
        CHKERR m.second->setParams(fe_ptr, gg);
        MoFEMFunctionReturnHot(0);
      }
    }
 
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Meta elastic material parameters not found for entity %lu",
            static_cast<unsigned long>(ent));
    MoFEMFunctionReturn(0);
  }

◆ updateState()

template<int DIM>

MoFEMErrorCode MatOps::MatMetaElastic< DIM >::updateState	(	int	tag,
		EntityHandle	entity,
		int	gg
	)

inlineoverridevirtual

Reimplemented from MatOps::PhysicalEquations.

Definition at line 99 of file MatMetaElastic.cpp.

                                                                            {
    return metaElasticAtPts.at(tag)->updateState(tag, entity, gg);
  }

Member Data Documentation

◆ metaElasticAtPts

template<int DIM>

std::map<int, boost::shared_ptr<PhysicalEquations> > MatOps::MatMetaElastic< DIM >::metaElasticAtPts

protected

Definition at line 111 of file MatMetaElastic.cpp.

                                                                     {
 
      {
 
          // First tag is always MetaElastic tag
          MatOps::MatOpsTagsRegistry::setTagName("MetaElastic"),
          createMatOpsPhysicalEquationsPtr<ELASTICITY::NEOHOOKEAN, modelType>(
              this->matOpsDataPtr,
              MatOps::MatOpsTagsRegistry::setTagName("MetaElastic"))
 
      },
 
      {
 
          MatOps::MatOpsTagsRegistry::setTagName("MooneyRivlinWriggersEq63"),
          createMatOpsPhysicalEquationsPtr<ELASTICITY::MOONEYRIVLINWRIGGERSEQ63,
                                          modelType>(
              this->matOpsDataPtr, MatOps::MatOpsTagsRegistry::setTagName(
                                      "MooneyRivlinWriggersEq63"))
 
      },
 
      {
 
          MatOps::MatOpsTagsRegistry::setTagName("StVenantKirchhoff"),
          createMatOpsPhysicalEquationsPtr<ELASTICITY::STVENANTKIRCHHOFF,
                                          modelType>(
              this->matOpsDataPtr,
              MatOps::MatOpsTagsRegistry::setTagName("StVenantKirchhoff"))
 
      },
 
  };