www.mooseframework.org
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
RadialReturnStressUpdate Class Referenceabstract

RadialReturnStressUpdate computes the radial return stress increment for an isotropic elastic-viscoplasticity model after interating on the difference between new and old trial stress increments. More...

#include <RadialReturnStressUpdate.h>

Inheritance diagram for RadialReturnStressUpdate:
[legend]

Public Member Functions

 RadialReturnStressUpdate (const InputParameters &parameters)
 
virtual void updateState (RankTwoTensor &strain_increment, RankTwoTensor &inelastic_strain_increment, const RankTwoTensor &rotation_increment, RankTwoTensor &stress_new, const RankTwoTensor &stress_old, const RankFourTensor &elasticity_tensor, const RankTwoTensor &elastic_strain_old, bool compute_full_tangent_operator, RankFourTensor &tangent_operator) override
 A radial return (J2) mapping method is performed with return mapping iterations. More...
 
virtual Real computeReferenceResidual (const Real effective_trial_stress, const Real scalar_effective_inelastic_strain) override
 Compute a reference quantity to be used for checking relative convergence. More...
 
virtual Real minimumPermissibleValue (const Real) const override
 Compute the minimum permissible value of the scalar. More...
 
virtual Real maximumPermissibleValue (const Real effective_trial_stress) const override
 Compute the maximum permissible value of the scalar. More...
 
virtual Real computeTimeStepLimit () override
 Compute the limiting value of the time step for this material. More...
 
bool requiresIsotropicTensor () override
 Does the model require the elasticity tensor to be isotropic? More...
 
void setQp (unsigned int qp)
 Sets the value of the global variable _qp for inheriting classes. More...
 
virtual void propagateQpStatefulProperties ()
 If updateState is not called during a timestep, this will be. More...
 
virtual TangentCalculationMethod getTangentCalculationMethod ()
 
void resetQpProperties () final
 Retained as empty methods to avoid a warning from Material.C in framework. These methods are unused in all inheriting classes and should not be overwritten. More...
 
void resetProperties () final
 

Protected Member Functions

virtual void initQpStatefulProperties () override
 
void propagateQpStatefulPropertiesRadialReturn ()
 Propagate the properties pertaining to this intermediate class. More...
 
virtual void computeStressInitialize (const Real, const RankFourTensor &)
 Perform any necessary initialization before return mapping iterations. More...
 
virtual Real computeStressDerivative (const Real, const Real)
 Calculate the derivative of the strain increment with respect to the updated stress. More...
 
virtual void computeStressFinalize (const RankTwoTensor &)
 Perform any necessary steps to finalize state after return mapping iterations. More...
 
void outputIterationSummary (std::stringstream *iter_output, const unsigned int total_it) override
 Output summary information for the convergence history of the model. More...
 
void returnMappingSolve (const Real effective_trial_stress, Real &scalar, const ConsoleStream &console)
 Perform the return mapping iterations. More...
 
virtual Real initialGuess (const Real)
 Compute an initial guess for the value of the scalar. More...
 
virtual Real computeResidual (const Real effective_trial_stress, const Real scalar)=0
 Compute the residual for a predicted value of the scalar. More...
 
virtual Real computeDerivative (const Real effective_trial_stress, const Real scalar)=0
 Compute the derivative of the residual as a function of the scalar variable. More...
 
virtual void iterationFinalize (Real)
 Finalize internal state variables for a model for a given iteration. More...
 
virtual void outputIterationStep (std::stringstream *iter_output, const unsigned int it, const Real effective_trial_stress, const Real scalar, const Real residual, const Real reference_residual)
 Output information for a single iteration step to build the convergence history of the model. More...
 

Protected Attributes

Real _three_shear_modulus
 3 * shear modulus More...
 
MaterialProperty< Real > & _effective_inelastic_strain
 
const MaterialProperty< Real > & _effective_inelastic_strain_old
 
