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

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

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

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

Public Member Functions
	OpAssembleTotalContactAreaImpl (boost::shared_ptr< CommonData > common_data_ptr, bool is_axisymmetric=false, boost::shared_ptr< Range > contact_range_ptr=nullptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Public Attributes
SurfaceDistanceFunction	surfaceDistanceFunction = surface_distance_function
GradSurfaceDistanceFunction	gradSurfaceDistanceFunction

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr
bool	isAxisymmetric
boost::shared_ptr< Range >	contactRange

Detailed Description

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

Definition at line 505 of file ContactOps.hpp.

Constructor & Destructor Documentation

OpAssembleTotalContactAreaImpl()

template<int DIM, typename BoundaryEleOp >

ContactOps::OpAssembleTotalContactAreaImpl< DIM, GAUSS, BoundaryEleOp >::OpAssembleTotalContactAreaImpl	(	boost::shared_ptr< CommonData >	common_data_ptr,
		bool	is_axisymmetric = false,
		boost::shared_ptr< Range >	contact_range_ptr = nullptr
	)

Definition at line 697 of file ContactOps.hpp.

    : BoundaryEleOp(NOSPACE, BoundaryEleOp::OPSPACE),
      commonDataPtr(common_data_ptr), isAxisymmetric(is_axisymmetric),
      contactRange(contact_range_ptr) {}

Member Function Documentation

doWork()

template<int DIM, typename BoundaryEleOp >

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

Definition at line 707 of file ContactOps.hpp.

                                              {
  MoFEMFunctionBegin;
 
  auto fe_type = BoundaryEleOp::getFEType();
 
  const auto fe_ent = BoundaryEleOp::getFEEntityHandle();
 
  if (contactRange->find(fe_ent) != contactRange->end()) {
    FTensor::Index<'i', DIM> i;
    FTensor::Index<'j', DIM> j;
    FTensor::Tensor1<double, 2> t_sum_a{0., 0.};
 
    auto t_w = BoundaryEleOp::getFTensor0IntegrationWeight();
    auto t_traction = getFTensor1FromMat<DIM>(commonDataPtr->contactTraction);
    auto t_coords = BoundaryEleOp::getFTensor1CoordsAtGaussPts();
 
    auto t_grad = getFTensor2FromMat<DIM, DIM>(commonDataPtr->contactDispGrad);
    auto t_normal_at_pts = BoundaryEleOp::getFTensor1NormalsAtGaussPts();
 
    const auto nb_gauss_pts = BoundaryEleOp::getGaussPts().size2();
    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);
 
    auto ts_time = BoundaryEleOp::getTStime();
    auto ts_time_step = BoundaryEleOp::getTStimeStep();
    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 t_sdf = getFTensor0FromVec(v_sdf);
    auto t_grad_sdf = getFTensor1FromMat<3>(m_grad_sdf);
    for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
      double jacobian = 1.;
      if (isAxisymmetric) {
        jacobian = 2. * M_PI * t_coords(0); // Axisymmetric Jacobian
      }
      auto tn = -t_traction(i) * t_grad_sdf(i);
      auto c = constrain(t_sdf, tn);
      double alpha = t_w * jacobian;
 
      FTensor::Tensor2<double, DIM, DIM> F;
      FTensor::Tensor2<double, DIM, DIM> invF;
      FTensor::Tensor1<double, DIM> t_normal_current;
 
      F(i, j) = t_grad(i, j) + kronecker_delta(i, j);
      auto det = determinantTensor(F);
      CHKERR invertTensor(F, det, invF);
      t_normal_current(i) = det * (invF(j, i) * t_normal_at_pts(j));
 
      alpha *= sqrt(t_normal_current(i) * t_normal_current(i));
 
      if (fe_type == MBTRI) {
        alpha /= 2;
      }
      if (c > 1e-12) {
        t_sum_a(0) += alpha; // real area
      }
      t_sum_a(1) += alpha; // Potential area
      ++t_w;
      ++t_traction;
      ++t_coords;
      ++t_sdf;
      ++t_grad_sdf;
 
      ++t_grad;
      ++t_normal_at_pts;
    }
    constexpr int ind[] = {AREA_ACTIVE, AREA_POTENTIAL};
    CHKERR VecSetValues(commonDataPtr->totalTractionAndArea, 2, ind, &t_sum_a(0),
                        ADD_VALUES);
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

template<int DIM, typename BoundaryEleOp >

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

private

Definition at line 518 of file ContactOps.hpp.

contactRange

template<int DIM, typename BoundaryEleOp >

boost::shared_ptr<Range> ContactOps::OpAssembleTotalContactAreaImpl< DIM, GAUSS, BoundaryEleOp >::contactRange

private

Definition at line 520 of file ContactOps.hpp.

gradSurfaceDistanceFunction

template<int DIM, typename BoundaryEleOp >

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

Initial value:

=
      grad_surface_distance_function

Definition at line 514 of file ContactOps.hpp.

isAxisymmetric

template<int DIM, typename BoundaryEleOp >

bool ContactOps::OpAssembleTotalContactAreaImpl< DIM, GAUSS, BoundaryEleOp >::isAxisymmetric

private

Definition at line 519 of file ContactOps.hpp.

surfaceDistanceFunction

template<int DIM, typename BoundaryEleOp >

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

Definition at line 513 of file ContactOps.hpp.

---

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

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