EshelbianPlasticity::OpCohesive_dJ_dkappa Struct Reference

Inheritance diagram for EshelbianPlasticity::OpCohesive_dJ_dkappa:

Collaboration diagram for EshelbianPlasticity::OpCohesive_dJ_dkappa:

Public Types
using	OP = OpCohesiveRhs
Public Types inherited from EshelbianPlasticity::OpCohesiveRhs
using	OP = OpBrokenBaseCohesive
Public Types inherited from EshelbianPlasticity::OpBrokenBaseCohesive
using	OP = FormsIntegrators< FaceElementForcesAndSourcesCore::UserDataOperator >::Assembly< A >::OpBrokenBase

Public Member Functions
MoFEMErrorCode	iNtegrate (EntData &data)
MoFEMErrorCode	aSsemble (EntData &data)
	OpCohesiveRhs (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< double > gc_ptr, boost::shared_ptr< VectorDouble > kappa_ptr, boost::shared_ptr< VectorDouble > kappa_delta_ptr, boost::shared_ptr< std::array< MatrixDouble, 2 > > lambda_ptr=nullptr, Tag dissipation_tags=Tag(), Tag grad_dissipation_tags=Tag(), SmartPetscObj< Vec > vec_dJ_dx=SmartPetscObj< Vec >(), boost::shared_ptr< Range > ents_ptr=nullptr)
Public Member Functions inherited from EshelbianPlasticity::OpCohesiveRhs
	OpCohesiveRhs (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< double > gc_ptr, boost::shared_ptr< VectorDouble > kappa_ptr, boost::shared_ptr< VectorDouble > kappa_delta_ptr, boost::shared_ptr< std::array< MatrixDouble, 2 > > lambda_ptr=nullptr, Tag dissipation_tags=Tag(), Tag grad_dissipation_tags=Tag(), SmartPetscObj< Vec > vec_dJ_dx=SmartPetscObj< Vec >(), boost::shared_ptr< Range > ents_ptr=nullptr)
MoFEMErrorCode	iNtegrate (EntData &data)
Public Member Functions inherited from EshelbianPlasticity::OpBrokenBaseCohesive
	OpBrokenBaseCohesive (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_flux_data_ptr, boost::shared_ptr< Range > ents_ptr=nullptr)
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)

Additional Inherited Members
Protected Attributes inherited from EshelbianPlasticity::OpCohesiveRhs
boost::shared_ptr< MatrixDouble >	uGammaPtr
boost::shared_ptr< double >	gcPtr
boost::shared_ptr< VectorDouble >	kappaPtr
boost::shared_ptr< VectorDouble >	kappaDeltaPtr
boost::shared_ptr< std::array< MatrixDouble, 2 > >	lambdaPtr
boost::shared_ptr< double >	totalDissipation
boost::shared_ptr< double >	tatalDissipationGrad
SmartPetscObj< Vec >	vec_dJdu
Tag	dissipationTag
Tag	gradDissipationTag
Protected Attributes inherited from EshelbianPlasticity::OpBrokenBaseCohesive
boost::shared_ptr< MatrixDouble >	fluxMatPtr
int	faceSense = 0

Detailed Description

Definition at line 571 of file EshelbianCohesive.cpp.

Member Typedef Documentation

OP

using EshelbianPlasticity::OpCohesive_dJ_dkappa::OP = OpCohesiveRhs

Definition at line 573 of file EshelbianCohesive.cpp.

Member Function Documentation

aSsemble()

MoFEMErrorCode EshelbianPlasticity::OpCohesive_dJ_dkappa::aSsemble ( EntData &  data )

inline

Definition at line 672 of file EshelbianCohesive.cpp.

                                         {
    MoFEMFunctionBegin;
    MoFEMFunctionReturn(0);
  }

iNtegrate()

MoFEMErrorCode EshelbianPlasticity::OpCohesive_dJ_dkappa::iNtegrate ( EntData &  data )

inline

