23 {
25
27
29
30 const size_t nb_integration_pts = row_data.getN().size1();
31 const size_t nb_row_base_functions = row_data.getN().size2();
32 auto t_w = getFTensor0IntegrationWeight();
33
34 auto t_row_base = row_data.getFTensor0N();
35
38 auto t_flow =
39 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->plasticFlow);
40 auto t_D =
41 getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*
commonDataPtr->mDPtr);
42
43 auto t_D_Op = getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*
mDPtr);
44
45 auto t_grad =
46 getFTensor2FromMat<SPACE_DIM, SPACE_DIM>(*(
commonDataPtr->mGradPtr));
47 auto t_logC_dC = getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM>(
49
52
53 for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
54
55 double alpha = getMeasure() * t_w;
56
61
63 t_D_Op,
66 t_flow_stress_dlogC(
i,
j,
k,
l) =
67 t_D(
i,
j,
m,
n) * t_diff_plastic_flow_dstrain(
m,
n,
k,
l);
68
70 t_diff_plastic_flow_stress_diff_symm;
71 t_diff_plastic_flow_stress_diff_symm(
i,
j,
k,
l) =
72 t_flow_stress_dlogC(
i,
j,
m,
n) * t_diff_symmetrize(
m,
n,
k,
l);
73
75 t_diff_plastic_flow_stress_dC;
76 t_diff_plastic_flow_stress_dC(
i,
j,
k,
l) =
77 t_flow_stress_dlogC(
i,
j,
m,
n) * (t_logC_dC(
m,
n,
k,
l) / 2);
78
80 t_diff_plastic_flow_stress_dgrad;
81 t_diff_plastic_flow_stress_dgrad(
i,
j,
k,
l) = 0;
88 t_diff_plastic_flow_stress_dgrad(ii, jj, kk, ll) +=
89 t_diff_plastic_flow_stress_dC(ii, jj, mm, nn) *
90 dC_dF(mm, nn, kk, ll);
91
92
93
94 size_t rr = 0;
96
97 auto t_mat =
99
100 const double c0 =
alpha * t_row_base * t_tau_dot;
101
102 auto t_col_diff_base = col_data.getFTensor1DiffN<
SPACE_DIM>(gg, 0);
104
106 t_mat(
i,
j,
l) -= c0 * (t_diff_plastic_flow_stress_dgrad(
i,
j,
l,
k) *
108
109 ++t_mat;
110 ++t_col_diff_base;
111 }
112
113 ++t_row_base;
114 }
115 for (; rr < nb_row_base_functions; ++rr)
116 ++t_row_base;
117
118 ++t_w;
119 ++t_f;
120 ++t_flow;
121 ++t_tau_dot;
122
123 ++t_logC_dC;
124 ++t_grad;
125 }
126
128}
#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.
FTensor::Index< 'j', SPACE_DIM > j
[Common data]
FTensor::Index< 'l', SPACE_DIM > l
FTensor::Index< 'k', SPACE_DIM > k
auto diff_tensor()
[Operators definitions]
auto diff_deviator(FTensor::Ddg< double, SPACE_DIM, SPACE_DIM > &&t_diff_stress)
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'm', SPACE_DIM > m
static FTensor::Tensor3< FTensor::PackPtr< double *, 2 >, 2, 2, 2 > get_mat_tensor_sym_dvector(size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
[Lambda functions]
auto diff_plastic_flow_dstress(long double f, FTensor::Tensor2_symmetric< T, SPACE_DIM > &t_flow, FTensor::Ddg< double, 3, SPACE_DIM > &&t_diff_deviator)
auto diff_plastic_flow_dstrain(FTensor::Ddg< T, SPACE_DIM, SPACE_DIM > &t_D, FTensor::Ddg< double, SPACE_DIM, SPACE_DIM > &&t_diff_plastic_flow_dstress)
FTensor::Index< 'n', SPACE_DIM > n
int nbRows
number of dofs on rows
MatrixDouble locMat
local entity block matrix
int nbCols
number if dof on column
Public Member Functions inherited from MoFEM::OpBaseImpl< A, EleOp >