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

CrystalPlasticityHCPDislocationSlipBeyerleinUpdate computes the dislocation forest evolution for the prismatic, pyramidal, and basal slip systems. More...

#include <CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h>

Inheritance diagram for CrystalPlasticityHCPDislocationSlipBeyerleinUpdate:
[legend]

Public Types

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

Public Member Functions

 CrystalPlasticityHCPDislocationSlipBeyerleinUpdate (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 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 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 PK2 stress, resolved shear stress, plastic deformation gradient, slip system resistances, etc. More...
 
virtual void setMaterialVectorSize () override
 
virtual void setInitialConstitutiveVariableValues () override
 This virtual method is called to set the constitutive internal state variables current value and the previous substep value to the old property value for the start of the stress convergence while loop. More...
 
virtual void setSubstepConstitutiveVariableValues () override
 This virtual method is called to set the current constitutive internal state variable value to that of the previous substep at the beginning of the next substep increment. More...
 
virtual void updateSubstepConstitutiveVariableValues () override
 Stores the current value of the constitutive internal state variables into a separate material property in case substepping is required, once the constitutive variables have passed convergence tolerances. 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 HCP structures proposed by Beyerlein and Tome. 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...
 
void calculateForestDislocationEvolutionIncrement ()
 Computes the evolution increment of the forest (glide) dislocations on each slip system, from equation 3.6 in Beyerlein and Tome IJP 24 (2008) 5, 867-895, and using the increment approach shown in Ardeljan et al. More...
 
void calculateSubstructureDensityEvolutionIncrement ()
 Computes the evolution of the total substructure dislocations, using the nomenclature introduced in Capolungo et al. More...
 
bool calculateForestDislocationDensity ()
 Calculate the current value of the incremented forest dislocation density on each slip system. More...
 
bool calculateSubstructureDislocationDensity ()
 Calculates the current value of the incremented substructure dislocations. More...
 
virtual void calculateSlipResistance () override
 Calculates the sum of the contribution of the initial slip resistance and the forest and substructure dislocations to the slip resistance on each slip system, equation 3.16 from Beyerlein and Tome IJP 24 (2008) 5, 867-895. More...
 
virtual void calculateGrainSizeResistance (DenseVector< Real > &lattice_resistance)
 Calculates the influence of the grain size on the slip system resistance, in the case where twin boundaries are not considered, following equation 3.21 from Beyerlein and Tome IJP 24 (2008) 5, 867-895. More...
 
virtual bool areConstitutiveStateVariablesConverged () override
 Determines if the state variables, e.g. More...
 
bool isSubstructureDislocationDensityConverged ()
 Determine if the substructure (debris) dislocation density values 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

const VariableValue_temperature
 Coupled temperature variable. More...
 
const unsigned int _slip_system_modes
 The number of the different slip systems types to consider in the simulation, which is dependent on the specific material. More...
 
const std::vector< unsigned int_number_slip_systems_per_mode
 The number of slip systems per type, stored as a vector. More...
 
const std::vector< Real_shear_modulus
 Forest hardening coefficients, Eqns. 3.19. More...
 
const std::vector< Real_substructure_rate_coefficient
 Substructure (debris) dislocation generation rate coefficient. More...
 
const Real _substructure_hardening_coefficient
 Substructure (debris) Taylor law hardening coefficient. More...
 
const std::vector< Real_hallpetch_like_coefficient
 Microscale Hall-Petch like coefficient, used to capture the effect of the grain size on the slip system resistance in the case without twins (or before twins have formed). More...
 
const Real _grain_size
 Value of the grain size, either average or associated with a particular block. 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
 
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
 Slip forest dislocation densities. More...
 
const MaterialProperty< std::vector< Real > > & _forest_dislocation_density_old
 
MaterialProperty< std::vector< Real > > & _forest_dislocation_increment
 
MaterialProperty< std::vector< Real > > & _forest_dislocations_removed_increment
 
const Real _initial_forest_dislocation_density
 
MaterialProperty< Real > & _total_substructure_density
 Sessile substructure dislocation density across all slip systems. More...
 
const MaterialProperty< Real > & _total_substructure_density_old
 
MaterialProperty< Real > & _total_substructure_density_increment
 
const Real _initial_substructure_density
 
const std::vector< Real_lattice_friction
 
const Real _reference_strain_rate
 Power-law slip rate calculation coefficients, from Wang et al IJP 49(2013)36-52. More...
 
const Real _rate_sensitivity_exponent
 
const std::vector< Real_burgers_vector
 Coefficients for slip dislocation evolution, Eqns 3.6 and 3.8. More...
 
const std::vector< Real_slip_generation_coefficient
 
const std::vector< Real_slip_activation_energy
 
const std::vector< Real_proportionality_factor
 
const Real _forest_interaction_coefficient
 
const Real _boltzman_constant
 
const Real _macro_applied_strain_rate
 
const Real _macro_reference_strain_rate
 
std::vector< Real_previous_substep_slip_resistance
 Stores the slip system resistance, dislocation densities from the previous substep. More...
 
std::vector< Real_previous_substep_forest_dislocations
 
Real _previous_substep_total_substructure_density
 
std::vector< Real_slip_resistance_before_update
 Caching current slip resistance, dislocation density values before final update. More...
 
std::vector< Real_forest_dislocations_before_update
 
Real _total_substructure_density_before_update
 
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

CrystalPlasticityHCPDislocationSlipBeyerleinUpdate computes the dislocation forest evolution for the prismatic, pyramidal, and basal slip systems.

Definition at line 21 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.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

◆ CrystalPlasticityHCPDislocationSlipBeyerleinUpdate()

CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::CrystalPlasticityHCPDislocationSlipBeyerleinUpdate ( const InputParameters parameters)

Definition at line 130 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

132 
133  _temperature(coupledValue("temperature")),
135  declareProperty<std::vector<Real>>(_base_name + "forest_dislocation_density")),
137  getMaterialPropertyOld<std::vector<Real>>(_base_name + "forest_dislocation_density")),
139  declareProperty<std::vector<Real>>(_base_name + "forest_dislocation_increment")),
141  declareProperty<std::vector<Real>>(_base_name + "_forest_dislocations_removed_increment")),
142  _initial_forest_dislocation_density(getParam<Real>("initial_forest_dislocation_density")),
143  _total_substructure_density(declareProperty<Real>(_base_name + "total_substructure_density")),
145  getMaterialPropertyOld<Real>(_base_name + "total_substructure_density")),
147  declareProperty<Real>(_base_name + "total_substructure_increment")),
148  _initial_substructure_density(getParam<Real>("initial_substructure_density")),
149 
150  _slip_system_modes(getParam<unsigned int>("slip_system_modes")),
152  getParam<std::vector<unsigned int>>("number_slip_systems_per_mode")),
153  _lattice_friction(getParam<std::vector<Real>>("lattice_friction_per_mode")),
154 
155  _reference_strain_rate(getParam<Real>("gamma_o")),
156  _rate_sensitivity_exponent(getParam<Real>("strain_rate_sensitivity_exponent")),
157 
158  _burgers_vector(getParam<std::vector<Real>>("burgers_vector_per_mode")),
160  getParam<std::vector<Real>>("slip_generation_coefficient_per_mode")),
162  getParam<std::vector<Real>>("normalized_slip_activiation_energy_per_mode")),
164  getParam<std::vector<Real>>("slip_energy_proportionality_factor_per_mode")),
165  _forest_interaction_coefficient(getParam<Real>("forest_interaction_parameter")),
166  _boltzman_constant(getParam<Real>("Boltzman_constant")),
167  _macro_applied_strain_rate(getParam<Real>("applied_strain_rate")),
168  _macro_reference_strain_rate(getParam<Real>("reference_macroscopic_strain_rate")),
169 
170  _shear_modulus(getParam<std::vector<Real>>("effective_shear_modulus_per_mode")),
172  getParam<std::vector<Real>>("substructure_rate_coefficient_per_mode")),
173  _substructure_hardening_coefficient(getParam<Real>("substructure_hardening_coefficient")),
174  _hallpetch_like_coefficient(getParam<std::vector<Real>>("Hall_Petch_like_constant_per_mode")),
175  _grain_size(getParam<Real>("grain_size")),
176 
177  // Twinning contributions, if used
178  _include_twinning_in_Lp(parameters.isParamValid("total_twin_volume_fraction")),
180  ? &getMaterialPropertyOld<Real>("total_twin_volume_fraction")
181  : nullptr)
182 {
183  // resize local caching vectors used for substepping
188 
189  // check that the number of slip systems is equal to the sum of the types of slip system
191  paramError("number_slip_systems_per_mode",
192  "The size the number of slip systems per mode is not equal to the number of slip "
193  "system types.");
194 
195  // Check that the number of slip mode dependent parameters is given matches the number of slip
196  // modes
197  if (_burgers_vector.size() != _slip_system_modes)
198  paramError("burgers_vector_per_mode",
199  "Please ensure that the size of burgers_vector_per_mode equals the value supplied "
200  "for slip_system_modes");
201 
203  paramError("slip_generation_coefficient_per_mode",
204  "Please ensure that the size of slip_generation_coefficient_per_mode equals the "
205  "value supplied for slip_system_modes");
206 
208  paramError("normalized_slip_activiation_energy_per_mode",
209  "Please ensure that the size of normalized_slip_activiation_energy_per_mode equals "
210  "the value supplied for slip_system_modes");
211 
213  paramError("slip_energy_proportionality_factor_per_mode",
214  "Please ensure that the size of slip_energy_proportionality_factor_per_mode equals "
215  "the value supplied for slip_system_modes");
216 
217  if (_shear_modulus.size() != _slip_system_modes)
218  paramError("effective_shear_modulus_per_mode",
219  "Please ensure that the size of effective_shear_modulus_per_mode equals the "
220  "value supplied for slip_system_modes");
221 
223  paramError("substructure_rate_coefficient_per_mode",
224  "Please ensure that the size of substructure_rate_coefficient_per_mode equals the "
225  "value supplied for slip_system_modes");
226 
228  paramError("Hall_Petch_like_constant_per_mode",
229  "Please ensure that the size of Hall_Petch_like_constant_per_mode equals the value "
230  "supplied for slip_system_modes");
231 
233  paramError("lattice_friction_per_mode",
234  "Please ensure that the size of lattice_friction_per_mode equals the value supplied "
235  "for slip_system_modes");
236 
237  unsigned int sum = 0;
238  for (const auto i : make_range(_slip_system_modes))
240  if (sum != _number_slip_systems)
241  paramError("slip_system_modes",
242  "The number of slip systems and the sum of the slip systems in each of the slip "
243  "system modes are not equal");
244 }
MaterialProperty< Real > & _total_substructure_density
Sessile substructure dislocation density across all slip systems.
const std::vector< Real > _shear_modulus
Forest hardening coefficients, Eqns. 3.19.
MaterialProperty< T > & declareProperty(const std::string &name)
CrystalPlasticityStressUpdateBase(const InputParameters &parameters)
const unsigned int _slip_system_modes
The number of the different slip systems types to consider in the simulation, which is dependent on t...
const MaterialProperty< T > & getMaterialPropertyOld(const std::string &name, MaterialData &material_data)
const std::string _base_name
Base name prepended to all material property names to allow for multi-material systems.
const std::vector< Real > _burgers_vector
Coefficients for slip dislocation evolution, Eqns 3.6 and 3.8.
const Real _reference_strain_rate
Power-law slip rate calculation coefficients, from Wang et al IJP 49(2013)36-52.
const MaterialProperty< Real > *const _twin_volume_fraction_total
User-defined material property name for the total volume fraction of twins in a twinning propagation ...
const std::vector< Real > _hallpetch_like_coefficient
Microscale Hall-Petch like coefficient, used to capture the effect of the grain size on the slip syst...
virtual const VariableValue & coupledValue(const std::string &var_name, unsigned int comp=0) const
std::vector< Real > _previous_substep_slip_resistance
Stores the slip system resistance, dislocation densities from the previous substep.
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
const std::vector< unsigned int > _number_slip_systems_per_mode
The number of slip systems per type, stored as a vector.
const T & getParam(const std::string &name) const
void paramError(const std::string &param, Args... args) const
std::vector< Real > _slip_resistance_before_update
Caching current slip resistance, dislocation density values before final update.
const Real _grain_size
Value of the grain size, either average or associated with a particular block.
IntRange< T > make_range(T beg, T end)
const InputParameters & parameters() const
const Real _substructure_hardening_coefficient
Substructure (debris) Taylor law hardening coefficient.
const bool _include_twinning_in_Lp
Flag to include the total twin volume fraction in the plastic velocity gradient calculation, per Kalidindi IJP (2001).
const std::vector< Real > _substructure_rate_coefficient
Substructure (debris) dislocation generation rate coefficient.
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.
bool isParamValid(const std::string &name) const

