1018 {
1019 int nb_points_avtive_on_elem = 0;
1020 int nb_points_on_elem = 0;
1021
1023 auto t_tau_dot = getFTensor0FromVec(
commonDataPtr->plasticTauDot);
1024 auto t_f = getFTensor0FromVec(
commonDataPtr->plasticSurface);
1025 auto t_plastic_strain_dot =
1026 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->plasticStrainDot);
1027 auto t_temp = getFTensor0FromVec(TPCommonDataPtr->temperature);
1028
1029 auto dt = this->getTStimeStep();
1030
1031 for (auto gg = 0; gg != nb_gauss_pts; ++gg) {
1033 const auto ww =
1034 w(eqiv, t_tau_dot, t_f,
1038 TPCommonDataPtr->temp0, t_temp),
1041
1042 ++nb_points_on_elem;
1043 if (sign_ww > 0) {
1044 ++nb_points_avtive_on_elem;
1045 }
1046
1047 ++t_tau;
1048 ++t_tau_dot;
1049 ++t_f;
1050 ++t_plastic_strain_dot;
1051 ++t_temp;
1052 }
1053
1059
1060 ++nb_elements;
1061 nb_points += nb_points_on_elem;
double constrian_sign(double x, double dt)
double w(double eqiv, double dot_tau, double f, double sigma_y, double sigma_Y)
auto equivalent_strain_dot(FTensor::Tensor2_symmetric< T, DIM > &t_plastic_strain_dot)
static std::array< int, 5 > activityData
boost::shared_ptr< CommonData > commonDataPtr
double iso_hardening(double tau, double H, double Qinf, double b_iso, double sigmaY)