EshlabianCohesive.hpp

Go to the documentation of this file.

/**
 * @file EshlabianCohesive.hpp
 * @author your name (you@domain.com)
 * @brief Add cohesive elements to Eshelbian plasticity module
 * @version 0.1
 * @date 2025-12-14
 *
 * @copyright Copyright (c) 2025
 *
 */
 
#ifndef MOFEM_ESHELBIAN_COHESIVE_HPP
#define MOFEM_ESHELBIAN_COHESIVE_HPP
 
namespace EshelbianPlasticity {
 
MoFEMErrorCode pushCohesiveOpsRhs(
    EshelbianCore &ep,
    ForcesAndSourcesCore::GaussHookFun set_integration_at_front_face,
    boost::shared_ptr<Range> interface_range_ptr,
    boost::ptr_deque<ForcesAndSourcesCore::UserDataOperator> &pip);
 
MoFEMErrorCode pushCohesiveOpsLhs(
    EshelbianCore &ep,
    ForcesAndSourcesCore::GaussHookFun set_integration_at_front_face,
    boost::shared_ptr<Range> interface_range_ptr,
    boost::ptr_deque<ForcesAndSourcesCore::UserDataOperator> &pip);
 
MoFEMErrorCode initializeCohesiveKappaField(EshelbianCore &ep);
 
struct CohesiveTAOCtx {
  CohesiveTAOCtx() = default;
  virtual SmartPetscObj<Vec> duplicateKappaVec() = 0;
  virtual SmartPetscObj<Vec> duplicateDissipationVec() = 0;
  virtual SmartPetscObj<Vec> duplicateGradientVec() = 0;
  virtual CommInterface::EntitiesPetscVector &getKappaVec() = 0;
};
 
boost::shared_ptr<CohesiveTAOCtx> createCohesiveTAOCtx(
    EshelbianCore *ep,
    ForcesAndSourcesCore::GaussHookFun set_integration_at_front_face,
    SmartPetscObj<TS> time_solver);
 
MoFEMErrorCode cohesiveEvaluateObjectiveAndGradient(Tao tao, Vec x,
                                                    PetscReal *f, Vec g,
                                                    void *ctx);
Tag get_kappa_tag(moab::Interface &moab);
Tag get_delta_kappa_tag(moab::Interface &moab);
 
} // namespace EshelbianPlasticity
 
#endif // MOFEM_ESHELBIAN_COHESIVE_HPP