https://mooseframework.inl.gov
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
ComputeSimoHughesJ2PlasticityStress Class Referenceabstract

#include <ComputeSimoHughesJ2PlasticityStress.h>

Inheritance diagram for ComputeSimoHughesJ2PlasticityStress:
[legend]

Public Types

typedef DerivativeMaterialPropertyNameInterface::SymbolName SymbolName
 

Public Member Functions

 ComputeSimoHughesJ2PlasticityStress (const InputParameters &parameters)
 
virtual void initialSetup () override
 
const GenericMaterialProperty< U, is_ad > & getDefaultMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< U, is_ad > & getDefaultMaterialPropertyByName (const std::string &name)
 
void validateDerivativeMaterialPropertyBase (const std::string &base)
 
const MaterialPropertyName derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
void validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
 
void validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
 
void validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
const MaterialPropertyName propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 
const MaterialPropertyName propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

virtual void initQpStatefulProperties () override
 
virtual void computeQpPK1Stress () override
 
void returnMappingSolve (const GenericReal< is_ad > &effective_trial_stress, GenericReal< is_ad > &scalar, const ConsoleStream &console)
 Perform the return mapping iterations. More...
 
virtual GenericReal< is_ad > minimumPermissibleValue (const GenericReal< is_ad > &effective_trial_stress) const
 Compute the minimum permissible value of the scalar. More...
 
virtual GenericReal< is_ad > maximumPermissibleValue (const GenericReal< is_ad > &effective_trial_stress) const
 Compute the maximum permissible value of the scalar. More...
 
virtual GenericReal< is_ad > initialGuess (const GenericReal< is_ad > &)
 Compute an initial guess for the value of the scalar. More...
 
virtual GenericReal< is_ad > computeResidual (const GenericReal< is_ad > &, const GenericReal< is_ad > &)=0
 Compute the residual for a predicted value of the scalar. More...
 
virtual GenericReal< is_ad > computeDerivative (const GenericReal< is_ad > &, const GenericReal< is_ad > &)=0
 Compute the derivative of the residual as a function of the scalar variable. More...
 
virtual GenericChainedReal< is_ad > computeResidualAndDerivative (const GenericReal< is_ad > &, const GenericChainedReal< is_ad > &)
 Compute the residual and the derivative for a predicted value of the scalar. More...
 
virtual Real computeReferenceResidual (const GenericReal< is_ad > &effective_trial_stress, const GenericReal< is_ad > &scalar)=0
 Compute a reference quantity to be used for checking relative convergence. More...
 
virtual void preStep (const GenericReal< is_ad > &, const GenericReal< is_ad > &, const GenericReal< is_ad > &)
 This method is called before taking a step in the return mapping algorithm. More...
 
virtual void iterationFinalize (const GenericReal< is_ad > &)
 Finalize internal state variables for a model for a given iteration. More...
 
virtual void outputIterationSummary (std::stringstream *iter_output, const unsigned int total_it)
 Output summary information for the convergence history of the model. More...
 
virtual void outputIterationStep (std::stringstream *iter_output, const GenericReal< is_ad > &effective_trial_stress, const GenericReal< is_ad > &scalar, const Real reference_residual)
 Output information for a single iteration step to build the convergence history of the model. More...
 
bool converged (const GenericReal< is_ad > &residual, const Real reference)
 Check to see whether the residual is within the convergence limits. More...
 
virtual Real computeReferenceResidual (const Real &effective_trial_stress, const Real &scalar) override
 The return mapping residual and derivative. More...
 
virtual Real computeResidual (const Real &effective_trial_stress, const Real &scalar) override
 
virtual Real computeDerivative (const Real &effective_trial_stress, const Real &scalar) override
 
virtual void preStep (const Real &scalar_old, const Real &residual, const Real &jacobian) override
 

Protected Attributes

const MaterialPropertyName _elasticity_tensor_name
 
const MaterialProperty< RankFourTensor > & _elasticity_tensor
 
const MaterialProperty< RankTwoTensor > & _F_old
 
const std::string _ep_name
 
MaterialProperty< Real > & _ep
 
const MaterialProperty< Real > & _ep_old
 
MaterialProperty< RankTwoTensor > & _be
 
const MaterialProperty< RankTwoTensor > & _be_old
 
MaterialProperty< RankTwoTensor > & _Np
 
MaterialBase * _flow_stress_material
 
const std::string _flow_stress_name
 
const MaterialProperty< Real > & _H
 
const MaterialProperty< Real > & _dH
 
const MaterialProperty< Real > & _d2H
 
bool _check_range
 Whether to check to see whether iterative solution is within admissible range, and set within that range if outside. More...
 
bool _line_search
 Whether to use line searches to improve convergence. More...
 
bool _bracket_solution
 Whether to save upper and lower bounds of root for scalar, and set solution to the midpoint between those bounds if outside them. More...
 

Private Attributes

RankFourTensor _d_be_d_F
 Helper (dummy) variables for iteratively updating the consistant tangent during return mapping. More...
 
RankFourTensor _d_n_d_be
 
RankTwoTensor _d_deltaep_d_betr
 
RankTwoTensor _d_R_d_betr
 
RankTwoTensor _d_J_d_betr
 

Detailed Description

Definition at line 20 of file ComputeSimoHughesJ2PlasticityStress.h.

Constructor & Destructor Documentation

◆ ComputeSimoHughesJ2PlasticityStress()

