ComputeMultipleCrystalPlasticityStress (used together with CrystalPlasticityStressUpdateBase) uses the multiplicative decomposition of the deformation gradient and solves the PK2 stress residual equation at the intermediate configuration to evolve the material state. More...

#include <ComputeMultipleCrystalPlasticityStress.h>

Inheritance diagram for ComputeMultipleCrystalPlasticityStress:

[legend]

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

enum	TEST_TYPE

typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

typedef DataFileName	DataFileParameterType

Public Member Functions
	ComputeMultipleCrystalPlasticityStress (const InputParameters &parameters)

virtual void	initialSetup () override

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

void	validateDerivativeMaterialPropertyBase (const std::string &base)

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

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)

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)

std::unordered_map< SubdomainID, std::vector< MaterialBase *> >	buildRequiredMaterials (bool allow_stateful=true)

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 ADVariableValue *	getADDefaultValue (const std::string &var_name) 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 ADVectorVariableValue *	getADDefaultVectorValue (const std::string &var_name) const

const ADVariableGradient &	getADDefaultGradient () const

const ADVectorVariableGradient &	getADDefaultVectorGradient () const

const ADVariableSecond &	getADDefaultSecond () const

const ADVectorVariableCurl &	getADDefaultCurl () const

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)

void	statefulPropertiesAllowed (bool)

bool	getMaterialPropertyCalled () const

const GenericMaterialProperty< T, is_ad > &	getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)

const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

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

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

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	TangentModuliType { TangentModuliType::EXACT, TangentModuliType::NONE }
	Type of tangent moduli calculation. More...

enum	LineSearchMethod { LineSearchMethod::CutHalf, LineSearchMethod::Bisection }
	strain formulation More...

enum	QP_Data_Type

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

virtual void	updateStress (RankTwoTensor &cauchy_stress, RankFourTensor &jacobian_mult)
	Updates the stress (PK2) at a quadrature point by calling constiutive relationship as defined in a child class Solves stress residual equation using Newton - Raphson: Updates slip system resistances iteratively. More...

virtual void	initQpStatefulProperties () override
	initializes the stateful properties such as PK2 stress, resolved shear stress, plastic deformation gradient, slip system resistances, etc. More...

void	calculateResidualAndJacobian ()
	Calls the residual and jacobian functions used in the stress update algorithm. More...

void	preSolveQp ()
	Reset the PK2 stress and the inverse deformation gradient to old values and provide an interface for inheriting classes to reset material properties. More...

void	solveQp ()
	Solve the stress and internal state variables (e.g. More...

void	postSolveQp (RankTwoTensor &stress_new, RankFourTensor &jacobian_mult)
	Save the final stress and internal variable values after the iterative solve. More...

void	solveStateVariables ()
	Solves the internal variables stress as a function of the slip specified by the constitutive model defined in the inheriting class. More...

void	solveStress ()
	solves for stress, updates plastic deformation gradient. More...

void	calculateResidual ()
	Calculate stress residual as the difference between the stored material property PK2 stress and the elastic PK2 stress calculated from the constitutively defined equivalent_slip_increment. More...

void	calculateJacobian ()
	Calculates the jacobian as ${J} = {I} - {C} {d{E}^e}{d{F}^e} {d{F}^e}{d{F}^P^{-1}} {d{F}^P^{-1}}{d{PK2}}$. More...

bool	lineSearchUpdate (const Real &rnorm_prev, const RankTwoTensor &dpk2)
	performs the line search update More...

void	calculateEigenstrainDeformationGrad ()
	Calculates the deformation gradient due to eigenstrain. 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 ()

void	markMatPropRequested (const std::string &)

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

bool	checkVar (const std::string &var_name, unsigned int comp=0, unsigned int comp_bound=0) const

void	validateExecutionerType (const std::string &name, const std::string &fn_name) 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 ()

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

std::vector< T >	coupledVectorHelper (const std::string &var_name, const Func &func) const

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)


void	calcTangentModuli (RankFourTensor &jacobian_mult)
	Calculates the tangent moduli for use as a preconditioner, using the elastic or elastic-plastic option as specified by the user. More...

void	elasticTangentModuli (RankFourTensor &jacobian_mult)

void	elastoPlasticTangentModuli (RankFourTensor &jacobian_mult)

Protected Attributes
const unsigned	_num_models
	number of plastic models More...

std::vector< CrystalPlasticityStressUpdateBase * >	_models
	The user supplied cyrstal plasticity consititutive models. More...

const unsigned	_num_eigenstrains
	number of eigenstrains More...

std::vector< ComputeCrystalPlasticityEigenstrainBase * >	_eigenstrains
	The user supplied cyrstal plasticity eigenstrains. More...

const std::string	_base_name
	optional parameter to define several mechanical systems on the same block, e.g. multiple phases More...

const MaterialProperty< RankFourTensor > &	_elasticity_tensor
	Elasticity tensor as defined by a separate class. More...

Real	_rtol
	Stress residual equation relative tolerance. More...

Real	_abs_tol
	Stress residual equation absolute tolerance. More...

RankTwoTensor	_residual_tensor
	Residual tensor. More...

RankFourTensor	_jacobian
	Jacobian tensor. More...

unsigned int	_maxiter
	Maximum number of iterations for stress update. More...

unsigned int	_maxiterg
	Maximum number of iterations for internal variable update. More...

enum ComputeMultipleCrystalPlasticityStress::TangentModuliType	_tan_mod_type

unsigned int	_max_substep_iter
	Maximum number of substep iterations. More...

Real	_substep_dt
	time step size during substepping More...

bool	_use_line_search
	Flag to activate line serach. More...

Real	_min_line_search_step_size
	Minimum line search step size. More...

Real	_line_search_tolerance
	Line search bisection method tolerance. More...

unsigned int	_line_search_max_iterations
	Line search bisection method maximum iteration number. More...

enum ComputeMultipleCrystalPlasticityStress::LineSearchMethod	_line_search_method

MaterialProperty< RankTwoTensor > &	_total_lagrangian_strain
	Lagrangian total strain measure for the entire crystal. More...

MaterialProperty< RankTwoTensor > &	_updated_rotation
	Tracks the rotation of the crystal during deformation Note: this rotation tensor is not applied to the crystal lattice. More...

const MaterialProperty< RankTwoTensor > &	_crysrot
	Crystal rotation in the original, or reference, configuration as defined by Euler angle arguments in the ComputeElasticityTensor classes. More...

const bool	_print_convergence_message
	Flag to print to console warning messages on stress, constitutive model convergence. More...

bool	_convergence_failed
	Flag to check whether convergence is achieved or if substepping is needed. More...

Real	_dfgrd_scale_factor
	Scales the substepping increment to obtain deformation gradient at a substep iteration. More...

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

MaterialProperty< 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 MaterialProperty< RankTwoTensor > &	_strain_increment
	Strain increment material property. More...

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

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

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

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

MaterialProperty< RankTwoTensor > &	_stress
	Stress material property. More...

