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
CrystalPlasticityKalidindiUpdate Class Reference

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

Inheritance diagram for CrystalPlasticityKalidindiUpdate:
[legend]

Public Types

enum  ConstantTypeEnum { ConstantTypeEnum::NONE, ConstantTypeEnum::ELEMENT, ConstantTypeEnum::SUBDOMAIN }
 
enum  TEST_TYPE
 
typedef DataFileName DataFileParameterType
 

Public Member Functions

 CrystalPlasticityKalidindiUpdate (const InputParameters &parameters)
 
void setQp (const unsigned int &qp)
 Sets the value of the global variable _qp for inheriting classes. More...
 
void setSubstepDt (const Real &substep_dt)
 Sets the value of the _substep_dt for inheriting classes. More...
 
virtual void setMaterialVectorSize ()
 
virtual void getSlipSystems ()
 A helper method to read in plane normal and direction vectors from a file and to normalize the vectors. More...
 
void transformHexagonalMillerBravaisSlipSystems (const MooseUtils::DelimitedFileReader &reader)
 A helper method to transform the Miller-Bravais 4-index notation for HCP crystals into a a 3-index Cartesian representation, using the convention a$_1$ = x of axis alignment in the basal plane. More...
 
void calculateFlowDirection (const RankTwoTensor &crysrot)
 Computes the Schmid tensor (m x n) for the original (reference) crystal lattice orientation for each glide slip system. More...
 
void calculateShearStress (const RankTwoTensor &pk2, const RankTwoTensor &inverse_eigenstrain_deformation_grad, const unsigned int &num_eigenstrains)
 Computes the shear stess for each slip system. More...
 
virtual void calculateTotalPlasticDeformationGradientDerivative (RankFourTensor &dfpinvdpk2, const RankTwoTensor &inverse_plastic_deformation_grad_old, const RankTwoTensor &inverse_eigenstrain_deformation_grad_old, const unsigned int &num_eigenstrains)
 Calculates the total value of ${d{F}^P^{-1}}{d{PK2}}$ and is intended to be an overwritten helper method for inheriting classes with multiple constitutive dislocation slip mechanisms, e.g. More...
 
void calculateSchmidTensor (const unsigned int &number_dislocation_systems, const std::vector< RealVectorValue > &plane_normal_vector, const std::vector< RealVectorValue > &direction_vector, std::vector< RankTwoTensor > &schmid_tensor, const RankTwoTensor &crysrot)
 A helper method to rotate the a direction and plane normal system set into the local crystal llatice orientation as defined by the crystal rotation tensor from the Elasticity tensor class. More...
 
void sortCrossSlipFamilies ()
 A helper method to sort the slip systems of a crystal into cross slip families based on common slip directions. More...
 
unsigned int identifyCrossSlipFamily (const unsigned int index)
 A helper method for inherting classes to identify to which cross slip family vector a particular slip system index belongs after the slip systems have been sorted. More...
 
virtual void calculateSlipResistance ()
 
virtual bool isConstitutiveStateVariableConverged (const std::vector< Real > &current_var, const std::vector< Real > &var_before_update, const std::vector< Real > &previous_substep_var, const Real &tolerance)
 Check if a typical state variable, e.g. More...
 
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
 
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
 
bool ghostable () const override final
 
virtual void resolveOptionalProperties () override
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty ()
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty ()
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty ()
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialPropertyByName (const std::string &prop_name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialPropertyByName (const std::string &prop_name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialPropertyByName (const std::string &prop_name)
 
const MaterialProperty< T > & getZeroMaterialProperty (Ts... args)
 
const MaterialProperty< T > & getZeroMaterialProperty (Ts... args)
 
const MaterialProperty< T > & getZeroMaterialProperty (Ts... args)
 
virtual void initStatefulProperties (unsigned int n_points)
 
virtual bool isInterfaceMaterial ()
 
virtual void computePropertiesAtQp (unsigned int qp)
 
const MaterialProperty< T > & getZeroMaterialPropertyByName (Ts... args)
 
virtual const std::set< std::string > & getRequestedItems () override
 
virtual const std::set< std::string > & getSuppliedItems () override
 
const std::set< unsigned int > & getSuppliedPropIDs ()
 
void checkStatefulSanity () const
 
std::set< OutputName > getOutputs ()
 
bool hasStatefulProperties () const
 
void setFaceInfo (const FaceInfo &fi)
 
void setActiveProperties (const std::unordered_set< unsigned int > &needed_props)
 
bool forceStatefulInit () const
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
MooseAppgetMooseApp () const
 
const std::string & type () const
 
virtual const std::string & name () const
 
std::string typeAndName () const
 
std::string errorPrefix (const std::string &error_type) const
 
void callMooseError (std::string msg, const bool with_prefix) const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
const InputParametersparameters () const
 
MooseObjectName uniqueName () const
 
const T & getParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const T * queryParam (const std::string &name) const
 
const T & getRenamedParam (const std::string &old_name, const std::string &new_name) const
 
getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &nm) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void mooseError (Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
const std::vector< SubdomainName > & blocks () const
 
unsigned int numBlocks () const
 
virtual const std::set< SubdomainID > & blockIDs () const
 
unsigned int blocksMaxDimension () const
 
bool hasBlocks (const SubdomainName &name) const
 
bool hasBlocks (const std::vector< SubdomainName > &names) const
 
bool hasBlocks (const std::set< SubdomainName > &names) const
 
bool hasBlocks (SubdomainID id) const
 
bool hasBlocks (const std::vector< SubdomainID > &ids) const
 
bool hasBlocks (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::vector< SubdomainID > &ids) const
 
bool hasBlockMaterialProperty (const std::string &prop_name)
 
const std::set< SubdomainID > & meshBlockIDs () const
 
virtual bool blockRestricted () const
 
virtual void checkVariable (const MooseVariableFieldBase &variable) const
 
virtual const std::set< BoundaryID > & boundaryIDs () const
 
const std::vector< BoundaryName > & boundaryNames () const
 
unsigned int numBoundaryIDs () const
 
bool hasBoundary (const BoundaryName &name) const
 
bool hasBoundary (const std::vector< BoundaryName > &names) const
 
bool hasBoundary (const BoundaryID &id) const
 
bool hasBoundary (const std::vector< BoundaryID > &ids, TEST_TYPE type=ALL) const
 
bool hasBoundary (const std::set< BoundaryID > &ids, TEST_TYPE type=ALL) const
 
bool isBoundarySubset (const std::set< BoundaryID > &ids) const
 
bool isBoundarySubset (const std::vector< BoundaryID > &ids) const
 
bool hasBoundaryMaterialProperty (const std::string &prop_name) const
 
virtual bool boundaryRestricted () const
 
const std::set< BoundaryID > & meshBoundaryIDs () const
 
virtual bool checkVariableBoundaryIntegrity () const
 
virtual void initialSetup ()
 
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 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
 
UserObjectName getUserObjectName (const std::string &param_name) const
 
const T & getUserObject (const std::string &param_name, bool is_dependency=true) const
 
const T & getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const UserObjectgetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectgetUserObjectBaseByName (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 ADVariableValuegetADDefaultValue (const std::string &var_name) const
 
const ADVectorVariableValuegetADDefaultVectorValue (const std::string &var_name) const
 
const ADVariableGradientgetADDefaultGradient () const
 
const ADVectorVariableGradientgetADDefaultVectorGradient () const
 
const ADVariableSecondgetADDefaultSecond () const
 
const ADVectorVariableCurlgetADDefaultCurl () const
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const std::string &prop_name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const std::string &prop_name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const std::string &prop_name)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const std::string &prop_name)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const std::string &prop_name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const std::string &prop_name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const std::string &prop_name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const std::string &prop_name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > > getBlockMaterialProperty (const MaterialPropertyName &name)
 
std::set< 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)
 
std::unordered_map< SubdomainID, std::vector< MaterialBase *> > buildRequiredMaterials (bool allow_stateful=true)
 
void statefulPropertiesAllowed (bool)
 
bool getMaterialPropertyCalled () const
 
const GenericMaterialProperty< T, is_ad > & getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)
 
const GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
MaterialProperty< T > & declarePropertyByName (const std::string &prop_name)
 
MaterialProperty< T > & declarePropertyByName (const std::string &prop_name)
 
MaterialProperty< T > & declareProperty (const std::string &name)
 
MaterialProperty< T > & declareProperty (const std::string &name)
 
ADMaterialProperty< T > & declareADPropertyByName (const std::string &prop_name)
 
ADMaterialProperty< T > & declareADPropertyByName (const std::string &prop_name)
 
ADMaterialProperty< T > & declareADProperty (const std::string &name)
 
ADMaterialProperty< T > & declareADProperty (const std::string &name)
 
auto & declareGenericProperty (const std::string &prop_name)
 
auto & declareGenericProperty (const std::string &prop_name)
 
GenericMaterialProperty< T, is_ad > & declareGenericPropertyByName (const std::string &prop_name)
 
GenericMaterialProperty< T, is_ad > & declareGenericPropertyByName (const std::string &prop_name)
 
const DistributiongetDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const T & getDistributionByName (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const T & getDistributionByName (const std::string &name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
const 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 hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
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
 
virtual void resetQpProperties () final
 Retained as empty methods to avoid a warning from Material.C in framework. These methods are unused in all inheriting classes and should not be overwritten. More...
 
virtual void resetProperties () final
 

Static Public Member Functions

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

Public Attributes

 ALL
 
 ANY
 
const ConsoleStream _console
 

Static Public Attributes

static constexpr PropertyValue::id_type default_property_id
 
static constexpr PropertyValue::id_type zero_property_id
 

Protected Types

enum  CrystalLatticeType { CrystalLatticeType::BCC, CrystalLatticeType::FCC, CrystalLatticeType::HCP }
 
enum  QP_Data_Type
 

Protected Member Functions

virtual void initQpStatefulProperties () override
 initializes the stateful properties such as stress, plastic deformation gradient, slip system resistances, etc. More...
 
virtual void setInitialConstitutiveVariableValues () override
 Sets the value of the current and previous substep iteration slip system resistance to the old value at the start of the PK2 stress convergence while loop. More...
 
virtual void setSubstepConstitutiveVariableValues () override
 Sets the current slip system resistance value to the previous substep value. More...
 
virtual void updateSubstepConstitutiveVariableValues () override
 Stores the current value of the slip system resistance into a separate material property in case substepping is needed. More...
 
virtual bool calculateSlipRate () override
 This virtual method is called to calculate the slip system slip increment based on the constitutive model defined in the child class. More...
 
virtual void calculateEquivalentSlipIncrement (RankTwoTensor &) override
 
virtual void calculateConstitutiveSlipDerivative (std::vector< Real > &dslip_dtau) override
 This virtual method is called to find the derivative of the slip increment with respect to the applied shear stress on the slip system based on the constiutive model defined in the child class. More...
 
virtual void cacheStateVariablesBeforeUpdate () override
 Finalizes the values of the state variables and slip system resistance for the current timestep after convergence has been reached. More...
 
virtual void calculateStateVariableEvolutionRateComponent () override
 Following the Constitutive model for slip system resistance as given in Kalidindi, S.R., C.A. More...
 
virtual bool updateStateVariables () override
 Finalizes the values of the state variables and slip system resistance for the current timestep after convergence has been reached. More...
 
virtual bool areConstitutiveStateVariablesConverged () override
 Determines if all the state variables have converged. More...
 
virtual void checkMaterialProperty (const std::string &name, const unsigned int state) override
 
virtual const MaterialDatamaterialData () const override
 
virtual MaterialDatamaterialData () override
 
virtual const QBase & qRule () const override
 
virtual void computeQpProperties ()
 
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 ()
 
virtual bool hasBlockMaterialPropertyHelper (const std::string &prop_name)
 
void initializeBlockRestrictable (const MooseObject *moose_object)
 
Moose::CoordinateSystemType getBlockCoordSystem ()
 
bool hasBoundaryMaterialPropertyHelper (const std::string &prop_name) const
 
bool isCoupledScalar (const std::string &var_name, unsigned int i=0) const
 
unsigned int coupledScalarComponents (const std::string &var_name) const
 
unsigned int coupledScalar (const std::string &var_name, unsigned int comp=0) const
 
libMesh::Order coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const
 
const 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
 
virtual void addUserObjectDependencyHelper (const UserObject &) const
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &) const
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &) const
 
T & declareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< T > declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const T & getRestartableData (const std::string &data_name) const
 
T & declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
T & declareRecoverableData (const std::string &data_name, Args &&... args)
 
T & declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
T & declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
 
std::string restartableName (const std::string &data_name) const
 
std::string deduceFunctorName (const std::string &name) const
 
const Moose::Functor< T > & getFunctor (const std::string &name)
 
const Moose::Functor< T > & getFunctor (const std::string &name, THREAD_ID tid)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem, THREAD_ID tid)
 
bool isFunctor (const std::string &name) const
 
bool isFunctor (const std::string &name, const SubProblem &subproblem) const
 
Moose::ElemArg makeElemArg (const Elem *elem, bool correct_skewnewss=false) const
 
void checkFunctorSupportsSideIntegration (const std::string &name, bool qp_integration)
 
void flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
 
InvalidSolutionID registerInvalidSolutionInternal (const std::string &message, const bool warning) const
 
virtual void coupledCallback (const std::string &, bool) const
 
virtual bool isCoupled (const std::string &var_name, unsigned int i=0) const
 
virtual bool isCoupledConstant (const std::string &var_name) const
 
unsigned int coupledComponents (const std::string &var_name) const
 
VariableName coupledName (const std::string &var_name, unsigned int comp=0) const
 
std::vector< VariableName > coupledNames (const std::string &var_name) const
 
virtual unsigned int coupled (const std::string &var_name, unsigned int comp=0) const
 
std::vector< unsigned 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 T & coupledNodalValue (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 T & coupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const 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 ()
 
bool checkVar (const std::string &var_name, unsigned int comp=0, unsigned int comp_bound=0) const
 
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)
 
const T * getVarHelper (const std::string &var_name, unsigned int comp) const
 
T * getVarHelper (const std::string &var_name, unsigned int comp)
 
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
 
void validateExecutionerType (const std::string &name, const std::string &fn_name) const
 
std::vector< T > coupledVectorHelper (const std::string &var_name, const Func &func) const
 
void markMatPropRequested (const std::string &)
 
MaterialPropertyName getMaterialPropertyName (const std::string &name) const
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 

Protected Attributes

std::vector< Real_hb
 Slip system interaction matrix used to calculate the hardening contributions from the self and latent slip systems, from Kalidindi et al (1992). More...
 
std::vector< Real_slip_resistance_increment
 Increment of increased resistance for each slip system. More...
 
std::vector< Real_previous_substep_slip_resistance
 Stores the values of the slip system resistance from the previous substep In classes which use dislocation densities, analogous dislocation density substep vectors will be required. More...
 
std::vector< Real_slip_resistance_before_update
 Caches the value of the current slip system resistance immediately prior to the update of the slip system resistance, and is used to calculate the the slip system resistance increment for the current substep (or step if only one substep is taken) for the convergence check tolerance comparison. More...
 
const bool _include_twinning_in_Lp
 Flag to include the total twin volume fraction in the plastic velocity gradient calculation, per Kalidindi IJP (2001). More...
 
const MaterialProperty< Real > *const _twin_volume_fraction_total
 User-defined material property name for the total volume fraction of twins in a twinning propagation constitutive model, when this class is used in conjunction with the twinning propagation model. More...
 
const std::string _base_name
 Base name prepended to all material property names to allow for multi-material systems. More...
 
enum CrystalPlasticityStressUpdateBase::CrystalLatticeType _crystal_lattice_type
 
const std::vector< Real_unit_cell_dimension
 
const unsigned int _number_slip_systems
 Maximum number of active slip systems for the crystalline material being modeled. More...
 
std::string _slip_sys_file_name
 File should contain slip plane normal and direction. More...
 
Real _rel_state_var_tol
 Internal variable update equation tolerance. More...
 
Real _slip_incr_tol
 Slip increment tolerance. More...
 
Real _resistance_tol
 Tolerance for change in slip system resistance over an increment. More...
 
Real _zero_tol
 Residual tolerance when variable value is zero. Default 1e-12. More...
 
MaterialProperty< std::vector< Real > > & _slip_increment
 Current slip increment material property. More...
 
MaterialProperty< std::vector< Real > > & _tau
 Resolved shear stress on each slip system. More...
 
const bool _print_convergence_message
 Flag to print to console warning messages on stress, constitutive model convergence. More...
 
Real _substep_dt
 Substepping time step value used within the inheriting constitutive models. More...
 
std::vector< std::vector< unsigned int > > _cross_slip_familes
 Sorted slip system indices into cross slip family groups. More...
 
bool _calculate_cross_slip
 Flag to run the cross slip calculations if cross slip numbers are specified. More...
 
 CURR
 
 PREV
 
bool _bnd
 
bool _neighbor
 
const MooseArray< Point > & _q_point
 
const QBase *const & _qrule
 
const MooseArray< Real > & _JxW
 
const Elem *const & _current_elem
 
const SubdomainID_current_subdomain_id
 
const unsigned int_current_side
 
const ConstantTypeEnum _constant_option
 
SubProblem_subproblem
 
FEProblemBase_fe_problem
 
THREAD_ID _tid
 
Assembly_assembly
 
unsigned int _qp
 
const MooseArray< Real > & _coord
 
const MooseArray< Point > & _normals
 
MooseMesh_mesh
 
const Moose::CoordinateSystemType_coord_sys
 
std::set< std::string > _requested_props
 
std::set< std::string > _supplied_props
 
std::set< unsigned int_supplied_prop_ids
 
std::unordered_set< unsigned int_active_prop_ids
 
const bool _compute
 
std::unordered_map< unsigned int, unsigned int_props_to_min_states
 
std::vector< unsigned int_displacements
 
bool _has_stateful_property
 
bool _overrides_init_stateful_props
 
const FaceInfo_face_info
 
const bool & _enabled
 
MooseApp_app
 
