Inheritance diagram for MatOps::MatMetaElastic< DIM, MODEL_TYPE >:

Collaboration diagram for MatOps::MatMetaElastic< DIM, MODEL_TYPE >:

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)

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

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

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 = MODEL_TYPE

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, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>
struct MatOps::MatMetaElastic< DIM, MODEL_TYPE >

Definition at line 16 of file MatMetaElastic.cpp.

Member Typedef Documentation

◆ A

template<int DIM, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

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

Definition at line 19 of file MatMetaElastic.cpp.

Member Function Documentation

◆ evaluateDerivatives()

template<int DIM, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

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

inlinevirtual

Reimplemented from MatOps::PhysicalEquations.

Definition at line 101 of file MatMetaElastic.cpp.

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

◆ evaluateVariable()

template<int DIM, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

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

inlinevirtual

Reimplemented from MatOps::PhysicalEquations.

Definition at line 97 of file MatMetaElastic.cpp.

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

◆ getMetaElasticMap()

template<int DIM, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

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

inline

Definition at line 109 of file MatMetaElastic.cpp.

                                                                       {
    return metaElasticAtPts;
  }

◆ getOptions()

template<int DIM, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

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

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 21 of file MatMetaElastic.cpp.

                                                                            {
    MoFEMFunctionBegin;
 
    // Activate UMAT only for an explicit library or material block.
    char umat_library[PETSC_MAX_PATH_LEN] = "";
    CHKERR PetscOptionsGetString(PETSC_NULLPTR, PETSC_NULLPTR,
                                 "-umat_library", umat_library,
                                 sizeof(umat_library), PETSC_NULLPTR);
 
    if (umat_library[0] ||
        (m_field_ptr &&
         !m_field_ptr->getInterface<MeshsetsManager>()
              ->getCubitMeshsetPtr(std::regex("MAT_UMAT(.*)"))
              .empty())) {
      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, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

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

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 89 of file MatMetaElastic.cpp.

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

◆ setParams()

template<int DIM, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

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

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 57 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, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

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

inlinevirtual

Reimplemented from MatOps::PhysicalEquations.

Definition at line 105 of file MatMetaElastic.cpp.

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

Member Data Documentation

◆ metaElasticAtPts

template<int DIM, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

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

protected

Definition at line 116 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"))
 
      },
 
  };

◆ modelType

template<int DIM, int MODEL_TYPE = (DIM == 3) ? MatOps::MODEL_3D : MatOps::MODEL_2D_PLANE_STRAIN>

constexpr int MatOps::MatMetaElastic< DIM, MODEL_TYPE >::modelType = MODEL_TYPE

inlinestaticconstexprprotected

Definition at line 114 of file MatMetaElastic.cpp.

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

src/materials/impl/MatMetaElastic.cpp

Public Types

Public Member Functions

Protected Attributes

Static Protected Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ A

Member Function Documentation

◆ evaluateDerivatives()

◆ evaluateVariable()

◆ getMetaElasticMap()

◆ getOptions()

◆ recordTape()

◆ setParams()

◆ updateState()

Member Data Documentation

◆ metaElasticAtPts

◆ modelType