MatGenericElastic.cpp

Go to the documentation of this file.

#include <MoFEM.hpp>
 
using namespace MoFEM;
 
#include "MatOps.hpp"
#include "MatGenericElastic.hpp"
 
namespace MatOps {
 
GenericElastic::GenericElastic(boost::shared_ptr<MatOpsData> mat_ops_data_ptr,
                               int tag)
    : PhysicalEquations(mat_ops_data_ptr, tag) {}
 
MoFEMErrorCode GenericElastic::getOptions(MoFEM::Interface *m_field_ptr) {
  MoFEMFunctionBegin;
  (void)m_field_ptr;
  MoFEMFunctionReturn(0);
}
 
MoFEMErrorCode GenericElastic::setParams(FEMethod *fe_ptr, int gg) {
  MoFEMFunctionBegin;
  (void)fe_ptr;
  (void)gg;
  MoFEMFunctionReturn(0);
}
 
MoFEMErrorCode GenericElastic::recordTape() {
  MoFEMFunctionBegin;
  matOpsDataPtr->insertCommonData("grad", MatrixDouble());
  matOpsDataPtr->insertCommonData("P", MatrixDouble());
  matOpsDataPtr->insertCommonData("P_dF", MatrixDouble());
  matOpsDataPtr->insertActiveData("F", MatrixDouble());
  matOpsDataPtr->insertDependentData("P", MatrixDouble());
  matOpsDataPtr->insertDependentDerivativesData("P_dF", MatrixDouble());
  MoFEMFunctionReturn(0);
}
 
template <int DIM> struct GenericElasticImpl : public GenericElastic {
  using GenericElastic::GenericElastic;
  ForcesAndSourcesCore::UserDataOperator *
  createOp(boost::shared_ptr<PhysicalEquations> physical_ptr, bool eval_stress,
           bool eval_tangent, bool update) {
    using EleOp = ForcesAndSourcesCore::UserDataOperator;
    return new OpEvalMatSimpleMaterialImpl<DIM, EleOp>(
        physical_ptr, eval_stress, eval_tangent, update);
  }
};
 
template <>
boost::shared_ptr<PhysicalEquations>
createMatOpsPhysicalEquationsPtr<GenericElastic, MODEL_3D>(
    boost::shared_ptr<MatOpsData> mat_ops_data_ptr, int tag) {
  return boost::make_shared<GenericElasticImpl<3>>(mat_ops_data_ptr, tag);
}
 
template <>
boost::shared_ptr<PhysicalEquations>
createMatOpsPhysicalEquationsPtr<GenericElastic, MODEL_2D_PLANE_STRAIN>(
    boost::shared_ptr<MatOpsData> mat_ops_data_ptr, int tag) {
  return boost::make_shared<GenericElasticImpl<2>>(mat_ops_data_ptr, tag);
}
 
} // namespace MatOps