EshelbianPlasticity::OpConstrainBoundaryHDivLhs_dU Struct Reference

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

Inheritance diagram for EshelbianPlasticity::OpConstrainBoundaryHDivLhs_dU:

Collaboration diagram for EshelbianPlasticity::OpConstrainBoundaryHDivLhs_dU:

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

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

Detailed Description

Examples

EshelbianPlasticity.cpp.

Definition at line 177 of file EshelbianContact.hpp.

Constructor & Destructor Documentation

OpConstrainBoundaryHDivLhs_dU()

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

Definition at line 773 of file EshelbianContact.cpp.

    : ContactOps::AssemblyBoundaryEleOp(row_field_name, col_field_name,
                                        ContactOps::BoundaryEleOp::OPROWCOL),
      commonDataPtr(common_data_ptr), contactTreePtr(contact_tree_ptr) {
  sYmm = false;
}

Member Function Documentation

iNtegrate()

MoFEMErrorCode EshelbianPlasticity::OpConstrainBoundaryHDivLhs_dU::iNtegrate	(	EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

Definition at line 784 of file EshelbianContact.cpp.

                                                                             {
  MoFEMFunctionBegin;
 
  using namespace ContactOps;
 
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
 
  auto nb_rows = row_data.getIndices().size();
  auto nb_cols = col_data.getIndices().size();
 
  auto &locMat = AssemblyBoundaryEleOp::locMat;
  locMat.resize(nb_rows, nb_cols, false);
  locMat.clear();
 
  if (nb_cols && nb_rows) {
 
    const size_t nb_gauss_pts = getGaussPts().size2();
 
    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;
    };
 
    constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
 
    auto face_data_vec_ptr =
        contactTreePtr->findFaceDataVecPtr(getFEEntityHandle());
    auto face_gauss_pts_it = face_data_vec_ptr->begin();
 
    auto t_row_base = row_data.getFTensor1N<3>();
    auto nb_face_functions = row_data.getN().size2() / 3;
    for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
 
      const auto alpha = t_w / 2.;
 
      size_t rr = 0;
      for (; rr != nb_rows / 3; ++rr) {
 
        auto row_base = alpha * (t_row_base(i) * t_material_normal(i));
 
        auto t_mat = getFTensor2FromArray<3, 3, 3>(locMat, 3 * rr);
        auto t_col_base = col_data.getFTensor0N(gg, 0);
 
        for (size_t cc = 0; cc != nb_cols / 3; ++cc) {
          const auto beta = row_base * t_col_base;
          t_mat(i, j) -= beta * t_kd(i, j);
          ++t_col_base;
          ++t_mat;
        }
 
        ++t_row_base;
      }
      for (; rr < nb_face_functions; ++rr)
        ++t_row_base;
 
      next();
    }
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

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

private

Definition at line 189 of file EshelbianContact.hpp.

contactTreePtr

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

private

Definition at line 190 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