MaterialProperty< RankTwoTensor > &	_elastic_strain
	Elastic strain material property. More...

const MaterialProperty< RankTwoTensor > &	_extra_stress
	Extra stress tensor. More...

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

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

	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

std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > >	_default_value

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


MaterialProperty< RankTwoTensor > &	_plastic_deformation_gradient
	Plastic deformation gradient RankTwoTensor for the crystal. More...

const MaterialProperty< RankTwoTensor > &	_plastic_deformation_gradient_old


MaterialProperty< RankTwoTensor > *	_eigenstrain_deformation_gradient
	Generalized eigenstrain deformation gradient RankTwoTensor for the crystal. More...

const MaterialProperty< RankTwoTensor > *	_eigenstrain_deformation_gradient_old


const MaterialProperty< RankTwoTensor > &	_deformation_gradient
	Total deformation gradient RankTwoTensor for the crystal. More...

const MaterialProperty< RankTwoTensor > &	_deformation_gradient_old


MaterialProperty< RankTwoTensor > &	_pk2
	Second Piola-Kirchoff stress measure. More...

const MaterialProperty< RankTwoTensor > &	_pk2_old


RankTwoTensor	_temporary_deformation_gradient
	Helper deformation gradient tensor variables used in iterative solve. More...

RankTwoTensor	_elastic_deformation_gradient

RankTwoTensor	_inverse_plastic_deformation_grad

RankTwoTensor	_inverse_plastic_deformation_grad_old

RankTwoTensor	_inverse_eigenstrain_deformation_grad


RankTwoTensor	_delta_deformation_gradient
	Used for substepping; Uniformly divides the increment in deformation gradient. More...

RankTwoTensor	_temporary_deformation_gradient_old

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Detailed Description

ComputeMultipleCrystalPlasticityStress (used together with CrystalPlasticityStressUpdateBase) uses the multiplicative decomposition of the deformation gradient and solves the PK2 stress residual equation at the intermediate configuration to evolve the material state.

The internal variables are updated using an iterative predictor-corrector algorithm. Backward Euler integration rule is used for the rate equations.

This material that is not called by MOOSE because of the compute=false flag set in the parameter list.

Definition at line 30 of file ComputeMultipleCrystalPlasticityStress.h.

Member Enumeration Documentation

◆ LineSearchMethod

enum ComputeMultipleCrystalPlasticityStress::LineSearchMethod

strongprotected

strain formulation

Enumerator
CutHalf
Bisection

Definition at line 175 of file ComputeMultipleCrystalPlasticityStress.h.

175 { CutHalf, Bisection } _line_search_method;

ComputeMultipleCrystalPlasticityStress::_line_search_method

enum ComputeMultipleCrystalPlasticityStress::LineSearchMethod _line_search_method

◆ TangentModuliType

enum ComputeMultipleCrystalPlasticityStress::TangentModuliType

strongprotected

Type of tangent moduli calculation.

Enumerator
EXACT
NONE

Definition at line 154 of file ComputeMultipleCrystalPlasticityStress.h.

154 { EXACT, NONE } _tan_mod_type;

EXACT

NONE

ComputeMultipleCrystalPlasticityStress::_tan_mod_type

enum ComputeMultipleCrystalPlasticityStress::TangentModuliType _tan_mod_type

Constructor & Destructor Documentation

◆ ComputeMultipleCrystalPlasticityStress()

ComputeMultipleCrystalPlasticityStress::ComputeMultipleCrystalPlasticityStress ( const InputParameters & parameters )

Definition at line 64 of file ComputeMultipleCrystalPlasticityStress.C.

   : ComputeFiniteStrainElasticStress(parameters),
     _num_models(getParam<std::vector<MaterialName>>("crystal_plasticity_models").size()),
     _num_eigenstrains(getParam<std::vector<MaterialName>>("eigenstrain_names").size()),
     _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
     _elasticity_tensor(getMaterialPropertyByName<RankFourTensor>(_base_name + "elasticity_tensor")),
     _rtol(getParam<Real>("rtol")),
     _abs_tol(getParam<Real>("abs_tol")),
     _maxiter(getParam<unsigned int>("maxiter")),
     _maxiterg(getParam<unsigned int>("maxiter_state_variable")),
     _tan_mod_type(getParam<MooseEnum>("tan_mod_type").getEnum<TangentModuliType>()),
     _max_substep_iter(getParam<unsigned int>("maximum_substep_iteration")),
     _use_line_search(getParam<bool>("use_line_search")),
     _min_line_search_step_size(getParam<Real>("min_line_search_step_size")),
     _line_search_tolerance(getParam<Real>("line_search_tol")),
     _line_search_max_iterations(getParam<unsigned int>("line_search_maxiter")),
     _line_search_method(getParam<MooseEnum>("line_search_method").getEnum<LineSearchMethod>()),
     _plastic_deformation_gradient(declareProperty<RankTwoTensor>("plastic_deformation_gradient")),
     _plastic_deformation_gradient_old(
         getMaterialPropertyOld<RankTwoTensor>("plastic_deformation_gradient")),
     _eigenstrain_deformation_gradient(
         _num_eigenstrains ? &declareProperty<RankTwoTensor>("eigenstrain_deformation_gradient")
                           : nullptr),
     _eigenstrain_deformation_gradient_old(
         _num_eigenstrains
             ? &getMaterialPropertyOld<RankTwoTensor>("eigenstrain_deformation_gradient")
             : nullptr),
     _deformation_gradient(getMaterialProperty<RankTwoTensor>(_base_name + "deformation_gradient")),
     _deformation_gradient_old(
         getMaterialPropertyOld<RankTwoTensor>(_base_name + "deformation_gradient")),
     _pk2(declareProperty<RankTwoTensor>("second_piola_kirchhoff_stress")),
     _pk2_old(getMaterialPropertyOld<RankTwoTensor>("second_piola_kirchhoff_stress")),
     _total_lagrangian_strain(
         declareProperty<RankTwoTensor>("total_lagrangian_strain")), // Lagrangian strain
     _updated_rotation(declareProperty<RankTwoTensor>("updated_rotation")),
     _crysrot(getMaterialProperty<RankTwoTensor>(
         _base_name + "crysrot")), // defined in the elasticity tensor classes for crystal plasticity
     _print_convergence_message(getParam<bool>("print_state_variable_convergence_error_messages"))
 {
   _convergence_failed = false;
 }

Member Function Documentation

◆ calcTangentModuli()

void ComputeMultipleCrystalPlasticityStress::calcTangentModuli ( RankFourTensor & jacobian_mult )

protected

Calculates the tangent moduli for use as a preconditioner, using the elastic or elastic-plastic option as specified by the user.

Definition at line 570 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by postSolveQp().

 {
   switch (_tan_mod_type)
   {
     case TangentModuliType::EXACT:
       elastoPlasticTangentModuli(jacobian_mult);
       break;
     default:
       elasticTangentModuli(jacobian_mult);
   }
 }

◆ calculateEigenstrainDeformationGrad()