Member Function Documentation

◆ areConstitutiveStateVariablesConverged()

bool CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::areConstitutiveStateVariablesConverged ( )
overrideprotectedvirtual

Determines if the state variables, e.g.

defect densities, have converged by comparing the change in the values over the iteration period.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 397 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

398 {
408  return true;
409  return false;
410 }
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
std::vector< Real > _previous_substep_slip_resistance
Stores the slip system resistance, dislocation densities from the previous substep.
Real _rel_state_var_tol
Internal variable update equation tolerance.
std::vector< Real > _slip_resistance_before_update
Caching current slip resistance, dislocation density values before final update.
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.
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.
bool isSubstructureDislocationDensityConverged()
Determine if the substructure (debris) dislocation density values have converged. ...

◆ cacheStateVariablesBeforeUpdate()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::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 438 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

439 {
443 }
MaterialProperty< Real > & _total_substructure_density
Sessile substructure dislocation density across all slip systems.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
std::vector< Real > _slip_resistance_before_update
Caching current slip resistance, dislocation density values before final update.
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.

◆ calculateConstitutiveSlipDerivative()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::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 383 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

385 {
386  for (const auto i : make_range(_number_slip_systems))
387  {
389  dslip_dtau[i] = 0.0;
390  else
391  dslip_dtau[i] = _slip_increment[_qp][i] /
392  (_rate_sensitivity_exponent * std::abs(_tau[_qp][i])) * _substep_dt;
393  }
394 }
bool absoluteFuzzyEqual(const T &var1, const T2 &var2, const T3 &tol=libMesh::TOLERANCE *libMesh::TOLERANCE)
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< Real > > & _tau
Resolved shear stress on each slip system.
IntRange< T > make_range(T beg, T end)

