OpContactTools::OpConstrainBoundaryLhs_dTraction Struct Reference

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

Inheritance diagram for OpContactTools::OpConstrainBoundaryLhs_dTraction:

Collaboration diagram for OpContactTools::OpConstrainBoundaryLhs_dTraction:

Public Member Functions
	OpConstrainBoundaryLhs_dTraction (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 63 of file ContactOperators.hpp.

Constructor & Destructor Documentation

OpConstrainBoundaryLhs_dTraction()

OpContactTools::OpConstrainBoundaryLhs_dTraction::OpConstrainBoundaryLhs_dTraction	(	const std::string	row_field_name,
		const std::string	col_field_name,
		boost::shared_ptr< CommonData >	common_data_ptr
	)

Definition at line 416 of file ContactOperators.hpp.

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

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpContactTools::OpConstrainBoundaryLhs_dTraction::iNtegrate	(	EntData &	row_data,
		EntData &	col_data
	)

Definition at line 423 of file ContactOperators.hpp.

                                          {
  MoFEMFunctionBegin;
 
  const size_t nb_gauss_pts = getGaussPts().size2();
  const size_t row_nb_dofs = row_data.getIndices().size();
  const size_t col_nb_dofs = col_data.getIndices().size();
 
  if (row_nb_dofs && col_nb_dofs) {
 
    auto t_normal = getFTensor1Normal();
    const double ls = sqrt(t_normal(i) * t_normal(i));
    t_normal(i) /= ls;
 
    auto t_disp = getFTensor1FromMat<3>(*(commonDataPtr->mDispPtr));
    auto t_traction =
        getFTensor1FromMat<3>(*(commonDataPtr->contactTractionPtr));
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_row_base = row_data.getFTensor1N<3>();
    size_t nb_face_functions = row_data.getN().size2() / 3;
 
    auto t_contact_normal =
        getFTensor1FromMat<3>(*(commonDataPtr->contactNormalPtr));
    auto t_gap = getFTensor0FromVec(*(commonDataPtr->contactGapPtr));
 
    for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
 
      const double alpha = t_w * getMeasure();
 
      if (getNormalsAtGaussPts().size1() == nb_gauss_pts) {
        VectorDouble n = getNormalsAtGaussPts(gg);
        auto t_n = getFTensor1FromPtr<3>(&*n.data().begin());
        t_normal(i) = t_n(i) / sqrt(t_n(j) * t_n(j));
      }
 
      // auto t_contact_normal = normal(t_coords, t_disp);
      Tensor2<double, 3, 3> t_P;
      t_P(i, j) = t_contact_normal(i) * t_contact_normal(j);
 
      Tensor2<double, 3, 3> t_Q;
      t_Q(i, j) = kronecker_delta(i, j) - t_P(i, j);
 
      const double diff_traction = diff_constraints_dtraction(
          t_gap, normal_traction(t_traction, t_contact_normal));
 
      size_t rr = 0;
      for (; rr != row_nb_dofs / 3; ++rr) {
        auto t_mat = getFTensor2FromArray<3, 3, 3>(locMat, 3 * rr);
        const double row_base = t_row_base(i) * t_normal(i);
 
        auto t_col_base = col_data.getFTensor1N<3>(gg, 0);
 
        for (size_t cc = 0; cc != col_nb_dofs / 3; ++cc) {
          const double col_base = t_col_base(i) * t_normal(i);
          const double beta = alpha * row_base * col_base;
 
          t_mat(i, j) += (beta * diff_traction) * t_P(i, j);
          t_mat(i, j) += beta * (*cache).cn_cont * t_Q(i, j);
 
          ++t_col_base;
          ++t_mat;
        }
 
        ++t_row_base;
      }
      for (; rr < nb_face_functions; ++rr)
        ++t_row_base;
 
      ++t_contact_normal;
      ++t_gap;
      ++t_disp;
      ++t_traction;
      ++t_coords;
      ++t_w;
    }
 
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> OpContactTools::OpConstrainBoundaryLhs_dTraction::commonDataPtr

private

Definition at line 71 of file ContactOperators.hpp.

---

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

• users_modules/multifield_plasticity/src/ContactOperators.hpp