EshelbianPlasticity::OpCohesiveLhs_dPdP Struct Reference

Inheritance diagram for EshelbianPlasticity::OpCohesiveLhs_dPdP:

Collaboration diagram for EshelbianPlasticity::OpCohesiveLhs_dPdP:

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 521 of file EshelbianCohesive.cpp.

Member Typedef Documentation

OP

using EshelbianPlasticity::OpCohesiveLhs_dPdP::OP = OpCohesiveRhs

Definition at line 523 of file EshelbianCohesive.cpp.

Member Function Documentation

aSsemble()

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

inline

Definition at line 606 of file EshelbianCohesive.cpp.

                                         {
    MoFEMFunctionBegin;
 
    if (!this->timeScalingFun.empty())
      this->locMat *= this->timeScalingFun(this->getFEMethod()->ts_t);
    if (!this->feScalingFun.empty())
      this->locMat *= this->feScalingFun(this->getFEMethod());
    // assemble local matrix
    CHKERR this->matSetValuesHook(this, data, data, this->locMat);
 
    MoFEMFunctionReturn(0);
  }

iNtegrate()

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

inline

Definition at line 526 of file EshelbianCohesive.cpp.

                                          {
    MoFEMFunctionBegin;
 
    auto get_sense_index = [this]() { return (faceSense == 1) ? 0 : 1; };
 
    auto nb_dofs = data.getIndices().size();
    FTENSOR_INDEX(3, i);
    FTENSOR_INDEX(3, j);
    FTENSOR_INDEX(3, J);
 
    int nb_integration_pts = OP::getGaussPts().size2();
 
    auto t_P = getFTensor2FromMat<3, 3>(*(OP::fluxMatPtr));
    auto t_kappa = getFTensor0FromVec<0>(*kappaPtr);
    auto t_delta_kappa = getFTensor0FromVec<0>(*kappaDeltaPtr);
    auto t_face_normal = getFTensor1NormalsAtGaussPts();
 
    int nb_base_functions = data.getN().size2() / 3;
    auto t_row_base_fun = data.getFTensor1N<3>();
    auto t_w = OP::getFTensor0IntegrationWeight();
 
    auto next = [&]() {
      ++t_P;
      ++t_kappa;
      ++t_delta_kappa;
      ++t_face_normal;
      ++t_w;
    };
 
    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 fracture = [this](auto &t_traction, auto &t_normalized_normal,
                             auto &t_kappa, auto &t_delta_kappa, auto gc) {
        double kappa = static_cast<double>(t_kappa + t_delta_kappa);
        auto alpha = GriffithCohesiveLaw::getAlpha(
            kappa, gc, OpGetParameters::min_stiffness);
        auto t_dgap = GriffithCohesiveLaw::calculateDiffGapDTraction(
            t_traction, t_normalized_normal, alpha, OpGetParameters::beta,
            true);
        FTENSOR_INDEX(3, i);
        FTENSOR_INDEX(3, j);
        FTensor::Tensor2<double, 3, 3> t_dgap_double;
        t_dgap_double(i, j) = -t_dgap(i, j) / 2.;
        return t_dgap_double;
      };
 
      auto assemble = [&](const FTensor::Tensor2<double, 3, 3> &&t_dgap) {
        int rr = 0;
        for (; rr != nb_dofs / SPACE_DIM; ++rr) {
          auto t_mat = getFTensor2FromArray<3, 3, 3>(OP::locMat, 3 * rr);
          double row_base = t_w * (t_row_base_fun(J) * t_normal(J));
          auto t_col_base_fun = data.getFTensor1N<3>(gg, 0);
          for (int cc = 0; cc != nb_dofs / SPACE_DIM; ++cc) {
            double col_base = t_col_base_fun(J) * t_normalized_normal(J);
            t_mat(i, j) += (row_base * col_base) * t_dgap(i, j);
            ++t_mat;
            ++t_col_base_fun;
          }
          ++t_row_base_fun;
        }
        for (; rr != nb_base_functions; ++rr)
          ++t_row_base_fun;
      };
 
      assemble(fracture(t_traction, t_normalized_normal, t_kappa, t_delta_kappa,
                        *gcPtr));
 
      next();
    }
 
    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 421 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