Finite Strain Crystal Plasticity

Crystal Plasticity base class: FCC system with power law flow rule implemented

Description

warning

The FiniteStrainCrystalPlasticity model is not actively developed. Use of the ComputeMultipleCrystalPlasticityStress system are recommend instead.

Constitutive models are used to calculate the plastic slip rate. In this crystal plasticity material the slip rate is modeled as a power law: $(1)var element = document.getElementById("moose-equation-4030cc98-2990-4163-97e4-03546c56699e");katex.render(" \\dot{\\gamma}^{\\alpha} = \\dot{\\gamma}_o \\left| \\frac{\\tau^{\\alpha}}{g^{\\alpha}} \\right|^{1/m} sign \\left( \\tau^{\\alpha} \\right)", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-2026fbad-a305-4532-ba62-2c602faa9c53");katex.render("\\dot{\\gamma}_o", element, {displayMode:false,throwOnError:false});$ is a reference slip rate, $var element = document.getElementById("moose-equation-33867484-d465-48cf-b832-708c6f1fac84");katex.render("\\tau^{\\alpha}", element, {displayMode:false,throwOnError:false});$ is the applied shear stress on each slip system $var element = document.getElementById("moose-equation-66a3424b-3bec-4738-b2b8-cb7231ecb6bf");katex.render("\\alpha", element, {displayMode:false,throwOnError:false});$ , $var element = document.getElementById("moose-equation-908cf8cf-6b1b-4aac-b5f9-5b2ab9f1a7c7");katex.render("g^{\\alpha}", element, {displayMode:false,throwOnError:false});$ is the slip system strength, or resistance to slip, and $var element = document.getElementById("moose-equation-b6ce73e8-69a6-4869-901b-6f53e7e39021");katex.render("m", element, {displayMode:false,throwOnError:false});$ is the strain rate sensitivity exponent. The strength of each slip system is solved with an iterative process as a function of the slip increment $(2)var element = document.getElementById("moose-equation-797d96b4-a386-4bcd-aaad-77b2db6d9ce0");katex.render(" g^{\\alpha} = g_o + \\Delta \\gamma^{\\alpha} q^{\\alpha \\beta} h _o \\left| 1 - \\frac{g^{\\alpha}}{g_{sat}} \\right|^a sign \\left( 1 - \\frac{g^{\\alpha}}{g_{sat}} \\right)", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-8e663be6-f949-4c01-8967-a440e22603cb");katex.render("q^{\\alpha \\beta}", element, {displayMode:false,throwOnError:false});$ is a hardening coefficient matrix that accounts for the different in self and latent hardening, Eq. (3), $var element = document.getElementById("moose-equation-408540d0-5b0f-4f69-8f22-20ee87ebaf35");katex.render("h_o", element, {displayMode:false,throwOnError:false});$ is an initial hardening term, $var element = document.getElementById("moose-equation-779a92a6-db61-41a0-8222-b7fba4f8b56b");katex.render("g_{sat}", element, {displayMode:false,throwOnError:false});$ is a constant saturated hardening value, and $var element = document.getElementById("moose-equation-75fb6915-9e9b-4534-9497-97e5623b911e");katex.render("a", element, {displayMode:false,throwOnError:false});$ is the hardening exponent (Kalidindi et al., 1992). The self and latent hardening of the crystal is defined for an FCC system as $(3)var element = document.getElementById("moose-equation-1c2263a8-8d79-41ec-8e2a-1e23a46c849e");katex.render(" q^{\\alpha \\beta} = \\begin{Bmatrix} 1.0 & q & q & q \\\\ q & 1.0 & q & q \\\\ q & q & 1.0 & q \\\\ q & q & q & 1.0 \\end{Bmatrix}", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-95566d9c-34aa-4492-b2df-55248916eb1d");katex.render("q", element, {displayMode:false,throwOnError:false});$ is a constant value of latent hardening among non-coplanar slip systems. In Eq. (3) the slip systems which share the same slip plane normal (e.g. $var element = document.getElementById("moose-equation-f6b4f81e-7f70-40a6-931f-1df59e086a42");katex.render("[\\bar{1}11]", element, {displayMode:false,throwOnError:false});$ ) are coplanar and grouped together with a latent hardening rate of unity (Kalidindi et al., 1992). Each matrix entry in Eq. (3) represents the interaction among two different coplanar slip system groups, that is a total of six slip systems (Kalidindi et al., 1992).

In this class, the crystal strain and stress response with response to the residual calculated by the second Piola-Kirchoff stress increment. In contrast, in FiniteStrainCPSlipRateRes the convergence of the crystal plasticity strain and stress response is determined with respect to the slip rate on each slip system of the crystal plasticity model.

Units Assumed in the Crystal Plasticity Materials

The simulation domain for crystal plasticity models is resolved on the order of individual crystal grains, and, as such, the mesh size is small. Although MOOSE itself is dimension agnostic, the crystal plasticity models are implemented in the mm-MPa-s unit system. This dimension system choice impacts the input files in the following manner:

Mesh dimensions should be constructed in mm
Elastic constant values (e.g. Young's modulus and shear modulus) are entered in MPa
Initial slip system strength values are entered in MPa
Simulation times are given in s
Strain rates and displacement loading rates are given in 1/s and mm/s, respectively

In physically based models, which maybe based on this class, initial densities of crystal defects (e.g. dislocations, point defects) should be given in 1/mm $var element = document.getElementById("moose-equation-86017caa-4882-46ef-9632-f0c68bac289a");katex.render("^2", element, {displayMode:false,throwOnError:false});$ or 1/mm $var element = document.getElementById("moose-equation-4b2a8b14-0590-4013-9bfd-d8b11f4ac622");katex.render("^3", element, {displayMode:false,throwOnError:false});$

Input Parameters

nssNumber of slip systems
C++ Type:int
Controllable:No
Description:Number of slip systems
slip_sys_file_nameName of the file containing the slip system
C++ Type:FileName
Controllable:No
Description:Name of the file containing the slip system

Required Parameters

abs_tol1e-06Constitutive stress residue absolute tolerance
Default:1e-06
C++ Type:double
Controllable:No
Description:Constitutive stress residue absolute tolerance
base_nameOptional parameter that allows the user to define multiple mechanics material systems on the same block, i.e. for multiple phases
C++ Type:std::string
Controllable:No
Description:Optional parameter that allows the user to define multiple mechanics material systems on the same block, i.e. for multiple phases
blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
boundaryThe list of boundaries (ids or names) from the mesh where this object applies
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The list of boundaries (ids or names) from the mesh where this object applies
computeTrueWhen false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.
Default:True
C++ Type:bool
Controllable:No
Description:When false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.
constant_onNONEWhen ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped
Default:NONE
C++ Type:MooseEnum
Options:NONE, ELEMENT, SUBDOMAIN
Controllable:No
Description:When ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped
declare_suffixAn optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Controllable:No
Description:An optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.
flowpropsParameters used in slip rate equations
C++ Type:std::vector<double>
Controllable:No
Description:Parameters used in slip rate equations
gen_random_stress_flagFalseFlag to generate random stress to perform time cutback on constitutive failure
Default:False
C++ Type:bool
Controllable:No
Description:Flag to generate random stress to perform time cutback on constitutive failure
gpropsInitial values of slip system resistances
C++ Type:std::vector<double>
Controllable:No
Description:Initial values of slip system resistances
gtol100Constitutive slip system resistance residual tolerance
Default:100
C++ Type:double
Controllable:No
Description:Constitutive slip system resistance residual tolerance
hpropsHardening properties
C++ Type:std::vector<double>
Controllable:No
Description:Hardening properties
input_random_scaling_varFalseFlag to input scaling variable: _Cijkl(0,0,0,0) when false
Default:False
C++ Type:bool
Controllable:No
Description:Flag to input scaling variable: _Cijkl(0,0,0,0) when false
intvar_read_typenoneRead from options for initial value of internal variables: Default from .i file
Default:none
C++ Type:MooseEnum
Options:slip_sys_file, slip_sys_res_file, none
Controllable:No
Description:Read from options for initial value of internal variables: Default from .i file
line_search_maxiter20Line search bisection method maximum number of iteration
Default:20
C++ Type:unsigned int
Controllable:No
Description:Line search bisection method maximum number of iteration
line_search_methodCUT_HALFThe method used in line search
Default:CUT_HALF
C++ Type:MooseEnum
Options:CUT_HALF, BISECTION
Controllable:No
Description:The method used in line search
line_search_tol0.5Line search bisection method tolerance
Default:0.5
C++ Type:double
Controllable:No
Description:Line search bisection method tolerance
maximum_substep_iteration1Maximum number of substep iteration
Default:1
C++ Type:unsigned int
Controllable:No
Description:Maximum number of substep iteration
maxiter100Maximum number of iterations for stress update
Default:100
C++ Type:unsigned int
Controllable:No
Description:Maximum number of iterations for stress update
maxitergss100Maximum number of iterations for slip system resistance update
Default:100
C++ Type:unsigned int
Controllable:No
Description:Maximum number of iterations for slip system resistance update
min_line_search_step_size0.01Minimum line search step size
Default:0.01
C++ Type:double
Controllable:No
Description:Minimum line search step size
num_slip_sys_flowrate_props2Number of flow rate properties for a slip system
Default:2
C++ Type:unsigned int
Controllable:No
Description:Number of flow rate properties for a slip system
num_slip_sys_props0Number of slip system specific properties provided in the file containing slip system normals and directions
Default:0
C++ Type:unsigned int
Controllable:No
Description:Number of slip system specific properties provided in the file containing slip system normals and directions
prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
random_scaling_var1e+09Random scaling variable: Large value can cause non-positive definiteness
Default:1e+09
C++ Type:double
Controllable:No
Description:Random scaling variable: Large value can cause non-positive definiteness
random_seed2000Random integer used to generate random stress when constitutive failure occurs
Default:2000
C++ Type:unsigned int
Controllable:No
Description:Random integer used to generate random stress when constitutive failure occurs
read_prop_user_objectThe ElementReadPropertyFile GeneralUserObject to read element specific property values from file
C++ Type:UserObjectName
Controllable:No
Description:The ElementReadPropertyFile GeneralUserObject to read element specific property values from file
rtol1e-06Constitutive stress residue relative tolerance
Default:1e-06
C++ Type:double
Controllable:No
Description:Constitutive stress residue relative tolerance
slip_incr_tol0.02Maximum allowable slip in an increment
Default:0.02
C++ Type:double
Controllable:No
Description:Maximum allowable slip in an increment
slip_sys_flow_prop_file_nameName of the file containing the values of slip rate equation parameters
C++ Type:FileName
Controllable:No
Description:Name of the file containing the values of slip rate equation parameters
slip_sys_hard_prop_file_nameName of the file containing the values of hardness evolution parameters
C++ Type:FileName
Controllable:No
Description:Name of the file containing the values of hardness evolution parameters
slip_sys_res_prop_file_nameName of the file containing the initial values of slip system resistances
C++ Type:FileName
Controllable:No
Description:Name of the file containing the initial values of slip system resistances
tan_mod_typenoneType of tangent moduli for preconditioner: default elastic
Default:none
C++ Type:MooseEnum
Options:exact, none
Controllable:No
Description:Type of tangent moduli for preconditioner: default elastic
use_line_searchFalseUse line search in constitutive update
Default:False
C++ Type:bool
Controllable:No
Description:Use line search in constitutive update

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:Yes
Description:Set the enabled status of the MooseObject.
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Controllable:No
Description:Determines whether this object is calculated using an implicit or explicit form
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Controllable:No
Description:The seed for the master random number generator

Advanced Parameters

output_propertiesList of material properties, from this material, to output (outputs must also be defined to an output type)
C++ Type:std::vector<std::string>
Controllable:No
Description:List of material properties, from this material, to output (outputs must also be defined to an output type)
outputsnone Vector of output names where you would like to restrict the output of variables(s) associated with this object
Default:none
C++ Type:std::vector<OutputName>
Controllable:No
Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object

Outputs Parameters

References

Surya R Kalidindi, Curt A Bronkhorst, and Lallit Anand. Crystallographic texture evolution in bulk deformation processing of fcc metals. Journal of the Mechanics and Physics of Solids, 40(3):537–569, 1992.

@article{kalidindi1992,
    author = "Kalidindi, Surya R and Bronkhorst, Curt A and Anand, Lallit",
    title = "Crystallographic texture evolution in bulk deformation processing of FCC metals",
    journal = "Journal of the Mechanics and Physics of Solids",
    volume = "40",
    number = "3",
    pages = "537--569",
    year = "1992",
    publisher = "Elsevier"
}