MatOps::MatHuHu< DIM > Struct Template Reference

Inheritance diagram for MatOps::MatHuHu< DIM >:

Collaboration diagram for MatOps::MatHuHu< DIM >:

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

Public Member Functions
ForcesAndSourcesCore::UserDataOperator *	createOp (boost::shared_ptr< PhysicalEquations > physical_ptr, bool eval_stress, bool eval_tangent, bool update) override
MoFEMErrorCode	getOptions (MoFEM::Interface *m_field_ptr=nullptr) override
MoFEMErrorCode	setParams (FEMethod *fe_ptr, int gg) override
MoFEMErrorCode	recordTape () 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)
	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	setParams (FEMethod *fe_ptr, int gg)=0
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)

Protected Attributes
double	K = 1e-6
std::vector< double >	defaultMaterialParameters = {K}

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

Detailed Description

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

Definition at line 21 of file MatHuHu.cpp.

Member Typedef Documentation

A

template<int DIM>

using MatOps::MatHuHu< DIM >::A = PhysicalEquations

Definition at line 23 of file MatHuHu.cpp.

Member Function Documentation

createOp() [1/3]

ForcesAndSourcesCore::UserDataOperator * MatOps::MatHuHu< 3 >::createOp ( boost::shared_ptr< PhysicalEquations >  physical_ptr, bool  eval_stress, bool  eval_tangent, bool  update  )

virtual

Implements MatOps::PhysicalEquations.

Definition at line 291 of file MatHuHu.cpp.

                                                                         {
  using EleOp = ForcesAndSourcesCore::UserDataOperator;
  return createOpImpl<3, EleOp>(physical_ptr, eval_stress, eval_tangent,
                                update);
}

createOp() [2/3]

ForcesAndSourcesCore::UserDataOperator * MatOps::MatHuHu< 2 >::createOp ( boost::shared_ptr< PhysicalEquations >  physical_ptr, bool  eval_stress, bool  eval_tangent, bool  update  )

virtual

Implements MatOps::PhysicalEquations.

Definition at line 300 of file MatHuHu.cpp.

                                                                         {
  using EleOp = ForcesAndSourcesCore::UserDataOperator;
  return createOpImpl<2, EleOp>(physical_ptr, eval_stress, eval_tangent,
                                update);
}

createOp() [3/3]

template<int DIM>

ForcesAndSourcesCore::UserDataOperator * MatOps::MatHuHu< DIM >::createOp ( boost::shared_ptr< PhysicalEquations >  physical_ptr, bool  eval_stress, bool  eval_tangent, bool  update  )

overridevirtual

Implements MatOps::PhysicalEquations.

getOptions()

template<int DIM>

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

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 28 of file MatHuHu.cpp.

                                                                            {
    MoFEMFunctionBegin;
    MOFEM_LOG_CHANNEL("WORLD");
 
    PetscOptionsBegin(PETSC_COMM_WORLD, "huhu_", "", "none");
    CHKERR PetscOptionsScalar("-K", "Bulk modulus K", "", K, &K, PETSC_NULLPTR);
    PetscOptionsEnd();
 
    MOFEM_TAG_AND_LOG("WORLD", Sev::inform, "Default HuHu parameters:")
        << " K = " << K;
    defaultMaterialParameters = {K};
 
    std::string block_name = "MAT_HUHU";
 
    for (auto &m :
 
         m_field_ptr->getInterface<MeshsetsManager>()->
 
         getCubitMeshsetPtr(
             std::regex((boost::format("%s(.*)") % block_name).str()))
 
    ) {
 
      std::vector<double> block_data;
      CHKERR m->getAttributes(block_data);
      if (block_data.size() < 1) {
        block_data.push_back(
            K); // Use default K if not provided in block attributes
      }
      auto get_block_ents = [&]() {
        Range ents;
        CHK_MOAB_THROW(m_field_ptr->get_moab().get_entities_by_handle(
                           m->meshset, ents, true),
                       "can not get block entities");
        return ents;
      };
 
      A::paramVecByRange.push_back({get_block_ents(), {block_data[0]}});
 
      MOFEM_TAG_AND_LOG("WORLD", Sev::inform, "MatBlock for HuHu")
          << *m << " K = " << block_data[0];
    }
 
    MoFEMFunctionReturn(0);
  };

PhysicalEquations() [1/2]

template<int DIM>

MatOps::PhysicalEquations::PhysicalEquations ( )

delete

PhysicalEquations() [2/2]

template<int DIM>

MatOps::PhysicalEquations::PhysicalEquations ( boost::shared_ptr< MatOpsData >  mat_ops_data_ptr, int  tag, boost::shared_ptr< std::map< int, Range > >  tag_vs_range_ptr = nullptr  )

inline

Definition at line 114 of file MatOps.hpp.

      : matOpsDataPtr(mat_ops_data_ptr), tAg(tag),
        tagVsRangePtr(tag_vs_range_ptr) {}

recordTape()

template<int DIM>

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

inlineoverridevirtual

Implements MatOps::PhysicalEquations.

Definition at line 88 of file MatHuHu.cpp.

                                       {
    MoFEMFunctionBegin;
 
    A::matOpsDataPtr->insertCommonData("grad", MatrixDouble(DIM * DIM, 1));
    A::matOpsDataPtr->insertCommonData("k", MatrixDouble(1, 1));
    A::matOpsDataPtr->insertCommonData("k_dF", MatrixDouble(DIM * DIM, 1));
 
    A::matOpsDataPtr->insertActiveData("F", MatrixDouble(DIM, DIM));
    A::matOpsDataPtr->insertDependentData("k", MatrixDouble(1, 1));
    A::matOpsDataPtr->insertDependentDerivativesData("k_dF",
                                                    MatrixDouble(DIM, DIM));
 
    auto t_F = getFTensor2FromPtr<DIM, DIM>(
        A::matOpsDataPtr->getActiveDataPtr("F")->data().data());
    auto k = (*A::matOpsDataPtr->getDependentDataPtr("k"))(0, 0);
 
    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) = t_kd(i, J);
 
    FTensor::Tensor2<adouble, DIM, DIM> ta_F;
    adouble ta_k;
 
    adouble det_aF;
    FTensor::Tensor2<adouble, DIM, DIM> ta_invF;
 
    trace_on(A::tAg);
    auto p_K = mkparam(K);
    ta_F(i, J) <<= t_F(i, J);
    ta_k = p_K;
    ta_k >>= k;
 
    trace_off();
 
    MoFEMFunctionReturn(0);
  }

setParams()

template<int DIM>

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

inlineoverride

Definition at line 74 of file MatHuHu.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;
      }
    }
    set_param_vec(A::tAg, defaultMaterialParameters.size(),
                  defaultMaterialParameters.data());
    return 0;
  }

Member Data Documentation

defaultMaterialParameters

template<int DIM>

std::vector<double> MatOps::MatHuHu< DIM >::defaultMaterialParameters = {K}

protected

Definition at line 132 of file MatHuHu.cpp.

{K};

K

template<int DIM>

double MatOps::MatHuHu< DIM >::K = 1e-6

protected

Definition at line 131 of file MatHuHu.cpp.

---

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

• src/materials/impl/MatHuHu.cpp