◆ calculateEquivalentSlipIncrement()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateEquivalentSlipIncrement ( RankTwoTensor equivalent_slip_increment)
overrideprotectedvirtual

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 369 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

371 {
373  {
374  for (const auto i : make_range(_number_slip_systems))
375  equivalent_slip_increment += (1.0 - (*_twin_volume_fraction_total)[_qp]) *
377  }
378  else // if no twinning volume fraction material property supplied, use base class
380 }
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.

◆ calculateForestDislocationDensity()

bool CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateForestDislocationDensity ( )
protected

Calculate the current value of the incremented forest dislocation density on each slip system.

Definition at line 601 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

Referenced by updateStateVariables().

602 {
603  for (const auto i : make_range(_number_slip_systems))
604  {
608  else
611 
612  if (_forest_dislocation_density[_qp][i] < 0.0)
613  return false;
614  }
615  return true;
616 }
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.
IntRange< T > make_range(T beg, T end)
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.

◆ calculateForestDislocationEvolutionIncrement()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateForestDislocationEvolutionIncrement ( )
protected

Computes the evolution increment of the forest (glide) dislocations on each slip system, from equation 3.6 in Beyerlein and Tome IJP 24 (2008) 5, 867-895, and using the increment approach shown in Ardeljan et al.

Journal of Mechanics and Physics of Solids 66 (2014) 16-31, equation 9.

Definition at line 453 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

Referenced by calculateStateVariableEvolutionRateComponent().

454 {
455  DenseVector<Real> k1_term(_number_slip_systems);
456  DenseVector<Real> k2_term(_number_slip_systems);
457 
458  const Real temperature_strain_term =
461 
462  // solve first for the coefficients, which depend on the given slip mode
463  unsigned int slip_mode = 0;
464  unsigned int counter_adjustment = 0;
465  for (const auto i : make_range(_number_slip_systems))
466  {
467  Real interaction_term = 0.0;
468  Real volume_term = 0.0;
469  if ((i - counter_adjustment) < _number_slip_systems_per_mode[slip_mode])
470  {
471  k1_term(i) = _slip_generation_coefficient[slip_mode];
472  interaction_term = _forest_interaction_coefficient * _burgers_vector[slip_mode] /
473  _slip_activation_energy[slip_mode];
474  volume_term =
475  _proportionality_factor[slip_mode] * Utility::pow<3>(_burgers_vector[slip_mode]);
476  }
477  else
478  {
479  counter_adjustment += _number_slip_systems_per_mode[slip_mode];
480  ++slip_mode;
481 
482  k1_term(i) = _slip_generation_coefficient[slip_mode];
483  interaction_term = _forest_interaction_coefficient * _burgers_vector[slip_mode] /
484  _slip_activation_energy[slip_mode];
485  volume_term =
486  _proportionality_factor[slip_mode] * Utility::pow<3>(_burgers_vector[slip_mode]);
487  }
488  k2_term(i) = interaction_term * k1_term(i) * (1.0 - temperature_strain_term / volume_term);
489  }
490 
491  for (const auto i : make_range(_number_slip_systems))
492  {
493  const Real abs_slip_increment = std::abs(_slip_increment[_qp][i]);
494  Real generated_dislocations = 0.0;
495 
496  if (_forest_dislocation_density[_qp][i] > 0.0)
497  generated_dislocations = k1_term(i) * std::sqrt(_forest_dislocation_density[_qp][i]) *
498  abs_slip_increment * _substep_dt;
499 
501  k2_term(i) * _forest_dislocation_density[_qp][i] * abs_slip_increment * _substep_dt;
502 
504  generated_dislocations - _forest_dislocations_removed_increment[_qp][i];
505  }
506 }
MaterialProperty< std::vector< Real > > & _slip_increment
Current slip increment material property.
const std::vector< Real > _burgers_vector
Coefficients for slip dislocation evolution, Eqns 3.6 and 3.8.
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. ...
const std::vector< unsigned int > _number_slip_systems_per_mode
The number of slip systems per type, stored as a vector.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
IntRange< T > make_range(T beg, T end)
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.

