https://mooseframework.inl.gov
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Types | Protected Member Functions | Protected Attributes | Static Protected Attributes | List of all members
ComputeMultipleCrystalPlasticityStress Class Reference

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  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_typegetElementID (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_typegetElementIDNeighbor (const std::string &id_parameter_name, unsigned int comp=0) const override
 
virtual const dof_id_typegetElementIDByName (const std::string &id_parameter_name) const override
 
virtual const dof_id_typegetElementIDNeighborByName (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)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false, bool no_dep=false)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false)
 
MaterialBasegetMaterialByName (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
 
 CreateMooseEnumClass (ConstantTypeEnum, NONE, ELEMENT, SUBDOMAIN)
 
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 (const 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 () const
 
void checkStatefulSanity () const
 
std::set< OutputName > getOutputs ()
 
bool hasStatefulProperties () const
 
bool hasRestoredProperties () 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 hasActiveProperties ()
 
bool forceStatefulInit () const
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
bool isKokkosObject () const
 
MooseAppgetMooseApp () const
 
const std::string & type () const
 
const std::string & name () const
 
std::string typeAndName () const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
MooseObjectName uniqueName () const
 
const InputParametersparameters () const
 
const hit::Node * getHitNode () const
 
bool hasBase () const
 
const std::string & getBase () const
 
const TgetParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const TqueryParam (const std::string &name) const
 
const TgetRenamedParam (const std::string &old_name, const std::string &new_name) const
 
T getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool haveParameter (const std::string &name) const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &name) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (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
 
std::string messagePrefix (const bool hit_prefix=true) const
 
std::string errorPrefix (const std::string &) const
 
void mooseError (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecatedNoTrace (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
void callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) 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 ExecFlagEnumgetExecuteOnEnum () 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 ADVariableValuegetADDefaultValue (const std::string &var_name) const
 
const FunctiongetFunction (const std::string &name) const
 
const FunctiongetFunctionByName (const FunctionName &name) const
 
bool hasFunction (const std::string &param_name) const
 
bool hasFunctionByName (const FunctionName &name) const
 
Moose::Kokkos::Function getKokkosFunction (const std::string &name) const
 
const TgetKokkosFunction (const std::string &name) const
 
Moose::Kokkos::Function getKokkosFunctionByName (const FunctionName &name) const
 
const TgetKokkosFunctionByName (const FunctionName &name) const
 
bool hasKokkosFunction (const std::string &param_name) const
 
bool hasKokkosFunctionByName (const FunctionName &name) const
 
UserObjectName getUserObjectName (const std::string &param_name) const
 
const TgetUserObject (const std::string &param_name, bool is_dependency=true) const
 
const TgetUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const UserObjectBasegetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectBasegetUserObjectBaseByName (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 VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOldByName (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_typegetAllElemIDs (unsigned int elem_id_index) const
 
std::set< dof_id_typegetElemIDsOnBlocks (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 ADVectorVariableValuegetADDefaultVectorValue (const std::string &var_name) const
 
const ADVariableGradientgetADDefaultGradient () const
 
const ADVectorVariableGradientgetADDefaultVectorGradient () const
 
const ADVariableSecondgetADDefaultSecond () const
 
const ADVectorVariableCurlgetADDefaultCurl () const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledVectorTagGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledVectorTagGradientsByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagNodalValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagDofValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledVectorTagGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledVectorTagGradients (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagNodalValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagDofValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledValues (const std::string &var_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradients (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValues (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValues (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledValuesOld (const std::string &var_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradientsOld (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValuesOld (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValuesOld (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradientsOlder (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledDots (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalDots (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledVectorTagGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledVectorTagGradientsByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagNodalValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagDofValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledVectorTagGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledVectorTagGradients (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagNodalValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagDofValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValues (const std::string &var_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradients (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValues (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValues (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValuesOld (const std::string &var_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradientsOld (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValuesOld (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValuesOld (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradientsOlder (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDots (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalDots (const std::string &var_name) const
 
Moose::Kokkos::Scalar< const RealkokkosCoupledDotDu (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosZeroValue () const
 
Moose::Kokkos::VariableGradient kokkosZeroGradient () const
 
Moose::Kokkos::VariableValue kokkosZeroNodalValue () const
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyByName (const std::string &prop_name_in)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOldByName (const std::string &prop_name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOlderByName (const std::string &prop_name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialProperty (const std::string &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOld (const std::string &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOlder (const std::string &name)
 
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > > getBlockMaterialProperty (const MaterialPropertyName &name)
 
std::pair< Moose::Kokkos::MaterialProperty< T, dimension >, std::set< SubdomainID > > getKokkosBlockMaterialProperty (const MaterialPropertyName &name)
 
std::set< SubdomainIDgetMaterialPropertyBlocks (const std::string &name)
 
std::vector< SubdomainName > getMaterialPropertyBlockNames (const std::string &name)
 
std::set< BoundaryIDgetMaterialPropertyBoundaryIDs (const std::string &name)
 
std::vector< BoundaryName > getMaterialPropertyBoundaryNames (const std::string &name)
 
void checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)
 
void statefulPropertiesAllowed (bool)
 
virtual 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 DistributiongetDistribution (const std::string &name) const
 
const TgetDistribution (const std::string &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const TgetDistribution (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const TgetDistributionByName (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const TgetDistributionByName (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 PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (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 hasKokkosMaterialProperty (const std::string &name)
 
bool hasKokkosMaterialProperty (const std::string &name)
 
bool hasKokkosMaterialPropertyByName (const std::string &name)
 
bool hasKokkosMaterialPropertyByName (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
 
PenetrationLocatorgetPenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
 
PenetrationLocatorgetQuadraturePenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
 
NearestNodeLocatorgetNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)
 
NearestNodeLocatorgetQuadratureNearestNodeLocator (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 void callMooseError (MooseApp *const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node *node, const bool show_trace=true)
 
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, NameFunc &&name_func)
 
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
 
 usingCombinedWarningSolutionWarnings
 
const ConsoleStream _console
 

Static Public Attributes

static const std::string type_param
 
static const std::string name_param
 
static const std::string unique_name_param
 
static const std::string app_param
 
static const std::string moose_base_param
 
static const std::string kokkos_object_param
 
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 MaterialDatamaterialData () const override
 
virtual MaterialDatamaterialData () override
 
virtual Moose::MaterialDataType materialDataType () override
 
virtual const QBase & qRule () const override
 
virtual void resetQpProperties ()
 
virtual const FEProblemBasemiProblem () const
 
virtual FEProblemBasemiProblem ()
 
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 &)
 
void flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
 
InvalidSolutionID registerInvalidSolutionInternal (const std::string &message, const bool warning) const
 
virtual bool hasBlockMaterialPropertyHelper (const std::string &prop_name)
 
void initializeBlockRestrictable (const MooseObject *moose_object)
 
void initializeKokkosBlockRestrictable ()
 
Moose::CoordinateSystemType getBlockCoordSystem ()
 
KOKKOS_FUNCTION dof_id_type numKokkosBlockElements () const
 
KOKKOS_FUNCTION dof_id_type numKokkosBlockNodes () const
 
KOKKOS_FUNCTION dof_id_type numKokkosBlockSides () const
 
KOKKOS_FUNCTION ContiguousElementID kokkosBlockElementID (Moose::Kokkos::ThreadID tid) const
 
KOKKOS_FUNCTION ContiguousElementID kokkosBlockNodeID (Moose::Kokkos::ThreadID tid) const
 
KOKKOS_FUNCTION auto kokkosBlockElementSideID (Moose::Kokkos::ThreadID tid) const
 
bool hasBoundaryMaterialPropertyHelper (const std::string &prop_name) const
 
KOKKOS_FUNCTION dof_id_type numKokkosBoundaryNodes () const
 
KOKKOS_FUNCTION dof_id_type numKokkosBoundarySides () const
 
KOKKOS_FUNCTION ContiguousNodeID kokkosBoundaryNodeID (Moose::Kokkos::ThreadID tid) const
 
KOKKOS_FUNCTION auto kokkosBoundaryElementSideID (Moose::Kokkos::ThreadID tid) 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 VariableValuecoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarValue (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 VariableValuecoupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalargetScalarVar (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 UserObjectBase &) const
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &) const
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &) const
 
TdeclareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< TdeclareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const TgetRestartableData (const std::string &data_name) const
 
TdeclareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
TdeclareRecoverableData (const std::string &data_name, Args &&... args)
 
TdeclareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
TdeclareRestartableDataWithObjectNameWithContext (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)
 
const Moose::Functor< T > & getFunctorByName (const std::string &name)
 
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 intcoupledIndices (const std::string &var_name) const
 
virtual const VariableValuecoupledValue (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 VariableValuecoupledValueLower (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledValues (const std::string &var_name) const
 
const ADVariableValueadCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVectorVariableValue *> adCoupledVectorValues (const std::string &var_name) const
 
virtual const VariableValuecoupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagValue (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 ArrayVariableValuecoupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableValuecoupledVectorTagArrayValue (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 VariableGradientcoupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableGradientcoupledVectorTagGradient (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 ArrayVariableGradientcoupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableGradientcoupledVectorTagArrayGradient (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 VariableValuecoupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
const ArrayVariableValuecoupledVectorTagArrayDofValue (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 VariableValuecoupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledMatrixTagValue (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 VectorVariableValuecoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ArrayVariableValue *> coupledArrayValues (const std::string &var_name) const
 
MooseWritableVariablewritableVariable (const std::string &var_name, unsigned int comp=0)
 
virtual VariableValuewritableCoupledValue (const std::string &var_name, unsigned int comp=0)
 
void checkWritableVar (MooseWritableVariable *var)
 
virtual const VariableValuecoupledValueOld (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 VariableValuecoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledValuesOlder (const std::string &var_name) const
 
virtual const VariableValuecoupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradients (const std::string &var_name) const
 
const ADVariableGradientadCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableGradientadCoupledGradientDot (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 ADVectorVariableGradientadCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableSecondadCoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableSecondadCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradientsOld (const std::string &var_name) const
 
virtual const VariableGradientcoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOlder (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableCurladCoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDiv (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledDots (const std::string &var_name) const
 
virtual const VariableValuecoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledDots (const std::string &var_name) const
 
const ADVariableValueadCoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledArrayDotDu (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
const Moose::ADType< T >::typeadCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValues (const std::string &var_name) const
 
virtual const VariableValuecoupledDofValuesOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValuesOld (const std::string &var_name) const
 
virtual const VariableValuecoupledDofValuesOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValuesOlder (const std::string &var_name) const
 
virtual const ArrayVariableValuecoupledArrayDofValues (const std::string &var_name, unsigned int comp=0) const
 
virtual const ADVariableValueadCoupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadZeroValue () const
 
const ADVariableGradientadZeroGradient () const
 
const ADVariableSecondadZeroSecond () 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 MooseVariableFieldBasegetFEVar (const std::string &var_name, unsigned int comp) const
 
const MooseVariableFieldBasegetFieldVar (const std::string &var_name, unsigned int comp) const
 
MooseVariableFieldBasegetFieldVar (const std::string &var_name, unsigned int comp)
 
std::vector< const MooseVariableFieldBase *> getFieldVars (const std::string &var_name) const
 
const TgetVarHelper (const std::string &var_name, unsigned int comp) const
 
TgetVarHelper (const std::string &var_name, unsigned int comp)
 
std::vector< const T *> getVarsHelper (const std::string &var_name) const
 
MooseVariablegetVar (const std::string &var_name, unsigned int comp)
 
const MooseVariablegetVar (const std::string &var_name, unsigned int comp) const
 
VectorMooseVariablegetVectorVar (const std::string &var_name, unsigned int comp)
 
const VectorMooseVariablegetVectorVar (const std::string &var_name, unsigned int comp) const
 
ArrayMooseVariablegetArrayVar (const std::string &var_name, unsigned int comp)
 
const ArrayMooseVariablegetArrayVar (const std::string &var_name, unsigned int comp) const
 
std::vector< TcoupledVectorHelper (const std::string &var_name, const Func &func) const
 
virtual void getKokkosMaterialPropertyHook (const std::string &, const unsigned int)
 
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 > & _updated_rotation_old
 
MaterialProperty< Real > & _misorientation
 Misorientation angle of the crystal during deformation. 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 MaterialPropertyName _declare_suffix
 
const bool & _enabled
 
MooseApp_app
 
Factory_factory
 
ActionFactory_action_factory
 
const std::string & _type
 
const std::string & _name
 
const InputParameters_pars
 
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< MooseVariableField< Real > *> _coupled_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 bool _is_kokkos_object
 
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

strain formulation

Enumerator
CutHalf 
Bisection 

Definition at line 175 of file ComputeMultipleCrystalPlasticityStress.h.

175 { CutHalf, Bisection } _line_search_method;
enum ComputeMultipleCrystalPlasticityStress::LineSearchMethod _line_search_method

◆ TangentModuliType

Type of tangent moduli calculation.

Enumerator
EXACT 
NONE 

Definition at line 154 of file ComputeMultipleCrystalPlasticityStress.h.

EXACT
NONE
enum ComputeMultipleCrystalPlasticityStress::TangentModuliType _tan_mod_type

Constructor & Destructor Documentation

◆ ComputeMultipleCrystalPlasticityStress()

ComputeMultipleCrystalPlasticityStress::ComputeMultipleCrystalPlasticityStress ( const InputParameters parameters)

Definition at line 64 of file ComputeMultipleCrystalPlasticityStress.C.

67  _num_models(getParam<std::vector<MaterialName>>("crystal_plasticity_models").size()),
68  _num_eigenstrains(getParam<std::vector<MaterialName>>("eigenstrain_names").size()),
69  _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
70  _elasticity_tensor(getMaterialPropertyByName<RankFourTensor>(_base_name + "elasticity_tensor")),
71  _rtol(getParam<Real>("rtol")),
72  _abs_tol(getParam<Real>("abs_tol")),
73  _maxiter(getParam<unsigned int>("maxiter")),
74  _maxiterg(getParam<unsigned int>("maxiter_state_variable")),
75  _tan_mod_type(getParam<MooseEnum>("tan_mod_type").getEnum<TangentModuliType>()),
76  _max_substep_iter(getParam<unsigned int>("maximum_substep_iteration")),
77  _use_line_search(getParam<bool>("use_line_search")),
78  _min_line_search_step_size(getParam<Real>("min_line_search_step_size")),
79  _line_search_tolerance(getParam<Real>("line_search_tol")),
80  _line_search_max_iterations(getParam<unsigned int>("line_search_maxiter")),
81  _line_search_method(getParam<MooseEnum>("line_search_method").getEnum<LineSearchMethod>()),
82  _plastic_deformation_gradient(declareProperty<RankTwoTensor>("plastic_deformation_gradient")),
84  getMaterialPropertyOld<RankTwoTensor>("plastic_deformation_gradient")),
86  _num_eigenstrains ? &declareProperty<RankTwoTensor>("eigenstrain_deformation_gradient")
87  : nullptr),
90  ? &getMaterialPropertyOld<RankTwoTensor>("eigenstrain_deformation_gradient")
91  : nullptr),
92  _deformation_gradient(getMaterialProperty<RankTwoTensor>(_base_name + "deformation_gradient")),
94  getMaterialPropertyOld<RankTwoTensor>(_base_name + "deformation_gradient")),
95  _pk2(declareProperty<RankTwoTensor>("second_piola_kirchhoff_stress")),
96  _pk2_old(getMaterialPropertyOld<RankTwoTensor>("second_piola_kirchhoff_stress")),
98  declareProperty<RankTwoTensor>("total_lagrangian_strain")), // Lagrangian strain
99  _updated_rotation(declareProperty<RankTwoTensor>("updated_rotation")),
100  _updated_rotation_old(getMaterialPropertyOld<RankTwoTensor>("updated_rotation")),
101  _misorientation(declareProperty<Real>("misorientation")),
102  _crysrot(getMaterialProperty<RankTwoTensor>(
103  _base_name + "crysrot")), // defined in the elasticity tensor classes for crystal plasticity
104  _print_convergence_message(getParam<bool>("print_state_variable_convergence_error_messages"))
105 {
106  _convergence_failed = false;
107 }
const MaterialProperty< RankTwoTensor > & _pk2_old
MaterialProperty< RankTwoTensor > & _total_lagrangian_strain
Lagrangian total strain measure for the entire crystal.
Real _line_search_tolerance
Line search bisection method tolerance.
bool _convergence_failed
Flag to check whether convergence is achieved or if substepping is needed.
const T & getParam(const std::string &name) const
unsigned int _maxiter
Maximum number of iterations for stress update.
const InputParameters & parameters() const
MaterialProperty< RankTwoTensor > * _eigenstrain_deformation_gradient
Generalized eigenstrain deformation gradient RankTwoTensor for the crystal.
const MaterialProperty< RankTwoTensor > & _crysrot
Crystal rotation in the original, or reference, configuration as defined by Euler angle arguments in ...
enum ComputeMultipleCrystalPlasticityStress::LineSearchMethod _line_search_method
const MaterialProperty< RankTwoTensor > & _deformation_gradient
Total deformation gradient RankTwoTensor for the crystal.
Real _min_line_search_step_size
Minimum line search step size.
MaterialProperty< RankTwoTensor > & _pk2
Second Piola-Kirchoff stress measure.
ComputeFiniteStrainElasticStress(const InputParameters &parameters)
Real _abs_tol
Stress residual equation absolute tolerance.
const MaterialProperty< RankTwoTensor > * _eigenstrain_deformation_gradient_old
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Elasticity tensor as defined by a separate class.
unsigned int _line_search_max_iterations
Line search bisection method maximum iteration number.
MaterialProperty< RankTwoTensor > & _updated_rotation
Tracks the rotation of the crystal during deformation Note: this rotation tensor is not applied to th...
unsigned int _max_substep_iter
Maximum number of substep iterations.
MaterialProperty< RankTwoTensor > & _plastic_deformation_gradient
Plastic deformation gradient RankTwoTensor for the crystal.
enum ComputeMultipleCrystalPlasticityStress::TangentModuliType _tan_mod_type
const MaterialProperty< RankTwoTensor > & _updated_rotation_old
unsigned int _maxiterg
Maximum number of iterations for internal variable update.
const MaterialProperty< RankTwoTensor > & _deformation_gradient_old
const MaterialProperty< RankTwoTensor > & _plastic_deformation_gradient_old
MaterialProperty< Real > & _misorientation
Misorientation angle of the crystal during deformation.
bool isParamValid(const std::string &name) const
const bool _print_convergence_message
Flag to print to console warning messages on stress, constitutive model convergence.
Real _rtol
Stress residual equation relative tolerance.
const std::string _base_name
optional parameter to define several mechanical systems on the same block, e.g. multiple phases ...

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 589 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by postSolveQp().

590 {
591  switch (_tan_mod_type)
592  {
594  elastoPlasticTangentModuli(jacobian_mult);
595  break;
596  default:
597  elasticTangentModuli(jacobian_mult);
598  }
599 }
void elastoPlasticTangentModuli(RankFourTensor &jacobian_mult)
enum ComputeMultipleCrystalPlasticityStress::TangentModuliType _tan_mod_type

◆ calculateEigenstrainDeformationGrad()

void ComputeMultipleCrystalPlasticityStress::calculateEigenstrainDeformationGrad ( )
protected

Calculates the deformation gradient due to eigenstrain.

Definition at line 730 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by updateStress().

731 {
734 
735  for (unsigned int i = 0; i < _num_eigenstrains; ++i)
736  {
737  _eigenstrains[i]->setSubstepDt(_substep_dt);
738  _eigenstrains[i]->computeQpProperties();
739  _inverse_eigenstrain_deformation_grad *= _eigenstrains[i]->getDeformationGradientInverse();
740  }
741  (*_eigenstrain_deformation_gradient)[_qp] = _inverse_eigenstrain_deformation_grad.inverse();
742 }
RankTwoTensorTempl< Real > inverse() const
void addIa(const Real &a)
std::vector< ComputeCrystalPlasticityEigenstrainBase * > _eigenstrains
The user supplied cyrstal plasticity eigenstrains.

◆ 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 554 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by calculateResidualAndJacobian().

555 {
556  // may not need to cache the dfpinvdpk2 here. need to double check
557  RankFourTensor dfedfpinv, deedfe, dfpinvdpk2, dfpinvdpk2_per_model;
558 
560 
561  for (const auto i : make_range(Moose::dim))
562  for (const auto j : make_range(Moose::dim))
563  for (const auto k : make_range(Moose::dim))
564  dfedfpinv(i, j, k, j) = ffeiginv(i, k);
565 
566  for (const auto i : make_range(Moose::dim))
567  for (const auto j : make_range(Moose::dim))
568  for (const auto k : make_range(Moose::dim))
569  {
570  deedfe(i, j, k, i) = deedfe(i, j, k, i) + _elastic_deformation_gradient(k, j) * 0.5;
571  deedfe(i, j, k, j) = deedfe(i, j, k, j) + _elastic_deformation_gradient(k, i) * 0.5;
572  }
573 
574  for (unsigned int i = 0; i < _num_models; ++i)
575  {
576  _models[i]->calculateTotalPlasticDeformationGradientDerivative(
577  dfpinvdpk2_per_model,
581  dfpinvdpk2 += dfpinvdpk2_per_model;
582  }
583 
584  _jacobian =
585  RankFourTensor::IdentityFour() - (_elasticity_tensor[_qp] * deedfe * dfedfpinv * dfpinvdpk2);
586 }
RankTwoTensor _temporary_deformation_gradient
Helper deformation gradient tensor variables used in iterative solve.
static RankFourTensorTempl< Real > IdentityFour()
static constexpr std::size_t dim
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Elasticity tensor as defined by a separate class.
std::vector< CrystalPlasticityStressUpdateBase * > _models
The user supplied cyrstal plasticity consititutive models.
IntRange< T > make_range(T beg, T end)
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
static const std::string k
Definition: NS.h:134

◆ 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 510 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by calculateResidualAndJacobian(), and lineSearchUpdate().

511 {
512  RankTwoTensor ce, elastic_strain, ce_pk2, equivalent_slip_increment_per_model,
513  equivalent_slip_increment, pk2_new;
514 
515  equivalent_slip_increment.zero();
516 
517  // calculate slip rate in order to compute F^{p-1}
518  for (unsigned int i = 0; i < _num_models; ++i)
519  {
520  equivalent_slip_increment_per_model.zero();
521 
522  // calculate shear stress with consideration of contribution from other physics
523  _models[i]->calculateShearStress(
525 
526  _convergence_failed = !_models[i]->calculateSlipRate();
527 
529  return;
530 
531  _models[i]->calculateEquivalentSlipIncrement(equivalent_slip_increment_per_model);
532  equivalent_slip_increment += equivalent_slip_increment_per_model;
533  }
534 
535  RankTwoTensor residual_equivalent_slip_increment =
536  RankTwoTensor::Identity() - equivalent_slip_increment;
538  _inverse_plastic_deformation_grad_old * residual_equivalent_slip_increment;
539 
543 
545 
546  elastic_strain = ce - RankTwoTensor::Identity();
547  elastic_strain *= 0.5;
548 
549  pk2_new = _elasticity_tensor[_qp] * elastic_strain;
550  _residual_tensor = _pk2[_qp] - pk2_new;
551 }
RankTwoTensor _temporary_deformation_gradient
Helper deformation gradient tensor variables used in iterative solve.
bool _convergence_failed
Flag to check whether convergence is achieved or if substepping is needed.
static RankTwoTensorTempl Identity()
MaterialProperty< RankTwoTensor > & _pk2
Second Piola-Kirchoff stress measure.
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Elasticity tensor as defined by a separate class.
RankTwoTensorTempl< Real > transpose() const
std::vector< CrystalPlasticityStressUpdateBase * > _models
The user supplied cyrstal plasticity consititutive models.

◆ calculateResidualAndJacobian()

void ComputeMultipleCrystalPlasticityStress::calculateResidualAndJacobian ( )
protected

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

Definition at line 501 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by solveStress().

502 {
505  return;
507 }
bool _convergence_failed
Flag to check whether convergence is achieved or if substepping is needed.
void calculateResidual()
Calculate stress residual as the difference between the stored material property PK2 stress and the e...
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{PK...

◆ computeQpProperties()

void ComputeGeneralStressBase::computeQpProperties ( )
overrideprotectedvirtualinherited

Reimplemented from DerivativeMaterialInterface< Material >.

Definition at line 44 of file ComputeGeneralStressBase.C.

45 {
47 
48  // Add in extra stress
50 }
const MaterialProperty< RankTwoTensor > & _extra_stress
Extra stress tensor.
MaterialProperty< RankTwoTensor > & _stress
Stress material property.
virtual void computeQpStress()=0
Compute the stress and store it in the _stress material property for the current quadrature point...

◆ 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 188 of file ComputeMultipleCrystalPlasticityStress.C.

189 {
190  for (unsigned int i = 0; i < _num_models; ++i)
191  {
192  _models[i]->setQp(_qp);
193  _models[i]->setMaterialVectorSize();
194  }
195 
196  for (unsigned int i = 0; i < _num_eigenstrains; ++i)
197  _eigenstrains[i]->setQp(_qp);
198 
199  updateStress(_stress[_qp], _Jacobian_mult[_qp]); // This is NOT the exact jacobian
200 }
MaterialProperty< RankFourTensor > & _Jacobian_mult
derivative of stress w.r.t. strain (_dstress_dstrain)
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 ch...
std::vector< CrystalPlasticityStressUpdateBase * > _models
The user supplied cyrstal plasticity consititutive models.
std::vector< ComputeCrystalPlasticityEigenstrainBase * > _eigenstrains
The user supplied cyrstal plasticity eigenstrains.
MaterialProperty< RankTwoTensor > & _stress
Stress material property.

◆ elasticTangentModuli()

void ComputeMultipleCrystalPlasticityStress::elasticTangentModuli ( RankFourTensor jacobian_mult)
protected

Definition at line 648 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by calcTangentModuli().

649 {
650  // update jacobian_mult
651  jacobian_mult = _elasticity_tensor[_qp];
652 }
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Elasticity tensor as defined by a separate class.

◆ elastoPlasticTangentModuli()

void ComputeMultipleCrystalPlasticityStress::elastoPlasticTangentModuli ( RankFourTensor jacobian_mult)
protected

Definition at line 602 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by calcTangentModuli().

603 {
604  RankFourTensor tan_mod;
605  RankTwoTensor pk2fet, fepk2, feiginvfpinv;
606  RankFourTensor deedfe, dsigdpk2dfe, dfedf;
607 
608  // Fill in the matrix stiffness material property
609  for (const auto i : make_range(Moose::dim))
610  for (const auto j : make_range(Moose::dim))
611  for (const auto k : make_range(Moose::dim))
612  {
613  deedfe(i, j, k, i) = deedfe(i, j, k, i) + _elastic_deformation_gradient(k, j) * 0.5;
614  deedfe(i, j, k, j) = deedfe(i, j, k, j) + _elastic_deformation_gradient(k, i) * 0.5;
615  }
616 
617  usingTensorIndices(i_, j_, k_, l_);
618  dsigdpk2dfe = _elastic_deformation_gradient.times<i_, k_, j_, l_>(_elastic_deformation_gradient) *
619  _elasticity_tensor[_qp] * deedfe;
620 
623 
624  for (const auto i : make_range(Moose::dim))
625  for (const auto j : make_range(Moose::dim))
626  for (const auto l : make_range(Moose::dim))
627  {
628  tan_mod(i, j, i, l) += pk2fet(l, j);
629  tan_mod(i, j, j, l) += fepk2(i, l);
630  }
631 
632  tan_mod += dsigdpk2dfe;
633 
634  const auto je = _elastic_deformation_gradient.det();
635  if (je > 0.0)
636  tan_mod /= je;
637 
639  for (const auto i : make_range(Moose::dim))
640  for (const auto j : make_range(Moose::dim))
641  for (const auto l : make_range(Moose::dim))
642  dfedf(i, j, i, l) = feiginvfpinv(l, j);
643 
644  jacobian_mult = tan_mod * dfedf;
645 }
static constexpr std::size_t dim
MaterialProperty< RankTwoTensor > & _pk2
Second Piola-Kirchoff stress measure.
const MaterialProperty< RankFourTensor > & _elasticity_tensor
Elasticity tensor as defined by a separate class.
RankTwoTensorTempl< Real > transpose() const
RankFourTensorTempl< Real > times(const RankTwoTensorTempl< Real > &b) const
IntRange< T > make_range(T beg, T end)
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
static const std::string k
Definition: NS.h:134

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

30 {
32  mooseError("hasGuaranteedMaterialProperty() needs to be called in initialSetup()");
33 
34  // Reference to MaterialWarehouse for testing and retrieving block ids
35  const auto & warehouse = _gc_feproblem->getMaterialWarehouse();
36 
37  // Complete set of ids that this object is active
38  const auto & ids = (_gc_block_restrict && _gc_block_restrict->blockRestricted())
41 
42  // Loop over each id for this object
43  for (const auto & id : ids)
44  {
45  // If block materials exist, look if any issue the required guarantee
46  if (warehouse.hasActiveBlockObjects(id))
47  {
48  const std::vector<std::shared_ptr<MaterialBase>> & mats = warehouse.getActiveBlockObjects(id);
49  for (const auto & mat : mats)
50  {
51  const auto & mat_props = mat->getSuppliedItems();
52  if (mat_props.count(prop_name))
53  {
54  auto guarantee_mat = dynamic_cast<GuaranteeProvider *>(mat.get());
55  if (guarantee_mat && !guarantee_mat->hasGuarantee(prop_name, guarantee))
56  {
57  // we found at least one material on the set of block we operate on
58  // that does _not_ provide the requested guarantee
59  return false;
60  }
61  }
62  }
63  }
64  }
65 
66  return true;
67 }
void mooseError(Args &&... args)
Add-on class that provides the functionality to issue guarantees for declared material properties...
BlockRestrictable *const _gc_block_restrict
Access block restrictions of the object with this interface.
virtual const std::set< SubdomainID > & blockIDs() const
virtual bool blockRestricted() const
const MaterialWarehouse & getMaterialWarehouse() const
FEProblemBase *const _gc_feproblem
Reference to the FEProblemBase class.
virtual MooseMesh & mesh() override
virtual bool startedInitialSetup()
const std::set< SubdomainID > & meshSubdomains() const

◆ initialSetup()

void ComputeMultipleCrystalPlasticityStress::initialSetup ( )
overridevirtual

Reimplemented from ComputeFiniteStrainElasticStress.

Definition at line 148 of file ComputeMultipleCrystalPlasticityStress.C.

149 {
150  // get crystal plasticity models
151  std::vector<MaterialName> model_names =
152  getParam<std::vector<MaterialName>>("crystal_plasticity_models");
153 
154  for (unsigned int i = 0; i < _num_models; ++i)
155  {
157  dynamic_cast<CrystalPlasticityStressUpdateBase *>(&getMaterialByName(model_names[i]));
158 
159  if (model)
160  {
161  _models.push_back(model);
162  // TODO: check to make sure that the material model is compatible with this class
163  }
164  else
165  mooseError("Model " + model_names[i] +
166  " is not compatible with ComputeMultipleCrystalPlasticityStress");
167  }
168 
169  // get crystal plasticity eigenstrains
170  std::vector<MaterialName> eigenstrain_names =
171  getParam<std::vector<MaterialName>>("eigenstrain_names");
172 
173  for (unsigned int i = 0; i < _num_eigenstrains; ++i)
174  {
177  &getMaterialByName(eigenstrain_names[i]));
178 
179  if (eigenstrain)
180  _eigenstrains.push_back(eigenstrain);
181  else
182  mooseError("Eigenstrain" + eigenstrain_names[i] +
183  " is not compatible with ComputeMultipleCrystalPlasticityStress");
184  }
185 }
ComputeCrystalPlasticityEigenstrainBase is the base class for computing eigenstrain tensors in crysta...
const PertinentGeochemicalSystem model(database, {"H2O", "H+", "HCO3-", "O2(aq)", "Ca++", ">(s)FeOH", "radius_neg1", "radius_neg1.5"}, {"Calcite"}, {}, {"Calcite_asdf"}, {"CH4(aq)"}, {">(s)FeOCa+"}, "O2(aq)", "e-")
MaterialBase & getMaterialByName(const std::string &name, bool no_warn=false, bool no_dep=false)
std::vector< CrystalPlasticityStressUpdateBase * > _models
The user supplied cyrstal plasticity consititutive models.
void mooseError(Args &&... args) const
std::vector< ComputeCrystalPlasticityEigenstrainBase * > _eigenstrains
The user supplied cyrstal plasticity eigenstrains.

◆ 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 110 of file ComputeMultipleCrystalPlasticityStress.C.

111 {
114 
115  if (_num_eigenstrains)
116  {
117  (*_eigenstrain_deformation_gradient)[_qp].zero();
118  (*_eigenstrain_deformation_gradient)[_qp].addIa(1.0);
119  }
120  else
121  {
122  // set to identity if no eigenstrain is added
125  }
126 
127  _pk2[_qp].zero();
128 
130 
131  _updated_rotation[_qp].zero();
132  _updated_rotation[_qp].addIa(1.0);
133 
134  for (unsigned int i = 0; i < _num_models; ++i)
135  {
136  _models[i]->setQp(_qp);
137  _models[i]->initQpStatefulProperties();
138  }
139 
140  for (unsigned int i = 0; i < _num_eigenstrains; ++i)
141  {
142  _eigenstrains[i]->setQp(_qp);
143  _eigenstrains[i]->initQpStatefulProperties();
144  }
145 }
MaterialProperty< RankTwoTensor > & _total_lagrangian_strain
Lagrangian total strain measure for the entire crystal.
void addIa(const Real &a)
MaterialProperty< RankTwoTensor > & _pk2
Second Piola-Kirchoff stress measure.
MaterialProperty< RankTwoTensor > & _updated_rotation
Tracks the rotation of the crystal during deformation Note: this rotation tensor is not applied to th...
MaterialProperty< RankTwoTensor > & _plastic_deformation_gradient
Plastic deformation gradient RankTwoTensor for the crystal.
std::vector< CrystalPlasticityStressUpdateBase * > _models
The user supplied cyrstal plasticity consititutive models.
std::vector< ComputeCrystalPlasticityEigenstrainBase * > _eigenstrains
The user supplied cyrstal plasticity eigenstrains.

◆ lineSearchUpdate()

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

performs the line search update

Definition at line 655 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by solveStress().

657 {
659  {
660  Real rnorm;
661  Real step = 1.0;
662 
663  do
664  {
665  _pk2[_qp] = _pk2[_qp] - step * dpk2;
666  step /= 2.0;
667  _pk2[_qp] = _pk2[_qp] + step * dpk2;
668 
670  rnorm = _residual_tensor.L2norm();
671  } while (rnorm > rnorm_prev && step > _min_line_search_step_size);
672 
673  // has norm improved or is the step still above minumum search step size?
674  return (rnorm <= rnorm_prev || step > _min_line_search_step_size);
675  }
677  {
678  unsigned int count = 0;
679  Real step_a = 0.0;
680  Real step_b = 1.0;
681  Real step = 1.0;
682  Real s_m = 1000.0;
683  Real rnorm = 1000.0;
684 
686  auto s_b = _residual_tensor.doubleContraction(dpk2);
687  const auto rnorm1 = _residual_tensor.L2norm();
688  _pk2[_qp] = _pk2[_qp] - dpk2;
690  auto s_a = _residual_tensor.doubleContraction(dpk2);
691  const auto rnorm0 = _residual_tensor.L2norm();
692  _pk2[_qp] = _pk2[_qp] + dpk2;
693 
694  if ((rnorm1 / rnorm0) < _line_search_tolerance || s_a * s_b > 0)
695  {
697  return true;
698  }
699 
700  while ((rnorm / rnorm0) > _line_search_tolerance && count < _line_search_max_iterations)
701  {
702  _pk2[_qp] = _pk2[_qp] - step * dpk2;
703  step = 0.5 * (step_b + step_a);
704  _pk2[_qp] = _pk2[_qp] + step * dpk2;
707  rnorm = _residual_tensor.L2norm();
708 
709  if (s_m * s_a < 0.0)
710  {
711  step_b = step;
712  s_b = s_m;
713  }
714  if (s_m * s_b < 0.0)
715  {
716  step_a = step;
717  s_a = s_m;
718  }
719  count++;
720  }
721 
722  // below tolerance and max iterations?
723  return ((rnorm / rnorm0) < _line_search_tolerance && count < _line_search_max_iterations);
724  }
725  else
726  mooseError("Line search method is not provided.");
727 }
Real _line_search_tolerance
Line search bisection method tolerance.
void calculateResidual()
Calculate stress residual as the difference between the stored material property PK2 stress and the e...
enum ComputeMultipleCrystalPlasticityStress::LineSearchMethod _line_search_method
Real _min_line_search_step_size
Minimum line search step size.
MaterialProperty< RankTwoTensor > & _pk2
Second Piola-Kirchoff stress measure.
Real doubleContraction(const RankTwoTensorTempl< Real > &a) const
unsigned int _line_search_max_iterations
Line search bisection method maximum iteration number.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
void mooseError(Args &&... args) const

◆ 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 303 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by updateStress().

305 {
306  cauchy_stress = _elastic_deformation_gradient * _pk2[_qp] *
308 
309  calcTangentModuli(jacobian_mult);
310 
315 
316  // Calculate crystal rotation to track separately
317  RankTwoTensor rot;
319  _updated_rotation[_qp] = rot * _crysrot[_qp];
320 
321  // Calculate the misorientation
322  auto _delta_misorientation_mat = _updated_rotation[_qp] * _updated_rotation_old[_qp].inverse();
323  auto cos_val = (_delta_misorientation_mat.trace() - 1.0) / 2.0;
324  // Mathematically, the values should lie in the range of [-1, 1]. However, there maybe numerical
325  // errors during calculation. Let's guard small numerical errors here and throw an error when the
326  // numerical error is too big, which indicates other issues.
327  if (cos_val > 1.0 || cos_val < -1.0)
328  {
329  if (MooseUtils::absoluteFuzzyEqual(cos_val, -1.0))
330  cos_val = -1.0;
331  else if (MooseUtils::absoluteFuzzyEqual(cos_val, 1.0))
332  cos_val = 1.0;
333  else
334  mooseError("Misorientation value is undefined. This should not happen.");
335  }
336  _misorientation[_qp] = std::acos(cos_val);
337 }
MaterialProperty< RankTwoTensor > & _total_lagrangian_strain
Lagrangian total strain measure for the entire crystal.
static RankTwoTensorTempl Identity()
const MaterialProperty< RankTwoTensor > & _crysrot
Crystal rotation in the original, or reference, configuration as defined by Euler angle arguments in ...
const MaterialProperty< RankTwoTensor > & _deformation_gradient
Total deformation gradient RankTwoTensor for the crystal.
MaterialProperty< RankTwoTensor > & _pk2
Second Piola-Kirchoff stress measure.
MaterialProperty< RankTwoTensor > & _updated_rotation
Tracks the rotation of the crystal during deformation Note: this rotation tensor is not applied to th...
RankTwoTensorTempl< Real > transpose() const
const MaterialProperty< RankTwoTensor > & _updated_rotation_old
MaterialProperty< Real > & _misorientation
Misorientation angle of the crystal during deformation.
void getRUDecompositionRotation(RankTwoTensorTempl< Real > &rot) const
void mooseError(Args &&... args) const
void calcTangentModuli(RankFourTensor &jacobian_mult)
Calculates the tangent moduli for use as a preconditioner, using the elastic or elastic-plastic optio...

◆ 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 271 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by updateStress().

272 {
273  for (unsigned int i = 0; i < _num_models; ++i)
274  _models[i]->setInitialConstitutiveVariableValues();
275 
276  _pk2[_qp] = _pk2_old[_qp];
278 }
const MaterialProperty< RankTwoTensor > & _pk2_old
MaterialProperty< RankTwoTensor > & _pk2
Second Piola-Kirchoff stress measure.
std::vector< CrystalPlasticityStressUpdateBase * > _models
The user supplied cyrstal plasticity consititutive models.
const MaterialProperty< RankTwoTensor > & _plastic_deformation_gradient_old

◆ 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 281 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by updateStress().

282 {
283  for (unsigned int i = 0; i < _num_models; ++i)
284  {
285  _models[i]->setSubstepConstitutiveVariableValues();
286  _models[i]->calculateSlipResistance();
287  }
288 
290 
293  return; // pop back up and take a smaller substep
294 
295  for (unsigned int i = 0; i < _num_models; ++i)
296  _models[i]->updateSubstepConstitutiveVariableValues();
297 
298  // save off the old F^{p} inverse now that have converged on the stress and state variables
300 }
void solveStateVariables()
Solves the internal variables stress as a function of the slip specified by the constitutive model de...
bool _convergence_failed
Flag to check whether convergence is achieved or if substepping is needed.
std::vector< CrystalPlasticityStressUpdateBase * > _models
The user supplied cyrstal plasticity consititutive models.

◆ 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 340 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by solveQp().

341 {
342  unsigned int iteration;
343  bool iter_flag = true;
344 
345  iteration = 0;
346  // Check for slip system resistance update tolerance
347  do
348  {
349  solveStress();
351  return;
352 
354  _inverse_plastic_deformation_grad.inverse(); // the postSoveStress
355 
356  // Update slip system resistance and state variable after the stress has been finalized
357  // We loop through all the models for each calculation
358  // in order to make sure that when coupling appears, the state variables are updated based on
359  // the same components
360  for (unsigned int i = 0; i < _num_models; ++i)
361  _models[i]->cacheStateVariablesBeforeUpdate();
362 
363  for (unsigned int i = 0; i < _num_models; ++i)
364  _models[i]->calculateStateVariableEvolutionRateComponent();
365 
366  for (unsigned int i = 0; i < _num_models; ++i)
367  if (!_models[i]->updateStateVariables())
368  _convergence_failed = true;
369 
370  for (unsigned int i = 0; i < _num_models; ++i)
371  _models[i]->calculateSlipResistance();
372 
374  return;
375 
376  for (unsigned int i = 0; i < _num_models; ++i)
377  {
378  // iter_flag = false, stop iteration if all values are converged and good to go
379  // iter_flag = true, continue iteration if any value is not converged
380  if (!_models[i]->areConstitutiveStateVariablesConverged())
381  {
382  iter_flag = true;
383  break;
384  }
385  else
386  iter_flag = false; // iter_flag = false, stop iteration only when all models returns true
387  }
388 
389  if (iter_flag)
390  {
392  mooseWarning("ComputeMultipleCrystalPlasticityStress: State variables (or the system "
393  "resistance) did not converge at element ",
394  _current_elem->id(),
395  " and qp ",
396  _qp,
397  "\n");
398  }
399  iteration++;
400  } while (iter_flag && iteration < _maxiterg);
401 
402  if (iteration == _maxiterg)
403  {
405  mooseWarning(
406  "ComputeMultipleCrystalPlasticityStress: Hardness Integration error. Reached the "
407  "maximum number of iterations to solve for the state variables at element ",
408  _current_elem->id(),
409  " and qp ",
410  _qp,
411  "\n");
412 
413  _convergence_failed = true;
414  }
415 }
RankTwoTensorTempl< Real > inverse() const
bool _convergence_failed
Flag to check whether convergence is achieved or if substepping is needed.
MaterialProperty< RankTwoTensor > & _plastic_deformation_gradient
Plastic deformation gradient RankTwoTensor for the crystal.
std::vector< CrystalPlasticityStressUpdateBase * > _models
The user supplied cyrstal plasticity consititutive models.
unsigned int _maxiterg
Maximum number of iterations for internal variable update.
void mooseWarning(Args &&... args) const
void solveStress()
solves for stress, updates plastic deformation gradient.
const bool _print_convergence_message
Flag to print to console warning messages on stress, constitutive model convergence.

◆ solveStress()

void ComputeMultipleCrystalPlasticityStress::solveStress ( )
protected

solves for stress, updates plastic deformation gradient.

Definition at line 418 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by solveStateVariables().

419 {
420  unsigned int iteration = 0;
421  RankTwoTensor dpk2;
422  Real rnorm, rnorm0, rnorm_prev;
423 
424  // Calculate stress residual
427  {
429  mooseWarning("ComputeMultipleCrystalPlasticityStress: the slip increment exceeds tolerance "
430  "at element ",
431  _current_elem->id(),
432  " and Gauss point ",
433  _qp);
434 
435  return;
436  }
437 
438  rnorm = _residual_tensor.L2norm();
439  rnorm0 = rnorm;
440 
441  // Check for stress residual tolerance; different from user object version which
442  // compares the absolute tolerance of only the original rnorm value
443  while (rnorm > _rtol * rnorm0 && rnorm > _abs_tol && iteration < _maxiter)
444  {
445  // Calculate stress increment
446  dpk2 = -_jacobian.invSymm() * _residual_tensor;
447  _pk2[_qp] = _pk2[_qp] + dpk2;
448 
450 
452  {
454  mooseWarning("ComputeMultipleCrystalPlasticityStress: the slip increment exceeds tolerance "
455  "at element ",
456  _current_elem->id(),
457  " and Gauss point ",
458  _qp);
459 
460  return;
461  }
462 
463  rnorm_prev = rnorm;
464  rnorm = _residual_tensor.L2norm();
465 
466  if (_use_line_search && rnorm > rnorm_prev && !lineSearchUpdate(rnorm_prev, dpk2))
467  {
469  mooseWarning("ComputeMultipleCrystalPlasticityStress: Failed with line search");
470 
471  _convergence_failed = true;
472  return;
473  }
474 
475  if (_use_line_search)
476  rnorm = _residual_tensor.L2norm();
477 
478  iteration++;
479  }
480 
481  if (iteration >= _maxiter)
482  {
484  mooseWarning("ComputeMultipleCrystalPlasticityStress: Stress Integration error rmax = ",
485  rnorm,
486  " and the tolerance is ",
487  _rtol * rnorm0,
488  " when the rnorm0 value is ",
489  rnorm0,
490  " for element ",
491  _current_elem->id(),
492  " and qp ",
493  _qp);
494 
495  _convergence_failed = true;
496  }
497 }
bool _convergence_failed
Flag to check whether convergence is achieved or if substepping is needed.
unsigned int _maxiter
Maximum number of iterations for stress update.
bool lineSearchUpdate(const Real &rnorm_prev, const RankTwoTensor &dpk2)
performs the line search update
MaterialProperty< RankTwoTensor > & _pk2
Second Piola-Kirchoff stress measure.
Real _abs_tol
Stress residual equation absolute tolerance.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
void mooseWarning(Args &&... args) const
void calculateResidualAndJacobian()
Calls the residual and jacobian functions used in the stress update algorithm.
const bool _print_convergence_message
Flag to print to console warning messages on stress, constitutive model convergence.
RankFourTensorTempl< T > invSymm() const
Real _rtol
Stress residual equation relative tolerance.

◆ 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 203 of file ComputeMultipleCrystalPlasticityStress.C.

Referenced by computeQpStress().

205 {
206  // Does not support face/boundary material property calculation
207  if (isBoundaryMaterial())
208  return;
209 
210  // Initialize substepping variables
211  unsigned int substep_iter = 1;
212  unsigned int num_substep = 1;
213 
217 
219 
220  // Loop through all models and calculate the schmid tensor for the current state of the crystal
221  // lattice
222  // Not sure if we should pass in the updated or the original rotation here
223  // If not, then we should not need to compute the flow direction every iteration here
224  for (unsigned int i = 0; i < _num_models; ++i)
225  _models[i]->calculateFlowDirection(_crysrot[_qp]);
226 
227  do
228  {
229  _convergence_failed = false;
230  preSolveQp();
231 
232  _substep_dt = _dt / num_substep;
233  for (unsigned int i = 0; i < _num_models; ++i)
234  _models[i]->setSubstepDt(_substep_dt);
235 
236  // calculate F^{eigen} only when we have eigenstrain
239  if (_num_eigenstrains)
241 
242  for (unsigned int istep = 0; istep < num_substep; ++istep)
243  {
245  (static_cast<Real>(istep) + 1) / num_substep * _delta_deformation_gradient;
247 
248  solveQp();
249 
251  {
253  mooseWarning(
254  "The crystal plasticity constitutive model has failed to converge. Increasing "
255  "the number of substeps.");
256 
257  substep_iter++;
258  num_substep *= 2;
259  break;
260  }
261  }
262 
263  if (substep_iter > _max_substep_iter && _convergence_failed)
264  mooseException("ComputeMultipleCrystalPlasticityStress: Constitutive failure");
265  } while (_convergence_failed);
266 
267  postSolveQp(cauchy_stress, jacobian_mult);
268 }
void postSolveQp(RankTwoTensor &stress_new, RankFourTensor &jacobian_mult)
Save the final stress and internal variable values after the iterative solve.
RankTwoTensor _temporary_deformation_gradient
Helper deformation gradient tensor variables used in iterative solve.
void calculateEigenstrainDeformationGrad()
Calculates the deformation gradient due to eigenstrain.
bool _convergence_failed
Flag to check whether convergence is achieved or if substepping is needed.
virtual bool isBoundaryMaterial() const override
const MaterialProperty< RankTwoTensor > & _crysrot
Crystal rotation in the original, or reference, configuration as defined by Euler angle arguments in ...
const MaterialProperty< RankTwoTensor > & _deformation_gradient
Total deformation gradient RankTwoTensor for the crystal.
void addIa(const Real &a)
unsigned int _max_substep_iter
Maximum number of substep iterations.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
std::vector< CrystalPlasticityStressUpdateBase * > _models
The user supplied cyrstal plasticity consititutive models.
const MaterialProperty< RankTwoTensor > & _deformation_gradient_old
void mooseWarning(Args &&... args) const
void solveQp()
Solve the stress and internal state variables (e.g.
const bool _print_convergence_message
Flag to print to console warning messages on stress, constitutive model convergence.
RankTwoTensor _delta_deformation_gradient
Used for substepping; Uniformly divides the increment in deformation gradient.
void preSolveQp()
Reset the PK2 stress and the inverse deformation gradient to old values and provide an interface for ...

◆ validParams()

InputParameters ComputeMultipleCrystalPlasticityStress::validParams ( )
static

Definition at line 20 of file ComputeMultipleCrystalPlasticityStress.C.

21 {
23 
24  params.addClassDescription(
25  "Crystal Plasticity base class: handles the Newton iteration over the stress residual and "
26  "calculates the Jacobian based on constitutive laws from multiple material classes "
27  "that are inherited from CrystalPlasticityStressUpdateBase");
28  params.addParam<std::string>(
29  "base_name",
30  "Optional parameter that allows the user to define multiple mechanics material systems on "
31  "the same block, i.e. for multiple phases");
32 
33  params.addRequiredParam<std::vector<MaterialName>>(
34  "crystal_plasticity_models",
35  "The material objects to use to calculate crystal plasticity stress and strains.");
36  params.addParam<std::vector<MaterialName>>(
37  "eigenstrain_names", {}, "The material objects to calculate eigenstrains.");
38  params.addParam<MooseEnum>("tan_mod_type",
39  MooseEnum("exact none", "none"),
40  "Type of tangent moduli for preconditioner: default elastic");
41  params.addParam<Real>("rtol", 1e-6, "Constitutive stress residual relative tolerance");
42  params.addParam<Real>("abs_tol", 1e-6, "Constitutive stress residual absolute tolerance");
43  params.addParam<unsigned int>("maxiter", 100, "Maximum number of iterations for stress update");
44  params.addParam<unsigned int>(
45  "maxiter_state_variable", 100, "Maximum number of iterations for state variable update");
46  params.addParam<unsigned int>(
47  "maximum_substep_iteration", 1, "Maximum number of substep iteration");
48  params.addParam<bool>("use_line_search", false, "Use line search in constitutive update");
49  params.addParam<Real>("min_line_search_step_size", 0.01, "Minimum line search step size");
50  params.addParam<Real>("line_search_tol", 0.5, "Line search bisection method tolerance");
51  params.addParam<unsigned int>(
52  "line_search_maxiter", 20, "Line search bisection method maximum number of iteration");
53  params.addParam<MooseEnum>("line_search_method",
54  MooseEnum("CUT_HALF BISECTION", "CUT_HALF"),
55  "The method used in line search");
56  params.addParam<bool>(
57  "print_state_variable_convergence_error_messages",
58  false,
59  "Whether or not to print warning messages from the crystal plasticity specific convergence "
60  "checks on the stress measure and general constitutive model quantinties.");
61  return params;
62 }
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
void addRequiredParam(const std::string &name, const std::string &doc_string)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
void addClassDescription(const std::string &doc_string)

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 230 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 216 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 233 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 238 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

◆ _elastic_strain

MaterialProperty<RankTwoTensor>& ComputeGeneralStressBase::_elastic_strain
protectedinherited

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

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

◆ _inverse_plastic_deformation_grad

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_inverse_plastic_deformation_grad
protected

◆ _inverse_plastic_deformation_grad_old

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_inverse_plastic_deformation_grad_old
protected

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

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

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

◆ _misorientation

MaterialProperty<Real>& ComputeMultipleCrystalPlasticityStress::_misorientation
protected

Misorientation angle of the crystal during deformation.

Definition at line 210 of file ComputeMultipleCrystalPlasticityStress.h.

Referenced by postSolveQp().

◆ _models

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

◆ _num_eigenstrains

const unsigned ComputeMultipleCrystalPlasticityStress::_num_eigenstrains
protected

◆ _num_models

const unsigned ComputeMultipleCrystalPlasticityStress::_num_models
protected

◆ _pk2

MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_pk2
protected

◆ _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 227 of file ComputeMultipleCrystalPlasticityStress.h.

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

◆ _residual_tensor

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_residual_tensor
protected

◆ _rotation_increment

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainElasticStress::_rotation_increment
protectedinherited

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

◆ _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 219 of file ComputeMultipleCrystalPlasticityStress.h.

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

◆ _temporary_deformation_gradient_old

RankTwoTensor ComputeMultipleCrystalPlasticityStress::_temporary_deformation_gradient_old
protected

Definition at line 234 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().

◆ _updated_rotation_old

const MaterialProperty<RankTwoTensor>& ComputeMultipleCrystalPlasticityStress::_updated_rotation_old
protected

Definition at line 205 of file ComputeMultipleCrystalPlasticityStress.h.

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