Real _max_inelastic_increment
 
const RankTwoTensor _identity_two
 Rank two identity tensor. More...
 
const RankFourTensor _identity_symmetric_four
 Rank four symmetric identity tensor. More...
 
const RankFourTensor _deviatoric_projection_four
 Rank four deviatoric projection tensor. More...
 
const std::string _base_name
 Name used as a prefix for all material properties related to the stress update model. More...
 
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...
 

Detailed Description

RadialReturnStressUpdate computes the radial return stress increment for an isotropic elastic-viscoplasticity model after interating on the difference between new and old trial stress increments.

This radial return mapping class acts as a base class for the radial return creep and plasticity classes / combinations. The stress increment computed by RadialReturnStressUpdate is used by ComputeMultipleInelasticStress which computes the elastic stress for finite strains. This return mapping class is acceptable for finite strains but not total strains. This class is based on the Elasto-viscoplasticity algorithm in F. Dunne and N. Petrinic's Introduction to Computational Plasticity (2004) Oxford University Press.

Definition at line 35 of file RadialReturnStressUpdate.h.

Constructor & Destructor Documentation

◆ RadialReturnStressUpdate()

RadialReturnStressUpdate::RadialReturnStressUpdate ( const InputParameters &  parameters)

Definition at line 36 of file RadialReturnStressUpdate.C.

37  : StressUpdateBase(parameters),
39  _effective_inelastic_strain(declareProperty<Real>(
40  _base_name + getParam<std::string>("effective_inelastic_strain_name"))),
41  _effective_inelastic_strain_old(getMaterialPropertyOld<Real>(
42  _base_name + getParam<std::string>("effective_inelastic_strain_name"))),
43  _max_inelastic_increment(parameters.get<Real>("max_inelastic_increment")),
44  _identity_two(RankTwoTensor::initIdentity),
45  _identity_symmetric_four(RankFourTensor::initIdentitySymmetricFour),
47  _identity_two.outerProduct(_identity_two) / 3.0)
48 {
49 }
const RankFourTensor _deviatoric_projection_four
Rank four deviatoric projection tensor.
SingleVariableReturnMappingSolution(const InputParameters &parameters)
const MaterialProperty< Real > & _effective_inelastic_strain_old
const RankFourTensor _identity_symmetric_four
Rank four symmetric identity tensor.
const RankTwoTensor _identity_two
Rank two identity tensor.
const std::string _base_name
Name used as a prefix for all material properties related to the stress update model.
MaterialProperty< Real > & _effective_inelastic_strain
StressUpdateBase(const InputParameters &parameters)

Member Function Documentation

◆ computeDerivative()

virtual Real SingleVariableReturnMappingSolution::computeDerivative ( const Real  effective_trial_stress,
const Real  scalar 
)
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 IsotropicPlasticityStressUpdate, PowerLawCreepStressUpdate, IsotropicPlasticity, CLSHPlasticModel, and PowerLawCreepModel.

Referenced by RadialReturnCreepStressUpdateBase::computeStressDerivative(), and SingleVariableReturnMappingSolution::internalSolve().

◆ computeReferenceResidual()

Real RadialReturnStressUpdate::computeReferenceResidual ( const Real  effective_trial_stress,
const Real  scalar 
)
overridevirtual

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

Implements SingleVariableReturnMappingSolution.

Definition at line 149 of file RadialReturnStressUpdate.C.

151 {
152  return effective_trial_stress / _three_shear_modulus - scalar_effective_inelastic_strain;
153 }
Real _three_shear_modulus
3 * shear modulus

◆ computeResidual()

virtual Real SingleVariableReturnMappingSolution::computeResidual ( const Real  effective_trial_stress,
const Real  scalar 
)
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 IsotropicPlasticityStressUpdate, PowerLawCreepStressUpdate, IsotropicPlasticity, CLSHPlasticModel, and PowerLawCreepModel.