◆ calculateGrainSizeResistance()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateGrainSizeResistance ( DenseVector< Real > &  lattice_resistance)
protectedvirtual

Calculates the influence of the grain size on the slip system resistance, in the case where twin boundaries are not considered, following equation 3.21 from Beyerlein and Tome IJP 24 (2008) 5, 867-895.

These values are added to the initial lattice friction vector. In the case where twins are considered, this method should be overwritten to account for the influence of the twin dislocations.

Definition at line 299 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

Referenced by calculateSlipResistance(), and initQpStatefulProperties().

301 {
302  unsigned int slip_mode = 0;
303  unsigned int counter_adjustment = 0;
304  for (const auto i : make_range(_number_slip_systems))
305  {
306  Real hallpetch_burgers_term = 0.0;
307  if ((i - counter_adjustment) < _number_slip_systems_per_mode[slip_mode])
308  hallpetch_burgers_term = _hallpetch_like_coefficient[slip_mode] * _shear_modulus[slip_mode] *
309  std::sqrt(_burgers_vector[slip_mode]);
310  else
311  {
312  counter_adjustment += _number_slip_systems_per_mode[slip_mode];
313  ++slip_mode;
314  hallpetch_burgers_term = _hallpetch_like_coefficient[slip_mode] * _shear_modulus[slip_mode] *
315  std::sqrt(_burgers_vector[slip_mode]);
316  }
317  lattice_resistance(i) += hallpetch_burgers_term / std::sqrt(_grain_size);
318  }
319 }
const std::vector< Real > _shear_modulus
Forest hardening coefficients, Eqns. 3.19.
const std::vector< Real > _burgers_vector
Coefficients for slip dislocation evolution, Eqns 3.6 and 3.8.
const std::vector< Real > _hallpetch_like_coefficient
Microscale Hall-Petch like coefficient, used to capture the effect of the grain size on the slip syst...
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
const std::vector< unsigned int > _number_slip_systems_per_mode
The number of slip systems per type, stored as a vector.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const Real _grain_size
Value of the grain size, either average or associated with a particular block.
IntRange< T > make_range(T beg, T end)

◆ 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 CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::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 343 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

344 {
345  for (const auto i : make_range(_number_slip_systems))
346  {
347  Real driving_force = std::abs(_tau[_qp][i] / _slip_resistance[_qp][i]);
348  if (driving_force < _zero_tol)
349  _slip_increment[_qp][i] = 0.0;
350  else
351  {
352  _slip_increment[_qp][i] =
354  if (_tau[_qp][i] < 0.0)
355  _slip_increment[_qp][i] *= -1.0;
356  }
357  if (std::abs(_slip_increment[_qp][i]) * _substep_dt > _slip_incr_tol)
358  {
360  mooseWarning("Maximum allowable slip increment exceeded ",
361  std::abs(_slip_increment[_qp][i]) * _substep_dt);
362  return false;
363  }
364  }
365  return true;
366 }
MaterialProperty< std::vector< Real > > & _slip_increment
Current slip increment material property.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
const Real _reference_strain_rate
Power-law slip rate calculation coefficients, from Wang et al IJP 49(2013)36-52.
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. ...
Real _zero_tol
Residual tolerance when variable value is zero. Default 1e-12.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
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()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateSlipResistance ( )
overrideprotectedvirtual

Calculates the sum of the contribution of the initial slip resistance and the forest and substructure dislocations to the slip resistance on each slip system, equation 3.16 from Beyerlein and Tome IJP 24 (2008) 5, 867-895.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 538 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

539 {
540  DenseVector<Real> forest_hardening(_number_slip_systems, 0.0);
541  DenseVector<Real> substructure_hardening(_number_slip_systems, 0.0);
542  DenseVector<Real> lattice_resistance(_number_slip_systems, 0.0);
543 
544  unsigned int slip_mode = 0;
545  unsigned int counter_adjustment = 0;
546  for (const auto i : make_range(_number_slip_systems))
547  {
548  Real burgers = 0.0;
549  Real shear_modulus = 0.0;
550  if ((i - counter_adjustment) < _number_slip_systems_per_mode[slip_mode])
551  {
552  burgers = _burgers_vector[slip_mode];
553  shear_modulus = _shear_modulus[slip_mode];
554  lattice_resistance(i) = _lattice_friction[slip_mode];
555  }
556  else
557  {
558  counter_adjustment += _number_slip_systems_per_mode[slip_mode];
559  ++slip_mode;
560  burgers = _burgers_vector[slip_mode];
561  shear_modulus = _shear_modulus[slip_mode];
562  lattice_resistance(i) = _lattice_friction[slip_mode];
563  }
564 
565  // forest dislocation hardening
566  if (_forest_dislocation_density[_qp][i] > 0.0)
567  forest_hardening(i) = _forest_interaction_coefficient * burgers * shear_modulus *
568  std::sqrt(_forest_dislocation_density[_qp][i]);
569  else
570  forest_hardening(i) = 0.0;
571 
572  // substructure dislocation hardening
573  if (_total_substructure_density[_qp] > 0.0)
574  {
575  const Real spacing_term = burgers * std::sqrt(_total_substructure_density[_qp]);
576  substructure_hardening(i) = _substructure_hardening_coefficient * shear_modulus *
577  spacing_term * std::log10(1.0 / spacing_term);
578  }
579  else
580  substructure_hardening(i) = 0.0;
581  }
582 
583  calculateGrainSizeResistance(lattice_resistance);
584 
585  // have the constant initial value, while it's not a function of temperature, sum
586  for (const auto i : make_range(_number_slip_systems))
587  _slip_resistance[_qp][i] =
588  lattice_resistance(i) + forest_hardening(i) + substructure_hardening(i);
589 }
MaterialProperty< Real > & _total_substructure_density
Sessile substructure dislocation density across all slip systems.
virtual void calculateGrainSizeResistance(DenseVector< Real > &lattice_resistance)
Calculates the influence of the grain size on the slip system resistance, in the case where twin boun...
const std::vector< Real > _shear_modulus
Forest hardening coefficients, Eqns. 3.19.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
const std::vector< Real > _burgers_vector
Coefficients for slip dislocation evolution, Eqns 3.6 and 3.8.
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
const std::vector< unsigned int > _number_slip_systems_per_mode
The number of slip systems per type, stored as a vector.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
IntRange< T > make_range(T beg, T end)
const Real _substructure_hardening_coefficient
Substructure (debris) Taylor law hardening coefficient.
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.