Definition at line 576 of file EshelbianCohesive.cpp.

                                          {
    MoFEMFunctionBegin;
 
    auto get_sense_index = [this]() { return (faceSense == 1) ? 0 : 1; };
 
    FTENSOR_INDEX(3, i);
    FTENSOR_INDEX(3, J);
    FTENSOR_INDEX(3, K);
 
    int nb_integration_pts = OP::getGaussPts().size2();
 
    auto t_P = getFTensor2FromMat<3, 3>(*(OP::fluxMatPtr));
    auto t_lambda = getFTensor2FromMat<3, 3>(lambdaPtr->at(get_sense_index()));
    auto t_kappa = getFTensor0FromVec<0>(*kappaPtr);
    auto t_delta_kappa = getFTensor0FromVec<0>(*kappaDeltaPtr);
    auto t_face_normal = getFTensor1NormalsAtGaussPts();
    auto t_w = OP::getFTensor0IntegrationWeight();
 
    auto next = [&]() {
      ++t_P;
      ++t_face_normal;
      ++t_kappa;
      ++t_delta_kappa;
      ++t_lambda;
      ++t_w;
    };
 
    double face_dissipation = 0.0;
    double face_grad_dissipation = 0.0;
    auto face_handle = OP::getFEEntityHandle();
    CHKERR getMoab().tag_get_data(dissipationTag, &face_handle, 1,
                                  &face_dissipation);
    CHKERR getMoab().tag_get_data(gradDissipationTag, &face_handle, 1,
                                  &face_grad_dissipation);
 
    for (auto gg = 0; gg != nb_integration_pts; ++gg) {
      FTensor::Tensor1<double, 3> t_normal;
      t_normal(J) = t_face_normal(J) * (faceSense / 2.);
      FTensor::Tensor1<double, 3> t_normalized_normal;
      t_normalized_normal(J) = t_normal(J);
      t_normalized_normal.normalize();
      FTensor::Tensor1<double, 3> t_traction;
      t_traction(i) = t_P(i, J) * t_normalized_normal(J);
 
      auto dJ_dkappa = [](auto &t_delta_kappa, auto &t_traction,
                          auto &t_normalized_normal, auto &t_kappa, auto gc) {
        double kappa = static_cast<double>(t_kappa);
        double delta_kappa = static_cast<double>(t_delta_kappa);
        double teff = GriffithCohesiveLaw::calculateEffectiveTraction(
            t_traction, t_normalized_normal, OpGetParameters::beta);
        double m = GriffithCohesiveLaw::calculateDissipationSurplus(
            delta_kappa, teff, kappa, gc);
        double m_grad =
            GriffithCohesiveLaw::calculateDissipationSurplusDiffKappa(
                delta_kappa, teff, kappa, gc);
        return boost::make_tuple(m, m_grad);
      };
 
      auto dr_kappa = [](auto &t_delta_kappa, auto &t_traction,
                         auto &t_normalized_normal, auto &t_kappa, auto gc) {
        double kappa = static_cast<double>(t_kappa);
        double delta_kappa = static_cast<double>(t_delta_kappa);
        double kappa_plus_delta = kappa + delta_kappa;
        double diff_alpha =
            GriffithCohesiveLaw::getDiffAlpha(kappa_plus_delta, gc);
        auto [teff, t_gap] = GriffithCohesiveLaw::calculateGap(
            t_traction, t_normalized_normal, diff_alpha, OpGetParameters::beta);
        FTENSOR_INDEX(3, i);
        FTensor::Tensor1<double, 3> t_gap_double;
        t_gap_double(i) = -t_gap(i);
        return t_gap_double;
      };
 
      auto [J, dJ] = dJ_dkappa(t_delta_kappa, t_traction, t_normalized_normal,
                               t_kappa, *gcPtr);
      face_dissipation += t_w * J * t_normal.l2();
      face_grad_dissipation += t_w * dJ * t_normal.l2();
 
      auto t_dr = dr_kappa(t_delta_kappa, t_traction, t_normalized_normal,
                           t_kappa, *gcPtr);
      FTensor::Tensor1<double, 3> t_l;
      t_l(i) = t_lambda(i, K) * t_normal(K);
      face_grad_dissipation -= t_w * t_l(i) * t_dr(i);
 
      next();
    }
 
    CHKERR getMoab().tag_set_data(dissipationTag, &face_handle, 1,
                                  &face_dissipation);
    CHKERR getMoab().tag_set_data(gradDissipationTag, &face_handle, 1,
                                  &face_grad_dissipation);
 
 
    MoFEMFunctionReturn(0);
  };

OpCohesiveRhs()

EshelbianPlasticity::OpCohesiveRhs::OpCohesiveRhs ( boost::shared_ptr< std::vector< BrokenBaseSideData > >  broken_base_side_data, boost::shared_ptr< double >  gc_ptr, boost::shared_ptr< VectorDouble >  kappa_ptr, boost::shared_ptr< VectorDouble >  kappa_delta_ptr, boost::shared_ptr< std::array< MatrixDouble, 2 > >  lambda_ptr = nullptr, Tag  dissipation_tags = Tag(), Tag  grad_dissipation_tags = Tag(), SmartPetscObj< Vec >  vec_dJ_dx = SmartPetscObj<Vec>(), boost::shared_ptr< Range >  ents_ptr = nullptr  )

inline

Definition at line 372 of file EshelbianCohesive.cpp.

      : OpBrokenBaseCohesive(broken_base_side_data, ents_ptr), gcPtr(gc_ptr),
        kappaPtr(kappa_ptr), kappaDeltaPtr(kappa_delta_ptr),
        lambdaPtr(lambda_ptr), dissipationTag(dissipation_tags),
        gradDissipationTag(grad_dissipation_tags), vec_dJdu(vec_dJ_dx) {}

---

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

• users_modules/eshelbian_plasticity/src/impl/EshelbianCohesive.cpp