ADComputeSmearedCrackingStress computes the stress for a finite strain material with smeared cracking. More...

#include <ADComputeSmearedCrackingStress.h>

Inheritance diagram for ADComputeSmearedCrackingStress:

Public Types
enum	ConstantTypeEnum { ConstantTypeEnum::NONE, ConstantTypeEnum::ELEMENT, ConstantTypeEnum::SUBDOMAIN }

enum	TEST_TYPE

typedef DataFileName	DataFileParameterType

Public Member Functions
	ADComputeSmearedCrackingStress (const InputParameters &parameters)

virtual void	initialSetup () override

virtual void	initQpStatefulProperties () override

virtual void	computeQpStress () override

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

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

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
	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	CrackedElasticityType { CrackedElasticityType::DIAGONAL, CrackedElasticityType::FULL }
	Enum defining the method used to adjust the elasticity tensor for cracking. More...

using	ADR2 = Moose::GenericType< R2, true >

using	ADR4 = Moose::GenericType< R4, true >

enum	QP_Data_Type

Protected Member Functions
void	updateLocalElasticityTensor ()
	Update the local elasticity tensor (_local_elasticity_tensor) due to the effects of cracking. More...

virtual void	updateCrackingStateAndStress ()
	Update all cracking-related state variables and the stress tensor due to cracking in all directions. More...

virtual unsigned int	getNumKnownCrackDirs () const
	Get the number of known crack directions. More...

void	computeCrackStrainAndOrientation (ADRealVectorValue &strain_in_crack_dir)
	Compute the crack strain in the crack coordinate system. More...

void	updateStressTensorForCracking (ADRankTwoTensor &tensor, const ADRealVectorValue &sigma)
	Updates the full stress tensor to account for the effect of cracking using the provided stresses in the crack directions. More...

bool	previouslyCracked ()
	Check to see whether there was cracking in any diretion in the previous time step. More...

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

virtual void	finiteStrainRotation ()
	Rotate _elastic_strain, _stress, and _inelastic_strain to the new configuration. More...

virtual void	updateQpState (ADRankTwoTensor &elastic_strain_increment, ADRankTwoTensor &combined_inelastic_strain_increment)
	Given the _strain_increment[_qp], iterate over all of the user-specified recompute materials in order to find an admissible stress (which is placed into _stress[_qp]) and set of inelastic strains. More...

virtual void	updateQpStateSingleModel (unsigned model_number, ADRankTwoTensor &elastic_strain_increment, ADRankTwoTensor &combined_inelastic_strain_increment)
	An optimised version of updateQpState that gets used when the number of plastic models is unity, or when we're cycling through models Given the _strain_increment[_qp], find an admissible stress (which is put into _stress[_qp]) and inelastic strain. More...

virtual void	computeAdmissibleState (unsigned model_number, ADRankTwoTensor &elastic_strain_increment, ADRankTwoTensor &inelastic_strain_increment)
	Given a trial stress (_stress[_qp]) and a strain increment (elastic_strain_increment) let the model_number model produce an admissible stress (gets placed back in _stress[_qp]), and decompose the strain increment into an elastic part (gets placed back into elastic_strain_increment) and an inelastic part (inelastic_strain_increment). More...

virtual void	computeQpProperties () override

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

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

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)

bool	hasGuaranteedMaterialProperty (const MaterialPropertyName &prop, Guarantee guarantee)

Protected Attributes
enum ADComputeSmearedCrackingStress::CrackedElasticityType	_cracked_elasticity_type

ADMaterialProperty< RealVectorValue > &	_crack_flags
	Vector of values going from 1 to 0 as crack damage accumulates. More...

ADRankFourTensor	_local_elasticity_tensor
	Variables used by multiple methods within the calculation for a single material point. More...

std::vector< ADSmearedCrackSofteningBase * >	_softening_models
	The user-supplied list of softening models to be used in the 3 crack directions. More...

std::vector< ADReal >	_local_elastic_vector
	Vector helper to update local elasticity tensor. More...

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

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

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

const unsigned	_num_models
	number of plastic models More...

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

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

MaterialProperty< Real > &	_material_timestep_limit

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

bool	_is_elasticity_tensor_guaranteed_isotropic
	is the elasticity tensor guaranteed to be isotropic? More...

DamageBaseTempl< true > *	_damage_model
	Pointer to the damage model. More...

RankTwoTensor	_undamaged_stress_old

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

const ADMaterialProperty< R4 > &	_elasticity_tensor
	Elasticity tensor material property. More...

const ADMaterialProperty< R2 > &	_strain_increment

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

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

const ADMaterialProperty< RankTwoTensor > &	_rotation_increment

const MaterialProperty< R2 > &	_stress_old
	The old stress tensor. More...

const MaterialProperty< R2 > &	_elastic_strain_old
	The old elastic strain is used to calculate the old stress in the case of variable elasticity tensors. More...

	usingComputeStressBaseMembers

const std::string	_base_name
	Base name of the material system. More...

const ADMaterialProperty< R2 > &	_mechanical_strain

ADMaterialProperty< R2 > &	_stress
	The stress tensor to be calculated. More...

ADMaterialProperty< R2 > &	_elastic_strain

std::vector< const MaterialProperty< R2 > * >	_extra_stresses
	Extra stress tensors. More...

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

	CURR

	PREV

bool	_bnd

bool	_neighbor

const MooseArray< Point > &	_q_point

const QBase *const &	_qrule

const MooseArray< Real > &	_JxW

const Elem *const &	_current_elem

const SubdomainID &	_current_subdomain_id

const unsigned int &	_current_side

const ConstantTypeEnum	_constant_option

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

THREAD_ID	_tid

Assembly &	_assembly

unsigned int	_qp

const MooseArray< Real > &	_coord

const MooseArray< Point > &	_normals

MooseMesh &	_mesh

const Moose::CoordinateSystemType &	_coord_sys

std::set< std::string >	_requested_props

std::set< std::string >	_supplied_props

std::set< unsigned int >	_supplied_prop_ids

std::unordered_set< unsigned int >	_active_prop_ids

const bool	_compute

std::unordered_map< unsigned int, unsigned int >	_props_to_min_states

std::vector< unsigned int >	_displacements

bool	_has_stateful_property

bool	_overrides_init_stateful_props

const FaceInfo *	_face_info

const bool &	_enabled

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

const MaterialData *	_blk_material_data

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

FEProblemBase &	_sc_fe_problem

const THREAD_ID	_sc_tid

const Real &	_real_zero

const VariableValue &	_scalar_zero

const Point &	_point_zero

const InputParameters &	_ti_params

FEProblemBase &	_ti_feproblem

bool	_is_implicit

Real &	_t

const Real &	_t_old

int &	_t_step

Real &	_dt

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


const ADVariableValue &	_cracking_stress
	Input parameters for smeared crack models. More...

std::vector< unsigned int >	_prescribed_crack_directions
	User-prescribed cracking directions. More...

const unsigned int	_max_cracks
	Maximum number of cracks permitted at a material point. More...

const ADReal	_cracking_neg_fraction
	Defines transition to changed stiffness during unloading. More...

const ADReal	_shear_retention_factor
	Controls the amount of shear retained. More...

const ADReal	_max_stress_correction
	Controls the maximum amount that the damaged elastic stress is corrected to folow the release model during a time step. More...


ADMaterialProperty< RealVectorValue > &	_crack_damage

const MaterialProperty< RealVectorValue > &	_crack_damage_old


ADMaterialProperty< RankTwoTensor > &	_crack_rotation

const MaterialProperty< RankTwoTensor > &	_crack_rotation_old


ADMaterialProperty< RealVectorValue > &	_crack_initiation_strain

const MaterialProperty< RealVectorValue > &	_crack_initiation_strain_old


ADMaterialProperty< RealVectorValue > &	_crack_max_strain

const MaterialProperty< RealVectorValue > &	_crack_max_strain_old


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

const Real	_relative_tolerance

const Real	_absolute_tolerance