ComputeSimoHughesJ2PlasticityStress::ComputeSimoHughesJ2PlasticityStress ( const InputParameters &  parameters)

Definition at line 27 of file ComputeSimoHughesJ2PlasticityStress.C.

29  : DerivativeMaterialInterface<ComputeLagrangianStressPK1>(parameters),
30  SingleVariableReturnMappingSolution(parameters),
31  _elasticity_tensor_name(_base_name + getParam<MaterialPropertyName>("elasticity_tensor")),
32  _elasticity_tensor(getMaterialProperty<RankFourTensor>(_elasticity_tensor_name)),
33  _F_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "deformation_gradient")),
34  _ep_name(_base_name + "effective_plastic_strain"),
35  _ep(declareProperty<Real>(_ep_name)),
36  _ep_old(getMaterialPropertyOldByName<Real>(_ep_name)),
37  _be(declareProperty<RankTwoTensor>(_base_name +
38  "volume_preserving_elastic_left_cauchy_green_strain")),
39  _be_old(getMaterialPropertyOldByName<RankTwoTensor>(
40  _base_name + "volume_preserving_elastic_left_cauchy_green_strain")),
41  _Np(declareProperty<RankTwoTensor>(_base_name + "flow_direction")),
42  _flow_stress_material(nullptr),
43  _flow_stress_name(_base_name + "flow_stress"),
44  _H(getMaterialPropertyByName<Real>(_flow_stress_name)),
45  _dH(getDefaultMaterialPropertyByName<Real, false>(
46  derivativePropertyName(_flow_stress_name, {_ep_name}))),
47  _d2H(getDefaultMaterialPropertyByName<Real, false>(
48  derivativePropertyName(_flow_stress_name, {_ep_name, _ep_name})))
49 {
50 }
ComputeSimoHughesJ2PlasticityStress::_elasticity_tensor
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Definition: ComputeSimoHughesJ2PlasticityStress.h:46
ComputeSimoHughesJ2PlasticityStress::_be_old
const MaterialProperty< RankTwoTensor > & _be_old
Definition: ComputeSimoHughesJ2PlasticityStress.h:53
ComputeSimoHughesJ2PlasticityStress::_F_old
const MaterialProperty< RankTwoTensor > & _F_old
Definition: ComputeSimoHughesJ2PlasticityStress.h:48
DerivativeMaterialInterface< ComputeLagrangianStressPK1 >::derivativePropertyName
const MaterialPropertyName derivativePropertyName(const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
SingleVariableReturnMappingSolution
SingleVariableReturnMappingSolutionTempl< false > SingleVariableReturnMappingSolution
Definition: SingleVariableReturnMappingSolution.h:283

Member Function Documentation

◆ computeDerivative() [1/2]

Real ComputeSimoHughesJ2PlasticityStress::computeDerivative ( const Real &  effective_trial_stress,
const Real &  scalar 
)
overrideprotectedvirtual

Definition at line 167 of file ComputeSimoHughesJ2PlasticityStress.C.

169 {
170  const Real G = ElasticityTensorTools::getIsotropicShearModulus(_elasticity_tensor[_qp]);
171 
172  // Update the flow stress
173  _ep[_qp] = _ep_old[_qp] + scalar;
174  _flow_stress_material->computePropertiesAtQp(_qp);
175 
176  return -G * _be[_qp].trace() - _dH[_qp];
177 }
ComputeSimoHughesJ2PlasticityStress::_elasticity_tensor
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Definition: ComputeSimoHughesJ2PlasticityStress.h:46
ElasticityTensorTools::getIsotropicShearModulus
T getIsotropicShearModulus(const RankFourTensorTempl< T > &elasticity_tensor)
Get the shear modulus for an isotropic elasticity tensor param elasticity_tensor the tensor (must be ...
Definition: ElasticityTensorTools.h:70
NS::G
static const std::string G
Definition: NS.h:170
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
MaterialBase::computePropertiesAtQp
virtual void computePropertiesAtQp(unsigned int qp)

◆ computeDerivative() [2/2]

template<bool is_ad>
virtual GenericReal<is_ad> SingleVariableReturnMappingSolutionTempl< is_ad >::computeDerivative ( const GenericReal< is_ad > &  ,
const GenericReal< is_ad > &   
)
protectedpure virtualinherited

Compute the derivative of the residual as a function of the scalar variable.

The residual should be in strain increment units for all models for consistency.

Parameters
effective_trial_stressEffective trial stress
scalarInelastic strain increment magnitude being solved for

Implemented in LAROMANCEStressUpdateBaseTempl< is_ad >, CompositePowerLawCreepStressUpdateTempl< is_ad >, IsotropicPlasticityStressUpdateTempl< is_ad >, PowerLawCreepStressUpdateTempl< is_ad >, PowerLawCreepTestTempl< is_ad >, and CombinedNonlinearHardeningPlasticityTempl< is_ad >.

◆ computeQpPK1Stress()

void ComputeSimoHughesJ2PlasticityStress::computeQpPK1Stress ( )
overrideprotectedvirtual

Definition at line 73 of file ComputeSimoHughesJ2PlasticityStress.C.

74 {
75  usingTensorIndices(i, j, k, l, m);
76  const Real G = ElasticityTensorTools::getIsotropicShearModulus(_elasticity_tensor[_qp]);
77  const Real K = ElasticityTensorTools::getIsotropicBulkModulus(_elasticity_tensor[_qp]);
78  const auto I = RankTwoTensor::Identity();
79  const auto Fit = _F[_qp].inverse().transpose();
80  const auto detJ = _F[_qp].det();
81 
82  // Update configuration
83  RankTwoTensor f = _inv_df[_qp].inverse();
84  RankTwoTensor f_bar = f / std::cbrt(f.det());
85 
86  // Elastic predictor
87  _be[_qp] = f_bar * _be_old[_qp] * f_bar.transpose();
88  RankTwoTensor s = G * _be[_qp].deviatoric();
89  _Np[_qp] = MooseUtils::absoluteFuzzyEqual(s.norm(), 0) ? std::sqrt(1. / 2.) * I
90  : std::sqrt(3. / 2.) * s / s.norm();
91  Real s_eff = s.doubleContraction(_Np[_qp]);
92 
93  // Compute the derivative of the strain before return mapping
94  if (_fe_problem.currentlyComputingJacobian())
95  _d_be_d_F = _F_old[_qp].inverse().times<l, m, i, j, k, m>(
96  (I.times<i, k, j, l>(f_bar * _be_old[_qp].transpose()) +
97  I.times<j, k, i, l>(f_bar * _be_old[_qp])) /
98  std::cbrt(f.det()) -
99  2. / 3. * _be[_qp].times<i, j, l, k>(_inv_df[_qp]));
100 
101  // Check for plastic loading and do return mapping
102  Real delta_ep = 0;
103  if (computeResidual(s_eff, 0) > 0)
104  {
105  // Initialize the derivative of the internal variable
106  if (_fe_problem.currentlyComputingJacobian())
107  {
108  _d_deltaep_d_betr.zero();
109  if (MooseUtils::absoluteFuzzyEqual(s.norm(), 0))
110  _d_n_d_be.zero();
111  else
112  _d_n_d_be = G / std::sqrt(6) / s.norm() *
113  (3 * I.times<i, k, j, l>(I) - 2 * _Np[_qp].times<i, j, k, l>(_Np[_qp]) -
114  I.times<i, j, k, l>(I));
115  }
116 
117  returnMappingSolve(s_eff, delta_ep, _console);
118 
119  // Correct the derivative of the strain after return mapping
120  if (_fe_problem.currentlyComputingJacobian())
121  _d_be_d_F -=
122  2. / 3. *
123  (_be[_qp].trace() * _Np[_qp].times<i, j, k, l>(_d_deltaep_d_betr) +
124  delta_ep * _Np[_qp].times<i, j, k, l>(I) + delta_ep * _be[_qp].trace() * _d_n_d_be) *
125  _d_be_d_F;
126  }
127 
128  // Update intermediate and current configurations
129  _ep[_qp] = _ep_old[_qp] + delta_ep;
130  _be[_qp] -= 2. / 3. * delta_ep * _be[_qp].trace() * _Np[_qp];
131  s = G * _be[_qp].deviatoric();
132  RankTwoTensor tau = (K * (detJ * detJ - 1) / 2) * I + s;
133  _pk1_stress[_qp] = tau * Fit;
134 
135  // Compute the consistent tangent, i.e. the derivative of the PK1 stress w.r.t. the deformation
136  // gradient.
137  if (_fe_problem.currentlyComputingJacobian())
138  {
139  RankFourTensor d_tau_d_F = K * detJ * detJ * I.times<i, j, k, l>(Fit) +
140  G * (_d_be_d_F - I.times<i, j, k, l>(I) * _d_be_d_F / 3);
141  _pk1_jacobian[_qp] = Fit.times<m, j, i, m, k, l>(d_tau_d_F) - Fit.times<k, j, i, l>(tau * Fit);
142  }
143 }
ComputeSimoHughesJ2PlasticityStress::_elasticity_tensor
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Definition: ComputeSimoHughesJ2PlasticityStress.h:46
ElasticityTensorTools::getIsotropicShearModulus
T getIsotropicShearModulus(const RankFourTensorTempl< T > &elasticity_tensor)
Get the shear modulus for an isotropic elasticity tensor param elasticity_tensor the tensor (must be ...
Definition: ElasticityTensorTools.h:70
NS::K
static const std::string K
Definition: NS.h:174
trace
Real trace(const RealTensor &A, const unsigned int &dim)
inverse
void inverse(const std::vector< std::vector< Real >> &m, std::vector< std::vector< Real >> &m_inv)
ComputeSimoHughesJ2PlasticityStress::computeResidual
virtual Real computeResidual(const Real &effective_trial_stress, const Real &scalar) override
Definition: ComputeSimoHughesJ2PlasticityStress.C:154
RankTwoTensorTempl< Real >::Identity
static RankTwoTensorTempl Identity()
ComputeSimoHughesJ2PlasticityStress::_be_old
const MaterialProperty< RankTwoTensor > & _be_old
Definition: ComputeSimoHughesJ2PlasticityStress.h:53
NS::G
static const std::string G
Definition: NS.h:170
ElasticityTensorTools::getIsotropicBulkModulus
T getIsotropicBulkModulus(const RankFourTensorTempl< T > &elasticity_tensor)
Get the bulk modulus for an isotropic elasticity tensor param elasticity_tensor the tensor (must be i...
Definition: ElasticityTensorTools.h:81
SingleVariableReturnMappingSolutionTempl::returnMappingSolve
void returnMappingSolve(const GenericReal< is_ad > &effective_trial_stress, GenericReal< is_ad > &scalar, const ConsoleStream &console)
Perform the return mapping iterations.
Definition: SingleVariableReturnMappingSolution.C:100
ComputeSimoHughesJ2PlasticityStress::_F_old
const MaterialProperty< RankTwoTensor > & _F_old
Definition: ComputeSimoHughesJ2PlasticityStress.h:48
f
Real f(Real x)
Test function for Brents method.
Definition: BrentsMethodTest.C:21
RankTwoTensorTempl< Real >::doubleContraction
Real doubleContraction(const RankTwoTensorTempl< Real > &a) const
RankTwoTensorTempl< Real >::transpose
RankTwoTensorTempl< Real > transpose() const
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
ComputeSimoHughesJ2PlasticityStress::_d_be_d_F
RankFourTensor _d_be_d_F
Helper (dummy) variables for iteratively updating the consistant tangent during return mapping...
Definition: ComputeSimoHughesJ2PlasticityStress.h:64
NS::k
static const std::string k
Definition: NS.h:134

◆ computeReferenceResidual() [1/2]

Real ComputeSimoHughesJ2PlasticityStress::computeReferenceResidual ( const Real &  effective_trial_stress,
const Real &  scalar 
)
overrideprotectedvirtual

The return mapping residual and derivative.

Definition at line 146 of file ComputeSimoHughesJ2PlasticityStress.C.

148 {
149  const Real G = ElasticityTensorTools::getIsotropicShearModulus(_elasticity_tensor[_qp]);
150  return effective_trial_stress - G * scalar * _be[_qp].trace();
151 }
ComputeSimoHughesJ2PlasticityStress::_elasticity_tensor
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Definition: ComputeSimoHughesJ2PlasticityStress.h:46
ElasticityTensorTools::getIsotropicShearModulus
T getIsotropicShearModulus(const RankFourTensorTempl< T > &elasticity_tensor)
Get the shear modulus for an isotropic elasticity tensor param elasticity_tensor the tensor (must be ...
Definition: ElasticityTensorTools.h:70
NS::G
static const std::string G
Definition: NS.h:170
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ computeReferenceResidual() [2/2]

template<bool is_ad>
virtual Real SingleVariableReturnMappingSolutionTempl< is_ad >::computeReferenceResidual ( const GenericReal< is_ad > &  effective_trial_stress,
const GenericReal< is_ad > &  scalar 
)
protectedpure virtualinherited

Compute a reference quantity to be used for checking relative convergence.

This should be in strain increment units for all models for consistency.

Parameters
effective_trial_stressEffective trial stress
scalarInelastic strain increment magnitude being solved for

Implemented in RadialReturnStressUpdateTempl< is_ad >.

◆ computeResidual() [1/2]

Real ComputeSimoHughesJ2PlasticityStress::computeResidual ( const Real &  effective_trial_stress,
const Real &  scalar 
)
overrideprotectedvirtual

Definition at line 154 of file ComputeSimoHughesJ2PlasticityStress.C.

Referenced by computeQpPK1Stress().

156 {
157  const Real G = ElasticityTensorTools::getIsotropicShearModulus(_elasticity_tensor[_qp]);
158 
159  // Update the flow stress
160  _ep[_qp] = _ep_old[_qp] + scalar;
161  _flow_stress_material->computePropertiesAtQp(_qp);
162 
163  return effective_trial_stress - G * scalar * _be[_qp].trace() - _H[_qp];
164 }
ComputeSimoHughesJ2PlasticityStress::_elasticity_tensor
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Definition: ComputeSimoHughesJ2PlasticityStress.h:46
ElasticityTensorTools::getIsotropicShearModulus
T getIsotropicShearModulus(const RankFourTensorTempl< T > &elasticity_tensor)
Get the shear modulus for an isotropic elasticity tensor param elasticity_tensor the tensor (must be ...
Definition: ElasticityTensorTools.h:70
NS::G
static const std::string G
Definition: NS.h:170
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
MaterialBase::computePropertiesAtQp
virtual void computePropertiesAtQp(unsigned int qp)

◆ computeResidual() [2/2]

template<bool is_ad>
virtual GenericReal<is_ad> SingleVariableReturnMappingSolutionTempl< is_ad >::computeResidual ( const GenericReal< is_ad > &  ,
const GenericReal< is_ad > &   
)
protectedpure virtualinherited

Compute the residual for a predicted value of the scalar.

This residual should be in strain increment units for all models for consistency.

Parameters
effective_trial_stressEffective trial stress
scalarInelastic strain increment magnitude being solved for

Implemented in LAROMANCEStressUpdateBaseTempl< is_ad >, IsotropicPlasticityStressUpdateTempl< is_ad >, CompositePowerLawCreepStressUpdateTempl< is_ad >, PowerLawCreepStressUpdateTempl< is_ad >, PowerLawCreepTestTempl< is_ad >, and CombinedNonlinearHardeningPlasticityTempl< is_ad >.

◆ computeResidualAndDerivative()

template<bool is_ad>
virtual GenericChainedReal<is_ad> SingleVariableReturnMappingSolutionTempl< is_ad >::computeResidualAndDerivative ( const GenericReal< is_ad > &  ,
const GenericChainedReal< is_ad > &   
)
inlineprotectedvirtualinherited

Compute the residual and the derivative for a predicted value of the scalar.

This residual should be in strain increment units for all models for consistency.

Parameters
effective_trial_stressEffective trial stress
scalarInelastic strain increment magnitude being solved for

Reimplemented in CompositePowerLawCreepStressUpdateTempl< is_ad >, and PowerLawCreepStressUpdateTempl< is_ad >.

Definition at line 97 of file SingleVariableReturnMappingSolution.h.

99  {
100  mooseError("computeResidualAndDerivative has to be implemented if "
101  "automatic_differentiation_return_mapping = true.");
102  return 0;
103  };
mooseError
void mooseError(Args &&... args)

◆ converged()

template<bool is_ad>
bool SingleVariableReturnMappingSolutionTempl< is_ad >::converged ( const GenericReal< is_ad > &  residual,
const Real  reference 
)
protectedinherited

Check to see whether the residual is within the convergence limits.

Parameters
residualCurrent value of the residual
referenceCurrent value of the reference quantity
Returns
Whether the model converged

Definition at line 327 of file SingleVariableReturnMappingSolution.C.

329 {
330  const Real residual = MetaPhysicL::raw_value(ad_residual);
331  return (std::abs(residual) <= _absolute_tolerance ||
332  std::abs(residual / reference) <= _relative_tolerance);
333 }
MetaPhysicL::raw_value
auto raw_value(const Eigen::Map< T > &in)
SingleVariableReturnMappingSolutionTempl::_relative_tolerance
Real _relative_tolerance
Relative convergence tolerance.
Definition: SingleVariableReturnMappingSolution.h:194
SingleVariableReturnMappingSolutionTempl::_absolute_tolerance
Real _absolute_tolerance
Absolute convergence tolerance.
Definition: SingleVariableReturnMappingSolution.h:197
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ initialGuess()

template<bool is_ad>
virtual GenericReal<is_ad> SingleVariableReturnMappingSolutionTempl< is_ad >::initialGuess ( const GenericReal< is_ad > &  )
inlineprotectedvirtualinherited

Compute an initial guess for the value of the scalar.

For some cases, an intellegent starting point can provide enhanced robustness in the Newton iterations. This is also an opportunity for classes that derive from this to perform initialization tasks.

Parameters
effective_trial_stressEffective trial stress

Reimplemented in PowerLawCreepTestTempl< is_ad >.

Definition at line 66 of file SingleVariableReturnMappingSolution.h.

67  {
68  return 0.0;
69  }

◆ initialSetup()

void ComputeSimoHughesJ2PlasticityStress::initialSetup ( )
overridevirtual

Definition at line 53 of file ComputeSimoHughesJ2PlasticityStress.C.

54 {
55  ComputeLagrangianStressBase::initialSetup();
56 
57  _flow_stress_material = &getMaterial("flow_stress_material");
58 
59  // Enforce isotropic elastic tensor
60  if (!hasGuaranteedMaterialProperty(_elasticity_tensor_name, Guarantee::ISOTROPIC))
61  mooseError("ComputeSimoHughesJ2PlasticityStress requires an isotropic elasticity tensor");
62 }
mooseError
void mooseError(Args &&... args)
ComputeLagrangianStressBase::initialSetup
virtual void initialSetup() override
Derive _large_kinematics from the strain calculator&#39;s LARGE_KINEMATICS guarantee. ...
Definition: ComputeLagrangianStressBase.C:59

◆ initQpStatefulProperties()

void ComputeSimoHughesJ2PlasticityStress::initQpStatefulProperties ( )
overrideprotectedvirtual

Definition at line 65 of file ComputeSimoHughesJ2PlasticityStress.C.

66 {
67  ComputeLagrangianStressPK1::initQpStatefulProperties();
68  _be[_qp].setToIdentity();
69  _ep[_qp] = 0;
70 }
ComputeLagrangianStressBase::initQpStatefulProperties
virtual void initQpStatefulProperties() override
Initialize everything with zeros.
Definition: ComputeLagrangianStressBase.C:81

◆ iterationFinalize()

template<bool is_ad>
virtual void SingleVariableReturnMappingSolutionTempl< is_ad >::iterationFinalize ( const GenericReal< is_ad > &  )
inlineprotectedvirtualinherited

Finalize internal state variables for a model for a given iteration.

Parameters
scalarInelastic strain increment magnitude being solved for

Reimplemented in IsotropicPlasticityStressUpdateTempl< is_ad >, and CombinedNonlinearHardeningPlasticityTempl< is_ad >.

Definition at line 128 of file SingleVariableReturnMappingSolution.h.

128 {}

◆ maximumPermissibleValue()

template<bool is_ad>
GenericReal< is_ad > SingleVariableReturnMappingSolutionTempl< is_ad >::maximumPermissibleValue ( const GenericReal< is_ad > &  effective_trial_stress) const
protectedvirtualinherited

Compute the maximum permissible value of the scalar.

For some models, the magnitude of this may be known.

Parameters
effective_trial_stressEffective trial stress

Reimplemented in RadialReturnStressUpdateTempl< is_ad >, LAROMANCEStressUpdateBaseTempl< is_ad >, and PowerLawCreepTestTempl< is_ad >.

Definition at line 92 of file SingleVariableReturnMappingSolution.C.

94 {
95  return std::numeric_limits<Real>::max();
96 }

◆ minimumPermissibleValue()

template<bool is_ad>
GenericReal< is_ad > SingleVariableReturnMappingSolutionTempl< is_ad >::minimumPermissibleValue ( const GenericReal< is_ad > &  effective_trial_stress) const
protectedvirtualinherited

Compute the minimum permissible value of the scalar.

For some models, the magnitude of this may be known.

Parameters
effective_trial_stressEffective trial stress

Reimplemented in RadialReturnStressUpdateTempl< is_ad >, and PowerLawCreepTestTempl< is_ad >.

Definition at line 84 of file SingleVariableReturnMappingSolution.C.

86 {
87  return std::numeric_limits<Real>::lowest();
88 }

◆ outputIterationStep()

template<bool is_ad>
void SingleVariableReturnMappingSolutionTempl< is_ad >::outputIterationStep ( std::stringstream *  iter_output,
const GenericReal< is_ad > &  effective_trial_stress,
const GenericReal< is_ad > &  scalar,
const Real  reference_residual 
)
protectedvirtualinherited

Output information for a single iteration step to build the convergence history of the model.

Parameters
iter_outputOutput stream
effective_trial_stressEffective trial stress
residualCurrent value of the residual
referenceCurrent value of the reference quantity

Reimplemented in LAROMANCEStressUpdateBaseTempl< is_ad >.

Definition at line 421 of file SingleVariableReturnMappingSolution.C.

Referenced by LAROMANCEStressUpdateBaseTempl< is_ad >::outputIterationStep().

426 {
427  if (iter_output)
428  {
429  const unsigned int it = _iteration;
430  const Real residual = MetaPhysicL::raw_value(_residual);
431 
432  *iter_output << " iteration=" << it
433  << " trial_stress=" << MetaPhysicL::raw_value(effective_trial_stress)
434  << " scalar=" << MetaPhysicL::raw_value(scalar) << " residual=" << residual
435  << " ref_res=" << reference_residual
436  << " rel_res=" << std::abs(residual) / reference_residual
437  << " rel_tol=" << _relative_tolerance << " abs_res=" << std::abs(residual)
438  << " abs_tol=" << _absolute_tolerance << '\n';
439  }
440 }
MetaPhysicL::raw_value
auto raw_value(const Eigen::Map< T > &in)
SingleVariableReturnMappingSolutionTempl::_relative_tolerance
Real _relative_tolerance
Relative convergence tolerance.
Definition: SingleVariableReturnMappingSolution.h:194
SingleVariableReturnMappingSolutionTempl::_absolute_tolerance
Real _absolute_tolerance
Absolute convergence tolerance.
Definition: SingleVariableReturnMappingSolution.h:197
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ outputIterationSummary()

template<bool is_ad>
void SingleVariableReturnMappingSolutionTempl< is_ad >::outputIterationSummary ( std::stringstream *  iter_output,
const unsigned int  total_it 
)
protectedvirtualinherited

Output summary information for the convergence history of the model.

Parameters
iter_outputOutput stream
total_itTotal iteration count

Reimplemented in RadialReturnStressUpdateTempl< is_ad >, LAROMANCEStressUpdateBaseTempl< is_ad >, and ADViscoplasticityStressUpdate.

Definition at line 444 of file SingleVariableReturnMappingSolution.C.

Referenced by ADViscoplasticityStressUpdate::outputIterationSummary(), LAROMANCEStressUpdateBaseTempl< is_ad >::outputIterationSummary(), and RadialReturnStressUpdateTempl< is_ad >::outputIterationSummary().

446 {
447  if (iter_output)
448  *iter_output << "In " << total_it << " iterations the residual went from "
449  << MetaPhysicL::raw_value(_initial_residual) << " to "
450  << MetaPhysicL::raw_value(_residual) << " in '" << _svrms_name << "'."
451  << std::endl;
452 }
MetaPhysicL::raw_value
auto raw_value(const Eigen::Map< T > &in)
SingleVariableReturnMappingSolutionTempl::_initial_residual
GenericReal< is_ad > _initial_residual
Residual values, kept as members to retain solver state for summary outputting.
Definition: SingleVariableReturnMappingSolution.h:217
SingleVariableReturnMappingSolutionTempl::_svrms_name
const std::string _svrms_name
MOOSE input name of the object performing the solve.
Definition: SingleVariableReturnMappingSolution.h:225

◆ preStep() [1/2]

void ComputeSimoHughesJ2PlasticityStress::preStep ( const Real &  scalar_old,
const Real &  residual,
const Real &  jacobian 
)
overrideprotectedvirtual

Definition at line 180 of file ComputeSimoHughesJ2PlasticityStress.C.

181 {
182  if (!_fe_problem.currentlyComputingJacobian())
183  return;
184 
185  const auto I = RankTwoTensor::Identity();
186  const Real G = ElasticityTensorTools::getIsotropicShearModulus(_elasticity_tensor[_qp]);
187 
188  // Update the flow stress
189  _ep[_qp] = _ep_old[_qp] + scalar;
190  _flow_stress_material->computePropertiesAtQp(_qp);
191 
192  _d_R_d_betr =
193  G * _Np[_qp] - G * scalar * I - (G * _be[_qp].trace() + _dH[_qp]) * _d_deltaep_d_betr;
194  _d_J_d_betr = -G * I - _d2H[_qp] * _d_deltaep_d_betr;
195  _d_deltaep_d_betr += -1 / J * _d_R_d_betr + R / J / J * _d_J_d_betr;
196 }
ComputeSimoHughesJ2PlasticityStress::_elasticity_tensor
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Definition: ComputeSimoHughesJ2PlasticityStress.h:46
ElasticityTensorTools::getIsotropicShearModulus
T getIsotropicShearModulus(const RankFourTensorTempl< T > &elasticity_tensor)
Get the shear modulus for an isotropic elasticity tensor param elasticity_tensor the tensor (must be ...
Definition: ElasticityTensorTools.h:70
R
const double R
Definition: GasLiquidMassTransferTest.C:21
RankTwoTensorTempl< Real >::Identity
static RankTwoTensorTempl Identity()
NS::G
static const std::string G
Definition: NS.h:170
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
MaterialBase::computePropertiesAtQp
virtual void computePropertiesAtQp(unsigned int qp)

◆ preStep() [2/2]

template<bool is_ad>
virtual void SingleVariableReturnMappingSolutionTempl< is_ad >::preStep ( const GenericReal< is_ad > &  ,
const GenericReal< is_ad > &  ,
const GenericReal< is_ad > &   
)
inlineprotectedvirtualinherited

This method is called before taking a step in the return mapping algorithm.

A typical use case is to accumulate the exact algorithmic tangent during return mapping.

Definition at line 118 of file SingleVariableReturnMappingSolution.h.

121  {
122  }

◆ returnMappingSolve()

template<bool is_ad>
void SingleVariableReturnMappingSolutionTempl< is_ad >::returnMappingSolve ( const GenericReal< is_ad > &  effective_trial_stress,
GenericReal< is_ad > &  scalar,
const ConsoleStream &  console 
)
protectedinherited

Perform the return mapping iterations.

Parameters
effective_trial_stressEffective trial stress
scalarInelastic strain increment magnitude being solved for
consoleConsole output

Definition at line 100 of file SingleVariableReturnMappingSolution.C.

Referenced by ADViscoplasticityStressUpdate::computeInelasticStrainIncrement(), and computeQpPK1Stress().

104 {
105  // construct the stringstream here only if the debug level is set to ALL
106  std::unique_ptr<std::stringstream> iter_output =
107  (_internal_solve_output_on == InternalSolveOutput::ALWAYS)
108  ? std::make_unique<std::stringstream>()
109  : nullptr;
110 
111  // do the internal solve and capture iteration info during the first round
112  // iff full history output is requested regardless of whether the solve failed or succeeded
113  auto solve_state =
114  internalSolve(effective_trial_stress,
115  scalar,
116  _internal_solve_full_iteration_history ? iter_output.get() : nullptr);
117  if (solve_state != SolveState::SUCCESS &&
118  _internal_solve_output_on != InternalSolveOutput::ALWAYS)
119  {
120  // output suppressed by user, throw immediately
121  if (_internal_solve_output_on == InternalSolveOutput::NEVER)
122  mooseException("");
123 
124  // user expects some kind of output, if necessary setup output stream now
125  if (!iter_output)
126  iter_output = std::make_unique<std::stringstream>();
127 
128  // add the appropriate error message to the output
129  switch (solve_state)
130  {
131  case SolveState::NAN_INF:
132  *iter_output << "Encountered inf or nan in material return mapping iterations.\n";
133  break;
134 
135  case SolveState::EXCEEDED_ITERATIONS:
136  *iter_output << "Exceeded maximum iterations in material return mapping iterations.\n";
137  break;
138 
139  default:
140  mooseError("Unhandled solver state");
141  }
142 
143  // if full history output is only requested for failed solves we have to repeat
144  // the solve a second time
145  if (_internal_solve_full_iteration_history)
146  internalSolve(effective_trial_stress, scalar, iter_output.get());
147 
148  // Append summary and throw exception
149  outputIterationSummary(iter_output.get(), _iteration);
150  mooseException(iter_output->str());
151  }
152 
153  if (_internal_solve_output_on == InternalSolveOutput::ALWAYS)
154  {
155  // the solve did not fail but the user requested debug output anyways
156  outputIterationSummary(iter_output.get(), _iteration);
157  console << iter_output->str() << std::flush;
158  }
159 }
mooseError
void mooseError(Args &&... args)
SingleVariableReturnMappingSolutionTempl::internalSolve
SolveState internalSolve(const GenericReal< is_ad > effective_trial_stress, GenericReal< is_ad > &scalar, std::stringstream *iter_output=nullptr)
Method called from within this class to perform the actual return mappping iterations.
Definition: SingleVariableReturnMappingSolution.C:163
SingleVariableReturnMappingSolutionTempl::_internal_solve_full_iteration_history
const bool _internal_solve_full_iteration_history
Whether to output iteration information all the time (regardless of whether iterations converge) ...
Definition: SingleVariableReturnMappingSolution.h:191
SingleVariableReturnMappingSolutionTempl::outputIterationSummary
virtual void outputIterationSummary(std::stringstream *iter_output, const unsigned int total_it)
Output summary information for the convergence history of the model.
Definition: SingleVariableReturnMappingSolution.C:444
SingleVariableReturnMappingSolutionTempl::_internal_solve_output_on
enum SingleVariableReturnMappingSolutionTempl::InternalSolveOutput _internal_solve_output_on

◆ validParams()

InputParameters ComputeSimoHughesJ2PlasticityStress::validParams ( )
static

Definition at line 15 of file ComputeSimoHughesJ2PlasticityStress.C.

16 {
17  InputParameters params = DerivativeMaterialInterface<ComputeLagrangianStressPK1>::validParams();
18  params += SingleVariableReturnMappingSolution::validParams();
19  params.addClassDescription("The Simo-Hughes style J2 plasticity.");
20  params.addParam<MaterialPropertyName>(
21  "elasticity_tensor", "elasticity_tensor", "The name of the elasticity tensor.");
22  params.addRequiredParam<MaterialName>("flow_stress_material",
23  "The material defining the flow stress");
24  return params;
25 }
InputParameters::addParam
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
InputParameters::addRequiredParam
void addRequiredParam(const std::string &name, const std::string &doc_string)
validParams
InputParameters validParams()
InputParameters::addClassDescription
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _be

MaterialProperty<RankTwoTensor>& ComputeSimoHughesJ2PlasticityStress::_be
protected

Definition at line 52 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeDerivative(), computeQpPK1Stress(), computeReferenceResidual(), computeResidual(), initQpStatefulProperties(), and preStep().

◆ _be_old

const MaterialProperty<RankTwoTensor>& ComputeSimoHughesJ2PlasticityStress::_be_old
protected

Definition at line 53 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeQpPK1Stress().

◆ _bracket_solution

template<bool is_ad>
bool SingleVariableReturnMappingSolutionTempl< is_ad >::_bracket_solution
protectedinherited

Whether to save upper and lower bounds of root for scalar, and set solution to the midpoint between those bounds if outside them.

Definition at line 145 of file SingleVariableReturnMappingSolution.h.

◆ _check_range

template<bool is_ad>
bool SingleVariableReturnMappingSolutionTempl< is_ad >::_check_range
protectedinherited

Whether to check to see whether iterative solution is within admissible range, and set within that range if outside.

Definition at line 138 of file SingleVariableReturnMappingSolution.h.

Referenced by ADViscoplasticityStressUpdate::ADViscoplasticityStressUpdate(), LAROMANCEStressUpdateBaseTempl< is_ad >::LAROMANCEStressUpdateBaseTempl(), and PowerLawCreepTestTempl< is_ad >::PowerLawCreepTestTempl().

◆ _d2H

const MaterialProperty<Real>& ComputeSimoHughesJ2PlasticityStress::_d2H
protected

Definition at line 60 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by preStep().

◆ _d_be_d_F

RankFourTensor ComputeSimoHughesJ2PlasticityStress::_d_be_d_F
private

Helper (dummy) variables for iteratively updating the consistant tangent during return mapping.

Definition at line 64 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeQpPK1Stress().

◆ _d_deltaep_d_betr

RankTwoTensor ComputeSimoHughesJ2PlasticityStress::_d_deltaep_d_betr
private

Definition at line 66 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeQpPK1Stress(), and preStep().

◆ _d_J_d_betr

RankTwoTensor ComputeSimoHughesJ2PlasticityStress::_d_J_d_betr
private

Definition at line 68 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by preStep().

◆ _d_n_d_be

RankFourTensor ComputeSimoHughesJ2PlasticityStress::_d_n_d_be
private

Definition at line 65 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeQpPK1Stress().

◆ _d_R_d_betr

RankTwoTensor ComputeSimoHughesJ2PlasticityStress::_d_R_d_betr
private

Definition at line 67 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by preStep().

◆ _dH

const MaterialProperty<Real>& ComputeSimoHughesJ2PlasticityStress::_dH
protected

Definition at line 59 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeDerivative(), and preStep().

◆ _elasticity_tensor

const MaterialProperty<RankFourTensor>& ComputeSimoHughesJ2PlasticityStress::_elasticity_tensor
protected

Definition at line 46 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeDerivative(), computeQpPK1Stress(), computeReferenceResidual(), computeResidual(), and preStep().

◆ _elasticity_tensor_name

const MaterialPropertyName ComputeSimoHughesJ2PlasticityStress::_elasticity_tensor_name
protected

Definition at line 45 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by initialSetup().

◆ _ep

MaterialProperty<Real>& ComputeSimoHughesJ2PlasticityStress::_ep
protected

Definition at line 50 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeDerivative(), computeQpPK1Stress(), computeResidual(), initQpStatefulProperties(), and preStep().

◆ _ep_name

const std::string ComputeSimoHughesJ2PlasticityStress::_ep_name
protected

Definition at line 49 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by ComputeSimoHughesJ2PlasticityStress().

◆ _ep_old

const MaterialProperty<Real>& ComputeSimoHughesJ2PlasticityStress::_ep_old
protected

Definition at line 51 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeDerivative(), computeQpPK1Stress(), computeResidual(), and preStep().

◆ _F_old

const MaterialProperty<RankTwoTensor>& ComputeSimoHughesJ2PlasticityStress::_F_old
protected

Definition at line 48 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeQpPK1Stress().

◆ _flow_stress_material

MaterialBase* ComputeSimoHughesJ2PlasticityStress::_flow_stress_material
protected

Definition at line 56 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeDerivative(), computeResidual(), initialSetup(), and preStep().

◆ _flow_stress_name

const std::string ComputeSimoHughesJ2PlasticityStress::_flow_stress_name
protected

Definition at line 57 of file ComputeSimoHughesJ2PlasticityStress.h.

◆ _H

const MaterialProperty<Real>& ComputeSimoHughesJ2PlasticityStress::_H
protected

Definition at line 58 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeResidual().

◆ _line_search

template<bool is_ad>
bool SingleVariableReturnMappingSolutionTempl< is_ad >::_line_search
protectedinherited

Whether to use line searches to improve convergence.

Definition at line 141 of file SingleVariableReturnMappingSolution.h.

◆ _Np

MaterialProperty<RankTwoTensor>& ComputeSimoHughesJ2PlasticityStress::_Np
protected

Definition at line 54 of file ComputeSimoHughesJ2PlasticityStress.h.

Referenced by computeQpPK1Stress(), and preStep().

The documentation for this class was generated from the following files:
Generated on Mon Aug 3 2026 21:33:02 for https://mooseframework.inl.gov by   doxygen 1.8.14