◆ calculateStateVariableEvolutionRateComponent()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateStateVariableEvolutionRateComponent ( )
overrideprotectedvirtual

Following the constitutive model for HCP structures proposed by Beyerlein and Tome.

A dislocation-based constitutive law for pure Zr including temperature effects. International Journal of Plasticity 24(2008)5, 867-895. This class includes only the contributions for forest dislocation evolution and the associated hardening contributions due to the forest and substructure dislocations. Contributions for twin dislocations, as given in the paper above, are included in a separate class.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 446 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

447 {
450 }
void calculateSubstructureDensityEvolutionIncrement()
Computes the evolution of the total substructure dislocations, using the nomenclature introduced in C...
void calculateForestDislocationEvolutionIncrement()
Computes the evolution increment of the forest (glide) dislocations on each slip system, from equation 3.6 in Beyerlein and Tome IJP 24 (2008) 5, 867-895, and using the increment approach shown in Ardeljan et al.

◆ calculateSubstructureDensityEvolutionIncrement()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateSubstructureDensityEvolutionIncrement ( )
protected

Computes the evolution of the total substructure dislocations, using the nomenclature introduced in Capolungo et al.

Materials Science and Engineering A 513-514 (2009) 42-51, and the equation form given in Zecevic et al. International Journal of Plasticity 84 (2016) 138-159.

Definition at line 509 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

Referenced by calculateStateVariableEvolutionRateComponent().

510 {
511  // calculate the generation coefficient, which depends on the slip mode
512  DenseVector<Real> generation_term(_number_slip_systems, 0.0);
513 
514  unsigned int slip_mode = 0;
515  unsigned int counter_adjustment = 0;
516  for (const auto i : make_range(_number_slip_systems))
517  {
518  if ((i - counter_adjustment) < _number_slip_systems_per_mode[slip_mode])
519  generation_term(i) = _substructure_rate_coefficient[slip_mode] * _burgers_vector[slip_mode];
520  else
521  {
522  counter_adjustment += _number_slip_systems_per_mode[slip_mode];
523  ++slip_mode;
524  generation_term(i) = _substructure_rate_coefficient[slip_mode] * _burgers_vector[slip_mode];
525  }
526  }
527 
528  // perform the summing calculation over all slip systems
530  const Real sqrt_substructures = std::sqrt(_total_substructure_density[_qp]);
531 
532  for (const auto i : make_range(_number_slip_systems))
534  generation_term(i) * sqrt_substructures * _forest_dislocations_removed_increment[_qp][i];
535 }
MaterialProperty< Real > & _total_substructure_density
Sessile substructure dislocation density across all slip systems.
const std::vector< Real > _burgers_vector
Coefficients for slip dislocation evolution, Eqns 3.6 and 3.8.
unsigned int _qp
const unsigned int _number_slip_systems
Maximum number of active slip systems for the crystalline material being modeled. ...
const std::vector< unsigned int > _number_slip_systems_per_mode
The number of slip systems per type, stored as a vector.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
IntRange< T > make_range(T beg, T end)
const std::vector< Real > _substructure_rate_coefficient
Substructure (debris) dislocation generation rate coefficient.

◆ calculateSubstructureDislocationDensity()

bool CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::calculateSubstructureDislocationDensity ( )
protected

Calculates the current value of the incremented substructure dislocations.

Definition at line 619 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

Referenced by updateStateVariables().

◆ 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 CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::initQpStatefulProperties ( )
overrideprotectedvirtual

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

This class is often overwritten by inherting classes.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 247 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

248 {
250 
251  // Resize constitutive-model specific material properties
253 
254  // Set constitutive-model specific initial values from parameters
255  const Real forest_density_per_system = _initial_forest_dislocation_density / _number_slip_systems;
256  for (const auto i : make_range(_number_slip_systems))
257  {
258  _forest_dislocation_density[_qp][i] = forest_density_per_system;
260  _slip_increment[_qp][i] = 0.0;
261  }
262 
263  // Set initial resistance from lattice friction, which is type dependent
264  DenseVector<Real> lattice_resistance(_number_slip_systems, 0.0);
265  unsigned int slip_mode = 0;
266  unsigned int counter_adjustment = 0;
267  for (const auto i : make_range(_number_slip_systems))
268  {
269  if ((i - counter_adjustment) < _number_slip_systems_per_mode[slip_mode])
270  lattice_resistance(i) = _lattice_friction[slip_mode];
271  else
272  {
273  counter_adjustment += _number_slip_systems_per_mode[slip_mode];
274  ++slip_mode;
275  lattice_resistance(i) = _lattice_friction[slip_mode];
276  }
277  }
278 
279  calculateGrainSizeResistance(lattice_resistance);
280 
281  for (const auto i : make_range(_number_slip_systems))
282  _slip_resistance[_qp][i] = lattice_resistance(i);
283 
286 }
MaterialProperty< Real > & _total_substructure_density
Sessile substructure dislocation density across all slip systems.
virtual void calculateGrainSizeResistance(DenseVector< Real > &lattice_resistance)
Calculates the influence of the grain size on the slip system resistance, in the case where twin boun...
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.
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. ...
const std::vector< unsigned int > _number_slip_systems_per_mode
The number of slip systems per type, stored as a vector.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
IntRange< T > make_range(T beg, T end)
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.

◆ 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(), CrystalPlasticityKalidindiUpdate::areConstitutiveStateVariablesConverged(), and 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)

◆ isSubstructureDislocationDensityConverged()

bool CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::isSubstructureDislocationDensityConverged ( )
protected

Determine if the substructure (debris) dislocation density values have converged.

Definition at line 413 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

Referenced by areConstitutiveStateVariablesConverged().

