OpCalculateTractionJump Struct Reference

#include "tutorials/vec-8_elasticity_adaptive/src/AdaptiveOrderRef.hpp"

Inheritance diagram for OpCalculateTractionJump:

Collaboration diagram for OpCalculateTractionJump:

Public Member Functions
	OpCalculateTractionJump (const std::string field_name, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< MatrixDouble > traction_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< MatrixDouble >	stressPtr
boost::shared_ptr< MatrixDouble >	tractionJumpPtr

Detailed Description

Definition at line 99 of file AdaptiveOrderRef.hpp.

Constructor & Destructor Documentation

OpCalculateTractionJump()

OpCalculateTractionJump::OpCalculateTractionJump ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  stress_ptr, boost::shared_ptr< MatrixDouble >  traction_ptr  )

inline

Definition at line 100 of file AdaptiveOrderRef.hpp.

      : SideEleOp(field_name, SideEleOp::OPROW, false), stressPtr(stress_ptr),
        tractionJumpPtr(traction_ptr) {
    std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
    doEntities[MBVERTEX] = true;
  }

Member Function Documentation

doWork()

MoFEMErrorCode OpCalculateTractionJump::doWork ( int  side, EntityType  type, EntData &  data  )

inline

Definition at line 109 of file AdaptiveOrderRef.hpp.

                                                                  {
    MoFEMFunctionBegin;
    FTensor::Index<'i', SPACE_DIM> i;
    FTensor::Index<'j', SPACE_DIM> j;
    auto nb_gauss_pts = getGaussPts().size2();
 
    // auto n_total_side_eles = getLoopSize();
    auto n_InLoop = getNinTheLoop();
    auto e_sense = getSkeletonSense();
 
    if (n_InLoop == 0) {
      tractionJumpPtr->resize(nb_gauss_pts, SPACE_DIM, false);
      tractionJumpPtr->clear();
    }
 
    auto t_normal = getFTensor1NormalsAtGaussPts();
    auto t_stress = getFTensor2SymmetricFromMat<SPACE_DIM>(*stressPtr);
    auto t_traction = getFTensor1FromMat<SPACE_DIM>(*tractionJumpPtr);
 
    for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
      FTensor::Tensor1<double, 3> t_unit_normal;
      t_unit_normal(i) = t_normal(i);
      t_unit_normal.normalize();
 
      t_traction(i) -= e_sense * (t_stress(i, j) * t_unit_normal(j));
 
      ++t_normal;
      ++t_stress;
      ++t_traction;
    }
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

stressPtr

boost::shared_ptr<MatrixDouble> OpCalculateTractionJump::stressPtr

private

Definition at line 143 of file AdaptiveOrderRef.hpp.

tractionJumpPtr

boost::shared_ptr<MatrixDouble> OpCalculateTractionJump::tractionJumpPtr

private

Definition at line 144 of file AdaptiveOrderRef.hpp.

---

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

• tutorials/vec-8_elasticity_adaptive/src/AdaptiveOrderRef.hpp