EshelbianPlasticity::OpConstrainBoundaryHDivRhs Struct Reference

#include <users_modules/eshelbian_plasticit/src/EshelbianContact.hpp>

Inheritance diagram for EshelbianPlasticity::OpConstrainBoundaryHDivRhs:

Collaboration diagram for EshelbianPlasticity::OpConstrainBoundaryHDivRhs:

Public Member Functions
	OpConstrainBoundaryHDivRhs (const std::string row_field_name, boost::shared_ptr< ContactOps::CommonData > common_data_ptr, boost::shared_ptr< ContactTree > contact_tree_ptr)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data)

Private Attributes
boost::shared_ptr< ContactOps::CommonData >	commonDataPtr
boost::shared_ptr< ContactTree >	contactTreePtr

Detailed Description

Examples

EshelbianPlasticity.cpp.

Definition at line 129 of file EshelbianContact.hpp.

Constructor & Destructor Documentation

OpConstrainBoundaryHDivRhs()

EshelbianPlasticity::OpConstrainBoundaryHDivRhs::OpConstrainBoundaryHDivRhs	(	const std::string	row_field_name,
		boost::shared_ptr< ContactOps::CommonData >	common_data_ptr,
		boost::shared_ptr< ContactTree >	contact_tree_ptr
	)

Definition at line 426 of file EshelbianContact.cpp.

    : ContactOps::AssemblyBoundaryEleOp(row_field_name, row_field_name,
                                        ContactOps::BoundaryEleOp::OPROW),
      commonDataPtr(common_data_ptr), contactTreePtr(contact_tree_ptr) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode EshelbianPlasticity::OpConstrainBoundaryHDivRhs::iNtegrate

(

EntitiesFieldData::EntData &

row_data

)

Definition at line 435 of file EshelbianContact.cpp.

                                                                    {
  MoFEMFunctionBegin;
 
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
  FTensor::Index<'l', 3> l;
 
  const size_t nb_gauss_pts = getGaussPts().size2();
 
#ifndef NDEBUG
  if (commonDataPtr->contactDisp.size2() != nb_gauss_pts) {
    SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Wrong number of integration pts %d != %d",
             commonDataPtr->contactDisp.size2(), nb_gauss_pts);
  }
#endif // !NDEBUG
 
  auto &nf = locF;
  locF.clear();
 
  auto t_w = getFTensor0IntegrationWeight();
  auto t_disp = getFTensor1FromMat<3>(commonDataPtr->contactDisp);
  auto t_traction = getFTensor1FromMat<3>(commonDataPtr->contactTraction);
  auto t_coords = getFTensor1CoordsAtGaussPts();
  auto t_material_normal = getFTensor1NormalsAtGaussPts();
 
  auto next = [&]() {
    ++t_w;
    ++t_disp;
    ++t_traction;
    ++t_coords;
    ++t_material_normal;
  };
 
  auto face_data_vec_ptr =
      contactTreePtr->findFaceDataVecPtr(getFEEntityHandle());
  auto face_gauss_pts_it = face_data_vec_ptr->begin();
 
  auto nb_base_functions = data.getN().size2() / 3;
  auto t_base = data.getFTensor1N<3>();
  for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
 
    auto t_nf = getFTensor1FromPtr<3>(&nf[0]);
    const double alpha = t_w / 2.;
 
    size_t bb = 0;
    for (; bb != nbRows / 3; ++bb) {
      const double beta = alpha * t_base(i) * t_material_normal(i);
      t_nf(i) -= beta * t_disp(i);
      ++t_nf;
      ++t_base;
    }
    for (; bb < nb_base_functions; ++bb)
      ++t_base;
 
    next();
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

boost::shared_ptr<ContactOps::CommonData> EshelbianPlasticity::OpConstrainBoundaryHDivRhs::commonDataPtr

private

Definition at line 139 of file EshelbianContact.hpp.

contactTreePtr

boost::shared_ptr<ContactTree> EshelbianPlasticity::OpConstrainBoundaryHDivRhs::contactTreePtr

private

Definition at line 140 of file EshelbianContact.hpp.

---

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

• users_modules/eshelbian_plasticit/src/EshelbianContact.hpp
• users_modules/eshelbian_plasticit/src/impl/EshelbianContact.cpp