OpContactTools::OpConstrainBoundaryRhs Struct Reference

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

Inheritance diagram for OpContactTools::OpConstrainBoundaryRhs:

Collaboration diagram for OpContactTools::OpConstrainBoundaryRhs:

Public Member Functions
	OpConstrainBoundaryRhs (const std::string field_name, boost::shared_ptr< CommonData > common_data_ptr)
MoFEMErrorCode	iNtegrate (EntData &data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

Definition at line 43 of file ContactOperators.hpp.

Constructor & Destructor Documentation

OpConstrainBoundaryRhs()

OpContactTools::OpConstrainBoundaryRhs::OpConstrainBoundaryRhs	(	const std::string	field_name,
		boost::shared_ptr< CommonData >	common_data_ptr
	)

Definition at line 125 of file ContactOperators.hpp.

    : BoundaryEleOpAssembly(field_name, field_name, BoundaryEleOp::OPROW),
      commonDataPtr(common_data_ptr) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode OpContactTools::OpConstrainBoundaryRhs::iNtegrate

(

EntData &

data

)

Definition at line 130 of file ContactOperators.hpp.

                                                              {
  MoFEMFunctionBegin;
 
  const size_t nb_gauss_pts = getGaussPts().size2();
  const size_t nb_dofs = data.getIndices().size();
 
  if (nb_dofs) {
 
    auto t_normal = getFTensor1Normal();
    const double ls = sqrt(t_normal(i) * t_normal(i));
    t_normal(i) /= ls;
 
    auto t_w = getFTensor0IntegrationWeight();
    auto t_disp = getFTensor1FromMat<3>(*(commonDataPtr->mDispPtr));
    auto t_traction =
        getFTensor1FromMat<3>(*(commonDataPtr->contactTractionPtr));
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    size_t nb_base_functions = data.getN().size2() / 3;
    auto t_base = data.getFTensor1N<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) {
 
      auto t_nf = getFTensor1FromArray<3, 3>(locF);
 
      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));
      }
      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 gap_0 = gap0(t_gap, t_disp, t_contact_normal);
      Tensor1<double, 3> t_rhs_constrains;
      t_rhs_constrains(i) =
          t_contact_normal(i) *
          constraint(gap_0, t_gap,
                     normal_traction(t_traction, t_contact_normal));
 
      Tensor1<double, 3> t_rhs_tangent_disp, t_rhs_tangent_traction;
      t_rhs_tangent_disp(i) = t_Q(i, j) * t_disp(j);
      t_rhs_tangent_traction(i) = (*cache).cn_cont * t_Q(i, j) * t_traction(j);
 
      size_t bb = 0;
      for (; bb != nb_dofs / 3; ++bb) {
        const double beta = alpha * (t_base(i) * t_normal(i));
 
        t_nf(i) -= beta * t_rhs_constrains(i);
        t_nf(i) -= beta * t_rhs_tangent_disp(i);
        t_nf(i) += beta * t_rhs_tangent_traction(i);
 
        ++t_nf;
        ++t_base;
      }
      for (; bb < nb_base_functions; ++bb)
        ++t_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::OpConstrainBoundaryRhs::commonDataPtr

private

Definition at line 49 of file ContactOperators.hpp.

---

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

• users_modules/multifield_plasticity/src/ContactOperators.hpp