914 {
918 };
922 int nb_integration_pts = getGaussPts().size2();
923 auto t_w = getFTensor0IntegrationWeight();
924 auto t_traction = getFTensor1FromMat<SPACE_DIM>(
dataAtPts->tractionAtPts);
925 auto t_coords = getFTensor1CoordsAtGaussPts();
926 auto t_spatial_disp = getFTensor1FromMat<SPACE_DIM>(
dataAtPts->wL2AtPts);
927
929
933
934 for (auto gg = 0; gg != nb_integration_pts; ++gg) {
935 loc_reaction_forces(
i) += (t_traction(
i)) * t_w * getMeasure();
936 t_coords_spatial(
i) = t_coords(
i) + t_spatial_disp(
i);
937
938 loc_moment_forces(
i) +=
939 (FTensor::levi_civita<double>(
i,
j,
k) * t_coords_spatial(
j)) *
940 t_traction(
k) * t_w * getMeasure();
941 ++t_coords;
942 ++t_spatial_disp;
943 ++t_w;
944 ++t_traction;
945 }
946
953
955}
#define FTENSOR_INDEX(DIM, I)
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#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
FTensor::Index< 'k', 3 > k