v0.16.0
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< DataAtIntegrationPtsdataAtPts
 data at integration pts
 

Detailed Description

Definition at line 299 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

◆ OpCalculateTractionFromSideEle()

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

Definition at line 300 of file EshelbianOperators.hpp.

302 : SideEleOp(NOSPACE, SideEleOp::OPSPACE), dataAtPts(data_ptr) {
303 std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
304 for (auto t = moab::CN::TypeDimensionMap[SPACE_DIM].first;
305 t <= moab::CN::TypeDimensionMap[SPACE_DIM].second; ++t)
306 doEntities[t] = true;
307 sYmm = false;
308 }
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 944 of file EshelbianOperators.cpp.

945 {
948
949 auto n_in_the_loop = getNinTheLoop();
950 auto loop_size = getLoopSize();
951 auto sense = getSkeletonSense();
952 auto nb_gauss_pts = getGaussPts().size2();
953 auto t_normal = getFTensor1NormalsAtGaussPts();
954
955 auto t_sigma = dataAtPts->getFTensorApproxP(getGaussPts().size2());
956 auto get_tracion =
958 dataAtPts->tractionAtPts, nb_gauss_pts);
959 if (!n_in_the_loop) {
960 dataAtPts->tractionAtPts.clear();
961 }
962
963 auto t_traction = get_tracion();
964 for (int gg = 0; gg != nb_gauss_pts; gg++) {
965 t_traction(i) +=
966 t_sigma(i, j) * sense * (t_normal(j) / t_normal.l2()) / loop_size;
967 ++t_traction;
968 ++t_sigma;
969 ++t_normal;
970 }
971
973}
#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 313 of file EshelbianOperators.hpp.


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