const std::string _type
 
const std::string _name
 
const InputParameters_pars
 
Factory_factory
 
ActionFactory_action_factory
 
const MaterialData_blk_material_data
 
const ExecFlagEnum_execute_enum
 
const ExecFlagType_current_execute_flag
 
FEProblemBase_sc_fe_problem
 
const THREAD_ID _sc_tid
 
const Real_real_zero
 
const VariableValue_scalar_zero
 
const Point & _point_zero
 
const InputParameters_ti_params
 
FEProblemBase_ti_feproblem
 
bool _is_implicit
 
Real_t
 
const Real_t_old
 
int_t_step
 
Real_dt
 
Real_dt_old
 
bool _is_transient
 
MooseApp_restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase_mci_feproblem
 
GeometricSearchData_geometric_search_data
 
bool _requires_geometric_search
 
const InputParameters_c_parameters
 
const std::string & _c_name
 
const std::string & _c_type
 
FEProblemBase_c_fe_problem
 
const SystemBase *const _c_sys
 
std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> > _coupled_vars
 
std::vector< MooseVariableFieldBase *> _coupled_moose_vars
 
std::vector< MooseVariable *> _coupled_standard_moose_vars
 
std::vector< VectorMooseVariable *> _coupled_vector_moose_vars
 
std::vector< ArrayMooseVariable *> _coupled_array_moose_vars
 
std::vector< MooseVariableFV< Real > *> _coupled_standard_fv_moose_vars
 
std::vector< MooseLinearVariableFV< Real > *> _coupled_standard_linear_fv_moose_vars
 
const std::unordered_map< std::string, std::string > & _new_to_deprecated_coupled_vars
 
bool _c_nodal
 
bool _c_is_implicit
 
const bool _c_allow_element_to_nodal_coupling
 
THREAD_ID _c_tid
 
std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > > _default_value
 
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADReal > > > _ad_default_value
 
std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > > _default_vector_value
 
std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > > _default_array_value
 
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > > _ad_default_vector_value
 
VariableValue _default_value_zero
 
VariableGradient _default_gradient
 
MooseArray< ADRealVectorValue_ad_default_gradient
 
MooseArray< ADRealTensorValue_ad_default_vector_gradient
 
VariableSecond _default_second
 
MooseArray< ADRealTensorValue_ad_default_second
 
MooseArray< ADRealVectorValue_ad_default_curl
 
const VariableValue_zero
 
const VariablePhiValue_phi_zero
 
const MooseArray< ADReal > & _ad_zero
 
const VariableGradient_grad_zero
 
const MooseArray< ADRealVectorValue > & _ad_grad_zero
 
const VariablePhiGradient_grad_phi_zero
 
const VariableSecond_second_zero
 
const MooseArray< ADRealTensorValue > & _ad_second_zero
 
const VariablePhiSecond_second_phi_zero
 
const VectorVariableValue_vector_zero
 
const VectorVariableCurl_vector_curl_zero
 
VectorVariableValue _default_vector_value_zero
 
VectorVariableGradient _default_vector_gradient
 
VectorVariableCurl _default_vector_curl
 
VectorVariableDivergence _default_div
 
ArrayVariableValue _default_array_value_zero
 
ArrayVariableGradient _default_array_gradient
 
bool _coupleable_neighbor
 
const InputParameters_mi_params
 
const std::string _mi_name
 
const MooseObjectName _mi_moose_object_name
 
FEProblemBase_mi_feproblem
 
SubProblem_mi_subproblem
 
const THREAD_ID _mi_tid
 
const Moose::MaterialDataType _material_data_type
 
MaterialData_material_data
 
bool _stateful_allowed
 
bool _get_material_property_called
 
std::vector< std::unique_ptr< PropertyValue > > _default_properties
 
std::unordered_set< unsigned int_material_property_dependencies
 
const MaterialPropertyName _get_suffix
 
const bool _use_interpolated_state
 
const Parallel::Communicator & _communicator
 
const Real _r
 Varibles used in the Kalidindi 1992 slip system resistance constiutive model. More...
 
const Real _h
 
const Real _tau_sat
 
const Real _gss_a
 
const Real _ao
 
const Real _xm
 
const Real _gss_initial
 
const Real _number_cross_slip_directions
 Parameters to characterize the cross slip behavior of the crystal. More...
 
const Real _number_cross_slip_planes
 
MaterialProperty< std::vector< Real > > & _slip_resistance
 Slip system resistance. More...
 
const MaterialProperty< std::vector< Real > > & _slip_resistance_old
 
std::vector< RealVectorValue_slip_direction
 Slip system direction and normal and associated Schmid tensors. More...
 
std::vector< RealVectorValue_slip_plane_normal
 
MaterialProperty< std::vector< RankTwoTensor > > & _flow_direction
 

Static Protected Attributes

static const std::string _interpolated_old
 
static const std::string _interpolated_older
 

Detailed Description

CrystalPlasticityKalidindiUpdate 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 interative predictor-corrector algorithm. Backward Euler integration rule is used for the rate equations.

Definition at line 24 of file CrystalPlasticityKalidindiUpdate.h.

Member Enumeration Documentation

◆ CrystalLatticeType

Enumerator
BCC 
FCC 
HCP 

Definition at line 192 of file CrystalPlasticityStressUpdateBase.h.

192 { BCC, FCC, HCP } _crystal_lattice_type;
enum CrystalPlasticityStressUpdateBase::CrystalLatticeType _crystal_lattice_type

Constructor & Destructor Documentation

◆ CrystalPlasticityKalidindiUpdate()

CrystalPlasticityKalidindiUpdate::CrystalPlasticityKalidindiUpdate ( const InputParameters parameters)

Definition at line 35 of file CrystalPlasticityKalidindiUpdate.C.

38  // Constitutive values
39  _r(getParam<Real>("r")),
40  _h(getParam<Real>("h")),
41  _tau_sat(getParam<Real>("t_sat")),
42  _gss_a(getParam<Real>("gss_a")),
43  _ao(getParam<Real>("ao")),
44  _xm(getParam<Real>("xm")),
45  _gss_initial(getParam<Real>("gss_initial")),
46 
47  // resize vectors used in the consititutive slip hardening
50 
51  // resize local caching vectors used for substepping
54 
55  // Twinning contributions, if used
56  _include_twinning_in_Lp(parameters.isParamValid("total_twin_volume_fraction")),
58  ? &getMaterialPropertyOld<Real>("total_twin_volume_fraction")
59  : nullptr)
60 {
61 }
const Real _r
Varibles used in the Kalidindi 1992 slip system resistance constiutive model.
CrystalPlasticityStressUpdateBase(const InputParameters &parameters)
std::vector< Real > _previous_substep_slip_resistance
Stores the values of the slip system resistance from the previous substep In classes which use disloc...
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
const MaterialProperty< Real > *const _twin_volume_fraction_total
User-defined material property name for the total volume fraction of twins in a twinning propagation ...
std::vector< Real > _slip_resistance_increment
Increment of increased resistance for each slip system.
const InputParameters & parameters() const
std::vector< Real > _hb
Slip system interaction matrix used to calculate the hardening contributions from the self and latent...
std::vector< Real > _slip_resistance_before_update
Caches the value of the current slip system resistance immediately prior to the update of the slip sy...
const bool _include_twinning_in_Lp
Flag to include the total twin volume fraction in the plastic velocity gradient calculation, per Kalidindi IJP (2001).
bool isParamValid(const std::string &name) const

Member Function Documentation

◆ areConstitutiveStateVariablesConverged()

bool CrystalPlasticityKalidindiUpdate::areConstitutiveStateVariablesConverged ( )
overrideprotectedvirtual

Determines if all the state variables have converged.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 141 of file CrystalPlasticityKalidindiUpdate.C.

142 {
147 }
std::vector< Real > _previous_substep_slip_resistance
Stores the values of the slip system resistance from the previous substep In classes which use disloc...
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
virtual bool isConstitutiveStateVariableConverged(const std::vector< Real > &current_var, const std::vector< Real > &var_before_update, const std::vector< Real > &previous_substep_var, const Real &tolerance)
Check if a typical state variable, e.g.
Real _resistance_tol
Tolerance for change in slip system resistance over an increment.
std::vector< Real > _slip_resistance_before_update
Caches the value of the current slip system resistance immediately prior to the update of the slip sy...

◆ cacheStateVariablesBeforeUpdate()

void CrystalPlasticityKalidindiUpdate::cacheStateVariablesBeforeUpdate ( )
overrideprotectedvirtual

Finalizes the values of the state variables and slip system resistance for the current timestep after convergence has been reached.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 157 of file CrystalPlasticityKalidindiUpdate.C.

158 {
160 }
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
std::vector< Real > _slip_resistance_before_update
Caches the value of the current slip system resistance immediately prior to the update of the slip sy...

◆ calculateConstitutiveSlipDerivative()

void CrystalPlasticityKalidindiUpdate::calculateConstitutiveSlipDerivative ( std::vector< Real > &  )
overrideprotectedvirtual

This virtual method is called to find the derivative of the slip increment with respect to the applied shear stress on the slip system based on the constiutive model defined in the child class.