414 {
415  bool converged_flag = true;
416 
417  Real substructure_diff =
419 
420  if (_previous_substep_total_substructure_density < _zero_tol && substructure_diff > _zero_tol)
421  converged_flag = false;
424  converged_flag = false;
425 
426  return converged_flag;
427 }
MaterialProperty< Real > & _total_substructure_density
Sessile substructure dislocation density across all slip systems.
unsigned int _qp
Real _zero_tol
Residual tolerance when variable value is zero. Default 1e-12.
Real _rel_state_var_tol
Internal variable update equation tolerance.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ 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 CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::setInitialConstitutiveVariableValues ( )
overrideprotectedvirtual

This virtual method is called to set the constitutive internal state variables current value and the previous substep value to the old property value for the start of the stress convergence while loop.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 322 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

323 {
326 
329 
332 }
MaterialProperty< Real > & _total_substructure_density
Sessile substructure dislocation density across all slip systems.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
std::vector< Real > _previous_substep_slip_resistance
Stores the slip system resistance, dislocation densities from the previous substep.
const MaterialProperty< std::vector< Real > > & _slip_resistance_old
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.

◆ setMaterialVectorSize()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::setMaterialVectorSize ( )
overrideprotectedvirtual

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 289 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

290 {
292 
293  // Resize non-stateful material properties
296 }
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. ...

◆ 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 CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::setSubstepConstitutiveVariableValues ( )
overrideprotectedvirtual

This virtual method is called to set the current constitutive internal state variable value to that of the previous substep at the beginning of the next substep increment.

In cases where only one substep is taken (or when the first) substep is taken, this method sets the current constitutive internal state variable value to the old value.

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 335 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

336 {
340 }
MaterialProperty< Real > & _total_substructure_density
Sessile substructure dislocation density across all slip systems.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
std::vector< Real > _previous_substep_slip_resistance
Stores the slip system resistance, dislocation densities from the previous substep.
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.

◆ 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 CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::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 592 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

593 {
595  return true;
596  else
597  return false;
598 }
bool calculateSubstructureDislocationDensity()
Calculates the current value of the incremented substructure dislocations.
bool calculateForestDislocationDensity()
Calculate the current value of the incremented forest dislocation density on each slip system...

◆ updateSubstepConstitutiveVariableValues()

void CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::updateSubstepConstitutiveVariableValues ( )
overrideprotectedvirtual

Stores the current value of the constitutive internal state variables into a separate material property in case substepping is required, once the constitutive variables have passed convergence tolerances.

This separate material property is used as the previous substep value in the associated setSubstepConstitutiveVariableValues method in the next substep (if taken).

Reimplemented from CrystalPlasticityStressUpdateBase.

Definition at line 430 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

431 {
435 }
MaterialProperty< Real > & _total_substructure_density
Sessile substructure dislocation density across all slip systems.
MaterialProperty< std::vector< Real > > & _slip_resistance
Slip system resistance.
unsigned int _qp
std::vector< Real > _previous_substep_slip_resistance
Stores the slip system resistance, dislocation densities from the previous substep.
MaterialProperty< std::vector< Real > > & _forest_dislocation_density
Slip forest dislocation densities.

◆ validParams()

InputParameters CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::validParams ( )
static

Definition at line 16 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.C.

17 {
19  params.addClassDescription("Two-term dislocation slip model for hexagonal close packed crystals "
20  "from Beyerline and Tome");
21 
22  params.set<MooseEnum>("crystal_lattice_type") = "HCP";
23  params.suppressParameter<MooseEnum>("crystal_lattice_type");
24 
25  params.addCoupledVar("temperature", "The name of the temperature variable");
27  "initial_forest_dislocation_density",
28  "initial_forest_dislocation_density>0",
29  "The initial density of the forest dislocations, in 1/mm^2, assumed "
30  "to be split evenly among all slip systems");
32  "initial_substructure_density",
33  "initial_substructure_density>0",
34  "The initial total density of the sessile dislocations, in 1/mm^2");
35 
36  params.addParam<unsigned int>(
37  "slip_system_modes",
38  1,
39  "Number of different types of slip systems in this HCP crystal, e.g. for a material with "
40  "basal<a>, prismatic<a>, and pyramidal<a> active slip systems, this number would be 3");
41  params.addParam<std::vector<unsigned int>>(
42  "number_slip_systems_per_mode",
43  std::vector<unsigned int>(),
44  "The number of slip systems per each slip system type. The sum of the entries of the vector "
45  "given here must equal the value given for the total number of slip systems.");
46  params.addParam<std::vector<Real>>(
47  "lattice_friction_per_mode",
48  std::vector<Real>(),
49  "Value of the lattice friction for each type of the slip system, units of MPa. The order "
50  "must be consistent with the number of slip systems per type vector.");
51 
52  params.addParam<std::vector<Real>>(
53  "effective_shear_modulus_per_mode",
54  std::vector<Real>(),
55  "Effective isotropic shear modulus value, mu, in MPa. The order "
56  "must be consistent with the number of slip systems per type vector.");
57 
58  params.addParam<std::vector<Real>>(
59  "burgers_vector_per_mode",
60  std::vector<Real>(),
61  "Value of the Burgers vector, b, for each type of the slip system, units of mm. The order "
62  "must "
63  "be consistent with the number of slip systems per type vector.");
64  params.addParam<std::vector<Real>>(
65  "slip_generation_coefficient_per_mode",
66  std::vector<Real>(),
67  "Slip dislocation generation coefficient value for each type of the slip system, k_1, units "
68  "of 1/mm. The order "
69  "must be consistent with the number of slip systems per type vector.");
70  params.addParam<std::vector<Real>>(
71  "normalized_slip_activiation_energy_per_mode",
72  std::vector<Real>(),
73  "Value of the slip dislocation attraction activation energy for each type of the slip "
74  "system, g, dimensionless. The order must be consistent with the number of slip systems per "
75  "type vector.");
76  params.addParam<std::vector<Real>>(
77  "slip_energy_proportionality_factor_per_mode",
78  std::vector<Real>(),
79  "Value of the the dislocation slip attraction energy proportionality factor for each type of "
80  "the slip system, D, units of MPa. The order must be consistent with the number of slip "
81  "systems "
82  "per type vector.");
83  params.addParam<std::vector<Real>>(
84  "substructure_rate_coefficient_per_mode",
85  std::vector<Real>(),
86  "Material-independent rate constant that accounts for locking of slip dislocations in "
87  "sessile substructure dislocation segments, q, dimensionless. This value is often determined "
88  "through dislocation dynamics calculations. The order must be consistent with the number of "
89  "slip systems per type vector.");
90 
91  params.addParam<Real>("gamma_o", 1.0e-3, "Reference strain rate on each slip system, in 1/s");
92  params.addParam<Real>("strain_rate_sensitivity_exponent",
93  0.05,
94  "The strain rate sensitivity exponent for the power law relationship of "
95  "resolved shear stress");
96  params.addParam<Real>("forest_interaction_parameter",
97  0.9,
98  "Forest dislocation interaction parameter, Chi, dimensionless.");
99  params.addParam<Real>("Boltzman_constant", 1.38065e-20, "Boltzman constant, in MPa-mm^3/K");
100  params.addRangeCheckedParam<Real>("applied_strain_rate",
101  1.0e-4,
102  "applied_strain_rate<=1.0e-3 & applied_strain_rate>=1.0e-5",
103  "Value of the applied macroscopic strain rate and should "
104  "correspond to the simulation loading conditions, in 1/s.");
105  params.addParam<Real>("reference_macroscopic_strain_rate",
106  1.0e7,
107  "Value of the reference macroscopic strain rate for the thermal "
108  "dislocation attraction, in 1/s.");
109  params.addParam<Real>("substructure_hardening_coefficient",
110  0.086,
111  "Value of the coefficient for the expanded Taylor hardening substructure "
112  "hardening relation, set to recover the Taylor hardening law for low "
113  "substructure densities, k_{sub}, dimensionless.");
114  params.addParam<std::vector<Real>>(
115  "Hall_Petch_like_constant_per_mode",
116  std::vector<Real>(),
117  "The microstructure Hall-Petch like coefficient value used to capture the influence of grain "
118  "size on slip system resistance in the absence of twin dislocations, dimensionless");
119  params.addRequiredRangeCheckedParam<Real>(
120  "grain_size", "grain_size>0", "Value of the crystal grain size, in mm");
121 
122  params.addParam<MaterialPropertyName>(
123  "total_twin_volume_fraction",
124  "Total twin volume fraction, if twinning is considered in the simulation");
125 
126  return params;
127 }
void addRequiredRangeCheckedParam(const std::string &name, const std::string &parsed_function, const std::string &doc_string)
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
T & set(const std::string &name, bool quiet_mode=false)
void suppressParameter(const std::string &name)
void addCoupledVar(const std::string &name, const std::string &doc_string)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
void addClassDescription(const std::string &doc_string)

