#include "users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp"
Definition at line 300 of file EshelbianOperators.hpp.
◆ OpCalculateTractionFromSideEl()
| OpCalculateTractionFromSideEl::OpCalculateTractionFromSideEl |
( |
boost::shared_ptr< DataAtIntegrationPts > |
data_ptr | ) |
|
|
inline |
Definition at line 301 of file EshelbianOperators.hpp.
304 std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
305 for (
auto t = moab::CN::TypeDimensionMap[
SPACE_DIM].first;
306 t <= moab::CN::TypeDimensionMap[
SPACE_DIM].second; ++
t)
307 doEntities[
t] =
true;
308 sYmm = false;
309 }
constexpr double t
plate stiffness
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
data at integration pts
EleOnSide::UserDataOperator SideEleOp
◆ doWork()
- Examples
- mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.
Definition at line 884 of file EshelbianOperators.cpp.
885 {
889
890 auto N_InLoop = getNinTheLoop();
891 auto sensee = getSkeletonSense();
892 auto nb_gauss_pts = getGaussPts().size2();
893 auto t_normal = getFTensor1NormalsAtGaussPts();
894
895 auto t_sigma_ptr =
896 getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(
dataAtPts->approxPAtPts);
897 if (N_InLoop == 0) {
900 }
901 auto t_traction = getFTensor1FromMat<SPACE_DIM>(
dataAtPts->tractionAtPts);
902
903 for (int gg = 0; gg != nb_gauss_pts; gg++) {
904 t_traction(
i) = t_sigma_ptr(
i,
j) * sensee * (t_normal(
j) / t_normal.l2());
905 ++t_traction;
906 ++t_sigma_ptr;
907 ++t_normal;
908 }
909
911}
#define FTENSOR_INDEX(DIM, I)
#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
◆ dataAtPts
The documentation for this struct was generated from the following files: