EshelbianPlasticity::OpCohesive_dJ_dP Struct Reference

Inheritance diagram for EshelbianPlasticity::OpCohesive_dJ_dP:

Collaboration diagram for EshelbianPlasticity::OpCohesive_dJ_dP:

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

Member Typedef Documentation

OP

using EshelbianPlasticity::OpCohesive_dJ_dP::OP = OpCohesiveRhs

Definition at line 733 of file EshelbianCohesive.cpp.

Member Function Documentation

aSsemble()

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

inline

Definition at line 811 of file EshelbianCohesive.cpp.

                                         {
    MoFEMFunctionBegin;
    return VecSetValues<AssemblyTypeSelector<A>>(vec_dJdu, data, OP::locF,
                                                 ADD_VALUES);
    MoFEMFunctionReturn(0);
  }

iNtegrate()

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

inline

Definition at line 736 of file EshelbianCohesive.cpp.

                                          {
    MoFEMFunctionBegin;
 
    FTENSOR_INDEX(3, i);
    FTENSOR_INDEX(3, J);
 
    int nb_integration_pts = OP::getGaussPts().size2();
    auto nb_dofs = data.getIndices().size();
 
    int nb_base_functions = data.getN().size2() / 3;
    auto t_row_base_fun = data.getFTensor1N<3>();
    auto t_w = OP::getFTensor0IntegrationWeight();
 
    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();
 
    auto next = [&]() {
      ++t_P;
      ++t_face_normal;
      ++t_kappa;
      ++t_delta_kappa;
      ++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 dJ_dtraction = [](auto &t_traction, auto &t_normalized_normal,
                             auto &t_delta_kappa, auto &t_kappa, auto gc) {
        double kappa = static_cast<double>(t_kappa);
        double delta_kappa = static_cast<double>(t_delta_kappa);
        FTENSOR_INDEX(3, i);
        FTensor::Tensor1<double, 3> t_traction_double;
        t_traction_double(i) = t_traction(i);
        auto t_dM =
            GriffithCohesiveLaw::calculateDissipationSurplusDiffTraction(
                t_traction_double, delta_kappa, kappa, t_normalized_normal, gc,
                OpGetParameters::beta, OpGetParameters::min_stiffness);
        FTensor::Tensor1<double, 3> t_dJ_double;
        t_dJ_double(i) = t_dM(i);
        return t_dJ_double;
      };
 
      auto assemble = [&](const FTensor::Tensor1<double, 3> &&t_dJ) {
 
        auto t_nf = getFTensor1FromPtr<3>(&*OP::locF.begin());
        int bb = 0;
        for (; bb != nb_dofs / SPACE_DIM; ++bb) {
          double row_base = t_w * (t_row_base_fun(J) * t_normal(J));
          t_nf(i) += row_base * t_dJ(i);
          ++t_nf;
          ++t_row_base_fun;
        }
        for (; bb != nb_base_functions; ++bb)
          ++t_row_base_fun;
      };
 
      assemble(dJ_dtraction(t_traction, t_normalized_normal, t_delta_kappa,
                            t_kappa, *gcPtr));
 
      // cerr << "AAAA "  << OP::locF << endl;
      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