ContactOps::OpConstrainBoundaryRhs Struct Reference

#include <users_modules/tutorials/adv-1/src/ContactOps.hpp>

Inheritance diagram for ContactOps::OpConstrainBoundaryRhs:

Collaboration diagram for ContactOps::OpConstrainBoundaryRhs:

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

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

Definition at line 171 of file ContactOps.hpp.

Constructor & Destructor Documentation

OpConstrainBoundaryRhs()

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

Examples

ContactOps.hpp.

Definition at line 273 of file ContactOps.hpp.

    : AssemblyBoundaryEleOp(field_name, field_name, DomainEleOp::OPROW),
      commonDataPtr(common_data_ptr) {}

Member Function Documentation

iNtegrate()

MoFEMErrorCode ContactOps::OpConstrainBoundaryRhs::iNtegrate

(

EntitiesFieldData::EntData &

data

)

Examples

ContactOps.hpp.

Definition at line 279 of file ContactOps.hpp.

                                                                {
  MoFEMFunctionBegin;
 
  const size_t nb_gauss_pts = AssemblyBoundaryEleOp::getGaussPts().size2();
 
  auto &nf = AssemblyBoundaryEleOp::locF;
 
  auto t_normal = getFTensor1Normal();
  t_normal(i) /= sqrt(t_normal(j) * t_normal(j));
 
  auto t_w = getFTensor0IntegrationWeight();
  auto t_disp = getFTensor1FromMat<SPACE_DIM>(commonDataPtr->contactDisp);
  auto t_traction =
      getFTensor1FromMat<SPACE_DIM>(commonDataPtr->contactTraction);
  auto t_coords = getFTensor1CoordsAtGaussPts();
 
  size_t nb_base_functions = data.getN().size2() / 3;
  auto t_base = data.getFTensor1N<3>();
  for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
 
    auto t_nf = getFTensor1FromPtr<SPACE_DIM>(&nf[0]);
    const double alpha = t_w * getMeasure();
 
    FTensor::Tensor1<double, 3> t_spatial_coords{0., 0., 0.};
    t_spatial_coords(i) = t_coords(i) + t_disp(i);
 
    auto sdf =
        surface_distance_function(getTStime(), t_spatial_coords, t_traction);
    auto t_grad_sdf = grad_surface_distance_function(
        getTStime(), t_spatial_coords, t_traction);
    auto tn = -t_traction(i) * t_grad_sdf(i);
    auto c = constrain(sdf, tn);
 
    FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM> t_cP;
    t_cP(i, j) = (c * t_grad_sdf(i)) * t_grad_sdf(j);
    FTensor::Tensor2<double, SPACE_DIM, SPACE_DIM> t_cQ;
    t_cQ(i, j) = kronecker_delta(i, j) - t_cP(i, j);
 
    FTensor::Tensor1<double, SPACE_DIM> t_rhs;
    t_rhs(i) =
 
        t_cQ(i, j) * (t_disp(j) - cn * t_traction(j))
 
        +
 
        t_cP(i, j) * t_disp(j) +
        c * (sdf * t_grad_sdf(i)); // add gap0 displacements
 
    size_t bb = 0;
    for (; bb != AssemblyBoundaryEleOp::nbRows / SPACE_DIM; ++bb) {
      const double beta = alpha * (t_base(i) * t_normal(i));
      t_nf(i) -= beta * t_rhs(i);
 
      ++t_nf;
      ++t_base;
    }
    for (; bb < nb_base_functions; ++bb)
      ++t_base;
 
    ++t_disp;
    ++t_traction;
    ++t_coords;
    ++t_w;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> ContactOps::OpConstrainBoundaryRhs::commonDataPtr

private

Examples

ContactOps.hpp.

Definition at line 177 of file ContactOps.hpp.

---

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

• users_modules/tutorials/adv-1/src/ContactOps.hpp