const bool	_internal_solve_full_iteration_history

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Detailed Description

ADComputeSmearedCrackingStress computes the stress for a finite strain material with smeared cracking.

Definition at line 21 of file ADComputeSmearedCrackingStress.h.

Member Typedef Documentation

◆ ADR2

template<typename R2 , typename R4 >

using ADComputeFiniteStrainElasticStressTempl< R2, R4 >::ADR2 = Moose::GenericType<R2, true>

protectedinherited

Definition at line 32 of file ADComputeFiniteStrainElasticStress.h.

◆ ADR4

template<typename R2 , typename R4 >

using ADComputeFiniteStrainElasticStressTempl< R2, R4 >::ADR4 = Moose::GenericType<R4, true>

protectedinherited

Definition at line 33 of file ADComputeFiniteStrainElasticStress.h.

Member Enumeration Documentation

◆ CrackedElasticityType

enum ADComputeSmearedCrackingStress::CrackedElasticityType

strongprotected

Enum defining the method used to adjust the elasticity tensor for cracking.

Enumerator
DIAGONAL
FULL

Definition at line 99 of file ADComputeSmearedCrackingStress.h.

99 { DIAGONAL, FULL } _cracked_elasticity_type;

ADComputeSmearedCrackingStress::_cracked_elasticity_type

enum ADComputeSmearedCrackingStress::CrackedElasticityType _cracked_elasticity_type

TangentCalculationMethod::FULL

DIAGONAL

Constructor & Destructor Documentation

◆ ADComputeSmearedCrackingStress()

ADComputeSmearedCrackingStress::ADComputeSmearedCrackingStress ( const InputParameters & parameters )

Definition at line 70 of file ADComputeSmearedCrackingStress.C.

   : ADComputeMultipleInelasticStress(parameters),
     _cracking_stress(adCoupledValue("cracking_stress")),
     _max_cracks(getParam<unsigned int>("max_cracks")),
     _cracking_neg_fraction(getParam<Real>("cracking_neg_fraction")),
     _shear_retention_factor(getParam<Real>("shear_retention_factor")),
     _max_stress_correction(getParam<Real>("max_stress_correction")),
     _cracked_elasticity_type(
         getParam<MooseEnum>("cracked_elasticity_type").getEnum<CrackedElasticityType>()),
     _crack_damage(declareADProperty<RealVectorValue>(_base_name + "crack_damage")),
     _crack_damage_old(getMaterialPropertyOld<RealVectorValue>(_base_name + "crack_damage")),
     _crack_flags(declareADProperty<RealVectorValue>(_base_name + "crack_flags")),
     _crack_rotation(declareADProperty<RankTwoTensor>(_base_name + "crack_rotation")),
     _crack_rotation_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "crack_rotation")),
     _crack_initiation_strain(
         declareADProperty<RealVectorValue>(_base_name + "crack_initiation_strain")),
     _crack_initiation_strain_old(
         getMaterialPropertyOld<RealVectorValue>(_base_name + "crack_initiation_strain")),
     _crack_max_strain(declareADProperty<RealVectorValue>(_base_name + "crack_max_strain")),
     _crack_max_strain_old(getMaterialPropertyOld<RealVectorValue>(_base_name + "crack_max_strain"))
 {
   MultiMooseEnum prescribed_crack_directions =
       getParam<MultiMooseEnum>("prescribed_crack_directions");
   if (prescribed_crack_directions.size() > 0)
   {
     if (prescribed_crack_directions.size() > 3)
       mooseError("A maximum of three crack directions may be specified");
     for (unsigned int i = 0; i < prescribed_crack_directions.size(); ++i)
     {
       for (unsigned int j = 0; j < i; ++j)
         if (prescribed_crack_directions[i] == prescribed_crack_directions[j])
           mooseError("Entries in 'prescribed_crack_directions' cannot be repeated");
       _prescribed_crack_directions.push_back(
           static_cast<unsigned int>(prescribed_crack_directions.get(i)));
     }
 
     // Fill in the last remaining direction if 2 are specified
     if (_prescribed_crack_directions.size() == 2)
     {
       std::set<unsigned int> available_dirs = {0, 1, 2};
       for (auto dir : _prescribed_crack_directions)
         if (available_dirs.erase(dir) != 1)
           mooseError("Invalid prescribed crack direction:" + Moose::stringify(dir));
       if (available_dirs.size() != 1)
         mooseError("Error in finding remaining available crack direction");
       _prescribed_crack_directions.push_back(*available_dirs.begin());
     }
   }
   if (!isParamSetByUser("cracked_elasticity_type"))
     paramWarning(
         "cracked_elasticity_type",
         "Defaulting to the legacy option of 'DIAGONAL', but the 'FULL' option is preferred");
 
   _local_elastic_vector.resize(9);
 }

Member Function Documentation

◆ computeAdmissibleState()

void ADComputeMultipleInelasticStress::computeAdmissibleState	(	unsigned	model_number,
		ADRankTwoTensor &	elastic_strain_increment,
		ADRankTwoTensor &	inelastic_strain_increment
	)

protectedvirtualinherited

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

Parameters

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

Definition at line 389 of file ADComputeMultipleInelasticStress.C.

Referenced by ADComputeMultipleInelasticStress::updateQpState(), and ADComputeMultipleInelasticStress::updateQpStateSingleModel().

 {
   // Properly update material properties (necessary if substepping is employed).
   _models[model_number]->resetIncrementalMaterialProperties();
 
   if (_damage_model)
     _models[model_number]->updateState(elastic_strain_increment,
                                        inelastic_strain_increment,
                                        _rotation_increment[_qp],
                                        _stress[_qp],
                                        _undamaged_stress_old,
                                        _elasticity_tensor[_qp],
                                        _elastic_strain_old[_qp]);
   else if (_models[model_number]->substeppingCapabilityEnabled() &&
            (_is_elasticity_tensor_guaranteed_isotropic || !_perform_finite_strain_rotations))
   {
     _models[model_number]->updateStateSubstep(elastic_strain_increment,
                                               inelastic_strain_increment,
                                               _rotation_increment[_qp],
                                               _stress[_qp],
                                               _stress_old[_qp],
                                               _elasticity_tensor[_qp],
                                               _elastic_strain_old[_qp]);
   }
   else
     _models[model_number]->updateState(elastic_strain_increment,
                                        inelastic_strain_increment,
                                        _rotation_increment[_qp],
                                        _stress[_qp],
                                        _stress_old[_qp],
                                        _elasticity_tensor[_qp],
                                        _elastic_strain_old[_qp]);
 }

◆ computeCrackStrainAndOrientation()

void ADComputeSmearedCrackingStress::computeCrackStrainAndOrientation ( ADRealVectorValue & strain_in_crack_dir )

protected

Compute the crack strain in the crack coordinate system.

Also computes the crack orientations, and stores in _crack_rotation.

Parameters

strain_in_crack_dir Computed strains in crack directions

Definition at line 484 of file ADComputeSmearedCrackingStress.C.

Referenced by updateCrackingStateAndStress().

 {
   // The rotation tensor is ordered such that directions for pre-existing cracks appear first
   // in the list of columns.  For example, if there is one existing crack, its direction is in the
   // first column in the rotation tensor.
   const unsigned int num_known_dirs = getNumKnownCrackDirs();
 
   if (num_known_dirs == 0)
   {
     std::vector<ADReal> eigval(3, 0.0);
     ADRankTwoTensor eigvec;
 
     _elastic_strain[_qp].symmetricEigenvaluesEigenvectors(eigval, eigvec);
 
     // If the elastic strain is beyond the cracking strain, save the eigen vectors as
     // the rotation tensor. Reverse their order so that the third principal strain
     // (most tensile) will correspond to the first crack.
     _crack_rotation[_qp].fillColumn(0, eigvec.column(2));
     _crack_rotation[_qp].fillColumn(1, eigvec.column(1));
     _crack_rotation[_qp].fillColumn(2, eigvec.column(0));
 
     strain_in_crack_dir(0) = eigval[2];
     strain_in_crack_dir(1) = eigval[1];
     strain_in_crack_dir(2) = eigval[0];
   }
   else if (num_known_dirs == 1)
   {
     // This is easily the most complicated case.
     // 1.  Rotate the elastic strain to the orientation associated with the known
     //     crack.
     // 2.  Extract the lower 2x2 block into a separate tensor.
     // 3.  Run the eigen solver on the result.
     // 4.  Update the rotation tensor to reflect the effect of the 2 eigenvectors.
 
     // 1.
     const ADRankTwoTensor & R = _crack_rotation[_qp];
     RankTwoTensor ePrime(raw_value(_elastic_strain[_qp]));
     ePrime.rotate(raw_value(R.transpose())); // elastic strain in crack coordinates
 
     // 2.
     ColumnMajorMatrix e2x2(2, 2);
     e2x2(0, 0) = ePrime(1, 1);
     e2x2(1, 0) = ePrime(2, 1);
     e2x2(0, 1) = ePrime(1, 2);
     e2x2(1, 1) = ePrime(2, 2);
 
     // 3.
     ColumnMajorMatrix e_val2x1(2, 1);
     ColumnMajorMatrix e_vec2x2(2, 2);
     e2x2.eigen(e_val2x1, e_vec2x2);
 
     // 4.
     ADRankTwoTensor eigvec(
         1.0, 0.0, 0.0, 0.0, e_vec2x2(0, 1), e_vec2x2(1, 1), 0.0, e_vec2x2(0, 0), e_vec2x2(1, 0));
 
     _crack_rotation[_qp] = _crack_rotation_old[_qp] * eigvec; // Roe implementation
 
     strain_in_crack_dir(0) = ePrime(0, 0);
     strain_in_crack_dir(1) = e_val2x1(1, 0);
     strain_in_crack_dir(2) = e_val2x1(0, 0);
   }
   else if (num_known_dirs == 2 || num_known_dirs == 3)
   {
     // Rotate to cracked orientation and pick off the strains in the rotated
     // coordinate directions.
     const ADRankTwoTensor & R = _crack_rotation[_qp];
     ADRankTwoTensor ePrime(_elastic_strain[_qp]);
     ePrime.rotate(R.transpose()); // elastic strain in crack coordinates
 
     strain_in_crack_dir(0) = ePrime(0, 0);
     strain_in_crack_dir(1) = ePrime(1, 1);
     strain_in_crack_dir(2) = ePrime(2, 2);
   }
   else
     mooseError("Invalid number of known crack directions");
 }

◆ computeQpProperties()

template<typename R2 >

void ADComputeStressBaseTempl< R2 >::computeQpProperties ( )

overrideprotectedvirtualinherited

Reimplemented from Material.

Reimplemented in ADComputeMultiplePorousInelasticStress.

Definition at line 59 of file ADComputeStressBase.C.

Referenced by ADComputeMultiplePorousInelasticStress::computeQpProperties().

 {
   computeQpStress();
 
   // Add in extra stress
   for (MooseIndex(_extra_stresses) i = 0; i < _extra_stresses.size(); ++i)
     _stress[_qp] += (*_extra_stresses[i])[_qp];
 }

◆ computeQpStress()

void ADComputeSmearedCrackingStress::computeQpStress ( )

overridevirtual

Reimplemented from ADComputeMultipleInelasticStress.

Definition at line 219 of file ADComputeSmearedCrackingStress.C.

 {
 
   if (!previouslyCracked())
     computeQpStressIntermediateConfiguration();
   else
   {
     _elastic_strain[_qp] = _elastic_strain_old[_qp] + _strain_increment[_qp];
 
     // Propagate behavior from the (now inactive) inelastic models
     _inelastic_strain[_qp] = _inelastic_strain_old[_qp];
     for (auto model : _models)
     {
       model->setQp(_qp);
       model->propagateQpStatefulProperties();
     }
 
     // Since the other inelastic models are inactive, they will not limit the time step
     _material_timestep_limit[_qp] = std::numeric_limits<Real>::max();
 
     // update _local_elasticity_tensor based on cracking state in previous time step
     updateLocalElasticityTensor();
 
     // Calculate stress in intermediate configuration
     _stress[_qp] = _local_elasticity_tensor * _elastic_strain[_qp];
   }
 
   // compute crack status and adjust stress
   updateCrackingStateAndStress();
 
   if (_perform_finite_strain_rotations)
   {
     finiteStrainRotation();
     _crack_rotation[_qp] = _rotation_increment[_qp] * _crack_rotation[_qp];
   }
 }

◆ computeQpStressIntermediateConfiguration()

void ADComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration ( )

protectedvirtualinherited

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

Definition at line 163 of file ADComputeMultipleInelasticStress.C.

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

 {
   ADRankTwoTensor elastic_strain_increment;
   ADRankTwoTensor combined_inelastic_strain_increment;
 
   if (_num_models == 0)
   {
     _elastic_strain[_qp] = _elastic_strain_old[_qp] + _strain_increment[_qp];
 
     // If the elasticity tensor values have changed and the tensor is isotropic,
     // use the old strain to calculate the old stress
     if (_is_elasticity_tensor_guaranteed_isotropic || !_perform_finite_strain_rotations)
       _stress[_qp] = _elasticity_tensor[_qp] * (_elastic_strain_old[_qp] + _strain_increment[_qp]);
     else
     {
       if (_damage_model)
         paramError(
             "damage_model",
             "Damage models cannot be used with inelastic models and elastic anisotropic behavior");
 
       ADRankFourTensor elasticity_tensor_rotated = _elasticity_tensor[_qp];
       elasticity_tensor_rotated.rotate(_rotation_total_old[_qp]);
 
       _stress[_qp] =
           elasticity_tensor_rotated * (_elastic_strain_old[_qp] + _strain_increment[_qp]);
 
       // Update current total rotation matrix to be used in next step
       _rotation_total[_qp] = _rotation_increment[_qp] * _rotation_total_old[_qp];
     }
   }
   else
   {
     if (_num_models == 1 || _cycle_models)
       updateQpStateSingleModel((_t_step - 1) % _num_models,
                                elastic_strain_increment,
                                combined_inelastic_strain_increment);
     else
       updateQpState(elastic_strain_increment, combined_inelastic_strain_increment);
 
     _elastic_strain[_qp] = _elastic_strain_old[_qp] + elastic_strain_increment;
     _inelastic_strain[_qp] = _inelastic_strain_old[_qp] + combined_inelastic_strain_increment;
   }
 }

◆ finiteStrainRotation()

void ADComputeMultipleInelasticStress::finiteStrainRotation ( )

protectedvirtualinherited

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

Definition at line 208 of file ADComputeMultipleInelasticStress.C.

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

 {
   _elastic_strain[_qp] =
       _rotation_increment[_qp] * _elastic_strain[_qp] * _rotation_increment[_qp].transpose();
   _stress[_qp] = _rotation_increment[_qp] * _stress[_qp] * _rotation_increment[_qp].transpose();
   _inelastic_strain[_qp] =
       _rotation_increment[_qp] * _inelastic_strain[_qp] * _rotation_increment[_qp].transpose();
 }

◆ getNumKnownCrackDirs()

unsigned int ADComputeSmearedCrackingStress::getNumKnownCrackDirs ( ) const

protectedvirtual

Get the number of known crack directions.

This includes cracks in prescribed directions (even if not yet active) and active cracks in other directions.

Returns: number of known crack directions

Definition at line 563 of file ADComputeSmearedCrackingStress.C.

Referenced by computeCrackStrainAndOrientation().

 {
   unsigned int num_known_dirs = 0;
   for (unsigned int i = 0; i < 3; ++i)
   {
     if (_crack_damage_old[_qp](i) > 0.0 || _prescribed_crack_directions.size() >= i + 1)
       ++num_known_dirs;
   }
   return num_known_dirs;
 }

◆ hasGuaranteedMaterialProperty()

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

