#include <LAROMANCEStressUpdateBase.h>

Inheritance diagram for LAROMANCEStressUpdateBaseTempl< is_ad >:

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
	LAROMANCEStressUpdateBaseTempl (const InputParameters &parameters)

virtual void	resetIncrementalMaterialProperties () override
	Reset material properties. More...

virtual void	storeIncrementalMaterialProperties (const unsigned int total_number_substeps) override
	Properly set up the incremental calculation storage of the stateful material properties in the inheriting classes. 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...

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 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 *	queryParam (const std::string &name) 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 (const std::set< 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	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)

void	addMooseVariableDependency (MooseVariableFieldBase *var)

void	addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)

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
const Elem *const &	_current_elem

Real &	_dt

	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	ROMInputTransform { ROMInputTransform::LINEAR, ROMInputTransform::LOG, ROMInputTransform::EXP }

enum	WindowFailure { WindowFailure::ERROR, WindowFailure::EXCEPTION, WindowFailure::WARN, WindowFailure::IGNORE, WindowFailure::DONOTHING, WindowFailure::USELIMIT, WindowFailure::EXTRAPOLATE }
	Enum to error, warn, ignore, or extrapolate if input is outside of window of applicability. More...

enum	QP_Data_Type

Protected Member Functions
virtual void	exportJSON ()

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	initialSetup () override

virtual void	setupUnitConversionFactors (const InputParameters &parameters)

virtual void	initQpStatefulProperties () override

virtual void	computeStressInitialize (const GenericReal< is_ad > &effective_trial_stress, const GenericRankFourTensor< is_ad > &elasticity_tensor) override
	Perform any necessary initialization before return mapping iterations. More...

virtual GenericReal< is_ad >	computeResidual (const GenericReal< is_ad > &effective_trial_stress, const GenericReal< is_ad > &scalar) override
	Compute the residual for a predicted value of the scalar. More...

virtual GenericReal< is_ad >	computeDerivative (const GenericReal< is_ad > &, const GenericReal< is_ad > &) override
	Compute the derivative of the residual as a function of the scalar variable. 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 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...

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	outputIterationStep (std::stringstream *iter_output, const GenericReal< is_ad > &effective_trial_stress, const GenericReal< is_ad > &scalar, const Real reference_residual) override
	Output information for a single iteration step to build the convergence history of the model. More...

void	precomputeROM (const unsigned out_index)
	Precompute the ROM strain rate information for all inputs except for strain. More...

GenericReal< is_ad >	computeROM (const unsigned int tile, const unsigned int partition, const unsigned out_index, const bool derivative=false)
	Computes the ROM calculated increment for a given output and tile. More...

bool	checkInTile (const unsigned int p, const unsigned int t) const
	Checks if the input combination is in a specific tile. More...

bool	areTilesNotIdentical (const unsigned int p, const unsigned int t, const unsigned int tt, const unsigned int in_index)
	Checks if two tile domains are equal. More...

GenericReal< is_ad >	normalizeInput (const GenericReal< is_ad > &input, const ROMInputTransform transform, const Real transform_coef, const std::vector< Real > &transformed_limits, const bool derivative=false)
	Convert the input variables into the form expected by the ROM Legendre polynomials to have a normalized space from [-1, 1] so that every variable has equal weight. More...

void	buildPolynomials (const unsigned int p, const GenericReal< is_ad > &rom_input, std::vector< GenericReal< is_ad >> &polynomial_inputs, const GenericReal< is_ad > &drom_input=0, const bool derivative=false)
	Assemble the array of Legendre polynomials to be multiplied by the ROM coefficients. More...

void	precomputeValues (const unsigned int p, const std::vector< Real > &coefs, const std::vector< std::vector< GenericReal< is_ad >>> &polynomial_inputs, std::vector< GenericReal< is_ad >> &precomputed)
	Arranges the calculated Legendre polynomials into the proper oder and multiplies the Legendre polynomials by the ROM coefficients to compute the predicted output values. More...

GenericReal< is_ad >	computeValues (const unsigned int p, const std::vector< GenericReal< is_ad >> &precomputed, const std::vector< std::vector< GenericReal< is_ad >>> &polynomial_inputs, const std::vector< GenericReal< is_ad >> &dpolynomial_inputs={}, const bool derivative=false)
	Arranges the calculated Legendre polynomials into the proper oder and multiplies the Legendre polynomials by the ROM coefficients to compute the predicted output values. More...

virtual GenericReal< is_ad >	convertOutput (const unsigned int p, const std::vector< Real > &old_input_values, const GenericReal< is_ad > &rom_output, const unsigned out_index, const GenericReal< is_ad > &drom_output=0.0, const bool derivative=false)
	Computes the output variable increments from the ROM predictions by bringing out of the normalized map to the actual physical values. More...

GenericReal< is_ad >	computePolynomial (const GenericReal< is_ad > &value, const unsigned int degree, const bool derivative=false)
	Calculate the value or derivative of Legendre polynomial up to 3rd order. More...

GenericReal< is_ad >	sigmoid (const Real lower, const Real upper, const GenericReal< is_ad > &val, const bool derivative=false)
	Calculate the sigmoid function weighting for the input based on the limits. More...

void	computeTileWeight (std::vector< std::vector< GenericReal< is_ad >>> &weights, GenericReal< is_ad > &input, const unsigned int in_index, const bool derivative=false)
	Compute the contribution (weight) of each tile in each partition, based on the input and tile boundaries (in terms of input domain). More...

virtual void	computePartitionWeights (std::vector< GenericReal< is_ad >> &weights, std::vector< GenericReal< is_ad >> &dweights_dstress)
	Compute the weight of the different partitions. More...

template<typename T >
void	convertValue (T &x, const ROMInputTransform transform, const Real coef, const bool derivative=false)
	Convert input based on the transform type. More...

std::vector< unsigned int >	getMakeFrameHelper (const unsigned int p) const

std::vector< std::vector< std::vector< std::vector< Real > > > >	getTransformedLimits (const unsigned int p, const std::vector< std::vector< std::vector< Real >>> limits)

virtual std::vector< std::vector< std::vector< std::vector< ROMInputTransform > > > >	getTransform ()

virtual std::vector< std::vector< std::vector< std::vector< Real > > > >	getTransformCoefs ()

virtual std::vector< std::vector< std::vector< std::vector< Real > > > >	getNormalizationLimits ()

virtual std::vector< std::vector< std::vector< std::vector< Real > > > >	getInputLimits ()

virtual std::vector< std::vector< std::vector< std::vector< Real > > > >	getCoefs ()

virtual std::vector< std::vector< unsigned int > >	getTilings ()

virtual std::vector< Real >	getStrainCutoff ()

void	checkJSONKey (const std::string &key)
	check if a JSON file was loaded and if the specified key exists More...

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

virtual void	propagateQpStatefulProperties () override
	If updateState is not called during a timestep, this will be. 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	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

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 GenericVectorVariableValue< is_ad > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp=0) const

const GenericVectorVariableValue< false > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const

const GenericVectorVariableValue< true > &	coupledGenericVectorValue (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

std::vector< const VectorVariableValue *>	coupledVectorValuesOld (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 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...

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 GenericVariableValue< is_ad > &	_temperature
	Coupled temperature variable. More...

const GenericMaterialProperty< Real, is_ad > *	_environmental
	Optionally coupled environmental factor. More...

std::vector< std::pair< WindowFailure, WindowFailure > >	_window_failure

const bool	_verbose
	Flag to output verbose infromation. More...

const Real	_max_cell_increment
	Maximum cell dislocation increment. More...

const Function *const	_cell_function
	Optional cell dislocation forcing function. More...

GenericReal< is_ad >	_cell_dislocation_increment
	Container for cell dislocation increment. More...

const Real	_max_wall_increment
	Maximum wall dislocation increment. More...

const Function *const	_wall_function
	Optional wall dislocation forcing function. More...

const Function *const	_stress_function
	Optiontal effective stress forcing function. More...

GenericReal< is_ad >	_wall_dislocation_increment
	Container for wall dislocation increment. More...

const unsigned int	_cell_input_index
	Index corresponding to the position for the dislocations with in the cell in the input vector. More...

const unsigned int	_wall_input_index
	Index corresponding to the position for the dislocations within the cell wall in the input vector. More...

const unsigned int	_stress_input_index
	Index corresponding to the position for the stress in the input vector. More...

const unsigned int	_old_strain_input_index
	Index corresponding to the position for the old strain in the input vector. More...

const unsigned int	_temperature_input_index
	Index corresponding to the position for the tempeature in the input vector. More...

const unsigned int	_environmental_input_index
	Index corresponding to the position for the environmental factor in the input vector. More...

const unsigned int	_cell_output_index
	Index corresponding to the position for cell dislocations increment in the output vector. More...

const unsigned int	_wall_output_index
	Index corresponding to the position for wall dislocations increment in the output vector. More...

const unsigned int	_strain_output_index
	Index corresponding to the position for strain increment in the output vector. More...

const Function *const	_creep_strain_old_forcing_function
	Optional old creep strain forcing function. More...

unsigned int	_num_partitions
	Number of partitions. More...

std::vector< unsigned int >	_num_tiles
	Number of ROM tiles per partition. More...

unsigned int	_num_inputs
	Number of inputs for the ROM data set. More...

unsigned int	_num_outputs
	Number of inputs to the ROM data set. More...

std::vector< unsigned int >	_degree
	Legendre polynomial degree for the ROM data set for each partition. More...

std::vector< unsigned int >	_num_coefs
	Total number of Legendre polynomial coefficients for the ROM data set in each parition. More...

std::vector< std::vector< std::vector< std::vector< ROMInputTransform > > > >	_transform
	Transform rules defined by the ROM data set for each partition. More...

std::vector< std::vector< std::vector< std::vector< Real > > > >	_transform_coefs
	Transform coefficients defined by the ROM data set for each partition. More...

std::vector< std::vector< std::vector< std::vector< Real > > > >	_input_limits
	Input limits defined by the ROM data set for each partition. More...

std::vector< std::vector< std::vector< std::vector< Real > > > >	_normalization_limits
	Normalization limits defined by the ROM data set for each partition. More...

std::vector< std::vector< std::vector< std::vector< Real > > > >	_coefs
	Coefficients used with Legendre polynomials defined by the ROM data set for each partition. More...

std::vector< std::vector< std::vector< std::vector< std::vector< Real > > > > >	_transformed_normalization_limits
	Limits transformed from readabile input to ROM readable limits for normalization. More...

std::vector< std::vector< unsigned int > >	_makeframe_helper
	Helper container defined by the ROM data set. More...

GenericMaterialProperty< Real, is_ad > &	_creep_rate
	Creep rate material property. More...

GenericMaterialProperty< Real, is_ad > &	_cell_rate
	Cell dislocations rate of change. More...

GenericMaterialProperty< Real, is_ad > &	_wall_rate
	Wall dislocations rate of change. More...

MaterialProperty< Real > &	_extrapolation
	Material property to hold smootherstep applied in order to extrapolate. More...

GenericMaterialProperty< Real, is_ad > &	_second_partition_weight
	Material property to store partition weight. More...

GenericReal< is_ad >	_derivative
	Container for derivative of creep increment with respect to strain. More...

std::vector< GenericReal< is_ad > >	_input_values
	Container for input values. More...

std::vector< Real >	_old_input_values
	Container for old input values. More...

std::vector< std::vector< std::vector< GenericReal< is_ad > > > >	_rom_inputs
	Container for converted rom_inputs. More...

std::vector< std::vector< std::vector< std::vector< GenericReal< is_ad > > > > >	_polynomial_inputs
	Container for ROM polynomial inputs. More...

std::vector< std::vector< std::vector< GenericReal< is_ad > > > >	_precomputed_vals
	Container for ROM precomputed values. More...

std::vector< std::pair< Real, Real > >	_global_limits
	Container for global limits. More...

std::vector< std::vector< GenericReal< is_ad > > >	_non_stress_weights
	Container for weights for each tile as computed for all input values beside stress. More...

std::vector< std::vector< GenericReal< is_ad > > >	_weights
	Container for weights for each tile as computed for all input values beside stress. More...

std::vector< GenericReal< is_ad > >	_partition_weights
	Container for weights for each tile as computed for all input values beside stress. More...

std::vector< GenericReal< is_ad > >	_dpartition_weight_dstress
	Container for d_parition_weights d_stress. More...

std::vector< std::vector< unsigned int > >	_tiling
	Container for tiling orientations. More...

std::vector< Real >	_cutoff
	Container for strain cutoff. More...

Real	_stress_ucf
	Unit conversion factors required to convert from the specified unit to MPa. More...

RankTwoTensor	_plastic_strain_increment
	Total plastic strain increment in step (summing substep contributions) More...

MaterialProperty< Real > &	_number_of_substeps
	Material property capturing number of substeps for output purposes (defaults to one if substepping isn't used) More...

std::vector< std::string >	_index_name
	index names for error output More...

nlohmann::json	_json
	JSON object constructed from the datafile. More...

	usingTransientInterfaceMembers

const std::string	_name

const MooseArray< Point > &	_q_point

unsigned int	_qp

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


GenericMaterialProperty< Real, is_ad > &	_cell_dislocations
	Material properties for cell (glissile) dislocation densities (1/m^2) More...

const MaterialProperty< Real > &	_cell_dislocations_old


GenericMaterialProperty< Real, is_ad > &	_wall_dislocations
	Material properties for wall (locked) dislocation densities (1/m^2) More...

const MaterialProperty< Real > &	_wall_dislocations_old


GenericMaterialProperty< Real, is_ad > &	_wall_dislocations_step
	Material properties accumulated at substeps. More...

GenericMaterialProperty< Real, is_ad > &	_cell_dislocations_step

Static Protected Attributes
static RankFourTensor	_identityTensor

static const std::string	_interpolated_old

static const std::string	_interpolated_older

Detailed Description

template<bool is_ad>
class LAROMANCEStressUpdateBaseTempl< is_ad >

Definition at line 16 of file LAROMANCEStressUpdateBase.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

◆ ROMInputTransform

template<bool is_ad>

enum LAROMANCEStressUpdateBaseTempl::ROMInputTransform

strongprotected

Enumerator
LINEAR
LOG
EXP

Definition at line 33 of file LAROMANCEStressUpdateBase.h.

   {
     LINEAR,
     LOG,
     EXP
   };

◆ 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
   };

◆ WindowFailure

template<bool is_ad>

enum LAROMANCEStressUpdateBaseTempl::WindowFailure

strongprotected

Enum to error, warn, ignore, or extrapolate if input is outside of window of applicability.

Enumerator
ERROR
EXCEPTION
WARN
IGNORE
DONOTHING
USELIMIT
EXTRAPOLATE

Definition at line 80 of file LAROMANCEStressUpdateBase.h.

   {
     ERROR,
     EXCEPTION,
     WARN,
     IGNORE,
     DONOTHING,
     USELIMIT,
     EXTRAPOLATE
   };

Constructor & Destructor Documentation

◆ LAROMANCEStressUpdateBaseTempl()

template<bool is_ad>

LAROMANCEStressUpdateBaseTempl< is_ad >::LAROMANCEStressUpdateBaseTempl ( const InputParameters & parameters )

Definition at line 166 of file LAROMANCEStressUpdateBase.C.

   : RadialReturnCreepStressUpdateBaseTempl<is_ad>(parameters),
     _temperature(this->template coupledGenericValue<is_ad>("temperature")),
     _environmental(
         this->isParamValid("environmental_factor")
             ? &this->template getGenericMaterialProperty<Real, is_ad>("environmental_factor")
             : nullptr),
     _window_failure(_environmental ? 6 : 5),
 
     _verbose(this->template getParam<bool>("verbose")),
     _cell_dislocations(
         this->template declareGenericProperty<Real, is_ad>(this->_base_name + "cell_dislocations")),
     _cell_dislocations_old(
         this->template getMaterialPropertyOld<Real>(this->_base_name + "cell_dislocations")),
     _max_cell_increment(this->template getParam<Real>("max_relative_cell_dislocation_increment")),
     _cell_function(this->isParamValid("cell_dislocation_density_forcing_function")
                        ? &this->getFunction("cell_dislocation_density_forcing_function")
                        : NULL),
     _cell_dislocation_increment(0.0),
     _wall_dislocations(
         this->template declareGenericProperty<Real, is_ad>(this->_base_name + "wall_dislocations")),
     _wall_dislocations_old(
         this->template getMaterialPropertyOld<Real>(this->_base_name + "wall_dislocations")),
     _max_wall_increment(this->template getParam<Real>("max_relative_wall_dislocation_increment")),
     _wall_function(this->isParamValid("wall_dislocation_density_forcing_function")
                        ? &this->getFunction("wall_dislocation_density_forcing_function")
                        : NULL),
     _stress_function(this->isParamValid("effective_stress_forcing_function")
                          ? &this->getFunction("effective_stress_forcing_function")
                          : NULL),
     _wall_dislocation_increment(0.0),
     _cell_input_index(0),
     _wall_input_index(1),
     _stress_input_index(2),
     _old_strain_input_index(3),
     _temperature_input_index(4),
     _environmental_input_index(5),
     _cell_output_index(0),
     _wall_output_index(1),
     _strain_output_index(2),
 
     _creep_strain_old_forcing_function(this->isParamValid("old_creep_strain_forcing_function")
                                            ? &this->getFunction("old_creep_strain_forcing_function")
                                            : NULL),
 
     _creep_rate(
         this->template declareGenericProperty<Real, is_ad>(this->_base_name + "creep_rate")),
     _cell_rate(this->template declareGenericProperty<Real, is_ad>(this->_base_name +
                                                                   "cell_dislocation_rate")),
     _wall_rate(this->template declareGenericProperty<Real, is_ad>(this->_base_name +
                                                                   "wall_dislocation_rate")),
     _extrapolation(this->template declareProperty<Real>("ROM_extrapolation")),
     _second_partition_weight(
         this->template declareGenericProperty<Real, is_ad>("partition_weight")),
 
     _derivative(0.0),
     _old_input_values(3),
     _wall_dislocations_step(this->template declareGenericProperty<Real, is_ad>(
         this->_base_name + "wall_dislocations_step")),
     _cell_dislocations_step(this->template declareGenericProperty<Real, is_ad>(
         this->_base_name + "cell_dislocations_step")),
     _plastic_strain_increment(),
     _number_of_substeps(
         this->template declareProperty<Real>(this->_base_name + "number_of_substeps")),
     _index_name(_window_failure.size())
 {
   this->_check_range = true; // this may not be necessary?
 
   _index_name[_cell_input_index] = "cell";
   _index_name[_wall_input_index] = "wall";
   _index_name[_stress_input_index] = "stress";
   _index_name[_old_strain_input_index] = "old strain";
   _index_name[_temperature_input_index] = "temperature";
 
   _window_failure[_cell_input_index].first =
       parameters.get<MooseEnum>("cell_input_window_low_failure").getEnum<WindowFailure>();
   _window_failure[_cell_input_index].second =
       parameters.get<MooseEnum>("cell_input_window_high_failure").getEnum<WindowFailure>();
   _window_failure[_wall_input_index].first =
       parameters.get<MooseEnum>("wall_input_window_low_failure").getEnum<WindowFailure>();
   _window_failure[_wall_input_index].second =
       parameters.get<MooseEnum>("wall_input_window_high_failure").getEnum<WindowFailure>();
   _window_failure[_stress_input_index].first =
       parameters.get<MooseEnum>("stress_input_window_low_failure").getEnum<WindowFailure>();
   _window_failure[_stress_input_index].second =
       parameters.get<MooseEnum>("stress_input_window_high_failure").getEnum<WindowFailure>();
   _window_failure[_old_strain_input_index].first =
       parameters.get<MooseEnum>("old_strain_input_window_low_failure").getEnum<WindowFailure>();
   _window_failure[_old_strain_input_index].second =
       parameters.get<MooseEnum>("old_strain_input_window_high_failure").getEnum<WindowFailure>();
   _window_failure[_temperature_input_index].first =
       parameters.get<MooseEnum>("temperature_input_window_low_failure").getEnum<WindowFailure>();
   _window_failure[_temperature_input_index].second =
       parameters.get<MooseEnum>("temperature_input_window_high_failure").getEnum<WindowFailure>();
   if (_environmental)
   {
     _index_name[_environmental_input_index] = "environmental";
     _window_failure[_environmental_input_index].first =
         parameters.get<MooseEnum>("environment_input_window_low_failure").getEnum<WindowFailure>();
     _window_failure[_environmental_input_index].second =
         parameters.get<MooseEnum>("environment_input_window_high_failure").getEnum<WindowFailure>();
   }
 
   // load JSON datafile
   if (this->isParamValid("model"))
   {
     const auto model_file_name = this->template getParam<DataFileName>("model");
     std::ifstream model_file(model_file_name.c_str());
     model_file >> _json;
   }
 
   setupUnitConversionFactors(parameters);
 }

Member Function Documentation

◆ areTilesNotIdentical()

template<bool is_ad>

bool LAROMANCEStressUpdateBaseTempl< is_ad >::areTilesNotIdentical	(	const unsigned int	p,
		const unsigned int	t,
		const unsigned int	tt,
		const unsigned int	in_index
	)

protected

Checks if two tile domains are equal.

Parameters

p	Partition index
t	Tile 1 index
tt	Tile 2 index
in_index	input index

Returns: integer 0 (false) or 1 (true)

Definition at line 761 of file LAROMANCEStressUpdateBase.C.

 {
   if (_input_limits[p][t][in_index][0] != _input_limits[p][tt][in_index][0] &&
       _input_limits[p][t][in_index][1] != _input_limits[p][tt][in_index][1])
     return true;
   else
     return false;
 }

◆ buildPolynomials()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::buildPolynomials	(	const unsigned int	p,
		const GenericReal< is_ad > &	rom_input,
		std::vector< GenericReal< is_ad >> &	polynomial_inputs,
		const GenericReal< is_ad > &	drom_input = `0`,
		const bool	derivative = `false`
	)

protected

Assemble the array of Legendre polynomials to be multiplied by the ROM coefficients.

Parameters

p	Partition index
rom_input	ROM input
polynomial_inputs	Vector of transformed Legendre polynomials
drom_input	Optional derivative of ROM input with respect to stress
derivative	Optional flag to return derivative of converted input with respect to stress.

Definition at line 1015 of file LAROMANCEStressUpdateBase.C.

 {
   for (unsigned int d = 0; d < _degree[p]; ++d)
   {
     polynomial_inputs[d] = computePolynomial(rom_input, d);
 
     if (derivative)
       polynomial_inputs[d] = drom_input * computePolynomial(rom_input, d, derivative);
   }
 }

◆ 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.

160 { return 1; }

◆ checkInTile()

template<bool is_ad>

bool LAROMANCEStressUpdateBaseTempl< is_ad >::checkInTile	(	const unsigned int	p,
		const unsigned int	t
	)		const

protected

Checks if the input combination is in a specific tile.

Parameters

p	Partition index
t	Tile index

Returns: bool if in tile

Definition at line 750 of file LAROMANCEStressUpdateBase.C.

 {
   for (unsigned int i = 0; i < _num_inputs; ++i)
     if (_input_values[i] < _input_limits[p][t][i][0] ||
         _input_values[i] > _input_limits[p][t][i][1])
       return false;
   return true;
 }

◆ checkJSONKey()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::checkJSONKey ( const std::string & key )

protected

check if a JSON file was loaded and if the specified key exists

Definition at line 1393 of file LAROMANCEStressUpdateBase.C.

Referenced by LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::getClassificationEll(), LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::getClassificationEta(), LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::getClassificationLuu(), LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::getClassificationMmean(), LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::getClassificationMscale(), LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::getClassificationVind(), LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::getClassificationXu(), LAROMANCEStressUpdateBaseTempl< is_ad >::getCoefs(), LAROMANCEStressUpdateBaseTempl< is_ad >::getInputLimits(), LAROMANCEStressUpdateBaseTempl< is_ad >::getStrainCutoff(), LAROMANCEStressUpdateBaseTempl< is_ad >::getTransform(), and LAROMANCEStressUpdateBaseTempl< is_ad >::getTransformCoefs().

 {
   if (!this->isParamValid("model"))
     this->paramError("model", "Specify a JSON data filename.");
 
   const auto model_file_name = this->template getParam<DataFileName>("model");
   if (_json.empty())
     this->paramError("model", "The specified JSON data file '", model_file_name, "' is empty.");
   if (!_json.contains(key))
     this->paramError(
         "model", "The key '", key, "' is missing from the JSON data file '", model_file_name, "'.");
 }

◆ computeDerivative()

template<bool is_ad>

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

inlineoverrideprotectedvirtual

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

Implements SingleVariableReturnMappingSolutionTempl< is_ad >.

Definition at line 59 of file LAROMANCEStressUpdateBase.h.

   {
     return _derivative;
   }

◆ computePartitionWeights()

template<bool is_ad>

virtual void LAROMANCEStressUpdateBaseTempl< is_ad >::computePartitionWeights	(	std::vector< GenericReal< is_ad >> &	weights,
		std::vector< GenericReal< is_ad >> &	dweights_dstress
	)

inlineprotectedvirtual

Compute the weight of the different partitions.

Parameters

weights	Weights for each partition
derivative	Optional flag to return derivative of the sigmoid w.r.t. the input

Reimplemented in LAROMANCEPartitionStressUpdateBaseTempl< is_ad >.

Definition at line 257 of file LAROMANCEStressUpdateBase.h.

   {
     if (_num_partitions != 1)
       mooseError("Internal error: If number of partitions is not one, then computePartitionWeights "
                  "must be defined");
     weights[0] = 1.0;
     dweights_dstress[0] = 0.0;
   }

◆ computePolynomial()

template<bool is_ad>

GenericReal< is_ad > LAROMANCEStressUpdateBaseTempl< is_ad >::computePolynomial	(	const GenericReal< is_ad > &	value,
		const unsigned int	degree,
		const bool	derivative = `false`
	)

protected

Calculate the value or derivative of Legendre polynomial up to 3rd order.

Parameters

value	Input to Legendre polynomial
degree	Degree of Legendre polynomial
derivative	Optional flag to return derivative of Legendre polynomial Legendre

Returns: Computed value from Legendre polynomial

Definition at line 1107 of file LAROMANCEStressUpdateBase.C.

 {
   if (degree == 0)
   {
     if (derivative)
       return 0.0;
     return 1.0;
   }
   else if (degree == 1)
   {
     if (derivative)
       return 1.0;
     return value;
   }
   else if (degree == 2)
   {
     if (derivative)
       return 3.0 * value;
     return 1.5 * Utility::pow<2>(value) - 0.5;
   }
   else
   {
     if (derivative)
       return 7.5 * Utility::pow<2>(value) - 1.5;
     return 2.5 * Utility::pow<3>(value) - 1.5 * value;
   }
 }

◆ 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()

template<bool is_ad>

GenericReal< is_ad > LAROMANCEStressUpdateBaseTempl< is_ad >::computeResidual	(	const GenericReal< is_ad > &	,
		const GenericReal< is_ad > &
	)

overrideprotectedvirtual

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

Implements SingleVariableReturnMappingSolutionTempl< is_ad >.

Definition at line 775 of file LAROMANCEStressUpdateBase.C.

 {
   mooseAssert(std::isfinite(MetaPhysicL::raw_value(effective_trial_stress)),
               "computeResidual: effective_trial_stress must be finite");
   mooseAssert(std::isfinite(MetaPhysicL::raw_value(scalar)),
               "computeResidual: scalar must be finite!");
   // Update new stress
   GenericReal<is_ad> trial_stress_mpa = _stress_function
                                             ? _stress_function->value(_t, _q_point[_qp])
                                             : effective_trial_stress * _stress_ucf;
   GenericReal<is_ad> dtrial_stress_dscalar = 0.0;
 
   // Update stress if strain is being applied, i.e. non-testing simulation
   if (this->_apply_strain)
   {
     trial_stress_mpa -= this->_three_shear_modulus * scalar * _stress_ucf;
     dtrial_stress_dscalar -= this->_three_shear_modulus * _stress_ucf;
   }
   _input_values[_stress_input_index] = trial_stress_mpa;
 
   // Update weights for each partition with new stress
   for (unsigned int p = 0; p < _num_partitions; ++p)
     _weights[p] = _non_stress_weights[p];
   computeTileWeight(_weights, _input_values[_stress_input_index], _stress_input_index);
   auto dweights_dstress = _non_stress_weights;
   computeTileWeight(
       dweights_dstress, _input_values[_stress_input_index], _stress_input_index, true);
 
   computePartitionWeights(_partition_weights, _dpartition_weight_dstress);
 
   // Save extrapolation as a material property in order quantify adequate tiling range
   _extrapolation[_qp] = 0.0;
   for (unsigned int p = 0; p < _num_partitions; ++p)
     for (unsigned int t = 0; t < _num_tiles[p]; ++t)
       _extrapolation[_qp] += MetaPhysicL::raw_value(_weights[p][t] * _partition_weights[p]);
 
   GenericReal<is_ad> total_rom_effective_strain_inc = 0.0;
   GenericReal<is_ad> dtotal_rom_effective_strain_inc_dstress = 0.0;
 
   // Run ROM if all values are within windows.
   for (unsigned int p = 0; p < _num_partitions; p++)
   {
     if (_partition_weights[p])
     {
       // compute weight normalizing factor
       GenericReal<is_ad> weight_normalizer = 0;
       unsigned int number_of_active_tiles = 0;
       for (unsigned int t = 0; t < _num_tiles[p]; ++t)
       {
         // count number of active tiles
         number_of_active_tiles += checkInTile(p, t);
         if (_weights[p][t])
         {
           // tile normalization factor (sum of tile weights)
           weight_normalizer += _weights[p][t];
         }
       }
 
       // normalize weights only when 3 tiles overlap
       if (number_of_active_tiles == 3)
         for (unsigned int t = 0; t < _num_tiles[p]; ++t)
         {
           _weights[p][t] /= weight_normalizer;
           dweights_dstress[p][t] /= weight_normalizer;
         }
 
       for (unsigned int t = 0; t < _num_tiles[p]; ++t)
         if (_weights[p][t])
         {
           const GenericReal<is_ad> rom = computeROM(t, p, _strain_output_index);
           if (rom == std::numeric_limits<float>::infinity())
             mooseError("In ", _name, ": Output for strain increment reaches infinity: ", rom);
 
           total_rom_effective_strain_inc += _partition_weights[p] * _weights[p][t] * rom;
 
           dtotal_rom_effective_strain_inc_dstress +=
               _partition_weights[p] * _weights[p][t] * computeROM(t, p, _strain_output_index, true);
           if (_dpartition_weight_dstress[p])
             dtotal_rom_effective_strain_inc_dstress +=
                 _dpartition_weight_dstress[p] * _weights[p][t] * rom;
           if (dweights_dstress[p][t])
             dtotal_rom_effective_strain_inc_dstress +=
                 _partition_weights[p] * dweights_dstress[p][t] * rom;
         }
     }
   }
 
   if (_verbose)
   {
     Moose::err << std::setprecision(9);
     GenericReal<is_ad> environmental = 0.0;
     if (_environmental)
       environmental = (*_environmental)[_qp];
     Moose::err << "Verbose information from " << _name << ": \n";
     Moose::err << " dt: " << _dt << "\n";
     Moose::err << " old cell disl: " << _old_input_values[_cell_output_index] << "\n";
     Moose::err << " old wall disl: " << _old_input_values[_wall_output_index] << "\n";
     Moose::err << " initial stress (MPa): "
                << MetaPhysicL::raw_value(effective_trial_stress) * _stress_ucf << "\n";
     Moose::err << " temperature: " << MetaPhysicL::raw_value(_temperature[_qp]) << "\n";
     Moose::err << " environmental factor: " << MetaPhysicL::raw_value(environmental) << "\n";
     Moose::err << " calculated scalar strain value: " << MetaPhysicL::raw_value(scalar) << "\n";
     Moose::err << " trial stress into rom (MPa): " << MetaPhysicL::raw_value(trial_stress_mpa)
                << "\n";
     Moose::err << " old effective strain: " << _old_input_values[_strain_output_index] << "\n";
     Moose::err << " extrapolation: " << MetaPhysicL::raw_value(_extrapolation[_qp]) << "\n";
     Moose::err << " partition 2 weight: " << MetaPhysicL::raw_value(_second_partition_weight[_qp])
                << "\n";
     Moose::err << " weights by tile, partition 1: ";
     for (unsigned int t = 0; t < _num_tiles[0]; ++t)
       Moose::err << " (" << t << ", " << MetaPhysicL::raw_value(_weights[0][t]) << ") ";
     Moose::err << "\n";
     if (_num_partitions > 1)
     {
       Moose::err << " weights by tile, partition 2: ";
       for (unsigned int t = 0; t < _num_tiles[1]; ++t)
         Moose::err << " (" << t << ", " << MetaPhysicL::raw_value(_weights[1][t]) << ") ";
     }
     Moose::err << "\n";
     Moose::err << " nonstress weights by tile, partition 1: ";
     for (unsigned int t = 0; t < _num_tiles[0]; ++t)
       Moose::err << " (" << t << ", " << MetaPhysicL::raw_value(_non_stress_weights[0][t]) << ") ";
     Moose::err << "\n";
     if (_num_partitions > 1)
     {
       Moose::err << " nonstress weights by tile, partition 2: ";
       for (unsigned int t = 0; t < _num_tiles[1]; ++t)
         Moose::err << " (" << t << ", " << MetaPhysicL::raw_value(_non_stress_weights[1][t])
                    << ") ";
     }
     Moose::err << "\n";
     Moose::err << " effective strain increment: "
                << MetaPhysicL::raw_value(total_rom_effective_strain_inc) << std::endl;
   }
 
   _creep_rate[_qp] = total_rom_effective_strain_inc / _dt;
   _derivative = dtotal_rom_effective_strain_inc_dstress * dtrial_stress_dscalar - 1.0;
 
   if (!this->_apply_strain)
   {
     if (_verbose)
       Moose::err << "    Strain not applied due to apply_strain input parameter!" << std::endl;
     _derivative = 1.0;
     return 0.0;
   }
   return total_rom_effective_strain_inc - scalar;
 }

◆ 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;
   };

◆ computeROM()

template<bool is_ad>

GenericReal< is_ad > LAROMANCEStressUpdateBaseTempl< is_ad >::computeROM	(	const unsigned int	tile,
		const unsigned int	partition,
		const unsigned	out_index,
		const bool	derivative = `false`
	)

protected

Computes the ROM calculated increment for a given output and tile.

Parameters

tile	Tile index
partition	Partition index
out_index	Output index
derivative	Optional flag to return derivative of ROM increment with respect to stress.

Returns: ROM computed increment

Definition at line 956 of file LAROMANCEStressUpdateBase.C.

 {
   // Update due to new stress
   _rom_inputs[p][t][_stress_input_index] =
       normalizeInput(_input_values[_stress_input_index],
                      _transform[p][t][out_index][_stress_input_index],
                      _transform_coefs[p][t][out_index][_stress_input_index],
                      _transformed_normalization_limits[p][t][out_index][_stress_input_index]);
 
   buildPolynomials(
       p, _rom_inputs[p][t][_stress_input_index], _polynomial_inputs[p][t][_stress_input_index]);
 
   // Compute ROM values
   const GenericReal<is_ad> rom_output =
       computeValues(p, _precomputed_vals[p][t], _polynomial_inputs[p][t]);
 
   // Return converted output if not derivative
   if (!derivative)
     return convertOutput(p, _old_input_values, rom_output, out_index);
 
   const GenericReal<is_ad> drom_input =
       normalizeInput(_input_values[_stress_input_index],
                      _transform[p][t][out_index][_stress_input_index],
                      _transform_coefs[p][t][out_index][_stress_input_index],
                      _transformed_normalization_limits[p][t][out_index][_stress_input_index],
                      derivative);
 
   std::vector<GenericReal<is_ad>> dpolynomial_inputs(_degree[p], 0.0);
   buildPolynomials(
       p, _rom_inputs[p][t][_stress_input_index], dpolynomial_inputs, drom_input, derivative);
 
   const GenericReal<is_ad> drom_output = computeValues(
       p, _precomputed_vals[p][t], _polynomial_inputs[p][t], dpolynomial_inputs, derivative);
   return convertOutput(p, _old_input_values, rom_output, out_index, drom_output, derivative);
 }

◆ computeStrainEnergyRateDensity()

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 LAROMANCEStressUpdateBaseTempl< is_ad >::computeStressFinalize ( const GenericRankTwoTensor< is_ad > & )

overrideprotectedvirtual

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

Parameters

inelasticStrainIncrement Inelastic strain increment

Reimplemented from RadialReturnCreepStressUpdateBaseTempl< is_ad >.

Definition at line 1182 of file LAROMANCEStressUpdateBase.C.

 {
   _cell_dislocation_increment = 0.0;
   _wall_dislocation_increment = 0.0;
   if (_input_values[_stress_input_index])
   {
     precomputeROM(_cell_output_index);
     for (unsigned int p = 0; p < _num_partitions; ++p)
       if (_partition_weights[p])
         for (unsigned int t = 0; t < _num_tiles[p]; ++t)
           if (_weights[p][t])
             _cell_dislocation_increment +=
                 _partition_weights[p] * _weights[p][t] * computeROM(t, p, _cell_output_index);
     precomputeROM(_wall_output_index);
     for (unsigned int p = 0; p < _num_partitions; ++p)
       if (_partition_weights[p])
         for (unsigned int t = 0; t < _num_tiles[p]; ++t)
           if (_weights[p][t])
             _wall_dislocation_increment +=
                 _partition_weights[p] * _weights[p][t] * computeROM(t, p, _wall_output_index);
   }
 
   _cell_rate[_qp] = _cell_dislocation_increment / _dt;
   _wall_rate[_qp] = _wall_dislocation_increment / _dt;
   _cell_dislocations[_qp] = _old_input_values[_cell_output_index] + _cell_dislocation_increment;
   _wall_dislocations[_qp] = _old_input_values[_wall_output_index] + _wall_dislocation_increment;
 
   // For (possibly) substepping.
   _plastic_strain_increment += MetaPhysicL::raw_value(plastic_strain_increment);
 
   // Prevent the ROM from calculating and proceeding with negative dislocations
   if (_apply_strain && (_cell_dislocations[_qp] < 0.0 || _wall_dislocations[_qp] < 0.0))
     mooseException("In ",
                    _name,
                    ": Negative disclocation density calculated for cell (old : ",
                    MetaPhysicL::raw_value(_old_input_values[_cell_output_index]),
                    " increment: ",
                    MetaPhysicL::raw_value(_cell_dislocation_increment),
                    " value: ",
                    MetaPhysicL::raw_value(_cell_dislocations[_qp]),
                    ") or wall (old : ",
                    MetaPhysicL::raw_value(_old_input_values[_wall_output_index]),
                    " increment: ",
                    MetaPhysicL::raw_value(_wall_dislocation_increment),
                    " value: ",
                    MetaPhysicL::raw_value(_wall_dislocations[_qp]),
                    ").");
 
   if (_verbose)
   {
     Moose::err << " Finalized ROM output\n";
     Moose::err << "  effective creep strain increment: "
                << std::sqrt(2.0 / 3.0 *
                             MetaPhysicL::raw_value(_plastic_strain_increment.doubleContraction(
                                 _plastic_strain_increment)))
                << "\n";
     Moose::err << "  total effective creep strain: "
                << std::sqrt(2.0 / 3.0 *
                             MetaPhysicL::raw_value(this->_creep_strain[_qp].doubleContraction(
                                 this->_creep_strain[_qp])))
                << "\n";
     Moose::err << "  creep rate: " << MetaPhysicL::raw_value(_creep_rate[_qp]) << "\n";
     Moose::err << "  cell dislocation rate: " << MetaPhysicL::raw_value(_cell_rate[_qp]) << "\n";
     Moose::err << "  wall dislocation rate: " << MetaPhysicL::raw_value(_wall_rate[_qp]) << "\n";
     Moose::err << "  new cell dislocations: " << MetaPhysicL::raw_value(_cell_dislocations[_qp])
                << "\n";
     Moose::err << "  new wall dislocations: " << MetaPhysicL::raw_value(_wall_dislocations[_qp])
                << "\n"
                << std::endl;
   }
 
   RadialReturnCreepStressUpdateBaseTempl<is_ad>::computeStressFinalize(
       MetaPhysicL::raw_value(_plastic_strain_increment));
 }

◆ computeStressInitialize()

template<bool is_ad>

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

overrideprotectedvirtual

Perform any necessary initialization before return mapping iterations.

Parameters

effective_trial_stress	Effective trial stress
elasticityTensor	Elasticity tensor

Reimplemented from RadialReturnStressUpdateTempl< is_ad >.

Definition at line 550 of file LAROMANCEStressUpdateBase.C.

 {
   RadialReturnCreepStressUpdateBaseTempl<is_ad>::computeStressInitialize(effective_trial_stress,
                                                                          elasticity_tensor);
   // Previous substep creep strain
   RankTwoTensor creep_strain_substep = this->_creep_strain_old[_qp] + _plastic_strain_increment;
 
   // Prepare old values
   _old_input_values[_cell_output_index] =
       _cell_function
           ? _cell_function->value(_t, _q_point[_qp])
           : (_cell_dislocations_old[_qp] + MetaPhysicL::raw_value(_cell_dislocations_step[_qp]));
   _old_input_values[_wall_output_index] =
       _wall_function
           ? _wall_function->value(_t, _q_point[_qp])
           : (_wall_dislocations_old[_qp] + MetaPhysicL::raw_value(_wall_dislocations_step[_qp]));
   _old_input_values[_strain_output_index] =
       _creep_strain_old_forcing_function
           ? _creep_strain_old_forcing_function->value(_t, _q_point[_qp])
           : MetaPhysicL::raw_value(
                 std::sqrt((creep_strain_substep).doubleContraction(creep_strain_substep) / 1.5));
 
   // Prepare input
   _input_values[_cell_input_index] = _old_input_values[_cell_output_index];
   _input_values[_wall_input_index] = _old_input_values[_wall_output_index];
   _input_values[_stress_input_index] = _stress_function ? _stress_function->value(_t, _q_point[_qp])
                                                         : effective_trial_stress * _stress_ucf;
   _input_values[_old_strain_input_index] = _old_input_values[_strain_output_index];
   _input_values[_temperature_input_index] = _temperature[_qp];
   if (_environmental)
     _input_values[_environmental_input_index] = (*_environmental)[_qp];
 
   // Determine tile weight and check if input is in range
   for (unsigned int p = 0; p < _num_partitions; ++p)
     std::fill(_non_stress_weights[p].begin(), _non_stress_weights[p].end(), 1.0);
   for (unsigned int i = 0; i < _num_inputs; i++)
     if (i != _stress_input_index)
       computeTileWeight(_non_stress_weights, _input_values[i], i);
 
   // Precompute transformed input and prebuild polynomials for inputs other than strain
   precomputeROM(_strain_output_index);
 }

◆ 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

◆ computeTileWeight()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::computeTileWeight	(	std::vector< std::vector< GenericReal< is_ad >>> &	weights,
		GenericReal< is_ad > &	input,
		const unsigned int	in_index,
		const bool	derivative = `false`
	)

protected

Compute the contribution (weight) of each tile in each partition, based on the input and tile boundaries (in terms of input domain).

Parameters

weights	Weights for each tile
input	Input value
in_index	Input index
derivative	Optional flag to return derivative of the sigmoid w.r.t. the input

Definition at line 597 of file LAROMANCEStressUpdateBase.C.

 {
   mooseAssert(std::isfinite(MetaPhysicL::raw_value(input)),
               "computeTileWeight input must be finite");
 
   for (unsigned int p = 0; p < _num_partitions; ++p)
   {
     for (unsigned int t = 0; t < _num_tiles[p]; ++t)
     {
       // If input is within a specfic tile's window of applicability
       if (input >= _input_limits[p][t][in_index][0] && input <= _input_limits[p][t][in_index][1])
       {
         // If there is not more than one tile for this input, set derivative of
         // weight to zero
         if (_tiling[p][in_index] == 1)
         {
           if (derivative)
             weights[p][t] = 0.0;
         }
         else
         {
           // Flag to ensure weights are applied only once
           bool overlap = false;
           for (unsigned int tt = 0; tt < _num_tiles[p]; ++tt)
           {
             if (!overlap && t != tt)
             {
               // If input is within another tile's window of applicability, check to see if inputs
               // place us in that tile and ensure the two tiles are different in the dimension of
               // interest
               if (areTilesNotIdentical(p, t, tt, in_index) && checkInTile(p, tt))
               {
                 overlap = true;
 
                 // If current tile is below the second tile's window of applicability
                 if (_input_limits[p][t][in_index][0] < _input_limits[p][tt][in_index][0] &&
                     _input_limits[p][t][in_index][1] > _input_limits[p][tt][in_index][0])
                 {
                   weights[p][t] *= sigmoid(_input_limits[p][tt][in_index][0],
                                            _input_limits[p][t][in_index][1],
                                            input,
                                            derivative);
                 }
                 // If current tile is above the second tile's window of applicability
                 else if (_input_limits[p][t][in_index][0] > _input_limits[p][tt][in_index][0] &&
                          _input_limits[p][t][in_index][0] < _input_limits[p][tt][in_index][1])
                 {
                   if (derivative)
                     weights[p][t] *= -sigmoid(_input_limits[p][t][in_index][0],
                                               _input_limits[p][tt][in_index][1],
                                               input,
                                               derivative);
                   else
                     weights[p][t] *= (1.0 - sigmoid(_input_limits[p][t][in_index][0],
                                                     _input_limits[p][tt][in_index][1],
                                                     input));
                 }
               }
             }
           }
 
           // If not overlapping, weight is not updated, and the derivative of tile weight is set to
           // zero
           if (!overlap && derivative)
             weights[p][t] = 0.0;
         }
       }
       // If input is below the lower tile limit
       else if (input < _input_limits[p][t][in_index][0])
       {
         // If the lower tile limit equals the lower global limit
         if (_input_limits[p][t][in_index][0] == _global_limits[in_index].first)
         {
           if (_window_failure[in_index].first == WindowFailure::EXTRAPOLATE)
           {
             if (derivative)
               weights[p][t] *= -sigmoid(0.0, _input_limits[p][t][in_index][0], input, derivative);
             else
               weights[p][t] *= (1.0 - sigmoid(0.0, _input_limits[p][t][in_index][0], input));
             input = _input_limits[p][t][in_index][0];
           }
           else if (_window_failure[in_index].first == WindowFailure::USELIMIT)
             input = _input_limits[p][t][in_index][0];
           else
           {
             weights[p][t] = 0.0;
             std::stringstream msg;
             msg << "In " << _name << ": " << _index_name[in_index]
                 << " input parameter with value (" << MetaPhysicL::raw_value(input)
                 << ") is out of lower global range (" << _input_limits[p][t][in_index][0] << ")";
 
             // Set input to edge of limit so it is found later in computePartitionWeights
             input = _input_limits[p][t][in_index][0];
 
             if (_window_failure[in_index].first == WindowFailure::WARN)
               mooseWarning(msg.str());
             else if (_window_failure[in_index].first == WindowFailure::ERROR)
               mooseError(msg.str());
             else if (_window_failure[in_index].first == WindowFailure::EXCEPTION)
               mooseException(msg.str());
             // if (_window_failure[in_index].first == WindowFailure::DONOTHING) <- nothing is done
           }
         }
         // if input below tile limit, update weight of tile to be zero
         else
           weights[p][t] = 0.0;
       }
       // If input is above the upper tile limit
       else if (input > _input_limits[p][t][in_index][1])
       {
         if (_input_limits[p][t][in_index][1] == _global_limits[in_index].second)
         {
           if (_window_failure[in_index].second == WindowFailure::EXTRAPOLATE)
             mooseError("Internal error. Extrapolate not appropriate for upper bound");
           else if (_window_failure[in_index].second == WindowFailure::USELIMIT)
             input = _input_limits[p][t][in_index][1];
           else
           {
             weights[p][t] = 0.0;
             std::stringstream msg;
             msg << "In " << _name << ": " << _index_name[in_index]
                 << " input parameter with value (" << MetaPhysicL::raw_value(input)
                 << ") is out of upper global range (" << _input_limits[p][t][in_index][1] << ")";
 
             // Set input to edge of limit so it is found later in computePartitionWeights
             input = _input_limits[p][t][in_index][1];
 
             if (_window_failure[in_index].second == WindowFailure::WARN)
               mooseWarning(msg.str());
             else if (_window_failure[in_index].second == WindowFailure::ERROR)
               mooseError(msg.str());
             else if (_window_failure[in_index].second == WindowFailure::EXCEPTION)
               mooseException(msg.str());
             // if (_window_failure[in_index].second == WindowFailure::DONOTHING) <- nothing is done
           }
         }
         // if input above tile limit, update weight of tile to be zero
         else
           weights[p][t] = 0.0;
       }
       // If input is outside window of applicability, weight is zero
       else
         mooseError("In ", _name, ": Internal error. Outside input limits, input=", input);
     }
   }
 }

◆ computeTimeStepLimit()

template<bool is_ad>

Real LAROMANCEStressUpdateBaseTempl< is_ad >::computeTimeStepLimit ( )

overrideprotectedvirtual

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

Returns: Limiting time step

Reimplemented from RadialReturnStressUpdateTempl< is_ad >.

Definition at line 1260 of file LAROMANCEStressUpdateBase.C.

 {
   Real limited_dt = RadialReturnCreepStressUpdateBaseTempl<is_ad>::computeTimeStepLimit();
 
   Real cell_strain_inc = std::abs(MetaPhysicL::raw_value(_cell_dislocation_increment));
   if (cell_strain_inc && _old_input_values[_cell_output_index])
     limited_dt = std::min(limited_dt,
                           _dt * _max_cell_increment * _old_input_values[_cell_output_index] /
                               cell_strain_inc);
   Real wall_strain_inc = std::abs(MetaPhysicL::raw_value(_wall_dislocation_increment));
   if (wall_strain_inc && _old_input_values[_wall_output_index])
     limited_dt = std::min(limited_dt,
                           _dt * _max_wall_increment * _old_input_values[_wall_output_index] /
                               wall_strain_inc);
 
   return limited_dt;
 }

◆ computeValues()

template<bool is_ad>

GenericReal< is_ad > LAROMANCEStressUpdateBaseTempl< is_ad >::computeValues	(	const unsigned int	p,
		const std::vector< GenericReal< is_ad >> &	precomputed,
		const std::vector< std::vector< GenericReal< is_ad >>> &	polynomial_inputs,
		const std::vector< GenericReal< is_ad >> &	dpolynomial_inputs = `{}`,
		const bool	derivative = `false`
	)

protected

Arranges the calculated Legendre polynomials into the proper oder and multiplies the Legendre polynomials by the ROM coefficients to compute the predicted output values.

This method only manipulates the stress input, with all others handled in precomputeValues

Parameters

p	Partition index
precomputed	Precomputed multiplication of polynomials
polynomial_inputs	Vector of Legendre polynomial transformation
dpolynomial_inputs	Vector of derivative of Legendre polynomial transformation with respect to stress
derivative	Optional flag to return derivative of converted computed values with respect to stress.

Returns: ROM output

Definition at line 1050 of file LAROMANCEStressUpdateBase.C.

 {
   GenericReal<is_ad> rom_output = 0.0;
   for (unsigned int c = 0; c < _num_coefs[p]; ++c)
   {
     if (!derivative)
       rom_output +=
           precomputed[c] *
           polynomial_inputs[_stress_input_index]
                            [_makeframe_helper[p][c + _num_coefs[p] * _stress_input_index]];
 
     else
       rom_output +=
           precomputed[c] *
           dpolynomial_inputs[_makeframe_helper[p][c + _num_coefs[p] * _stress_input_index]];
   }
   return rom_output;
 }

◆ 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);
 }

◆ convertOutput()

template<bool is_ad>

GenericReal< is_ad > LAROMANCEStressUpdateBaseTempl< is_ad >::convertOutput	(	const unsigned int	p,
		const std::vector< Real > &	old_input_values,
		const GenericReal< is_ad > &	rom_output,
		const unsigned	out_index,
		const GenericReal< is_ad > &	drom_output = `0.0`,
		const bool	derivative = `false`
	)

protectedvirtual

Computes the output variable increments from the ROM predictions by bringing out of the normalized map to the actual physical values.

Parameters

p	Partition index
old_input_values	Previous timestep values of ROM inputs
rom_output	Outputs from ROM
out_index	Output index
drom_output	Derivative of output with respect to stress
derivative	Optional flag to return derivative of output with respect to stress.

Returns: Converted ROM output

Definition at line 1076 of file LAROMANCEStressUpdateBase.C.

 {
   if (out_index == _strain_output_index)
   {
     if (derivative)
       return std::exp(rom_output) * _dt * drom_output;
     else
       return std::exp(rom_output) * _dt;
   }
 
   if (derivative)
     return 0.0;
 
   GenericReal<is_ad> expout = std::exp(rom_output);
   mooseAssert(expout > 0.0, "ROM calculated strain increment must be positive");
 
   if (expout > _cutoff[p])
     expout -= _cutoff[p];
   else
     expout = -_cutoff[p] * _cutoff[p] / expout + _cutoff[p];
 
   return -expout * old_input_values[out_index] * _dt;
 }

◆ convertValue()

template<bool is_ad>

template<typename T >

void LAROMANCEStressUpdateBaseTempl< is_ad >::convertValue	(	T &	x,
		const ROMInputTransform	transform,
		const Real	coef,
		const bool	derivative = `false`
	)

inlineprotected

Convert input based on the transform type.

Parameters

x	Input value
transform	Enum declaring transform to be performed
coef	Coefficient applied during transformation
derivative	Optional flag to return derivative of the sigmoid w.r.t. the input

Definition at line 275 of file LAROMANCEStressUpdateBase.h.

   {
     if (transform == ROMInputTransform::EXP)
     {
       mooseAssert(coef != 0, "Coefficient must not be zero.");
       if (derivative)
         x = std::exp(x / coef) / coef;
       else
         x = std::exp(x / coef);
     }
     else if (transform == ROMInputTransform::LOG)
     {
       mooseAssert(x + coef > 0, "Sum must be greater than 0.");
       if (derivative)
         x = 1.0 / (x + coef);
       else
         x = std::log(x + coef);
     }
     else if (transform == ROMInputTransform::LINEAR)
     {
       mooseAssert(coef == 0.0, "Coefficient cannot be supplied with linear transformation");
       if (derivative)
         x = 1.0;
     }
   }

◆ coupledGenericValue() [1/3]

template<bool is_ad>

const GenericVariableValue< is_ad > & Material::coupledGenericValue

protected

◆ coupledGenericValue() [2/3]

template<bool is_ad>

const GenericVariableValue< false > & Material::coupledGenericValue

protected

◆ coupledGenericValue() [3/3]

template<bool is_ad>

const GenericVariableValue< true > & Material::coupledGenericValue

protected

◆ 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;
   }

◆ exportJSON()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::exportJSON ( )

protectedvirtual

Reimplemented in LAROMANCEPartitionStressUpdateBaseTempl< is_ad >.

Definition at line 283 of file LAROMANCEStressUpdateBase.C.

Referenced by LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::exportJSON().

 {
   _json["strain_cutoff"] = getStrainCutoff();
   _json["transform"] = getTransform();
   _json["transform_coefs"] = getTransformCoefs();
   _json["input_limits"] = getInputLimits();
   _json["normalization_limits"] = getNormalizationLimits();
   _json["coefs"] = getCoefs();
   _json["tiling"] = getTilings();
   _json["cutoff"] = getStrainCutoff();
 }

◆ getCoefs()

template<bool is_ad>

virtual std::vector<std::vector<std::vector<std::vector<Real> > > > LAROMANCEStressUpdateBaseTempl< is_ad >::getCoefs ( )

inlineprotectedvirtual

Reimplemented in LAROMANCE3TileTestTempl< is_ad >, and SS316HLAROMANCEStressUpdateTestTempl< is_ad >.

Definition at line 423 of file LAROMANCEStressUpdateBase.h.

   {
     checkJSONKey("coefs");
     return _json["coefs"].template get<std::vector<std::vector<std::vector<std::vector<Real>>>>>();
   }

◆ getInputLimits()

template<bool is_ad>

virtual std::vector<std::vector<std::vector<std::vector<Real> > > > LAROMANCEStressUpdateBaseTempl< is_ad >::getInputLimits ( )

inlineprotectedvirtual

Reimplemented in LAROMANCE3TileTestTempl< is_ad >, and SS316HLAROMANCEStressUpdateTestTempl< is_ad >.

Definition at line 411 of file LAROMANCEStressUpdateBase.h.

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

   {
     checkJSONKey("input_limits");
     return _json["input_limits"]
         .template get<std::vector<std::vector<std::vector<std::vector<Real>>>>>();
   }

◆ getMakeFrameHelper()

template<bool is_ad>

std::vector< unsigned int > LAROMANCEStressUpdateBaseTempl< is_ad >::getMakeFrameHelper ( const unsigned int p ) const

protected

Definition at line 1168 of file LAROMANCEStressUpdateBase.C.

 {
   std::vector<unsigned int> v(_num_coefs[p] * _num_inputs);
 
   for (unsigned int numcoeffs = 0; numcoeffs < _num_coefs[p]; ++numcoeffs)
     for (unsigned int invar = 0; invar < _num_inputs; ++invar)
       v[numcoeffs + _num_coefs[p] * invar] =
           numcoeffs / MathUtils::pow(_degree[p], invar) % _degree[p];
 
   return v;
 }

◆ getNormalizationLimits()

template<bool is_ad>

virtual std::vector<std::vector<std::vector<std::vector<Real> > > > LAROMANCEStressUpdateBaseTempl< is_ad >::getNormalizationLimits ( )

inlineprotectedvirtual

Reimplemented in LAROMANCE3TileTestTempl< is_ad >.

Definition at line 391 of file LAROMANCEStressUpdateBase.h.

   {
     if (_json.contains("normalization_limits"))
       return _json["normalization_limits"]
           .template get<std::vector<std::vector<std::vector<std::vector<Real>>>>>();
 
     return getInputLimits();
   }

◆ getStrainCutoff()

template<bool is_ad>

virtual std::vector<Real> LAROMANCEStressUpdateBaseTempl< is_ad >::getStrainCutoff ( )

inlineprotectedvirtual

Reimplemented in LAROMANCE3TileTestTempl< is_ad >, and SS316HLAROMANCEStressUpdateTestTempl< is_ad >.

Definition at line 449 of file LAROMANCEStressUpdateBase.h.

   {
     checkJSONKey("cutoff");
     return _json["cutoff"].template get<std::vector<Real>>();
   }

◆ 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;
   }

◆ getTilings()

template<bool is_ad>

virtual std::vector<std::vector<unsigned int> > LAROMANCEStressUpdateBaseTempl< is_ad >::getTilings ( )

inlineprotectedvirtual

Reimplemented in LAROMANCE3TileTestTempl< is_ad >.

Definition at line 434 of file LAROMANCEStressUpdateBase.h.

   {
     if (_json.contains("tiling"))
       return _json["tiling"].template get<std::vector<std::vector<unsigned int>>>();
 
     if (_environmental)
       return {{1, 1, 1, 1, 1, 1}};
     return {{1, 1, 1, 1, 1}};
   };

◆ getTransform()

template<bool is_ad>

virtual std::vector<std::vector<std::vector<std::vector<ROMInputTransform> > > > LAROMANCEStressUpdateBaseTempl< is_ad >::getTransform ( )

inlineprotectedvirtual

Reimplemented in LAROMANCE3TileTestTempl< is_ad >, and SS316HLAROMANCEStressUpdateTestTempl< is_ad >.

Definition at line 349 of file LAROMANCEStressUpdateBase.h.

   {
     checkJSONKey("transform");
     return _json["transform"]
         .template get<std::vector<std::vector<std::vector<std::vector<ROMInputTransform>>>>>();
   }

◆ getTransformCoefs()

template<bool is_ad>

virtual std::vector<std::vector<std::vector<std::vector<Real> > > > LAROMANCEStressUpdateBaseTempl< is_ad >::getTransformCoefs ( )

inlineprotectedvirtual

Reimplemented in LAROMANCE3TileTestTempl< is_ad >, and SS316HLAROMANCEStressUpdateTestTempl< is_ad >.

Definition at line 372 of file LAROMANCEStressUpdateBase.h.

   {
     checkJSONKey("transform_coefs");
     return _json["transform_coefs"]
         .template get<std::vector<std::vector<std::vector<std::vector<Real>>>>>();
   }

◆ getTransformedLimits()

template<bool is_ad>

std::vector< std::vector< std::vector< std::vector< Real > > > > LAROMANCEStressUpdateBaseTempl< is_ad >::getTransformedLimits	(	const unsigned int	p,
		const std::vector< std::vector< std::vector< Real >>>	limits
	)

protected

Definition at line 1139 of file LAROMANCEStressUpdateBase.C.

 {
   std::vector<std::vector<std::vector<std::vector<Real>>>> transformed_limits(
       _num_tiles[p],
       std::vector<std::vector<std::vector<Real>>>(
           _num_outputs, std::vector<std::vector<Real>>(_num_inputs, std::vector<Real>(3))));
   for (unsigned int t = 0; t < _num_tiles[p]; ++t)
   {
     for (unsigned int o = 0; o < _num_outputs; ++o)
     {
       for (unsigned int i = 0; i < _num_inputs; ++i)
       {
         for (unsigned int k = 0; k < 2; ++k)
         {
           transformed_limits[t][o][i][k] = limits[t][i][k];
           convertValue(
               transformed_limits[t][o][i][k], _transform[p][t][o][i], _transform_coefs[p][t][o][i]);
         }
         transformed_limits[t][o][i][2] =
             (transformed_limits[t][o][i][1] - transformed_limits[t][o][i][0]) / 2.0;
       }
     }
   }
   return transformed_limits;
 }

◆ 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;
   }

◆ initialSetup()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::initialSetup ( )

overrideprotectedvirtual

Reimplemented from Material.

Reimplemented in LAROMANCEPartitionStressUpdateBaseTempl< is_ad >.

Definition at line 315 of file LAROMANCEStressUpdateBase.C.

