468 const size_t nb_integration_pts = row_data.getN().size1();
469 const size_t nb_row_base_functions = row_data.getN().size2();
470 auto t_w = getFTensor0IntegrationWeight();
472 auto t_row_base = row_data.getFTensor0N();
477 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->plasticFlow);
480 getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*
commonDataPtr->mDPtr);
482 auto t_D_Op = getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*
mDPtr);
486 for (
size_t gg = 0; gg != nb_integration_pts; ++gg) {
487 double alpha = getMeasure() * t_w;
488 double beta =
alpha * getTSa();
493 const double c0 =
alpha * t_row_base * t_tau_dot;
494 const double c1 = beta * t_row_base;
503 auto t_col_base = col_data.getFTensor0N(gg, 0);
508 (t_D(
i,
j,
m,
n) * (c1 * t_diff_plastic_strain(
m,
n,
k,
l) +
509 c0 * t_diff_plastic_flow_dstrain(
m,
n,
k,
l)));
517 for (; rr < nb_row_base_functions; ++rr)
#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 auto get_mat_tensor_sym_dtensor_sym(size_t rr, MatrixDouble &mat, FTensor::Number< 2 >)
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
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