ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp > Struct Template Reference

#include "tutorials/adv-1/src/ContactOps.hpp"

Inheritance diagram for ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >:

Collaboration diagram for ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >:

Public Member Functions
	OpEvaluateSDFImpl (boost::shared_ptr< CommonData > common_data_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr
SurfaceDistanceFunction	surfaceDistanceFunction = surface_distance_function
GradSurfaceDistanceFunction	gradSurfaceDistanceFunction
HessSurfaceDistanceFunction	hessSurfaceDistanceFunction

Detailed Description

template<int DIM, typename BoundaryEleOp>
struct ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >

Definition at line 479 of file ContactOps.hpp.

Constructor & Destructor Documentation

OpEvaluateSDFImpl()

template<int DIM, typename BoundaryEleOp >

ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >::OpEvaluateSDFImpl

(

boost::shared_ptr< CommonData >

common_data_ptr

)

Definition at line 745 of file ContactOps.hpp.

    : BoundaryEleOp(NOSPACE, BoundaryEleOp::OPSPACE),
      commonDataPtr(common_data_ptr) {}

Member Function Documentation

doWork()

template<int DIM, typename BoundaryEleOp >

MoFEMErrorCode ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Definition at line 752 of file ContactOps.hpp.

                                                                    {
  MoFEMFunctionBegin;
 
  const auto nb_gauss_pts = BoundaryEleOp::getGaussPts().size2();
  auto &sdf_vec = commonDataPtr->sdfVals;
  auto &grad_mat = commonDataPtr->gradsSdf;
  auto &hess_mat = commonDataPtr->hessSdf;
  auto &constraint_vec = commonDataPtr->constraintVals;
  auto &contactTraction_mat = commonDataPtr->contactTraction;
 
  sdf_vec.resize(nb_gauss_pts, false);
  grad_mat.resize(DIM, nb_gauss_pts, false);
  hess_mat.resize((DIM * (DIM + 1)) / 2, nb_gauss_pts, false);
  constraint_vec.resize(nb_gauss_pts, false);
 
  auto t_traction = getFTensor1FromMat<DIM>(contactTraction_mat);
 
  auto t_sdf = getFTensor0FromVec(sdf_vec);
  auto t_grad_sdf = getFTensor1FromMat<DIM>(grad_mat);
  auto t_hess_sdf = getFTensor2SymmetricFromMat<DIM>(hess_mat);
  auto t_constraint = getFTensor0FromVec(constraint_vec);
 
  auto t_disp = getFTensor1FromMat<DIM>(commonDataPtr->contactDisp);
  auto t_coords = BoundaryEleOp::getFTensor1CoordsAtGaussPts();
  auto t_normal_at_pts = BoundaryEleOp::getFTensor1NormalsAtGaussPts();
 
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
 
  auto ts_time = BoundaryEleOp::getTStime();
  auto ts_time_step = BoundaryEleOp::getTStimeStep();
 
  auto m_spatial_coords = get_spatial_coords(
      BoundaryEleOp::getFTensor1CoordsAtGaussPts(),
      getFTensor1FromMat<DIM>(commonDataPtr->contactDisp), nb_gauss_pts);
  auto m_normals_at_pts = get_normalize_normals(
      BoundaryEleOp::getFTensor1NormalsAtGaussPts(), nb_gauss_pts);
 
  // placeholder to pass boundary block id to python
  int block_id = 0;
 
  auto v_sdf =
      surfaceDistanceFunction(ts_time_step, ts_time, nb_gauss_pts,
                              m_spatial_coords, m_normals_at_pts, block_id);
 
  auto m_grad_sdf =
      gradSurfaceDistanceFunction(ts_time_step, ts_time, nb_gauss_pts,
                                  m_spatial_coords, m_normals_at_pts, block_id);
 
  auto m_hess_sdf =
      hessSurfaceDistanceFunction(ts_time_step, ts_time, nb_gauss_pts,
                                  m_spatial_coords, m_normals_at_pts, block_id);
 
  auto t_sdf_v = getFTensor0FromVec(v_sdf);
  auto t_grad_sdf_v = getFTensor1FromMat<3>(m_grad_sdf);
  auto t_hess_sdf_v = getFTensor2SymmetricFromMat<3>(m_hess_sdf);
 
  auto next = [&]() {
    ++t_sdf;
    ++t_sdf_v;
    ++t_grad_sdf;
    ++t_grad_sdf_v;
    ++t_hess_sdf;
    ++t_hess_sdf_v;
    ++t_disp;
    ++t_traction;
    ++t_constraint;
  };
 
  for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
 
    auto tn = -t_traction(i) * t_grad_sdf_v(i);
    auto c = constrain(t_sdf_v, tn);
 
    t_sdf = t_sdf_v;
    t_grad_sdf(i) = t_grad_sdf_v(i);
    t_hess_sdf(i, j) = t_hess_sdf_v(i, j);
    t_constraint = c;
 
    next();
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

template<int DIM, typename BoundaryEleOp >

boost::shared_ptr<CommonData> ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >::commonDataPtr

private

Definition at line 484 of file ContactOps.hpp.

gradSurfaceDistanceFunction

template<int DIM, typename BoundaryEleOp >

GradSurfaceDistanceFunction ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >::gradSurfaceDistanceFunction

private

Initial value:

=
      grad_surface_distance_function

Definition at line 487 of file ContactOps.hpp.

hessSurfaceDistanceFunction

template<int DIM, typename BoundaryEleOp >

HessSurfaceDistanceFunction ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >::hessSurfaceDistanceFunction

private

Initial value:

=
      hess_surface_distance_function

Definition at line 489 of file ContactOps.hpp.

surfaceDistanceFunction

template<int DIM, typename BoundaryEleOp >

SurfaceDistanceFunction ContactOps::OpEvaluateSDFImpl< DIM, GAUSS, BoundaryEleOp >::surfaceDistanceFunction = surface_distance_function

private

Definition at line 486 of file ContactOps.hpp.

---

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

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