OpRotatingFrameTools::OpCalculateJumpOnSkeleton Struct Reference

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

Inheritance diagram for OpRotatingFrameTools::OpCalculateJumpOnSkeleton:

Collaboration diagram for OpRotatingFrameTools::OpCalculateJumpOnSkeleton:

Public Member Functions
	OpCalculateJumpOnSkeleton (const std::string field_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< DomainSideEle > side_op_fe)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr
boost::shared_ptr< DomainSideEle >	sideOpFe

Detailed Description

Definition at line 90 of file RotatingFrameOperators.hpp.

Constructor & Destructor Documentation

OpCalculateJumpOnSkeleton()

OpRotatingFrameTools::OpCalculateJumpOnSkeleton::OpCalculateJumpOnSkeleton	(	const std::string	field_name,
		boost::shared_ptr< CommonData >	common_data_ptr,
		boost::shared_ptr< DomainSideEle >	side_op_fe
	)

Definition at line 383 of file RotatingFrameOperators.cpp.

    : SkeletonEleOp(field_name, field_name, SkeletonEleOp::OPROW),
      commonDataPtr(common_data_ptr), sideOpFe(side_op_fe) {
  std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
  doEntities[MBTRI] = doEntities[MBQUAD] = true;
}

Member Function Documentation

doWork()

MoFEMErrorCode OpRotatingFrameTools::OpCalculateJumpOnSkeleton::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Definition at line 392 of file RotatingFrameOperators.cpp.

                                                                {
  MoFEMFunctionBegin;
  // if (skipEnts.find(getFEEntityHandle()) != skipEnts.end())
  //   MoFEMFunctionReturnHot(0);
  if ((type == MBTRI || type == MBQUAD) && side == 0) {
    const size_t nb_dofs = data.getIndices().size();
    // if (nb_dofs == 0)
    //   MoFEMFunctionReturnHot(0);
 
    const size_t nb_integration_pts = data.getN().size1();
    const size_t nb_base_functions = data.getN().size2();
    commonDataPtr->plasticStrainSideMap.clear();
    CHKERR loopSideVolumes("dFE", *sideOpFe);
    const int check_size = commonDataPtr->plasticStrainSideMap.size();
 
    if (commonDataPtr->plasticStrainSideMap.size() == 1)
      if (static_cast<int>(
              commonDataPtr->plasticStrainSideMap.begin()->second.size2()) !=
          nb_integration_pts)
        SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
                "wrong number of integration points");
 
    auto &plastic_mat_l = commonDataPtr->plasticStrainSideMap.at(-1);
    auto &plastic_mat_r = commonDataPtr->plasticStrainSideMap.at(1);
    auto t_plastic_strain_l = getFTensor2SymmetricFromMat<3>(plastic_mat_l);
    auto t_plastic_strain_r = getFTensor2SymmetricFromMat<3>(plastic_mat_r);
    auto &tau_vec_l = commonDataPtr->plasticTauSideMap.at(-1);
    auto &tau_vec_r = commonDataPtr->plasticTauSideMap.at(1);
    auto t_tau_l = getFTensor0FromVec(tau_vec_l);
    auto t_tau_r = getFTensor0FromVec(tau_vec_r);
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_base = data.getFTensor0N();
    auto t_coords = getFTensor1CoordsAtGaussPts();
    commonDataPtr->plasticTauJumpPtr->resize(nb_integration_pts, false);
    commonDataPtr->plasticStrainJumpPtr->resize(6, nb_integration_pts, false);
    auto t_plastic_strain_jump =
        getFTensor2SymmetricFromMat<3>(*commonDataPtr->plasticStrainJumpPtr);
    auto t_tau_jump = getFTensor0FromVec(*commonDataPtr->plasticTauJumpPtr);
 
    for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
      double alpha = getMeasure() * t_w;
 
      // Tensor2_symmetric<PackPtr<double *, 6>, 3> t_nf{&nf[0], &nf[1],
      // &nf[2],
      //                                                 &nf[3], &nf[4],
      //                                                 &nf[5]};
      // this is just work in progress
      t_plastic_strain_jump(i, j) =
          t_plastic_strain_l(i, j) - t_plastic_strain_r(i, j);
      t_tau_jump = t_tau_l - t_tau_r;
 
      // size_t bb = 0;
      // for (; bb != nb_dofs / 6; ++bb) {
      //   t_nf(i, j) += alpha * t_base;
      //   ++t_base;
      //   ++t_nf;
      // }
 
      // for (; bb < nb_base_functions; ++bb)
      //   ++t_base;
 
      ++t_plastic_strain_jump;
      ++t_tau_jump;
      ++t_plastic_strain_l;
      ++t_plastic_strain_r;
      ++t_tau_l;
      ++t_tau_r;
      ++t_w;
    }
 
    // CHKERR VecSetValues(getTSf(), data, nf.data(), ADD_VALUES);
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> OpRotatingFrameTools::OpCalculateJumpOnSkeleton::commonDataPtr

private

Definition at line 97 of file RotatingFrameOperators.hpp.

sideOpFe

boost::shared_ptr<DomainSideEle> OpRotatingFrameTools::OpCalculateJumpOnSkeleton::sideOpFe

private

Definition at line 98 of file RotatingFrameOperators.hpp.

---

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

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