MatUmat.hpp

Go to the documentation of this file.

/**
 * @file MatUmat.hpp
 * @brief 
 * @version 0.1
 * @date 2026-04-02
 * 
 * @copyright Copyright (c) 2026
 * 
 */
 
#ifndef MAT_OPS_UMAT_HPP
#define MAT_OPS_UMAT_HPP
 
#include "MatOps.hpp"
 
namespace MatOps {
 
struct UMAT {
    UMAT() = delete;
};
struct UMatOpData : public MatOpsData::UserData {
  UMatOpData(ForcesAndSourcesCore::UserDataOperator *raw_op_ptr)
      : rawOpPtr(raw_op_ptr) {}
  ForcesAndSourcesCore::UserDataOperator *rawOpPtr;
};
 
template <int DIM, int MODEL_TYPE = (DIM == 3) ? MODEL_3D : MODEL_2D_PLANE_STRAIN>
struct MatUmatImpl;
 
template <int DIM, int MODEL_TYPE = (DIM == 3) ? MODEL_3D : MODEL_2D_PLANE_STRAIN>
struct StrainHandler {
    virtual MoFEMErrorCode runUmat(MatUmatImpl<DIM, MODEL_TYPE> &matUmatImpl,
                                 EntityHandle entity, int gg) = 0;
  virtual MoFEMErrorCode
    convertStress(MatUmatImpl<DIM, MODEL_TYPE> &matUmatImpl,
                boost::shared_ptr<MatOpsData> mat_ops_data_ptr) = 0;
  virtual MoFEMErrorCode
    convertTangent(MatUmatImpl<DIM, MODEL_TYPE> &matUmatImpl,
                 boost::shared_ptr<MatOpsData> mat_ops_data_ptr) = 0;
  virtual ~StrainHandler() = default;
};
 
template <>
boost::shared_ptr<PhysicalEquations>
createMatOpsPhysicalEquationsPtr<UMAT, MODEL_3D>(
    boost::shared_ptr<MatOpsData> mat_ops_data_ptr, int tag);
 
template <>
boost::shared_ptr<PhysicalEquations>
createMatOpsPhysicalEquationsPtr<UMAT, MODEL_2D_PLANE_STRAIN>(
    boost::shared_ptr<MatOpsData> mat_ops_data_ptr, int tag);
 
template <>
boost::shared_ptr<PhysicalEquations>
createMatOpsPhysicalEquationsPtr<UMAT, MODEL_2D_PLANE_STRESS>(
    boost::shared_ptr<MatOpsData> mat_ops_data_ptr, int tag);
 
template <>
boost::shared_ptr<PhysicalEquations>
createMatOpsPhysicalEquationsPtr<UMAT, MODEL_AXISYMMETRIC>(
    boost::shared_ptr<MatOpsData> mat_ops_data_ptr, int tag);
 
} // namespace MatOps
 
#endif // MAT_OPS_UMAT_HPP