ADMultiplePowerLawCreepStressUpdate Class Reference

This class uses the stress update material in a radial return isotropic creep model. More...

#include <ADMultiplePowerLawCreepStressUpdate.h>

Inheritance diagram for ADMultiplePowerLawCreepStressUpdate:

Public Types
enum	SubsteppingType { SubsteppingType::NONE, SubsteppingType::ERROR_BASED, SubsteppingType::INCREMENT_BASED }
using	GR2 = GenericRankTwoTensor< is_ad >
using	GSR2 = Moose::GenericType< RankTwoTensor, is_ad >
using	GR4 = GenericRankFourTensor< is_ad >
using	GSR4 = Moose::GenericType< RankFourTensor, is_ad >
enum	ConstantTypeEnum { ConstantTypeEnum::NONE, ConstantTypeEnum::ELEMENT, ConstantTypeEnum::SUBDOMAIN }
enum	TEST_TYPE
typedef DataFileName	DataFileParameterType

Public Member Functions
	ADMultiplePowerLawCreepStressUpdate (const InputParameters &parameters)
virtual Real	computeStrainEnergyRateDensity (const ADMaterialProperty< RankTwoTensor > &stress, const ADMaterialProperty< RankTwoTensor > &strain_rate) override
virtual bool	substeppingCapabilityEnabled () override
	Does the model include the infrastructure for substep decomposition of the elastic strain initially used to calculate the trial stress guess Inheriting classes which wish to use the substepping capability should overwrite this method and set it to return true. More...
virtual void	updateState (GenericRankTwoTensor< is_ad > &strain_increment, GenericRankTwoTensor< is_ad > &inelastic_strain_increment, const GenericRankTwoTensor< is_ad > &rotation_increment, GenericRankTwoTensor< is_ad > &stress_new, const RankTwoTensor &stress_old, const GenericRankFourTensor< is_ad > &elasticity_tensor, const RankTwoTensor &elastic_strain_old, bool compute_full_tangent_operator=false, RankFourTensor &tangent_operator=StressUpdateBaseTempl< is_ad >::_identityTensor) override
	A radial return (J2) mapping method is performed with return mapping iterations. More...
virtual void	updateState (GR2 &strain_increment, GR2 &inelastic_strain_increment, const GR2 &rotation_increment, GR2 &stress_new, const RankTwoTensor &stress_old, const GR4 &elasticity_tensor, const RankTwoTensor &elastic_strain_old, bool compute_full_tangent_operator=false, RankFourTensor &tangent_operator=StressUpdateBaseTempl< is_ad >::_identityTensor)
	Given a strain increment that results in a trial stress, perform some procedure (such as an iterative return-mapping process) to produce an admissible stress, an elastic strain increment and an inelastic strain increment. More...
virtual void	updateStateSubstepInternal (GenericRankTwoTensor< is_ad > &, GenericRankTwoTensor< is_ad > &, const GenericRankTwoTensor< is_ad > &, GenericRankTwoTensor< is_ad > &, const RankTwoTensor &, const GenericRankFourTensor< is_ad > &, const RankTwoTensor &, unsigned int total_number_substeps, bool compute_full_tangent_operator=false, RankFourTensor &tangent_operator=StressUpdateBaseTempl< is_ad >::_identityTensor)
virtual void	updateStateSubstep (GenericRankTwoTensor< is_ad > &, GenericRankTwoTensor< is_ad > &, const GenericRankTwoTensor< is_ad > &, GenericRankTwoTensor< is_ad > &, const RankTwoTensor &, const GenericRankFourTensor< is_ad > &, const RankTwoTensor &, bool compute_full_tangent_operator=false, RankFourTensor &tangent_operator=StressUpdateBaseTempl< is_ad >::_identityTensor) override
	Similar to the updateState function, this method updates the strain and stress for one substep. More...
virtual void	updateStateSubstep (GR2 &, GR2 &, const GR2 &, GR2 &, const RankTwoTensor &, const GR4 &, const RankTwoTensor &, bool compute_full_tangent_operator=false, RankFourTensor &tangent_operator=StressUpdateBaseTempl< is_ad >::_identityTensor)
	Similar to the updateState function, this method updates the strain and stress for one substep. More...
virtual Real	computeReferenceResidual (const GenericReal< is_ad > &effective_trial_stress, const GenericReal< is_ad > &scalar_effective_inelastic_strain) override
	Compute a reference quantity to be used for checking relative convergence. More...
virtual GenericReal< is_ad >	minimumPermissibleValue (const GenericReal< is_ad > &) const override
	Compute the minimum permissible value of the scalar. More...
virtual GenericReal< is_ad >	maximumPermissibleValue (const GenericReal< is_ad > &effective_trial_stress) const override
	Compute the maximum permissible value of the scalar. More...
virtual Real	computeTimeStepLimit () override
	Compute the limiting value of the time step for this material. More...
bool	requiresIsotropicTensor () override
	Does the model require the elasticity tensor to be isotropic? More...
bool	isIsotropic () override
	Radial return mapped models should be isotropic by default! More...
virtual int	calculateNumberSubsteps (const GenericRankTwoTensor< is_ad > &strain_increment) override
	If substepping is enabled, calculate the number of substeps as a function of the elastic strain increment guess and the maximum inelastic strain increment ratio based on a user-specified tolerance. More...
virtual int	calculateNumberSubsteps (const GR2 &)
	Given the elastic strain increment compute the number of substeps required to bring a substepped trial stress guess distance from the yield surface into the tolerance specified in the individual child class. More...
virtual bool	substeppingCapabilityRequested () override
	Has the user requested usage of (possibly) implemented substepping capability for inelastic models. More...
const GenericReal< is_ad > &	effectiveInelasticStrainIncrement () const
	Current value of scalar inelastic strain. More...
void	updateEffectiveInelasticStrainIncrement (const GenericReal< is_ad > &eisi)
void	updateEffectiveInelasticStrain (const GenericReal< is_ad > &increment)
void	computeTangentOperator (Real effective_trial_stress, const RankTwoTensor &stress_new, RankFourTensor &tangent_operator)
	Calculate the tangent_operator. More...
template<>
void	computeTangentOperator (Real effective_trial_stress, const RankTwoTensor &stress_new, RankFourTensor &tangent_operator)
template<>
void	computeTangentOperator (Real effective_trial_stress, const RankTwoTensor &stress_new, RankFourTensor &tangent_operator)
void	setQp (unsigned int qp)
	Sets the value of the global variable _qp for inheriting classes. More...
virtual void	storeIncrementalMaterialProperties (const unsigned int)
	Properly set up the incremental calculation storage of the stateful material properties in the inheriting classes. More...
virtual void	resetIncrementalMaterialProperties ()
	Reset material properties. More...
virtual Real	computeStrainEnergyRateDensity (const GenericMaterialProperty< RankTwoTensor, is_ad > &, const GenericMaterialProperty< RankTwoTensor, is_ad > &)
	Compute the strain energy rate density for this inelastic model for the current step. More...
virtual const dof_id_type &	getElementID (const std::string &id_parameter_name, unsigned int comp=0) const override
dof_id_type	getElementID (const Elem *elem, unsigned int elem_id_index) const
virtual const dof_id_type &	getElementIDNeighbor (const std::string &id_parameter_name, unsigned int comp=0) const override
virtual const dof_id_type &	getElementIDByName (const std::string &id_parameter_name) const override
virtual const dof_id_type &	getElementIDNeighborByName (const std::string &id_parameter_name) const override
virtual void	computeProperties () override
MaterialBase &	getMaterial (const std::string &name)
MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false, bool no_dep=false)
MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)
MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)
virtual bool	isBoundaryMaterial () const override
virtual const std::unordered_set< unsigned int > &	getMatPropDependencies () const override
virtual void	subdomainSetup () override
bool	ghostable () const override final
virtual void	resolveOptionalProperties () override
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)
const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)
const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)
const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)
virtual void	initStatefulProperties (unsigned int n_points)
virtual bool	isInterfaceMaterial ()
virtual void	computePropertiesAtQp (unsigned int qp)
const MaterialProperty< T > &	getZeroMaterialPropertyByName (Ts... args)
virtual const std::set< std::string > &	getRequestedItems () override
virtual const std::set< std::string > &	getSuppliedItems () override
const std::set< unsigned int > &	getSuppliedPropIDs ()
void	checkStatefulSanity () const
std::set< OutputName >	getOutputs ()
bool	hasStatefulProperties () const
void	setFaceInfo (const FaceInfo &fi)
void	setActiveProperties (const std::unordered_set< unsigned int > &needed_props)
bool	forceStatefulInit () const
virtual bool	enabled () const
std::shared_ptr< MooseObject >	getSharedPtr ()
std::shared_ptr< const MooseObject >	getSharedPtr () const
MooseApp &	getMooseApp () const
const std::string &	type () const
virtual const std::string &	name () const
std::string	typeAndName () const
std::string	errorPrefix (const std::string &error_type) const
void	callMooseError (std::string msg, const bool with_prefix) const
MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const
const InputParameters &	parameters () const
MooseObjectName	uniqueName () const
const T &	getParam (const std::string &name) const
std::vector< std::pair< T1, T2 > >	getParam (const std::string &param1, const std::string &param2) const
const T &	getRenamedParam (const std::string &old_name, const std::string &new_name) const
T	getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
bool	isParamValid (const std::string &name) const
bool	isParamSetByUser (const std::string &nm) const
void	paramError (const std::string &param, Args... args) const
void	paramWarning (const std::string &param, Args... args) const
void	paramInfo (const std::string &param, Args... args) const
void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
void	mooseError (Args &&... args) const
void	mooseErrorNonPrefixed (Args &&... args) const
void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
void	mooseWarning (Args &&... args) const
void	mooseWarningNonPrefixed (Args &&... args) const
void	mooseDeprecated (Args &&... args) const
void	mooseInfo (Args &&... args) const
std::string	getDataFileName (const std::string &param) const
std::string	getDataFileNameByName (const std::string &relative_path) const
std::string	getDataFilePath (const std::string &relative_path) const
const std::vector< SubdomainName > &	blocks () const
unsigned int	numBlocks () const
virtual const std::set< SubdomainID > &	blockIDs () const
unsigned int	blocksMaxDimension () const
bool	hasBlocks (const SubdomainName &name) const
bool	hasBlocks (const std::vector< SubdomainName > &names) const
bool	hasBlocks (SubdomainID id) const
bool	hasBlocks (const std::vector< SubdomainID > &ids) const
bool	hasBlocks (const std::set< SubdomainID > &ids) const
bool	isBlockSubset (const std::set< SubdomainID > &ids) const
bool	isBlockSubset (const std::vector< SubdomainID > &ids) const
bool	hasBlockMaterialProperty (const std::string &prop_name)
const std::set< SubdomainID > &	meshBlockIDs () const
virtual bool	blockRestricted () const
virtual void	checkVariable (const MooseVariableFieldBase &variable) const
virtual const std::set< BoundaryID > &	boundaryIDs () const
const std::vector< BoundaryName > &	boundaryNames () const
unsigned int	numBoundaryIDs () const
bool	hasBoundary (const BoundaryName &name) const
bool	hasBoundary (const std::vector< BoundaryName > &names) const
bool	hasBoundary (const BoundaryID &id) const
bool	hasBoundary (const std::vector< BoundaryID > &ids, TEST_TYPE type=ALL) const
bool	hasBoundary (const std::set< BoundaryID > &ids, TEST_TYPE type=ALL) const
bool	isBoundarySubset (const std::set< BoundaryID > &ids) const
bool	isBoundarySubset (const std::vector< BoundaryID > &ids) const
bool	hasBoundaryMaterialProperty (const std::string &prop_name) const
virtual bool	boundaryRestricted () const
const std::set< BoundaryID > &	meshBoundaryIDs () const
virtual bool	checkVariableBoundaryIntegrity () const
virtual void	initialSetup ()
virtual void	timestepSetup ()
virtual void	jacobianSetup ()
virtual void	residualSetup ()
virtual void	customSetup (const ExecFlagType &)
const ExecFlagEnum &	getExecuteOnEnum () const
const std::set< MooseVariableFieldBase *> &	getMooseVariableDependencies () const
std::set< MooseVariableFieldBase *>	checkAllVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_omit={})
std::set< MooseVariableFieldBase *>	checkVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_check)
const std::vector< MooseVariableScalar *> &	getCoupledMooseScalarVars ()
const std::set< TagID > &	getScalarVariableCoupleableVectorTags () const
const std::set< TagID > &	getScalarVariableCoupleableMatrixTags () const
const Function &	getFunction (const std::string &name) const
const Function &	getFunctionByName (const FunctionName &name) const
bool	hasFunction (const std::string &param_name) const
bool	hasFunctionByName (const FunctionName &name) const
UserObjectName	getUserObjectName (const std::string &param_name) const
const T &	getUserObject (const std::string &param_name, bool is_dependency=true) const
const T &	getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
const UserObject &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const
const UserObject &	getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const
bool	isImplicit ()
Moose::StateArg	determineState () const
bool	isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
bool	hasPostprocessor (const std::string &param_name, const unsigned int index=0) const
bool	hasPostprocessorByName (const PostprocessorName &name) const
std::size_t	coupledPostprocessors (const std::string &param_name) const
const PostprocessorName &	getPostprocessorName (const std::string &param_name, const unsigned int index=0) const
const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
bool	hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
bool	hasVectorPostprocessor (const std::string &param_name) const
bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name) const
const VectorPostprocessorName &	getVectorPostprocessorName (const std::string &param_name) const
virtual void	meshChanged ()
void	buildOutputHideVariableList (std::set< std::string > variable_names)
void	setRandomResetFrequency (ExecFlagType exec_flag)
unsigned long	getRandomLong () const
Real	getRandomReal () const
unsigned int	getSeed (std::size_t id)
unsigned int	getMasterSeed () const
bool	isNodal () const
ExecFlagType	getResetOnTime () const
void	setRandomDataPointer (RandomData *random_data)
virtual unsigned int	getElementIDIndex (const std::string &id_parameter_name, unsigned int comp=0) const
virtual unsigned int	getElementIDIndexByName (const std::string &id_name) const
bool	hasElementID (const std::string &id_name) const
dof_id_type	maxElementID (unsigned int elem_id_index) const
dof_id_type	minElementID (unsigned int elem_id_index) const
bool	areElemIDsIdentical (const std::string &id_name1, const std::string &id_name2) const
std::unordered_map< dof_id_type, std::set< dof_id_type > >	getElemIDMapping (const std::string &id_name1, const std::string &id_name2) const
std::set< dof_id_type >	getAllElemIDs (unsigned int elem_id_index) const
std::set< dof_id_type >	getElemIDsOnBlocks (unsigned int elem_id_index, const std::set< SubdomainID > &blks) const
const std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> > &	getCoupledVars () const
const std::vector< MooseVariableFieldBase *> &	getCoupledMooseVars () const
const std::vector< MooseVariable *> &	getCoupledStandardMooseVars () const
const std::vector< VectorMooseVariable *> &	getCoupledVectorMooseVars () const
const std::vector< ArrayMooseVariable *> &	getCoupledArrayMooseVars () const
void	addFEVariableCoupleableVectorTag (TagID tag)
void	addFEVariableCoupleableMatrixTag (TagID tag)
std::set< TagID > &	getFEVariableCoupleableVectorTags ()
const std::set< TagID > &	getFEVariableCoupleableVectorTags () const
std::set< TagID > &	getFEVariableCoupleableMatrixTags ()
const std::set< TagID > &	getFEVariableCoupleableMatrixTags () const
auto &	getWritableCoupledVariables () const
bool	hasWritableCoupledVariables () const
const ADVariableValue *	getADDefaultValue (const std::string &var_name) const
const ADVectorVariableValue *	getADDefaultVectorValue (const std::string &var_name) const
const ADVariableGradient &	getADDefaultGradient () const
const ADVectorVariableGradient &	getADDefaultVectorGradient () const
const ADVariableSecond &	getADDefaultSecond () const
const ADVectorVariableCurl &	getADDefaultCurl () const
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)
const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name, MaterialData &material_data)
const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)
const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)
const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name, MaterialData &material_data)
const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)
const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)
const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)
const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)
const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const std::string &name, const unsigned int state=0)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const std::string &name, const unsigned int state=0)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialPropertyByName (const std::string &prop_name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialPropertyByName (const std::string &prop_name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)
const ADMaterialProperty< T > &	getADMaterialPropertyByName (const std::string &prop_name)
const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)
const ADMaterialProperty< T > &	getADMaterialPropertyByName (const std::string &prop_name)
const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)
const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)
const MaterialProperty< T > &	getMaterialPropertyOldByName (const std::string &prop_name)
const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)
const MaterialProperty< T > &	getMaterialPropertyOldByName (const std::string &prop_name)
const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)
const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)
const MaterialProperty< T > &	getMaterialPropertyOlderByName (const std::string &prop_name)
const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)
const MaterialProperty< T > &	getMaterialPropertyOlderByName (const std::string &prop_name)
const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > >	getBlockMaterialProperty (const MaterialPropertyName &name)
std::set< SubdomainID >	getMaterialPropertyBlocks (const std::string &name)
std::vector< SubdomainName >	getMaterialPropertyBlockNames (const std::string &name)
std::set< BoundaryID >	getMaterialPropertyBoundaryIDs (const std::string &name)
std::vector< BoundaryName >	getMaterialPropertyBoundaryNames (const std::string &name)
void	checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)
std::unordered_map< SubdomainID, std::vector< MaterialBase *> >	buildRequiredMaterials (bool allow_stateful=true)
void	statefulPropertiesAllowed (bool)
bool	getMaterialPropertyCalled () const
const GenericMaterialProperty< T, is_ad > &	getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)
const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)
const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)
const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)
const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)
const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)
const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)
MaterialProperty< T > &	declarePropertyByName (const std::string &prop_name)
MaterialProperty< T > &	declarePropertyByName (const std::string &prop_name)
MaterialProperty< T > &	declareProperty (const std::string &name)
MaterialProperty< T > &	declareProperty (const std::string &name)
ADMaterialProperty< T > &	declareADPropertyByName (const std::string &prop_name)
ADMaterialProperty< T > &	declareADPropertyByName (const std::string &prop_name)
ADMaterialProperty< T > &	declareADProperty (const std::string &name)
ADMaterialProperty< T > &	declareADProperty (const std::string &name)
auto &	declareGenericProperty (const std::string &prop_name)
auto &	declareGenericProperty (const std::string &prop_name)
GenericMaterialProperty< T, is_ad > &	declareGenericPropertyByName (const std::string &prop_name)
GenericMaterialProperty< T, is_ad > &	declareGenericPropertyByName (const std::string &prop_name)
const Distribution &	getDistribution (const std::string &name) const
const T &	getDistribution (const std::string &name) const
const Distribution &	getDistribution (const std::string &name) const
const T &	getDistribution (const std::string &name) const
const Distribution &	getDistributionByName (const DistributionName &name) const
const T &	getDistributionByName (const std::string &name) const
const Distribution &	getDistributionByName (const DistributionName &name) const
const T &	getDistributionByName (const std::string &name) const
bool	hasUserObject (const std::string &param_name) const
bool	hasUserObject (const std::string &param_name) const
bool	hasUserObject (const std::string &param_name) const
bool	hasUserObject (const std::string &param_name) const
bool	hasUserObjectByName (const UserObjectName &object_name) const
bool	hasUserObjectByName (const UserObjectName &object_name) const
bool	hasUserObjectByName (const UserObjectName &object_name) const
bool	hasUserObjectByName (const UserObjectName &object_name) const
const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const
virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const
const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const
const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const
const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const
const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const
bool	isVectorPostprocessorDistributed (const std::string &param_name) const
bool	isVectorPostprocessorDistributed (const std::string &param_name) const
bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
bool	hasMaterialProperty (const std::string &name)
bool	hasMaterialProperty (const std::string &name)
bool	hasMaterialPropertyByName (const std::string &name)
bool	hasMaterialPropertyByName (const std::string &name)
bool	hasADMaterialProperty (const std::string &name)
bool	hasADMaterialProperty (const std::string &name)
bool	hasADMaterialPropertyByName (const std::string &name)
bool	hasADMaterialPropertyByName (const std::string &name)
bool	hasGenericMaterialProperty (const std::string &name)
bool	hasGenericMaterialProperty (const std::string &name)
bool	hasGenericMaterialPropertyByName (const std::string &name)
bool	hasGenericMaterialPropertyByName (const std::string &name)
PenetrationLocator &	getPenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
PenetrationLocator &	getQuadraturePenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
NearestNodeLocator &	getNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)
NearestNodeLocator &	getQuadratureNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)
bool	requiresGeometricSearch () const
const Parallel::Communicator &	comm () const
processor_id_type	n_processors () const
processor_id_type	processor_id () const
void	resetQpProperties () final
	Retained as empty methods to avoid a warning from Material.C in framework. These methods are unused in all inheriting classes and should not be overwritten. More...
void	resetProperties () final

Static Public Member Functions
static InputParameters	validParams ()
static std::deque< MaterialBase *>	buildRequiredMaterials (const Consumers &mat_consumers, const std::vector< std::shared_ptr< MaterialBase >> &mats, const bool allow_stateful)
static bool	restricted (const std::set< BoundaryID > &ids)
static void	sort (typename std::vector< T > &vector)
static void	sortDFS (typename std::vector< T > &vector)
static void	cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)
static std::string	deduceFunctorName (const std::string &name, const InputParameters &params)

Public Attributes
unsigned int	_qp
const Elem *const &	_current_elem
Real &	_dt
const MooseArray< Point > &	_q_point
	ALL
	ANY
const ConsoleStream	_console

Static Public Attributes
static constexpr PropertyValue::id_type	default_property_id
static constexpr PropertyValue::id_type	zero_property_id

Protected Types
enum	QP_Data_Type

Protected Member Functions
virtual void	computeStressInitialize (const ADReal &effective_trial_stress, const ADRankFourTensor &elasticity_tensor) override
virtual ADReal	computeResidual (const ADReal &effective_trial_stress, const ADReal &scalar) override
virtual ADReal	computeDerivative (const ADReal &effective_trial_stress, const ADReal &scalar) override
std::size_t	stressIndex (const ADReal &effective_trial_stress)
virtual void	initQpStatefulProperties () override
virtual void	propagateQpStatefulProperties () override
	If updateState is not called during a timestep, this will be. More...
virtual void	computeStressFinalize (const GenericRankTwoTensor< is_ad > &plastic_strain_increment) override
	Perform any necessary steps to finalize state after return mapping iterations. More...
virtual Real	computeStressDerivative (const Real effective_trial_stress, const Real scalar) override
	This method returns the derivative of the creep strain with respect to the von mises stress. More...
template<>
Real	computeStressDerivative (const Real effective_trial_stress, const Real scalar)
	Calculate the derivative of the strain increment with respect to the updated stress. More...
virtual TangentCalculationMethod	getTangentCalculationMethod () override
void	propagateQpStatefulPropertiesRadialReturn ()
	Propagate the properties pertaining to this intermediate class. More...
virtual void	computeStressInitialize (const GenericReal< is_ad > &effective_trial_stress, const GenericRankFourTensor< is_ad > &elasticity_tensor)
	Perform any necessary initialization before return mapping iterations. More...
void	outputIterationSummary (std::stringstream *iter_output, const unsigned int total_it) override
	Output summary information for the convergence history of the model. More...
virtual void	checkMaterialProperty (const std::string &name, const unsigned int state) override
virtual const MaterialData &	materialData () const override
virtual MaterialData &	materialData () override
virtual const QBase &	qRule () const override
virtual void	computeQpProperties ()
virtual const FEProblemBase &	miProblem () const
virtual FEProblemBase &	miProblem ()
bool	isPropertyActive (const unsigned int prop_id) const
void	registerPropName (const std::string &prop_name, bool is_get, const unsigned int state)
void	checkExecutionStage ()
void	checkExecutionStage ()
virtual bool	hasBlockMaterialPropertyHelper (const std::string &prop_name)
void	initializeBlockRestrictable (const MooseObject *moose_object)
Moose::CoordinateSystemType	getBlockCoordSystem ()
bool	hasBoundaryMaterialPropertyHelper (const std::string &prop_name) const
void	addMooseVariableDependency (MooseVariableFieldBase *var)
void	addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)
bool	isCoupledScalar (const std::string &var_name, unsigned int i=0) const
unsigned int	coupledScalarComponents (const std::string &var_name) const
unsigned int	coupledScalar (const std::string &var_name, unsigned int comp=0) const
libMesh::Order	coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarValue (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< is_ad > &	coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
const VariableValue &	coupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
const VariableValue &	coupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
const VariableValue &	coupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDot (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
const MooseVariableScalar *	getScalarVar (const std::string &var_name, unsigned int comp) const
virtual void	addUserObjectDependencyHelper (const UserObject &) const
virtual void	addPostprocessorDependencyHelper (const PostprocessorName &) const
virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &) const
T &	declareRestartableData (const std::string &data_name, Args &&... args)
ManagedValue< T >	declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
const T &	getRestartableData (const std::string &data_name) const
T &	declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
T &	declareRecoverableData (const std::string &data_name, Args &&... args)
T &	declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
T &	declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
std::string	restartableName (const std::string &data_name) const
std::string	deduceFunctorName (const std::string &name) const
const Moose::Functor< T > &	getFunctor (const std::string &name)
const Moose::Functor< T > &	getFunctor (const std::string &name, THREAD_ID tid)
const Moose::Functor< T > &	getFunctor (const std::string &name, SubProblem &subproblem)
const Moose::Functor< T > &	getFunctor (const std::string &name, SubProblem &subproblem, THREAD_ID tid)
bool	isFunctor (const std::string &name) const
bool	isFunctor (const std::string &name, const SubProblem &subproblem) const
Moose::ElemArg	makeElemArg (const Elem *elem, bool correct_skewnewss=false) const
void	checkFunctorSupportsSideIntegration (const std::string &name, bool qp_integration)
void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) const
virtual void	coupledCallback (const std::string &, bool) const
virtual bool	isCoupled (const std::string &var_name, unsigned int i=0) const
virtual bool	isCoupledConstant (const std::string &var_name) const
unsigned int	coupledComponents (const std::string &var_name) const
VariableName	coupledName (const std::string &var_name, unsigned int comp=0) const
std::vector< VariableName >	coupledNames (const std::string &var_name) const
virtual unsigned int	coupled (const std::string &var_name, unsigned int comp=0) const
std::vector< unsigned int >	coupledIndices (const std::string &var_name) const
virtual const VariableValue &	coupledValue (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledValues (const std::string &var_name) const
std::vector< const VectorVariableValue *>	coupledVectorValues (const std::string &var_name) const
const GenericVariableValue< is_ad > &	coupledGenericValue (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericValue (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericValue (const std::string &var_name, unsigned int comp) const
std::vector< const GenericVariableValue< is_ad > *>	coupledGenericValues (const std::string &var_name) const
std::vector< const GenericVariableValue< false > *>	coupledGenericValues (const std::string &var_name) const
std::vector< const GenericVariableValue< true > *>	coupledGenericValues (const std::string &var_name) const
const GenericVariableValue< is_ad > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< is_ad > &	coupledGenericDot (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericDot (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericDot (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< is_ad > &	coupledGenericDotDot (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
virtual const VariableValue &	coupledValueLower (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adCoupledValue (const std::string &var_name, unsigned int comp=0) const
std::vector< const ADVariableValue *>	adCoupledValues (const std::string &var_name) const
const ADVariableValue &	adCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const
const ADVectorVariableValue &	adCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
std::vector< const ADVectorVariableValue *>	adCoupledVectorValues (const std::string &var_name) const
virtual const VariableValue &	coupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const VariableValue &	coupledVectorTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const VariableValue *>	coupledVectorTagValues (const std::string &var_names, TagID tag) const
std::vector< const VariableValue *>	coupledVectorTagValues (const std::string &var_names, const std::string &tag_name) const
virtual const ArrayVariableValue &	coupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const ArrayVariableValue &	coupledVectorTagArrayValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const ArrayVariableValue *>	coupledVectorTagArrayValues (const std::string &var_names, TagID tag) const
std::vector< const ArrayVariableValue *>	coupledVectorTagArrayValues (const std::string &var_names, const std::string &tag_name) const
virtual const VariableGradient &	coupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const VariableGradient &	coupledVectorTagGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const VariableGradient *>	coupledVectorTagGradients (const std::string &var_names, TagID tag) const
std::vector< const VariableGradient *>	coupledVectorTagGradients (const std::string &var_names, const std::string &tag_name) const
virtual const ArrayVariableGradient &	coupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const ArrayVariableGradient &	coupledVectorTagArrayGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const ArrayVariableGradient *>	coupledVectorTagArrayGradients (const std::string &var_names, TagID tag) const
std::vector< const ArrayVariableGradient *>	coupledVectorTagArrayGradients (const std::string &var_names, const std::string &tag_name) const
virtual const VariableValue &	coupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const
virtual const VariableValue &	coupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
const ArrayVariableValue &	coupledVectorTagArrayDofValue (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledVectorTagDofValues (const std::string &var_names, TagID tag) const
std::vector< const VariableValue *>	coupledVectorTagDofValues (const std::string &var_names, const std::string &tag_name) const
virtual const VariableValue &	coupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const VariableValue &	coupledMatrixTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const VariableValue *>	coupledMatrixTagValues (const std::string &var_names, TagID tag) const
std::vector< const VariableValue *>	coupledMatrixTagValues (const std::string &var_names, const std::string &tag_name) const
virtual const VectorVariableValue &	coupledVectorValue (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayValue (const std::string &var_name, unsigned int comp=0) const
std::vector< const ArrayVariableValue *>	coupledArrayValues (const std::string &var_name) const
MooseWritableVariable &	writableVariable (const std::string &var_name, unsigned int comp=0)
virtual VariableValue &	writableCoupledValue (const std::string &var_name, unsigned int comp=0)
void	checkWritableVar (MooseWritableVariable *var)
virtual const VariableValue &	coupledValueOld (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledValuesOld (const std::string &var_name) const
virtual const VariableValue &	coupledValueOlder (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledValuesOlder (const std::string &var_name) const
virtual const VariableValue &	coupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradient (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableGradient *>	coupledGradients (const std::string &var_name) const
const ADVariableGradient &	adCoupledGradient (const std::string &var_name, unsigned int comp=0) const
const ADVariableGradient &	adCoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
std::vector< const ADVariableGradient *>	adCoupledGradients (const std::string &var_name) const
const GenericVariableGradient< is_ad > &	coupledGenericGradient (const std::string &var_name, unsigned int comp=0) const
const GenericVariableGradient< false > &	coupledGenericGradient (const std::string &var_name, unsigned int comp) const
const GenericVariableGradient< true > &	coupledGenericGradient (const std::string &var_name, unsigned int comp) const
std::vector< const GenericVariableGradient< is_ad > *>	coupledGenericGradients (const std::string &var_name) const
std::vector< const GenericVariableGradient< false > *>	coupledGenericGradients (const std::string &var_name) const
std::vector< const GenericVariableGradient< true > *>	coupledGenericGradients (const std::string &var_name) const
const ADVectorVariableGradient &	adCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
const ADVariableSecond &	adCoupledSecond (const std::string &var_name, unsigned int comp=0) const
const ADVectorVariableSecond &	adCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradientOld (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableGradient *>	coupledGradientsOld (const std::string &var_name) const
virtual const VariableGradient &	coupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradientDot (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableGradient &	coupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableGradient &	coupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableGradient &	coupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableGradient &	coupledArrayGradient (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableGradient &	coupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableGradient &	coupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableGradient &	coupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableCurl &	coupledCurl (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableCurl &	coupledCurlOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableCurl &	coupledCurlOlder (const std::string &var_name, unsigned int comp=0) const
const ADVectorVariableCurl &	adCoupledCurl (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableDivergence &	coupledDiv (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableDivergence &	coupledDivOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableDivergence &	coupledDivOlder (const std::string &var_name, unsigned int comp=0) const
virtual const VariableSecond &	coupledSecond (const std::string &var_name, unsigned int comp=0) const
virtual const VariableSecond &	coupledSecondOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableSecond &	coupledSecondOlder (const std::string &var_name, unsigned int comp=0) const
virtual const VariableSecond &	coupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDot (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledDots (const std::string &var_name) const
virtual const VariableValue &	coupledDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDotDotOld (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adCoupledDot (const std::string &var_name, unsigned int comp=0) const
std::vector< const ADVariableValue *>	adCoupledDots (const std::string &var_name) const
const ADVariableValue &	adCoupledDotDot (const std::string &var_name, unsigned int comp=0) const
const ADVectorVariableValue &	adCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayDot (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDotDu (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDotDotDu (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledArrayDotDu (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalValue (const std::string &var_name, unsigned int comp=0) const
const Moose::ADType< T >::type &	adCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalDot (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledNodalDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledNodalDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledNodalDotDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDofValues (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledAllDofValues (const std::string &var_name) const
virtual const VariableValue &	coupledDofValuesOld (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledAllDofValuesOld (const std::string &var_name) const
virtual const VariableValue &	coupledDofValuesOlder (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledAllDofValuesOlder (const std::string &var_name) const
virtual const ArrayVariableValue &	coupledArrayDofValues (const std::string &var_name, unsigned int comp=0) const
virtual const ADVariableValue &	adCoupledDofValues (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adZeroValue () const
const ADVariableGradient &	adZeroGradient () const
const ADVariableSecond &	adZeroSecond () const
const GenericVariableValue< is_ad > &	genericZeroValue ()
const GenericVariableValue< false > &	genericZeroValue ()
const GenericVariableValue< true > &	genericZeroValue ()
const GenericVariableGradient< is_ad > &	genericZeroGradient ()
const GenericVariableGradient< false > &	genericZeroGradient ()
const GenericVariableGradient< true > &	genericZeroGradient ()
const GenericVariableSecond< is_ad > &	genericZeroSecond ()
const GenericVariableSecond< false > &	genericZeroSecond ()
const GenericVariableSecond< true > &	genericZeroSecond ()
bool	checkVar (const std::string &var_name, unsigned int comp=0, unsigned int comp_bound=0) const
const MooseVariableFieldBase *	getFEVar (const std::string &var_name, unsigned int comp) const
const MooseVariableFieldBase *	getFieldVar (const std::string &var_name, unsigned int comp) const
MooseVariableFieldBase *	getFieldVar (const std::string &var_name, unsigned int comp)
const T *	getVarHelper (const std::string &var_name, unsigned int comp) const
T *	getVarHelper (const std::string &var_name, unsigned int comp)
MooseVariable *	getVar (const std::string &var_name, unsigned int comp)
const MooseVariable *	getVar (const std::string &var_name, unsigned int comp) const
VectorMooseVariable *	getVectorVar (const std::string &var_name, unsigned int comp)
const VectorMooseVariable *	getVectorVar (const std::string &var_name, unsigned int comp) const
ArrayMooseVariable *	getArrayVar (const std::string &var_name, unsigned int comp)
const ArrayMooseVariable *	getArrayVar (const std::string &var_name, unsigned int comp) const
void	validateExecutionerType (const std::string &name, const std::string &fn_name) const
std::vector< T >	coupledVectorHelper (const std::string &var_name, const Func &func) const
void	markMatPropRequested (const std::string &)
MaterialPropertyName	getMaterialPropertyName (const std::string &name) const
const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)
const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)
const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)
const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)
const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)
const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)
void	returnMappingSolve (const GenericReal< is_ad > &effective_trial_stress, GenericReal< is_ad > &scalar, const ConsoleStream &console)
	Perform the return mapping iterations. More...
virtual GenericReal< is_ad >	initialGuess (const GenericReal< is_ad > &)
	Compute an initial guess for the value of the scalar. More...
virtual GenericReal< is_ad >	computeResidual (const GenericReal< is_ad > &, const GenericReal< is_ad > &)=0
	Compute the residual for a predicted value of the scalar. More...
virtual GenericReal< is_ad >	computeDerivative (const GenericReal< is_ad > &, const GenericReal< is_ad > &)=0
	Compute the derivative of the residual as a function of the scalar variable. More...
virtual GenericChainedReal< is_ad >	computeResidualAndDerivative (const GenericReal< is_ad > &, const GenericChainedReal< is_ad > &)
	Compute the residual and the derivative for a predicted value of the scalar. More...
virtual void	preStep (const GenericReal< is_ad > &, const GenericReal< is_ad > &, const GenericReal< is_ad > &)
	This method is called before taking a step in the return mapping algorithm. More...
virtual void	iterationFinalize (const GenericReal< is_ad > &)
	Finalize internal state variables for a model for a given iteration. More...
virtual void	outputIterationStep (std::stringstream *iter_output, const GenericReal< is_ad > &effective_trial_stress, const GenericReal< is_ad > &scalar, const Real reference_residual)
	Output information for a single iteration step to build the convergence history of the model. More...
bool	converged (const GenericReal< is_ad > &residual, const Real reference)
	Check to see whether the residual is within the convergence limits. More...

Protected Attributes
const ADVariableValue *const	_temperature
	Temperature variable value. More...
const std::vector< Real >	_coefficient
	Leading coefficient vector. More...
const std::vector< Real >	_n_exponent
	Exponent on the effective stress vector. More...
const std::vector< Real >	_m_exponent
	Exponent on time vector. More...
const std::vector< Real >	_stress_thresholds
	Stress thresholds vector. More...
const std::vector< Real >	_activation_energy
	Activation energy for exp term vector. More...
const Real	_gas_constant
	Gas constant for exp term. More...
const Real	_start_time
	Simulation start time. More...
ADReal	_exponential
	Exponential calculated from activation, gas constant, and temperature. More...
Real	_exp_time
	Exponential calculated from current time. More...
unsigned int	_stress_index
	Quadrature-point value pointing to the right power law parameters. More...
const unsigned int	_number_of_models
	Total number of models provided by the user. More...
GenericMaterialProperty< RankTwoTensor, is_ad > &	_creep_strain
	Creep strain material property. More...
const MaterialProperty< RankTwoTensor > &	_creep_strain_old
GenericReal< is_ad >	_three_shear_modulus
	3 * shear modulus More...
GenericMaterialProperty< Real, is_ad > &	_effective_inelastic_strain
const MaterialProperty< Real > &	_effective_inelastic_strain_old
GenericReal< is_ad >	_effective_inelastic_strain_increment
	Stores the scalar effective inelastic strain increment from Newton iteration. More...
Real	_max_inelastic_increment
	Maximum allowable scalar inelastic strain increment, used to control the timestep size in conjunction with a user object. More...
const Real	_substep_tolerance
	Used to calculate the number of substeps taken in the radial return algorithm, when substepping is enabled, based on the elastic strain increment ratio to the maximum inelastic increment. More...
const RankTwoTensor	_identity_two
	Rank two identity tensor. More...
const RankFourTensor	_identity_symmetric_four
	Rank four symmetric identity tensor. More...
const RankFourTensor	_deviatoric_projection_four
	Rank four deviatoric projection tensor. More...
const bool	_apply_strain
	Debugging option to enable specifying instead of calculating strain. More...
SubsteppingType	_use_substepping
	Whether user has requested the use of substepping technique to improve convergence [make const later]. More...
const bool	_adaptive_substepping
	Use adaptive substepping, cutting substep sizes until convergence is achieved. More...
const unsigned int	_maximum_number_substeps
	Maximum number of substeps. If the calculation results in a larger number, cut overall time step. More...
Real	_dt_original
	original timestep (to be restored after substepping is completed) More...
const std::string	_base_name
	Name used as a prefix for all material properties related to the stress update model. More...
	CURR
	PREV
bool	_bnd
bool	_neighbor
const QBase *const &	_qrule
const MooseArray< Real > &	_JxW
const SubdomainID &	_current_subdomain_id
const unsigned int &	_current_side
const ConstantTypeEnum	_constant_option
SubProblem &	_subproblem
FEProblemBase &	_fe_problem
THREAD_ID	_tid
Assembly &	_assembly
const MooseArray< Real > &	_coord
const MooseArray< Point > &	_normals
MooseMesh &	_mesh
const Moose::CoordinateSystemType &	_coord_sys
std::set< std::string >	_requested_props
std::set< std::string >	_supplied_props
std::set< unsigned int >	_supplied_prop_ids
std::unordered_set< unsigned int >	_active_prop_ids
const bool	_compute
std::unordered_map< unsigned int, unsigned int >	_props_to_min_states
std::vector< unsigned int >	_displacements
bool	_has_stateful_property
bool	_overrides_init_stateful_props
const FaceInfo *	_face_info
const bool &	_enabled
MooseApp &	_app
const std::string	_type
const std::string	_name
const InputParameters &	_pars
Factory &	_factory
ActionFactory &	_action_factory
const MaterialData *	_blk_material_data
const ExecFlagEnum &	_execute_enum
const ExecFlagType &	_current_execute_flag
FEProblemBase &	_sc_fe_problem
const THREAD_ID	_sc_tid
const Real &	_real_zero
const VariableValue &	_scalar_zero
const Point &	_point_zero
const InputParameters &	_ti_params
FEProblemBase &	_ti_feproblem
bool	_is_implicit
Real &	_t
const Real &	_t_old
int &	_t_step
Real &	_dt_old
bool	_is_transient
MooseApp &	_restartable_app
const std::string	_restartable_system_name
const THREAD_ID	_restartable_tid
const bool	_restartable_read_only
FEProblemBase &	_mci_feproblem
GeometricSearchData &	_geometric_search_data
bool	_requires_geometric_search
const InputParameters &	_c_parameters
const std::string &	_c_name
const std::string &	_c_type
FEProblemBase &	_c_fe_problem
const SystemBase *const	_c_sys
std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> >	_coupled_vars
std::vector< MooseVariableFieldBase *>	_coupled_moose_vars
std::vector< MooseVariable *>	_coupled_standard_moose_vars
std::vector< VectorMooseVariable *>	_coupled_vector_moose_vars
std::vector< ArrayMooseVariable *>	_coupled_array_moose_vars
std::vector< MooseVariableFV< Real > *>	_coupled_standard_fv_moose_vars
std::vector< MooseLinearVariableFV< Real > *>	_coupled_standard_linear_fv_moose_vars
const std::unordered_map< std::string, std::string > &	_new_to_deprecated_coupled_vars
bool	_c_nodal
bool	_c_is_implicit
const bool	_c_allow_element_to_nodal_coupling
THREAD_ID	_c_tid
std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > >	_default_value
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADReal > > >	_ad_default_value
std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > >	_default_vector_value
std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > >	_default_array_value
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > >	_ad_default_vector_value
VariableValue	_default_value_zero
VariableGradient	_default_gradient
MooseArray< ADRealVectorValue >	_ad_default_gradient
MooseArray< ADRealTensorValue >	_ad_default_vector_gradient
VariableSecond	_default_second
MooseArray< ADRealTensorValue >	_ad_default_second
MooseArray< ADRealVectorValue >	_ad_default_curl
const VariableValue &	_zero
const VariablePhiValue &	_phi_zero
const MooseArray< ADReal > &	_ad_zero
const VariableGradient &	_grad_zero
const MooseArray< ADRealVectorValue > &	_ad_grad_zero
const VariablePhiGradient &	_grad_phi_zero
const VariableSecond &	_second_zero
const MooseArray< ADRealTensorValue > &	_ad_second_zero
const VariablePhiSecond &	_second_phi_zero
const VectorVariableValue &	_vector_zero
const VectorVariableCurl &	_vector_curl_zero
VectorVariableValue	_default_vector_value_zero
VectorVariableGradient	_default_vector_gradient
VectorVariableCurl	_default_vector_curl
VectorVariableDivergence	_default_div
ArrayVariableValue	_default_array_value_zero
ArrayVariableGradient	_default_array_gradient
bool	_coupleable_neighbor
const InputParameters &	_mi_params
const std::string	_mi_name
const MooseObjectName	_mi_moose_object_name
FEProblemBase &	_mi_feproblem
SubProblem &	_mi_subproblem
const THREAD_ID	_mi_tid
const Moose::MaterialDataType	_material_data_type
MaterialData &	_material_data
bool	_stateful_allowed
bool	_get_material_property_called
std::vector< std::unique_ptr< PropertyValue > >	_default_properties
std::unordered_set< unsigned int >	_material_property_dependencies
const MaterialPropertyName	_get_suffix
const bool	_use_interpolated_state
const Parallel::Communicator &	_communicator
bool	_check_range
	Whether to check to see whether iterative solution is within admissible range, and set within that range if outside. More...
bool	_line_search
	Whether to use line searches to improve convergence. More...
bool	_bracket_solution
	Whether to save upper and lower bounds of root for scalar, and set solution to the midpoint between those bounds if outside them. More...

Static Protected Attributes
static RankFourTensor	_identityTensor
static const std::string	_interpolated_old
static const std::string	_interpolated_older

Detailed Description

This class uses the stress update material in a radial return isotropic creep model.

This class is one of the basic radial return constitutive models; more complex constitutive models combine creep and plasticity.

This class inherits from RadialReturnCreepStressUpdateBase and must be used in conjunction with ComputeMultipleInelasticStress. It calculates creep based on stress, temperature, and time effects. Particularly, it allows the user to provide multiple power law parameters which are chosen within this object depending on the von Mises stress level at the quadrature point.

Definition at line 25 of file ADMultiplePowerLawCreepStressUpdate.h.

Member Typedef Documentation

GR2

using StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::GR2 = GenericRankTwoTensor<is_ad>

inherited

Definition at line 57 of file StressUpdateBase.h.

GR4

using StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::GR4 = GenericRankFourTensor<is_ad>

inherited

Definition at line 59 of file StressUpdateBase.h.

GSR2

using StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::GSR2 = Moose::GenericType<RankTwoTensor , is_ad>

inherited

Definition at line 58 of file StressUpdateBase.h.

GSR4

using StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::GSR4 = Moose::GenericType<RankFourTensor , is_ad>

inherited

Definition at line 60 of file StressUpdateBase.h.

Member Enumeration Documentation

SubsteppingType

template<bool is_ad>

enum RadialReturnStressUpdateTempl::SubsteppingType

stronginherited

Enumerator
NONE
ERROR_BASED
INCREMENT_BASED

Definition at line 43 of file RadialReturnStressUpdate.h.

  {
    NONE,
    ERROR_BASED,
    INCREMENT_BASED
  };

Constructor & Destructor Documentation

ADMultiplePowerLawCreepStressUpdate()

ADMultiplePowerLawCreepStressUpdate::ADMultiplePowerLawCreepStressUpdate

(

const InputParameters &

parameters

)

Definition at line 40 of file ADMultiplePowerLawCreepStressUpdate.C.

  : ADRadialReturnCreepStressUpdateBase(parameters),
    _temperature(isParamValid("temperature") ? &adCoupledValue("temperature") : nullptr),
    _coefficient(getParam<std::vector<Real>>("coefficient")),
    _n_exponent(getParam<std::vector<Real>>("n_exponent")),
    _m_exponent(getParam<std::vector<Real>>("m_exponent")),
    _stress_thresholds(getParam<std::vector<Real>>("stress_thresholds")),
    _activation_energy(getParam<std::vector<Real>>("activation_energy")),
    _gas_constant(getParam<Real>("gas_constant")),
    _start_time(getParam<Real>("start_time")),
    _exponential(1),
    _exp_time(0),
    _stress_index(0),
    _number_of_models(_m_exponent.size())
{

  if (!std::is_sorted(_stress_thresholds.begin(), _stress_thresholds.end()))
    paramError("stress_thresholds",
               "Stress thresholds input must be provided in increasing ordered");

  if (_coefficient.size() != _n_exponent.size() || _coefficient.size() != _m_exponent.size() ||
      _coefficient.size() != _activation_energy.size())
    paramError(
        "n_exponent",
        "Inputs to ADMultiplePowerLawCreepStressUpdate creep models must have the same size");

  if (_number_of_models != _stress_thresholds.size() - 1)
    paramError("stress_thresholds",
               "The number of creep models must be the number of stress thresholds minute. That "
               "is, one creep model will be valid between two stress thresholds");

  for (unsigned int i = 0; i < _number_of_models; i++)
    if (_start_time < this->_app.getStartTime() && (std::trunc(_m_exponent[i]) != _m_exponent[i]))
      paramError("start_time",
                 "Start time must be equal to or greater than the Executioner start_time if a "
                 "non-integer m_exponent is used");
}

Member Function Documentation

calculateNumberSubsteps() [1/2]

template<bool is_ad>

int RadialReturnStressUpdateTempl< is_ad >::calculateNumberSubsteps ( const GenericRankTwoTensor< is_ad > &  strain_increment )

overridevirtualinherited

If substepping is enabled, calculate the number of substeps as a function of the elastic strain increment guess and the maximum inelastic strain increment ratio based on a user-specified tolerance.

Parameters

strain_increment

When called, this is the elastic strain guess

Returns

The number of substeps required

Definition at line 150 of file RadialReturnStressUpdate.C.

{
  // compute an effective elastic strain measure
  const GenericReal<is_ad> contracted_elastic_strain =
      strain_increment.doubleContraction(strain_increment);
  const Real effective_elastic_strain =
      std::sqrt(3.0 / 2.0 * MetaPhysicL::raw_value(contracted_elastic_strain));

  if (MooseUtils::absoluteFuzzyEqual(effective_elastic_strain, 0.0))
    return 1;

  if (_use_substepping == SubsteppingType::INCREMENT_BASED)
  {
    const Real ratio = effective_elastic_strain / _max_inelastic_increment;

    if (ratio > _substep_tolerance)
      return std::ceil(ratio / _substep_tolerance);
    return 1;
  }

  if (_use_substepping == SubsteppingType::ERROR_BASED)
  {
    const Real accurate_time_step_ratio = _substep_tolerance / effective_elastic_strain;

    if (accurate_time_step_ratio < 1.0)
      return std::ceil(1.0 / accurate_time_step_ratio);
    return 1;
  }

  mooseError("calculateNumberSubsteps should not have been called. Notify a developer.");
}

calculateNumberSubsteps() [2/2]

virtual int StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::calculateNumberSubsteps ( const GR2 &  )

inlinevirtualinherited

Given the elastic strain increment compute the number of substeps required to bring a substepped trial stress guess distance from the yield surface into the tolerance specified in the individual child class.

Definition at line 160 of file StressUpdateBase.h.

{ return 1; }

computeDerivative() [1/2]

ADReal ADMultiplePowerLawCreepStressUpdate::computeDerivative ( const ADReal &  effective_trial_stress, const ADReal &  scalar  )

overrideprotectedvirtual

Definition at line 106 of file ADMultiplePowerLawCreepStressUpdate.C.

{
  const ADReal stress_delta = effective_trial_stress - _three_shear_modulus * scalar;
  const ADReal creep_rate_derivative =
      -_coefficient[_stress_index] * _three_shear_modulus * _n_exponent[_stress_index] *
      std::pow(stress_delta, _n_exponent[_stress_index] - 1.0) * _exponential * _exp_time;
  return creep_rate_derivative * _dt - 1.0;
}

computeDerivative() [2/2]

template<bool is_ad>

virtual GenericReal<is_ad> SingleVariableReturnMappingSolutionTempl< is_ad >::computeDerivative ( const GenericReal< is_ad > &  , const GenericReal< is_ad > &    )

protectedpure virtualinherited

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

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

Parameters

effective_trial_stress	Effective trial stress
scalar	Inelastic strain increment magnitude being solved for

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

computeReferenceResidual()

template<bool is_ad>

Real RadialReturnStressUpdateTempl< is_ad >::computeReferenceResidual ( const GenericReal< is_ad > &  effective_trial_stress, const GenericReal< is_ad > &  scalar  )

overridevirtualinherited

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

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

Parameters

effective_trial_stress	Effective trial stress
scalar	Inelastic strain increment magnitude being solved for

Implements SingleVariableReturnMappingSolutionTempl< is_ad >.

Definition at line 468 of file RadialReturnStressUpdate.C.

{
  return MetaPhysicL::raw_value(effective_trial_stress / _three_shear_modulus) -
         MetaPhysicL::raw_value(scalar_effective_inelastic_strain);
}

computeResidual() [1/2]

ADReal ADMultiplePowerLawCreepStressUpdate::computeResidual ( const ADReal &  effective_trial_stress, const ADReal &  scalar  )

overrideprotectedvirtual

Definition at line 95 of file ADMultiplePowerLawCreepStressUpdate.C.

{
  const ADReal stress_delta = effective_trial_stress - _three_shear_modulus * scalar;
  const ADReal creep_rate = _coefficient[_stress_index] *
                            std::pow(stress_delta, _n_exponent[_stress_index]) * _exponential *
                            _exp_time;
  return creep_rate * _dt - scalar;
}

computeResidual() [2/2]

template<bool is_ad>

virtual GenericReal<is_ad> SingleVariableReturnMappingSolutionTempl< is_ad >::computeResidual ( const GenericReal< is_ad > &  , const GenericReal< is_ad > &    )

protectedpure virtualinherited

Compute the residual for a predicted value of the scalar.

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

Parameters

effective_trial_stress	Effective trial stress
scalar	Inelastic strain increment magnitude being solved for

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

computeResidualAndDerivative()

template<bool is_ad>

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

inlineprotectedvirtualinherited

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

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

Parameters

effective_trial_stress	Effective trial stress
scalar	Inelastic strain increment magnitude being solved for

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

Definition at line 97 of file SingleVariableReturnMappingSolution.h.

  {
    mooseError("computeResidualAndDerivative has to be implemented if "
               "automatic_differentiation_return_mapping = true.");
    return 0;
  };

computeStrainEnergyRateDensity() [1/2]

Real ADMultiplePowerLawCreepStressUpdate::computeStrainEnergyRateDensity ( const ADMaterialProperty< RankTwoTensor > &  stress, const ADMaterialProperty< RankTwoTensor > &  strain_rate  )

overridevirtual

Definition at line 117 of file ADMultiplePowerLawCreepStressUpdate.C.

{
  if (_n_exponent[_stress_index] <= 1)
    return 0.0;

  Real creep_factor = _n_exponent[_stress_index] / (_n_exponent[_stress_index] + 1);

  return MetaPhysicL::raw_value(creep_factor * stress[_qp].doubleContraction((strain_rate)[_qp]));
}

computeStrainEnergyRateDensity() [2/2]

virtual Real StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::computeStrainEnergyRateDensity ( const GenericMaterialProperty< RankTwoTensor, is_ad > &  , const GenericMaterialProperty< RankTwoTensor, is_ad > &    )

inlinevirtualinherited

Compute the strain energy rate density for this inelastic model for the current step.

Parameters

stress	The stress tensor at the end of the step
strain_rate	The strain rate at the end of the step

Returns

The computed strain energy rate density

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

Definition at line 180 of file StressUpdateBase.h.

  {
    mooseError(
        "The computation of strain energy rate density needs to be implemented by a child class");
    return 0.0;
  }

computeStressDerivative() [1/2]

template<bool is_ad>

Real RadialReturnCreepStressUpdateBaseTempl< is_ad >::computeStressDerivative ( const Real  effective_trial_stress, const Real  scalar  )

overrideprotectedvirtualinherited

This method returns the derivative of the creep strain with respect to the von mises stress.

It assumes the stress delta (von mises stress used to determine the creep rate) is calculated as: effective_trial_stress - _three_shear_modulus * scalar

Reimplemented from RadialReturnStressUpdateTempl< is_ad >.

Definition at line 52 of file RadialReturnCreepStressUpdateBase.C.

{
  mooseError("computeStressDerivative called: no stress derivative computation is needed for AD");
}

computeStressDerivative() [2/2]

template<>

Real RadialReturnCreepStressUpdateBaseTempl< false >::computeStressDerivative ( const Real  , const Real    )

protectedvirtualinherited

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

Parameters

effective_trial_stress	Effective trial stress
scalar	Inelastic strain increment magnitude being solved for

Reimplemented from RadialReturnStressUpdateTempl< is_ad >.

Definition at line 60 of file RadialReturnCreepStressUpdateBase.C.

{
  return -(computeDerivative(effective_trial_stress, scalar) + 1.0) / this->_three_shear_modulus;
}

computeStressFinalize()

template<bool is_ad>

void RadialReturnCreepStressUpdateBaseTempl< is_ad >::computeStressFinalize ( const GenericRankTwoTensor< is_ad > &  )

overrideprotectedvirtualinherited

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

Parameters

inelasticStrainIncrement

Inelastic strain increment

Reimplemented from RadialReturnStressUpdateTempl< is_ad >.

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

Definition at line 68 of file RadialReturnCreepStressUpdateBase.C.

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

{
  _creep_strain[_qp] = _creep_strain_old[_qp] + plastic_strain_increment;
}

computeStressInitialize() [1/2]

void ADMultiplePowerLawCreepStressUpdate::computeStressInitialize ( const ADReal &  effective_trial_stress, const ADRankFourTensor &  elasticity_tensor  )

overrideprotectedvirtual

Definition at line 80 of file ADMultiplePowerLawCreepStressUpdate.C.

{
  ADRadialReturnStressUpdate::computeStressInitialize(effective_trial_stress, elasticity_tensor);

  _stress_index = stressIndex(effective_trial_stress);

  if (_temperature)
    _exponential =
        std::exp(-_activation_energy[_stress_index] / (_gas_constant * (*_temperature)[_qp]));

  _exp_time = std::pow(_t - _start_time, _m_exponent[_stress_index]);
}

computeStressInitialize() [2/2]

template<bool is_ad>

void RadialReturnStressUpdateTempl< is_ad >::computeStressInitialize ( const GenericReal< is_ad > &  effective_trial_stress, const GenericRankFourTensor< is_ad > &  elasticity_tensor  )

protectedvirtualinherited

Perform any necessary initialization before return mapping iterations.

Parameters

effective_trial_stress	Effective trial stress
elasticityTensor	Elasticity tensor

Reimplemented in LAROMANCEStressUpdateBaseTempl< is_ad >, IsotropicPlasticityStressUpdateTempl< is_ad >, IsotropicPowerLawHardeningStressUpdateTempl< is_ad >, CompositePowerLawCreepStressUpdateTempl< is_ad >, PowerLawCreepStressUpdateTempl< is_ad >, TemperatureDependentHardeningStressUpdateTempl< is_ad >, and CombinedNonlinearHardeningPlasticityTempl< is_ad >.

Definition at line 133 of file RadialReturnStressUpdate.C.

Referenced by CombinedNonlinearHardeningPlasticityTempl< is_ad >::computeStressInitialize(), TemperatureDependentHardeningStressUpdateTempl< is_ad >::computeStressInitialize(), computeStressInitialize(), PowerLawCreepStressUpdateTempl< is_ad >::computeStressInitialize(), CompositePowerLawCreepStressUpdateTempl< is_ad >::computeStressInitialize(), IsotropicPlasticityStressUpdateTempl< is_ad >::computeStressInitialize(), IsotropicPowerLawHardeningStressUpdateTempl< is_ad >::computeStressInitialize(), and LAROMANCEStressUpdateBaseTempl< is_ad >::computeStressInitialize().

{
  // Set the value of 3 * shear modulus for use as a reference residual value
  _three_shear_modulus = 3.0 * ElasticityTensorTools::getIsotropicShearModulus(elasticity_tensor);
}

computeTangentOperator() [1/3]

template<bool is_ad>

void RadialReturnStressUpdateTempl< is_ad >::computeTangentOperator ( Real  effective_trial_stress, const RankTwoTensor &  stress_new, RankFourTensor &  tangent_operator  )

inherited

Calculate the tangent_operator.

Definition at line 185 of file RadialReturnStressUpdate.C.

Referenced by ComputeCreepPlasticityStress::computeTangentOperators().

{
  mooseError("computeTangentOperator called: no tangent computation is needed for AD");
}

computeTangentOperator() [2/3]

template<>

void RadialReturnStressUpdateTempl< false >::computeTangentOperator ( Real  effective_trial_stress, const RankTwoTensor &  stress_new, RankFourTensor &  tangent_operator  )

inherited

Definition at line 194 of file RadialReturnStressUpdate.C.

{
  if (getTangentCalculationMethod() == TangentCalculationMethod::PARTIAL)
  {
    if (MooseUtils::absoluteFuzzyEqual(_effective_inelastic_strain_increment, 0.0))
      tangent_operator.zero();
    else
    {
      // mu = _three_shear_modulus / 3.0;
      // norm_dev_stress = ||s_n+1||
      // effective_trial_stress = von mises trial stress = std::sqrt(3.0 / 2.0) * ||s_n+1^trial||
      // scalar_effective_inelastic_strain = Delta epsilon^cr_n+1
      // deriv = derivative of scalar_effective_inelastic_strain w.r.t. von mises stress
      // deriv = std::sqrt(3.0 / 2.0) partial Delta epsilon^cr_n+1n over partial ||s_n+1^trial||

      mooseAssert(_three_shear_modulus != 0.0, "Shear modulus is zero");

      const RankTwoTensor deviatoric_stress = stress_new.deviatoric();
      const Real norm_dev_stress_squared = deviatoric_stress.doubleContraction(deviatoric_stress);
      if (MooseUtils::absoluteFuzzyEqual(norm_dev_stress_squared, 0.0))
      {
        tangent_operator.zero();
        return;
      }
      const Real norm_dev_stress = std::sqrt(norm_dev_stress_squared);

      const RankTwoTensor flow_direction = deviatoric_stress / norm_dev_stress;
      const RankFourTensor flow_direction_dyad = flow_direction.outerProduct(flow_direction);
      const Real deriv =
          computeStressDerivative(effective_trial_stress, _effective_inelastic_strain_increment);
      const Real scalar_one = _three_shear_modulus * _effective_inelastic_strain_increment /
                              std::sqrt(1.5) / norm_dev_stress;

      tangent_operator = scalar_one * _deviatoric_projection_four +
                         (_three_shear_modulus * deriv - scalar_one) * flow_direction_dyad;
    }
  }
}

computeTangentOperator() [3/3]

template<>

void RadialReturnStressUpdateTempl< false >::computeTangentOperator ( Real  effective_trial_stress, const RankTwoTensor &  stress_new, RankFourTensor &  tangent_operator  )

inherited

computeTimeStepLimit()

template<bool is_ad>

Real RadialReturnStressUpdateTempl< is_ad >::computeTimeStepLimit ( )

overridevirtualinherited

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

Returns

Limiting time step

Reimplemented from StressUpdateBaseTempl< is_ad >.

Reimplemented in LAROMANCEStressUpdateBaseTempl< is_ad >.

Definition at line 486 of file RadialReturnStressUpdate.C.

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

{
  const Real scalar_inelastic_strain_incr =
      std::abs(MetaPhysicL::raw_value(_effective_inelastic_strain[_qp]) -
               _effective_inelastic_strain_old[_qp]);
  if (!scalar_inelastic_strain_incr)
    return std::numeric_limits<Real>::max();

  return _dt * _max_inelastic_increment / scalar_inelastic_strain_incr;
}

converged()

template<bool is_ad>

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

protectedinherited

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

Parameters

residual	Current value of the residual
reference	Current value of the reference quantity

Returns

Whether the model converged

Definition at line 326 of file SingleVariableReturnMappingSolution.C.

{
  const Real residual = MetaPhysicL::raw_value(ad_residual);
  return (std::abs(residual) <= _absolute_tolerance ||
          std::abs(residual / reference) <= _relative_tolerance);
}

effectiveInelasticStrainIncrement()

template<bool is_ad>

const GenericReal<is_ad>& RadialReturnStressUpdateTempl< is_ad >::effectiveInelasticStrainIncrement ( ) const

inlineinherited

Current value of scalar inelastic strain.

Definition at line 151 of file RadialReturnStressUpdate.h.

Referenced by ComputeCreepPlasticityStress::updateQpState().

  {
    return _effective_inelastic_strain_increment;
  }

getTangentCalculationMethod()

template<bool is_ad>

virtual TangentCalculationMethod RadialReturnCreepStressUpdateBaseTempl< is_ad >::getTangentCalculationMethod ( )

inlineoverrideprotectedvirtualinherited

Reimplemented from StressUpdateBaseTempl< is_ad >.

Definition at line 48 of file RadialReturnCreepStressUpdateBase.h.

  {
    return TangentCalculationMethod::PARTIAL;
  }

initialGuess()

template<bool is_ad>

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

inlineprotectedvirtualinherited

Compute an initial guess for the value of the scalar.

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

Parameters

effective_trial_stress

Effective trial stress

Reimplemented in PowerLawCreepTestTempl< is_ad >.

Definition at line 66 of file SingleVariableReturnMappingSolution.h.

  {
    return 0.0;
  }

initQpStatefulProperties()

template<bool is_ad>

void RadialReturnCreepStressUpdateBaseTempl< is_ad >::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Reimplemented from RadialReturnStressUpdateTempl< is_ad >.

Reimplemented in LAROMANCEStressUpdateBaseTempl< is_ad >.

Definition at line 34 of file RadialReturnCreepStressUpdateBase.C.

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

{
  _creep_strain[_qp].zero();

  RadialReturnStressUpdateTempl<is_ad>::initQpStatefulProperties();
}

isIsotropic()

template<bool is_ad>

bool RadialReturnStressUpdateTempl< is_ad >::isIsotropic ( )

inlineoverridevirtualinherited

Radial return mapped models should be isotropic by default!

Reimplemented from StressUpdateBaseTempl< is_ad >.

Definition at line 128 of file RadialReturnStressUpdate.h.

{ return true; };

iterationFinalize()

template<bool is_ad>

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

inlineprotectedvirtualinherited

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

Parameters

scalar

Inelastic strain increment magnitude being solved for

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

Definition at line 128 of file SingleVariableReturnMappingSolution.h.

{}

maximumPermissibleValue()

template<bool is_ad>

GenericReal< is_ad > RadialReturnStressUpdateTempl< is_ad >::maximumPermissibleValue ( const GenericReal< is_ad > &  effective_trial_stress ) const

overridevirtualinherited

Compute the maximum permissible value of the scalar.

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

Parameters

effective_trial_stress

Effective trial stress

Reimplemented from SingleVariableReturnMappingSolutionTempl< is_ad >.

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

Definition at line 478 of file RadialReturnStressUpdate.C.

{
  return effective_trial_stress / _three_shear_modulus;
}

minimumPermissibleValue()

template<bool is_ad>

virtual GenericReal<is_ad> RadialReturnStressUpdateTempl< is_ad >::minimumPermissibleValue ( const GenericReal< is_ad > &  effective_trial_stress ) const

inlineoverridevirtualinherited

Compute the minimum permissible value of the scalar.

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

Parameters

effective_trial_stress

Effective trial stress

Reimplemented from SingleVariableReturnMappingSolutionTempl< is_ad >.

Reimplemented in PowerLawCreepTestTempl< is_ad >.

Definition at line 106 of file RadialReturnStressUpdate.h.

  {
    return 0.0;
  }

outputIterationStep()

template<bool is_ad>

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

protectedvirtualinherited

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

Parameters

iter_output	Output stream
effective_trial_stress	Effective trial stress
residual	Current value of the residual
reference	Current value of the reference quantity

Reimplemented in LAROMANCEStressUpdateBaseTempl< is_ad >.

Definition at line 419 of file SingleVariableReturnMappingSolution.C.

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

{
  if (iter_output)
  {
    const unsigned int it = _iteration;
    const Real residual = MetaPhysicL::raw_value(_residual);

    *iter_output << " iteration=" << it
                 << " trial_stress=" << MetaPhysicL::raw_value(effective_trial_stress)
                 << " scalar=" << MetaPhysicL::raw_value(scalar) << " residual=" << residual
                 << " ref_res=" << reference_residual
                 << " rel_res=" << std::abs(residual) / reference_residual
                 << " rel_tol=" << _relative_tolerance << " abs_res=" << std::abs(residual)
                 << " abs_tol=" << _absolute_tolerance << '\n';
  }
}

outputIterationSummary()

template<bool is_ad>

void RadialReturnStressUpdateTempl< is_ad >::outputIterationSummary ( std::stringstream *  iter_output, const unsigned int  total_it  )

overrideprotectedvirtualinherited

Output summary information for the convergence history of the model.

Parameters

iter_output	Output stream
total_it	Total iteration count

Reimplemented from SingleVariableReturnMappingSolutionTempl< is_ad >.

Definition at line 499 of file RadialReturnStressUpdate.C.

{
  if (iter_output)
  {
    *iter_output << "At element " << _current_elem->id() << " _qp=" << _qp << " Coordinates "
                 << _q_point[_qp] << " block=" << _current_elem->subdomain_id() << '\n';
  }
  SingleVariableReturnMappingSolutionTempl<is_ad>::outputIterationSummary(iter_output, total_it);
}

preStep()

template<bool is_ad>

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

inlineprotectedvirtualinherited

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

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

Definition at line 118 of file SingleVariableReturnMappingSolution.h.

  {
  }

propagateQpStatefulProperties()

template<bool is_ad>

void RadialReturnCreepStressUpdateBaseTempl< is_ad >::propagateQpStatefulProperties ( )

overrideprotectedvirtualinherited

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

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

Reimplemented from StressUpdateBaseTempl< is_ad >.

Definition at line 43 of file RadialReturnCreepStressUpdateBase.C.

{
  _creep_strain[_qp] = _creep_strain_old[_qp];

  propagateQpStatefulPropertiesRadialReturn();
}

propagateQpStatefulPropertiesRadialReturn()

template<bool is_ad>

void RadialReturnStressUpdateTempl< is_ad >::propagateQpStatefulPropertiesRadialReturn ( )

protectedinherited

Propagate the properties pertaining to this intermediate class.

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

Definition at line 143 of file RadialReturnStressUpdate.C.

Referenced by CombinedNonlinearHardeningPlasticityTempl< is_ad >::propagateQpStatefulProperties(), and IsotropicPlasticityStressUpdateTempl< is_ad >::propagateQpStatefulProperties().

{
  _effective_inelastic_strain[_qp] = _effective_inelastic_strain_old[_qp];
}

requiresIsotropicTensor()

template<bool is_ad>

bool RadialReturnStressUpdateTempl< is_ad >::requiresIsotropicTensor ( )

inlineoverridevirtualinherited

Does the model require the elasticity tensor to be isotropic?

Implements StressUpdateBaseTempl< is_ad >.

Definition at line 123 of file RadialReturnStressUpdate.h.

{ return true; }

resetIncrementalMaterialProperties()

virtual void StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::resetIncrementalMaterialProperties ( )

inlinevirtualinherited

Reset material properties.

Useful for substepping with inelastic models.

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

Definition at line 172 of file StressUpdateBase.h.

{};

resetProperties()

void StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::resetProperties ( )

inlinefinalvirtualinherited

Reimplemented from Material.

Definition at line 138 of file StressUpdateBase.h.

{}

resetQpProperties()

void StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::resetQpProperties ( )

inlinefinalvirtualinherited

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

Reimplemented from Material.

Definition at line 137 of file StressUpdateBase.h.

{}

returnMappingSolve()

template<bool is_ad>

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

protectedinherited

Perform the return mapping iterations.

Parameters

effective_trial_stress	Effective trial stress
scalar	Inelastic strain increment magnitude being solved for
console	Console output

Definition at line 100 of file SingleVariableReturnMappingSolution.C.

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

{
  // construct the stringstream here only if the debug level is set to ALL
  std::unique_ptr<std::stringstream> iter_output =
      (_internal_solve_output_on == InternalSolveOutput::ALWAYS)
          ? std::make_unique<std::stringstream>()
          : nullptr;

  // do the internal solve and capture iteration info during the first round
  // iff full history output is requested regardless of whether the solve failed or succeeded
  auto solve_state =
      internalSolve(effective_trial_stress,
                    scalar,
                    _internal_solve_full_iteration_history ? iter_output.get() : nullptr);
  if (solve_state != SolveState::SUCCESS &&
      _internal_solve_output_on != InternalSolveOutput::ALWAYS)
  {
    // output suppressed by user, throw immediately
    if (_internal_solve_output_on == InternalSolveOutput::NEVER)
      mooseException("");

    // user expects some kind of output, if necessary setup output stream now
    if (!iter_output)
      iter_output = std::make_unique<std::stringstream>();

    // add the appropriate error message to the output
    switch (solve_state)
    {
      case SolveState::NAN_INF:
        *iter_output << "Encountered inf or nan in material return mapping iterations.\n";
        break;

      case SolveState::EXCEEDED_ITERATIONS:
        *iter_output << "Exceeded maximum iterations in material return mapping iterations.\n";
        break;

      default:
        mooseError("Unhandled solver state");
    }

    // if full history output is only requested for failed solves we have to repeat
    // the solve a second time
    if (_internal_solve_full_iteration_history)
      internalSolve(effective_trial_stress, scalar, iter_output.get());

    // Append summary and throw exception
    outputIterationSummary(iter_output.get(), _iteration);
    mooseException(iter_output->str());
  }

  if (_internal_solve_output_on == InternalSolveOutput::ALWAYS)
  {
    // the solve did not fail but the user requested debug output anyways
    outputIterationSummary(iter_output.get(), _iteration);
    console << iter_output->str() << std::flush;
  }
}

setQp()

void StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::setQp ( unsigned int  qp )

inherited

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

Definition at line 46 of file StressUpdateBase.C.

{
  _qp = qp;
}

storeIncrementalMaterialProperties()

virtual void StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::storeIncrementalMaterialProperties ( const unsigned int  )

inlinevirtualinherited

Properly set up the incremental calculation storage of the stateful material properties in the inheriting classes.

Reimplemented in LAROMANCEStressUpdateBaseTempl< is_ad >.

Definition at line 167 of file StressUpdateBase.h.

{};

stressIndex()

std::size_t ADMultiplePowerLawCreepStressUpdate::stressIndex ( const ADReal &  effective_trial_stress )

protected

Definition at line 136 of file ADMultiplePowerLawCreepStressUpdate.C.

Referenced by computeStressInitialize().

{
  // Check the correct model for *this* radial return
  std::size_t i = 0;
  while (i < _number_of_models - 1 && effective_trial_stress > _stress_thresholds[i + 1])
    ++i;

  return i;
}

substeppingCapabilityEnabled()

bool ADMultiplePowerLawCreepStressUpdate::substeppingCapabilityEnabled ( )

overridevirtual

Does the model include the infrastructure for substep decomposition of the elastic strain initially used to calculate the trial stress guess Inheriting classes which wish to use the substepping capability should overwrite this method and set it to return true.

Reimplemented from StressUpdateBaseTempl< is_ad >.

Definition at line 130 of file ADMultiplePowerLawCreepStressUpdate.C.

{
  return this->_use_substepping != ADRadialReturnStressUpdate::SubsteppingType::NONE;
}

substeppingCapabilityRequested()

template<bool is_ad>

virtual bool RadialReturnStressUpdateTempl< is_ad >::substeppingCapabilityRequested ( )

inlineoverridevirtualinherited

Has the user requested usage of (possibly) implemented substepping capability for inelastic models.

Parent classes set this to false, but RadialReturn inelastic models have the ability to implement substepping.

Reimplemented from StressUpdateBaseTempl< is_ad >.

Definition at line 145 of file RadialReturnStressUpdate.h.

  {
    return _use_substepping != SubsteppingType::NONE;
  }

updateEffectiveInelasticStrain()

template<bool is_ad>

void RadialReturnStressUpdateTempl< is_ad >::updateEffectiveInelasticStrain ( const GenericReal< is_ad > &  increment )

inlineinherited

Definition at line 161 of file RadialReturnStressUpdate.h.

Referenced by ComputeCreepPlasticityStress::finalizeConstitutiveModels().

  {
    _effective_inelastic_strain[_qp] = _effective_inelastic_strain_old[_qp] + increment;
  }

updateEffectiveInelasticStrainIncrement()

template<bool is_ad>

void RadialReturnStressUpdateTempl< is_ad >::updateEffectiveInelasticStrainIncrement ( const GenericReal< is_ad > &  eisi )

inlineinherited

Definition at line 156 of file RadialReturnStressUpdate.h.

Referenced by ComputeCreepPlasticityStress::finalizeConstitutiveModels().

  {
    _effective_inelastic_strain_increment = eisi;
  }

updateState() [1/2]

template<bool is_ad>

void RadialReturnStressUpdateTempl< is_ad >::updateState ( GenericRankTwoTensor< is_ad > &  strain_increment, GenericRankTwoTensor< is_ad > &  inelastic_strain_increment, const GenericRankTwoTensor< is_ad > &  rotation_increment, GenericRankTwoTensor< is_ad > &  stress_new, const RankTwoTensor &  stress_old, const GenericRankFourTensor< is_ad > &  elasticity_tensor, const RankTwoTensor &  elastic_strain_old, bool  compute_full_tangent_operator = false, RankFourTensor &  tangent_operator = StressUpdateBaseTempl<is_ad>::_identityTensor  )

overridevirtualinherited

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

Parameters

strain_increment	Sum of elastic and inelastic strain increments
inelastic_strain_increment	Inelastic strain increment calculated by this class
rotation	increment Not used by this class
stress_new	New trial stress from pure elastic calculation
stress_old	Old state of stress
elasticity_tensor	Rank 4 C_{ijkl}, must be isotropic
elastic_strain_old	Old state of total elastic strain
compute_full_tangent_operator	Flag currently unused by this class
tangent_operator	Currently a copy of the elasticity tensor in this class

Reimplemented in RadialReturnBackstressStressUpdateBaseTempl< is_ad >.

Definition at line 237 of file RadialReturnStressUpdate.C.

Referenced by RadialReturnBackstressStressUpdateBaseTempl< is_ad >::updateState().

{

  // compute the deviatoric trial stress and trial strain from the current intermediate
  // configuration
  GenericRankTwoTensor<is_ad> deviatoric_trial_stress = stress_new.deviatoric();

  // compute the effective trial stress
  GenericReal<is_ad> dev_trial_stress_squared =
      deviatoric_trial_stress.doubleContraction(deviatoric_trial_stress);
  GenericReal<is_ad> effective_trial_stress = MetaPhysicL::raw_value(dev_trial_stress_squared)
                                                  ? std::sqrt(3.0 / 2.0 * dev_trial_stress_squared)
                                                  : 0.0;

  computeStressInitialize(effective_trial_stress, elasticity_tensor);

  mooseAssert(
      _three_shear_modulus != 0.0,
      "Shear modulus is zero. Ensure that the base class computeStressInitialize() is called.");

  // Use Newton iteration to determine the scalar effective inelastic strain increment
  _effective_inelastic_strain_increment = 0.0;
  if (!MooseUtils::absoluteFuzzyEqual(effective_trial_stress, 0.0))
  {
    this->returnMappingSolve(
        effective_trial_stress, _effective_inelastic_strain_increment, this->_console);
    if (_effective_inelastic_strain_increment != 0.0)
      inelastic_strain_increment =
          deviatoric_trial_stress *
          (1.5 * _effective_inelastic_strain_increment / effective_trial_stress);
    else
      inelastic_strain_increment.zero();
  }
  else
    inelastic_strain_increment.zero();

  if (_apply_strain)
  {
    strain_increment -= inelastic_strain_increment;
    updateEffectiveInelasticStrain(_effective_inelastic_strain_increment);

    // Use the old elastic strain here because we require tensors used by this class
    // to be isotropic and this method natively allows for changing in time
    // elasticity tensors
    stress_new = elasticity_tensor * (strain_increment + elastic_strain_old);
  }

  computeStressFinalize(inelastic_strain_increment);

  if constexpr (!is_ad)
  {
    if (compute_full_tangent_operator)
      computeTangentOperator(effective_trial_stress, stress_new, tangent_operator);
  }
  else
  {
    libmesh_ignore(compute_full_tangent_operator);
    libmesh_ignore(tangent_operator);
  }
}

updateState() [2/2]

void StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::updateState ( GR2 &  strain_increment, GR2 &  inelastic_strain_increment, const GR2 &  rotation_increment, GR2 &  stress_new, const RankTwoTensor &  stress_old, const GR4 &  elasticity_tensor, const RankTwoTensor &  elastic_strain_old, bool  compute_full_tangent_operator = false, RankFourTensor &  tangent_operator = StressUpdateBaseTempl<is_ad>::_identityTensor  )

virtualinherited

Given a strain increment that results in a trial stress, perform some procedure (such as an iterative return-mapping process) to produce an admissible stress, an elastic strain increment and an inelastic strain increment.

If _fe_problem.currentlyComputingJacobian() = true, then updateState also computes d(stress)/d(strain) (or some approximation to it).

This method is called by ComputeMultipleInelasticStress. This method is pure virutal: all inheriting classes must overwrite this method.

Parameters

strain_increment	Upon input: the strain increment. Upon output: the elastic strain increment
inelastic_strain_increment	The inelastic_strain resulting from the interative procedure
rotation_increment	The finite-strain rotation increment
stress_new	Upon input: the trial stress that results from applying strain_increment as an elastic strain. Upon output: the admissible stress
stress_old	The old value of stress
elasticity_tensor	The elasticity tensor
compute_full_tangent_operator	The calling routine would like the full consistent tangent operator to be placed in tangent_operator, if possible. This is irrelevant if _fe_problem.currentlyComputingJacobian() = false
tangent_operator	d(stress)/d(strain), or some approximation to it If compute_full_tangent_operator=false, then tangent_operator=elasticity_tensor is an appropriate choice. tangent_operator is only computed if _fe_problem.currentlyComputingJacobian() = true

Definition at line 68 of file StressUpdateBase.C.

{
  mooseError("updateState called: it needs to be implemented by your inelastic model");
}

updateStateSubstep() [1/2]

template<bool is_ad>

void RadialReturnStressUpdateTempl< is_ad >::updateStateSubstep ( GenericRankTwoTensor< is_ad > &  strain_increment, GenericRankTwoTensor< is_ad > &  inelastic_strain_increment, const GenericRankTwoTensor< is_ad > &  rotation_increment, GenericRankTwoTensor< is_ad > &  stress_new, const RankTwoTensor &  stress_old, const GenericRankFourTensor< is_ad > &  elasticity_tensor, const RankTwoTensor &  elastic_strain_old, bool  compute_full_tangent_operator = false, RankFourTensor &  tangent_operator = StressUpdateBaseTempl<is_ad>::_identityTensor  )

overridevirtualinherited

Similar to the updateState function, this method updates the strain and stress for one substep.

Definition at line 412 of file RadialReturnStressUpdate.C.

{
  unsigned int num_substeps = calculateNumberSubsteps(strain_increment);
  _dt_original = _dt;
  while (true)
  {
    try
    {
      updateStateSubstepInternal(strain_increment,
                                 inelastic_strain_increment,
                                 rotation_increment,
                                 stress_new,
                                 stress_old,
                                 elasticity_tensor,
                                 elastic_strain_old,
                                 num_substeps,
                                 compute_full_tangent_operator,
                                 tangent_operator);
    }
    catch (MooseException & e)
    {
      // if we are not using adaptive substepping we just rethrow the exception
      if (!_adaptive_substepping)
        throw e;

      // otherwise we double the number of substeps and try again
      num_substeps *= 2;
      if (num_substeps <= _maximum_number_substeps)
        continue;

      // too meany substeps, break out of the loop
      break;
    }

    // updateStateSubstepInternal was successful (didn't throw)
    _dt = _dt_original;
    return;
  }

  // recover the original timestep
  _dt = _dt_original;

  mooseException("Adaptive substepping failed. Maximum number of substeps exceeded.");
}

updateStateSubstep() [2/2]

void StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::updateStateSubstep ( GR2 &  , GR2 &  , const GR2 &  , GR2 &  , const RankTwoTensor &  , const GR4 &  , const RankTwoTensor &  , bool  compute_full_tangent_operator = false, RankFourTensor &  tangent_operator = StressUpdateBaseTempl<is_ad>::_identityTensor  )

virtualinherited

Similar to the updateState function, this method updates the strain and stress for one substep.

Definition at line 84 of file StressUpdateBase.C.

{
  this->template paramError(
      "use_substep",
      "updateStateSubstep called: it needs to be implemented by your inelastic model");
}

updateStateSubstepInternal()

template<bool is_ad>

void RadialReturnStressUpdateTempl< is_ad >::updateStateSubstepInternal ( GenericRankTwoTensor< is_ad > &  strain_increment, GenericRankTwoTensor< is_ad > &  inelastic_strain_increment, const GenericRankTwoTensor< is_ad > &  rotation_increment, GenericRankTwoTensor< is_ad > &  stress_new, const RankTwoTensor &  stress_old, const GenericRankFourTensor< is_ad > &  elasticity_tensor, const RankTwoTensor &  elastic_strain_old, unsigned int  total_number_substeps, bool  compute_full_tangent_operator = false, RankFourTensor &  tangent_operator = StressUpdateBaseTempl<is_ad>::_identityTensor  )

virtualinherited

Definition at line 309 of file RadialReturnStressUpdate.C.

{
  // if only one substep is needed, then call the original update state method
  if (total_number_substeps == 1)
  {
    updateState(strain_increment,
                inelastic_strain_increment,
                rotation_increment,
                stress_new,
                stress_old,
                elasticity_tensor,
                elastic_strain_old,
                compute_full_tangent_operator,
                tangent_operator);

    this->storeIncrementalMaterialProperties(total_number_substeps);
    return;
  }

  if (total_number_substeps > _maximum_number_substeps)
    mooseException("The number of substeps computed exceeds 'maximum_number_substeps'.");

  // cut the original timestep
  _dt = _dt_original / total_number_substeps;

  // initialize the inputs
  const GenericRankTwoTensor<is_ad> strain_increment_per_step =
      strain_increment / total_number_substeps;
  GenericRankTwoTensor<is_ad> sub_stress_new = elasticity_tensor * elastic_strain_old;
  GenericRankTwoTensor<is_ad> sub_elastic_strain_old = elastic_strain_old;

  // clear the original inputs
  MathUtils::mooseSetToZero(strain_increment);
  MathUtils::mooseSetToZero(inelastic_strain_increment);
  MathUtils::mooseSetToZero(stress_new);

  GenericReal<is_ad> sub_effective_inelastic_strain_increment = 0.0;
  GenericRankTwoTensor<is_ad> sub_inelastic_strain_increment = inelastic_strain_increment;

  for (unsigned int step = 0; step < total_number_substeps; ++step)
  {
    // set up input for this substep
    GenericRankTwoTensor<is_ad> sub_strain_increment = strain_increment_per_step;
    sub_stress_new += elasticity_tensor * sub_strain_increment;

    Real effective_sub_stress_new;
    if constexpr (!is_ad)
    {
      // compute effective_sub_stress_new
      const RankTwoTensor deviatoric_sub_stress_new = sub_stress_new.deviatoric();
      const Real dev_sub_stress_new_squared =
          deviatoric_sub_stress_new.doubleContraction(deviatoric_sub_stress_new);
      effective_sub_stress_new = std::sqrt(3.0 / 2.0 * dev_sub_stress_new_squared);
    }
    else
      libmesh_ignore(effective_sub_stress_new);

    // update stress and strain based on the strain increment
    updateState(sub_strain_increment,
                sub_inelastic_strain_increment,
                rotation_increment, // not used in updateState
                sub_stress_new,
                stress_old, // not used in updateState
                elasticity_tensor,
                elastic_strain_old,
                false);
    // do not compute tangent until the end of this substep (or not at all for is_ad == true)

    // update strain and stress
    strain_increment += sub_strain_increment;
    inelastic_strain_increment += sub_inelastic_strain_increment;
    sub_elastic_strain_old += sub_strain_increment;
    sub_stress_new = elasticity_tensor * sub_elastic_strain_old;

    // accumulate scalar_effective_inelastic_strain
    sub_effective_inelastic_strain_increment += _effective_inelastic_strain_increment;

    if constexpr (!is_ad)
      computeTangentOperator(effective_sub_stress_new, sub_stress_new, tangent_operator);

    // store incremental material properties for this step
    this->storeIncrementalMaterialProperties(total_number_substeps);
  }

  // update stress
  stress_new = sub_stress_new;

  // update effective inelastic strain
  updateEffectiveInelasticStrain(sub_effective_inelastic_strain_increment);
}

validParams()

InputParameters ADMultiplePowerLawCreepStressUpdate::validParams ( )

static

Definition at line 16 of file ADMultiplePowerLawCreepStressUpdate.C.

{
  InputParameters params = ADRadialReturnCreepStressUpdateBase::validParams();
  params.addClassDescription(
      "This class uses the stress update material in a radial return isotropic power law creep "
      "model. This class can be used in conjunction with other creep and plasticity materials "
      "for more complex simulations.");

  // Linear strain hardening parameters
  params.addCoupledVar("temperature", "Coupled temperature");
  params.addRequiredParam<std::vector<Real>>("coefficient",
                                             "Leading coefficient in power-law equation");
  params.addRequiredParam<std::vector<Real>>("n_exponent",
                                             "Exponent on effective stress in power-law equation");
  params.addRequiredParam<std::vector<Real>>("m_exponent",
                                             "Exponent on time in power-law equation");
  params.addRequiredParam<std::vector<Real>>("stress_thresholds",
                                             "Stress intervals to switch creep behavior");
  params.addRequiredParam<std::vector<Real>>("activation_energy", "Activation energy");
  params.addParam<Real>("gas_constant", 8.3143, "Universal gas constant");
  params.addParam<Real>("start_time", 0.0, "Start time (if not zero)");
  return params;
}

Member Data Documentation

_activation_energy

const std::vector<Real> ADMultiplePowerLawCreepStressUpdate::_activation_energy

protected

Activation energy for exp term vector.

Definition at line 64 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by ADMultiplePowerLawCreepStressUpdate(), and computeStressInitialize().

_adaptive_substepping

template<bool is_ad>

const bool RadialReturnStressUpdateTempl< is_ad >::_adaptive_substepping

protectedinherited

Use adaptive substepping, cutting substep sizes until convergence is achieved.

Definition at line 260 of file RadialReturnStressUpdate.h.

Referenced by RadialReturnStressUpdateTempl< is_ad >::RadialReturnStressUpdateTempl().

_apply_strain

template<bool is_ad>

const bool RadialReturnStressUpdateTempl< is_ad >::_apply_strain

protectedinherited

Debugging option to enable specifying instead of calculating strain.

Definition at line 254 of file RadialReturnStressUpdate.h.

_base_name

const std::string StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::_base_name

protectedinherited

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

Definition at line 191 of file StressUpdateBase.h.

_bracket_solution

template<bool is_ad>

bool SingleVariableReturnMappingSolutionTempl< is_ad >::_bracket_solution

protectedinherited

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

Definition at line 145 of file SingleVariableReturnMappingSolution.h.

_check_range

template<bool is_ad>

bool SingleVariableReturnMappingSolutionTempl< is_ad >::_check_range

protectedinherited

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

Definition at line 138 of file SingleVariableReturnMappingSolution.h.

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

_coefficient

const std::vector<Real> ADMultiplePowerLawCreepStressUpdate::_coefficient

protected

Leading coefficient vector.

Definition at line 52 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by ADMultiplePowerLawCreepStressUpdate(), computeDerivative(), and computeResidual().

_creep_strain

template<bool is_ad>

GenericMaterialProperty<RankTwoTensor, is_ad>& RadialReturnCreepStressUpdateBaseTempl< is_ad >::_creep_strain

protectedinherited

Creep strain material property.

Definition at line 54 of file RadialReturnCreepStressUpdateBase.h.

_creep_strain_old

template<bool is_ad>

const MaterialProperty<RankTwoTensor>& RadialReturnCreepStressUpdateBaseTempl< is_ad >::_creep_strain_old

protectedinherited

Definition at line 55 of file RadialReturnCreepStressUpdateBase.h.

_current_elem

template<bool is_ad>

const Elem *const & Material::_current_elem

inherited

_deviatoric_projection_four

template<bool is_ad>

const RankFourTensor RadialReturnStressUpdateTempl< is_ad >::_deviatoric_projection_four

protectedinherited

Rank four deviatoric projection tensor.

Definition at line 251 of file RadialReturnStressUpdate.h.

_dt

template<bool is_ad>

Real & Material::_dt

inherited

_dt_original

template<bool is_ad>

Real RadialReturnStressUpdateTempl< is_ad >::_dt_original

protectedinherited

original timestep (to be restored after substepping is completed)

Definition at line 266 of file RadialReturnStressUpdate.h.

_effective_inelastic_strain

template<bool is_ad>

GenericMaterialProperty<Real, is_ad>& RadialReturnStressUpdateTempl< is_ad >::_effective_inelastic_strain

protectedinherited

Definition at line 219 of file RadialReturnStressUpdate.h.

Referenced by RadialReturnStressUpdateTempl< is_ad >::updateEffectiveInelasticStrain().

_effective_inelastic_strain_increment

template<bool is_ad>

GenericReal<is_ad> RadialReturnStressUpdateTempl< is_ad >::_effective_inelastic_strain_increment

protectedinherited

Stores the scalar effective inelastic strain increment from Newton iteration.

Definition at line 223 of file RadialReturnStressUpdate.h.

Referenced by RadialReturnStressUpdateTempl< is_ad >::effectiveInelasticStrainIncrement(), and RadialReturnStressUpdateTempl< is_ad >::updateEffectiveInelasticStrainIncrement().

_effective_inelastic_strain_old

template<bool is_ad>

const MaterialProperty<Real>& RadialReturnStressUpdateTempl< is_ad >::_effective_inelastic_strain_old

protectedinherited

Definition at line 220 of file RadialReturnStressUpdate.h.

Referenced by RadialReturnStressUpdateTempl< is_ad >::updateEffectiveInelasticStrain().

_exp_time

Real ADMultiplePowerLawCreepStressUpdate::_exp_time

protected

Exponential calculated from current time.

Definition at line 76 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by computeDerivative(), computeResidual(), and computeStressInitialize().

_exponential

ADReal ADMultiplePowerLawCreepStressUpdate::_exponential

protected

Exponential calculated from activation, gas constant, and temperature.

Definition at line 73 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by computeDerivative(), computeResidual(), and computeStressInitialize().

_gas_constant

const Real ADMultiplePowerLawCreepStressUpdate::_gas_constant

protected

Gas constant for exp term.

Definition at line 67 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by computeStressInitialize().

_identity_symmetric_four

template<bool is_ad>

const RankFourTensor RadialReturnStressUpdateTempl< is_ad >::_identity_symmetric_four

protectedinherited

Rank four symmetric identity tensor.

Definition at line 246 of file RadialReturnStressUpdate.h.

_identity_two

template<bool is_ad>

const RankTwoTensor RadialReturnStressUpdateTempl< is_ad >::_identity_two

protectedinherited

Rank two identity tensor.

Definition at line 241 of file RadialReturnStressUpdate.h.

_identityTensor

RankFourTensor StressUpdateBaseTempl< is_ad, RankTwoTensor , RankFourTensor >::_identityTensor

staticprotectedinherited

Definition at line 193 of file StressUpdateBase.h.

_line_search

template<bool is_ad>

bool SingleVariableReturnMappingSolutionTempl< is_ad >::_line_search

protectedinherited

Whether to use line searches to improve convergence.

Definition at line 141 of file SingleVariableReturnMappingSolution.h.

_m_exponent

const std::vector<Real> ADMultiplePowerLawCreepStressUpdate::_m_exponent

protected

Exponent on time vector.

Definition at line 58 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by ADMultiplePowerLawCreepStressUpdate(), and computeStressInitialize().

_max_inelastic_increment

template<bool is_ad>

Real RadialReturnStressUpdateTempl< is_ad >::_max_inelastic_increment

protectedinherited

Maximum allowable scalar inelastic strain increment, used to control the timestep size in conjunction with a user object.

Definition at line 229 of file RadialReturnStressUpdate.h.

_maximum_number_substeps

template<bool is_ad>

const unsigned int RadialReturnStressUpdateTempl< is_ad >::_maximum_number_substeps

protectedinherited

Maximum number of substeps. If the calculation results in a larger number, cut overall time step.

Definition at line 263 of file RadialReturnStressUpdate.h.

_n_exponent

const std::vector<Real> ADMultiplePowerLawCreepStressUpdate::_n_exponent

protected

Exponent on the effective stress vector.

Definition at line 55 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by ADMultiplePowerLawCreepStressUpdate(), computeDerivative(), computeResidual(), and computeStrainEnergyRateDensity().

_number_of_models

const unsigned int ADMultiplePowerLawCreepStressUpdate::_number_of_models

protected

Total number of models provided by the user.

Definition at line 82 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by ADMultiplePowerLawCreepStressUpdate().

_q_point

template<bool is_ad>

const MooseArray< Point > & Material::_q_point

inherited

_qp

template<bool is_ad>

unsigned int Material::_qp

inherited

_start_time

const Real ADMultiplePowerLawCreepStressUpdate::_start_time

protected

Simulation start time.

Definition at line 70 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by computeStressInitialize().

_stress_index

unsigned int ADMultiplePowerLawCreepStressUpdate::_stress_index

protected

Quadrature-point value pointing to the right power law parameters.

Definition at line 79 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by computeDerivative(), computeResidual(), computeStrainEnergyRateDensity(), and computeStressInitialize().

_stress_thresholds

const std::vector<Real> ADMultiplePowerLawCreepStressUpdate::_stress_thresholds

protected

Stress thresholds vector.

Definition at line 61 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by ADMultiplePowerLawCreepStressUpdate(), and stressIndex().

_substep_tolerance

template<bool is_ad>

const Real RadialReturnStressUpdateTempl< is_ad >::_substep_tolerance

protectedinherited

Used to calculate the number of substeps taken in the radial return algorithm, when substepping is enabled, based on the elastic strain increment ratio to the maximum inelastic increment.

Definition at line 236 of file RadialReturnStressUpdate.h.

_temperature

const ADVariableValue* const ADMultiplePowerLawCreepStressUpdate::_temperature

protected

Temperature variable value.

Definition at line 49 of file ADMultiplePowerLawCreepStressUpdate.h.

Referenced by computeStressInitialize().

_three_shear_modulus

template<bool is_ad>

GenericReal<is_ad> RadialReturnStressUpdateTempl< is_ad >::_three_shear_modulus

protectedinherited

3 * shear modulus

Definition at line 217 of file RadialReturnStressUpdate.h.

Referenced by computeDerivative(), and computeResidual().

_use_substepping

template<bool is_ad>

SubsteppingType RadialReturnStressUpdateTempl< is_ad >::_use_substepping

protectedinherited

Whether user has requested the use of substepping technique to improve convergence [make const later].

Definition at line 257 of file RadialReturnStressUpdate.h.

Referenced by RadialReturnStressUpdateTempl< is_ad >::RadialReturnStressUpdateTempl(), substeppingCapabilityEnabled(), and RadialReturnStressUpdateTempl< is_ad >::substeppingCapabilityRequested().

---

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

• solid_mechanics/include/materials/ADMultiplePowerLawCreepStressUpdate.h
• solid_mechanics/src/materials/ADMultiplePowerLawCreepStressUpdate.C