This method must be overwritten in the child class.

Implements CrystalPlasticityStressUpdateBase.

Definition at line 126 of file CrystalPlasticityKalidindiUpdate.C.

128 {
129  for (const auto i : make_range(_number_slip_systems))
130  {
132  dslip_dtau[i] = 0.0;
133  else
134  dslip_dtau[i] = _ao / _xm *
135  std::pow(std::abs(_tau[_qp][i] / _slip_resistance[_qp][i]), 1.0 / _xm - 1.0) /
136  _slip_resistance[_qp][i];
137  }
138 }
bool absoluteFuzzyEqual(const T &var1, const T2 &var2, const T3 &tol=libMesh::TOLERANCE *libMesh::TOLERANCE)
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
MaterialProperty< std::vector< Real > > & _tau
Resolved shear stress on each slip system.
IntRange< T > make_range(T beg, T end)
MooseUnits pow(const MooseUnits &, int)

◆ calculateEquivalentSlipIncrement()

void CrystalPlasticityKalidindiUpdate::calculateEquivalentSlipIncrement ( RankTwoTensor equivalent_slip_increment)
overrideprotectedvirtual

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 112 of file CrystalPlasticityKalidindiUpdate.C.

114 {
116  {
117  for (const auto i : make_range(_number_slip_systems))
118  equivalent_slip_increment += (1.0 - (*_twin_volume_fraction_total)[_qp]) *
120  }
121  else // if no twinning volume fraction material property supplied, use base class
123 }
virtual void calculateEquivalentSlipIncrement(RankTwoTensor &)
MaterialProperty< std::vector< Real > > & _slip_increment
Current slip increment material property.
Real _substep_dt
Substepping time step value used within the inheriting constitutive models.
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
MaterialProperty< std::vector< RankTwoTensor > > & _flow_direction
IntRange< T > make_range(T beg, T end)
const bool _include_twinning_in_Lp
Flag to include the total twin volume fraction in the plastic velocity gradient calculation, per Kalidindi IJP (2001).

◆ calculateFlowDirection()

void CrystalPlasticityStressUpdateBase::calculateFlowDirection ( const RankTwoTensor crysrot)
inherited

Computes the Schmid tensor (m x n) for the original (reference) crystal lattice orientation for each glide slip system.

Definition at line 369 of file CrystalPlasticityStressUpdateBase.C.

370 {
373 }
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
MaterialProperty< std::vector< RankTwoTensor > > & _flow_direction
void calculateSchmidTensor(const unsigned int &number_dislocation_systems, const std::vector< RealVectorValue > &plane_normal_vector, const std::vector< RealVectorValue > &direction_vector, std::vector< RankTwoTensor > &schmid_tensor, const RankTwoTensor &crysrot)
A helper method to rotate the a direction and plane normal system set into the local crystal llatice ...
std::vector< RealVectorValue > _slip_direction
Slip system direction and normal and associated Schmid tensors.

◆ calculateSchmidTensor()

void CrystalPlasticityStressUpdateBase::calculateSchmidTensor ( const unsigned int number_dislocation_systems,
const std::vector< RealVectorValue > &  plane_normal_vector,
const std::vector< RealVectorValue > &  direction_vector,
std::vector< RankTwoTensor > &  schmid_tensor,
const RankTwoTensor crysrot 
)
inherited

A helper method to rotate the a direction and plane normal system set into the local crystal llatice orientation as defined by the crystal rotation tensor from the Elasticity tensor class.

Definition at line 376 of file CrystalPlasticityStressUpdateBase.C.

Referenced by CrystalPlasticityStressUpdateBase::calculateFlowDirection().

382 {
383  std::vector<RealVectorValue> local_direction_vector, local_plane_normal;
384  local_direction_vector.resize(number_slip_systems);
385  local_plane_normal.resize(number_slip_systems);
386 
387  // Update slip direction and normal with crystal orientation
388  for (const auto i : make_range(_number_slip_systems))
389  {
390  local_direction_vector[i].zero();
391  local_plane_normal[i].zero();
392 
393  for (const auto j : make_range(LIBMESH_DIM))
394  for (const auto k : make_range(LIBMESH_DIM))
395  {
396  local_direction_vector[i](j) =
397  local_direction_vector[i](j) + crysrot(j, k) * direction_vector[i](k);
398 
399  local_plane_normal[i](j) =
400  local_plane_normal[i](j) + crysrot(j, k) * plane_normal_vector[i](k);
401  }
402 
403  // Calculate Schmid tensor
404  for (const auto j : make_range(LIBMESH_DIM))
405  for (const auto k : make_range(LIBMESH_DIM))
406  {
407  schmid_tensor[i](j, k) = local_direction_vector[i](j) * local_plane_normal[i](k);
408  }
409  }
410 }
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
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:130

◆ calculateShearStress()

void CrystalPlasticityStressUpdateBase::calculateShearStress ( const RankTwoTensor pk2,
const RankTwoTensor inverse_eigenstrain_deformation_grad,
const unsigned int num_eigenstrains 
)
inherited

Computes the shear stess for each slip system.

Definition at line 413 of file CrystalPlasticityStressUpdateBase.C.

417 {
418  if (!num_eigenstrains)
419  {
420  for (const auto i : make_range(_number_slip_systems))
422 
423  return;
424  }
425 
426  RankTwoTensor eigenstrain_deformation_grad = inverse_eigenstrain_deformation_grad.inverse();
427  for (const auto i : make_range(_number_slip_systems))
428  {
429  // compute PK2_hat using deformation gradient
430  RankTwoTensor pk2_hat = eigenstrain_deformation_grad.det() *
431  eigenstrain_deformation_grad.transpose() * pk2 *
432  inverse_eigenstrain_deformation_grad.transpose();
433  _tau[_qp][i] = pk2_hat.doubleContraction(_flow_direction[_qp][i]);
434  }
435 }
RankTwoTensorTempl< Real > inverse() const
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
MaterialProperty< std::vector< RankTwoTensor > > & _flow_direction
Real doubleContraction(const RankTwoTensorTempl< Real > &a) const
RankTwoTensorTempl< Real > transpose() const
MaterialProperty< std::vector< Real > > & _tau
Resolved shear stress on each slip system.
IntRange< T > make_range(T beg, T end)

◆ calculateSlipRate()

bool CrystalPlasticityKalidindiUpdate::calculateSlipRate ( )
overrideprotectedvirtual

This virtual method is called to calculate the slip system slip increment based on the constitutive model defined in the child class.

This method must be overwritten in the child class.

Implements CrystalPlasticityStressUpdateBase.

Definition at line 90 of file CrystalPlasticityKalidindiUpdate.C.

91 {
92  for (const auto i : make_range(_number_slip_systems))
93  {
94  _slip_increment[_qp][i] =
95  _ao * std::pow(std::abs(_tau[_qp][i] / _slip_resistance[_qp][i]), 1.0 / _xm);
96  if (_tau[_qp][i] < 0.0)
97  _slip_increment[_qp][i] *= -1.0;
98 
99  if (std::abs(_slip_increment[_qp][i]) * _substep_dt > _slip_incr_tol)
100  {
102  mooseWarning("Maximum allowable slip increment exceeded ",
103  std::abs(_slip_increment[_qp][i]) * _substep_dt);
104 
105  return false;
106  }
107  }
108  return true;
109 }
MaterialProperty< std::vector< Real > > & _slip_increment
Current slip increment material property.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
void mooseWarning(Args &&... args) const
Real _substep_dt
Substepping time step value used within the inheriting constitutive models.
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
MaterialProperty< std::vector< Real > > & _tau
Resolved shear stress on each slip system.
IntRange< T > make_range(T beg, T end)
const bool _print_convergence_message
Flag to print to console warning messages on stress, constitutive model convergence.
MooseUnits pow(const MooseUnits &, int)

◆ calculateSlipResistance()

virtual void CrystalPlasticityStressUpdateBase::calculateSlipResistance ( )
inlinevirtualinherited

◆ calculateStateVariableEvolutionRateComponent()

void CrystalPlasticityKalidindiUpdate::calculateStateVariableEvolutionRateComponent ( )
overrideprotectedvirtual

Following the Constitutive model for slip system resistance as given in Kalidindi, S.R., C.A.

Bronkhorst, and L. Anand. Crystallographic texture evolution in bulk deformation processing of FCC metals. Journal of the Mechanics and Physics of Solids 40, no. 3 (1992): 537-569. Eqns 40 - 43. The slip system resistant increment is calculated as $ g = | q^{ } h^{} |$ and a convergence check is performed on the slip system resistance increment

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 163 of file CrystalPlasticityKalidindiUpdate.C.