Referenced by SingleVariableReturnMappingSolution::internalSolve().

◆ computeStressDerivative()

virtual Real RadialReturnStressUpdate::computeStressDerivative ( const Real  ,
const Real   
)
inlineprotectedvirtual

Calculate the derivative of the strain increment with respect to the updated stress.

Parameters
effective_trial_stressEffective trial stress
scalarInelastic strain increment magnitude being solved for

Reimplemented in RadialReturnCreepStressUpdateBase.

Definition at line 112 of file RadialReturnStressUpdate.h.

Referenced by updateState().

113  {
114  return 0.0;
115  }

◆ computeStressFinalize()

virtual void RadialReturnStressUpdate::computeStressFinalize ( const RankTwoTensor &  )
inlineprotectedvirtual

Perform any necessary steps to finalize state after return mapping iterations.

Parameters
inelasticStrainIncrementInelastic strain increment

Reimplemented in IsotropicPlasticityStressUpdate, and RadialReturnCreepStressUpdateBase.

Definition at line 121 of file RadialReturnStressUpdate.h.

Referenced by updateState().

121 {}

◆ computeStressInitialize()

virtual void RadialReturnStressUpdate::computeStressInitialize ( const Real  ,
const RankFourTensor &   
)
inlineprotectedvirtual

Perform any necessary initialization before return mapping iterations.

Parameters
effective_trial_stressEffective trial stress
elasticityTensorElasticity tensor

Reimplemented in IsotropicPlasticityStressUpdate, IsotropicPowerLawHardeningStressUpdate, PowerLawCreepStressUpdate, and TemperatureDependentHardeningStressUpdate.

Definition at line 102 of file RadialReturnStressUpdate.h.

Referenced by updateState().

104  {
105  }

◆ computeTimeStepLimit()

Real RadialReturnStressUpdate::computeTimeStepLimit ( )
overridevirtual

Compute the limiting value of the time step for this material.

Returns
Limiting time step

Reimplemented from StressUpdateBase.

Definition at line 162 of file RadialReturnStressUpdate.C.

163 {
164  Real scalar_inelastic_strain_incr;
165 
166  scalar_inelastic_strain_incr =
168  if (MooseUtils::absoluteFuzzyEqual(scalar_inelastic_strain_incr, 0.0))
169  return std::numeric_limits<Real>::max();
170 
171  return _dt * _max_inelastic_increment / scalar_inelastic_strain_incr;
172 }
const MaterialProperty< Real > & _effective_inelastic_strain_old
MaterialProperty< Real > & _effective_inelastic_strain

◆ getTangentCalculationMethod()

virtual TangentCalculationMethod StressUpdateBase::getTangentCalculationMethod ( )
inlinevirtualinherited

◆ initialGuess()

virtual Real SingleVariableReturnMappingSolution::initialGuess ( const Real  )
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

Definition at line 63 of file SingleVariableReturnMappingSolution.h.

Referenced by SingleVariableReturnMappingSolution::internalSolve().

63 { return 0.0; }

◆ initQpStatefulProperties()

void RadialReturnStressUpdate::initQpStatefulProperties ( )
overrideprotectedvirtual

◆ iterationFinalize()

virtual void SingleVariableReturnMappingSolution::iterationFinalize ( Real  )
inlineprotectedvirtualinherited

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

Parameters
scalarInelastic strain increment magnitude being solved for

Reimplemented in IsotropicPlasticityStressUpdate, IsotropicPlasticity, and CLSHPlasticModel.

Definition at line 93 of file SingleVariableReturnMappingSolution.h.

Referenced by SingleVariableReturnMappingSolution::internalSolve().

93 {}

◆ maximumPermissibleValue()

Real RadialReturnStressUpdate::maximumPermissibleValue ( const Real  effective_trial_stress) const
overridevirtual

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 from SingleVariableReturnMappingSolution.

Definition at line 156 of file RadialReturnStressUpdate.C.

