OpPlasticTools::OpCalculateConstraintLhs_dEP Struct Reference

#include <users_modules/multifield_plasticity/src/PlasticOperators.hpp>

Inheritance diagram for OpPlasticTools::OpCalculateConstraintLhs_dEP:

Collaboration diagram for OpPlasticTools::OpCalculateConstraintLhs_dEP:

Public Member Functions
	OpCalculateConstraintLhs_dEP (const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< CommonData > common_data_ptr)
MoFEMErrorCode	iNtegrate (EntData &row_data, EntData &col_data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

Definition at line 185 of file PlasticOperators.hpp.

Constructor & Destructor Documentation

OpCalculateConstraintLhs_dEP()

OpPlasticTools::OpCalculateConstraintLhs_dEP::OpCalculateConstraintLhs_dEP	(	const std::string	row_field_name,
		const std::string	col_field_name,
		boost::shared_ptr< CommonData >	common_data_ptr
	)

Definition at line 931 of file PlasticOperators.cpp.

    : DomainEleOpAssembly(row_field_name, col_field_name,
                          DomainEleOp::OPROWCOL),
      commonDataPtr(common_data_ptr) {
  sYmm = false;
}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpPlasticTools::OpCalculateConstraintLhs_dEP::iNtegrate	(	EntData &	row_data,
		EntData &	col_data
	)

Definition at line 940 of file PlasticOperators.cpp.

                                                                          {
  MoFEMFunctionBegin;
 
  const size_t nb_row_dofs = row_data.getIndices().size();
  const size_t nb_col_dofs = col_data.getIndices().size();
  if (nb_row_dofs && nb_col_dofs) {
 
    auto get_dt = [&]() {
      double dt;
      CHKERR TSGetTimeStep(getFEMethod()->ts, &dt);
      dt = 1;
      return dt;
    };
    const auto dt = get_dt();
 
    const size_t nb_integration_pts = row_data.getN().size1();
    const size_t nb_row_base_functions = row_data.getN().size2();
    auto t_w = getFTensor0IntegrationWeight();
 
    auto get_row_base = [&]() {
      if (commonDataPtr->isDualBase) {
        double *base_ptr = &*commonDataPtr->dualBaseMat.data().begin();
        return Tensor0<PackPtr<double *, 1>>(base_ptr);
      } else {
        return row_data.getFTensor0N();
      }
    };
    auto t_row_base = get_row_base();
 
    auto t_f = getFTensor0FromVec(*(commonDataPtr->plasticSurfacePtr));
    auto t_tau = getFTensor0FromVec(*(commonDataPtr->plasticTauPtr));
    auto t_tau_dot = getFTensor0FromVec(*(commonDataPtr->plasticTauDotPtr));
    auto t_flow =
        getFTensor2SymmetricFromMat<3>(*(commonDataPtr->plasticFlowPtr));
 
    auto t_D = getFTensor4DdgFromMat<3, 3, 0>(*commonDataPtr->mtD);
    auto t_D_Deviator =
        getFTensor4DdgFromMat<3, 3, 0>(*commonDataPtr->mtD_Deviator);
 
    auto t_diff_symmetrize = diff_symmetrize();
    Index<'L', 6> L;
    auto t_L = symm_L_tensor();
 
    for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
      double alpha = dt * getMeasure() * t_w * (*cache).scale_constraint;
 
      auto mat_ptr = locMat.data().begin();
      auto t_diff_constrain_dstrain = diff_constrain_kin_hard_dstrain(
          t_D_Deviator,
          diff_constrain_dstress(
              diff_constrain_df(t_tau_dot, t_f, hardening(t_tau)), t_flow));
 
      Tensor1<PackPtr<double *, 6>, 6> t_mat(&locMat(0, 0), &locMat(0, 1),
                                             &locMat(0, 2), &locMat(0, 3),
                                             &locMat(0, 4), &locMat(0, 5));
 
      size_t rr = 0;
      for (; rr != nb_row_dofs; ++rr) {
 
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        for (size_t cc = 0; cc != nb_col_dofs / 6; cc++) {
 
          t_mat(L) -= (alpha * t_row_base * t_col_base) *
                      (t_diff_constrain_dstrain(i, j) * t_L(i, j, L));
 
          ++t_mat;
          ++t_col_base;
        }
 
        ++t_row_base;
      }
      for (; rr != nb_row_base_functions; ++rr)
        ++t_row_base;
 
      ++t_f;
      ++t_tau;
      ++t_tau_dot;
      ++t_flow;
      ++t_w;
    }
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> OpPlasticTools::OpCalculateConstraintLhs_dEP::commonDataPtr

private

Definition at line 192 of file PlasticOperators.hpp.

---

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

• users_modules/multifield_plasticity/src/PlasticOperators.hpp
• users_modules/multifield_plasticity/src/impl/PlasticOperators.cpp