void ComputeMultipleCrystalPlasticityStress::calculateEigenstrainDeformationGrad ( )

protected

Calculates the deformation gradient due to eigenstrain.

Definition at line 711 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by updateStress().

 {
   _inverse_eigenstrain_deformation_grad.zero();
   _inverse_eigenstrain_deformation_grad.addIa(1.0);
 
   for (unsigned int i = 0; i < _num_eigenstrains; ++i)
   {
     _eigenstrains[i]->setSubstepDt(_substep_dt);
     _eigenstrains[i]->computeQpProperties();
     _inverse_eigenstrain_deformation_grad *= _eigenstrains[i]->getDeformationGradientInverse();
   }
   (*_eigenstrain_deformation_gradient)[_qp] = _inverse_eigenstrain_deformation_grad.inverse();
 }

◆ calculateJacobian()

void ComputeMultipleCrystalPlasticityStress::calculateJacobian ( )

protected

Calculates the jacobian as ${J} = {I} - {C} {d{E}^e}{d{F}^e} {d{F}^e}{d{F}^P^{-1}} {d{F}^P^{-1}}{d{PK2}}$.

Definition at line 535 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by calculateResidualAndJacobian().

 {
   // may not need to cache the dfpinvdpk2 here. need to double check
   RankFourTensor dfedfpinv, deedfe, dfpinvdpk2, dfpinvdpk2_per_model;
 
   RankTwoTensor ffeiginv = _temporary_deformation_gradient * _inverse_eigenstrain_deformation_grad;
 
   for (const auto i : make_range(Moose::dim))
     for (const auto j : make_range(Moose::dim))
       for (const auto k : make_range(Moose::dim))
         dfedfpinv(i, j, k, j) = ffeiginv(i, k);
 
   for (const auto i : make_range(Moose::dim))
     for (const auto j : make_range(Moose::dim))
       for (const auto k : make_range(Moose::dim))
       {
         deedfe(i, j, k, i) = deedfe(i, j, k, i) + _elastic_deformation_gradient(k, j) * 0.5;
         deedfe(i, j, k, j) = deedfe(i, j, k, j) + _elastic_deformation_gradient(k, i) * 0.5;
       }
 
   for (unsigned int i = 0; i < _num_models; ++i)
   {
     _models[i]->calculateTotalPlasticDeformationGradientDerivative(
         dfpinvdpk2_per_model,
         _inverse_plastic_deformation_grad_old,
         _inverse_eigenstrain_deformation_grad,
         _num_eigenstrains);
     dfpinvdpk2 += dfpinvdpk2_per_model;
   }
 
   _jacobian =
       RankFourTensor::IdentityFour() - (_elasticity_tensor[_qp] * deedfe * dfedfpinv * dfpinvdpk2);
 }

◆ calculateResidual()

void ComputeMultipleCrystalPlasticityStress::calculateResidual ( )

protected

Calculate stress residual as the difference between the stored material property PK2 stress and the elastic PK2 stress calculated from the constitutively defined equivalent_slip_increment.

The equivalent_slip_increment is passed in as an input arguement.

Definition at line 491 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by calculateResidualAndJacobian(), and lineSearchUpdate().

 {
   RankTwoTensor ce, elastic_strain, ce_pk2, equivalent_slip_increment_per_model,
       equivalent_slip_increment, pk2_new;
 
   equivalent_slip_increment.zero();
 
   // calculate slip rate in order to compute F^{p-1}
   for (unsigned int i = 0; i < _num_models; ++i)
   {
     equivalent_slip_increment_per_model.zero();
 
     // calculate shear stress with consideration of contribution from other physics
     _models[i]->calculateShearStress(
         _pk2[_qp], _inverse_eigenstrain_deformation_grad, _num_eigenstrains);
 
     _convergence_failed = !_models[i]->calculateSlipRate();
 
     if (_convergence_failed)
       return;
 
     _models[i]->calculateEquivalentSlipIncrement(equivalent_slip_increment_per_model);
     equivalent_slip_increment += equivalent_slip_increment_per_model;
   }
 
   RankTwoTensor residual_equivalent_slip_increment =
       RankTwoTensor::Identity() - equivalent_slip_increment;
   _inverse_plastic_deformation_grad =
       _inverse_plastic_deformation_grad_old * residual_equivalent_slip_increment;
 
   _elastic_deformation_gradient = _temporary_deformation_gradient *
                                   _inverse_eigenstrain_deformation_grad *
                                   _inverse_plastic_deformation_grad;
 
   ce = _elastic_deformation_gradient.transpose() * _elastic_deformation_gradient;
 
   elastic_strain = ce - RankTwoTensor::Identity();
   elastic_strain *= 0.5;
 
   pk2_new = _elasticity_tensor[_qp] * elastic_strain;
   _residual_tensor = _pk2[_qp] - pk2_new;
 }

◆ calculateResidualAndJacobian()

void ComputeMultipleCrystalPlasticityStress::calculateResidualAndJacobian ( )

protected

Calls the residual and jacobian functions used in the stress update algorithm.

Definition at line 482 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by solveStress().

 {
   calculateResidual();
   if (_convergence_failed)
     return;
   calculateJacobian();
 }

◆ computeQpProperties()

void ComputeGeneralStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

Reimplemented from DerivativeMaterialInterface< Material >.

Definition at line 44 of file ComputeGeneralStressBase.C.

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

◆ computeQpStress()

void ComputeMultipleCrystalPlasticityStress::computeQpStress ( )

overrideprotectedvirtual

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

Reimplemented from ComputeFiniteStrainElasticStress.

Definition at line 186 of file ComputeMultipleCrystalPlasticityStress.C.

 {
   for (unsigned int i = 0; i < _num_models; ++i)
   {
     _models[i]->setQp(_qp);
     _models[i]->setMaterialVectorSize();
   }
 
   for (unsigned int i = 0; i < _num_eigenstrains; ++i)
     _eigenstrains[i]->setQp(_qp);
 
   updateStress(_stress[_qp], _Jacobian_mult[_qp]); // This is NOT the exact jacobian
 }

◆ elasticTangentModuli()

void ComputeMultipleCrystalPlasticityStress::elasticTangentModuli ( RankFourTensor & jacobian_mult )

protected

Definition at line 629 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by calcTangentModuli().

 {
   // update jacobian_mult
   jacobian_mult = _elasticity_tensor[_qp];
 }

◆ elastoPlasticTangentModuli()

void ComputeMultipleCrystalPlasticityStress::elastoPlasticTangentModuli ( RankFourTensor & jacobian_mult )

protected