164 {
165  for (const auto i : make_range(_number_slip_systems))
166  {
167  // Clear out increment from the previous iteration
169 
170  _hb[i] = _h * std::pow(std::abs(1.0 - _slip_resistance[_qp][i] / _tau_sat), _gss_a);
171  const Real hsign = 1.0 - _slip_resistance[_qp][i] / _tau_sat;
172  if (hsign < 0.0)
173  _hb[i] *= -1.0;
174  }
175 
176  for (const auto i : make_range(_number_slip_systems))
177  {
178  for (const auto j : make_range(_number_slip_systems))
179  {
180  unsigned int iplane, jplane;
181  iplane = i / 3;
182  jplane = j / 3;
183 
184  if (iplane == jplane) // self vs. latent hardening
186  std::abs(_slip_increment[_qp][j]) * _hb[j]; // q_{ab} = 1.0 for self hardening
187  else
189  std::abs(_slip_increment[_qp][j]) * _r * _hb[j]; // latent hardenign
190  }
191  }
192 }
const Real _r
Varibles used in the Kalidindi 1992 slip system resistance constiutive model.
MaterialProperty< std::vector< Real > > & _slip_increment
Current slip increment material property.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
std::vector< Real > _slip_resistance_increment
Increment of increased resistance for each slip system.
IntRange< T > make_range(T beg, T end)
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
std::vector< Real > _hb
Slip system interaction matrix used to calculate the hardening contributions from the self and latent...
MooseUnits pow(const MooseUnits &, int)

◆ calculateTotalPlasticDeformationGradientDerivative()

void CrystalPlasticityStressUpdateBase::calculateTotalPlasticDeformationGradientDerivative ( RankFourTensor dfpinvdpk2,
const RankTwoTensor inverse_plastic_deformation_grad_old,
const RankTwoTensor inverse_eigenstrain_deformation_grad_old,
const unsigned int num_eigenstrains 
)
virtualinherited

Calculates the total value of ${d{F}^P^{-1}}{d{PK2}}$ and is intended to be an overwritten helper method for inheriting classes with multiple constitutive dislocation slip mechanisms, e.g.

glide and twinning, $ {d{F}^P^{-1}}{d{PK2}_i}$

Definition at line 438 of file CrystalPlasticityStressUpdateBase.C.

443 {
444  std::vector<Real> dslip_dtau(_number_slip_systems, 0.0);
445  std::vector<RankTwoTensor> dtaudpk2(_number_slip_systems);
446  std::vector<RankTwoTensor> dfpinvdslip(_number_slip_systems);
447 
449 
450  for (const auto j : make_range(_number_slip_systems))
451  {
452  if (num_eigenstrains)
453  {
454  RankTwoTensor eigenstrain_deformation_grad_old =
455  inverse_eigenstrain_deformation_grad_old.inverse();
456  dtaudpk2[j] = eigenstrain_deformation_grad_old.det() * eigenstrain_deformation_grad_old *
457  _flow_direction[_qp][j] * inverse_eigenstrain_deformation_grad_old;
458  }
459  else
460  dtaudpk2[j] = _flow_direction[_qp][j];
461  dfpinvdslip[j] = -inverse_plastic_deformation_grad_old * _flow_direction[_qp][j];
462  dfpinvdpk2 += (dfpinvdslip[j] * dslip_dtau[j] * _substep_dt).outerProduct(dtaudpk2[j]);
463  }
464 }
RankTwoTensorTempl< Real > inverse() const
Real _substep_dt
Substepping time step value used within the inheriting constitutive models.
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
virtual void calculateConstitutiveSlipDerivative(std::vector< Real > &)=0
This virtual method is called to find the derivative of the slip increment with respect to the applie...
MaterialProperty< std::vector< RankTwoTensor > > & _flow_direction
IntRange< T > make_range(T beg, T end)
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")

◆ getSlipSystems()

void CrystalPlasticityStressUpdateBase::getSlipSystems ( )
virtualinherited

A helper method to read in plane normal and direction vectors from a file and to normalize the vectors.

This method is abstracted to allow for reuse in inheriting classes with multiple plane normal and direction vector pairs.

Definition at line 141 of file CrystalPlasticityStressUpdateBase.C.

Referenced by CrystalPlasticityStressUpdateBase::CrystalPlasticityStressUpdateBase().

142 {
143  bool orthonormal_error = false;
144 
145  // read in the slip system data from auxiliary text file
148  _reader.read();
149 
150  // check the size of the input
151  if (_reader.getData().size() != _number_slip_systems)
152  paramError(
153  "number_slip_systems",
154  "The number of rows in the slip system file should match the number of slip system.");
155 
156  for (const auto i : make_range(_number_slip_systems))
157  {
158  // initialize to zero
159  _slip_direction[i].zero();
160  _slip_plane_normal[i].zero();
161  }
162 
167  {
168  for (const auto i : make_range(_number_slip_systems))
169  {
170  // directly grab the raw data and scale it by the unit cell dimension
171  for (const auto j : index_range(_reader.getData(i)))
172  {
173  if (j < LIBMESH_DIM)
174  _slip_plane_normal[i](j) = _reader.getData(i)[j] / _unit_cell_dimension[j];
175  else
176  _slip_direction[i](j - LIBMESH_DIM) =
177  _reader.getData(i)[j] * _unit_cell_dimension[j - LIBMESH_DIM];
178  }
179  }
180  }
181 
182  for (const auto i : make_range(_number_slip_systems))
183  {
184  // normalize
185  _slip_plane_normal[i] /= _slip_plane_normal[i].norm();
186  _slip_direction[i] /= _slip_direction[i].norm();
187 
189  {
190  const auto magnitude = _slip_plane_normal[i] * _slip_direction[i];
191  if (std::abs(magnitude) > libMesh::TOLERANCE)
192  {
193  orthonormal_error = true;
194  break;
195  }
196  }
197  }
198 
199  if (orthonormal_error)
200  mooseError("CrystalPlasticityStressUpdateBase Error: The slip system file contains a slip "
201  "direction and plane normal pair that are not orthonormal in the Cartesian "
202  "coordinate system.");
203 }
std::string _slip_sys_file_name
File should contain slip plane normal and direction.
void transformHexagonalMillerBravaisSlipSystems(const MooseUtils::DelimitedFileReader &reader)
A helper method to transform the Miller-Bravais 4-index notation for HCP crystals into a a 3-index Ca...
static constexpr Real TOLERANCE
enum CrystalPlasticityStressUpdateBase::CrystalLatticeType _crystal_lattice_type
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
void paramError(const std::string &param, Args... args) const
IntRange< T > make_range(T beg, T end)
void mooseError(Args &&... args) const
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
auto index_range(const T &sizable)
std::vector< RealVectorValue > _slip_direction
Slip system direction and normal and associated Schmid tensors.

◆ identifyCrossSlipFamily()

unsigned int CrystalPlasticityStressUpdateBase::identifyCrossSlipFamily ( const unsigned int  index)
inherited

A helper method for inherting classes to identify to which cross slip family vector a particular slip system index belongs after the slip systems have been sorted.

Returns the integer value of the identified cross slip system family for the outer vector created in sortCrossSlipFamilies.

Definition at line 357 of file CrystalPlasticityStressUpdateBase.C.

358 {
359  for (unsigned int i = 0; i < _number_cross_slip_directions; ++i)
360  for (unsigned int j = 0; j < _number_cross_slip_planes; ++j)
361  if (_cross_slip_familes[i][j] == index)
362  return i;
363 
364  // Should never reach this statement
365  mooseError("The supplied slip system index is not among the slip system families sorted.");
366 }
std::vector< std::vector< unsigned int > > _cross_slip_familes
Sorted slip system indices into cross slip family groups.
const Real _number_cross_slip_directions
Parameters to characterize the cross slip behavior of the crystal.
void mooseError(Args &&... args) const
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")

◆ initQpStatefulProperties()

void CrystalPlasticityKalidindiUpdate::initQpStatefulProperties ( )
overrideprotectedvirtual

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

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 64 of file CrystalPlasticityKalidindiUpdate.C.

65 {
67  for (const auto i : make_range(_number_slip_systems))
68  {
70  _slip_increment[_qp][i] = 0.0;
71  }
72 }
virtual void initQpStatefulProperties() override
initializes the stateful properties such as PK2 stress, resolved shear stress, plastic deformation gr...
MaterialProperty< std::vector< Real > > & _slip_increment
Current slip increment material property.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
IntRange< T > make_range(T beg, T end)

◆ isConstitutiveStateVariableConverged()

bool CrystalPlasticityStressUpdateBase::isConstitutiveStateVariableConverged ( const std::vector< Real > &  current_var,
const std::vector< Real > &  var_before_update,
const std::vector< Real > &  previous_substep_var,
const Real tolerance 
)
virtualinherited

Check if a typical state variable, e.g.

defect density, has converged by comparing the change in the values over the iteration period.

Definition at line 488 of file CrystalPlasticityStressUpdateBase.C.

Referenced by CrystalPlasticityTwinningKalidindiUpdate::areConstitutiveStateVariablesConverged(), areConstitutiveStateVariablesConverged(), and CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::areConstitutiveStateVariablesConverged().