157 {
158  return effective_trial_stress / _three_shear_modulus;
159 }
Real _three_shear_modulus
3 * shear modulus

◆ minimumPermissibleValue()

virtual Real RadialReturnStressUpdate::minimumPermissibleValue ( const Real  effective_trial_stress) const
inlineoverridevirtual

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 from SingleVariableReturnMappingSolution.

Definition at line 66 of file RadialReturnStressUpdate.h.

67  {
68  return 0.0;
69  }

◆ outputIterationStep()

void SingleVariableReturnMappingSolution::outputIterationStep ( std::stringstream *  iter_output,
const unsigned int  it,
const Real  effective_trial_stress,
const Real  scalar,
const Real  residual,
const Real  reference_residual 
)
protectedvirtualinherited

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

Parameters
iter_outputOutput stream
itCurrent iteration count
effective_trial_stressEffective trial stress
scalarInelastic strain increment magnitude being solved for
residualCurrent value of the residual
referenceCurrent value of the reference quantity

Definition at line 327 of file SingleVariableReturnMappingSolution.C.

Referenced by SingleVariableReturnMappingSolution::internalSolve().

333 {
334  if (iter_output)
335  {
336  *iter_output << " iteration=" << it << " trial_stress=" << effective_trial_stress
337  << " scalar=" << scalar << " residual=" << residual
338  << " ref_res=" << reference_residual
339  << " rel_res=" << std::abs(residual) / reference_residual
340  << " rel_tol=" << _relative_tolerance << " abs_res=" << std::abs(residual)
341  << " abs_tol=" << _absolute_tolerance << '\n';
342  }
343 }
Real _relative_tolerance
Relative convergence tolerance.
Real _absolute_tolerance
Absolute convergence tolerance.

◆ outputIterationSummary()

void RadialReturnStressUpdate::outputIterationSummary ( std::stringstream *  iter_output,
const unsigned int  total_it 
)
overrideprotectedvirtual

Output summary information for the convergence history of the model.

Parameters
iter_outputOutput stream
total_itTotal iteration count

Reimplemented from SingleVariableReturnMappingSolution.

Definition at line 175 of file RadialReturnStressUpdate.C.

177 {
178  if (iter_output)
179  {
180  *iter_output << "At element " << _current_elem->id() << " _qp=" << _qp << " Coordinates "
181  << _q_point[_qp] << " block=" << _current_elem->subdomain_id() << '\n';
182  }
184 }
virtual void outputIterationSummary(std::stringstream *iter_output, const unsigned int total_it)
Output summary information for the convergence history of the model.

◆ propagateQpStatefulProperties()

void StressUpdateBase::propagateQpStatefulProperties ( )
virtualinherited

If updateState is not called during a timestep, this will be.

This method allows derived classes to set internal parameters from their Old values, for instance

Reimplemented in MultiParameterPlasticityStressUpdate, LinearViscoelasticStressUpdate, IsotropicPlasticityStressUpdate, and RadialReturnCreepStressUpdateBase.

Definition at line 48 of file StressUpdateBase.C.

49 {
50  mooseError(
51  "propagateQpStatefulProperties called: it needs to be implemented by your inelastic model");
52 }

◆ propagateQpStatefulPropertiesRadialReturn()

void RadialReturnStressUpdate::propagateQpStatefulPropertiesRadialReturn ( )
protected

Propagate the properties pertaining to this intermediate class.

This is intended to be called from propagateQpStatefulProperties() in classes that inherit from this one. This is intentionally named uniquely because almost all models that derive from this class have their own stateful properties, and this forces them to define their own implementations of propagateQpStatefulProperties().

Definition at line 58 of file RadialReturnStressUpdate.C.

Referenced by RadialReturnCreepStressUpdateBase::propagateQpStatefulProperties(), and IsotropicPlasticityStressUpdate::propagateQpStatefulProperties().