Definition at line 583 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by calcTangentModuli().

 {
   RankFourTensor tan_mod;
   RankTwoTensor pk2fet, fepk2, feiginvfpinv;
   RankFourTensor deedfe, dsigdpk2dfe, dfedf;
 
   // Fill in the matrix stiffness material property
   for (const auto i : make_range(Moose::dim))
     for (const auto j : make_range(Moose::dim))
       for (const auto k : make_range(Moose::dim))
       {
         deedfe(i, j, k, i) = deedfe(i, j, k, i) + _elastic_deformation_gradient(k, j) * 0.5;
         deedfe(i, j, k, j) = deedfe(i, j, k, j) + _elastic_deformation_gradient(k, i) * 0.5;
       }
 
   usingTensorIndices(i_, j_, k_, l_);
   dsigdpk2dfe = _elastic_deformation_gradient.times<i_, k_, j_, l_>(_elastic_deformation_gradient) *
                 _elasticity_tensor[_qp] * deedfe;
 
   pk2fet = _pk2[_qp] * _elastic_deformation_gradient.transpose();
   fepk2 = _elastic_deformation_gradient * _pk2[_qp];
 
   for (const auto i : make_range(Moose::dim))
     for (const auto j : make_range(Moose::dim))
       for (const auto l : make_range(Moose::dim))
       {
         tan_mod(i, j, i, l) += pk2fet(l, j);
         tan_mod(i, j, j, l) += fepk2(i, l);
       }
 
   tan_mod += dsigdpk2dfe;
 
   const auto je = _elastic_deformation_gradient.det();
   if (je > 0.0)
     tan_mod /= je;
 
   feiginvfpinv = _inverse_eigenstrain_deformation_grad * _inverse_plastic_deformation_grad;
   for (const auto i : make_range(Moose::dim))
     for (const auto j : make_range(Moose::dim))
       for (const auto l : make_range(Moose::dim))
         dfedf(i, j, i, l) = feiginvfpinv(l, j);
 
   jacobian_mult = tan_mod * dfedf;
 }

◆ hasGuaranteedMaterialProperty()

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

protectedinherited

Definition at line 28 of file GuaranteeConsumer.C.

Referenced by ComputeFiniteStrainElasticStress::computeQpStress(), ADComputeSmearedCrackingStress::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 ComputeMultipleCrystalPlasticityStress::initialSetup ( )

overridevirtual

Reimplemented from ComputeFiniteStrainElasticStress.

Definition at line 146 of file ComputeMultipleCrystalPlasticityStress.C.

 {
   // get crystal plasticity models
   std::vector<MaterialName> model_names =
       getParam<std::vector<MaterialName>>("crystal_plasticity_models");
 
   for (unsigned int i = 0; i < _num_models; ++i)
   {
     CrystalPlasticityStressUpdateBase * model =
         dynamic_cast<CrystalPlasticityStressUpdateBase *>(&getMaterialByName(model_names[i]));
 
     if (model)
     {
       _models.push_back(model);
       // TODO: check to make sure that the material model is compatible with this class
     }
     else
       mooseError("Model " + model_names[i] +
                  " is not compatible with ComputeMultipleCrystalPlasticityStress");
   }
 
   // get crystal plasticity eigenstrains
   std::vector<MaterialName> eigenstrain_names =
       getParam<std::vector<MaterialName>>("eigenstrain_names");
 
   for (unsigned int i = 0; i < _num_eigenstrains; ++i)
   {
     ComputeCrystalPlasticityEigenstrainBase * eigenstrain =
         dynamic_cast<ComputeCrystalPlasticityEigenstrainBase *>(
             &getMaterialByName(eigenstrain_names[i]));
 
     if (eigenstrain)
       _eigenstrains.push_back(eigenstrain);
     else
       mooseError("Eigenstrain" + eigenstrain_names[i] +
                  " is not compatible with ComputeMultipleCrystalPlasticityStress");
   }
 }

◆ initQpStatefulProperties()

void ComputeMultipleCrystalPlasticityStress::initQpStatefulProperties ( )

overrideprotectedvirtual

initializes the stateful properties such as PK2 stress, resolved shear stress, plastic deformation gradient, slip system resistances, etc.

This class is often overwritten by inherting classes.

Reimplemented from ComputeFiniteStrainElasticStress.

Definition at line 108 of file ComputeMultipleCrystalPlasticityStress.C.

 {
   _plastic_deformation_gradient[_qp].zero();
   _plastic_deformation_gradient[_qp].addIa(1.0);
 
   if (_num_eigenstrains)
   {
     (*_eigenstrain_deformation_gradient)[_qp].zero();
     (*_eigenstrain_deformation_gradient)[_qp].addIa(1.0);
   }
   else
   {
     // set to identity if no eigenstrain is added
     _inverse_eigenstrain_deformation_grad.zero();
     _inverse_eigenstrain_deformation_grad.addIa(1.0);
   }
 
   _pk2[_qp].zero();
 
   _total_lagrangian_strain[_qp].zero();
 
   _updated_rotation[_qp].zero();
   _updated_rotation[_qp].addIa(1.0);
 
   for (unsigned int i = 0; i < _num_models; ++i)
   {
     _models[i]->setQp(_qp);
     _models[i]->initQpStatefulProperties();
   }
 
   for (unsigned int i = 0; i < _num_eigenstrains; ++i)
   {
     _eigenstrains[i]->setQp(_qp);
     _eigenstrains[i]->initQpStatefulProperties();
   }
 }

◆ lineSearchUpdate()

bool ComputeMultipleCrystalPlasticityStress::lineSearchUpdate	(	const Real &	rnorm_prev,
		const RankTwoTensor &	dpk2
	)

protected

performs the line search update

Definition at line 636 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by solveStress().

 {
   if (_line_search_method == LineSearchMethod::CutHalf)
   {
     Real rnorm;
     Real step = 1.0;
 
     do
     {
       _pk2[_qp] = _pk2[_qp] - step * dpk2;
       step /= 2.0;
       _pk2[_qp] = _pk2[_qp] + step * dpk2;
 
       calculateResidual();
       rnorm = _residual_tensor.L2norm();
     } while (rnorm > rnorm_prev && step > _min_line_search_step_size);
 
     // has norm improved or is the step still above minumum search step size?
     return (rnorm <= rnorm_prev || step > _min_line_search_step_size);
   }
   else if (_line_search_method == LineSearchMethod::Bisection)
   {
     unsigned int count = 0;
     Real step_a = 0.0;
     Real step_b = 1.0;
     Real step = 1.0;
     Real s_m = 1000.0;
     Real rnorm = 1000.0;
 
     calculateResidual();
     auto s_b = _residual_tensor.doubleContraction(dpk2);
     const auto rnorm1 = _residual_tensor.L2norm();
     _pk2[_qp] = _pk2[_qp] - dpk2;
     calculateResidual();
     auto s_a = _residual_tensor.doubleContraction(dpk2);
     const auto rnorm0 = _residual_tensor.L2norm();
     _pk2[_qp] = _pk2[_qp] + dpk2;
 
     if ((rnorm1 / rnorm0) < _line_search_tolerance || s_a * s_b > 0)
     {
       calculateResidual();
       return true;
     }
 
     while ((rnorm / rnorm0) > _line_search_tolerance && count < _line_search_max_iterations)
     {
       _pk2[_qp] = _pk2[_qp] - step * dpk2;
       step = 0.5 * (step_b + step_a);
       _pk2[_qp] = _pk2[_qp] + step * dpk2;
       calculateResidual();
       s_m = _residual_tensor.doubleContraction(dpk2);
       rnorm = _residual_tensor.L2norm();
 
       if (s_m * s_a < 0.0)
       {
         step_b = step;
         s_b = s_m;
       }
       if (s_m * s_b < 0.0)
       {
         step_a = step;
         s_a = s_m;
       }
       count++;
     }
 
     // below tolerance and max iterations?
     return ((rnorm / rnorm0) < _line_search_tolerance && count < _line_search_max_iterations);
   }
   else
     mooseError("Line search method is not provided.");
 }

