102 {
104
105 auto t_flow =
106 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->plasticFlow);
107
108 auto t_plastic_strain_dot =
109 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->plasticStrainDot);
111
112 auto t_D =
113 getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*
commonDataPtr->mDPtr);
115
116 const size_t nb_integration_pts = data.getN().size1();
117 const size_t nb_base_functions = data.getN().size2();
118 auto t_w = getFTensor0IntegrationWeight();
119 auto t_base = data.getFTensor0N();
120
122
123 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
124 double alpha = getMeasure() * t_w;
125
127
128 const double tau_dot = t_tau_dot;
129
131 t_flow_stress_diff(
i,
j) = t_D(
i,
j,
k,
l) * (t_plastic_strain_dot(
k,
l) -
132 t_tau_dot * t_flow(
k,
l));
133
134 size_t bb = 0;
136 t_nf(
i,
j) += (
alpha * t_base) * t_flow_stress_diff(
i,
j);
137 ++t_base;
138 ++t_nf;
139 }
140 for (; bb < nb_base_functions; ++bb)
141 ++t_base;
142
143 ++t_flow;
144 ++t_plastic_strain_dot;
145 ++t_tau_dot;
146 ++t_w;
147 }
148
150}
#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()
static auto getFTensor0FromVec(ublas::vector< T, A > &data)
Get tensor rank 0 (scalar) form data vector.
static FTensor::Tensor2< FTensor::PackPtr< double *, 3 >, 2, 2 > get_nf(VectorDouble &nf, FTensor::Number< 2 >)
FTensor::Index< 'j', SPACE_DIM > j
[Common data]
FTensor::Index< 'l', SPACE_DIM > l
FTensor::Index< 'k', SPACE_DIM > k
FTensor::Index< 'i', SPACE_DIM > i
VectorDouble locF
local entity vector
int nbRows
number of dofs on rows
Public Member Functions inherited from MoFEM::OpBaseImpl< A, EleOp >