protectedinherited

Definition at line 28 of file GuaranteeConsumer.C.

Referenced by ComputeFiniteStrainElasticStress::computeQpStress(), initialSetup(), ComputeSmearedCrackingStress::initialSetup(), ComputeLinearElasticPFFractureStress::initialSetup(), ComputeSimoHughesJ2PlasticityStress::initialSetup(), CriticalTimeStep::initialSetup(), ComputeStVenantKirchhoffStress::initialSetup(), ComputeMultipleInelasticStressBase::initialSetup(), and ADComputeMultipleInelasticStress::initialSetup().

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

◆ initialSetup()

void ADComputeSmearedCrackingStress::initialSetup ( )

overridevirtual

Reimplemented from ADComputeMultipleInelasticStress.

Definition at line 190 of file ADComputeSmearedCrackingStress.C.

 {
   ADComputeMultipleInelasticStress::initialSetup();
 
   if (!hasGuaranteedMaterialProperty(_elasticity_tensor_name, Guarantee::ISOTROPIC))
     mooseError("ADComputeSmearedCrackingStress requires that the elasticity tensor be "
                "guaranteed isotropic");
 
   std::vector<MaterialName> soft_matls = getParam<std::vector<MaterialName>>("softening_models");
   for (auto soft_matl : soft_matls)
   {
     ADSmearedCrackSofteningBase * scsb =
         dynamic_cast<ADSmearedCrackSofteningBase *>(&getMaterialByName(soft_matl));
     if (scsb)
       _softening_models.push_back(scsb);
     else
       paramError("softening_models", "Model " + soft_matl + " is not a softening model");
   }
   if (_softening_models.size() == 1)
   {
     // Reuse the same model in all 3 directions
     _softening_models.push_back(_softening_models[0]);
     _softening_models.push_back(_softening_models[0]);
   }
   else if (_softening_models.size() != 3)
     paramError("softening_models", "Either 1 or 3 softening models must be specified");
 }

◆ initQpStatefulProperties()

void ADComputeSmearedCrackingStress::initQpStatefulProperties ( )

overridevirtual

Reimplemented from ADComputeMultipleInelasticStress.

Definition at line 127 of file ADComputeSmearedCrackingStress.C.

 {
   ADComputeMultipleInelasticStress::initQpStatefulProperties();
 
   _crack_damage[_qp] = 0.0;
 
   _crack_initiation_strain[_qp] = 0.0;
   _crack_max_strain[_qp](0) = 0.0;
 
   switch (_prescribed_crack_directions.size())
   {
     case 0:
     {
       _crack_rotation[_qp] = ADRankTwoTensor::Identity();
       break;
     }
     case 1:
     {
       _crack_rotation[_qp].zero();
       switch (_prescribed_crack_directions[0])
       {
         case 0:
         {
           _crack_rotation[_qp](0, 0) = 1.0;
           _crack_rotation[_qp](1, 1) = 1.0;
           _crack_rotation[_qp](2, 2) = 1.0;
           break;
         }
         case 1:
         {
           _crack_rotation[_qp](1, 0) = 1.0;
           _crack_rotation[_qp](0, 1) = 1.0;
           _crack_rotation[_qp](2, 2) = 1.0;
           break;
         }
         case 2:
         {
           _crack_rotation[_qp](2, 0) = 1.0;
           _crack_rotation[_qp](0, 1) = 1.0;
           _crack_rotation[_qp](1, 2) = 1.0;
           break;
         }
       }
       break;
     }
     case 2:
     {
       mooseError("Number of prescribed crack directions cannot be 2");
       break;
     }
     case 3:
     {
       for (unsigned int i = 0; i < _prescribed_crack_directions.size(); ++i)
       {
         ADRealVectorValue crack_dir_vec;
         crack_dir_vec(_prescribed_crack_directions[i]) = 1.0;
         _crack_rotation[_qp].fillColumn(i, crack_dir_vec);
       }
     }
   }
 }

◆ previouslyCracked()

bool ADComputeSmearedCrackingStress::previouslyCracked ( )

protected

Check to see whether there was cracking in any diretion in the previous time step.

Returns: true if cracked, false if not cracked

Definition at line 601 of file ADComputeSmearedCrackingStress.C.

Referenced by computeQpStress().

 {
   for (unsigned int i = 0; i < 3; ++i)
     if (_crack_damage_old[_qp](i) > 0.0)
       return true;
   return false;
 }

◆ updateCrackingStateAndStress()

void ADComputeSmearedCrackingStress::updateCrackingStateAndStress ( )

protectedvirtual

Update all cracking-related state variables and the stress tensor due to cracking in all directions.

Definition at line 373 of file ADComputeSmearedCrackingStress.C.

Referenced by computeQpStress().

 {
   const ADReal youngs_modulus =
       ElasticityTensorTools::getIsotropicYoungsModulus(_elasticity_tensor[_qp]);
 
   ADReal cracking_stress = _cracking_stress[_qp];
 
   if (cracking_stress > 0)
   {
     // Initializing crack states
     _crack_rotation[_qp] = _crack_rotation_old[_qp];
 
     for (unsigned i = 0; i < 3; ++i)
     {
       _crack_max_strain[_qp](i) = _crack_max_strain_old[_qp](i);
       _crack_initiation_strain[_qp](i) = _crack_initiation_strain_old[_qp](i);
       _crack_damage[_qp](i) = _crack_damage_old[_qp](i);
     }
 
     // Compute crack orientations: updated _crack_rotation[_qp] based on current strain
     ADRealVectorValue strain_in_crack_dir;
     computeCrackStrainAndOrientation(strain_in_crack_dir);
 
     for (unsigned i = 0; i < 3; ++i)
     {
       if (strain_in_crack_dir(i) > _crack_max_strain[_qp](i))
         _crack_max_strain[_qp](i) = strain_in_crack_dir(i);
     }
 
     // Check for new cracks.
     // Rotate stress to cracked orientation.
     const ADRankTwoTensor & R = _crack_rotation[_qp];
     ADRankTwoTensor sigmaPrime(_stress[_qp]);
     sigmaPrime.rotate(R.transpose()); // stress in crack coordinates
 
     unsigned int num_cracks = 0;
     for (unsigned int i = 0; i < 3; ++i)
     {
       if (_crack_damage_old[_qp](i) > 0.0)
         ++num_cracks;
     }
 
     bool cracked(false);
     ADRealVectorValue sigma;
     mooseAssert(_softening_models.size() == 3, "Must have 3 softening models");
     for (unsigned int i = 0; i < 3; ++i)
     {
       sigma(i) = sigmaPrime(i, i);
 
       ADReal stiffness_ratio = 1.0 - _crack_damage[_qp](i);
 
       const bool pre_existing_crack = (_crack_damage_old[_qp](i) > 0.0);
       const bool met_stress_criterion = (sigma(i) > cracking_stress);
       const bool loading_existing_crack = (strain_in_crack_dir(i) >= _crack_max_strain[_qp](i));
       const bool allowed_to_crack = (pre_existing_crack || num_cracks < _max_cracks);
       bool new_crack = false;
 
       cracked |= pre_existing_crack;
 
       // Adjustments for newly created cracks
       if (met_stress_criterion && !pre_existing_crack && allowed_to_crack)
       {
         new_crack = true;
         ++num_cracks;
 
         // Assume Poisson's ratio drops to zero for this direction.  Stiffness is then Young's
         // modulus.
         _crack_initiation_strain[_qp](i) = cracking_stress / youngs_modulus;
 
         if (_crack_max_strain[_qp](i) < _crack_initiation_strain[_qp](i))
           _crack_max_strain[_qp](i) = _crack_initiation_strain[_qp](i);
       }
 
       // Update stress and stiffness ratio according to specified crack release model
       if (new_crack || (pre_existing_crack && loading_existing_crack))
       {
         cracked = true;
         _softening_models[i]->computeCrackingRelease(sigma(i),
                                                      stiffness_ratio,
                                                      strain_in_crack_dir(i),
                                                      _crack_initiation_strain[_qp](i),
                                                      _crack_max_strain[_qp](i),
                                                      cracking_stress,
                                                      youngs_modulus);
         _crack_damage[_qp](i) = 1.0 - stiffness_ratio;
       }
 
       else if (cracked && _cracking_neg_fraction > 0 &&
                _crack_initiation_strain[_qp](i) * _cracking_neg_fraction > strain_in_crack_dir(i) &&
                -_crack_initiation_strain[_qp](i) * _cracking_neg_fraction < strain_in_crack_dir(i))
       {
         const ADReal etr = _cracking_neg_fraction * _crack_initiation_strain[_qp](i);
         const ADReal Eo = cracking_stress / _crack_initiation_strain[_qp](i);
         const ADReal Ec = Eo * (1.0 - _crack_damage_old[_qp](i));
         const ADReal a = (Ec - Eo) / (4.0 * etr);
         const ADReal b = 0.5 * (Ec + Eo);
         const ADReal c = 0.25 * (Ec - Eo) * etr;
         sigma(i) = (a * strain_in_crack_dir(i) + b) * strain_in_crack_dir(i) + c;
       }
     }
 
     if (cracked)
       updateStressTensorForCracking(_stress[_qp], sigma);
   }
 
   _crack_flags[_qp](0) = 1.0 - _crack_damage[_qp](2);
   _crack_flags[_qp](1) = 1.0 - _crack_damage[_qp](1);
   _crack_flags[_qp](2) = 1.0 - _crack_damage[_qp](0);
 }

