ComputeMultipleInelasticStress Class Reference

ComputeMultipleInelasticStress computes the stress, the consistent tangent operator (or an approximation to it), and a decomposition of the strain into elastic and inelastic parts. More...

#include <ComputeMultipleInelasticStress.h>

Inheritance diagram for ComputeMultipleInelasticStress:

Public Member Functions
	ComputeMultipleInelasticStress (const InputParameters &parameters)
virtual void	initialSetup () override

Static Public Member Functions
static InputParameters	validParams ()

Protected Types
enum	TangentOperatorEnum { TangentOperatorEnum::elastic, TangentOperatorEnum::nonlinear }
	what sort of Tangent operator to calculate More...

Protected Member Functions
virtual void	initQpStatefulProperties () override
virtual void	computeQpStress () override
	Compute the stress and store it in the _stress material property for the current quadrature point. More...
virtual void	computeQpStressIntermediateConfiguration ()
	Compute the stress for the current QP, but do not rotate tensors from the intermediate configuration to the new configuration. More...
virtual void	finiteStrainRotation (const bool force_elasticity_rotation=false)
	Rotate _elastic_strain, _stress, _inelastic_strain, and _Jacobian_mult to the new configuration. More...
virtual void	updateQpState (RankTwoTensor &elastic_strain_increment, RankTwoTensor &combined_inelastic_strain_increment)
	Given the _strain_increment[_qp], iterate over all of the user-specified recompute materials in order to find an admissible stress (which is placed into _stress[_qp]) and set of inelastic strains, as well as the tangent operator (which is placed into _Jacobian_mult[_qp]). More...
virtual void	updateQpStateSingleModel (unsigned model_number, RankTwoTensor &elastic_strain_increment, RankTwoTensor &combined_inelastic_strain_increment)
	An optimised version of updateQpState that gets used when the number of plastic models is unity, or when we're cycling through models Given the _strain_increment[_qp], find an admissible stress (which is put into _stress[_qp]) and inelastic strain, as well as the tangent operator (which is placed into _Jacobian_mult[_qp]) More...
virtual void	computeQpJacobianMult ()
	Using _elasticity_tensor[_qp] and the consistent tangent operators, _consistent_tangent_operator[...] computed by the inelastic models, compute _Jacobian_mult[_qp]. More...
virtual void	computeAdmissibleState (unsigned model_number, RankTwoTensor &elastic_strain_increment, RankTwoTensor &inelastic_strain_increment, RankFourTensor &consistent_tangent_operator)
	Given a trial stress (_stress[_qp]) and a strain increment (elastic_strain_increment) let the model_number model produce an admissible stress (gets placed back in _stress[_qp]), and decompose the strain increment into an elastic part (gets placed back into elastic_strain_increment) and an inelastic part (inelastic_strain_increment), as well as computing the consistent_tangent_operator. More...
virtual void	computeQpProperties () override
bool	hasGuaranteedMaterialProperty (const MaterialPropertyName &prop, Guarantee guarantee)

Protected Attributes
const bool	_perform_finite_strain_rotations
	after updateQpState, rotate the stress, elastic_strain, inelastic_strain and Jacobian_mult using _rotation_increment More...
MaterialProperty< RankTwoTensor > &	_inelastic_strain
	The sum of the inelastic strains that come from the plastic models. More...
const MaterialProperty< RankTwoTensor > &	_inelastic_strain_old
	old value of inelastic strain More...
enum ComputeMultipleInelasticStress::TangentOperatorEnum	_tangent_operator_type
const unsigned	_num_models
	number of plastic models More...
std::vector< bool >	_tangent_computation_flag
	Flags to compute tangent during updateState call. More...
TangentCalculationMethod	_tangent_calculation_method
	Calculation method for the tangent modulus. More...
const std::vector< Real >	_inelastic_weights
	_inelastic_strain = sum_i (_inelastic_weights_i * inelastic_strain_from_model_i) More...
std::vector< RankFourTensor >	_consistent_tangent_operator
	the consistent tangent operators computed by each plastic model More...
const bool	_cycle_models
	whether to cycle through the models, using only one model per timestep More...
MaterialProperty< Real > &	_matl_timestep_limit
const RankFourTensor	_identity_symmetric_four
	Rank four symmetric identity tensor. More...
std::vector< StressUpdateBase * >	_models
	The user supplied list of inelastic models to use in the simulation. More...
bool	_is_elasticity_tensor_guaranteed_isotropic
	is the elasticity tensor guaranteed to be isotropic? More...
bool	_all_models_isotropic
	are all inelastic models inherently isotropic? (not the case for e.g. weak plane plasticity models) More...
DamageBase *	_damage_model
	Pointer to the damage model. More...
RankTwoTensor	_undamaged_stress_old
const std::string	_elasticity_tensor_name
	Name of the elasticity tensor material property. More...
const MaterialProperty< RankFourTensor > &	_elasticity_tensor
	Elasticity tensor material property. More...
const MaterialProperty< RankTwoTensor > &	_rotation_increment
	Rotation increment material property. More...
const MaterialProperty< RankTwoTensor > &	_stress_old
	Old state of the stress tensor material property. More...
const std::string	_base_name
	Base name prepended to all material property names to allow for multi-material systems. More...
const MaterialProperty< RankTwoTensor > &	_mechanical_strain
	Mechanical strain material property. More...
MaterialProperty< RankTwoTensor > &	_stress
	Stress material property. More...
MaterialProperty< RankTwoTensor > &	_elastic_strain
	Elastic strain material property. More...
const MaterialProperty< RankTwoTensor > &	_extra_stress
	Extra stress tensor. More...
std::vector< const Function * >	_initial_stress_fcn
	initial stress components More...
MaterialProperty< RankFourTensor > &	_Jacobian_mult
	derivative of stress w.r.t. strain (_dstress_dstrain) More...
const unsigned int	_max_iterations
	Input parameters associated with the recompute iteration to return the stress state to the yield surface. More...
const Real	_relative_tolerance
const Real	_absolute_tolerance
const bool	_internal_solve_full_iteration_history
const MaterialProperty< RankTwoTensor > &	_elastic_strain_old
	Strain tensors. More...
const MaterialProperty< RankTwoTensor > &	_strain_increment

Private Attributes
const InputParameters &	_gc_params
	Parameters of the object with this interface. More...
FEProblemBase *const	_gc_feproblem
	Reference to the FEProblemBase class. More...
BlockRestrictable *const	_gc_block_restrict
	Access block restrictions of the object with this interface. More...