493 {
494  // sometimes the state variable size may not equal to the number of slip systems
495  unsigned int sz = current_var.size();
496  mooseAssert(current_var.size() == sz, "Current variable size does not match");
497  mooseAssert(var_before_update.size() == sz, "Variable before update size does not match");
498  mooseAssert(previous_substep_var.size() == sz, "Previous substep variable size does not match");
499 
500  bool is_converged = true;
501 
502  Real diff_val = 0.0;
503  Real abs_prev_substep_val = 0.0;
504  for (const auto i : make_range(sz))
505  {
506  diff_val = std::abs(var_before_update[i] - current_var[i]);
507  abs_prev_substep_val = std::abs(previous_substep_var[i]);
508 
509  // set to false if the state variable is not converged
510  if (abs_prev_substep_val < _zero_tol && diff_val > _zero_tol)
511  is_converged = false;
512  else if (abs_prev_substep_val > _zero_tol && diff_val > tolerance * abs_prev_substep_val)
513  is_converged = false;
514  }
515  return is_converged;
516 }
Real _zero_tol
Residual tolerance when variable value is zero. Default 1e-12.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
IntRange< T > make_range(T beg, T end)

◆ resetProperties()

virtual void CrystalPlasticityStressUpdateBase::resetProperties ( )
inlinefinalvirtualinherited

Reimplemented from Material.

Definition at line 32 of file CrystalPlasticityStressUpdateBase.h.

32 {}

◆ resetQpProperties()

virtual void CrystalPlasticityStressUpdateBase::resetQpProperties ( )
inlinefinalvirtualinherited

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

Reimplemented from Material.

Definition at line 31 of file CrystalPlasticityStressUpdateBase.h.

31 {}

◆ setInitialConstitutiveVariableValues()

void CrystalPlasticityKalidindiUpdate::setInitialConstitutiveVariableValues ( )
overrideprotectedvirtual

Sets the value of the current and previous substep iteration slip system resistance to the old value at the start of the PK2 stress convergence while loop.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 75 of file CrystalPlasticityKalidindiUpdate.C.

76 {
77  // Would also set old dislocation densities here if included in this model
80 }
std::vector< Real > _previous_substep_slip_resistance
Stores the values of the slip system resistance from the previous substep In classes which use disloc...
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
const MaterialProperty< std::vector< Real > > & _slip_resistance_old

◆ setMaterialVectorSize()

void CrystalPlasticityStressUpdateBase::setMaterialVectorSize ( )
virtualinherited

Reimplemented in CrystalPlasticityTwinningKalidindiUpdate, and CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.

Definition at line 126 of file CrystalPlasticityStressUpdateBase.C.