Member Data Documentation

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

◆ _boltzman_constant

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_boltzman_constant
protected

◆ _burgers_vector

const std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_burgers_vector
protected

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

◆ _forest_dislocation_density

MaterialProperty<std::vector<Real> >& CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_forest_dislocation_density
protected

◆ _forest_dislocation_density_old

const MaterialProperty<std::vector<Real> >& CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_forest_dislocation_density_old
protected

◆ _forest_dislocation_increment

MaterialProperty<std::vector<Real> >& CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_forest_dislocation_increment
protected

◆ _forest_dislocations_before_update

std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_forest_dislocations_before_update
protected

◆ _forest_dislocations_removed_increment

MaterialProperty<std::vector<Real> >& CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_forest_dislocations_removed_increment
protected

◆ _forest_interaction_coefficient

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_forest_interaction_coefficient
protected

◆ _grain_size

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_grain_size
protected

Value of the grain size, either average or associated with a particular block.

Definition at line 182 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by calculateGrainSizeResistance().

◆ _hallpetch_like_coefficient

const std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_hallpetch_like_coefficient
protected

Microscale Hall-Petch like coefficient, used to capture the effect of the grain size on the slip system resistance in the case without twins (or before twins have formed).

Should not be confused with the engineering scale Hall-Petch coefficient.

Definition at line 179 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by calculateGrainSizeResistance(), and CrystalPlasticityHCPDislocationSlipBeyerleinUpdate().

◆ _include_twinning_in_Lp

const bool CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_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 200 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by calculateEquivalentSlipIncrement().

◆ _initial_forest_dislocation_density

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_initial_forest_dislocation_density
protected

◆ _initial_substructure_density

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_initial_substructure_density
protected

◆ _lattice_friction

const std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_lattice_friction
protected

Initial slip resistance vectors, with values arranaged as per each slip system mode

