Crystal Plasticity State Var Rate Component Voce

Phenomenological Voce constitutive model state variable evolution rate component base class.

Description

This UserObject is meant to be used within the user object base crystal plasticity framework as a Voce hardening rule.

$(1)var element = document.getElementById("moose-equation-a153840c-39e2-40dd-9158-9869e5ced519");katex.render(" h^\\alpha = h_0^\\alpha\\left(1-\\frac{\\tau^{\\alpha}_{CRSS}-\\tau_0^\\alpha}{\\tau_{Sat}^\\alpha-\\tau_0^\\alpha} \\right)^{m^\\alpha}", element, {displayMode:true,throwOnError:false});$

with

$(2)var element = document.getElementById("moose-equation-ec4f217f-bebc-48f6-8c1e-33c5d4dcf1dc");katex.render(" \\dot\\tau^\\alpha_{CRSS} = \\sum_{\\beta=1}^N \\lvert \\dot\\gamma^\\beta \\rvert h^\\beta q^{\\alpha \\beta}", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-be4d6a9d-0e24-4915-951f-2f2ce082af3c");katex.render("\\tau^\\alpha_{CRSS}", element, {displayMode:false,throwOnError:false});$ is the current critical resolved shear stress, $var element = document.getElementById("moose-equation-9cca5032-27c3-4ede-a46d-e4687da0b631");katex.render("\\tau^{\\alpha}_{0}", element, {displayMode:false,throwOnError:false});$ is initial critical resolved shear stress, $var element = document.getElementById("moose-equation-27bdd19d-62d9-455d-ab2a-8cef885743cc");katex.render("\\tau_{Sat}^\\alpha", element, {displayMode:false,throwOnError:false});$ is the saturation resolved shear stress, $var element = document.getElementById("moose-equation-99c6aef3-17b0-48e4-b3d8-701abcc01736");katex.render("q^{\\alpha \\beta}", element, {displayMode:false,throwOnError:false});$ is the matrix containing self and latent hardening coefficients; $var element = document.getElementById("moose-equation-c2a5bca5-0025-4471-a1cf-3c38f32ca41d");katex.render("h_0^\\alpha", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-f09137cc-6dee-4bb3-b293-b462366b4dd9");katex.render("m^\\alpha", element, {displayMode:false,throwOnError:false});$ are hardening parameters.

Like for the others crystal plasticity routines values of all the parameters can be specified by groups of slip systems.

This UseroOject assumes that slip systems are provided in a well defined order and are grouped by planes (see modules/tensor_mechanics/test/tests/cp_user_object/input_slip_sys_bcc48.txt for BCC crystals and modules/tensor_mechanics/test/tests/cp_user_object/input_slip_sys.txt for FCC)

Seven variables need to be specified and one value is required for each group:

groups in which groups of slip systems are listed
h0_group_values $var element = document.getElementById("moose-equation-4bbb733c-485e-4938-bba8-9be7b881edeb");katex.render("h_0", element, {displayMode:false,throwOnError:false});$ value for each slip system group
tau0_group_values $var element = document.getElementById("moose-equation-26f915bd-6f1e-4fab-b5fd-84b1bc721829");katex.render("\\tau_0", element, {displayMode:false,throwOnError:false});$ value for each slip system group
tauSat_group_values = $var element = document.getElementById("moose-equation-264de294-cfaf-424c-9a85-8b31c7de5fa4");katex.render("\\tau_{Sat}", element, {displayMode:false,throwOnError:false});$ value for each slip system group
hardeningExponent_group_values $var element = document.getElementById("moose-equation-80a4642b-a2f3-4025-b862-0e8541e214a8");katex.render("m", element, {displayMode:false,throwOnError:false});$ value for each slip system group
selfHardening_group_values $var element = document.getElementById("moose-equation-b3258ff0-fd75-4724-8678-024701a924b0");katex.render("q^{\\alpha\\alpha}", element, {displayMode:false,throwOnError:false});$
coplanarHardening_group_values $var element = document.getElementById("moose-equation-35c61c38-653b-4352-9f0a-fc5f3f9bf7c7");katex.render("q^{\\alpha\\beta}", element, {displayMode:false,throwOnError:false});$ for co-planar slip systems in the same group
GroupGroup_Hardening_group_values $var element = document.getElementById("moose-equation-c11559a5-4241-4ca7-a348-511fd2213869");katex.render("q^{\\alpha\\beta}", element, {displayMode:false,throwOnError:false});$ This parameter requires NxN values with N being the number of groups. Values are listed as $var element = document.getElementById("moose-equation-cefe9994-26ac-4012-b685-d1e711b688e4");katex.render("value_ij", element, {displayMode:false,throwOnError:false});$ (e.g. ij=11,12,21,22) with $var element = document.getElementById("moose-equation-daa85a3e-6f9b-4c09-8242-eb70bc23ceb6");katex.render("i", element, {displayMode:false,throwOnError:false});$ being the actual group and $var element = document.getElementById("moose-equation-2d4c0ad8-bec4-4de1-97b4-f20d92736024");katex.render("j", element, {displayMode:false,throwOnError:false});$ the secondary group. For $var element = document.getElementById("moose-equation-8c9be16a-2cb1-49a5-a7eb-0a43708b8ea4");katex.render("i=j", element, {displayMode:false,throwOnError:false});$ the value represents the latent hardening coefficient between one slip system and all the non co-planar ones belonging the the same group. For $var element = document.getElementById("moose-equation-202f94ca-2787-4378-8121-9d3f39b55f13");katex.render("i\\neq j", element, {displayMode:false,throwOnError:false});$ , the value represents the latent hardening coefficient between all the slip systems belonging to group $var element = document.getElementById("moose-equation-2408243f-99d8-44bb-904b-05cb41cf88fb");katex.render("i", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-3443b896-8c9d-4369-a74e-4b6026f90596");katex.render("j", element, {displayMode:false,throwOnError:false});$ . Note that, usually, the matrix associated with GroupGroup_Hardening_group_values is symmetric.

Example Input File Syntax

[./state_var_evol_rate_comp_voce]
  type = CrystalPlasticityStateVarRateComponentVoce
  variable_size = 48
  crystal_lattice_type = 'BCC'
  groups = '0 12 24 48'
  h0_group_values = '1 2 3'
  tau0_group_values = '50 51 52'
  tauSat_group_values = '70 81 92'
  hardeningExponent_group_values = '1 2 3'
  selfHardening_group_values = '4 5 6'
  coplanarHardening_group_values = '7 8 9'
  GroupGroup_Hardening_group_values = '10 20 30
                                         40 50 60
                                         70 80 90'
  uo_slip_rate_name = slip_rate_gss
  uo_state_var_name = state_var_gss
[../]

In this example illustrate a BCC in which slip systems have been grouped by slip system class.

Input Parameters

variable_sizeThe userobject's variable size.
C++ Type:unsigned int
Controllable:No
Description:The userobject's variable size.

Required Parameters

GroupGroup_Hardening_group_valuesThe group-to-group latent hardening coefficient q_abThis is a NxN vector i.e. '1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0' means non-coplanar slip systems in gr_11,22,33= 1.0, 5.0 and 9.0 respectively.latent hardening between for gr_12,13 = 2.0 3.0 respectively
C++ Type:std::vector<double>
Controllable:No
Description:The group-to-group latent hardening coefficient q_abThis is a NxN vector i.e. '1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0' means non-coplanar slip systems in gr_11,22,33= 1.0, 5.0 and 9.0 respectively.latent hardening between for gr_12,13 = 2.0 3.0 respectively
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
coplanarHardening_group_valuesThe coplanar latent hardening coefficient q_abcorresponding to each group i.e. '1.0 2.0 3.0' means 0-11 = 1.0, 12-23 = 2.0 and 24-48 = 3.0
C++ Type:std::vector<double>
Controllable:No
Description:The coplanar latent hardening coefficient q_abcorresponding to each group i.e. '1.0 2.0 3.0' means 0-11 = 1.0, 12-23 = 2.0 and 24-48 = 3.0
crystal_lattice_typeFCCType of crystal lattyce structure output
Default:FCC
C++ Type:MooseEnum
Options:FCC, BCC
Controllable:No
Description:Type of crystal lattyce structure output
groupsTo group the initial values on different slip systems 'format: [start end)', i.e.'0 12 24 48' groups 0-11, 12-23 and 24-48
C++ Type:std::vector<unsigned int>
Controllable:No
Description:To group the initial values on different slip systems 'format: [start end)', i.e.'0 12 24 48' groups 0-11, 12-23 and 24-48
h0_group_valuesh0 hardening constant for each group i.e. '0.0 1.0 2.0' means 0-11 = 0.0, 12-23 = 1.0 and 24-48 = 2.0
C++ Type:std::vector<double>
Controllable:No
Description:h0 hardening constant for each group i.e. '0.0 1.0 2.0' means 0-11 = 0.0, 12-23 = 1.0 and 24-48 = 2.0
hardeningExponent_group_valuesThe hardening exponent mcorresponding to each group i.e. '1.0 2.0 3.0' means 0-11 = 1.0, 12-23 = 2.0 and 24-48 = 3.0
C++ Type:std::vector<double>
Controllable:No
Description:The hardening exponent mcorresponding to each group i.e. '1.0 2.0 3.0' means 0-11 = 1.0, 12-23 = 2.0 and 24-48 = 3.0
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.
selfHardening_group_valuesThe self hardening coefficient q_aacorresponding to each group i.e. '1.0 2.0 3.0' means 0-11 = 1.0, 12-23 = 2.0 and 24-48 = 3.0 usually these are all 1.
C++ Type:std::vector<double>
Controllable:No
Description:The self hardening coefficient q_aacorresponding to each group i.e. '1.0 2.0 3.0' means 0-11 = 1.0, 12-23 = 2.0 and 24-48 = 3.0 usually these are all 1.
tau0_group_valuesThe initial critical resolved shear stresscorresponding to each group i.e. '100.0 110.0 120.0' means 0-11 = 100.0, 12-23 = 110.0 and 24-48 = 120.0
C++ Type:std::vector<double>
Controllable:No
Description:The initial critical resolved shear stresscorresponding to each group i.e. '100.0 110.0 120.0' means 0-11 = 100.0, 12-23 = 110.0 and 24-48 = 120.0
tauSat_group_valuesThe saturation resolved shear stresscorresponding to each group i.e. '150.0 170.0 180.0' means 0-11 = 150.0, 12-23 = 170.0 and 24-48 = 180.0
C++ Type:std::vector<double>
Controllable:No
Description:The saturation resolved shear stresscorresponding to each group i.e. '150.0 170.0 180.0' means 0-11 = 150.0, 12-23 = 170.0 and 24-48 = 180.0
uo_slip_rate_nameName of slip rate property: Same as slip rate user object specified in input file.
C++ Type:std::string
Controllable:No
Description:Name of slip rate property: Same as slip rate user object specified in input file.
uo_state_var_nameName of state variable property: Same as state variable user object specified in input file.
C++ Type:std::string
Controllable:No
Description:Name of state variable property: Same as state variable user object specified in input file.

Optional Parameters

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Controllable:No
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
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.
execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
Default:0
C++ Type:int
Controllable:No
Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREAUX
force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup
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
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

(moose/modules/tensor_mechanics/test/tests/crystal_plasticity/user_object_based/user_object_Voce_BCC.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  elem_type = QUAD4
  displacements = 'disp_x disp_y'
  nx = 2
  ny = 2
[]

[GlobalParams]
  volumetric_locking_correction = true
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
[]

[Kernels]
  [./TensorMechanics]
    displacements = 'disp_x disp_y'
    use_displaced_mesh = true
  [../]
[]

[AuxVariables]
  [./stress_yy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./e_yy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./fp_yy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./gss]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[Functions]
  [./tdisp]
    type = ParsedFunction
    expression = 0.01*t
  [../]
[]

[UserObjects]
  [./prop_read]
    type = PropertyReadFile
    prop_file_name = 'euler_ang_file.txt'
    # Enter file data as prop#1, prop#2, .., prop#nprop
    nprop = 3
    read_type = element
  [../]
[]

[AuxKernels]
  [./stress_yy]
    type = RankTwoAux
    variable = stress_yy
    rank_two_tensor = stress
    index_j = 1
    index_i = 1
    execute_on = timestep_end
  [../]
  [./e_yy]
    type = RankTwoAux
    variable = e_yy
    rank_two_tensor = lage
    index_j = 1
    index_i = 1
    execute_on = timestep_end
  [../]
  [./fp_yy]
    type = RankTwoAux
    variable = fp_yy
    rank_two_tensor = fp
    index_j = 1
    index_i = 1
    execute_on = timestep_end
  [../]
  [./gss]
    type = MaterialStdVectorAux
    variable = gss
    property = state_var_gss
    index = 0
    execute_on = timestep_end
  [../]
[]

[BCs]
  [./fix_x]
    type = DirichletBC
    variable = disp_x
    boundary = 'left'
    value = 0
  [../]
  [./fix_y]
    type = DirichletBC
    variable = disp_y
    boundary = 'bottom'
    value = 0
  [../]
  [./tdisp]
    type = FunctionDirichletBC
    variable = disp_y
    boundary = top
    function = tdisp
  [../]
[]

[UserObjects]
  [./slip_rate_gss]
    type = CrystalPlasticitySlipRateGSS
    variable_size = 48
    slip_sys_file_name = input_slip_sys_bcc48.txt
    num_slip_sys_flowrate_props = 2
    flowprops = '1 12 0.001 0.1 13 24 0.001 0.1 25 48 0.001 0.1'
    uo_state_var_name = state_var_gss
  [../]
  [./slip_resistance_gss]
    type = CrystalPlasticitySlipResistanceGSS
    variable_size = 48
    uo_state_var_name = state_var_gss
  [../]
  [./state_var_gss]
    type = CrystalPlasticityStateVariable
    variable_size = 48
    groups = '0 12 24 48'
    group_values =  '50 51 52'
    uo_state_var_evol_rate_comp_name = state_var_evol_rate_comp_voce
    scale_factor = 1.0
  [../]
  [./state_var_evol_rate_comp_voce]
    type = CrystalPlasticityStateVarRateComponentVoce
    variable_size = 48
    crystal_lattice_type = 'BCC'
    groups = '0 12 24 48'
    h0_group_values = '1 2 3'
    tau0_group_values = '50 51 52'
    tauSat_group_values = '70 81 92'
    hardeningExponent_group_values = '1 2 3'
    selfHardening_group_values ='4 5 6'
    coplanarHardening_group_values='7 8 9'
    GroupGroup_Hardening_group_values = '10 20 30
                                         40 50 60
                                         70 80 90'
    uo_slip_rate_name = slip_rate_gss
    uo_state_var_name = state_var_gss
  [../]
[]

[Materials]
  [./crysp]
    type = FiniteStrainUObasedCP
    stol = 1e-2
    tan_mod_type = exact
    uo_slip_rates = 'slip_rate_gss'
    uo_slip_resistances = 'slip_resistance_gss'
    uo_state_vars = 'state_var_gss'
    uo_state_var_evol_rate_comps = 'state_var_evol_rate_comp_voce'
  [../]
  [./strain]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
  [../]
  [./elasticity_tensor]
    type = ComputeElasticityTensorCP
    C_ijkl = '1.684e5 1.214e5 1.214e5 1.684e5 1.214e5 1.684e5 0.754e5 0.754e5 0.754e5'
    fill_method = symmetric9
    read_prop_user_object = prop_read
  [../]
[]

[Postprocessors]
  [./stress_yy]
    type = ElementAverageValue
    variable = stress_yy
  [../]
  [./e_yy]
    type = ElementAverageValue
    variable = e_yy
  [../]
  [./fp_yy]
    type = ElementAverageValue
    variable = fp_yy
  [../]
  [./gss]
    type = ElementAverageValue
    variable = gss
  [../]
[]

[Preconditioning]
  [./smp]
    type = SMP
    full = true
  [../]
[]

[Executioner]
  type = Transient
  dt = 0.01
  solve_type = 'PJFNK'

  petsc_options_iname = -pc_hypre_type
  petsc_options_value = boomerang
  nl_abs_tol = 1e-10
  nl_rel_step_tol = 1e-10
  dtmax = 10.0
  nl_rel_tol = 1e-10
  end_time = 1
  dtmin = 0.01
  num_steps = 10
  nl_abs_step_tol = 1e-10
[]

[Outputs]
  exodus = true
[]