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

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

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

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

Public Member Functions
	OpAssembleTotalContactTractionImpl (boost::shared_ptr< CommonData > common_data_ptr, double scale=1, bool is_axisymmetric=false)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr
const double	scaleTraction
bool	isAxisymmetric

Detailed Description

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

Definition at line 446 of file ContactOps.hpp.

Constructor & Destructor Documentation

OpAssembleTotalContactTractionImpl()

template<int DIM, typename BoundaryEleOp >

ContactOps::OpAssembleTotalContactTractionImpl< DIM, GAUSS, BoundaryEleOp >::OpAssembleTotalContactTractionImpl	(	boost::shared_ptr< CommonData >	common_data_ptr,
		double	scale = 1,
		bool	is_axisymmetric = false
	)

Definition at line 609 of file ContactOps.hpp.

    : BoundaryEleOp(NOSPACE, BoundaryEleOp::OPSPACE),
      commonDataPtr(common_data_ptr), scaleTraction(scale),
      isAxisymmetric(is_axisymmetric) {}

Member Function Documentation

doWork()

template<int DIM, typename BoundaryEleOp >

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

Definition at line 619 of file ContactOps.hpp.

                                              {
  MoFEMFunctionBegin;
 
  FTensor::Index<'i', DIM> i;
  FTensor::Tensor1<double, 3> t_sum_t{0., 0., 0.};
 
  auto t_w = BoundaryEleOp::getFTensor0IntegrationWeight();
  auto t_traction = getFTensor1FromMat<DIM>(commonDataPtr->contactTraction);
  auto t_coords = BoundaryEleOp::getFTensor1CoordsAtGaussPts();
 
  const auto nb_gauss_pts = BoundaryEleOp::getGaussPts().size2();
  for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
    double jacobian = 1.;
    if (isAxisymmetric) {
      jacobian = 2. * M_PI * t_coords(0);
    }
    const double alpha = t_w * jacobian * BoundaryEleOp::getMeasure();
    t_sum_t(i) += alpha * t_traction(i);
    ++t_w;
    ++t_traction;
    ++t_coords;
  }
 
  t_sum_t(i) *= scaleTraction;
 
  constexpr int ind[] = {0, 1, 2};
  CHKERR VecSetValues(commonDataPtr->totalTraction, 3, ind, &t_sum_t(0),
                      ADD_VALUES);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

template<int DIM, typename BoundaryEleOp >

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

private

Definition at line 454 of file ContactOps.hpp.

isAxisymmetric

template<int DIM, typename BoundaryEleOp >

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

private

Definition at line 456 of file ContactOps.hpp.

scaleTraction

template<int DIM, typename BoundaryEleOp >

const double ContactOps::OpAssembleTotalContactTractionImpl< DIM, GAUSS, BoundaryEleOp >::scaleTraction

private

Definition at line 455 of file ContactOps.hpp.

---

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

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