Detailed Description

ComputeMultipleInelasticStress computes the stress, the consistent tangent operator (or an approximation to it), and a decomposition of the strain into elastic and inelastic parts.

By default finite strains are assumed.

The elastic strain is calculated by subtracting the computed inelastic strain increment tensor from the mechanical strain tensor. Mechanical strain is considered as the sum of the elastic and inelastic (plastic, creep, ect) strains.

This material is used to call the recompute iterative materials of a number of specified inelastic models that inherit from StressUpdateBase. It iterates over the specified inelastic models until the change in stress is within a user-specified tolerance, in order to produce the stress, the consistent tangent operator and the elastic and inelastic strains for the time increment.

Definition at line 38 of file ComputeMultipleInelasticStress.h.

Member Enumeration Documentation

TangentOperatorEnum

enum ComputeMultipleInelasticStress::TangentOperatorEnum

strongprotected

what sort of Tangent operator to calculate

Enumerator
elastic
nonlinear

Definition at line 143 of file ComputeMultipleInelasticStress.h.

{ elastic, nonlinear } _tangent_operator_type;

Constructor & Destructor Documentation

ComputeMultipleInelasticStress()

ComputeMultipleInelasticStress::ComputeMultipleInelasticStress

(

const InputParameters &

parameters

)

Definition at line 77 of file ComputeMultipleInelasticStress.C.

  : ComputeFiniteStrainElasticStress(parameters),
    _max_iterations(parameters.get<unsigned int>("max_iterations")),
    _relative_tolerance(parameters.get<Real>("relative_tolerance")),
    _absolute_tolerance(parameters.get<Real>("absolute_tolerance")),
    _internal_solve_full_iteration_history(getParam<bool>("internal_solve_full_iteration_history")),
    _perform_finite_strain_rotations(getParam<bool>("perform_finite_strain_rotations")),
    _elastic_strain_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "elastic_strain")),
    _strain_increment(getMaterialProperty<RankTwoTensor>(_base_name + "strain_increment")),
    _inelastic_strain(declareProperty<RankTwoTensor>(_base_name + "combined_inelastic_strain")),
    _inelastic_strain_old(
        getMaterialPropertyOld<RankTwoTensor>(_base_name + "combined_inelastic_strain")),
    _tangent_operator_type(getParam<MooseEnum>("tangent_operator").getEnum<TangentOperatorEnum>()),
    _num_models(getParam<std::vector<MaterialName>>("inelastic_models").size()),
    _tangent_computation_flag(_num_models, false),
    _tangent_calculation_method(TangentCalculationMethod::ELASTIC),
    _inelastic_weights(isParamValid("combined_inelastic_strain_weights")
                           ? getParam<std::vector<Real>>("combined_inelastic_strain_weights")
                           : std::vector<Real>(_num_models, true)),
    _consistent_tangent_operator(_num_models),
    _cycle_models(getParam<bool>("cycle_models")),
    _matl_timestep_limit(declareProperty<Real>("matl_timestep_limit")),
    _identity_symmetric_four(RankFourTensor::initIdentitySymmetricFour),
    _all_models_isotropic(true),
    _damage_model(isParamValid("damage_model")
                      ? dynamic_cast<DamageBase *>(&getMaterial("damage_model"))
                      : nullptr)
{
  if (_inelastic_weights.size() != _num_models)
    mooseError(
        "ComputeMultipleInelasticStress: combined_inelastic_strain_weights must contain the same "
        "number of entries as inelastic_models ",
        _inelastic_weights.size(),
        " ",
        _num_models);
}

Member Function Documentation

computeAdmissibleState()

void ComputeMultipleInelasticStress::computeAdmissibleState ( unsigned  model_number, RankTwoTensor &  elastic_strain_increment, RankTwoTensor &  inelastic_strain_increment, RankFourTensor &  consistent_tangent_operator  )

protectedvirtual

Given a trial stress (_stress[_qp]) and a strain increment (elastic_strain_increment) let the model_number model produce an admissible stress (gets placed back in _stress[_qp]), and decompose the strain increment into an elastic part (gets placed back into elastic_strain_increment) and an inelastic part (inelastic_strain_increment), as well as computing the consistent_tangent_operator.

Parameters

model_number	The inelastic model to use
elastic_strain_increment	Upon input, this is the strain increment. Upon output, it is the elastic part of the strain increment
inelastic_strain_increment	The inelastic strain increment corresponding to the supplied strain increment
consistent_tangent_operator	The consistent tangent operator

Reimplemented in ComputeMultipleInelasticCosseratStress.

Definition at line 466 of file ComputeMultipleInelasticStress.C.

{
  const bool jac = _fe_problem.currentlyComputingJacobian();
  if (_damage_model)
    _models[model_number]->updateState(elastic_strain_increment,
                                       inelastic_strain_increment,
                                       _rotation_increment[_qp],
                                       _stress[_qp],
                                       _undamaged_stress_old,
                                       _elasticity_tensor[_qp],
                                       _elastic_strain_old[_qp],
                                       (jac && _tangent_computation_flag[model_number]),
                                       consistent_tangent_operator);
  else
    _models[model_number]->updateState(elastic_strain_increment,
                                       inelastic_strain_increment,
                                       _rotation_increment[_qp],
                                       _stress[_qp],
                                       _stress_old[_qp],
                                       _elasticity_tensor[_qp],
                                       _elastic_strain_old[_qp],
                                       (jac && _tangent_computation_flag[model_number]),
                                       consistent_tangent_operator);
 
  if (jac && !_tangent_computation_flag[model_number])
  {
    if (_tangent_calculation_method == TangentCalculationMethod::PARTIAL)
      consistent_tangent_operator.zero();
    else
      consistent_tangent_operator = _elasticity_tensor[_qp];
  }
}

Referenced by ComputeMultipleInelasticCosseratStress::computeAdmissibleState(), updateQpState(), and updateQpStateSingleModel().

computeQpJacobianMult()

void ComputeMultipleInelasticStress::computeQpJacobianMult ( )

protectedvirtual

Using _elasticity_tensor[_qp] and the consistent tangent operators, _consistent_tangent_operator[...] computed by the inelastic models, compute _Jacobian_mult[_qp].

Reimplemented in ComputeMultipleInelasticCosseratStress.

Definition at line 393 of file ComputeMultipleInelasticStress.C.