◆ postSolveQp()

void ComputeMultipleCrystalPlasticityStress::postSolveQp	(	RankTwoTensor &	stress_new,
		RankFourTensor &	jacobian_mult
	)

protected

Save the final stress and internal variable values after the iterative solve.

Definition at line 301 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by updateStress().

 {
   cauchy_stress = _elastic_deformation_gradient * _pk2[_qp] *
                   _elastic_deformation_gradient.transpose() / _elastic_deformation_gradient.det();
 
   calcTangentModuli(jacobian_mult);
 
   _total_lagrangian_strain[_qp] =
       _deformation_gradient[_qp].transpose() * _deformation_gradient[_qp] -
       RankTwoTensor::Identity();
   _total_lagrangian_strain[_qp] = _total_lagrangian_strain[_qp] * 0.5;
 
   // Calculate crystal rotation to track separately
   RankTwoTensor rot;
   _elastic_deformation_gradient.getRUDecompositionRotation(rot);
   _updated_rotation[_qp] = rot * _crysrot[_qp];
 }

◆ preSolveQp()

void ComputeMultipleCrystalPlasticityStress::preSolveQp ( )

protected

Reset the PK2 stress and the inverse deformation gradient to old values and provide an interface for inheriting classes to reset material properties.

Definition at line 269 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by updateStress().

 {
   for (unsigned int i = 0; i < _num_models; ++i)
     _models[i]->setInitialConstitutiveVariableValues();
 
   _pk2[_qp] = _pk2_old[_qp];
   _inverse_plastic_deformation_grad_old = _plastic_deformation_gradient_old[_qp].inverse();
 }

◆ solveQp()

void ComputeMultipleCrystalPlasticityStress::solveQp ( )

protected

Solve the stress and internal state variables (e.g.

slip increment, slip system resistance) at each qp points

Definition at line 279 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by updateStress().

 {
   for (unsigned int i = 0; i < _num_models; ++i)
   {
     _models[i]->setSubstepConstitutiveVariableValues();
     _models[i]->calculateSlipResistance();
   }
 
   _inverse_plastic_deformation_grad = _inverse_plastic_deformation_grad_old;
 
   solveStateVariables();
   if (_convergence_failed)
     return; // pop back up and take a smaller substep
 
   for (unsigned int i = 0; i < _num_models; ++i)
     _models[i]->updateSubstepConstitutiveVariableValues();
 
   // save off the old F^{p} inverse now that have converged on the stress and state variables
   _inverse_plastic_deformation_grad_old = _inverse_plastic_deformation_grad;
 }

◆ solveStateVariables()

void ComputeMultipleCrystalPlasticityStress::solveStateVariables ( )

protected

Solves the internal variables stress as a function of the slip specified by the constitutive model defined in the inheriting class.

Definition at line 321 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by solveQp().

 {
   unsigned int iteration;
   bool iter_flag = true;
 
   iteration = 0;
   // Check for slip system resistance update tolerance
   do
   {
     solveStress();
     if (_convergence_failed)
       return;
 
     _plastic_deformation_gradient[_qp] =
         _inverse_plastic_deformation_grad.inverse(); // the postSoveStress
 
     // Update slip system resistance and state variable after the stress has been finalized
     // We loop through all the models for each calculation
     // in order to make sure that when coupling appears, the state variables are updated based on
     // the same components
     for (unsigned int i = 0; i < _num_models; ++i)
       _models[i]->cacheStateVariablesBeforeUpdate();
 
     for (unsigned int i = 0; i < _num_models; ++i)
       _models[i]->calculateStateVariableEvolutionRateComponent();
 
     for (unsigned int i = 0; i < _num_models; ++i)
       if (!_models[i]->updateStateVariables())
         _convergence_failed = true;
 
     for (unsigned int i = 0; i < _num_models; ++i)
       _models[i]->calculateSlipResistance();
 
     if (_convergence_failed)
       return;
 
     for (unsigned int i = 0; i < _num_models; ++i)
     {
       // iter_flag = false, stop iteration if all values are converged and good to go
       // iter_flag = true, continue iteration if any value is not converged
       if (!_models[i]->areConstitutiveStateVariablesConverged())
       {
         iter_flag = true;
         break;
       }
       else
         iter_flag = false; // iter_flag = false, stop iteration only when all models returns true
     }
 
     if (iter_flag)
     {
       if (_print_convergence_message)
         mooseWarning("ComputeMultipleCrystalPlasticityStress: State variables (or the system "
                      "resistance) did not converge at element ",
                      _current_elem->id(),
                      " and qp ",
                      _qp,
                      "\n");
     }
     iteration++;
   } while (iter_flag && iteration < _maxiterg);
 
   if (iteration == _maxiterg)
   {
     if (_print_convergence_message)
       mooseWarning(
           "ComputeMultipleCrystalPlasticityStress: Hardness Integration error. Reached the "
           "maximum number of iterations to solve for the state variables at element ",
           _current_elem->id(),
           " and qp ",
           _qp,
           "\n");
 
     _convergence_failed = true;
   }
 }

◆ solveStress()

void ComputeMultipleCrystalPlasticityStress::solveStress ( )

protected

solves for stress, updates plastic deformation gradient.

Definition at line 399 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by solveStateVariables().

 {
   unsigned int iteration = 0;
   RankTwoTensor dpk2;
   Real rnorm, rnorm0, rnorm_prev;
 
   // Calculate stress residual
   calculateResidualAndJacobian();
   if (_convergence_failed)
   {
     if (_print_convergence_message)
       mooseWarning("ComputeMultipleCrystalPlasticityStress: the slip increment exceeds tolerance "
                    "at element ",
                    _current_elem->id(),
                    " and Gauss point ",
                    _qp);
 
     return;
   }
 
   rnorm = _residual_tensor.L2norm();
   rnorm0 = rnorm;
 
   // Check for stress residual tolerance; different from user object version which
   // compares the absolute tolerance of only the original rnorm value
   while (rnorm > _rtol * rnorm0 && rnorm > _abs_tol && iteration < _maxiter)
   {
     // Calculate stress increment
     dpk2 = -_jacobian.invSymm() * _residual_tensor;
     _pk2[_qp] = _pk2[_qp] + dpk2;
 
     calculateResidualAndJacobian();
 
     if (_convergence_failed)
     {
       if (_print_convergence_message)
         mooseWarning("ComputeMultipleCrystalPlasticityStress: the slip increment exceeds tolerance "
                      "at element ",
                      _current_elem->id(),
                      " and Gauss point ",
                      _qp);
 
       return;
     }
 
     rnorm_prev = rnorm;
     rnorm = _residual_tensor.L2norm();
 
     if (_use_line_search && rnorm > rnorm_prev && !lineSearchUpdate(rnorm_prev, dpk2))
     {
       if (_print_convergence_message)
         mooseWarning("ComputeMultipleCrystalPlasticityStress: Failed with line search");
 
       _convergence_failed = true;
       return;
     }
 
     if (_use_line_search)
       rnorm = _residual_tensor.L2norm();
 
     iteration++;
   }
 
   if (iteration >= _maxiter)
   {
     if (_print_convergence_message)
       mooseWarning("ComputeMultipleCrystalPlasticityStress: Stress Integration error rmax = ",
                    rnorm,
                    " and the tolerance is ",
                    _rtol * rnorm0,
                    " when the rnorm0 value is ",
                    rnorm0,
                    " for element ",
                    _current_elem->id(),
                    " and qp ",
                    _qp);
 
     _convergence_failed = true;
   }
 }