◆ updateLocalElasticityTensor()

void ADComputeSmearedCrackingStress::updateLocalElasticityTensor ( )

protected

Update the local elasticity tensor (_local_elasticity_tensor) due to the effects of cracking.

Definition at line 257 of file ADComputeSmearedCrackingStress.C.

Referenced by computeQpStress().

 {
   const ADReal youngs_modulus =
       ElasticityTensorTools::getIsotropicYoungsModulus(_elasticity_tensor[_qp]);
 
   bool cracking_locally_active = false;
 
   const ADReal cracking_stress = _cracking_stress[_qp];
 
   if (cracking_stress > 0)
   {
     ADRealVectorValue stiffness_ratio_local(1.0, 1.0, 1.0);
     const ADRankTwoTensor & R = _crack_rotation_old[_qp];
     ADRankTwoTensor ePrime(_elastic_strain_old[_qp]);
     ePrime.rotate(R.transpose());
 
     for (unsigned int i = 0; i < 3; ++i)
     {
       // Update elasticity tensor based on crack status of the end of last time step
       if (_crack_damage_old[_qp](i) > 0.0)
       {
         if (_cracking_neg_fraction == 0.0 && MooseUtils::absoluteFuzzyLessThan(ePrime(i, i), 0.0))
           stiffness_ratio_local(i) = 1.0;
         else if (_cracking_neg_fraction > 0.0 &&
                  ePrime(i, i) < _crack_initiation_strain_old[_qp](i) * _cracking_neg_fraction &&
                  ePrime(i, i) > -_crack_initiation_strain_old[_qp](i) * _cracking_neg_fraction)
         {
           const ADReal etr = _cracking_neg_fraction * _crack_initiation_strain_old[_qp](i);
           const ADReal Eo = cracking_stress / _crack_initiation_strain_old[_qp](i);
           const ADReal Ec = Eo * (1.0 - _crack_damage_old[_qp](i));
           const ADReal a = (Ec - Eo) / (4 * etr);
           const ADReal b = (Ec + Eo) / 2;
           // Compute the ratio of the current transition stiffness to the original stiffness
           stiffness_ratio_local(i) = (2.0 * a * etr + b) / Eo;
           cracking_locally_active = true;
         }
         else
         {
           stiffness_ratio_local(i) = (1.0 - _crack_damage_old[_qp](i));
           cracking_locally_active = true;
         }
       }
     }
 
     if (cracking_locally_active)
     {
       // Update the elasticity tensor in the crack coordinate system
       if (_cracked_elasticity_type == CrackedElasticityType::DIAGONAL)
       {
         const bool c0_coupled = MooseUtils::absoluteFuzzyEqual(stiffness_ratio_local(0), 1.0);
         const bool c1_coupled = MooseUtils::absoluteFuzzyEqual(stiffness_ratio_local(1), 1.0);
         const bool c2_coupled = MooseUtils::absoluteFuzzyEqual(stiffness_ratio_local(2), 1.0);
 
         const ADReal c01 = (c0_coupled && c1_coupled ? 1.0 : 0.0);
         const ADReal c02 = (c0_coupled && c2_coupled ? 1.0 : 0.0);
         const ADReal c12 = (c1_coupled && c2_coupled ? 1.0 : 0.0);
 
         const ADReal c01_shear_retention =
             (c0_coupled && c1_coupled ? 1.0 : _shear_retention_factor);
         const ADReal c02_shear_retention =
             (c0_coupled && c2_coupled ? 1.0 : _shear_retention_factor);
         const ADReal c12_shear_retention =
             (c1_coupled && c2_coupled ? 1.0 : _shear_retention_factor);
 
         _local_elastic_vector[0] = (c0_coupled ? _elasticity_tensor[_qp](0, 0, 0, 0)
                                                : stiffness_ratio_local(0) * youngs_modulus);
         _local_elastic_vector[1] = _elasticity_tensor[_qp](0, 0, 1, 1) * c01;
         _local_elastic_vector[2] = _elasticity_tensor[_qp](0, 0, 2, 2) * c02;
         _local_elastic_vector[3] = (c1_coupled ? _elasticity_tensor[_qp](1, 1, 1, 1)
                                                : stiffness_ratio_local(1) * youngs_modulus);
         _local_elastic_vector[4] = _elasticity_tensor[_qp](1, 1, 2, 2) * c12;
         _local_elastic_vector[5] = (c2_coupled ? _elasticity_tensor[_qp](2, 2, 2, 2)
                                                : stiffness_ratio_local(2) * youngs_modulus);
         _local_elastic_vector[6] = _elasticity_tensor[_qp](1, 2, 1, 2) * c12_shear_retention;
         _local_elastic_vector[7] = _elasticity_tensor[_qp](0, 2, 0, 2) * c02_shear_retention;
         _local_elastic_vector[8] = _elasticity_tensor[_qp](0, 1, 0, 1) * c01_shear_retention;
       }
       else // _cracked_elasticity_type == CrackedElasticityType::FULL
       {
         const ADReal & c0 = stiffness_ratio_local(0);
         const ADReal & c1 = stiffness_ratio_local(1);
         const ADReal & c2 = stiffness_ratio_local(2);
 
         const ADReal c01 = c0 * c1;
         const ADReal c02 = c0 * c2;
         const ADReal c12 = c1 * c2;
 
         const ADReal c01_shear_retention = std::max(c01, _shear_retention_factor);
         const ADReal c02_shear_retention = std::max(c02, _shear_retention_factor);
         const ADReal c12_shear_retention = std::max(c12, _shear_retention_factor);
 
         _local_elastic_vector[0] = _elasticity_tensor[_qp](0, 0, 0, 0) * c0;
         _local_elastic_vector[1] = _elasticity_tensor[_qp](0, 0, 1, 1) * c01;
         _local_elastic_vector[2] = _elasticity_tensor[_qp](0, 0, 2, 2) * c02;
         _local_elastic_vector[3] = _elasticity_tensor[_qp](1, 1, 1, 1) * c1;
         _local_elastic_vector[4] = _elasticity_tensor[_qp](1, 1, 2, 2) * c12;
         _local_elastic_vector[5] = _elasticity_tensor[_qp](2, 2, 2, 2) * c2;
         _local_elastic_vector[6] = _elasticity_tensor[_qp](1, 2, 1, 2) * c12_shear_retention;
         _local_elastic_vector[7] = _elasticity_tensor[_qp](0, 2, 0, 2) * c02_shear_retention;
         _local_elastic_vector[8] = _elasticity_tensor[_qp](0, 1, 0, 1) * c01_shear_retention;
       }
 
       // Filling with 9 components is sufficient because these are the only nonzero entries
       // for isotropic or orthotropic materials.
       _local_elasticity_tensor.fillFromInputVector(_local_elastic_vector,
                                                    ADRankFourTensor::symmetric9);
 
       // Rotate the modified elasticity tensor back into global coordinates
       _local_elasticity_tensor.rotate(R);
     }
   }
   if (!cracking_locally_active)
     _local_elasticity_tensor = _elasticity_tensor[_qp];
 }