{
  if (_tangent_calculation_method == TangentCalculationMethod::ELASTIC)
    _Jacobian_mult[_qp] = _elasticity_tensor[_qp];
  else if (_tangent_calculation_method == TangentCalculationMethod::PARTIAL)
  {
    RankFourTensor A = _identity_symmetric_four;
    for (unsigned i_rmm = 0; i_rmm < _num_models; ++i_rmm)
      A += _consistent_tangent_operator[i_rmm];
    mooseAssert(A.isSymmetric(), "Tangent operator isn't symmetric");
    _Jacobian_mult[_qp] = A.invSymm() * _elasticity_tensor[_qp];
  }
  else
  {
    const RankFourTensor E_inv = _elasticity_tensor[_qp].invSymm();
    _Jacobian_mult[_qp] = _consistent_tangent_operator[0];
    for (unsigned i_rmm = 1; i_rmm < _num_models; ++i_rmm)
      _Jacobian_mult[_qp] = _consistent_tangent_operator[i_rmm] * E_inv * _Jacobian_mult[_qp];
  }
}

Referenced by updateQpState().

computeQpProperties()

void ComputeStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

Definition at line 50 of file ComputeStressBase.C.

{
  computeQpStress();
 
  // Add in extra stress
  _stress[_qp] += _extra_stress[_qp];
}

computeQpStress()

void ComputeMultipleInelasticStress::computeQpStress ( )

overrideprotectedvirtual

Compute the stress and store it in the _stress material property for the current quadrature point.

Reimplemented from ComputeFiniteStrainElasticStress.

Reimplemented in ComputeMultipleInelasticCosseratStress, and ComputeSmearedCrackingStress.

Definition at line 187 of file ComputeMultipleInelasticStress.C.

{
  if (_damage_model)
  {
    _undamaged_stress_old = _stress_old[_qp];
    _damage_model->setQp(_qp);
    _damage_model->computeUndamagedOldStress(_undamaged_stress_old);
  }
  computeQpStressIntermediateConfiguration();
 
  if (_damage_model)
  {
    _damage_model->setQp(_qp);
    _damage_model->updateDamage();
    _damage_model->updateStressForDamage(_stress[_qp]);
    _damage_model->finiteStrainRotation(_rotation_increment[_qp]);
    _damage_model->updateJacobianMultForDamage(_Jacobian_mult[_qp]);
 
    const Real damage_timestep_limit = _damage_model->computeTimeStepLimit();
    if (_matl_timestep_limit[_qp] > damage_timestep_limit)
      _matl_timestep_limit[_qp] = damage_timestep_limit;
  }
 
  if (_perform_finite_strain_rotations)
    finiteStrainRotation();
}

Referenced by ComputeMultipleInelasticCosseratStress::computeQpStress().

computeQpStressIntermediateConfiguration()

void ComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration ( )

protectedvirtual

Compute the stress for the current QP, but do not rotate tensors from the intermediate configuration to the new configuration.

Definition at line 215 of file ComputeMultipleInelasticStress.C.

{
  RankTwoTensor elastic_strain_increment;
  RankTwoTensor combined_inelastic_strain_increment;
 
  if (_num_models == 0)
  {
    _elastic_strain[_qp] = _elastic_strain_old[_qp] + _strain_increment[_qp];
 
    // If the elasticity tensor values have changed and the tensor is isotropic,
    // use the old strain to calculate the old stress
    if (_is_elasticity_tensor_guaranteed_isotropic || !_perform_finite_strain_rotations)
      _stress[_qp] = _elasticity_tensor[_qp] * (_elastic_strain_old[_qp] + _strain_increment[_qp]);
    else
    {
      if (_damage_model)
        _stress[_qp] = _undamaged_stress_old + _elasticity_tensor[_qp] * _strain_increment[_qp];
      else
        _stress[_qp] = _stress_old[_qp] + _elasticity_tensor[_qp] * _strain_increment[_qp];
    }
    if (_fe_problem.currentlyComputingJacobian())
      _Jacobian_mult[_qp] = _elasticity_tensor[_qp];
 
    _matl_timestep_limit[_qp] = std::numeric_limits<Real>::max();
  }
  else
  {
    if (_num_models == 1 || _cycle_models)
      updateQpStateSingleModel((_t_step - 1) % _num_models,
                               elastic_strain_increment,
                               combined_inelastic_strain_increment);
    else
      updateQpState(elastic_strain_increment, combined_inelastic_strain_increment);
 
    _elastic_strain[_qp] = _elastic_strain_old[_qp] + elastic_strain_increment;
    _inelastic_strain[_qp] = _inelastic_strain_old[_qp] + combined_inelastic_strain_increment;
  }
}

Referenced by ComputeSmearedCrackingStress::computeQpStress(), and computeQpStress().

finiteStrainRotation()

void ComputeMultipleInelasticStress::finiteStrainRotation ( const bool  force_elasticity_rotation = false )

protectedvirtual

Rotate _elastic_strain, _stress, _inelastic_strain, and _Jacobian_mult to the new configuration.

Parameters

force_elasticity_rotation

Force the elasticity tensor to be rotated, even if it is not deemed necessary.

Definition at line 255 of file ComputeMultipleInelasticStress.C.

{
  _elastic_strain[_qp] =
      _rotation_increment[_qp] * _elastic_strain[_qp] * _rotation_increment[_qp].transpose();
  _stress[_qp] = _rotation_increment[_qp] * _stress[_qp] * _rotation_increment[_qp].transpose();
  _inelastic_strain[_qp] =
      _rotation_increment[_qp] * _inelastic_strain[_qp] * _rotation_increment[_qp].transpose();
 
  if (force_elasticity_rotation ||
      !(_is_elasticity_tensor_guaranteed_isotropic &&
        (_all_models_isotropic ||
         _tangent_calculation_method == TangentCalculationMethod::ELASTIC || _num_models == 0)))
    _Jacobian_mult[_qp].rotate(_rotation_increment[_qp]);
}

Referenced by ComputeSmearedCrackingStress::computeQpStress(), and computeQpStress().

hasGuaranteedMaterialProperty()

bool GuaranteeConsumer::hasGuaranteedMaterialProperty ( const MaterialPropertyName &  prop, Guarantee  guarantee  )

protectedinherited

Definition at line 28 of file GuaranteeConsumer.C.

{
  if (!_gc_feproblem->startedInitialSetup())
    mooseError("hasGuaranteedMaterialProperty() needs to be called in initialSetup()");
 
  // Reference to MaterialWarehouse for testing and retrieving block ids
  const auto & warehouse = _gc_feproblem->getMaterialWarehouse();
 
  // Complete set of ids that this object is active
  const auto & ids = (_gc_block_restrict && _gc_block_restrict->blockRestricted())
                         ? _gc_block_restrict->blockIDs()
                         : _gc_feproblem->mesh().meshSubdomains();
 
  // Loop over each id for this object
  for (const auto & id : ids)
  {
    // If block materials exist, look if any issue the required guarantee
    if (warehouse.hasActiveBlockObjects(id))
    {
      const std::vector<std::shared_ptr<MaterialBase>> & mats = warehouse.getActiveBlockObjects(id);
      for (const auto & mat : mats)
      {
        const auto & mat_props = mat->getSuppliedItems();
        if (mat_props.count(prop_name))
        {
          auto guarantee_mat = dynamic_cast<GuaranteeProvider *>(mat.get());
          if (guarantee_mat && !guarantee_mat->hasGuarantee(prop_name, guarantee))
          {
            // we found at least one material on the set of block we operate on
            // that does _not_ provide the requested guarantee
            return false;
          }
        }
      }
    }
  }
 
  return true;
}

Referenced by ComputeFiniteStrainElasticStress::initialSetup(), ComputeSmearedCrackingStress::initialSetup(), ComputeLinearElasticPFFractureStress::initialSetup(), CriticalTimeStep::initialSetup(), and initialSetup().

initialSetup()

void ComputeMultipleInelasticStress::initialSetup ( )

overridevirtual

Reimplemented in ComputeSmearedCrackingStress.

Definition at line 122 of file ComputeMultipleInelasticStress.C.

{
  _is_elasticity_tensor_guaranteed_isotropic =
      hasGuaranteedMaterialProperty(_elasticity_tensor_name, Guarantee::ISOTROPIC);
 
  std::vector<MaterialName> models = getParam<std::vector<MaterialName>>("inelastic_models");
 
  for (unsigned int i = 0; i < _num_models; ++i)
  {
    StressUpdateBase * rrr = dynamic_cast<StressUpdateBase *>(&getMaterialByName(models[i]));
 
    if (rrr)
    {
      _models.push_back(rrr);
      _all_models_isotropic = _all_models_isotropic && rrr->isIsotropic();
      if (rrr->requiresIsotropicTensor() && !_is_elasticity_tensor_guaranteed_isotropic)
        mooseError("Model " + models[i] +
                   " requires an isotropic elasticity tensor, but the one supplied is not "
                   "guaranteed isotropic");
    }
    else
      mooseError("Model " + models[i] + " is not compatible with ComputeMultipleInelasticStress");
  }
 
  // Check if tangent calculation methods are consistent. If all models have
  // TangentOperatorEnum::ELASTIC or tangent_operator is set by the user as elasic, then the tangent
  // is never calculated: J_tot = C. If PARTIAL and NONE models are present, utilize PARTIAL
  // formulation: J_tot = (I + J_1 + ... J_N)^-1 C. If FULL and NONE models are present, utilize
  // FULL formulation: J_tot = J_1 * C^-1 * J_2 * C^-1 * ... J_N * C. If PARTIAL and FULL models are
  // present, error out.
 
  if (_tangent_operator_type != TangentOperatorEnum::elastic)
  {
    bool full_present = false;
    bool partial_present = false;
    for (unsigned int i = 0; i < _num_models; ++i)
    {
      if (_models[i]->getTangentCalculationMethod() == TangentCalculationMethod::FULL)
      {
        full_present = true;
        _tangent_computation_flag[i] = true;
        _tangent_calculation_method = TangentCalculationMethod::FULL;
      }
      else if (_models[i]->getTangentCalculationMethod() == TangentCalculationMethod::PARTIAL)
      {
        partial_present = true;
        _tangent_computation_flag[i] = true;
        _tangent_calculation_method = TangentCalculationMethod::PARTIAL;
      }
    }
    if (full_present && partial_present)
      mooseError("In ",
                 _name,
                 ": Models that calculate the full tangent operator and the partial tangent "
                 "operator are being combined. Either set tangent_operator to elastic, implement "
                 "the corrent tangent formulations, or use different models.");
  }
 
  if (isParamValid("damage_model") && !_damage_model)
    paramError("damage_model",
               "Damage Model " + _damage_model->name() +
                   " is not compatible with ComputeMultipleInelasticStress");
}

Referenced by ComputeSmearedCrackingStress::initialSetup().

initQpStatefulProperties()

void ComputeMultipleInelasticStress::initQpStatefulProperties ( )

overrideprotectedvirtual

Reimplemented from ComputeStressBase.

Reimplemented in ComputeMultipleInelasticCosseratStress, and ComputeSmearedCrackingStress.

Definition at line 115 of file ComputeMultipleInelasticStress.C.

{
  ComputeStressBase::initQpStatefulProperties();
  _inelastic_strain[_qp].zero();
}

Referenced by ComputeSmearedCrackingStress::initQpStatefulProperties(), and ComputeMultipleInelasticCosseratStress::initQpStatefulProperties().

updateQpState()

void ComputeMultipleInelasticStress::updateQpState ( RankTwoTensor &  elastic_strain_increment, RankTwoTensor &  combined_inelastic_strain_increment  )

protectedvirtual

Given the _strain_increment[_qp], iterate over all of the user-specified recompute materials in order to find an admissible stress (which is placed into _stress[_qp]) and set of inelastic strains, as well as the tangent operator (which is placed into _Jacobian_mult[_qp]).

Parameters

elastic_strain_increment	The elastic part of _strain_increment[_qp] after the iterative process has converged
combined_inelastic_strain_increment	The inelastic part of _strain_increment[_qp] after the iterative process has converged. This is a weighted sum of all the inelastic strains computed by all the plastic models during the Picard iterative scheme. The weights are dictated by the user using _inelastic_weights

Definition at line 271 of file ComputeMultipleInelasticStress.C.

{
  if (_internal_solve_full_iteration_history == true)
  {
    _console << std::endl
             << "iteration output for ComputeMultipleInelasticStress solve:"
             << " time=" << _t << " int_pt=" << _qp << std::endl;
  }
  Real l2norm_delta_stress;
  Real first_l2norm_delta_stress = 1.0;
  unsigned int counter = 0;
 
  std::vector<RankTwoTensor> inelastic_strain_increment;
  inelastic_strain_increment.resize(_num_models);
 
  for (unsigned i_rmm = 0; i_rmm < _models.size(); ++i_rmm)
    inelastic_strain_increment[i_rmm].zero();
 
  RankTwoTensor stress_max, stress_min;
 
  do
  {
    for (unsigned i_rmm = 0; i_rmm < _num_models; ++i_rmm)
    {
      _models[i_rmm]->setQp(_qp);
 
      // initially assume the strain is completely elastic
      elastic_strain_increment = _strain_increment[_qp];
      // and subtract off all inelastic strain increments calculated so far
      // except the one that we're about to calculate
      for (unsigned j_rmm = 0; j_rmm < _num_models; ++j_rmm)
        if (i_rmm != j_rmm)
          elastic_strain_increment -= inelastic_strain_increment[j_rmm];
 
      // form the trial stress, with the check for changed elasticity constants
      if (_is_elasticity_tensor_guaranteed_isotropic || !_perform_finite_strain_rotations)
        _stress[_qp] =
            _elasticity_tensor[_qp] * (_elastic_strain_old[_qp] + elastic_strain_increment);
      else
      {
        if (_damage_model)
          _stress[_qp] = _undamaged_stress_old + _elasticity_tensor[_qp] * elastic_strain_increment;
        else
          _stress[_qp] = _stress_old[_qp] + _elasticity_tensor[_qp] * elastic_strain_increment;
      }
 
      // given a trial stress (_stress[_qp]) and a strain increment (elastic_strain_increment)
      // let the i^th model produce an admissible stress (as _stress[_qp]), and decompose
      // the strain increment into an elastic part (elastic_strain_increment) and an
      // inelastic part (inelastic_strain_increment[i_rmm])
      computeAdmissibleState(i_rmm,
                             elastic_strain_increment,
                             inelastic_strain_increment[i_rmm],
                             _consistent_tangent_operator[i_rmm]);
 
      if (i_rmm == 0)
      {
        stress_max = _stress[_qp];
        stress_min = _stress[_qp];
      }
      else
      {
        for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
        {
          for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
          {
            if (_stress[_qp](i, j) > stress_max(i, j))
              stress_max(i, j) = _stress[_qp](i, j);
            else if (stress_min(i, j) > _stress[_qp](i, j))
              stress_min(i, j) = _stress[_qp](i, j);
          }
        }
      }
    }
 
    // now check convergence in the stress:
    // once the change in stress is within tolerance after each recompute material
    // consider the stress to be converged
    l2norm_delta_stress = (stress_max - stress_min).L2norm();
    if (counter == 0 && l2norm_delta_stress > 0.0)
      first_l2norm_delta_stress = l2norm_delta_stress;
 
    if (_internal_solve_full_iteration_history == true)
    {
      _console << "stress iteration number = " << counter << "\n"
               << " relative l2 norm delta stress = "
               << (0 == first_l2norm_delta_stress ? 0
                                                  : l2norm_delta_stress / first_l2norm_delta_stress)
               << "\n"
               << " stress convergence relative tolerance = " << _relative_tolerance << "\n"
               << " absolute l2 norm delta stress = " << l2norm_delta_stress << "\n"
               << " stress convergence absolute tolerance = " << _absolute_tolerance << std::endl;
    }
    ++counter;
  } while (counter < _max_iterations && l2norm_delta_stress > _absolute_tolerance &&
           (l2norm_delta_stress / first_l2norm_delta_stress) > _relative_tolerance &&
           _num_models != 1);
 
  if (counter == _max_iterations && l2norm_delta_stress > _absolute_tolerance &&
      (l2norm_delta_stress / first_l2norm_delta_stress) > _relative_tolerance)
    throw MooseException("Max stress iteration hit during ComputeMultipleInelasticStress solve!");
 
  combined_inelastic_strain_increment.zero();
  for (unsigned i_rmm = 0; i_rmm < _num_models; ++i_rmm)
    combined_inelastic_strain_increment +=
        _inelastic_weights[i_rmm] * inelastic_strain_increment[i_rmm];
 
  if (_fe_problem.currentlyComputingJacobian())
    computeQpJacobianMult();
 
  _matl_timestep_limit[_qp] = 0.0;
  for (unsigned i_rmm = 0; i_rmm < _num_models; ++i_rmm)
    _matl_timestep_limit[_qp] += 1.0 / _models[i_rmm]->computeTimeStepLimit();
 
  if (MooseUtils::absoluteFuzzyEqual(_matl_timestep_limit[_qp], 0.0))
    _matl_timestep_limit[_qp] = std::numeric_limits<Real>::max();
  else
    _matl_timestep_limit[_qp] = 1.0 / _matl_timestep_limit[_qp];
}

Referenced by computeQpStressIntermediateConfiguration().

updateQpStateSingleModel()

void ComputeMultipleInelasticStress::updateQpStateSingleModel ( unsigned  model_number, RankTwoTensor &  elastic_strain_increment, RankTwoTensor &  combined_inelastic_strain_increment  )

protectedvirtual

An optimised version of updateQpState that gets used when the number of plastic models is unity, or when we're cycling through models Given the _strain_increment[_qp], find an admissible stress (which is put into _stress[_qp]) and inelastic strain, as well as the tangent operator (which is placed into _Jacobian_mult[_qp])

Parameters

model_number	Use this model number
elastic_strain_increment	The elastic part of _strain_increment[_qp]
combined_inelastic_strain_increment	The inelastic part of _strain_increment[_qp]

Definition at line 415 of file ComputeMultipleInelasticStress.C.

{
  for (auto model : _models)
    model->setQp(_qp);
 
  elastic_strain_increment = _strain_increment[_qp];
 
  // If the elasticity tensor values have changed and the tensor is isotropic,
  // use the old strain to calculate the old stress
  if (_is_elasticity_tensor_guaranteed_isotropic || !_perform_finite_strain_rotations)
    _stress[_qp] = _elasticity_tensor[_qp] * (_elastic_strain_old[_qp] + elastic_strain_increment);
  else
  {
    if (_damage_model)
      _stress[_qp] = _undamaged_stress_old + _elasticity_tensor[_qp] * elastic_strain_increment;
    else
      _stress[_qp] = _stress_old[_qp] + _elasticity_tensor[_qp] * elastic_strain_increment;
  }
 
  computeAdmissibleState(model_number,
                         elastic_strain_increment,
                         combined_inelastic_strain_increment,
                         _consistent_tangent_operator[0]);
 
  if (_fe_problem.currentlyComputingJacobian())
  {
    if (_tangent_calculation_method == TangentCalculationMethod::ELASTIC)
      _Jacobian_mult[_qp] = _elasticity_tensor[_qp];
    else if (_tangent_calculation_method == TangentCalculationMethod::PARTIAL)
    {
      RankFourTensor A = _identity_symmetric_four + _consistent_tangent_operator[0];
      mooseAssert(A.isSymmetric(), "Tangent operator isn't symmetric");
      _Jacobian_mult[_qp] = A.invSymm() * _elasticity_tensor[_qp];
    }
    else
      _Jacobian_mult[_qp] = _consistent_tangent_operator[0];
  }
 
  _matl_timestep_limit[_qp] = _models[0]->computeTimeStepLimit();
 
  /* propagate internal variables, etc, to this timestep for those inelastic models where
   * "updateState" is not called */
  for (unsigned i_rmm = 0; i_rmm < _num_models; ++i_rmm)
    if (i_rmm != model_number)
      _models[i_rmm]->propagateQpStatefulProperties();
}