59 {
61 }
const MaterialProperty< Real > & _effective_inelastic_strain_old
MaterialProperty< Real > & _effective_inelastic_strain

◆ requiresIsotropicTensor()

bool RadialReturnStressUpdate::requiresIsotropicTensor ( )
inlineoverridevirtual

Does the model require the elasticity tensor to be isotropic?

Implements StressUpdateBase.

Definition at line 82 of file RadialReturnStressUpdate.h.

82 { return true; }

◆ resetProperties()

void StressUpdateBase::resetProperties ( )
inlinefinalinherited

Definition at line 117 of file StressUpdateBase.h.

117 {}

◆ resetQpProperties()

void StressUpdateBase::resetQpProperties ( )
inlinefinalinherited

Retained as empty methods to avoid a warning from Material.C in framework. These methods are unused in all inheriting classes and should not be overwritten.

Definition at line 116 of file StressUpdateBase.h.

116 {}

◆ returnMappingSolve()

void SingleVariableReturnMappingSolution::returnMappingSolve ( const Real  effective_trial_stress,
Real &  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 91 of file SingleVariableReturnMappingSolution.C.

Referenced by ReturnMappingModel::computeStress(), and updateState().

94 {
95  // construct the stringstream here only if the debug level is set to ALL
96  std::stringstream * iter_output =
97  (_internal_solve_output_on == InternalSolveOutput::ALWAYS) ? new std::stringstream : nullptr;
98 
99  // do the internal solve and capture iteration info during the first round
100  // iff full history output is requested regardless of whether the solve failed or succeeded
101  auto solve_state = internalSolve(effective_trial_stress,
102  scalar,
103  _internal_solve_full_iteration_history ? iter_output : nullptr);
104  if (solve_state != SolveState::SUCCESS &&
106  {
107  // output suppressed by user, throw immediately
109  throw MooseException("");
110 
111  // user expects some kind of output, if necessary setup output stream now
112  if (!iter_output)
113  iter_output = new std::stringstream;
114 
115  // add the appropriate error message to the output
116  switch (solve_state)
117  {
118  case SolveState::NAN_INF:
119  *iter_output << "Encountered inf or nan in material return mapping iterations.\n";
120  break;
121 
123  *iter_output << "Exceeded maximum iterations in material return mapping iterations.\n";
124  break;
125 
126  default:
127  mooseError("Unhandled solver state");
128  }
129 
130  // if full history output is only requested for failed solves we have to repeat
131  // the solve a second time
133  internalSolve(effective_trial_stress, scalar, iter_output);
134 
135  // Append summary and throw exception
136  outputIterationSummary(iter_output, _iteration);
137  throw MooseException(iter_output->str());
138  }
139 
141  {
142  // the solve did not fail but the user requested debug output anyways
143  outputIterationSummary(iter_output, _iteration);
144  console << iter_output->str();
145  }
146 }
virtual void outputIterationSummary(std::stringstream *iter_output, const unsigned int total_it)
Output summary information for the convergence history of the model.
enum SingleVariableReturnMappingSolution::InternalSolveOutput _internal_solve_output_on
SolveState internalSolve(const Real effective_trial_stress, Real &scalar, std::stringstream *iter_output=nullptr)
Method called from within this class to perform the actual return mappping iterations.
const bool _internal_solve_full_iteration_history
Whether to output iteration information all the time (regardless of whether iterations converge) ...

◆ setQp()

void StressUpdateBase::setQp ( unsigned int  qp)
inherited

Sets the value of the global variable _qp for inheriting classes.

Definition at line 42 of file StressUpdateBase.C.

43 {
44  _qp = qp;
45 }

◆ updateState()

void RadialReturnStressUpdate::updateState ( RankTwoTensor &  strain_increment,
RankTwoTensor &  inelastic_strain_increment,
const RankTwoTensor &  rotation_increment,
RankTwoTensor &  stress_new,
const RankTwoTensor &  stress_old,
const RankFourTensor &  elasticity_tensor,
const RankTwoTensor &  elastic_strain_old,
bool  compute_full_tangent_operator,
RankFourTensor &  tangent_operator 
)
overridevirtual

A radial return (J2) mapping method is performed with return mapping iterations.

Parameters
strain_incrementSum of elastic and inelastic strain increments
inelastic_strain_incrementInelastic strain increment calculated by this class
rotationincrement Not used by this class
stress_newNew trial stress from pure elastic calculation
stress_oldOld state of stress
elasticity_tensorRank 4 C_{ijkl}, must be isotropic
elastic_strain_oldOld state of total elastic strain
compute_full_tangent_operatorFlag currently unused by this class
tangent_operatorCurrently a copy of the elasticity tensor in this class

Implements StressUpdateBase.

Definition at line 64 of file RadialReturnStressUpdate.C.

73 {
74  // compute the deviatoric trial stress and trial strain from the current intermediate
75  // configuration
76  RankTwoTensor deviatoric_trial_stress = stress_new.deviatoric();
77 
78  // compute the effective trial stress
79  Real dev_trial_stress_squared =
80  deviatoric_trial_stress.doubleContraction(deviatoric_trial_stress);
81  Real effective_trial_stress = std::sqrt(3.0 / 2.0 * dev_trial_stress_squared);
82 
83  // Set the value of 3 * shear modulus for use as a reference residual value
85 
86  computeStressInitialize(effective_trial_stress, elasticity_tensor);
87 
88  // Use Newton iteration to determine the scalar effective inelastic strain increment
89  Real scalar_effective_inelastic_strain = 0.0;
90  if (!MooseUtils::absoluteFuzzyEqual(effective_trial_stress, 0.0))
91  {
92  returnMappingSolve(effective_trial_stress, scalar_effective_inelastic_strain, _console);
93  if (scalar_effective_inelastic_strain != 0.0)
94  inelastic_strain_increment =
95  deviatoric_trial_stress *
96  (1.5 * scalar_effective_inelastic_strain / effective_trial_stress);
97  else
98  inelastic_strain_increment.zero();
99  }
100  else
101  inelastic_strain_increment.zero();
102 
103  strain_increment -= inelastic_strain_increment;
105  _effective_inelastic_strain_old[_qp] + scalar_effective_inelastic_strain;
106 
107  // Use the old elastic strain here because we require tensors used by this class
108  // to be isotropic and this method natively allows for changing in time
109  // elasticity tensors
110  stress_new = elasticity_tensor * (strain_increment + elastic_strain_old);
111 
112  computeStressFinalize(inelastic_strain_increment);
113 
114  if (compute_full_tangent_operator &&
116  {
117  if (MooseUtils::absoluteFuzzyEqual(scalar_effective_inelastic_strain, 0.0))
118  tangent_operator.zero();
119  else
120  {
121  // mu = _three_shear_modulus / 3.0;
122  // norm_dev_stress = ||s_n+1||
123  // effective_trial_stress = von mises trial stress = std::sqrt(3.0 / 2.0) * ||s_n+1^trial||
124  // scalar_effective_inelastic_strain = Delta epsilon^cr_n+1
125  // deriv = derivative of scalar_effective_inelastic_strain w.r.t. von mises stress
126  // deriv = std::sqrt(3.0 / 2.0) partial Delta epsilon^cr_n+1n over partial ||s_n+1^trial||
127 
128  mooseAssert(_three_shear_modulus != 0.0, "Shear modulus is zero");
129 
130  const RankTwoTensor deviatoric_stress = stress_new.deviatoric();
131  const Real norm_dev_stress =
132  std::sqrt(deviatoric_stress.doubleContraction(deviatoric_stress));
133  mooseAssert(norm_dev_stress != 0.0, "Norm of the deviatoric is zero");
134 
135  const RankTwoTensor flow_direction = deviatoric_stress / norm_dev_stress;
136  const RankFourTensor flow_direction_dyad = flow_direction.outerProduct(flow_direction);
137  const Real deriv =
138  computeStressDerivative(effective_trial_stress, scalar_effective_inelastic_strain);
139  const Real scalar_one = _three_shear_modulus * scalar_effective_inelastic_strain /
140  std::sqrt(1.5) / norm_dev_stress;
141 
142  tangent_operator = scalar_one * _deviatoric_projection_four +
143  (_three_shear_modulus * deriv - scalar_one) * flow_direction_dyad;
144  }
145  }
146 }
const RankFourTensor _deviatoric_projection_four
Rank four deviatoric projection tensor.
virtual void computeStressInitialize(const Real, const RankFourTensor &)
Perform any necessary initialization before return mapping iterations.
const MaterialProperty< Real > & _effective_inelastic_strain_old
virtual TangentCalculationMethod getTangentCalculationMethod()
virtual void computeStressFinalize(const RankTwoTensor &)
Perform any necessary steps to finalize state after return mapping iterations.
void returnMappingSolve(const Real effective_trial_stress, Real &scalar, const ConsoleStream &console)
Perform the return mapping iterations.
virtual Real computeStressDerivative(const Real, const Real)
Calculate the derivative of the strain increment with respect to the updated stress.
MaterialProperty< Real > & _effective_inelastic_strain
Real getIsotropicShearModulus(const RankFourTensor &elasticity_tensor)
Get the shear modulus for an isotropic elasticity tensor param elasticity_tensor the tensor (must be ...
Real _three_shear_modulus
3 * shear modulus

Member Data Documentation

◆ _base_name

const std::string StressUpdateBase::_base_name
protectedinherited

Name used as a prefix for all material properties related to the stress update model.

Definition at line 122 of file StressUpdateBase.h.

◆ _bracket_solution

bool SingleVariableReturnMappingSolution::_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 126 of file SingleVariableReturnMappingSolution.h.

Referenced by SingleVariableReturnMappingSolution::internalSolve().

◆ _check_range

bool SingleVariableReturnMappingSolution::_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 119 of file SingleVariableReturnMappingSolution.h.

Referenced by SingleVariableReturnMappingSolution::internalSolve().

◆ _deviatoric_projection_four

const RankFourTensor RadialReturnStressUpdate::_deviatoric_projection_four
protected

Rank four deviatoric projection tensor.

Definition at line 146 of file RadialReturnStressUpdate.h.

Referenced by updateState().

◆ _effective_inelastic_strain

MaterialProperty<Real>& RadialReturnStressUpdate::_effective_inelastic_strain
protected

◆ _effective_inelastic_strain_old

const MaterialProperty<Real>& RadialReturnStressUpdate::_effective_inelastic_strain_old
protected

◆ _identity_symmetric_four

const RankFourTensor RadialReturnStressUpdate::_identity_symmetric_four
protected

Rank four symmetric identity tensor.

Definition at line 141 of file RadialReturnStressUpdate.h.

◆ _identity_two

const RankTwoTensor RadialReturnStressUpdate::_identity_two
protected

Rank two identity tensor.

Definition at line 136 of file RadialReturnStressUpdate.h.

◆ _line_search

bool SingleVariableReturnMappingSolution::_line_search
protectedinherited

Whether to use line searches to improve convergence.

Definition at line 122 of file SingleVariableReturnMappingSolution.h.

Referenced by SingleVariableReturnMappingSolution::internalSolve().

◆ _max_inelastic_increment

Real RadialReturnStressUpdate::_max_inelastic_increment
protected

Definition at line 131 of file RadialReturnStressUpdate.h.

Referenced by computeTimeStepLimit().

◆ _three_shear_modulus

Real RadialReturnStressUpdate::_three_shear_modulus
protected

The documentation for this class was generated from the following files: