v0.16.3
Loading...
Searching...
No Matches
Public Member Functions | Private Attributes | List of all members
OpCalculateTractionFromSideEle Struct Reference

#include "users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp"

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

Public Member Functions

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

Private Attributes

boost::shared_ptr< DataAtIntegrationPts > dataAtPts
 data at integration pts
 

Detailed Description

Definition at line 285 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

◆ OpCalculateTractionFromSideEle()

OpCalculateTractionFromSideEle::OpCalculateTractionFromSideEle ( boost::shared_ptr< DataAtIntegrationPts >  data_ptr)
inline

Definition at line 286 of file EshelbianOperators.hpp.

288 : SideEleOp(NOSPACE, SideEleOp::OPSPACE), dataAtPts(data_ptr) {
289 std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
290 for (auto t = moab::CN::TypeDimensionMap[SPACE_DIM].first;
291 t <= moab::CN::TypeDimensionMap[SPACE_DIM].second; ++t)
292 doEntities[t] = true;
293 sYmm = false;
294 }
constexpr int SPACE_DIM
@ NOSPACE
Definition definitions.h:83
constexpr double t
plate stiffness
Definition plate.cpp:58
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
data at integration pts
EleOnSide::UserDataOperator SideEleOp

Member Function Documentation

◆ doWork()

MoFEMErrorCode OpCalculateTractionFromSideEle::doWork ( int  side,
EntityType  type,
EntData &  data 
)
Examples
/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 996 of file EshelbianOperators.cpp.

997 {
1000
1001 auto n_in_the_loop = getNinTheLoop();
1002 auto loop_size = getLoopSize();
1003 auto sense = getSkeletonSense();
1004 auto nb_gauss_pts = getGaussPts().size2();
1005 auto t_normal = getFTensor1NormalsAtGaussPts();
1006
1007 auto t_sigma = dataAtPts->getFTensorApproxP(getGaussPts().size2());
1008 auto get_tracion =
1010 dataAtPts->tractionAtPts, nb_gauss_pts);
1011 if (!n_in_the_loop) {
1012 dataAtPts->tractionAtPts.clear();
1013 }
1014
1015 auto t_traction = get_tracion();
1016 for (int gg = 0; gg != nb_gauss_pts; gg++) {
1017 t_traction(i) +=
1018 t_sigma(i, j) * sense * (t_normal(j) / t_normal.l2()) / loop_size;
1019 ++t_traction;
1020 ++t_sigma;
1021 ++t_normal;
1022 }
1023
1025}
#define FTENSOR_INDEXES(DIM,...)
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
DataLayoutTraits< DataLayout::GaussByCoeffs > DL
Definition MatHuHu.hpp:33
decltype(GetFTensor1FromMatImpl< Tensor_Dim, S, DL, M >::get(std::declval< M & >(), 0, 0)) GetFTensor1FromMatType

Member Data Documentation

◆ dataAtPts

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

data at integration pts

Definition at line 299 of file EshelbianOperators.hpp.


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