v0.14.0
Public Member Functions | Private Attributes | List of all members
OpCalculateTractionFromSideEl Struct Reference

#include <users_modules/eshelbian_plasticit/src/EshelbianOperators.hpp>

Inheritance diagram for OpCalculateTractionFromSideEl:
[legend]
Collaboration diagram for OpCalculateTractionFromSideEl:
[legend]

Public Member Functions

 OpCalculateTractionFromSideEl (boost::shared_ptr< DataAtIntegrationPts > data_ptr)
 
MoFEMErrorCode doWork (int side, EntityType type, EntData &data)
 

Private Attributes

boost::shared_ptr< DataAtIntegrationPtsdataAtPts
 data at integration pts More...
 

Detailed Description

Definition at line 187 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

◆ OpCalculateTractionFromSideEl()

OpCalculateTractionFromSideEl::OpCalculateTractionFromSideEl ( boost::shared_ptr< DataAtIntegrationPts data_ptr)
inline

Definition at line 188 of file EshelbianOperators.hpp.

190  : SideEleOp(NOSPACE, SideEleOp::OPSPACE), dataAtPts(data_ptr) {
191  std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
192  for (auto t = moab::CN::TypeDimensionMap[SPACE_DIM].first;
193  t <= moab::CN::TypeDimensionMap[SPACE_DIM].second; ++t)
194  doEntities[t] = true;
195  sYmm = false;
196  }

Member Function Documentation

◆ doWork()

MoFEMErrorCode OpCalculateTractionFromSideEl::doWork ( int  side,
EntityType  type,
EntData data 
)
Examples
EshelbianOperators.cpp.

Definition at line 793 of file EshelbianOperators.cpp.

794  {
798 
799  auto n_total_side_eles = getLoopSize();
800  auto N_InLoop = getNinTheLoop();
801  auto sensee = getSkeletonSense();
802  auto nb_gauss_pts = getGaussPts().size2();
803  auto t_normal = getFTensor1NormalsAtGaussPts();
804 
805  auto t_sigma_ptr =
806  getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(dataAtPts->approxPAtPts);
807  if (N_InLoop == 0) {
808  dataAtPts->tractionAtPts.resize(SPACE_DIM, nb_gauss_pts, false);
809  dataAtPts->tractionAtPts.clear();
810  }
811  auto t_traction = getFTensor1FromMat<SPACE_DIM>(dataAtPts->tractionAtPts);
812 
813  for (int gg = 0; gg != nb_gauss_pts; gg++) {
814  t_traction(i) = t_sigma_ptr(i, j) * sensee * (t_normal(j) / t_normal.l2());
815  ++t_traction;
816  ++t_sigma_ptr;
817  ++t_normal;
818  }
819 
821 }

Member Data Documentation

◆ dataAtPts

boost::shared_ptr<DataAtIntegrationPts> OpCalculateTractionFromSideEl::dataAtPts
private

data at integration pts

Definition at line 201 of file EshelbianOperators.hpp.


The documentation for this struct was generated from the following files:
NOSPACE
@ NOSPACE
Definition: definitions.h:83
FTENSOR_INDEX
#define FTENSOR_INDEX(DIM, I)
Definition: Templates.hpp:2013
SideEleOp
EleOnSide::UserDataOperator SideEleOp
Definition: test_broken_space.cpp:42
SPACE_DIM
constexpr int SPACE_DIM
Definition: child_and_parent.cpp:16
t
constexpr double t
plate stiffness
Definition: plate.cpp:58
OpCalculateTractionFromSideEl::dataAtPts
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
data at integration pts
Definition: EshelbianOperators.hpp:201
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:429
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:359