ContactOps::OpConstrainBoundaryLhs_dU Struct Reference

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

Inheritance diagram for ContactOps::OpConstrainBoundaryLhs_dU:

Collaboration diagram for ContactOps::OpConstrainBoundaryLhs_dU:

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

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

Definition at line 180 of file ContactOps.hpp.

Constructor & Destructor Documentation

OpConstrainBoundaryLhs_dU()

ContactOps::OpConstrainBoundaryLhs_dU::OpConstrainBoundaryLhs_dU	(	const std::string	row_field_name,
		const std::string	col_field_name,
		boost::shared_ptr< CommonData >	common_data_ptr
	)

Examples

ContactOps.hpp.

Definition at line 347 of file ContactOps.hpp.

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

Member Function Documentation

iNtegrate()

MoFEMErrorCode ContactOps::OpConstrainBoundaryLhs_dU::iNtegrate	(	EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

Examples

ContactOps.hpp.

Definition at line 357 of file ContactOps.hpp.

                                                                         {
  MoFEMFunctionBegin;
 
  const size_t nb_gauss_pts = getGaussPts().size2();
  auto &locMat = AssemblyBoundaryEleOp::locMat;
 
  auto t_normal = getFTensor1Normal();
  t_normal(i) /= sqrt(t_normal(j) * t_normal(j));
 
  auto t_disp = getFTensor1FromMat<SPACE_DIM>(commonDataPtr->contactDisp);
  auto t_traction =
      getFTensor1FromMat<SPACE_DIM>(commonDataPtr->contactTraction);
  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;
 
  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
 
  for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
 
    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::Tensor2<double, 3, 3> t_res_dU;
    t_res_dU(i, j) =
        kronecker_delta(i, j) + t_cP(i, j);
 
    if(c>0) {
      auto t_hess_sdf = hess_surface_distance_function(
          getTStime(), t_spatial_coords, t_traction);
      t_res_dU(i, j) +=
          (c * cn) * (t_hess_sdf(i, j) * (t_grad_sdf(k) * t_traction(k)) +
                      t_grad_sdf(i) * t_hess_sdf(k, j) * t_traction(k)) +
          c * sdf * t_hess_sdf(i, j);
    }
 
    size_t rr = 0;
    for (; rr != AssemblyBoundaryEleOp::nbRows / SPACE_DIM; ++rr) {
 
      auto t_mat = getFTensor2FromArray<SPACE_DIM, SPACE_DIM, SPACE_DIM>(
          locMat, SPACE_DIM * rr);
 
      const double row_base = t_row_base(i) * t_normal(i);
 
      auto t_col_base = col_data.getFTensor0N(gg, 0);
      for (size_t cc = 0; cc != AssemblyBoundaryEleOp::nbCols / SPACE_DIM;
           ++cc) {
        const double beta = alpha * row_base * t_col_base;
 
        t_mat(i, j) -= beta * t_res_dU(i, j);
 
        ++t_col_base;
        ++t_mat;
      }
 
      ++t_row_base;
    }
    for (; rr < nb_face_functions; ++rr)
      ++t_row_base;
 
    ++t_disp;
    ++t_traction;
    ++t_coords;
    ++t_w;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

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

private

Examples

ContactOps.hpp.

Definition at line 188 of file ContactOps.hpp.

---

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

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