Referenced by CrystalPlasticityStressUpdateBase::initQpStatefulProperties(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::setMaterialVectorSize(), and CrystalPlasticityTwinningKalidindiUpdate::setMaterialVectorSize().

127 {
130  for (const auto i : make_range(_number_slip_systems))
131  {
132  _flow_direction[_qp][i].zero();
133  _tau[_qp][i] = 0.0;
134  }
135 
138 }
MaterialProperty< std::vector< Real > > & _slip_increment
Current slip increment material property.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
virtual void resize(const std::size_t size) override final
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
MaterialProperty< std::vector< RankTwoTensor > > & _flow_direction
MaterialProperty< std::vector< Real > > & _tau
Resolved shear stress on each slip system.
IntRange< T > make_range(T beg, T end)

◆ setQp()

void CrystalPlasticityStressUpdateBase::setQp ( const unsigned int qp)
inherited

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

Definition at line 476 of file CrystalPlasticityStressUpdateBase.C.

477 {
478  _qp = qp;
479 }
unsigned int _qp

◆ setSubstepConstitutiveVariableValues()

void CrystalPlasticityKalidindiUpdate::setSubstepConstitutiveVariableValues ( )
overrideprotectedvirtual

Sets the current slip system resistance value to the previous substep value.

In cases where only one substep is taken (or when the first) substep is taken, this method just sets the current value to the old slip system resistance value again.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 83 of file CrystalPlasticityKalidindiUpdate.C.

84 {
85  // Would also set substepped dislocation densities here if included in this model
87 }
std::vector< Real > _previous_substep_slip_resistance
Stores the values of the slip system resistance from the previous substep In classes which use disloc...
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp

◆ setSubstepDt()

void CrystalPlasticityStressUpdateBase::setSubstepDt ( const Real substep_dt)
inherited

Sets the value of the _substep_dt for inheriting classes.

Definition at line 482 of file CrystalPlasticityStressUpdateBase.C.

483 {
484  _substep_dt = substep_dt;
485 }
Real _substep_dt
Substepping time step value used within the inheriting constitutive models.

◆ sortCrossSlipFamilies()

void CrystalPlasticityStressUpdateBase::sortCrossSlipFamilies ( )
inherited

A helper method to sort the slip systems of a crystal into cross slip families based on common slip directions.

This method determines if slip directions are parallel, and stores the index of the slip systems from getSlipSystems (the same index is used for the applied shear stress and the internal state variables) in a vector of vectors, where the outer vector separates the individual slip system families and the inner vector stories the indices of the slip systems within a single family. This vector of vectors can then be used in the inheriting classes to calculate cross slip interactions. The values of number_cross_slip_directions and number_cross_slip_planes must be set to use this sorting method.

Definition at line 290 of file CrystalPlasticityStressUpdateBase.C.

Referenced by CrystalPlasticityStressUpdateBase::CrystalPlasticityStressUpdateBase().

291 {
293  {
294  _cross_slip_familes.resize(0);
295  return;
296  }
297 
298  // If cross slip does occur, then set up the system of vectors for the families
300  // and set the first index of each inner vector
301  for (unsigned int i = 0; i < _number_cross_slip_directions; ++i)
302  _cross_slip_familes[i].resize(1);
303 
304  // Sort the index of the slip system based vectors into separte families
305  unsigned int family_counter = 1;
306  _cross_slip_familes[0][0] = 0;
307 
308  for (unsigned int i = 1; i < _number_slip_systems; ++i)
309  {
310  for (unsigned int j = 0; j < family_counter; ++j)
311  {
312  // check to see if the slip system direction i matches any of the existing slip directions
313  // First calculate the dot product
314  Real dot_product = 0.0;
315  for (const auto k : make_range(Moose::dim))
316  {
317  unsigned int check_family_index = _cross_slip_familes[j][0];
318  dot_product += std::abs(_slip_direction[check_family_index](k) - _slip_direction[i](k));
319  }
320  // Then check if the dot product is one, if yes, add to family and break
321  if (MooseUtils::absoluteFuzzyEqual(dot_product, 0.0))
322  {
323  _cross_slip_familes[j].push_back(i);
325  mooseError(
326  "Exceeded the number of cross slip planes allowed in a single cross slip family");
327 
328  break; // exit the loop over the exisiting cross slip families and move to the next slip
329  // direction
330  }
331  // The slip system in question does not belong to an existing family
332  else if (j == (family_counter - 1) && !MooseUtils::absoluteFuzzyEqual(dot_product, 0.0))
333  {
334  if (family_counter > _number_cross_slip_directions)
335  mooseError("Exceeds the number of cross slip directions specified for this material");
336 
337  _cross_slip_familes[family_counter][0] = i;
338  family_counter++;
339  break;
340  }
341  }
342  }
343 
345  {
346  mooseWarning("Checking the slip system ordering now:");
347  for (unsigned int i = 0; i < _number_cross_slip_directions; ++i)
348  {
349  Moose::out << "In cross slip family " << i << std::endl;
350  for (unsigned int j = 0; j < _number_cross_slip_planes; ++j)
351  Moose::out << " is the slip direction number " << _cross_slip_familes[i][j] << std::endl;
352  }
353  }
354 }
bool absoluteFuzzyEqual(const T &var1, const T2 &var2, const T3 &tol=libMesh::TOLERANCE *libMesh::TOLERANCE)
std::vector< std::vector< unsigned int > > _cross_slip_familes
Sorted slip system indices into cross slip family groups.
static constexpr std::size_t dim
void mooseWarning(Args &&... args) const
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
const Real _number_cross_slip_directions
Parameters to characterize the cross slip behavior of the crystal.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
IntRange< T > make_range(T beg, T end)
void mooseError(Args &&... args) const
const bool _print_convergence_message
Flag to print to console warning messages on stress, constitutive model convergence.
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
static const std::string k
Definition: NS.h:130
std::vector< RealVectorValue > _slip_direction
Slip system direction and normal and associated Schmid tensors.

◆ transformHexagonalMillerBravaisSlipSystems()

void CrystalPlasticityStressUpdateBase::transformHexagonalMillerBravaisSlipSystems ( const MooseUtils::DelimitedFileReader reader)
inherited

A helper method to transform the Miller-Bravais 4-index notation for HCP crystals into a a 3-index Cartesian representation, using the convention a$_1$ = x of axis alignment in the basal plane.

Definition at line 206 of file CrystalPlasticityStressUpdateBase.C.

Referenced by CrystalPlasticityStressUpdateBase::getSlipSystems().

208 {
209  const unsigned int miller_bravais_indices = 4;
210  RealVectorValue temporary_slip_direction, temporary_slip_plane;
211  // temporary_slip_plane.resize(LIBMESH_DIM);
212  // temporary_slip_direction.resize(LIBMESH_DIM);
213 
216  mooseError("CrystalPlasticityStressUpdateBase Error: The specified unit cell dimensions are "
217  "not consistent with expectations for "
218  "HCP crystal hexagonal lattices.");
219  else if (reader.getData(0).size() != miller_bravais_indices * 2)
220  mooseError("CrystalPlasticityStressUpdateBase Error: The number of entries in the first row of "
221  "the slip system file is not consistent with the expectations for the 4-index "
222  "Miller-Bravais assumption for HCP crystals. This file should represent both the "
223  "slip plane normal and the slip direction with 4-indices each.");
224 
225  // set up the tranformation matrices
226  RankTwoTensor transform_matrix;
227  transform_matrix.zero();
228  transform_matrix(0, 0) = 1.0 / _unit_cell_dimension[0];
229  transform_matrix(1, 0) = 1.0 / (_unit_cell_dimension[0] * std::sqrt(3.0));
230  transform_matrix(1, 1) = 2.0 / (_unit_cell_dimension[0] * std::sqrt(3.0));
231  transform_matrix(2, 2) = 1.0 / (_unit_cell_dimension[2]);
232 
233  for (const auto i : make_range(_number_slip_systems))
234  {
235  // read in raw data from file and store in the temporary vectors
236  for (const auto j : index_range(reader.getData(i)))
237  {
238  // Check that the slip plane normal indices of the basal plane sum to zero for consistency
239  Real basal_pl_sum = 0.0;
240  for (const auto k : make_range(LIBMESH_DIM))
241  basal_pl_sum += reader.getData(i)[k];
242 
243  if (basal_pl_sum > _zero_tol)
244  mooseError(
245  "CrystalPlasticityStressUpdateBase Error: The specified HCP basal plane Miller-Bravais "
246  "indices do not sum to zero. Check the values supplied in the associated text file.");
247 
248  // Check that the slip direction indices of the basal plane sum to zero for consistency
249  Real basal_dir_sum = 0.0;
250  for (const auto k : make_range(miller_bravais_indices, miller_bravais_indices + LIBMESH_DIM))
251  basal_dir_sum += reader.getData(i)[k];
252 
253  if (basal_dir_sum > _zero_tol)
254  mooseError("CrystalPlasticityStressUpdateBase Error: The specified HCP slip direction "
255  "Miller-Bravais indices in the basal plane (U, V, and T) do not sum to zero "
256  "within the user specified tolerance (try loosing zero_tol if using the default "
257  "value). Check the values supplied in the associated text file.");
258 
259  if (j < miller_bravais_indices)
260  {
261  // Planes are directly copied over, per a_1 = x convention used here:
262  // Store the first two indices for the basal plane, (h and k), and drop
263  // the redundant third basal plane index (i)
264  if (j < 2)
265  temporary_slip_plane(j) = reader.getData(i)[j];
266  // Store the c-axis index as the third entry in the orthorombic index convention
267  else if (j == 3)
268  temporary_slip_plane(j - 1) = reader.getData(i)[j];
269  }
270  else
271  {
272  const auto direction_j = j - miller_bravais_indices;
273  // Store the first two indices for the slip direction in the basal plane,
274  //(U, V), and drop the redundant third basal plane index (T)
275  if (direction_j < 2)
276  temporary_slip_direction(direction_j) = reader.getData(i)[j];
277  // Store the c-axis index as the third entry in the orthorombic index convention
278  else if (direction_j == 3)
279  temporary_slip_direction(direction_j - 1) = reader.getData(i)[j];
280  }
281  }
282 
283  // perform transformation calculation
284  _slip_direction[i] = transform_matrix * temporary_slip_direction;
285  _slip_plane_normal[i] = transform_matrix * temporary_slip_plane;
286  }
287 }
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
Real _zero_tol
Residual tolerance when variable value is zero. Default 1e-12.
const std::vector< std::vector< T > > & getData() const
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
IntRange< T > make_range(T beg, T end)
void mooseError(Args &&... args) const
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
static const std::string k
Definition: NS.h:130
auto index_range(const T &sizable)
std::vector< RealVectorValue > _slip_direction
Slip system direction and normal and associated Schmid tensors.

◆ updateStateVariables()

bool CrystalPlasticityKalidindiUpdate::updateStateVariables ( )
overrideprotectedvirtual

Finalizes the values of the state variables and slip system resistance for the current timestep after convergence has been reached.

Implements CrystalPlasticityStressUpdateBase.

Definition at line 195 of file CrystalPlasticityKalidindiUpdate.C.

196 {
197  // Now perform the check to see if the slip system should be updated
198  for (const auto i : make_range(_number_slip_systems))
199  {
203  else
204  _slip_resistance[_qp][i] =
206 
207  if (_slip_resistance[_qp][i] < 0.0)
208  return false;
209  }
210  return true;
211 }
std::vector< Real > _previous_substep_slip_resistance
Stores the values of the slip system resistance from the previous substep In classes which use disloc...
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
Real _substep_dt
Substepping time step value used within the inheriting constitutive models.
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
Real _zero_tol
Residual tolerance when variable value is zero. Default 1e-12.
std::vector< Real > _slip_resistance_increment
Increment of increased resistance for each slip system.
IntRange< T > make_range(T beg, T end)

◆ updateSubstepConstitutiveVariableValues()

void CrystalPlasticityKalidindiUpdate::updateSubstepConstitutiveVariableValues ( )
overrideprotectedvirtual

Stores the current value of the slip system resistance into a separate material property in case substepping is needed.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 150 of file CrystalPlasticityKalidindiUpdate.C.

151 {
152  // Would also set substepped dislocation densities here if included in this model
154 }
std::vector< Real > _previous_substep_slip_resistance
Stores the values of the slip system resistance from the previous substep In classes which use disloc...
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp

◆ validParams()

InputParameters CrystalPlasticityKalidindiUpdate::validParams ( )
static

Definition at line 16 of file CrystalPlasticityKalidindiUpdate.C.

17 {
19  params.addClassDescription("Kalidindi version of homogeneous crystal plasticity.");
20  params.addParam<Real>("r", 1.0, "Latent hardening coefficient");
21  params.addParam<Real>("h", 541.5, "hardening constants");
22  params.addParam<Real>("t_sat", 109.8, "saturated slip system strength");
23  params.addParam<Real>("gss_a", 2.5, "coefficient for hardening");
24  params.addParam<Real>("ao", 0.001, "slip rate coefficient");
25  params.addParam<Real>("xm", 0.1, "exponent for slip rate");
26  params.addParam<Real>("gss_initial", 60.8, "initial lattice friction strength of the material");
27 
28  params.addParam<MaterialPropertyName>(
29  "total_twin_volume_fraction",
30  "Total twin volume fraction, if twinning is considered in the simulation");
31 
32  return params;
33 }
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, 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

◆ _ao

const Real CrystalPlasticityKalidindiUpdate::_ao
protected

◆ _base_name

const std::string CrystalPlasticityStressUpdateBase::_base_name
protectedinherited

Base name prepended to all material property names to allow for multi-material systems.

Definition at line 190 of file CrystalPlasticityStressUpdateBase.h.

◆ _calculate_cross_slip

bool CrystalPlasticityStressUpdateBase::_calculate_cross_slip
protectedinherited

Flag to run the cross slip calculations if cross slip numbers are specified.

Definition at line 243 of file CrystalPlasticityStressUpdateBase.h.

Referenced by CrystalPlasticityStressUpdateBase::CrystalPlasticityStressUpdateBase().

◆ _cross_slip_familes

std::vector<std::vector<unsigned int> > CrystalPlasticityStressUpdateBase::_cross_slip_familes
protectedinherited

Sorted slip system indices into cross slip family groups.

Definition at line 240 of file CrystalPlasticityStressUpdateBase.h.

Referenced by CrystalPlasticityStressUpdateBase::identifyCrossSlipFamily(), and CrystalPlasticityStressUpdateBase::sortCrossSlipFamilies().

◆ _crystal_lattice_type

enum CrystalPlasticityStressUpdateBase::CrystalLatticeType CrystalPlasticityStressUpdateBase::_crystal_lattice_type
protectedinherited

◆ _flow_direction

MaterialProperty<std::vector<RankTwoTensor> >& CrystalPlasticityStressUpdateBase::_flow_direction
protectedinherited

◆ _gss_a

const Real CrystalPlasticityKalidindiUpdate::_gss_a
protected

◆ _gss_initial

const Real CrystalPlasticityKalidindiUpdate::_gss_initial
protected

Definition at line 99 of file CrystalPlasticityKalidindiUpdate.h.

Referenced by initQpStatefulProperties().

◆ _h

const Real CrystalPlasticityKalidindiUpdate::_h
protected

◆ _hb

std::vector<Real> CrystalPlasticityKalidindiUpdate::_hb
protected

Slip system interaction matrix used to calculate the hardening contributions from the self and latent slip systems, from Kalidindi et al (1992).

Definition at line 106 of file CrystalPlasticityKalidindiUpdate.h.

Referenced by calculateStateVariableEvolutionRateComponent().

◆ _include_twinning_in_Lp

const bool CrystalPlasticityKalidindiUpdate::_include_twinning_in_Lp
protected

Flag to include the total twin volume fraction in the plastic velocity gradient calculation, per Kalidindi IJP (2001).

Definition at line 132 of file CrystalPlasticityKalidindiUpdate.h.

Referenced by calculateEquivalentSlipIncrement().

◆ _number_cross_slip_directions

const Real CrystalPlasticityStressUpdateBase::_number_cross_slip_directions
protectedinherited

Parameters to characterize the cross slip behavior of the crystal.

Definition at line 203 of file CrystalPlasticityStressUpdateBase.h.

Referenced by CrystalPlasticityStressUpdateBase::identifyCrossSlipFamily(), and CrystalPlasticityStressUpdateBase::sortCrossSlipFamilies().

◆ _number_cross_slip_planes

const Real CrystalPlasticityStressUpdateBase::_number_cross_slip_planes
protectedinherited

◆ _number_slip_systems

const unsigned int CrystalPlasticityStressUpdateBase::_number_slip_systems
protectedinherited

Maximum number of active slip systems for the crystalline material being modeled.

Definition at line 197 of file CrystalPlasticityStressUpdateBase.h.

Referenced by CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateConstitutiveSlipDerivative(), CrystalPlasticityTwinningKalidindiUpdate::calculateConstitutiveSlipDerivative(), calculateConstitutiveSlipDerivative(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateEquivalentSlipIncrement(), calculateEquivalentSlipIncrement(), CrystalPlasticityStressUpdateBase::calculateEquivalentSlipIncrement(), CrystalPlasticityStressUpdateBase::calculateFlowDirection(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateForestDislocationDensity(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateForestDislocationEvolutionIncrement(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateGrainSizeResistance(), CrystalPlasticityStressUpdateBase::calculateSchmidTensor(), CrystalPlasticityStressUpdateBase::calculateShearStress(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateSlipRate(), CrystalPlasticityTwinningKalidindiUpdate::calculateSlipRate(), calculateSlipRate(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateSlipResistance(), CrystalPlasticityTwinningKalidindiUpdate::calculateStateVariableEvolutionRateComponent(), calculateStateVariableEvolutionRateComponent(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateSubstructureDensityEvolutionIncrement(), CrystalPlasticityStressUpdateBase::calculateTotalPlasticDeformationGradientDerivative(), CrystalPlasticityTwinningKalidindiUpdate::calculateTwinResistance(), CrystalPlasticityTwinningKalidindiUpdate::calculateTwinVolumeFraction(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::CrystalPlasticityHCPDislocationSlipBeyerleinUpdate(), CrystalPlasticityStressUpdateBase::getSlipSystems(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::initQpStatefulProperties(), CrystalPlasticityTwinningKalidindiUpdate::initQpStatefulProperties(), initQpStatefulProperties(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::setMaterialVectorSize(), CrystalPlasticityTwinningKalidindiUpdate::setMaterialVectorSize(), CrystalPlasticityStressUpdateBase::setMaterialVectorSize(), CrystalPlasticityStressUpdateBase::sortCrossSlipFamilies(), CrystalPlasticityStressUpdateBase::transformHexagonalMillerBravaisSlipSystems(), and updateStateVariables().

◆ _previous_substep_slip_resistance

std::vector<Real> CrystalPlasticityKalidindiUpdate::_previous_substep_slip_resistance
protected

Stores the values of the slip system resistance from the previous substep In classes which use dislocation densities, analogous dislocation density substep vectors will be required.

Definition at line 116 of file CrystalPlasticityKalidindiUpdate.h.

Referenced by areConstitutiveStateVariablesConverged(), setInitialConstitutiveVariableValues(), setSubstepConstitutiveVariableValues(), updateStateVariables(), and updateSubstepConstitutiveVariableValues().

◆ _print_convergence_message

const bool CrystalPlasticityStressUpdateBase::_print_convergence_message
protectedinherited

◆ _r

const Real CrystalPlasticityKalidindiUpdate::_r
protected

Varibles used in the Kalidindi 1992 slip system resistance constiutive model.

Definition at line 93 of file CrystalPlasticityKalidindiUpdate.h.

Referenced by calculateStateVariableEvolutionRateComponent().

◆ _rel_state_var_tol

Real CrystalPlasticityStressUpdateBase::_rel_state_var_tol
protectedinherited

◆ _resistance_tol

Real CrystalPlasticityStressUpdateBase::_resistance_tol
protectedinherited

◆ _slip_direction

std::vector<RealVectorValue> CrystalPlasticityStressUpdateBase::_slip_direction
protectedinherited

◆ _slip_incr_tol

Real CrystalPlasticityStressUpdateBase::_slip_incr_tol
protectedinherited

◆ _slip_increment

MaterialProperty<std::vector<Real> >& CrystalPlasticityStressUpdateBase::_slip_increment
protectedinherited

◆ _slip_plane_normal

std::vector<RealVectorValue> CrystalPlasticityStressUpdateBase::_slip_plane_normal
protectedinherited

◆ _slip_resistance

MaterialProperty<std::vector<Real> >& CrystalPlasticityStressUpdateBase::_slip_resistance
protectedinherited

Slip system resistance.

Definition at line 217 of file CrystalPlasticityStressUpdateBase.h.

Referenced by CrystalPlasticityTwinningKalidindiUpdate::areConstitutiveStateVariablesConverged(), areConstitutiveStateVariablesConverged(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::areConstitutiveStateVariablesConverged(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::cacheStateVariablesBeforeUpdate(), CrystalPlasticityTwinningKalidindiUpdate::cacheStateVariablesBeforeUpdate(), cacheStateVariablesBeforeUpdate(), calculateConstitutiveSlipDerivative(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateSlipRate(), CrystalPlasticityTwinningKalidindiUpdate::calculateSlipRate(), calculateSlipRate(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateSlipResistance(), calculateStateVariableEvolutionRateComponent(), CrystalPlasticityTwinningKalidindiUpdate::calculateTwinResistance(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::initQpStatefulProperties(), CrystalPlasticityTwinningKalidindiUpdate::initQpStatefulProperties(), initQpStatefulProperties(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::setInitialConstitutiveVariableValues(), CrystalPlasticityTwinningKalidindiUpdate::setInitialConstitutiveVariableValues(), setInitialConstitutiveVariableValues(), CrystalPlasticityStressUpdateBase::setMaterialVectorSize(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::setSubstepConstitutiveVariableValues(), CrystalPlasticityTwinningKalidindiUpdate::setSubstepConstitutiveVariableValues(), setSubstepConstitutiveVariableValues(), updateStateVariables(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::updateSubstepConstitutiveVariableValues(), CrystalPlasticityTwinningKalidindiUpdate::updateSubstepConstitutiveVariableValues(), and updateSubstepConstitutiveVariableValues().

◆ _slip_resistance_before_update

std::vector<Real> CrystalPlasticityKalidindiUpdate::_slip_resistance_before_update
protected

Caches the value of the current slip system resistance immediately prior to the update of the slip system resistance, and is used to calculate the the slip system resistance increment for the current substep (or step if only one substep is taken) for the convergence check tolerance comparison.

In classes which use dislocation densities, analogous dislocation density caching vectors will also be required.

Definition at line 126 of file CrystalPlasticityKalidindiUpdate.h.

Referenced by areConstitutiveStateVariablesConverged(), and cacheStateVariablesBeforeUpdate().

◆ _slip_resistance_increment

std::vector<Real> CrystalPlasticityKalidindiUpdate::_slip_resistance_increment
protected

Increment of increased resistance for each slip system.

Definition at line 109 of file CrystalPlasticityKalidindiUpdate.h.

Referenced by calculateStateVariableEvolutionRateComponent(), and updateStateVariables().

◆ _slip_resistance_old

const MaterialProperty<std::vector<Real> >& CrystalPlasticityStressUpdateBase::_slip_resistance_old
protectedinherited

◆ _slip_sys_file_name

std::string CrystalPlasticityStressUpdateBase::_slip_sys_file_name
protectedinherited

File should contain slip plane normal and direction.

Definition at line 200 of file CrystalPlasticityStressUpdateBase.h.

Referenced by CrystalPlasticityStressUpdateBase::getSlipSystems().

◆ _substep_dt

Real CrystalPlasticityStressUpdateBase::_substep_dt
protectedinherited

◆ _tau

MaterialProperty<std::vector<Real> >& CrystalPlasticityStressUpdateBase::_tau
protectedinherited

◆ _tau_sat

const Real CrystalPlasticityKalidindiUpdate::_tau_sat
protected

◆ _twin_volume_fraction_total

const MaterialProperty<Real>* const CrystalPlasticityKalidindiUpdate::_twin_volume_fraction_total
protected

User-defined material property name for the total volume fraction of twins in a twinning propagation constitutive model, when this class is used in conjunction with the twinning propagation model.

Note that this value is the OLD material property and thus lags the current value by a single timestep.

Definition at line 141 of file CrystalPlasticityKalidindiUpdate.h.

◆ _unit_cell_dimension

const std::vector<Real> CrystalPlasticityStressUpdateBase::_unit_cell_dimension
protectedinherited

◆ _xm

const Real CrystalPlasticityKalidindiUpdate::_xm
protected

◆ _zero_tol

Real CrystalPlasticityStressUpdateBase::_zero_tol
protectedinherited

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