◆ updateStress()

void ComputeMultipleCrystalPlasticityStress::updateStress	(	RankTwoTensor &	cauchy_stress,
		RankFourTensor &	jacobian_mult
	)

protectedvirtual

Updates the stress (PK2) at a quadrature point by calling constiutive relationship as defined in a child class Solves stress residual equation using Newton - Raphson: Updates slip system resistances iteratively.

Definition at line 201 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by computeQpStress().

 {
   // Does not support face/boundary material property calculation
   if (isBoundaryMaterial())
     return;
 
   // Initialize substepping variables
   unsigned int substep_iter = 1;
   unsigned int num_substep = 1;
 
   _temporary_deformation_gradient_old = _deformation_gradient_old[_qp];
   if (_temporary_deformation_gradient_old.det() == 0)
     _temporary_deformation_gradient_old.addIa(1.0);
 
   _delta_deformation_gradient = _deformation_gradient[_qp] - _temporary_deformation_gradient_old;
 
   // Loop through all models and calculate the schmid tensor for the current state of the crystal
   // lattice
   // Not sure if we should pass in the updated or the original rotation here
   // If not, then we should not need to compute the flow direction every iteration here
   for (unsigned int i = 0; i < _num_models; ++i)
     _models[i]->calculateFlowDirection(_crysrot[_qp]);
 
   do
   {
     _convergence_failed = false;
     preSolveQp();
 
     _substep_dt = _dt / num_substep;
     for (unsigned int i = 0; i < _num_models; ++i)
       _models[i]->setSubstepDt(_substep_dt);
 
     // calculate F^{eigen} only when we have eigenstrain
     _inverse_eigenstrain_deformation_grad.zero();
     _inverse_eigenstrain_deformation_grad.addIa(1.0);
     if (_num_eigenstrains)
       calculateEigenstrainDeformationGrad();
 
     for (unsigned int istep = 0; istep < num_substep; ++istep)
     {
       _temporary_deformation_gradient =
           (static_cast<Real>(istep) + 1) / num_substep * _delta_deformation_gradient;
       _temporary_deformation_gradient += _temporary_deformation_gradient_old;
 
       solveQp();
 
       if (_convergence_failed)
       {
         if (_print_convergence_message)
           mooseWarning(
               "The crystal plasticity constitutive model has failed to converge. Increasing "
               "the number of substeps.");
 
         substep_iter++;
         num_substep *= 2;
         break;
       }
     }
 
     if (substep_iter > _max_substep_iter && _convergence_failed)
       mooseException("ComputeMultipleCrystalPlasticityStress: Constitutive failure");
   } while (_convergence_failed);
 
   postSolveQp(cauchy_stress, jacobian_mult);
 }

◆ validParams()

InputParameters ComputeMultipleCrystalPlasticityStress::validParams ( )

static

Definition at line 20 of file ComputeMultipleCrystalPlasticityStress.C.

 {
   InputParameters params = ComputeFiniteStrainElasticStress::validParams();
 
   params.addClassDescription(
       "Crystal Plasticity base class: handles the Newton iteration over the stress residual and "
       "calculates the Jacobian based on constitutive laws from multiple material classes "
       "that are inherited from CrystalPlasticityStressUpdateBase");
   params.addParam<std::string>(
       "base_name",
       "Optional parameter that allows the user to define multiple mechanics material systems on "
       "the same block, i.e. for multiple phases");
 
   params.addRequiredParam<std::vector<MaterialName>>(
       "crystal_plasticity_models",
       "The material objects to use to calculate crystal plasticity stress and strains.");
   params.addParam<std::vector<MaterialName>>(
       "eigenstrain_names", {}, "The material objects to calculate eigenstrains.");
   params.addParam<MooseEnum>("tan_mod_type",
                              MooseEnum("exact none", "none"),
                              "Type of tangent moduli for preconditioner: default elastic");
   params.addParam<Real>("rtol", 1e-6, "Constitutive stress residual relative tolerance");
   params.addParam<Real>("abs_tol", 1e-6, "Constitutive stress residual absolute tolerance");
   params.addParam<unsigned int>("maxiter", 100, "Maximum number of iterations for stress update");
   params.addParam<unsigned int>(
       "maxiter_state_variable", 100, "Maximum number of iterations for state variable update");
   params.addParam<unsigned int>(
       "maximum_substep_iteration", 1, "Maximum number of substep iteration");
   params.addParam<bool>("use_line_search", false, "Use line search in constitutive update");
   params.addParam<Real>("min_line_search_step_size", 0.01, "Minimum line search step size");
   params.addParam<Real>("line_search_tol", 0.5, "Line search bisection method tolerance");
   params.addParam<unsigned int>(
       "line_search_maxiter", 20, "Line search bisection method maximum number of iteration");
   params.addParam<MooseEnum>("line_search_method",
                              MooseEnum("CUT_HALF BISECTION", "CUT_HALF"),
                              "The method used in line search");
   params.addParam<bool>(
       "print_state_variable_convergence_error_messages",
       false,
       "Whether or not to print warning messages from the crystal plasticity specific convergence "
       "checks on the stress measure and general constitutive model quantinties.");
   return params;
 }

Member Data Documentation

◆ _abs_tol

Real ComputeMultipleCrystalPlasticityStress::_abs_tol

protected

Stress residual equation absolute tolerance.

Definition at line 141 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by solveStress().

◆ _base_name

const std::string ComputeMultipleCrystalPlasticityStress::_base_name

protected

optional parameter to define several mechanical systems on the same block, e.g. multiple phases

Definition at line 133 of file ComputeMultipleCrystalPlasticityStress.h.

◆ _convergence_failed

bool ComputeMultipleCrystalPlasticityStress::_convergence_failed

protected

Flag to check whether convergence is achieved or if substepping is needed.

Definition at line 224 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateResidual(), calculateResidualAndJacobian(), ComputeMultipleCrystalPlasticityStress(), solveQp(), solveStateVariables(), solveStress(), and updateStress().

