#include "users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp"
Definition at line 296 of file EshelbianOperators.hpp.
◆ OpCalculateTractionFromSideEl()
| OpCalculateTractionFromSideEl::OpCalculateTractionFromSideEl |
( |
boost::shared_ptr< DataAtIntegrationPts > |
data_ptr | ) |
|
|
inline |
Definition at line 297 of file EshelbianOperators.hpp.
300 std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
301 for (
auto t = moab::CN::TypeDimensionMap[
SPACE_DIM].first;
302 t <= moab::CN::TypeDimensionMap[
SPACE_DIM].second; ++
t)
303 doEntities[
t] =
true;
304 sYmm = false;
305 }
constexpr int SPACE_DIM
[Define dimension]
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 850 of file EshelbianOperators.cpp.
851 {
855
856 auto N_InLoop = getNinTheLoop();
857 auto sensee = getSkeletonSense();
858 auto nb_gauss_pts = getGaussPts().size2();
859 auto t_normal = getFTensor1NormalsAtGaussPts();
860
861 auto t_sigma_ptr =
862 getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(
dataAtPts->approxPAtPts);
863 if (N_InLoop == 0) {
866 }
867 auto t_traction = getFTensor1FromMat<SPACE_DIM>(
dataAtPts->tractionAtPts);
868
869 for (int gg = 0; gg != nb_gauss_pts; gg++) {
870 t_traction(
i) = t_sigma_ptr(
i,
j) * sensee * (t_normal(
j) / t_normal.l2());
871 ++t_traction;
872 ++t_sigma_ptr;
873 ++t_normal;
874 }
875
877}
#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: