ComputeCreepPlasticityStress 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 <ComputeCreepPlasticityStress.h>

Inheritance diagram for ComputeCreepPlasticityStress:

[legend]

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ComputeCreepPlasticityStress (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 Types
enum class	TangentOperatorEnum { elastic , nonlinear }
	what sort of Tangent operator to calculate More...

Protected Member Functions
virtual std::vector< MaterialName >	getInelasticModelNames () override

virtual void	updateQpState (RankTwoTensor &elastic_strain_increment, RankTwoTensor &combined_inelastic_strain_increment) override
	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]).

void	computeStress (const RankTwoTensor &elastic_strain_increment, RankTwoTensor &stress)

void	computeTangentOperators ()

void	computeInelasticStrainIncrements (Real effective_trial_stress, const RankTwoTensor &deviatoric_trial_stress)

void	finalizeConstitutiveModels ()

virtual void	initQpStatefulProperties () override

virtual void	computeQpStress () override
	Compute the stress and store it in the _stress material property for the current quadrature point.

virtual void	computeQpStressIntermediateConfiguration ()
	Compute the stress for the current QP, but do not rotate tensors from the intermediate configuration to the new configuration.

virtual void	finiteStrainRotation (const bool force_elasticity_rotation=false)
	Rotate _elastic_strain, _stress, _inelastic_strain, and _Jacobian_mult to the new configuration.

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])

virtual void	computeQpJacobianMult ()
	Using _elasticity_tensor[_qp] and the consistent tangent operators, _consistent_tangent_operator[...] computed by the inelastic models, compute _Jacobian_mult[_qp].

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.

virtual void	computeQpProperties () override

bool	hasGuaranteedMaterialProperty (const MaterialPropertyName &prop, Guarantee guarantee)

Protected Attributes
PowerLawCreepStressUpdate *	_creep_model

IsotropicPlasticityStressUpdate *	_plasticity_model

Real	_effective_creep_strain_increment

Real	_effective_plastic_strain_increment

std::array< RankTwoTensor, 2 >	_inelastic_strain_increment

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

MaterialProperty< RankTwoTensor > &	_inelastic_strain
	The sum of the inelastic strains that come from the plastic models.

const MaterialProperty< RankTwoTensor > &	_inelastic_strain_old
	old value of inelastic strain

enum ComputeMultipleInelasticStressBase::TangentOperatorEnum	_tangent_operator_type

unsigned	_num_models
	number of plastic models

std::vector< bool >	_tangent_computation_flag
	Flags to compute tangent during updateState call.

TangentCalculationMethod	_tangent_calculation_method
	Calculation method for the tangent modulus.

std::vector< Real >	_inelastic_weights
	_inelastic_strain = sum_i (_inelastic_weights_i * inelastic_strain_from_model_i)

std::vector< RankFourTensor >	_consistent_tangent_operator
	the consistent tangent operators computed by each plastic model

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

MaterialProperty< Real > &	_material_timestep_limit

const RankFourTensor	_identity_symmetric_four
	Rank four symmetric identity tensor.

std::vector< StressUpdateBase * >	_models
	The user supplied list of inelastic models to use in the simulation.

bool	_is_elasticity_tensor_guaranteed_isotropic
	is the elasticity tensor guaranteed to be isotropic?

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

DamageBaseTempl< false > *	_damage_model
	Pointer to the damage model.

RankTwoTensor	_undamaged_stress_old

const std::string	_elasticity_tensor_name
	Name of the elasticity tensor material property.

const MaterialProperty< RankFourTensor > &	_elasticity_tensor
	Elasticity tensor material property.

MaterialProperty< RankTwoTensor > &	_rotation_total
	Rotation up to current step "n" to compute anisotropic elasticity tensor.

const MaterialProperty< RankTwoTensor > &	_rotation_total_old
	Rotation up to "n - 1" (previous) step to compute anisotropic elasticity tensor.

const MaterialProperty< RankTwoTensor > &	_rotation_increment
	Rotation increment material property.

const MaterialProperty< RankTwoTensor > &	_stress_old
	Old state of the stress tensor material property.

const std::string	_base_name
	Base name prepended to all material property names to allow for multi-material systems.

const MaterialProperty< RankTwoTensor > &	_mechanical_strain
	Mechanical strain material property.

MaterialProperty< RankTwoTensor > &	_stress
	Stress material property.

MaterialProperty< RankTwoTensor > &	_elastic_strain
	Elastic strain material property.

const MaterialProperty< RankTwoTensor > &	_extra_stress
	Extra stress tensor.

std::vector< const Function * >	_initial_stress_fcn
	initial stress components

MaterialProperty< RankFourTensor > &	_Jacobian_mult
	derivative of stress w.r.t. strain (_dstress_dstrain)


const unsigned int	_max_iterations
	Input parameters associated with the recompute iteration to return the stress state to the yield surface.

const Real	_relative_tolerance

const Real	_absolute_tolerance

const bool	_internal_solve_full_iteration_history


const MaterialProperty< RankTwoTensor > &	_elastic_strain_old
	Strain tensors.

const MaterialProperty< RankTwoTensor > &	_strain_increment

Private Member Functions
bool	haveMaterialProperty (const std::string &prop_name)

std::vector< VariableName >	buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)

void	validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)

bool	isNotObjectVariable (const VariableName &name)

Private Attributes
FEProblemBase &	_dmi_fe_problem

const InputParameters &	_gc_params
	Parameters of the object with this interface.

FEProblemBase *const	_gc_feproblem
	Reference to the FEProblemBase class.

BlockRestrictable *const	_gc_block_restrict
	Access block restrictions of the object with this interface.

Detailed Description

ComputeCreepPlasticityStress 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 specifically to combine creep and plasticity. The solve happens for both materials simultaneously rather than with a staggered approach.

Definition at line 29 of file ComputeCreepPlasticityStress.h.

Member Enumeration Documentation

◆ TangentOperatorEnum

enum class ComputeMultipleInelasticStressBase::TangentOperatorEnum

strongprotectedinherited

what sort of Tangent operator to calculate

Enumerator
elastic
nonlinear

Definition at line 139 of file ComputeMultipleInelasticStressBase.h.

139{ elastic, nonlinear } _tangent_operator_type;

ComputeMultipleInelasticStressBase::TangentOperatorEnum::nonlinear

@ nonlinear

ComputeMultipleInelasticStressBase::TangentOperatorEnum::elastic

@ elastic

ComputeMultipleInelasticStressBase::_tangent_operator_type

enum ComputeMultipleInelasticStressBase::TangentOperatorEnum _tangent_operator_type

Constructor & Destructor Documentation

◆ ComputeCreepPlasticityStress()

ComputeCreepPlasticityStress::ComputeCreepPlasticityStress ( const InputParameters & parameters )

Definition at line 34 of file ComputeCreepPlasticityStress.C.

  : ComputeMultipleInelasticStressBase(parameters)
{
}

Member Function Documentation

◆ computeAdmissibleState()

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

protectedvirtualinherited

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 356 of file ComputeMultipleInelasticStressBase.C.

{
  // Reset properties to the beginning of the time step (necessary if substepping is employed).
  _models[model_number]->resetIncrementalMaterialProperties();
 
  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 if (_models[model_number]->substeppingCapabilityEnabled() &&
           (_is_elasticity_tensor_guaranteed_isotropic || !_perform_finite_strain_rotations))
    _models[model_number]->updateStateSubstep(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);
  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(), ComputeMultipleInelasticStress::updateQpState(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

◆ computeInelasticStrainIncrements()

void ComputeCreepPlasticityStress::computeInelasticStrainIncrements	(	Real	effective_trial_stress,
		const RankTwoTensor &	deviatoric_trial_stress
	)

protected

Definition at line 361 of file ComputeCreepPlasticityStress.C.

{
  if (!MooseUtils::absoluteFuzzyEqual(effective_trial_stress, 0.0))
  {
    if (_effective_creep_strain_increment != 0.0)
      _inelastic_strain_increment[0] =
          deviatoric_trial_stress *
          (1.5 * _effective_creep_strain_increment / effective_trial_stress);
    else
      _inelastic_strain_increment[0].zero();
 
    if (_effective_plastic_strain_increment != 0.0)
      _inelastic_strain_increment[1] =
          deviatoric_trial_stress *
          (1.5 * _effective_plastic_strain_increment / effective_trial_stress);
    else
      _inelastic_strain_increment[1].zero();
  }
  else
  {
    _inelastic_strain_increment[0].zero();
    _inelastic_strain_increment[1].zero();
  }
}

Referenced by updateQpState().

◆ computeQpJacobianMult()

void ComputeMultipleInelasticStressBase::computeQpJacobianMult ( )

protectedvirtualinherited

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 283 of file ComputeMultipleInelasticStressBase.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 (const auto i_rmm : make_range(_num_models))
      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 (const auto i_rmm : make_range(1u, _num_models))
      _Jacobian_mult[_qp] = _consistent_tangent_operator[i_rmm] * E_inv * _Jacobian_mult[_qp];
  }
}

Referenced by updateQpState(), and ComputeMultipleInelasticStress::updateQpState().

◆ computeQpProperties()

void ComputeGeneralStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

Definition at line 44 of file ComputeGeneralStressBase.C.

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

◆ computeQpStress()

void ComputeMultipleInelasticStressBase::computeQpStress ( )

overrideprotectedvirtualinherited

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 199 of file ComputeMultipleInelasticStressBase.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 (_material_timestep_limit[_qp] > damage_timestep_limit)
      _material_timestep_limit[_qp] = damage_timestep_limit;
  }
 
  if (_perform_finite_strain_rotations)
    finiteStrainRotation();
}

◆ computeQpStressIntermediateConfiguration()

void ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration ( )

protectedvirtualinherited

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

Definition at line 227 of file ComputeMultipleInelasticStressBase.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];
 
    _material_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 ComputeMultipleInelasticStressBase::computeQpStress(), and ComputeSmearedCrackingStress::computeQpStress().

◆ computeStress()

void ComputeCreepPlasticityStress::computeStress	(	const RankTwoTensor &	elastic_strain_increment,
		RankTwoTensor &	stress
	)

protected

Definition at line 322 of file ComputeCreepPlasticityStress.C.

{
  // form the stress, with the check for changed elasticity constants
  if (_is_elasticity_tensor_guaranteed_isotropic || !_perform_finite_strain_rotations)
    stress = _elasticity_tensor[_qp] * (_elastic_strain_old[_qp] + elastic_strain_increment);
  else
  {
    if (_damage_model)
      stress = _undamaged_stress_old + _elasticity_tensor[_qp] * elastic_strain_increment;
    else
      stress = _stress_old[_qp] + _elasticity_tensor[_qp] * elastic_strain_increment;
  }
}

Referenced by computeTangentOperators(), and updateQpState().

◆ computeTangentOperators()

void ComputeCreepPlasticityStress::computeTangentOperators ( )

protected

Definition at line 338 of file ComputeCreepPlasticityStress.C.

{
  auto elastic_strain_inc = _strain_increment[_qp] - _inelastic_strain_increment[1];
 
  RankTwoTensor stress;
  computeStress(elastic_strain_inc, stress);
  RankTwoTensor deviatoric_stress = stress.deviatoric();
  Real effective_stress = std::sqrt(1.5 * deviatoric_stress.doubleContraction(deviatoric_stress));
 
  _creep_model->computeTangentOperator(
      effective_stress, _stress[_qp], _consistent_tangent_operator[0]);
 
  elastic_strain_inc = _strain_increment[_qp] - _inelastic_strain_increment[0];
 
  computeStress(elastic_strain_inc, stress);
  deviatoric_stress = stress.deviatoric();
  effective_stress = std::sqrt(1.5 * deviatoric_stress.doubleContraction(deviatoric_stress));
 
  _plasticity_model->computeTangentOperator(
      effective_stress, _stress[_qp], _consistent_tangent_operator[1]);
}

Referenced by updateQpState().

◆ finalizeConstitutiveModels()

void ComputeCreepPlasticityStress::finalizeConstitutiveModels ( )

protected

Definition at line 388 of file ComputeCreepPlasticityStress.C.

{
  _creep_model->updateEffectiveInelasticStrainIncrement(_effective_creep_strain_increment);
  _plasticity_model->updateEffectiveInelasticStrainIncrement(_effective_plastic_strain_increment);
  _creep_model->updateEffectiveInelasticStrain(_effective_creep_strain_increment);
  _plasticity_model->updateEffectiveInelasticStrain(_effective_plastic_strain_increment);
  _creep_model->resetIncrementalMaterialProperties();
  _creep_model->computeStressFinalize(_inelastic_strain_increment[0]);
  _plasticity_model->computeStressFinalize(_inelastic_strain_increment[1]);
}

Referenced by updateQpState().

◆ finiteStrainRotation()

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

protectedvirtualinherited

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 267 of file ComputeMultipleInelasticStressBase.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 ComputeMultipleInelasticStressBase::computeQpStress(), and ComputeSmearedCrackingStress::computeQpStress().

◆ getInelasticModelNames()

std::vector< MaterialName > ComputeCreepPlasticityStress::getInelasticModelNames ( )

overrideprotectedvirtual

Implements ComputeMultipleInelasticStressBase.

Definition at line 40 of file ComputeCreepPlasticityStress.C.

{
  std::vector<MaterialName> names = {getParam<MaterialName>("creep_model"),
                                     getParam<MaterialName>("plasticity_model")};
  return names;
}

◆ 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::computeQpStress(), LagrangianStressDivergenceBase::initialSetup(), ADComputeMultipleInelasticStress::initialSetup(), ADComputeSmearedCrackingStress::initialSetup(), ComputeLinearElasticPFFractureStress::initialSetup(), ComputeMultipleInelasticStressBase::initialSetup(), ComputeSmearedCrackingStress::initialSetup(), ComputeLagrangianStressBase::initialSetup(), ComputeStVenantKirchhoffStress::initialSetup(), and CriticalTimeStep::initialSetup().

◆ initialSetup()

void ComputeCreepPlasticityStress::initialSetup ( )

overridevirtual

Reimplemented from ComputeMultipleInelasticStressBase.

Definition at line 48 of file ComputeCreepPlasticityStress.C.

{
  ComputeMultipleInelasticStressBase::initialSetup();
 
  if (_models.size() != 2)
    mooseError("Error in ComputeCreepPlasticityStress: two models are required.");
 
  _creep_model = dynamic_cast<PowerLawCreepStressUpdate *>(_models[0]);
  if (!_creep_model)
    mooseError("Model " + getParam<MaterialName>("creep_model") +
               " is not a compatible creep model in ComputeCreepPlasticityStress.");
 
  _plasticity_model = dynamic_cast<IsotropicPlasticityStressUpdate *>(_models[1]);
  if (!_plasticity_model)
    mooseError("Model " + getParam<MaterialName>("plasticity_model") +
               " is not a compatible plasticity model in ComputeCreepPlasticityStress.");
}

◆ initQpStatefulProperties()

void ComputeMultipleInelasticStressBase::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Reimplemented from ComputeFiniteStrainElasticStress.

Reimplemented in ComputeMultipleInelasticCosseratStress, and ComputeSmearedCrackingStress.

Definition at line 93 of file ComputeMultipleInelasticStressBase.C.

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

◆ updateQpState()

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

overrideprotectedvirtual

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.

Implements ComputeMultipleInelasticStressBase.

Definition at line 67 of file ComputeCreepPlasticityStress.C.

{
  if (_internal_solve_full_iteration_history == true)
  {
    _console << std::endl
             << "iteration output for ComputeCreepPlasticityStress solve:"
             << " time=" << _t << " int_pt=" << _qp << std::endl;
  }
 
  for (const auto i_rmm : index_range(_models))
    _inelastic_strain_increment[i_rmm].zero();
 
  elastic_strain_increment = _strain_increment[_qp];
  computeStress(elastic_strain_increment, _stress[_qp]);
 
  const RankTwoTensor trial_stress = _stress[_qp];
  const RankTwoTensor deviatoric_trial_stress = trial_stress.deviatoric();
  const Real effective_trial_stress =
      std::sqrt(1.5 * deviatoric_trial_stress.doubleContraction(deviatoric_trial_stress));
 
  //
  // First, compute creep response assuming no plasticity
  //
  _creep_model->setQp(_qp);
  computeAdmissibleState(
      0, elastic_strain_increment, _inelastic_strain_increment[0], _consistent_tangent_operator[0]);
  _effective_creep_strain_increment = _creep_model->effectiveInelasticStrainIncrement();
 
  //
  // Now check the plasticity model.
  // If no yielding, we are done.
  //
  _plasticity_model->setQp(_qp);
  computeAdmissibleState(
      1, elastic_strain_increment, _inelastic_strain_increment[1], _consistent_tangent_operator[1]);
  _effective_plastic_strain_increment = _plasticity_model->effectiveInelasticStrainIncrement();
  const Real yield_condition = _plasticity_model->yieldCondition();
 
  if (yield_condition > 0.0) // yielding even with maximum creep strain
  {
    const Real max_effective_creep_strain_increment = _effective_creep_strain_increment;
 
    //
    // Compute plastic response given no creep
    //
    elastic_strain_increment = _strain_increment[_qp];
    _stress[_qp] = trial_stress;
    computeAdmissibleState(1,
                           elastic_strain_increment,
                           _inelastic_strain_increment[1],
                           _consistent_tangent_operator[1]);
    const Real max_effective_plastic_strain_increment =
        _plasticity_model->effectiveInelasticStrainIncrement();
 
    //
    // Reset plasticity model
    //
    _plasticity_model->computeStressInitialize(effective_trial_stress, _elasticity_tensor[_qp]);
 
    //
    // The residual for the creep law is (creep rate)*dt - (creep strain increment)
    // We want the creep rate calculation to depend on both creep and plasticity.
    // Since we send in the total scalar inelastic strain (creep and plasticity), we need to correct
    // the second term by adding the plastic strain.
    //
    Real creep_residual = _creep_model->computeResidual(effective_trial_stress,
                                                        _effective_creep_strain_increment +
                                                            _effective_plastic_strain_increment);
    creep_residual += _effective_plastic_strain_increment;
 
    //
    // We want the plasticity resdiual to be (effective_trial_stress - r - yield_stress)/3G - (total
    // inelastic increment)
    // The standard residual for plasticity is (effective_trial_stress - r - yield_stress)/3G -
    // (plasticity strain increment)
    // Since we send in the plastic inelastic strain (for calculation of r), we must subtract the
    // creep strain to correct the final term in our desired residual
    //
    Real plasticity_residual = _plasticity_model->computeResidual(
        effective_trial_stress, _effective_plastic_strain_increment);
    plasticity_residual -= _effective_creep_strain_increment;
 
    unsigned int counter = 0;
    Real residual = 0;
    const Real threeG =
        3.0 * ElasticityTensorTools::getIsotropicShearModulus(_elasticity_tensor[_qp]);
    const Real reference_residual = effective_trial_stress / threeG;
 
    do
    {
      //
      // Solve Ax=b where x is the vector of creep and plastic inelastic strains
      //
      // x1 => creep
      // x2 => plasticity
      //
      //
      // A11 (creep,creep)
      //
      Real A11 = _creep_model->computeDerivative(effective_trial_stress,
                                                 _effective_creep_strain_increment +
                                                     _effective_plastic_strain_increment);
      //
      // A12 (creep,plasticity)
      //
      Real A12 = A11 + 1.0;
      //
      // A21 (plasticity,creep)
      //
      Real A21 = -1.0;
      //
      // A22 (plasticity,plasticity)
      //
      Real A22 = _plasticity_model->computeDerivative(effective_trial_stress,
                                                      _effective_plastic_strain_increment);
 
      Real rhs1 = creep_residual;
      Real rhs2 = plasticity_residual;
 
      //
      // Solve
      //
      // A11*a + A21 = 0
      // A11*a = -A21
      // a = -A21/A11
      const Real a = -A21 / A11;
      A21 = 0;
      A22 += a * A12;
      rhs2 += a * rhs1;
 
      Real x2 = rhs2 / A22;
      Real x1 = (rhs1 - A12 * x2) / A11;
 
      while (_effective_creep_strain_increment - x1 < 0)
        x1 *= 0.5;
      while (_effective_plastic_strain_increment - x2 < 0)
        x2 *= 0.5;
 
      while (_effective_creep_strain_increment - x1 > max_effective_creep_strain_increment)
        x1 *= 0.5;
      while (_effective_plastic_strain_increment - x2 > max_effective_plastic_strain_increment)
        x2 *= 0.5;
 
      // This check is to avoid a fpe in the creep law.
      // Maybe it is better to check for this condition in the creep law
      if (effective_trial_stress < threeG * (_effective_creep_strain_increment - x1 +
                                             _effective_plastic_strain_increment - x2))
      {
        const Real sum =
            _effective_creep_strain_increment - x1 + _effective_plastic_strain_increment - x2;
        const Real factor = 0.9 * effective_trial_stress / (threeG * sum);
        const Real cNew = (_effective_creep_strain_increment - x1) * factor;
        const Real pNew = (_effective_plastic_strain_increment - x2) * factor;
        x1 = _effective_creep_strain_increment - cNew;
        x2 = _effective_plastic_strain_increment - pNew;
      }
 
      _effective_creep_strain_increment -= x1;
      _effective_plastic_strain_increment -= x2;
 
      //
      // The residual for the creep law is (creep rate)*dt - (creep strain increment)
      //
      creep_residual = _creep_model->computeResidual(effective_trial_stress,
                                                     _effective_creep_strain_increment +
                                                         _effective_plastic_strain_increment);
      creep_residual += _effective_plastic_strain_increment;
      //
      // The residual for plasticity is (effective_trial_stress - r - yield_stress)/3G - scalar
      //
      plasticity_residual = _plasticity_model->computeResidual(effective_trial_stress,
                                                               _effective_plastic_strain_increment);
      plasticity_residual -= _effective_creep_strain_increment;
 
      residual =
          std::sqrt(creep_residual * creep_residual + plasticity_residual * plasticity_residual);
 
      if (_internal_solve_full_iteration_history == true)
      {
        _console << "stress iteration number = " << counter << "\n relative residual = "
                 << (0 == reference_residual ? 0 : residual / reference_residual)
                 << "\n stress convergence relative tolerance = " << _relative_tolerance
                 << "\n absolute residual = " << residual
                 << "\n creep residual = " << creep_residual
                 << "\n plasticity residual = " << plasticity_residual
                 << "\n creep iteration increment = " << x1
                 << "\n plasticity iteration increment = " << x2
                 << "\n creep scalar strain = " << _effective_creep_strain_increment
                 << "\n plasticity scalar strain = " << _effective_plastic_strain_increment
                 << "\n max creep scalar strain = " << max_effective_creep_strain_increment
                 << "\n max plasticity scalar strain = " << max_effective_plastic_strain_increment
                 << "\n stress convergence absolute tolerance = " << _absolute_tolerance
                 << std::endl;
      }
      ++counter;
    } while (counter < _max_iterations && residual > _absolute_tolerance &&
             (residual / reference_residual) > _relative_tolerance);
 
    if (counter == _max_iterations && residual > _absolute_tolerance &&
        (residual / reference_residual) > _relative_tolerance)
    {
      std::stringstream msg;
      msg << "\n relative residual = "
          << (0 == reference_residual ? 0 : residual / reference_residual)
          << "\n stress convergence relative tolerance = " << _relative_tolerance
          << "\n absolute residual = " << residual << "\n creep residual = " << creep_residual
          << "\n plasticity residual = " << plasticity_residual
          << "\n creep scalar strain = " << _effective_creep_strain_increment
          << "\n plasticity scalar strain = " << _effective_plastic_strain_increment
          << "\n max creep scalar strain = " << max_effective_creep_strain_increment
          << "\n max plasticity scalar strain = " << max_effective_plastic_strain_increment
          << "\n stress convergence absolute tolerance = " << _absolute_tolerance << std::endl;
 
      throw MooseException("Max stress iteration hit during ComputeCreepPlasticityStress solve! " +
                           _name + msg.str());
    }
  }
  computeInelasticStrainIncrements(effective_trial_stress, deviatoric_trial_stress);
 
  finalizeConstitutiveModels();
 
  combined_inelastic_strain_increment.zero();
  for (const auto i_rmm : make_range(_num_models))
    combined_inelastic_strain_increment += _inelastic_strain_increment[i_rmm];
 
  if (yield_condition > 0.0)
  {
    elastic_strain_increment = _strain_increment[_qp] - combined_inelastic_strain_increment;
    _stress[_qp] = _elasticity_tensor[_qp] * (elastic_strain_increment + _elastic_strain_old[_qp]);
 
    computeTangentOperators();
  }
  else
  {
    // The tangent for creep was called in the initial creep solve.
 
    // Set the tangent for plasticity to zero
    _consistent_tangent_operator[1].zero();
  }
 
  if (_fe_problem.currentlyComputingJacobian())
    computeQpJacobianMult();
 
  _material_timestep_limit[_qp] = 0.0;
  for (const auto i_rmm : make_range(_num_models))
    _material_timestep_limit[_qp] += 1.0 / _models[i_rmm]->computeTimeStepLimit();
 
  if (MooseUtils::absoluteFuzzyEqual(_material_timestep_limit[_qp], 0.0))
    _material_timestep_limit[_qp] = std::numeric_limits<Real>::max();
  else
    _material_timestep_limit[_qp] = 1.0 / _material_timestep_limit[_qp];
}

◆ updateQpStateSingleModel()

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

protectedvirtualinherited

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 305 of file ComputeMultipleInelasticStressBase.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];
  }
 
  _material_timestep_limit[_qp] = _models[0]->computeTimeStepLimit();
 
  /* propagate internal variables, etc, to this timestep for those inelastic models where
   * "updateState" is not called */
  for (const auto i_rmm : make_range(_num_models))
    if (i_rmm != model_number)
      _models[i_rmm]->propagateQpStatefulProperties();
}

Referenced by ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration().

◆ validParams()

InputParameters ComputeCreepPlasticityStress::validParams ( )

static

Definition at line 20 of file ComputeCreepPlasticityStress.C.

{
  InputParameters params = ComputeMultipleInelasticStressBase::validParams();
  params.addClassDescription("Compute state (stress and internal parameters such as inelastic "
                             "strains and internal parameters) using an Newton process for one "
                             "creep and one plasticity model");
  params.addRequiredParam<MaterialName>("creep_model",
                                        "Creep model that derives from PowerLawCreepStressUpdate.");
  params.addRequiredParam<MaterialName>(
      "plasticity_model", "Plasticity model that derives from IsotropicPlasticityStressUpdate.");
 
  return params;
}

Member Data Documentation

◆ _absolute_tolerance

const Real ComputeMultipleInelasticStressBase::_absolute_tolerance

protectedinherited

Definition at line 120 of file ComputeMultipleInelasticStressBase.h.

Referenced by updateQpState(), and ComputeMultipleInelasticStress::updateQpState().

◆ _all_models_isotropic

bool ComputeMultipleInelasticStressBase::_all_models_isotropic

protectedinherited

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

Definition at line 179 of file ComputeMultipleInelasticStressBase.h.

Referenced by ComputeMultipleInelasticStressBase::finiteStrainRotation(), and ComputeMultipleInelasticStressBase::initialSetup().

◆ _base_name

const std::string ComputeGeneralStressBase::_base_name

protectedinherited

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

Definition at line 43 of file ComputeGeneralStressBase.h.

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

◆ _consistent_tangent_operator

std::vector<RankFourTensor> ComputeMultipleInelasticStressBase::_consistent_tangent_operator

protectedinherited

the consistent tangent operators computed by each plastic model

Definition at line 154 of file ComputeMultipleInelasticStressBase.h.

Referenced by ComputeMultipleInelasticStressBase::computeQpJacobianMult(), ComputeMultipleInelasticCosseratStress::computeQpJacobianMult(), computeTangentOperators(), ComputeMultipleInelasticStressBase::initialSetup(), updateQpState(), ComputeMultipleInelasticStress::updateQpState(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

◆ _creep_model

PowerLawCreepStressUpdate* ComputeCreepPlasticityStress::_creep_model

protected

Definition at line 60 of file ComputeCreepPlasticityStress.h.

Referenced by computeTangentOperators(), finalizeConstitutiveModels(), initialSetup(), and updateQpState().

◆ _cycle_models

const bool ComputeMultipleInelasticStressBase::_cycle_models

protectedinherited

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

Definition at line 157 of file ComputeMultipleInelasticStressBase.h.

Referenced by ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration().

◆ _damage_model

DamageBaseTempl<false>* ComputeMultipleInelasticStressBase::_damage_model

protectedinherited

Pointer to the damage model.

Definition at line 182 of file ComputeMultipleInelasticStressBase.h.

Referenced by ComputeMultipleInelasticStressBase::computeAdmissibleState(), ComputeMultipleInelasticStressBase::computeQpStress(), ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration(), computeStress(), ComputeMultipleInelasticStressBase::initialSetup(), ComputeMultipleInelasticStress::updateQpState(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

◆ _effective_creep_strain_increment

Real ComputeCreepPlasticityStress::_effective_creep_strain_increment

protected

Definition at line 63 of file ComputeCreepPlasticityStress.h.

Referenced by computeInelasticStrainIncrements(), finalizeConstitutiveModels(), and updateQpState().

◆ _effective_plastic_strain_increment

Real ComputeCreepPlasticityStress::_effective_plastic_strain_increment

protected

Definition at line 64 of file ComputeCreepPlasticityStress.h.

Referenced by computeInelasticStrainIncrements(), finalizeConstitutiveModels(), and updateQpState().

◆ _elastic_strain

MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_elastic_strain

protectedinherited

Elastic strain material property.

Definition at line 50 of file ComputeGeneralStressBase.h.

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

◆ _elastic_strain_old

const MaterialProperty<RankTwoTensor>& ComputeMultipleInelasticStressBase::_elastic_strain_old

protectedinherited

Strain tensors.

Definition at line 128 of file ComputeMultipleInelasticStressBase.h.

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

◆ _elasticity_tensor

const MaterialProperty<RankFourTensor>& ComputeFiniteStrainElasticStress::_elasticity_tensor

protectedinherited

Elasticity tensor material property.

Definition at line 35 of file ComputeFiniteStrainElasticStress.h.

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

◆ _elasticity_tensor_name

const std::string ComputeFiniteStrainElasticStress::_elasticity_tensor_name

protectedinherited

Name of the elasticity tensor material property.

Definition at line 33 of file ComputeFiniteStrainElasticStress.h.

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

◆ _extra_stress

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_extra_stress

protectedinherited

Extra stress tensor.

Definition at line 53 of file ComputeGeneralStressBase.h.

Referenced by ComputeGeneralStressBase::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 ComputeMultipleInelasticStressBase::_identity_symmetric_four

protectedinherited

Rank four symmetric identity tensor.

Definition at line 164 of file ComputeMultipleInelasticStressBase.h.

Referenced by ComputeMultipleInelasticStressBase::computeQpJacobianMult(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

◆ _inelastic_strain

MaterialProperty<RankTwoTensor>& ComputeMultipleInelasticStressBase::_inelastic_strain

protectedinherited

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

Definition at line 133 of file ComputeMultipleInelasticStressBase.h.

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

◆ _inelastic_strain_increment

std::array<RankTwoTensor, 2> ComputeCreepPlasticityStress::_inelastic_strain_increment

protected

Definition at line 66 of file ComputeCreepPlasticityStress.h.

Referenced by computeInelasticStrainIncrements(), computeTangentOperators(), finalizeConstitutiveModels(), and updateQpState().

◆ _inelastic_strain_old

const MaterialProperty<RankTwoTensor>& ComputeMultipleInelasticStressBase::_inelastic_strain_old

protectedinherited

old value of inelastic strain

Definition at line 136 of file ComputeMultipleInelasticStressBase.h.

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

◆ _inelastic_weights

std::vector<Real> ComputeMultipleInelasticStressBase::_inelastic_weights

protectedinherited

_inelastic_strain = sum_i (_inelastic_weights_i * inelastic_strain_from_model_i)

Definition at line 151 of file ComputeMultipleInelasticStressBase.h.

Referenced by ComputeMultipleInelasticStressBase::initialSetup(), and ComputeMultipleInelasticStress::updateQpState().

◆ _initial_stress_fcn

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

protectedinherited

initial stress components

Definition at line 56 of file ComputeGeneralStressBase.h.

◆ _internal_solve_full_iteration_history

const bool ComputeMultipleInelasticStressBase::_internal_solve_full_iteration_history

protectedinherited

Definition at line 121 of file ComputeMultipleInelasticStressBase.h.

Referenced by updateQpState(), and ComputeMultipleInelasticStress::updateQpState().

◆ _is_elasticity_tensor_guaranteed_isotropic

bool ComputeMultipleInelasticStressBase::_is_elasticity_tensor_guaranteed_isotropic

protectedinherited

is the elasticity tensor guaranteed to be isotropic?

Definition at line 176 of file ComputeMultipleInelasticStressBase.h.

Referenced by ComputeMultipleInelasticStressBase::computeAdmissibleState(), ComputeMultipleInelasticStressBase::computeQpStressIntermediateConfiguration(), computeStress(), ComputeMultipleInelasticStressBase::finiteStrainRotation(), ComputeMultipleInelasticStressBase::initialSetup(), ComputeMultipleInelasticStress::updateQpState(), and ComputeMultipleInelasticStressBase::updateQpStateSingleModel().

◆ _Jacobian_mult

MaterialProperty<RankFourTensor>& ComputeGeneralStressBase::_Jacobian_mult

protectedinherited

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

Definition at line 59 of file ComputeGeneralStressBase.h.

◆ _material_timestep_limit

MaterialProperty<Real>& ComputeMultipleInelasticStressBase::_material_timestep_limit

protectedinherited

Definition at line 159 of file ComputeMultipleInelasticStressBase.h.

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

◆ _max_iterations

const unsigned int ComputeMultipleInelasticStressBase::_max_iterations

protectedinherited

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

Definition at line 118 of file ComputeMultipleInelasticStressBase.h.

Referenced by updateQpState(), and ComputeMultipleInelasticStress::updateQpState().

◆ _mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_mechanical_strain

protectedinherited

Mechanical strain material property.

Definition at line 46 of file ComputeGeneralStressBase.h.

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

◆ _models

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

protectedinherited

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 173 of file ComputeMultipleInelasticStressBase.h.

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