Referenced by computeQpStressIntermediateConfiguration().

validParams()

InputParameters ComputeMultipleInelasticStress::validParams ( )

static

Definition at line 21 of file ComputeMultipleInelasticStress.C.

{
  InputParameters params = ComputeFiniteStrainElasticStress::validParams();
  params.addClassDescription("Compute state (stress and internal parameters such as plastic "
                             "strains and internal parameters) using an iterative process.  "
                             "Combinations of creep models and plastic models may be used.");
  params.addParam<unsigned int>("max_iterations",
                                30,
                                "Maximum number of the stress update "
                                "iterations over the stress change after all "
                                "update materials are called");
  params.addParam<Real>("relative_tolerance",
                        1e-5,
                        "Relative convergence tolerance for the stress "
                        "update iterations over the stress change "
                        "after all update materials are called");
  params.addParam<Real>("absolute_tolerance",
                        1e-5,
                        "Absolute convergence tolerance for the stress "
                        "update iterations over the stress change "
                        "after all update materials are called");
  params.addParam<bool>(
      "internal_solve_full_iteration_history",
      false,
      "Set to true to output stress update iteration information over the stress change");
  params.addParam<bool>("perform_finite_strain_rotations",
                        true,
                        "Tensors are correctly rotated in "
                        "finite-strain simulations.  For "
                        "optimal performance you can set "
                        "this to 'false' if you are only "
                        "ever using small strains");
  MooseEnum tangent_operator("elastic nonlinear", "nonlinear");
  params.addParam<MooseEnum>(
      "tangent_operator",
      tangent_operator,
      "Type of tangent operator to return.  'elastic': return the "
      "elasticity tensor.  'nonlinear': return the full, general consistent tangent "
      "operator.");
  params.addRequiredParam<std::vector<MaterialName>>(
      "inelastic_models",
      "The material objects to use to calculate stress and inelastic strains. "
      "Note: specify creep models first and plasticity models second.");
  params.addParam<std::vector<Real>>("combined_inelastic_strain_weights",
                                     "The combined_inelastic_strain Material Property is a "
                                     "weighted sum of the model inelastic strains.  This parameter "
                                     "is a vector of weights, of the same length as "
                                     "inelastic_models.  Default = '1 1 ... 1'.  This "
                                     "parameter is set to 1 if the number of models = 1");
  params.addParam<bool>(
      "cycle_models", false, "At timestep N use only inelastic model N % num_models.");
  params.addParam<MaterialName>("damage_model", "Name of the damage model");
 
  return params;
}

Referenced by ComputeSmearedCrackingStress::validParams(), and ComputeMultipleInelasticCosseratStress::validParams().

Member Data Documentation

_absolute_tolerance

const Real ComputeMultipleInelasticStress::_absolute_tolerance

protected

Definition at line 124 of file ComputeMultipleInelasticStress.h.

Referenced by updateQpState().

_all_models_isotropic

bool ComputeMultipleInelasticStress::_all_models_isotropic

protected

are all inelastic models inherently isotropic? (not the case for e.g. weak plane plasticity models)

Definition at line 183 of file ComputeMultipleInelasticStress.h.

Referenced by finiteStrainRotation(), and initialSetup().

_base_name

const std::string ComputeStressBase::_base_name

protectedinherited

Base name prepended to all material property names to allow for multi-material systems.

Definition at line 45 of file ComputeStressBase.h.

Referenced by ComputeLinearElasticStress::initialSetup(), and ComputeCosseratLinearElasticStress::initialSetup().

_consistent_tangent_operator

std::vector<RankFourTensor> ComputeMultipleInelasticStress::_consistent_tangent_operator

protected

the consistent tangent operators computed by each plastic model

Definition at line 158 of file ComputeMultipleInelasticStress.h.

Referenced by ComputeMultipleInelasticCosseratStress::computeQpJacobianMult(), computeQpJacobianMult(), updateQpState(), and updateQpStateSingleModel().

_cycle_models

const bool ComputeMultipleInelasticStress::_cycle_models

protected

whether to cycle through the models, using only one model per timestep

Definition at line 161 of file ComputeMultipleInelasticStress.h.

Referenced by computeQpStressIntermediateConfiguration().

_damage_model

DamageBase* ComputeMultipleInelasticStress::_damage_model

protected

Pointer to the damage model.

Definition at line 186 of file ComputeMultipleInelasticStress.h.

Referenced by computeAdmissibleState(), computeQpStress(), computeQpStressIntermediateConfiguration(), initialSetup(), updateQpState(), and updateQpStateSingleModel().

_elastic_strain

MaterialProperty<RankTwoTensor>& ComputeStressBase::_elastic_strain

protectedinherited

Elastic strain material property.

Definition at line 52 of file ComputeStressBase.h.

Referenced by ComputeSmearedCrackingStress::computeCrackStrainAndOrientation(), ComputeLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeMultiPlasticityStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), computeQpStressIntermediateConfiguration(), finiteStrainRotation(), and ComputeStressBase::initQpStatefulProperties().

_elastic_strain_old

const MaterialProperty<RankTwoTensor>& ComputeMultipleInelasticStress::_elastic_strain_old

protected

Strain tensors.

Definition at line 132 of file ComputeMultipleInelasticStress.h.

Referenced by computeAdmissibleState(), ComputeSmearedCrackingStress::computeQpStress(), computeQpStressIntermediateConfiguration(), ComputeSmearedCrackingStress::updateLocalElasticityTensor(), updateQpState(), and updateQpStateSingleModel().

_elasticity_tensor