◆ updateQpState()

void ADComputeMultipleInelasticStress::updateQpState	(	ADRankTwoTensor &	elastic_strain_increment,
		ADRankTwoTensor &	combined_inelastic_strain_increment
	)

protectedvirtualinherited

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

Parameters

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

Definition at line 218 of file ADComputeMultipleInelasticStress.C.

Referenced by ADComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration().

 {
   if (_internal_solve_full_iteration_history == true)
   {
     _console << std::endl
              << "iteration output for ADComputeMultipleInelasticStress solve:"
              << " time=" << _t << " int_pt=" << _qp << std::endl;
   }
   Real l2norm_delta_stress;
   Real first_l2norm_delta_stress = 1.0;
   unsigned int counter = 0;
 
   std::vector<ADRankTwoTensor> inelastic_strain_increment;
   inelastic_strain_increment.resize(_num_models);
 
   for (unsigned i_rmm = 0; i_rmm < _models.size(); ++i_rmm)
     inelastic_strain_increment[i_rmm].zero();
 
   ADRankTwoTensor stress_max, stress_min;
 
   do
   {
     for (unsigned i_rmm = 0; i_rmm < _num_models; ++i_rmm)
     {
       _models[i_rmm]->setQp(_qp);
 
       // initially assume the strain is completely elastic
       elastic_strain_increment = _strain_increment[_qp];
       // and subtract off all inelastic strain increments calculated so far
       // except the one that we're about to calculate
       for (unsigned j_rmm = 0; j_rmm < _num_models; ++j_rmm)
         if (i_rmm != j_rmm)
           elastic_strain_increment -= inelastic_strain_increment[j_rmm];
 
       // form the trial stress, with the check for changed elasticity constants
       if (_is_elasticity_tensor_guaranteed_isotropic || !_perform_finite_strain_rotations)
         _stress[_qp] =
             _elasticity_tensor[_qp] * (_elastic_strain_old[_qp] + elastic_strain_increment);
       else
       {
         if (_damage_model)
           paramError("damage_model",
                      "Damage models cannot be used with inelastic models and elastic anisotropic "
                      "behavior");
 
         ADRankFourTensor elasticity_tensor_rotated = _elasticity_tensor[_qp];
         elasticity_tensor_rotated.rotate(_rotation_total_old[_qp]);
 
         _stress[_qp] =
             elasticity_tensor_rotated * (_elastic_strain_old[_qp] + elastic_strain_increment);
 
         // Update current total rotation matrix to be used in next step
         _rotation_total[_qp] = _rotation_increment[_qp] * _rotation_total_old[_qp];
       }
       // given a trial stress (_stress[_qp]) and a strain increment (elastic_strain_increment)
       // let the i^th model produce an admissible stress (as _stress[_qp]), and decompose
       // the strain increment into an elastic part (elastic_strain_increment) and an
       // inelastic part (inelastic_strain_increment[i_rmm])
       computeAdmissibleState(i_rmm, elastic_strain_increment, inelastic_strain_increment[i_rmm]);
 
       if (i_rmm == 0)
       {
         stress_max = _stress[_qp];
         stress_min = _stress[_qp];
       }
       else
       {
         for (const auto i : make_range(Moose::dim))
         {
           for (const auto j : make_range(Moose::dim))
           {
             if (_stress[_qp](i, j) > stress_max(i, j))
               stress_max(i, j) = _stress[_qp](i, j);
             else if (stress_min(i, j) > _stress[_qp](i, j))
               stress_min(i, j) = _stress[_qp](i, j);
           }
         }
       }
     }
 
     // now check convergence in the stress:
     // once the change in stress is within tolerance after each recompute material
     // consider the stress to be converged
     l2norm_delta_stress = MetaPhysicL::raw_value((stress_max - stress_min).L2norm());
     if (counter == 0 && l2norm_delta_stress > 0.0)
       first_l2norm_delta_stress = l2norm_delta_stress;
 
     if (_internal_solve_full_iteration_history == true)
     {
       _console << "stress iteration number = " << counter << "\n"
                << " relative l2 norm delta stress = "
                << (0 == first_l2norm_delta_stress ? 0
                                                   : l2norm_delta_stress / first_l2norm_delta_stress)
                << "\n"
                << " stress convergence relative tolerance = " << _relative_tolerance << "\n"
                << " absolute l2 norm delta stress = " << l2norm_delta_stress << "\n"
                << " stress convergence absolute tolerance = " << _absolute_tolerance << std::endl;
     }
     ++counter;
 
   } while (counter < _max_iterations && l2norm_delta_stress > _absolute_tolerance &&
            (l2norm_delta_stress / first_l2norm_delta_stress) > _relative_tolerance &&
            _num_models != 1);
 
   if (counter == _max_iterations && l2norm_delta_stress > _absolute_tolerance &&
       (l2norm_delta_stress / first_l2norm_delta_stress) > _relative_tolerance)
     mooseException(
         "In ", _name, ": Max stress iteration hit during ADComputeMultipleInelasticStress solve!");
 
   combined_inelastic_strain_increment.zero();
   for (unsigned i_rmm = 0; i_rmm < _num_models; ++i_rmm)
     combined_inelastic_strain_increment +=
         _inelastic_weights[i_rmm] * inelastic_strain_increment[i_rmm];
 
   _material_timestep_limit[_qp] = 0.0;
   for (unsigned i_rmm = 0; i_rmm < _num_models; ++i_rmm)
     _material_timestep_limit[_qp] += 1.0 / _models[i_rmm]->computeTimeStepLimit();
 
   if (MooseUtils::absoluteFuzzyEqual(_material_timestep_limit[_qp], 0.0))
     _material_timestep_limit[_qp] = std::numeric_limits<Real>::max();
   else
     _material_timestep_limit[_qp] = 1.0 / _material_timestep_limit[_qp];
 }

◆ updateQpStateSingleModel()

void ADComputeMultipleInelasticStress::updateQpStateSingleModel	(	unsigned	model_number,
		ADRankTwoTensor &	elastic_strain_increment,
		ADRankTwoTensor &	combined_inelastic_strain_increment
	)

protectedvirtualinherited

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

Parameters

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

Definition at line 345 of file ADComputeMultipleInelasticStress.C.

Referenced by ADComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration().

 {
   for (auto model : _models)
     model->setQp(_qp);
 
   elastic_strain_increment = _strain_increment[_qp];
 
   // If the elasticity tensor values have changed and the tensor is isotropic,
   // use the old strain to calculate the old stress
   if (_is_elasticity_tensor_guaranteed_isotropic || !_perform_finite_strain_rotations)
     _stress[_qp] = _elasticity_tensor[_qp] * (_elastic_strain_old[_qp] + elastic_strain_increment);
   else
   {
     if (_damage_model)
       paramError(
           "damage_model",
           "Damage models cannot be used with inelastic models and elastic anisotropic behavior");
 
     ADRankFourTensor elasticity_tensor_rotated = _elasticity_tensor[_qp];
     elasticity_tensor_rotated.rotate(_rotation_total_old[_qp]);
 
     _stress[_qp] =
         elasticity_tensor_rotated * (_elastic_strain_old[_qp] + elastic_strain_increment);
 
     // Update current total rotation matrix to be used in next step
     _rotation_total[_qp] = _rotation_increment[_qp] * _rotation_total_old[_qp];
   }
 
   computeAdmissibleState(
       model_number, elastic_strain_increment, combined_inelastic_strain_increment);
 
   _material_timestep_limit[_qp] = _models[0]->computeTimeStepLimit();
 
   /* propagate internal variables, etc, to this timestep for those inelastic models where
    * "updateState" is not called */
   for (unsigned i_rmm = 0; i_rmm < _num_models; ++i_rmm)
     if (i_rmm != model_number)
       _models[i_rmm]->propagateQpStatefulProperties();
 }

◆ updateStressTensorForCracking()

void ADComputeSmearedCrackingStress::updateStressTensorForCracking	(	ADRankTwoTensor &	tensor,
		const ADRealVectorValue &	sigma
	)

protected

Updates the full stress tensor to account for the effect of cracking using the provided stresses in the crack directions.

The tensor is rotated into the crack coordinates, modified, and rotated back.

Parameters

tensor	Stress tensor to be updated
sigma	Vector of stresses in crack directions

Definition at line 575 of file ADComputeSmearedCrackingStress.C.

Referenced by updateCrackingStateAndStress().

 {
   // Get transformation matrix
   const ADRankTwoTensor & R = _crack_rotation[_qp];
   // Rotate to crack frame
   tensor.rotate(R.transpose());
 
   // Reset stress if cracked
   for (unsigned int i = 0; i < 3; ++i)
     if (_crack_damage[_qp](i) > 0.0)
     {
       const ADReal stress_correction_ratio = (tensor(i, i) - sigma(i)) / tensor(i, i);
       if (stress_correction_ratio > _max_stress_correction)
         tensor(i, i) *= (1.0 - _max_stress_correction);
       else if (stress_correction_ratio < -_max_stress_correction)
         tensor(i, i) *= (1.0 + _max_stress_correction);
       else
         tensor(i, i) = sigma(i);
     }
 
   // Rotate back to global frame
   tensor.rotate(R);
 }

◆ validParams()

InputParameters ADComputeSmearedCrackingStress::validParams ( )

static

Definition at line 18 of file ADComputeSmearedCrackingStress.C.

 {
   InputParameters params = ADComputeMultipleInelasticStress::validParams();
   params.addClassDescription(
       "Compute stress using a fixed smeared cracking model. Uses automatic differentiation");
   params.addRequiredParam<std::vector<MaterialName>>(
       "softening_models",
       "The material objects used to compute softening behavior for loading a crack."
       "Either 1 or 3 models must be specified. If a single model is specified, it is"
       "used for all directions. If 3 models are specified, they will be used for the"
       "3 crack directions in sequence");
   params.addRequiredCoupledVar(
       "cracking_stress",
       "The stress threshold beyond which cracking occurs. Negative values prevent cracking.");
   MultiMooseEnum direction("x y z");
   params.addParam<MultiMooseEnum>(
       "prescribed_crack_directions", direction, "Prescribed directions of first cracks");
   params.addParam<unsigned int>(
       "max_cracks", 3, "The maximum number of cracks allowed at a material point.");
   params.addRangeCheckedParam<Real>("cracking_neg_fraction",
                                     0,
                                     "cracking_neg_fraction <= 1 & cracking_neg_fraction >= 0",
                                     "The fraction of the cracking strain at which "
                                     "a transition begins during decreasing "
                                     "strain to the original stiffness.");
   params.addParam<Real>(
       "max_stress_correction",
       1.0,
       "Maximum permitted correction to the predicted stress as a ratio of the "
       "stress change to the predicted stress from the previous step's damage level. "
       "Values less than 1 will improve robustness, but not be as accurate.");
 
   params.addRangeCheckedParam<Real>(
       "shear_retention_factor",
       0.0,
       "shear_retention_factor>=0 & shear_retention_factor<=1.0",
       "Fraction of original shear stiffness to be retained after cracking");
   params.set<std::vector<MaterialName>>("inelastic_models") = {};
 
   MooseEnum crackedElasticityType("DIAGONAL FULL", "DIAGONAL");
   crackedElasticityType.addDocumentation(
       "DIAGONAL", "Zero out terms coupling with directions orthogonal to a crack (legacy)");
   crackedElasticityType.addDocumentation(
       "FULL", "Consistently scale all entries based on damage (recommended)");
   params.addParam<MooseEnum>(
       "cracked_elasticity_type",
       crackedElasticityType,
       "Method to modify the local elasticity tensor to account for cracking");
 
   return params;
 }

Member Data Documentation

◆ _absolute_tolerance

const Real ADComputeMultipleInelasticStress::_absolute_tolerance

protectedinherited

Definition at line 105 of file ADComputeMultipleInelasticStress.h.

Referenced by ADComputeMultipleInelasticStress::updateQpState().

◆ _base_name

template<typename R2 >

const std::string ADComputeStressBaseTempl< R2 >::_base_name

protectedinherited

Base name of the material system.

Definition at line 43 of file ADComputeStressBase.h.

◆ _crack_damage

ADMaterialProperty<RealVectorValue>& ADComputeSmearedCrackingStress::_crack_damage

protected

Definition at line 102 of file ADComputeSmearedCrackingStress.h.

Referenced by initQpStatefulProperties(), updateCrackingStateAndStress(), and updateStressTensorForCracking().

◆ _crack_damage_old

const MaterialProperty<RealVectorValue>& ADComputeSmearedCrackingStress::_crack_damage_old

protected

Definition at line 103 of file ADComputeSmearedCrackingStress.h.

Referenced by getNumKnownCrackDirs(), previouslyCracked(), updateCrackingStateAndStress(), and updateLocalElasticityTensor().

◆ _crack_flags

ADMaterialProperty<RealVectorValue>& ADComputeSmearedCrackingStress::_crack_flags

protected

Vector of values going from 1 to 0 as crack damage accumulates.

Legacy property for backward compatibility – remove in the future.

Definition at line 108 of file ADComputeSmearedCrackingStress.h.

Referenced by updateCrackingStateAndStress().

◆ _crack_initiation_strain

ADMaterialProperty<RealVectorValue>& ADComputeSmearedCrackingStress::_crack_initiation_strain

protected

Definition at line 116 of file ADComputeSmearedCrackingStress.h.

Referenced by initQpStatefulProperties(), and updateCrackingStateAndStress().

◆ _crack_initiation_strain_old

const MaterialProperty<RealVectorValue>& ADComputeSmearedCrackingStress::_crack_initiation_strain_old

protected

Definition at line 117 of file ADComputeSmearedCrackingStress.h.

Referenced by updateCrackingStateAndStress(), and updateLocalElasticityTensor().

◆ _crack_max_strain

ADMaterialProperty<RealVectorValue>& ADComputeSmearedCrackingStress::_crack_max_strain

protected

Definition at line 121 of file ADComputeSmearedCrackingStress.h.

Referenced by initQpStatefulProperties(), and updateCrackingStateAndStress().

◆ _crack_max_strain_old

const MaterialProperty<RealVectorValue>& ADComputeSmearedCrackingStress::_crack_max_strain_old

protected

Definition at line 122 of file ADComputeSmearedCrackingStress.h.

Referenced by updateCrackingStateAndStress().

◆ _crack_rotation

ADMaterialProperty<RankTwoTensor>& ADComputeSmearedCrackingStress::_crack_rotation

protected

Definition at line 111 of file ADComputeSmearedCrackingStress.h.

Referenced by computeCrackStrainAndOrientation(), computeQpStress(), initQpStatefulProperties(), updateCrackingStateAndStress(), and updateStressTensorForCracking().

◆ _crack_rotation_old

const MaterialProperty<RankTwoTensor>& ADComputeSmearedCrackingStress::_crack_rotation_old

protected

Definition at line 112 of file ADComputeSmearedCrackingStress.h.

Referenced by computeCrackStrainAndOrientation(), updateCrackingStateAndStress(), and updateLocalElasticityTensor().

◆ _cracked_elasticity_type

enum ADComputeSmearedCrackingStress::CrackedElasticityType ADComputeSmearedCrackingStress::_cracked_elasticity_type

protected

Referenced by updateLocalElasticityTensor().

◆ _cracking_neg_fraction

const ADReal ADComputeSmearedCrackingStress::_cracking_neg_fraction

protected

Defines transition to changed stiffness during unloading.

Definition at line 88 of file ADComputeSmearedCrackingStress.h.

Referenced by updateCrackingStateAndStress(), and updateLocalElasticityTensor().

◆ _cracking_stress

const ADVariableValue& ADComputeSmearedCrackingStress::_cracking_stress

protected

Input parameters for smeared crack models.

Threshold at which cracking initiates if tensile stress exceeds it

Definition at line 79 of file ADComputeSmearedCrackingStress.h.

Referenced by updateCrackingStateAndStress(), and updateLocalElasticityTensor().

◆ _cycle_models

const bool ADComputeMultipleInelasticStress::_cycle_models

protectedinherited

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

Definition at line 125 of file ADComputeMultipleInelasticStress.h.

Referenced by ADComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration().

◆ _damage_model

DamageBaseTempl<true>* ADComputeMultipleInelasticStress::_damage_model

protectedinherited

Pointer to the damage model.

Definition at line 142 of file ADComputeMultipleInelasticStress.h.

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

◆ _elastic_strain

template<typename R2 >

ADMaterialProperty<R2>& ADComputeStressBaseTempl< R2 >::_elastic_strain

protectedinherited

Definition at line 49 of file ADComputeStressBase.h.

Referenced by computeCrackStrainAndOrientation(), computeQpStress(), ADComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration(), and ADComputeMultipleInelasticStress::finiteStrainRotation().

◆ _elastic_strain_old

template<typename R2 , typename R4 >

const MaterialProperty<R2>& ADComputeFiniteStrainElasticStressTempl< R2, R4 >::_elastic_strain_old

protectedinherited

The old elastic strain is used to calculate the old stress in the case of variable elasticity tensors.

Definition at line 56 of file ADComputeFiniteStrainElasticStress.h.

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

◆ _elasticity_tensor

template<typename R2 , typename R4 >

const ADMaterialProperty<R4>& ADComputeFiniteStrainElasticStressTempl< R2, R4 >::_elasticity_tensor

protectedinherited

Elasticity tensor material property.

Definition at line 40 of file ADComputeFiniteStrainElasticStress.h.

Referenced by ADComputeMultipleInelasticStress::computeAdmissibleState(), ADComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration(), updateCrackingStateAndStress(), updateLocalElasticityTensor(), ADComputeMultipleInelasticStress::updateQpState(), and ADComputeMultipleInelasticStress::updateQpStateSingleModel().

◆ _elasticity_tensor_name

template<typename R2 , typename R4 >

const std::string ADComputeFiniteStrainElasticStressTempl< R2, R4 >::_elasticity_tensor_name

protectedinherited

Name of the elasticity tensor material property.

Definition at line 38 of file ADComputeFiniteStrainElasticStress.h.

Referenced by initialSetup(), and ADComputeMultipleInelasticStress::initialSetup().

◆ _extra_stresses

template<typename R2 >

std::vector<const MaterialProperty<R2> *> ADComputeStressBaseTempl< R2 >::_extra_stresses

protectedinherited

Extra stress tensors.

Definition at line 52 of file ADComputeStressBase.h.

Referenced by ADComputeStressBaseTempl< R2 >::ADComputeStressBaseTempl().

◆ _inelastic_strain

ADMaterialProperty<RankTwoTensor>& ADComputeMultipleInelasticStress::_inelastic_strain

protectedinherited

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

Definition at line 113 of file ADComputeMultipleInelasticStress.h.

Referenced by ADComputeMultiplePorousInelasticStress::computeQpProperties(), computeQpStress(), ADComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration(), ADComputeMultipleInelasticStress::finiteStrainRotation(), and ADComputeMultipleInelasticStress::initQpStatefulProperties().

◆ _inelastic_strain_old

const MaterialProperty<RankTwoTensor>& ADComputeMultipleInelasticStress::_inelastic_strain_old

protectedinherited

old value of inelastic strain

Definition at line 116 of file ADComputeMultipleInelasticStress.h.

Referenced by ADComputeMultiplePorousInelasticStress::computeQpProperties(), computeQpStress(), and ADComputeMultipleInelasticStress::computeQpStressIntermediateConfiguration().

◆ _inelastic_weights

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

protectedinherited

_inelastic_strain = sum_i (_inelastic_weights_i * inelastic_strain_from_model_i)

Definition at line 122 of file ADComputeMultipleInelasticStress.h.

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

◆ _initial_stress_fcn

template<typename R2 >

std::vector<const Function *> ADComputeStressBaseTempl< R2 >::_initial_stress_fcn

protectedinherited

initial stress components

Definition at line 55 of file ADComputeStressBase.h.

◆ _internal_solve_full_iteration_history

const bool ADComputeMultipleInelasticStress::_internal_solve_full_iteration_history

protectedinherited

Definition at line 106 of file ADComputeMultipleInelasticStress.h.

Referenced by ADComputeMultipleInelasticStress::updateQpState().

◆ _is_elasticity_tensor_guaranteed_isotropic

bool ADComputeMultipleInelasticStress::_is_elasticity_tensor_guaranteed_isotropic

protectedinherited

is the elasticity tensor guaranteed to be isotropic?

Definition at line 139 of file ADComputeMultipleInelasticStress.h.

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

◆ _local_elastic_vector

std::vector<ADReal> ADComputeSmearedCrackingStress::_local_elastic_vector

protected

Vector helper to update local elasticity tensor.

Definition at line 132 of file ADComputeSmearedCrackingStress.h.

Referenced by ADComputeSmearedCrackingStress(), and updateLocalElasticityTensor().

◆ _local_elasticity_tensor

ADRankFourTensor ADComputeSmearedCrackingStress::_local_elasticity_tensor

protected

Variables used by multiple methods within the calculation for a single material point.

Definition at line 126 of file ADComputeSmearedCrackingStress.h.

Referenced by computeQpStress(), and updateLocalElasticityTensor().

◆ _material_timestep_limit

MaterialProperty<Real>& ADComputeMultipleInelasticStress::_material_timestep_limit

protectedinherited

Definition at line 127 of file ADComputeMultipleInelasticStress.h.

Referenced by computeQpStress(), ADComputeMultipleInelasticStress::computeQpStress(), ADComputeMultipleInelasticStress::updateQpState(), and ADComputeMultipleInelasticStress::updateQpStateSingleModel().

◆ _max_cracks

const unsigned int ADComputeSmearedCrackingStress::_max_cracks

protected

Maximum number of cracks permitted at a material point.

Definition at line 85 of file ADComputeSmearedCrackingStress.h.

Referenced by updateCrackingStateAndStress().

◆ _max_iterations

const unsigned int ADComputeMultipleInelasticStress::_max_iterations

protectedinherited

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

Definition at line 103 of file ADComputeMultipleInelasticStress.h.

Referenced by ADComputeMultipleInelasticStress::updateQpState().

◆ _max_stress_correction

const ADReal ADComputeSmearedCrackingStress::_max_stress_correction

protected

Controls the maximum amount that the damaged elastic stress is corrected to folow the release model during a time step.

Definition at line 95 of file ADComputeSmearedCrackingStress.h.

Referenced by updateStressTensorForCracking().

◆ _mechanical_strain

template<typename R2 >

const ADMaterialProperty<R2>& ADComputeStressBaseTempl< R2 >::_mechanical_strain

protectedinherited

Definition at line 45 of file ADComputeStressBase.h.

Referenced by ADComputeStrainIncrementBasedStress::computeQpStress().

◆ _models

std::vector<ADStressUpdateBase *> ADComputeMultipleInelasticStress::_models

protectedinherited

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

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

Definition at line 136 of file ADComputeMultipleInelasticStress.h.

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