Definition at line 146 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by calculateSlipResistance(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate(), and initQpStatefulProperties().

◆ _macro_applied_strain_rate

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_macro_applied_strain_rate
protected

◆ _macro_reference_strain_rate

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_macro_reference_strain_rate
protected

◆ _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 calculateConstitutiveSlipDerivative(), CrystalPlasticityTwinningKalidindiUpdate::calculateConstitutiveSlipDerivative(), CrystalPlasticityKalidindiUpdate::calculateConstitutiveSlipDerivative(), calculateEquivalentSlipIncrement(), CrystalPlasticityKalidindiUpdate::calculateEquivalentSlipIncrement(), CrystalPlasticityStressUpdateBase::calculateEquivalentSlipIncrement(), CrystalPlasticityStressUpdateBase::calculateFlowDirection(), calculateForestDislocationDensity(), calculateForestDislocationEvolutionIncrement(), calculateGrainSizeResistance(), CrystalPlasticityStressUpdateBase::calculateSchmidTensor(), CrystalPlasticityStressUpdateBase::calculateShearStress(), calculateSlipRate(), CrystalPlasticityTwinningKalidindiUpdate::calculateSlipRate(), CrystalPlasticityKalidindiUpdate::calculateSlipRate(), calculateSlipResistance(), CrystalPlasticityTwinningKalidindiUpdate::calculateStateVariableEvolutionRateComponent(), CrystalPlasticityKalidindiUpdate::calculateStateVariableEvolutionRateComponent(), calculateSubstructureDensityEvolutionIncrement(), CrystalPlasticityStressUpdateBase::calculateTotalPlasticDeformationGradientDerivative(), CrystalPlasticityTwinningKalidindiUpdate::calculateTwinResistance(), CrystalPlasticityTwinningKalidindiUpdate::calculateTwinVolumeFraction(), CrystalPlasticityHCPDislocationSlipBeyerleinUpdate(), CrystalPlasticityStressUpdateBase::getSlipSystems(), initQpStatefulProperties(), CrystalPlasticityTwinningKalidindiUpdate::initQpStatefulProperties(), CrystalPlasticityKalidindiUpdate::initQpStatefulProperties(), setMaterialVectorSize(), CrystalPlasticityTwinningKalidindiUpdate::setMaterialVectorSize(), CrystalPlasticityStressUpdateBase::setMaterialVectorSize(), CrystalPlasticityStressUpdateBase::sortCrossSlipFamilies(), CrystalPlasticityStressUpdateBase::transformHexagonalMillerBravaisSlipSystems(), and CrystalPlasticityKalidindiUpdate::updateStateVariables().

◆ _number_slip_systems_per_mode

const std::vector<unsigned int> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_number_slip_systems_per_mode
protected

◆ _previous_substep_forest_dislocations

std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_previous_substep_forest_dislocations
protected

◆ _previous_substep_slip_resistance

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

◆ _previous_substep_total_substructure_density

Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_previous_substep_total_substructure_density
protected

◆ _print_convergence_message

const bool CrystalPlasticityStressUpdateBase::_print_convergence_message
protectedinherited

◆ _proportionality_factor

const std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_proportionality_factor
protected

◆ _rate_sensitivity_exponent

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_rate_sensitivity_exponent
protected

◆ _reference_strain_rate

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_reference_strain_rate
protected

Power-law slip rate calculation coefficients, from Wang et al IJP 49(2013)36-52.

Definition at line 149 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by calculateSlipRate().

◆ _rel_state_var_tol

Real CrystalPlasticityStressUpdateBase::_rel_state_var_tol
protectedinherited

◆ _resistance_tol

Real CrystalPlasticityStressUpdateBase::_resistance_tol
protectedinherited

◆ _shear_modulus

const std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_shear_modulus
protected

◆ _slip_activation_energy

const std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_slip_activation_energy
protected

◆ _slip_direction

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

◆ _slip_generation_coefficient

const std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_slip_generation_coefficient
protected

◆ _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(), CrystalPlasticityKalidindiUpdate::areConstitutiveStateVariablesConverged(), areConstitutiveStateVariablesConverged(), cacheStateVariablesBeforeUpdate(), CrystalPlasticityTwinningKalidindiUpdate::cacheStateVariablesBeforeUpdate(), CrystalPlasticityKalidindiUpdate::cacheStateVariablesBeforeUpdate(), CrystalPlasticityKalidindiUpdate::calculateConstitutiveSlipDerivative(), calculateSlipRate(), CrystalPlasticityTwinningKalidindiUpdate::calculateSlipRate(), CrystalPlasticityKalidindiUpdate::calculateSlipRate(), calculateSlipResistance(), CrystalPlasticityKalidindiUpdate::calculateStateVariableEvolutionRateComponent(), CrystalPlasticityTwinningKalidindiUpdate::calculateTwinResistance(), initQpStatefulProperties(), CrystalPlasticityTwinningKalidindiUpdate::initQpStatefulProperties(), CrystalPlasticityKalidindiUpdate::initQpStatefulProperties(), setInitialConstitutiveVariableValues(), CrystalPlasticityTwinningKalidindiUpdate::setInitialConstitutiveVariableValues(), CrystalPlasticityKalidindiUpdate::setInitialConstitutiveVariableValues(), CrystalPlasticityStressUpdateBase::setMaterialVectorSize(), setSubstepConstitutiveVariableValues(), CrystalPlasticityTwinningKalidindiUpdate::setSubstepConstitutiveVariableValues(), CrystalPlasticityKalidindiUpdate::setSubstepConstitutiveVariableValues(), CrystalPlasticityKalidindiUpdate::updateStateVariables(), updateSubstepConstitutiveVariableValues(), CrystalPlasticityTwinningKalidindiUpdate::updateSubstepConstitutiveVariableValues(), and CrystalPlasticityKalidindiUpdate::updateSubstepConstitutiveVariableValues().

◆ _slip_resistance_before_update

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

Caching current slip resistance, dislocation density values before final update.

Definition at line 191 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by areConstitutiveStateVariablesConverged(), cacheStateVariablesBeforeUpdate(), and CrystalPlasticityHCPDislocationSlipBeyerleinUpdate().

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

◆ _slip_system_modes

const unsigned int CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_slip_system_modes
protected

The number of the different slip systems types to consider in the simulation, which is dependent on the specific material.

Types of slip systems, e.g. prismatic, with two different burgers vectors, e.g. and <c + a>, must be considered as two separate slip system types.

Definition at line 139 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by CrystalPlasticityHCPDislocationSlipBeyerleinUpdate().

◆ _substep_dt

Real CrystalPlasticityStressUpdateBase::_substep_dt
protectedinherited

◆ _substructure_hardening_coefficient

const Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_substructure_hardening_coefficient
protected

Substructure (debris) Taylor law hardening coefficient.

Definition at line 171 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by calculateSlipResistance().

◆ _substructure_rate_coefficient

const std::vector<Real> CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_substructure_rate_coefficient
protected

Substructure (debris) dislocation generation rate coefficient.

Definition at line 168 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by calculateSubstructureDensityEvolutionIncrement(), and CrystalPlasticityHCPDislocationSlipBeyerleinUpdate().

◆ _tau

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

◆ _temperature

const VariableValue& CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_temperature
protected

Coupled temperature variable.

Definition at line 116 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

Referenced by calculateForestDislocationEvolutionIncrement().

◆ _total_substructure_density

MaterialProperty<Real>& CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_total_substructure_density
protected

◆ _total_substructure_density_before_update

Real CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_total_substructure_density_before_update
protected

◆ _total_substructure_density_increment

MaterialProperty<Real>& CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_total_substructure_density_increment
protected

◆ _total_substructure_density_old

const MaterialProperty<Real>& CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_total_substructure_density_old
protected

◆ _twin_volume_fraction_total

const MaterialProperty<Real>* const CrystalPlasticityHCPDislocationSlipBeyerleinUpdate::_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 209 of file CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.h.

◆ _unit_cell_dimension

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

◆ _zero_tol

Real CrystalPlasticityStressUpdateBase::_zero_tol
protectedinherited

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