const MaterialProperty<RankFourTensor>& ComputeFiniteStrainElasticStress::_elasticity_tensor

protectedinherited

Elasticity tensor material property.

Definition at line 39 of file ComputeFiniteStrainElasticStress.h.

Referenced by computeAdmissibleState(), ComputeMultipleInelasticCosseratStress::computeQpJacobianMult(), computeQpJacobianMult(), ComputeFiniteStrainElasticStress::computeQpStress(), computeQpStressIntermediateConfiguration(), ComputeSmearedCrackingStress::updateCrackingStateAndStress(), ComputeSmearedCrackingStress::updateLocalElasticityTensor(), updateQpState(), and updateQpStateSingleModel().

_elasticity_tensor_name

const std::string ComputeFiniteStrainElasticStress::_elasticity_tensor_name

protectedinherited

Name of the elasticity tensor material property.

Definition at line 37 of file ComputeFiniteStrainElasticStress.h.

Referenced by ComputeFiniteStrainElasticStress::initialSetup(), ComputeSmearedCrackingStress::initialSetup(), and initialSetup().

_extra_stress

const MaterialProperty<RankTwoTensor>& ComputeStressBase::_extra_stress

protectedinherited

Extra stress tensor.

Definition at line 55 of file ComputeStressBase.h.

Referenced by ComputeStressBase::computeQpProperties().

_gc_block_restrict

BlockRestrictable* const GuaranteeConsumer::_gc_block_restrict

privateinherited

Access block restrictions of the object with this interface.

Definition at line 41 of file GuaranteeConsumer.h.

Referenced by GuaranteeConsumer::hasGuaranteedMaterialProperty().

_gc_feproblem

FEProblemBase* const GuaranteeConsumer::_gc_feproblem

privateinherited

Reference to the FEProblemBase class.

Definition at line 38 of file GuaranteeConsumer.h.

Referenced by GuaranteeConsumer::hasGuaranteedMaterialProperty().

_gc_params

const InputParameters& GuaranteeConsumer::_gc_params

privateinherited

Parameters of the object with this interface.

Definition at line 35 of file GuaranteeConsumer.h.

_identity_symmetric_four

const RankFourTensor ComputeMultipleInelasticStress::_identity_symmetric_four

protected

Rank four symmetric identity tensor.

Definition at line 168 of file ComputeMultipleInelasticStress.h.

Referenced by computeQpJacobianMult(), and updateQpStateSingleModel().

_inelastic_strain

MaterialProperty<RankTwoTensor>& ComputeMultipleInelasticStress::_inelastic_strain

protected

The sum of the inelastic strains that come from the plastic models.

Definition at line 137 of file ComputeMultipleInelasticStress.h.

Referenced by ComputeSmearedCrackingStress::computeQpStress(), computeQpStressIntermediateConfiguration(), finiteStrainRotation(), and initQpStatefulProperties().

_inelastic_strain_old

const MaterialProperty<RankTwoTensor>& ComputeMultipleInelasticStress::_inelastic_strain_old

protected

old value of inelastic strain

Definition at line 140 of file ComputeMultipleInelasticStress.h.

Referenced by ComputeSmearedCrackingStress::computeQpStress(), and computeQpStressIntermediateConfiguration().

_inelastic_weights

const std::vector<Real> ComputeMultipleInelasticStress::_inelastic_weights

protected

_inelastic_strain = sum_i (_inelastic_weights_i * inelastic_strain_from_model_i)

Definition at line 155 of file ComputeMultipleInelasticStress.h.

Referenced by ComputeMultipleInelasticStress(), and updateQpState().

_initial_stress_fcn

std::vector<const Function *> ComputeStressBase::_initial_stress_fcn

protectedinherited

initial stress components

Definition at line 58 of file ComputeStressBase.h.

_internal_solve_full_iteration_history

const bool ComputeMultipleInelasticStress::_internal_solve_full_iteration_history

protected

Definition at line 125 of file ComputeMultipleInelasticStress.h.

Referenced by updateQpState().

_is_elasticity_tensor_guaranteed_isotropic

bool ComputeMultipleInelasticStress::_is_elasticity_tensor_guaranteed_isotropic

protected

is the elasticity tensor guaranteed to be isotropic?

Definition at line 180 of file ComputeMultipleInelasticStress.h.

Referenced by computeQpStressIntermediateConfiguration(), finiteStrainRotation(), initialSetup(), updateQpState(), and updateQpStateSingleModel().

_Jacobian_mult

MaterialProperty<RankFourTensor>& ComputeStressBase::_Jacobian_mult

protectedinherited

derivative of stress w.r.t. strain (_dstress_dstrain)

Definition at line 61 of file ComputeStressBase.h.

Referenced by ComputeStrainIncrementBasedStress::computeQpJacobian(), FiniteStrainHyperElasticViscoPlastic::computeQpJacobian(), ComputeMultipleInelasticCosseratStress::computeQpJacobianMult(), computeQpJacobianMult(), ComputeLinearElasticStress::computeQpStress(), ComputeDamageStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeLinearElasticPFFractureStress::computeQpStress(), ComputeMultiPlasticityStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), computeQpStress(), computeQpStressIntermediateConfiguration(), ComputeLinearElasticPFFractureStress::computeStrainSpectral(), ComputeLinearElasticPFFractureStress::computeStrainVolDev(), ComputeLinearElasticPFFractureStress::computeStressSpectral(), FiniteStrainUObasedCP::elasticTangentModuli(), FiniteStrainUObasedCP::elastoPlasticTangentModuli(), finiteStrainRotation(), ComputeMultiPlasticityStress::postReturnMap(), FiniteStrainCrystalPlasticity::postSolveQp(), FiniteStrainCrystalPlasticity::preSolveQp(), and updateQpStateSingleModel().

_matl_timestep_limit

MaterialProperty<Real>& ComputeMultipleInelasticStress::_matl_timestep_limit

protected

Definition at line 163 of file ComputeMultipleInelasticStress.h.

Referenced by ComputeSmearedCrackingStress::computeQpStress(), computeQpStress(), computeQpStressIntermediateConfiguration(), updateQpState(), and updateQpStateSingleModel().

_max_iterations

const unsigned int ComputeMultipleInelasticStress::_max_iterations

protected

Input parameters associated with the recompute iteration to return the stress state to the yield surface.

Definition at line 122 of file ComputeMultipleInelasticStress.h.

Referenced by updateQpState().

_mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeStressBase::_mechanical_strain

protectedinherited

Mechanical strain material property.