Referenced by LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::initialSetup().

 {
   // export models that are compiled in
   if (this->isParamValid("export_model"))
   {
     exportJSON();
     std::ofstream out(this->template getParam<FileName>("export_model").c_str());
     out << _json;
   }
 
   // Pull in relevant ROM information and do sanity checks
   _transform = getTransform();
   _transform_coefs = getTransformCoefs();
   _input_limits = getInputLimits();
   _normalization_limits = getNormalizationLimits();
   _coefs = getCoefs();
   _tiling = getTilings();
   _cutoff = getStrainCutoff();
   // resize containers to be filled later based on partition dimension
   // and immediately run some sanity checks
   _num_partitions = _transform.size();
   if (_num_partitions < 1 || _num_partitions > 2)
     mooseError(
         "In ", _name, ": First dimension of getTransform must be either size one or size two");
   if (_transform[0].size() < 1)
     mooseError("In ", _name, ": Transform is not the correct shape");
   _num_outputs = _transform[0][0].size();
   if (_num_outputs != 3)
     mooseError("In ",
                _name,
                ": ",
                _num_outputs,
                " outputs detected. Three and only three outputs are currently supported.");
 
   _num_inputs = _transform[0][0][0].size();
   if (_num_inputs != 5 && _num_inputs != 6)
     mooseError("In ",
                _name,
                ": ",
                _num_inputs,
                " inputs detected. Only five or six inputs currently supported.");
   if (_num_inputs == 6 && !_environmental)
     this->paramError(
         "environmental_factor",
         "Number of ROM inputs indicate environmental factor is required to be coupled.");
   if (_num_inputs != 6 && _environmental)
     this->paramError("environmental_factor",
                      "Number of ROM inputs indicate environmental factor is not implemented, but "
                      "environmental factor coupled.");
   _num_tiles.resize(_num_partitions);
   _num_coefs.resize(_num_partitions);
   _degree.resize(_num_partitions);
   _precomputed_vals.resize(_num_partitions);
   _rom_inputs.resize(_num_partitions);
   _polynomial_inputs.resize(_num_partitions);
   _non_stress_weights.resize(_num_partitions);
   _weights.resize(_num_partitions);
   _partition_weights.resize(_num_partitions);
   _dpartition_weight_dstress.resize(_num_partitions);
   _transformed_normalization_limits.resize(_num_partitions);
   _makeframe_helper.resize(_num_partitions);
   _global_limits.resize(_num_inputs);
   // temporarily fill global limits with extreme numerical values, to later update
   for (unsigned int i = 0; i < _num_inputs; ++i)
     _global_limits[i] = {std::numeric_limits<Real>::max(), 0.0};
 
   // start loop over partitions to perform sanity checks, set global limits,
   // and print global configurations
   if (_transform.size() != _num_partitions || _transform_coefs.size() != _num_partitions ||
       _input_limits.size() != _num_partitions || _normalization_limits.size() != _num_partitions ||
       _coefs.size() != _num_partitions || _tiling.size() != _num_partitions ||
       _cutoff.size() != _num_partitions)
     mooseError(
         "In ", _name, ": one of the ROM inputs does not have the correct number of partitions");
 
   for (unsigned int p = 0; p < _num_partitions; ++p)
   {
     _num_tiles[p] = _transform[p].size();
     if (!_num_tiles[p])
       mooseError("In ", _name, ": No tiles detected. Double check your ROM input");
 
     bool correct_shape = true;
     if (_transform[p].size() != _num_tiles[p] || _transform_coefs[p].size() != _num_tiles[p] ||
         _input_limits[p].size() != _num_tiles[p] ||
         _normalization_limits[p].size() != _num_tiles[p] || _coefs[p].size() != _num_tiles[p])
       correct_shape = false;
     if (_tiling[p].size() != _num_inputs)
       correct_shape = false;
     if (_coefs[p][0].size() == 0)
       correct_shape = false;
     _num_coefs[p] = _coefs[p][0][0].size();
 
     // start loop over tiles to perform sanity checks.
     for (unsigned int t = 0; t < _num_tiles[p]; ++t)
     {
       if (_transform[p][t].size() != _num_outputs ||
           _transform_coefs[p][t].size() != _num_outputs || _coefs[p][t].size() != _num_outputs)
         correct_shape = false;
       for (unsigned int o = 0; o < _num_outputs; ++o)
         if (_transform[p][t][o].size() != _num_inputs ||
             _transform_coefs[p][t][o].size() != _num_inputs ||
             _coefs[p][t][o].size() != _num_coefs[p])
           correct_shape = false;
       if (_input_limits[p][t].size() != _num_inputs ||
           _normalization_limits[p][t].size() != _num_inputs)
         correct_shape = false;
       for (unsigned int i = 0; i < _num_inputs; ++i)
         if (_input_limits[p][t][i].size() != 2 || _normalization_limits[p][t][i].size() != 2)
           correct_shape = false;
     }
 
     if (!correct_shape)
       mooseError("In ", _name, ": ROM data is not the right shape.");
 
     _degree[p] = std::pow(_num_coefs[p], 1.0 / _num_inputs);
     if (!_degree[p] || _degree[p] > 4)
       mooseError("In ", _name, ": degree must be 1, 2, 3 or 4.");
 
     // Check input limits and find global limits of all partitions. Note that this will return the
     // extremes of the model! If the model is not flat along one input limit, then global limits
     // will not truely capture the mutli-dimensionality of the problem. Consequently, the user may
     // find input values that result in errors in the partition weight computation.
     for (unsigned int i = 0; i < _num_inputs; ++i)
     {
       for (unsigned int t = 0; t < _num_tiles[p]; ++t)
       {
         if (_input_limits[p][t][i][0] >= _input_limits[p][t][i][1])
           mooseError("In ", _name, ": Input limits are ordered incorrectly");
         _global_limits[i].first = std::min(_global_limits[i].first, _input_limits[p][t][i][0]);
         _global_limits[i].second = std::max(_global_limits[i].second, _input_limits[p][t][i][1]);
       }
     }
 
     // Precompute helper containers
     _transformed_normalization_limits[p] = getTransformedLimits(p, _normalization_limits[p]);
     _makeframe_helper[p] = getMakeFrameHelper(p);
 
     // Prepare containers for each partition
     _precomputed_vals[p].resize(_num_tiles[p], std::vector<GenericReal<is_ad>>(_num_coefs[p]));
     _rom_inputs[p].resize(_num_tiles[p], std::vector<GenericReal<is_ad>>(_num_inputs));
     _polynomial_inputs[p].resize(_num_tiles[p],
                                  std::vector<std::vector<GenericReal<is_ad>>>(
                                      _num_inputs, std::vector<GenericReal<is_ad>>(_degree[p])));
     _non_stress_weights[p].resize(_num_tiles[p]);
     _weights[p].resize(_num_tiles[p], 0);
   }
   // Prepare containers independent of partition
   _input_values.resize(_num_inputs);
 
   if (_verbose)
   {
     Moose::err << "ROM model info: " << _name << "\n";
     Moose::err << " number of tiles, partition 1: " << _num_tiles[0] << "\n";
     if (_num_partitions > 1)
       Moose::err << " number of tiles, partition 2: " << _num_tiles[1] << "\n";
     Moose::err << " number of outputs: " << _num_outputs << "\n";
     Moose::err << " number of inputs: " << _num_inputs << "\n";
     Moose::err << " degree (max Legendre degree + constant), partition 1: " << _degree[0] << "\n";
     if (_num_partitions > 1)
       Moose::err << " degree (max Legendre degree + constant), partition 2: " << _degree[1] << "\n";
     Moose::err << " number of coefficients, partition 1: " << _num_coefs[0] << "\n";
     if (_num_partitions > 1)
       Moose::err << " number of coefficients, partition 2: " << _num_coefs[1] << "\n";
     Moose::err << " Global limits:\n  cell dislocations ("
                << _global_limits[_cell_input_index].first << " - "
                << _global_limits[_cell_input_index].second << ")\n";
     Moose::err << "  wall dislocations (" << _global_limits[_wall_input_index].first << " - "
                << _global_limits[_wall_input_index].second << ")\n";
     Moose::err << "  Stress (" << _global_limits[_stress_input_index].first << " - "
                << _global_limits[_stress_input_index].second << ")\n";
     Moose::err << "  Old strain (" << _global_limits[_old_strain_input_index].first << " - "
                << _global_limits[_old_strain_input_index].second << ")\n";
     Moose::err << "  Temperature (" << _global_limits[_temperature_input_index].first << " - "
                << _global_limits[_temperature_input_index].second << ")\n";
     if (_environmental)
       Moose::err << "  Environmental factor (" << _global_limits[_environmental_input_index].first
                  << " - " << _global_limits[_environmental_input_index].second << ")\n";
     Moose::err << std::endl;
   }
 }

◆ initQpStatefulProperties()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::initQpStatefulProperties ( )

overrideprotectedvirtual

Reimplemented from RadialReturnCreepStressUpdateBaseTempl< is_ad >.

Definition at line 498 of file LAROMANCEStressUpdateBase.C.

 {
   MooseRandom rng;
   rng.seed(0, this->template getParam<unsigned int>("seed"));
 
   _cell_dislocations[_qp] = rng.randNormal(
       this->template getParam<Real>("initial_cell_dislocation_density"),
       this->template getParam<Real>("initial_cell_dislocation_density") *
           this->template getParam<Real>("cell_dislocations_normal_distribution_width"));
   _wall_dislocations[_qp] = rng.randNormal(
       this->template getParam<Real>("initial_wall_dislocation_density"),
       this->template getParam<Real>("initial_wall_dislocation_density") *
           this->template getParam<Real>("wall_dislocations_normal_distribution_width"));
 
   RadialReturnCreepStressUpdateBaseTempl<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.

128 { 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.

128 {}

◆ maximumPermissibleValue()

template<bool is_ad>

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

overrideprotectedvirtual

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 RadialReturnStressUpdateTempl< is_ad >.

Definition at line 517 of file LAROMANCEStressUpdateBase.C.

 {
   // Make maximum allowed scalar a little bit less than the deformation that would reduce the
   // trial stress to zero. This prevents negative trial stresses.
   return effective_trial_stress / this->_three_shear_modulus * 0.999999;
 }

◆ 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;
   }

◆ normalizeInput()

template<bool is_ad>

GenericReal< is_ad > LAROMANCEStressUpdateBaseTempl< is_ad >::normalizeInput	(	const GenericReal< is_ad > &	input,
		const ROMInputTransform	transform,
		const Real	transform_coef,
		const std::vector< Real > &	transformed_limits,
		const bool	derivative = `false`
	)

protected

Convert the input variables into the form expected by the ROM Legendre polynomials to have a normalized space from [-1, 1] so that every variable has equal weight.

Parameters

input	Input value
transform	ROMInputTransform enum indicating how the input is to be transformed
transform_coef	Transform coefficient for the given input
transformed_limits	Transformed limits for the given input
derivative	Optional flag to return derivative of converted input with respect to stress.

Returns: Converted input

Definition at line 997 of file LAROMANCEStressUpdateBase.C.

 {
   GenericReal<is_ad> x = input;
   convertValue(x, transform, transform_coef, derivative);
 
   // transformed_limits[2] = transformed_limits[1] - transformed_limits[0]
   if (derivative)
     return x / transformed_limits[2];
   else
     return (x - transformed_limits[0]) / transformed_limits[2] - 1.0;
 }

◆ outputIterationStep()

template<bool is_ad>

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

overrideprotectedvirtual

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 from SingleVariableReturnMappingSolutionTempl< is_ad >.

Definition at line 1377 of file LAROMANCEStressUpdateBase.C.

 {
   SingleVariableReturnMappingSolutionTempl<is_ad>::outputIterationStep(
       iter_output, effective_trial_stress, scalar, reference_residual);
   if (iter_output)
     *iter_output << " derivative: "
                  << MetaPhysicL::raw_value(computeDerivative(effective_trial_stress, scalar))
                  << std::endl;
 }

◆ outputIterationSummary()

template<bool is_ad>

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

overrideprotectedvirtual

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 1331 of file LAROMANCEStressUpdateBase.C.

 {
   if (iter_output)
   {
     *iter_output << "At element " << this->_current_elem->id() << " _qp=" << _qp << " Coordinates "
                  << _q_point[_qp] << " block=" << this->_current_elem->subdomain_id() << '\n';
     *iter_output << " dt " << _dt << " old cell disl: " << _old_input_values[_cell_output_index]
                  << " old wall disl: " << _old_input_values[_wall_output_index]
                  << " old effective strain: " << _old_input_values[_strain_output_index] << "\n";
 
     *iter_output << " temp: " << MetaPhysicL::raw_value(_temperature[_qp]) << " environmental: "
                  << (_environmental ? MetaPhysicL::raw_value((*_environmental)[_qp]) : 0.0)
                  << " trial stress into rom (MPa): "
                  << MetaPhysicL::raw_value(_input_values[_stress_input_index])
                  << " cell: " << MetaPhysicL::raw_value(_input_values[_cell_input_index])
                  << " wall: " << MetaPhysicL::raw_value(_input_values[_wall_input_index])
                  << " old strain: "
                  << MetaPhysicL::raw_value(_input_values[_old_strain_input_index]) << "\n";
     *iter_output << "  partition 2 weight: "
                  << MetaPhysicL::raw_value(_second_partition_weight[_qp]) << "\n";
     *iter_output << "  weights by tile, partition 1: ";
     for (unsigned int t = 0; t < _num_tiles[0]; ++t)
       *iter_output << " (" << t << ", " << MetaPhysicL::raw_value(_weights[0][t]) << ") ";
     *iter_output << "\n";
     *iter_output << "  weights by tile, partition 2: ";
     for (unsigned int t = 0; t < _num_tiles[1]; ++t)
       *iter_output << " (" << t << ", " << MetaPhysicL::raw_value(_weights[1][t]) << ") ";
     *iter_output << "\n";
     *iter_output << "  nonstress weights by tile, partition 1: ";
     for (unsigned int t = 0; t < _num_tiles[0]; ++t)
       *iter_output << " (" << t << ", " << MetaPhysicL::raw_value(_non_stress_weights[0][t])
                    << ") ";
     *iter_output << "\n";
     *iter_output << "  nonstress weights by tile, partition 2: ";
     for (unsigned int t = 0; t < _num_tiles[1]; ++t)
       *iter_output << " (" << t << ", " << MetaPhysicL::raw_value(_non_stress_weights[1][t])
                    << ") ";
     *iter_output << "\n";
   }
 
   SingleVariableReturnMappingSolutionTempl<is_ad>::outputIterationSummary(iter_output, total_it);
 }

◆ precomputeROM()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::precomputeROM ( const unsigned out_index )

protected

Precompute the ROM strain rate information for all inputs except for strain.

Strain will be computed in the radial return algorithm several times, while the remainder of the inputs remain constant.

Parameters

out_index Output index

Definition at line 926 of file LAROMANCEStressUpdateBase.C.

 {
   for (unsigned int p = 0; p < _num_partitions; ++p)
   {
     for (unsigned int t = 0; t < _num_tiles[p]; ++t)
     {
       // Only precompute for tiles that don't have zero weight
       if (_non_stress_weights[p][t])
       {
         for (unsigned int i = 0; i < _num_inputs; ++i)
         {
           if (i != _stress_input_index)
           {
             _rom_inputs[p][t][i] =
                 normalizeInput(_input_values[i],
                                _transform[p][t][out_index][i],
                                _transform_coefs[p][t][out_index][i],
                                _transformed_normalization_limits[p][t][out_index][i]);
             buildPolynomials(p, _rom_inputs[p][t][i], _polynomial_inputs[p][t][i]);
           }
         }
         precomputeValues(
             p, _coefs[p][t][out_index], _polynomial_inputs[p][t], _precomputed_vals[p][t]);
       }
     }
   }
 }

◆ precomputeValues()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::precomputeValues	(	const unsigned int	p,
		const std::vector< Real > &	coefs,
		const std::vector< std::vector< GenericReal< is_ad >>> &	polynomial_inputs,
		std::vector< GenericReal< is_ad >> &	precomputed
	)

protected

Arranges the calculated Legendre polynomials into the proper oder and multiplies the Legendre polynomials by the ROM coefficients to compute the predicted output values.

This method works with all inputs besides stress, while stress is handled by computeValues

Parameters

p	Partition index
coefs	Legendre polynomial coefficients
polynomial_inputs	Vector of transformed Legendre polynomial
precomputed	Vector that holds the precomputed ROM values

Definition at line 1033 of file LAROMANCEStressUpdateBase.C.

 {
   for (unsigned int c = 0; c < _num_coefs[p]; ++c)
   {
     precomputed[c] = coefs[c];
     for (unsigned int i = 0; i < _num_inputs; ++i)
       if (i != _stress_input_index)
         precomputed[c] *= polynomial_inputs[i][_makeframe_helper[p][c + _num_coefs[p] * i]];
   }
 }

◆ 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.

121 {

122 }

◆ 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.

123 { return true; }

◆ resetIncrementalMaterialProperties()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::resetIncrementalMaterialProperties ( )

overridevirtual

Reset material properties.

Useful for substepping with inelastic models.

Reimplemented from StressUpdateBaseTempl< is_ad >.

Definition at line 527 of file LAROMANCEStressUpdateBase.C.

 {
   _cell_dislocation_increment = 0.0;
   _wall_dislocation_increment = 0.0;
 
   _plastic_strain_increment.zero();
 
   _wall_dislocations_step[_qp] = 0.0;
   _cell_dislocations_step[_qp] = 0.0;
 }

◆ resetProperties()

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

inlinefinalvirtualinherited

Reimplemented from Material.

Definition at line 138 of file StressUpdateBase.h.

138 {}

◆ 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.

137 {}

◆ 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;
 }

◆ setupUnitConversionFactors()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::setupUnitConversionFactors ( const InputParameters & parameters )

protectedvirtual

Definition at line 304 of file LAROMANCEStressUpdateBase.C.

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

 {
   // Stress unit conversion factor
   const MooseUnits stress_unit_to("MPa");
   const MooseUnits stress_unit_from(parameters.get<std::string>("stress_unit"));
   _stress_ucf = stress_unit_to.convert(1, stress_unit_from);
 }

◆ sigmoid()

template<bool is_ad>

GenericReal< is_ad > LAROMANCEStressUpdateBaseTempl< is_ad >::sigmoid	(	const Real	lower,
		const Real	upper,
		const GenericReal< is_ad > &	val,
		const bool	derivative = `false`
	)

protected

Calculate the sigmoid function weighting for the input based on the limits.

Parameters

lower	Lower limit
upper	Upper limit
val	Input value
derivative	Optional flag to return derivative of the sigmoid w.r.t. the input

Returns: weight

Definition at line 1280 of file LAROMANCEStressUpdateBase.C.

 {
   mooseAssert(std::isfinite(MetaPhysicL::raw_value(val)), "Sigmoid value should must be infinite");
   mooseAssert(MetaPhysicL::raw_value(val) >= lower, "Input value must be greater than lower limit");
   mooseAssert(MetaPhysicL::raw_value(val) <= upper, "Input value must be greater than upper limit");
 
   // normalize val between 0 and 1, then shift between -1 and 1
   const GenericReal<is_ad> x = 2.0 * (val - lower) / (upper - lower) - 1.0;
 
   if (x == -1.0)
   {
     if (derivative)
       return 0.0;
     else
       return 1.0;
   }
   else if (x == 1.0)
   {
     if (derivative)
       return 0.0;
     else
       return 0.0;
   }
   else if (x < 1.0 && x > -1.0)
   {
     const GenericReal<is_ad> plus = std::exp(-2.0 / (1.0 + x));
     const GenericReal<is_ad> minus = std::exp(-2.0 / (1.0 - x));
 
     if (!derivative)
       return 1.0 - plus / (plus + minus);
 
     const GenericReal<is_ad> dplus_dx = plus * 2.0 / Utility::pow<2>(1.0 + x);
     const GenericReal<is_ad> dminus_dx = -minus * 2.0 / Utility::pow<2>(1.0 - x);
 
     return (plus * dminus_dx - dplus_dx * minus) / Utility::pow<2>(plus + minus) * 2.0 /
            (upper - lower);
   }
   else
     mooseError("Internal error: Sigmoid, value: x is out of bounds. val=",
                val,
                " low=",
                lower,
                " high=",
                upper);
 }

◆ storeIncrementalMaterialProperties()

template<bool is_ad>

void LAROMANCEStressUpdateBaseTempl< is_ad >::storeIncrementalMaterialProperties ( const unsigned int )

overridevirtual

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

Reimplemented from StressUpdateBaseTempl< is_ad >.

Definition at line 540 of file LAROMANCEStressUpdateBase.C.

 {
   _wall_dislocations_step[_qp] += _wall_dislocation_increment;
   _cell_dislocations_step[_qp] += _cell_dislocation_increment;
   _number_of_substeps[_qp] = total_number_of_substeps;
 }

◆ substeppingCapabilityEnabled()

template<bool is_ad>

bool LAROMANCEStressUpdateBaseTempl< is_ad >::substeppingCapabilityEnabled ( )

overrideprotectedvirtual

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 297 of file LAROMANCEStressUpdateBase.C.

 {
   return this->_use_substepping != RadialReturnStressUpdateTempl<is_ad>::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->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()

template<bool is_ad>

InputParameters LAROMANCEStressUpdateBaseTempl< is_ad >::validParams ( )

static

Definition at line 26 of file LAROMANCEStressUpdateBase.C.

Referenced by LAROMANCEPartitionStressUpdateBaseTempl< is_ad >::validParams(), and SS316HLAROMANCEStressUpdateTestTempl< is_ad >::validParams().

 {
   InputParameters params = RadialReturnCreepStressUpdateBaseTempl<is_ad>::validParams();
   params.addClassDescription("Base class to calculate the effective creep strain based on the "
                              "rates predicted by a material  specific Los Alamos Reduced Order "
                              "Model derived from a Visco-Plastic Self Consistent calculations.");
 
   params.addRequiredCoupledVar("temperature", "The coupled temperature (K)");
   params.addParam<MaterialPropertyName>("environmental_factor",
                                         "Optional coupled environmental factor");
 
   MooseEnum error_lower_limit_behavior("ERROR EXCEPTION WARN IGNORE DONTHING USELIMIT",
                                        "EXCEPTION");
   // Only allow ERROR and EXCEPTION on upper bounds
   MooseEnum error_upper_limit_behavior("ERROR EXCEPTION", "EXCEPTION");
   params.addParam<MooseEnum>(
       "cell_input_window_low_failure",
       error_lower_limit_behavior,
       "What to do if cell dislocation concentration is outside the lower global "
       "window of applicability.");
   params.addParam<MooseEnum>(
       "cell_input_window_high_failure",
       error_upper_limit_behavior,
       "What to do if cell dislocation concentration is outside the upper global "
       "window of applicability.");
   params.addParam<MooseEnum>("wall_input_window_low_failure",
                              error_lower_limit_behavior,
                              "What to do if wall dislocation concentration is outside the "
                              "lower global window of applicability.");
   params.addParam<MooseEnum>("wall_input_window_high_failure",
                              error_upper_limit_behavior,
                              "What to do if wall dislocation concentration is outside the "
                              "upper global window of applicability.");
   params.addParam<MooseEnum>(
       "old_strain_input_window_low_failure",
       error_lower_limit_behavior,
       "What to do if old strain is outside the lower global window of applicability.");
   params.addParam<MooseEnum>(
       "old_strain_input_window_high_failure",
       error_upper_limit_behavior,
       "What to do if old strain is outside the upper global window of applicability.");
 
   MooseEnum extrapolated_lower_limit_behavior(
       "ERROR EXCEPTION WARN IGNORE DONOTHING USELIMIT EXTRAPOLATE", "EXTRAPOLATE");
   params.addParam<MooseEnum>(
       "stress_input_window_low_failure",
       extrapolated_lower_limit_behavior,
       "What to do if stress is outside the lower global window of applicability.");
   params.addParam<MooseEnum>(
       "stress_input_window_high_failure",
       error_upper_limit_behavior,
       "What to do if stress is outside the upper global window of applicability.");
   params.addParam<MooseEnum>(
       "temperature_input_window_low_failure",
       extrapolated_lower_limit_behavior,
       "What to do if temperature is outside the lower global window of applicability.");
   params.addParam<MooseEnum>(
       "temperature_input_window_high_failure",
       error_upper_limit_behavior,
       "What to do if temperature is outside the upper global window of applicability.");
   params.addParam<MooseEnum>(
       "environment_input_window_low_failure",
       extrapolated_lower_limit_behavior,
       "What to do if environmental factor is outside the lower global window of applicability.");
   params.addParam<MooseEnum>(
       "environment_input_window_high_failure",
       error_upper_limit_behavior,
       "What to do if environmental factor is outside the upper global window of applicability.");
 
   params.addRequiredRangeCheckedParam<Real>(
       "initial_cell_dislocation_density",
       "initial_cell_dislocation_density >= 0.0",
       "Initial density of cell (glissile) dislocations (1/m^2)");
   params.addRangeCheckedParam<Real>(
       "cell_dislocations_normal_distribution_width",
       0.0,
       "cell_dislocations_normal_distribution_width >= 0.0",
       "Width of the normal distribution to assign to the initial cell dislocation value. This is "
       "given as a fraction of the initial_cell_dislocation_density.");
   params.addRangeCheckedParam<Real>(
       "max_relative_cell_dislocation_increment",
       0.5,
       "max_relative_cell_dislocation_increment > 0.0",
       "Maximum increment of density of cell (glissile) dislocations.");
 
   params.addRequiredRangeCheckedParam<Real>(
       "initial_wall_dislocation_density",
       "initial_wall_dislocation_density >= 0.0",
       "wall (locked) dislocation density initial value (1/m^2).");
   params.addRangeCheckedParam<Real>(
       "wall_dislocations_normal_distribution_width",
       0.0,
       "wall_dislocations_normal_distribution_width >= 0.0",
       "Width of the normal distribution to assign to the initial wall dislocation value. This is "
       "given as a fraction of the initial_wall_dislocation_density.");
   params.addRangeCheckedParam<Real>(
       "max_relative_wall_dislocation_increment",
       0.5,
       "max_relative_wall_dislocation_increment > 0.0",
       "Maximum increment of wall (locked) dislocation density initial value (1/m^2).");
 
   params.addParam<bool>("verbose", false, "Flag to output verbose information.");
 
   params.addParam<FunctionName>(
       "cell_dislocation_density_forcing_function",
       "Optional forcing function for wall dislocation. If provided, the wall dislocation "
       "density will be reset to the function value at the beginning of the timestep. Used for "
       "testing purposes only.");
   params.addParam<FunctionName>(
       "wall_dislocation_density_forcing_function",
       "Optional forcing function for wall dislocation. If provided, the wall dislocation "
       "density will be reset to the function value at the beginning of the timestep. Used for "
       "testing purposes only.");
   params.addParam<FunctionName>(
       "old_creep_strain_forcing_function",
       "Optional forcing function for the creep strain from the previous timestep. If provided, "
       "the old creep strain will be reset to the function value at the beginning of the "
       "timestep. Used for testing purposes only.");
   params.addParam<FunctionName>(
       "effective_stress_forcing_function",
       "Optional forcing function for the effective stress. If provided, the effective stress will "
       "be reset to the function value at the beginning of the timestep. Used for testing purposes "
       "only.");
 
   params.addParam<unsigned int>("seed", 0, "Random number generator seed");
   params.addParam<std::string>("stress_unit", "Pa", "unit of stress");
 
   // use std::string here to avoid automatic absolute path expansion
   params.addParam<DataFileName>("model", "LaRomance model JSON datafile");
   params.addParam<FileName>("export_model", "Write LaRomance model to JSON datafile");
 
   params.addParamNamesToGroup(
       "cell_dislocation_density_forcing_function wall_dislocation_density_forcing_function "
       "old_creep_strain_forcing_function effective_stress_forcing_function seed stress_unit",
       "Advanced");
 
   return params;
 }

Member Data Documentation

◆ _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.

◆ _cell_dislocation_increment

template<bool is_ad>

GenericReal<is_ad> LAROMANCEStressUpdateBaseTempl< is_ad >::_cell_dislocation_increment

protected

Container for cell dislocation increment.

Definition at line 482 of file LAROMANCEStressUpdateBase.h.

◆ _cell_dislocations

template<bool is_ad>

GenericMaterialProperty<Real, is_ad>& LAROMANCEStressUpdateBaseTempl< is_ad >::_cell_dislocations

protected

Material properties for cell (glissile) dislocation densities (1/m^2)

Definition at line 471 of file LAROMANCEStressUpdateBase.h.

◆ _cell_dislocations_old

template<bool is_ad>

const MaterialProperty<Real>& LAROMANCEStressUpdateBaseTempl< is_ad >::_cell_dislocations_old

protected

Definition at line 472 of file LAROMANCEStressUpdateBase.h.

◆ _cell_dislocations_step

template<bool is_ad>

GenericMaterialProperty<Real, is_ad>& LAROMANCEStressUpdateBaseTempl< is_ad >::_cell_dislocations_step

protected

Definition at line 630 of file LAROMANCEStressUpdateBase.h.

◆ _cell_function

template<bool is_ad>

const Function* const LAROMANCEStressUpdateBaseTempl< is_ad >::_cell_function

protected

Optional cell dislocation forcing function.

Definition at line 479 of file LAROMANCEStressUpdateBase.h.

◆ _cell_input_index

template<bool is_ad>

const unsigned int LAROMANCEStressUpdateBaseTempl< is_ad >::_cell_input_index

protected

Index corresponding to the position for the dislocations with in the cell in the input vector.

Definition at line 502 of file LAROMANCEStressUpdateBase.h.

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

◆ _cell_output_index

template<bool is_ad>

const unsigned int LAROMANCEStressUpdateBaseTempl< is_ad >::_cell_output_index

protected

Index corresponding to the position for cell dislocations increment in the output vector.

Definition at line 520 of file LAROMANCEStressUpdateBase.h.

◆ _cell_rate

template<bool is_ad>

GenericMaterialProperty<Real, is_ad>& LAROMANCEStressUpdateBaseTempl< is_ad >::_cell_rate

protected

Cell dislocations rate of change.

Definition at line 575 of file LAROMANCEStressUpdateBase.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().

◆ _coefs

template<bool is_ad>

std::vector<std::vector<std::vector<std::vector<Real> > > > LAROMANCEStressUpdateBaseTempl< is_ad >::_coefs

protected

Coefficients used with Legendre polynomials defined by the ROM data set for each partition.

Definition at line 562 of file LAROMANCEStressUpdateBase.h.

◆ _creep_rate

template<bool is_ad>

GenericMaterialProperty<Real, is_ad>& LAROMANCEStressUpdateBaseTempl< is_ad >::_creep_rate

protected

Creep rate material property.

Definition at line 572 of file LAROMANCEStressUpdateBase.h.

◆ _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.

◆ _creep_strain_old_forcing_function

template<bool is_ad>

const Function* const LAROMANCEStressUpdateBaseTempl< is_ad >::_creep_strain_old_forcing_function

protected

Optional old creep strain forcing function.

template<bool is_ad>

GenericReal<is_ad> LAROMANCEStressUpdateBaseTempl< is_ad >::_derivative

protected

Container for derivative of creep increment with respect to strain.

Definition at line 587 of file LAROMANCEStressUpdateBase.h.

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

◆ _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.

◆ _dpartition_weight_dstress

template<bool is_ad>

std::vector<GenericReal<is_ad> > LAROMANCEStressUpdateBaseTempl< is_ad >::_dpartition_weight_dstress

protected

Container for d_parition_weights d_stress.

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().

template<bool is_ad>

MaterialProperty<Real>& LAROMANCEStressUpdateBaseTempl< is_ad >::_extrapolation

protected

Material property to hold smootherstep applied in order to extrapolate.

Definition at line 581 of file LAROMANCEStressUpdateBase.h.

◆ _global_limits

template<bool is_ad>

std::vector<std::pair<Real, Real> > LAROMANCEStressUpdateBaseTempl< is_ad >::_global_limits

protected

Container for global limits.

Definition at line 605 of file LAROMANCEStressUpdateBase.h.

◆ _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.

◆ _index_name

template<bool is_ad>

template<bool is_ad>

unsigned int LAROMANCEStressUpdateBaseTempl< is_ad >::_num_partitions

protected

Number of partitions.

Definition at line 532 of file LAROMANCEStressUpdateBase.h.

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

◆ _num_tiles

template<bool is_ad>

std::vector<unsigned int> LAROMANCEStressUpdateBaseTempl< is_ad >::_num_tiles

protected

Number of ROM tiles per partition.

Definition at line 535 of file LAROMANCEStressUpdateBase.h.

◆ _number_of_substeps

template<bool is_ad>

MaterialProperty<Real>& LAROMANCEStressUpdateBaseTempl< is_ad >::_number_of_substeps

protected

Material property capturing number of substeps for output purposes (defaults to one if substepping isn't used)

Definition at line 637 of file LAROMANCEStressUpdateBase.h.

◆ _old_input_values

template<bool is_ad>

std::vector<Real> LAROMANCEStressUpdateBaseTempl< is_ad >::_old_input_values

protected

Container for old input values.

Definition at line 593 of file LAROMANCEStressUpdateBase.h.

◆ _old_strain_input_index

template<bool is_ad>

const unsigned int LAROMANCEStressUpdateBaseTempl< is_ad >::_old_strain_input_index

protected

Index corresponding to the position for the old strain in the input vector.

Definition at line 511 of file LAROMANCEStressUpdateBase.h.

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

◆ _partition_weights

template<bool is_ad>

std::vector<GenericReal<is_ad> > LAROMANCEStressUpdateBaseTempl< is_ad >::_partition_weights

protected

Container for weights for each tile as computed for all input values beside stress.

Definition at line 614 of file LAROMANCEStressUpdateBase.h.

◆ _plastic_strain_increment

template<bool is_ad>

RankTwoTensor LAROMANCEStressUpdateBaseTempl< is_ad >::_plastic_strain_increment

protected

Total plastic strain increment in step (summing substep contributions)

Definition at line 634 of file LAROMANCEStressUpdateBase.h.

◆ _polynomial_inputs

template<bool is_ad>

std::vector<std::vector<std::vector<std::vector<GenericReal<is_ad> > > > > LAROMANCEStressUpdateBaseTempl< is_ad >::_polynomial_inputs

protected

Container for ROM polynomial inputs.

Definition at line 599 of file LAROMANCEStressUpdateBase.h.

◆ _precomputed_vals

template<bool is_ad>

std::vector<std::vector<std::vector<GenericReal<is_ad> > > > LAROMANCEStressUpdateBaseTempl< is_ad >::_precomputed_vals

protected

Container for ROM precomputed values.

Definition at line 602 of file LAROMANCEStressUpdateBase.h.

◆ _q_point

template<bool is_ad>

const MooseArray< Point > & Material::_q_point

template<bool is_ad>

const unsigned int LAROMANCEStressUpdateBaseTempl< is_ad >::_stress_input_index

protected

Index corresponding to the position for the stress in the input vector.

Definition at line 508 of file LAROMANCEStressUpdateBase.h.

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

◆ _stress_ucf

template<bool is_ad>

Real LAROMANCEStressUpdateBaseTempl< is_ad >::_stress_ucf

protected

Unit conversion factors required to convert from the specified unit to MPa.

Definition at line 626 of file LAROMANCEStressUpdateBase.h.

◆ _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.

template<bool is_ad>

std::vector<std::vector<unsigned int> > LAROMANCEStressUpdateBaseTempl< is_ad >::_tiling

protected

Container for tiling orientations.

Definition at line 620 of file LAROMANCEStressUpdateBase.h.

◆ _transform

template<bool is_ad>

std::vector<std::vector<std::vector<std::vector<ROMInputTransform> > > > LAROMANCEStressUpdateBaseTempl< is_ad >::_transform

protected

Transform rules defined by the ROM data set for each partition.

Definition at line 550 of file LAROMANCEStressUpdateBase.h.

◆ _transform_coefs

template<bool is_ad>

std::vector<std::vector<std::vector<std::vector<Real> > > > LAROMANCEStressUpdateBaseTempl< is_ad >::_transform_coefs

protected

Transform coefficients defined by the ROM data set for each partition.

Definition at line 553 of file LAROMANCEStressUpdateBase.h.

◆ _transformed_normalization_limits

template<bool is_ad>

std::vector<std::vector<std::vector<std::vector<std::vector<Real> > > > > LAROMANCEStressUpdateBaseTempl< is_ad >::_transformed_normalization_limits

protected

Limits transformed from readabile input to ROM readable limits for normalization.

Definition at line 566 of file LAROMANCEStressUpdateBase.h.

◆ _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(), ADMultiplePowerLawCreepStressUpdate::substeppingCapabilityEnabled(), and RadialReturnStressUpdateTempl< is_ad >::substeppingCapabilityRequested().

◆ _verbose

template<bool is_ad>

const bool LAROMANCEStressUpdateBaseTempl< is_ad >::_verbose

protected

Flag to output verbose infromation.

Definition at line 468 of file LAROMANCEStressUpdateBase.h.

◆ _wall_dislocation_increment

template<bool is_ad>

GenericReal<is_ad> LAROMANCEStressUpdateBaseTempl< is_ad >::_wall_dislocation_increment

protected

Container for wall dislocation increment.

Definition at line 499 of file LAROMANCEStressUpdateBase.h.

◆ _wall_dislocations

template<bool is_ad>

GenericMaterialProperty<Real, is_ad>& LAROMANCEStressUpdateBaseTempl< is_ad >::_wall_dislocations

protected

Material properties for wall (locked) dislocation densities (1/m^2)

Definition at line 485 of file LAROMANCEStressUpdateBase.h.

template<bool is_ad>

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Types

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

template<bool is_ad> class LAROMANCEStressUpdateBaseTempl< is_ad >

Member Typedef Documentation

◆ GR2

◆ GR4

◆ GSR2

◆ GSR4

Member Enumeration Documentation

◆ ROMInputTransform

◆ SubsteppingType

◆ WindowFailure

Constructor & Destructor Documentation

◆ LAROMANCEStressUpdateBaseTempl()

Member Function Documentation

◆ areTilesNotIdentical()

◆ buildPolynomials()

◆ calculateNumberSubsteps() [1/2]

◆ calculateNumberSubsteps() [2/2]

◆ checkInTile()

◆ checkJSONKey()

◆ computeDerivative()

◆ computePartitionWeights()

◆ computePolynomial()

◆ computeReferenceResidual()

◆ computeResidual()

◆ computeResidualAndDerivative()

◆ computeROM()

◆ computeStrainEnergyRateDensity()

◆ computeStressDerivative() [1/2]

◆ computeStressDerivative() [2/2]

◆ computeStressFinalize()

◆ computeStressInitialize()

◆ computeTangentOperator() [1/3]

◆ computeTangentOperator() [2/3]

◆ computeTangentOperator() [3/3]

◆ computeTileWeight()

◆ computeTimeStepLimit()

◆ computeValues()

◆ converged()

◆ convertOutput()

◆ convertValue()

◆ coupledGenericValue() [1/3]

◆ coupledGenericValue() [2/3]

◆ coupledGenericValue() [3/3]

◆ effectiveInelasticStrainIncrement()

◆ exportJSON()

◆ getCoefs()

◆ getInputLimits()

◆ getMakeFrameHelper()

◆ getNormalizationLimits()

◆ getStrainCutoff()

◆ getTangentCalculationMethod()

◆ getTilings()

◆ getTransform()

◆ getTransformCoefs()

◆ getTransformedLimits()

◆ initialGuess()

◆ initialSetup()

◆ initQpStatefulProperties()

◆ isIsotropic()

◆ iterationFinalize()

◆ maximumPermissibleValue()

◆ minimumPermissibleValue()

◆ normalizeInput()

◆ outputIterationStep()

◆ outputIterationSummary()

◆ precomputeROM()

◆ precomputeValues()

◆ preStep()

◆ propagateQpStatefulProperties()

◆ propagateQpStatefulPropertiesRadialReturn()

template<bool is_ad>
class LAROMANCEStressUpdateBaseTempl< is_ad >