◆ _crysrot

const MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_crysrot

protected

Crystal rotation in the original, or reference, configuration as defined by Euler angle arguments in the ComputeElasticityTensor classes.

Definition at line 210 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by postSolveQp(), and updateStress().

◆ _deformation_gradient

const MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_deformation_gradient

protected

Total deformation gradient RankTwoTensor for the crystal.

Definition at line 188 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by postSolveQp(), and updateStress().

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_deformation_gradient_old

protected

Definition at line 189 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by updateStress().

◆ _delta_deformation_gradient

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_delta_deformation_gradient

protected

Used for substepping; Uniformly divides the increment in deformation gradient.

Definition at line 227 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by updateStress().

◆ _dfgrd_scale_factor

Real ComputeMultipleCrystalPlasticityStress::_dfgrd_scale_factor

protected

Scales the substepping increment to obtain deformation gradient at a substep iteration.

Definition at line 232 of file ComputeMultipleCrystalPlasticityStress.h.

◆ _eigenstrain_deformation_gradient

MaterialProperty<RankTwoTensor>* ComputeMultipleCrystalPlasticityStress::_eigenstrain_deformation_gradient

protected

Generalized eigenstrain deformation gradient RankTwoTensor for the crystal.

Definition at line 183 of file ComputeMultipleCrystalPlasticityStress.h.

◆ _eigenstrain_deformation_gradient_old

const MaterialProperty<RankTwoTensor>* ComputeMultipleCrystalPlasticityStress::_eigenstrain_deformation_gradient_old

protected

Definition at line 184 of file ComputeMultipleCrystalPlasticityStress.h.

◆ _eigenstrains

std::vector<ComputeCrystalPlasticityEigenstrainBase *> ComputeMultipleCrystalPlasticityStress::_eigenstrains

protected

The user supplied cyrstal plasticity eigenstrains.

Definition at line 130 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateEigenstrainDeformationGrad(), computeQpStress(), initialSetup(), and initQpStatefulProperties().

◆ _elastic_deformation_gradient

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_elastic_deformation_gradient

protected

Definition at line 214 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateJacobian(), calculateResidual(), elastoPlasticTangentModuli(), and postSolveQp().

◆ _elastic_strain

MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_elastic_strain

protectedinherited

Elastic strain material property.

Definition at line 50 of file ComputeGeneralStressBase.h.

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

◆ _elastic_strain_old

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainElasticStress::_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 51 of file ComputeFiniteStrainElasticStress.h.

Referenced by ComputeFiniteStrainElasticStress::computeQpStress().

◆ _elasticity_tensor

const MaterialProperty<RankFourTensor>& ComputeMultipleCrystalPlasticityStress::_elasticity_tensor

protected

Elasticity tensor as defined by a separate class.

Definition at line 136 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateJacobian(), calculateResidual(), elasticTangentModuli(), and elastoPlasticTangentModuli().

◆ _elasticity_tensor_name

const std::string ComputeFiniteStrainElasticStress::_elasticity_tensor_name

protectedinherited

Name of the elasticity tensor material property.

Definition at line 33 of file ComputeFiniteStrainElasticStress.h.

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

◆ _extra_stress

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_extra_stress

protectedinherited

Extra stress tensor.

Definition at line 53 of file ComputeGeneralStressBase.h.

Referenced by ComputeGeneralStressBase::computeQpProperties().

◆ _initial_stress_fcn

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

protectedinherited

initial stress components

Definition at line 56 of file ComputeGeneralStressBase.h.

◆ _inverse_eigenstrain_deformation_grad

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_inverse_eigenstrain_deformation_grad

protected

Definition at line 217 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateEigenstrainDeformationGrad(), calculateJacobian(), calculateResidual(), elastoPlasticTangentModuli(), initQpStatefulProperties(), and updateStress().

◆ _inverse_plastic_deformation_grad

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_inverse_plastic_deformation_grad

protected

Definition at line 215 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateResidual(), elastoPlasticTangentModuli(), solveQp(), and solveStateVariables().

◆ _inverse_plastic_deformation_grad_old

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_inverse_plastic_deformation_grad_old

protected

Definition at line 216 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateJacobian(), calculateResidual(), preSolveQp(), and solveQp().

◆ _jacobian

RankFourTensor ComputeMultipleCrystalPlasticityStress::_jacobian

protected

Jacobian tensor.

Definition at line 146 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateJacobian(), and solveStress().

◆ _Jacobian_mult

MaterialProperty<RankFourTensor>& ComputeGeneralStressBase::_Jacobian_mult

protectedinherited

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

Definition at line 59 of file ComputeGeneralStressBase.h.

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

◆ _line_search_max_iterations

unsigned int ComputeMultipleCrystalPlasticityStress::_line_search_max_iterations

protected

Line search bisection method maximum iteration number.

Definition at line 172 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by lineSearchUpdate().

◆ _line_search_method

enum ComputeMultipleCrystalPlasticityStress::LineSearchMethod ComputeMultipleCrystalPlasticityStress::_line_search_method

protected

Referenced by lineSearchUpdate().

◆ _line_search_tolerance

Real ComputeMultipleCrystalPlasticityStress::_line_search_tolerance

protected

Line search bisection method tolerance.

Definition at line 169 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by lineSearchUpdate().

◆ _max_substep_iter

unsigned int ComputeMultipleCrystalPlasticityStress::_max_substep_iter

protected

Maximum number of substep iterations.

Definition at line 157 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by updateStress().

◆ _maxiter

unsigned int ComputeMultipleCrystalPlasticityStress::_maxiter

protected

Maximum number of iterations for stress update.

Definition at line 149 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by solveStress().

◆ _maxiterg

unsigned int ComputeMultipleCrystalPlasticityStress::_maxiterg

protected

Maximum number of iterations for internal variable update.

Definition at line 151 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by solveStateVariables().

◆ _mechanical_strain

const MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_mechanical_strain

protectedinherited

Mechanical strain material property.

Definition at line 46 of file ComputeGeneralStressBase.h.

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

◆ _min_line_search_step_size

Real ComputeMultipleCrystalPlasticityStress::_min_line_search_step_size

protected

Minimum line search step size.

Definition at line 166 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by lineSearchUpdate().

◆ _models

std::vector<CrystalPlasticityStressUpdateBase *> ComputeMultipleCrystalPlasticityStress::_models

protected

The user supplied cyrstal plasticity consititutive models.

Definition at line 124 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateJacobian(), calculateResidual(), computeQpStress(), initialSetup(), initQpStatefulProperties(), preSolveQp(), solveQp(), solveStateVariables(), and updateStress().

◆ _num_eigenstrains

const unsigned ComputeMultipleCrystalPlasticityStress::_num_eigenstrains

protected

number of eigenstrains

Definition at line 127 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateEigenstrainDeformationGrad(), calculateJacobian(), calculateResidual(), computeQpStress(), initialSetup(), initQpStatefulProperties(), and updateStress().

◆ _num_models

const unsigned ComputeMultipleCrystalPlasticityStress::_num_models

protected

number of plastic models

Definition at line 121 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateJacobian(), calculateResidual(), computeQpStress(), initialSetup(), initQpStatefulProperties(), preSolveQp(), solveQp(), solveStateVariables(), and updateStress().

◆ _pk2

MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_pk2

protected

Second Piola-Kirchoff stress measure.

Definition at line 193 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateResidual(), elastoPlasticTangentModuli(), initQpStatefulProperties(), lineSearchUpdate(), postSolveQp(), preSolveQp(), and solveStress().

◆ _pk2_old

const MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_pk2_old

protected

Definition at line 194 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by preSolveQp().

◆ _plastic_deformation_gradient

MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_plastic_deformation_gradient

protected

Plastic deformation gradient RankTwoTensor for the crystal.

Definition at line 178 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by initQpStatefulProperties(), and solveStateVariables().

◆ _plastic_deformation_gradient_old

const MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_plastic_deformation_gradient_old

protected

Definition at line 179 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by preSolveQp().

◆ _print_convergence_message

const bool ComputeMultipleCrystalPlasticityStress::_print_convergence_message

protected

Flag to print to console warning messages on stress, constitutive model convergence.

Definition at line 221 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by solveStateVariables(), solveStress(), and updateStress().

◆ _residual_tensor

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_residual_tensor

protected

Residual tensor.

Definition at line 144 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateResidual(), lineSearchUpdate(), and solveStress().

◆ _rotation_increment

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainElasticStress::_rotation_increment

protectedinherited

Rotation increment material property.

Definition at line 43 of file ComputeFiniteStrainElasticStress.h.

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

◆ _rotation_total

MaterialProperty<RankTwoTensor>& ComputeFiniteStrainElasticStress::_rotation_total

protectedinherited

Rotation up to current step "n" to compute anisotropic elasticity tensor.

Definition at line 37 of file ComputeFiniteStrainElasticStress.h.

Referenced by ComputeFiniteStrainElasticStress::computeQpStress(), and ComputeFiniteStrainElasticStress::initQpStatefulProperties().

◆ _rotation_total_old

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainElasticStress::_rotation_total_old

protectedinherited

Rotation up to "n - 1" (previous) step to compute anisotropic elasticity tensor.

Definition at line 39 of file ComputeFiniteStrainElasticStress.h.

Referenced by ComputeFiniteStrainElasticStress::computeQpStress().

◆ _rtol

Real ComputeMultipleCrystalPlasticityStress::_rtol

protected

Stress residual equation relative tolerance.

Definition at line 139 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by solveStress().

◆ _strain_increment

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainElasticStress::_strain_increment

protectedinherited

Strain increment material property.

Definition at line 41 of file ComputeFiniteStrainElasticStress.h.

Referenced by ComputeFiniteStrainElasticStress::computeQpStress().

◆ _stress

MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_stress

protectedinherited

Stress material property.

Definition at line 48 of file ComputeGeneralStressBase.h.

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

◆ _stress_old

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainElasticStress::_stress_old

protectedinherited

Old state of the stress tensor material property.

Definition at line 45 of file ComputeFiniteStrainElasticStress.h.

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

◆ _substep_dt

Real ComputeMultipleCrystalPlasticityStress::_substep_dt

protected

time step size during substepping

Definition at line 160 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateEigenstrainDeformationGrad(), and updateStress().

◆ _tan_mod_type

enum ComputeMultipleCrystalPlasticityStress::TangentModuliType ComputeMultipleCrystalPlasticityStress::_tan_mod_type

protected

Referenced by calcTangentModuli().

◆ _temporary_deformation_gradient

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_temporary_deformation_gradient

protected

Helper deformation gradient tensor variables used in iterative solve.

Definition at line 213 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by calculateJacobian(), calculateResidual(), and updateStress().

◆ _temporary_deformation_gradient_old

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_temporary_deformation_gradient_old

protected

Definition at line 228 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by updateStress().

◆ _total_lagrangian_strain

MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_total_lagrangian_strain

protected

Lagrangian total strain measure for the entire crystal.

Definition at line 198 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

◆ _updated_rotation

MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_updated_rotation

protected

Tracks the rotation of the crystal during deformation Note: this rotation tensor is not applied to the crystal lattice.

Definition at line 204 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by initQpStatefulProperties(), and postSolveQp().

◆ _use_line_search

bool ComputeMultipleCrystalPlasticityStress::_use_line_search

protected

Flag to activate line serach.

Definition at line 163 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by solveStress().

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

solid_mechanics/include/materials/crystal_plasticity/ComputeMultipleCrystalPlasticityStress.h
solid_mechanics/src/materials/crystal_plasticity/ComputeMultipleCrystalPlasticityStress.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Types

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ LineSearchMethod

◆ TangentModuliType

Constructor & Destructor Documentation

◆ ComputeMultipleCrystalPlasticityStress()

Member Function Documentation

◆ calcTangentModuli()

◆ calculateEigenstrainDeformationGrad()

◆ calculateJacobian()

◆ calculateResidual()

◆ calculateResidualAndJacobian()

◆ computeQpProperties()

◆ computeQpStress()

◆ elasticTangentModuli()

◆ elastoPlasticTangentModuli()

◆ hasGuaranteedMaterialProperty()

◆ initialSetup()

◆ initQpStatefulProperties()

◆ lineSearchUpdate()

◆ postSolveQp()

◆ preSolveQp()

◆ solveQp()

◆ solveStateVariables()

◆ solveStress()

◆ updateStress()

◆ validParams()

Member Data Documentation

◆ _abs_tol

◆ _base_name

◆ _convergence_failed

◆ _crysrot

◆ _deformation_gradient

◆ _deformation_gradient_old

◆ _delta_deformation_gradient

◆ _dfgrd_scale_factor

◆ _eigenstrain_deformation_gradient

◆ _eigenstrain_deformation_gradient_old

◆ _eigenstrains

◆ _elastic_deformation_gradient

◆ _elastic_strain

◆ _elastic_strain_old

◆ _elasticity_tensor

◆ _elasticity_tensor_name

◆ _extra_stress

◆ _initial_stress_fcn

◆ _inverse_eigenstrain_deformation_grad

◆ _inverse_plastic_deformation_grad

◆ _inverse_plastic_deformation_grad_old

◆ _jacobian

◆ _Jacobian_mult

◆ _line_search_max_iterations

◆ _line_search_method

◆ _line_search_tolerance

◆ _max_substep_iter

◆ _maxiter

◆ _maxiterg

◆ _mechanical_strain

◆ _min_line_search_step_size

◆ _models

◆ _num_eigenstrains

◆ _num_models

◆ _pk2

◆ _pk2_old

◆ _plastic_deformation_gradient

◆ _plastic_deformation_gradient_old

◆ _print_convergence_message

◆ _residual_tensor

◆ _rotation_increment

◆ _rotation_total

◆ _rotation_total_old