Definition at line 48 of file ComputeStressBase.h.

Referenced by ComputeLinearElasticStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeLinearElasticPFFractureStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), ComputeLinearElasticPFFractureStress::computeStrainSpectral(), ComputeLinearElasticPFFractureStress::computeStrainVolDev(), and ComputeLinearElasticPFFractureStress::computeStressSpectral().

_models

std::vector<StressUpdateBase *> ComputeMultipleInelasticStress::_models

protected

The user supplied list of inelastic models to use in the simulation.

Users should take care to list creep models first and plasticity models last to allow for the case when a creep model relaxes the stress state inside of the yield surface in an iteration.

Definition at line 177 of file ComputeMultipleInelasticStress.h.

Referenced by computeAdmissibleState(), ComputeSmearedCrackingStress::computeQpStress(), initialSetup(), updateQpState(), and updateQpStateSingleModel().

_num_models

const unsigned ComputeMultipleInelasticStress::_num_models

protected

number of plastic models

Definition at line 146 of file ComputeMultipleInelasticStress.h.

Referenced by ComputeMultipleInelasticStress(), ComputeMultipleInelasticCosseratStress::computeQpJacobianMult(), computeQpJacobianMult(), computeQpStressIntermediateConfiguration(), finiteStrainRotation(), initialSetup(), updateQpState(), and updateQpStateSingleModel().

_perform_finite_strain_rotations

const bool ComputeMultipleInelasticStress::_perform_finite_strain_rotations

protected

after updateQpState, rotate the stress, elastic_strain, inelastic_strain and Jacobian_mult using _rotation_increment

Definition at line 129 of file ComputeMultipleInelasticStress.h.

Referenced by ComputeSmearedCrackingStress::computeQpStress(), ComputeMultipleInelasticCosseratStress::computeQpStress(), computeQpStress(), computeQpStressIntermediateConfiguration(), updateQpState(), and updateQpStateSingleModel().

_relative_tolerance

const Real ComputeMultipleInelasticStress::_relative_tolerance

protected

Definition at line 123 of file ComputeMultipleInelasticStress.h.

Referenced by updateQpState().

_rotation_increment

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainElasticStress::_rotation_increment

protectedinherited

Rotation increment material property.

Definition at line 43 of file ComputeFiniteStrainElasticStress.h.

Referenced by computeAdmissibleState(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeDamageStress::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), ComputeMultipleInelasticCosseratStress::computeQpStress(), computeQpStress(), and finiteStrainRotation().

_strain_increment

const MaterialProperty<RankTwoTensor>& ComputeMultipleInelasticStress::_strain_increment

protected

Definition at line 133 of file ComputeMultipleInelasticStress.h.

Referenced by ComputeSmearedCrackingStress::computeQpStress(), computeQpStressIntermediateConfiguration(), updateQpState(), and updateQpStateSingleModel().

_stress

MaterialProperty<RankTwoTensor>& ComputeStressBase::_stress

protectedinherited

Stress material property.

Definition at line 50 of file ComputeStressBase.h.

Referenced by ComputeMultipleInelasticCosseratStress::computeAdmissibleState(), computeAdmissibleState(), ComputeStressBase::computeQpProperties(), ComputeStrainIncrementBasedStress::computeQpStress(), ComputeLinearElasticStress::computeQpStress(), ComputeDamageStress::computeQpStress(), ComputeFiniteStrainElasticStress::computeQpStress(), ComputeCosseratLinearElasticStress::computeQpStress(), ComputeSmearedCrackingStress::computeQpStress(), ComputeLinearElasticPFFractureStress::computeQpStress(), FiniteStrainPlasticMaterial::computeQpStress(), ComputeMultiPlasticityStress::computeQpStress(), ComputeLinearViscoelasticStress::computeQpStress(), computeQpStress(), computeQpStressIntermediateConfiguration(), ComputeLinearElasticPFFractureStress::computeStrainSpectral(), ComputeLinearElasticPFFractureStress::computeStrainVolDev(), ComputeLinearElasticPFFractureStress::computeStressSpectral(), finiteStrainRotation(), ComputeStressBase::initQpStatefulProperties(), FiniteStrainCrystalPlasticity::initQpStatefulProperties(), FiniteStrainUObasedCP::initQpStatefulProperties(), FiniteStrainHyperElasticViscoPlastic::initQpStatefulProperties(), ComputeMultiPlasticityStress::postReturnMap(), FiniteStrainUObasedCP::postSolveQp(), FiniteStrainHyperElasticViscoPlastic::postSolveQp(), FiniteStrainCrystalPlasticity::postSolveQp(), ComputeSmearedCrackingStress::updateCrackingStateAndStress(), updateQpState(), and updateQpStateSingleModel().

_stress_old

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainElasticStress::_stress_old

protectedinherited

Old state of the stress tensor material property.

Definition at line 45 of file ComputeFiniteStrainElasticStress.h.

Referenced by computeAdmissibleState(), computeQpStress(), computeQpStressIntermediateConfiguration(), updateQpState(), and updateQpStateSingleModel().

_tangent_calculation_method

TangentCalculationMethod ComputeMultipleInelasticStress::_tangent_calculation_method

protected

Calculation method for the tangent modulus.

Definition at line 152 of file ComputeMultipleInelasticStress.h.

Referenced by computeAdmissibleState(), computeQpJacobianMult(), finiteStrainRotation(), initialSetup(), and updateQpStateSingleModel().

_tangent_computation_flag

std::vector<bool> ComputeMultipleInelasticStress::_tangent_computation_flag

protected

Flags to compute tangent during updateState call.

Definition at line 149 of file ComputeMultipleInelasticStress.h.

Referenced by computeAdmissibleState(), and initialSetup().

_tangent_operator_type

enum ComputeMultipleInelasticStress::TangentOperatorEnum ComputeMultipleInelasticStress::_tangent_operator_type

protected

Referenced by ComputeMultipleInelasticCosseratStress::computeQpJacobianMult(), and initialSetup().

_undamaged_stress_old

RankTwoTensor ComputeMultipleInelasticStress::_undamaged_stress_old

protected

Definition at line 188 of file ComputeMultipleInelasticStress.h.

Referenced by computeAdmissibleState(), computeQpStress(), computeQpStressIntermediateConfiguration(), updateQpState(), and updateQpStateSingleModel().

---

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

• tensor_mechanics/include/materials/ComputeMultipleInelasticStress.h
• tensor_mechanics/src/materials/ComputeMultipleInelasticStress.C