Basal-prismatic slip systems Crystal Plasticity #27763

p-antonioni · 2024-05-31T19:35:17Z

p-antonioni
May 31, 2024

Hi,

I have a doubt on crystal plasticity considering multiple slip systems. I have 2 material phases of which the first is BCC and the second one is HCP. I should simulate sliding on both the basal and prismatic slip planes of the alpha (hcp) phase using CrystalPlasticityKalidindiUpdate which have different hardening paremeters. At this point, I've tried modelling as two different slip systems as two distinct phases using base_name but I get errors but maybe I'm not using it correctly. Do you have any suggestions? Should I modify the CrystalPlasticityKalidindiUpdate or is there another way to do it?

Thanks,
Paolo

GiudGiud · 2024-05-31T22:51:14Z

GiudGiud
May 31, 2024
Collaborator

@laagesen

0 replies

laagesen · 2024-06-04T17:03:36Z

laagesen
Jun 4, 2024

I would ask for help from @sapitts or @parikshitbajpai or @dewenyushu on this one for crystal plasticity

2 replies

parikshitbajpai Jun 4, 2024

I would ask for help from @sapitts or @parikshitbajpai or @dewenyushu on this one for crystal plasticity

I think you meant @vprithiv

laagesen Jun 4, 2024

yes sorry about that!

vprithiv · 2024-06-05T17:07:28Z

vprithiv
Jun 5, 2024

@p-antonioni

One way to approach this problem is to define two different blocks for those two phases and specify them in the material blocks, this potentially should work without changing anything in the code

Just an example below (make sure to use the actual objects)

[elasticity_phase1]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1
poissons_ratio = 0.1
block = 'phase1'
[]
[elasticity_phase2]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 2
poissons_ratio = 0.2
block = 'phase2'
[]

6 replies

vprithiv Jun 5, 2024

If it is a very similar issue, please post in this thread.

Emmanuel-324 Jun 5, 2024

@vprithiv Here is my input script.

[GlobalParams]
  displacements = 'disp_x disp_y'
[]

[Mesh]
  [file]
      type = FileMeshGenerator
      file = my_kks_noload_twophase_out.e-s180
      use_for_exodus_restart = true
  []
[]
[Variables]
  [eta0]
    initial_from_file_var = eta1
  []

  [eta1]
    initial_from_file_var = eta3
  []
[]
[AuxVariables]
  [pk2]
    order = CONSTANT
    family = MONOMIAL
  []
  [fp_xx]
    order = CONSTANT
    family = MONOMIAL
  []
  [rotout]
    order = CONSTANT
    family = MONOMIAL
  []
  [e_xx]
    order = CONSTANT
    family = MONOMIAL
  []
  [gss]
    order = CONSTANT
    family = MONOMIAL
  []
  [slip_increment]
   order = CONSTANT
   family = MONOMIAL
  []
[]

[Physics/SolidMechanics/QuasiStatic/all]
  strain = FINITE
  add_variables = true
  generate_output = stress_xx
[]

[AuxKernels]
  [fp_xx]
    type = RankTwoAux
    variable = fp_xx
    rank_two_tensor = plastic_deformation_gradient
    index_j = 2
    index_i = 2
    execute_on = timestep_end
  []
  [pk2]
   type = RankTwoAux
   variable = pk2
   rank_two_tensor = second_piola_kirchhoff_stress
   index_j = 2
   index_i = 2
   execute_on = timestep_end
  []
  [e_xx]
    type = RankTwoAux
    variable = e_xx
    rank_two_tensor = total_lagrangian_strain
    index_j = 2
    index_i = 2
    execute_on = timestep_end
  []
  [eta0_gss]
    type = MaterialStdVectorAux
    variable = gss
    property = slip_resistance
    index = 0
    execute_on = timestep_end
    block = 0
  []
  [eta1_gss]
    type = MaterialStdVectorAux
    variable = gss
    property = slip_resistance
    index = 0
    execute_on = timestep_end
  []
  [slip_inc]
   type = MaterialStdVectorAux
   variable = slip_increment
   property = slip_increment
   index = 0
   execute_on = timestep_end
  []
[]

[BCs]
  [symmy]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [symmx]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  
  [tdisp]
    type = FunctionDirichletBC
    variable = disp_y
    boundary = right
    function = '0.1*t'
  []
[]

[Materials]
  [elasticity_tensor_eta0]
    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
    base_name = eta0
  []
  [stress_eta0]
    type = ComputeMultipleCrystalPlasticityStress
    crystal_plasticity_models = 'trial_xtalpl_eta0'
    tan_mod_type = exact
    maximum_substep_iteration = 1
    base_name = eta0
  []
  [trial_xtalpl_eta0]
    type = CrystalPlasticityKalidindiUpdate
    number_slip_systems = 12
    slip_sys_file_name = input_slip_sys.txt
    base_name = eta0
  []

  [elasticity_tensor_eta1]
    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
    base_name = eta1
  []
  [stress_eta1]
    type = ComputeMultipleCrystalPlasticityStress
    crystal_plasticity_models = 'trial_xtalpl_eta1'
    tan_mod_type = exact
    maximum_substep_iteration = 1
    base_name = eta1
  []
  [trial_xtalpl_eta1]
    type = CrystalPlasticityKalidindiUpdate
    number_slip_systems = 12
    slip_sys_file_name = input_slip_sys.txt
    base_name = eta1
  []

   # Switching functions for each phase
  [h0]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta0
    all_etas = 'eta0 eta1'
    h_name = h0
  []

  # Generate the global stress from the phase stresses
  [global_stress]
    type = MultiPhaseStressMaterial
    phase_base = 'eta0 eta1'
    h = 'h0     h1'
  []
[]

[Kernels]
  [eta0_dt]
      type = TimeDerivative
      variable = eta0
  []
  [eta1_dt]
      type = TimeDerivative
      variable = eta1
  []
[]

[Postprocessors]
  [stress_xx]
    type = ElementAverageValue
    variable = stress_xx
    block = 'ANY_BLOCK_ID 0'
  []
  [pk2]
   type = ElementAverageValue
   variable = pk2
  []
  [fp_xx]
    type = ElementAverageValue
    variable = fp_xx
  []
  [e_xx]
    type = ElementAverageValue
    variable = e_xx
  []
  [gss]
    type = ElementAverageValue
    variable = gss
  []
  [slip_increment]
   type = ElementAverageValue
   variable = slip_increment
  []
[]

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

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

  petsc_options_iname = '-pc_type -pc_asm_overlap -sub_pc_type -ksp_type -ksp_gmres_restart'
  petsc_options_value = ' asm      2              lu            gmres     200'
  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
[]

Below is the output error

The following error occurred in the Problem 'MOOSE Problem' of type FEProblem.check_inp(moose_srv)
Material property 'slip_resistance', requested by 'eta0_gss' is not defined on block 0check_inp(moose_srv)
Material property 'slip_resistance', requested by 'eta1_gss' is not defined on block 0check_inp(moose_srv)
Material property 'h1', requested by 'global_stress' is not defined on block 0check_inp(moose_srv)
Material property 'slip_increment', requested by 'slip_inc' is not defined on block 0check_inp(moose_srv)
Material property 'eta0_mechanical_strain', requested by 'stress_eta0' is not defined on block 0check_inp(moose_srv)
Material property 'eta0_strain_increment', requested by 'stress_eta0' is not defined on block 0check_inp(moose_srv)
Material property 'eta0_rotation_increment', requested by 'stress_eta0' is not defined on block 0check_inp(moose_srv)
Material property 'eta0_deformation_gradient', requested by 'stress_eta0' is not defined on block 0check_inp(moose_srv)
Material property 'eta1_mechanical_strain', requested by 'stress_eta1' is not defined on block 0check_inp(moose_srv)
Material property 'eta1_strain_increment', requested by 'stress_eta1' is not defined on block 0check_inp(moose_srv)
Material property 'eta1_rotation_increment', requested by 'stress_eta1' is not defined on block 0check_inp(moose_srv)
Material property 'eta1_deformation_gradient', requested by 'stress_eta1' is not defined on block 0check_inp(moose_srv)

Emmanuel-324 Jun 5, 2024

To give some context I was using the input script below at first which was working but it seems it is a deprecated code now so I have been trying almost all the crystal plasticity kernels/input script in the stress_update_method repository but I get similar material property errors.

@vprithiv this is just a reminder of this thread. Thanks

[Mesh]
  [file]
    type = FileMeshGenerator
    file = kks_multiphase_largescale_large_out.e-s183
    use_for_exodus_restart = true
  []
[]

[Variables]
  [disp_x]
  []
  [disp_y]
  []

  # order parameter 0
  [eta0]
    initial_from_file_var = eta1
  []
  # order parameter 1
  [eta1]
    initial_from_file_var = eta2
  []
  [eta2]
    initial_from_file_var = eta3
  []
[]

[AuxVariables]
  [vonmises]
    order = CONSTANT
    family = MONOMIAL
  []
  [stress_xx]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
[]

[AuxKernels]
  [vonmises]
    type = RankTwoScalarAux
    rank_two_tensor = stress
    variable = vonmises
    scalar_type = VonMisesStress
    execute_on = timestep_end
    #   block = 0
  []
  [stress_xx]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_j = 0
    index_i = 0
    execute_on = timestep_end
    block = 0
  []
[]

[BCs]
  [symmy]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [symmx]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  [tdisp]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = right
    function = '0.1*t'
  []
[]

[Materials]
  [elasticity_tensor_phase0]
    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
    base_name = phase0
    block = 0
  []
  [crysp_phase0]
    type = FiniteStrainCPSlipRateRes
    slip_sys_file_name = input_slip_sys.txt
    nss = 12
    num_slip_sys_flowrate_props = 2 #Number of properties in a slip system
    flowprops = '1 4 1e-6 0.05 5 8 1e-6 0.05 9 12 1e-6 0.05'
    hprops = '1.0 215 78.6 200.0 1.75'
    gprops = '1 4 78.6 5 8 78.6 9 12 78.6'
    slip_incr_tol = 1
    maximum_substep_iteration = 12
    rtol = 1e-9
    abs_tol = 1e-10
    base_name = phase0
  []
  [strain_phase0]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase0
    #    eigenstrain_names = eigenstrain2
  []

  [elasticity_tensor_phase1]
    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
    euler_angle_1 = 0.0
    euler_angle_2 = 45.0
    euler_angle_3 = 0.0
    base_name = phase1
  []
  [crysp_phase1]
    type = FiniteStrainCPSlipRateRes
    slip_sys_file_name = input_slip_sys.txt
    nss = 12
    num_slip_sys_flowrate_props = 2 #Number of properties in a slip system
    flowprops = '1 4 1e-6 0.05 5 8 1e-6 0.05 9 12 1e-6 0.05'
    hprops = '1.0 250 118.6 230.0 1.75'
    gprops = '1 4 118.6 5 8 118.6 9 12 118.6'
    slip_incr_tol = 1
    maximum_substep_iteration = 12
    rtol = 1e-9
    abs_tol = 1e-10
    base_name = phase1
  []
  [strain_phase1]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase1
  []

  [elasticity_tensor_phase2]
    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
    euler_angle_1 = 0.0
    euler_angle_2 = 45.0
    euler_angle_3 = 30.0
    base_name = phase2
  []
  [crysp_phase2]
    type = FiniteStrainCPSlipRateRes
    slip_sys_file_name = input_slip_sys.txt
    nss = 12
    num_slip_sys_flowrate_props = 2 #Number of properties in a slip system
    flowprops = '1 4 1e-6 0.05 5 8 1e-6 0.05 9 12 1e-6 0.05'
    hprops = '1.0 250 118.6 230.0 1.75'
    gprops = '1 4 118.6 5 8 118.6 9 12 118.6'
    slip_incr_tol = 1
    maximum_substep_iteration = 12
    rtol = 1e-9
    abs_tol = 1e-10
    base_name = phase2
  []
  [strain_phase2]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase2
  []

  # Switching functions for each phase
  [h0]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta0
    all_etas = 'eta0 eta1 eta2'
    h_name = h0
  []
  [h1]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta1
    all_etas = 'eta0 eta1 eta2'
    h_name = h1
  []
  [h2]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta2
    all_etas = 'eta0 eta1 eta2'
    h_name = h2
  []

  # Generate the global stress from the phase stresses
  [global_stress]
    type = MultiPhaseStressMaterial
    phase_base = 'phase0 phase1 phase2'
    h_name  =   'h0     h1     h2'
  []
[]

[Kernels]
  [eta0_dt]
    type = TimeDerivative
    variable = eta0
  []
  [eta1_dt]
    type = TimeDerivative
    variable = eta1
  []
  [eta2_dt]
    type = TimeDerivative
    variable = eta2
  []

  [TensorMechanics]
    displacements = 'disp_x disp_y'
  []
[]

[Postprocessors]
  [vonmises]
    type = ElementAverageValue
    variable = vonmises
    block = 'ANY_BLOCK_ID 0'
  []
  [stress_xx]
    type = ElementAverageValue
    variable = stress_xx
    block = 'ANY_BLOCK_ID 0'
  []
[]

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

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -snes_type'
  petsc_options_value = 'lu            superlu_dist          vinewtonrsls'

  #  petsc_options_iname = '-pc_type -pc_asm_overlap -sub_pc_type -ksp_type -ksp_gmres_restart'
  #  petsc_options_value = ' asm      2              lu            gmres     200'
  l_max_its = 20
  nl_max_its = 10
  nl_rel_tol = 1.0e-8
  nl_abs_tol = 1.0e-9

  end_time = 300

  [TimeStepper]
    type = IterationAdaptiveDT
    dt = 5e-3
    cutback_factor = 0.75
    growth_factor = 1.2
    optimal_iterations = 20
  []
[]

[Outputs]
  exodus = true
  csv = true
[]

p-antonioni Jun 10, 2024
Author

Hi @vprithiv and thanks,

Sorry for the late response.
I have done an attempt to introduce another slip mechanism in the alpha phase in this manner:
I defined the alpha phase blocks and have multiple 'CrystalPlasticityKalidindiUpdate' for each slip system family, which I then passed to 'ComputeMultipleCrystalPlasticityStress'. It appears to work as it compiles, but I received the following error:

ComputeMultipleCrystalPlasticityStress: Constitutive failure
Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0

Printing state variable error messages returns:

Maximum allowable slip increment exceeded

An attempt to simulate to reduce the timestep was performed but it fails to converge.
Do you think is it simply due to tolerance settings or material properties (some properties are not properly set) in the simulation or it could be an issue of how i define the problem?
I put h = 0 only because i wanted to exclude hardening, i don't know if it is a smart way for doing it.
Last question, if I understood correctly the stresses of the two slip systems are separated. If so, do I have to calculate the global stress in a similar way as @Emmanuel-324 did?
The following lines I am referring to:

  [global_stress]
    type = MultiPhaseStressMaterial
    phase_base = 'phase0 phase1 phase2'
    h_name  =   'h0     h1     h2'
  []

This is the alpha phase definition which is used as my input file:

  [stress_alpha]
    type = ComputeMultipleCrystalPlasticityStress
    crystal_plasticity_models = 'alpha_xtalpl_basal alpha_xtalpl_prism'
    tan_mod_type = exact
    maximum_substep_iteration = 10
    block = '2 4 6 8 10 12 14 16 18 20 22 24 26'
  []

  [alpha_xtalpl_basal]
    type = CrystalPlasticityKalidindiUpdate
    number_slip_systems = 3
    slip_sys_file_name = input_slip_sys_hcp_basal.txt
    crystal_lattice_type = HCP
    ao = 0.001
    gss_a = 2.5 # Hardening coefficient
    gss_initial = 150 # Initial Lattice friction strength of the material
    h = 0 #hardening constants
    number_cross_slip_directions = 0 
    number_cross_slip_planes = 0 
    r = 1 
    t_sat = 250
    unit_cell_dimension = '2.934e-7 2.934e-7 4.657e-7'
    xm = 0.1
    resistance_tol = 1.0e-2
    block = '2 4 6 8 10 12 14 16 18 20 22 24 26'
    base_name = basal
  []
  [alpha_xtalpl_prism]
    type = CrystalPlasticityKalidindiUpdate
    number_slip_systems = 3
    slip_sys_file_name = input_slip_sys_hcp_prism.txt
    crystal_lattice_type = HCP
    ao = 0.001
    gss_a = 2.5 
    gss_initial = 150
    h = 0 #hardening constants
    number_cross_slip_directions = 0
    number_cross_slip_planes = 0
    r = 1 
    t_sat = 250 
    unit_cell_dimension = '2.934e-7 2.934e-7 4.657e-7''
    xm = 0.1 
    resistance_tol = 1.0e-2
    #block = '2 4 6 8 10 12 14 16 18 20 22 24 26'
    base_name = prism
  []

Emmanuel-324 Jun 13, 2024

@vprithiv, my model works now. I'll reach out if I need further input. Thanks

sapitts · 2024-06-24T18:20:22Z

sapitts
Jun 24, 2024
Collaborator

@p-antonioni in case you have not already resolved this issue, there is a specific HCP class available in MOOSE, see https://mooseframework.inl.gov/source/materials/crystal_plasticity/CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.html I would recommend using this HCP specific class over modifications to the Kalidindi class if you are still having issues.

7 replies

Emmanuel-324 Jun 24, 2024

Thanks @sapitts . I want to find papers that have used the CrystalPlasticityKalidindiUpdate class. Do you have any suggestions?

Emmanuel-324 Jun 24, 2024

Lastly, my script below works but does not output the right strain. I mean I get only zeros for phase0_e_xx and phase1_e_xx

[Mesh]
  [file]
    type = FileMeshGenerator
    file = my_kks_noload_twophase_out.e-s180
    use_for_exodus_restart = true
  []
[]

[Variables]
  [disp_x]
  []
  [disp_y]
  []

  # order parameter 0
  [eta0]
    initial_from_file_var = eta1
  []
  # order parameter 1
  [eta1]
    initial_from_file_var = eta3
  []
[]

[AuxVariables]
  [phase0_pk2]
    order = CONSTANT
    family = MONOMIAL
  []
  [phase0_fp_xx]
    order = CONSTANT
    family = MONOMIAL
  []
  [phase0_e_xx]
    order = CONSTANT
    family = MONOMIAL
  []
  [phase1_pk2]
    order = CONSTANT
    family = MONOMIAL
  []
  [phase1_fp_xx]
    order = CONSTANT
    family = MONOMIAL
  []
  [phase1_e_xx]
    order = CONSTANT
    family = MONOMIAL
  []

  [vonmises]
    order = CONSTANT
    family = MONOMIAL
  []
  [stress_xx]
    order = CONSTANT
    family = MONOMIAL
  []

  [phase0_gss]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [phase0_slip_increment]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [phase1_gss]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
  [phase1_slip_increment]
    order = CONSTANT
    family = MONOMIAL
    block = 0
  []
[]

[AuxKernels]
  [phase0_pk2]
    type = RankTwoAux
    variable = phase0_pk2
    rank_two_tensor = phase0_second_piola_kirchhoff_stress
    index_j = 2
    index_i = 2
    execute_on = timestep_end
  []
  [phase0_fp_xx]
    type = RankTwoAux
    variable = phase0_fp_xx
    rank_two_tensor = phase0_plastic_deformation_gradient
    index_j = 2
    index_i = 2
    execute_on = timestep_end
  []
  [phase0_e_xx]
    type = RankTwoAux
    variable = phase0_e_xx
    rank_two_tensor = phase0_total_lagrangian_strain
    index_j = 2
    index_i = 2
    execute_on = timestep_end
  []
  [phase1_pk2]
    type = RankTwoAux
    variable = phase1_pk2
    rank_two_tensor = phase1_second_piola_kirchhoff_stress
    index_j = 2
    index_i = 2
    execute_on = timestep_end
  []
  [phase1_fp_xx]
    type = RankTwoAux
    variable = phase1_fp_xx
    rank_two_tensor = phase1_plastic_deformation_gradient
    index_j = 2
    index_i = 2
    execute_on = timestep_end
  []
  [phase1_e_xx]
    type = RankTwoAux
    variable = phase1_e_xx
    rank_two_tensor = phase1_total_lagrangian_strain
    index_j = 2
    index_i = 2
    execute_on = timestep_end
  []

  [vonmises]
    type = RankTwoScalarAux
    rank_two_tensor = stress
    variable = vonmises
    scalar_type = VonMisesStress
    execute_on = timestep_end
    #   block = 0
  []

  [stress_xx]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_j = 0
    index_i = 0
    execute_on = timestep_end
  []

  [gss_phase0]
    type = MaterialStdVectorAux
    variable = phase0_gss
    property = phase0_slip_resistance
    index = 0
    execute_on = timestep_end
  []
  [slip_inc_phase0]
    type = MaterialStdVectorAux
    variable = phase0_slip_increment
    property = phase0_slip_increment
    index = 0
    execute_on = timestep_end
  []
  [gss_phase1]
    type = MaterialStdVectorAux
    variable = phase1_gss
    property = phase1_slip_resistance
    index = 0
    execute_on = timestep_end
  []
  [slip_inc_phase1]
    type = MaterialStdVectorAux
    variable = phase1_slip_increment
    property = phase1_slip_increment
    index = 0
    execute_on = timestep_end
  []
[]

[BCs]
  [symmy]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [symmx]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  [tdisp]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = right
    function = '0.1*t'
  []
[]

[Materials]
  [elasticity_tensor_phase0]
    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
    base_name = phase0
  []
  [stress_phase0]
    type = ComputeMultipleCrystalPlasticityStress
    crystal_plasticity_models = 'trial_xtalpl_phase0'
    rtol = 1e-06
    tan_mod_type = exact
    base_name = phase0
  []
  [trial_xtalpl_phase0]
    type = CrystalPlasticityKalidindiUpdate
    number_slip_systems = 12
    slip_sys_file_name = input_slip_sys.txt
    crystal_lattice_type = FCC
    r = 1.0              
    h = 541.5            
    t_sat = 109.8        
    gss_a = 2.5          
    ao = 0.001           
    xm = 0.1             
    gss_initial = 60.8   
    base_name = phase0
  []
  [strain_phase0]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase0
    #    eigenstrain_names = eigenstrain2
  []

  [elasticity_tensor_phase1]
    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
    euler_angle_1 = 0.0
    euler_angle_2 = 45.0
    euler_angle_3 = 0.0
    base_name = phase1
  []
  [stress_phase1]
    type = ComputeMultipleCrystalPlasticityStress
    crystal_plasticity_models = 'trial_xtalpl_phase1'
    rtol = 1e-06
    tan_mod_type = exact
    base_name = phase1
  []
  [trial_xtalpl_phase1]
    type = CrystalPlasticityKalidindiUpdate
    number_slip_systems = 12
    slip_sys_file_name = input_slip_sys.txt
    crystal_lattice_type = FCC
    r = 1.0
    h = 541.5
    t_sat = 109.8
    gss_a = 2.5
    ao = 0.001
    xm = 0.1
    gss_initial = 60.8
    base_name = phase1
  []

  [strain_phase1]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y'
    base_name = phase1
  []

  # Switching functions for each phase
  [h0]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta0
    all_etas = 'eta0 eta1'
    h_name = h0
  []
  [h1]
    type = SwitchingFunctionMultiPhaseMaterial
    phase_etas = eta1
    all_etas = 'eta0 eta1'
    h_name = h1
  []

  # Generate the global stress from the phase stresses
  [global_stress]
    type = MultiPhaseStressMaterial
    phase_base = 'phase0 phase1'
    h = 'h0  h1'
  []
[]

[Kernels]
  [eta0_dt]
    type = TimeDerivative
    variable = eta0
  []
  [eta1_dt]
    type = TimeDerivative
    variable = eta1
  []
  [TensorMechanics]
    displacements = 'disp_x disp_y'
  []
[]

[Postprocessors]
  [phase0_pk2]
    type = ElementAverageValue
    variable = phase0_pk2
  []
  [phase0_fp_xx]
    type = ElementAverageValue
    variable = phase0_fp_xx
  []
  [phase0_e_xx]
    type = ElementAverageValue
    variable = phase0_e_xx
  []
  [phase1_pk2]
    type = ElementAverageValue
    variable = phase1_pk2
  []
  [phase1_fp_xx]
    type = ElementAverageValue
    variable = phase1_fp_xx
  []
  [phase1_e_xx]
    type = ElementAverageValue
    variable = phase1_e_xx
  []
  [vonmises]
    type = ElementAverageValue
    variable = vonmises
    block = 'ANY_BLOCK_ID 0'
  []

  [stress_xx]
    type = ElementAverageValue
    variable = stress_xx
  []
  [phase0_gss]
    type = ElementAverageValue
    variable = phase0_gss
  []
  [phase0_slip_increment]
    type = ElementAverageValue
    variable = phase0_slip_increment
  []
  [phase1_gss]
    type = ElementAverageValue
    variable = phase1_gss
  []
  [phase1_slip_increment]
    type = ElementAverageValue
    variable = phase1_slip_increment
  []
[]

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

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -snes_type'
  petsc_options_value = 'lu            superlu_dist          vinewtonrsls'

  #  petsc_options_iname = '-pc_type -pc_asm_overlap -sub_pc_type -ksp_type -ksp_gmres_restart'
  #  petsc_options_value = ' asm      2              lu            gmres     200'
  l_max_its = 20
  nl_max_its = 10
  nl_rel_tol = 1.0e-8
  nl_abs_tol = 1.0e-9

  end_time = 30

  [TimeStepper]
    type = IterationAdaptiveDT
    dt = 5e-3
    cutback_factor = 0.75
    growth_factor = 1.2
    optimal_iterations = 20
  []
[]

[Outputs]
  exodus = true
  print_linear_residuals = true
  [table]
    type = CSV
    execute_on = timestep_end
  []
  [checkpoint]
    type = Checkpoint
    num_files = 10
    interval = 200
  []
[]

sapitts Jun 24, 2024
Collaborator

@Emmanuel-324 if this last post doesn't refer to the original poster's question, this post needs to be a separate discussion issue, or posted on a discussion thread you have already started for this specific issue. Please follow the guidelines given by @GiudGiud

Emmanuel-324 Jun 24, 2024

Okay, thanks. Sorry about the inconvenience. Any suggestions about the papers?

p-antonioni Jun 26, 2024
Author

@p-antonioni in case you have not already resolved this issue, there is a specific HCP class available in MOOSE, see https://mooseframework.inl.gov/source/materials/crystal_plasticity/CrystalPlasticityHCPDislocationSlipBeyerleinUpdate.html I would recommend using this HCP specific class over modifications to the Kalidindi class if you are still having issues.

Hello @sapitts, thanks for your help and the suggestions.
I looked at the CrystalPlasticityHCPDislocationSlipBeyerleinUpdate documentation, but since I have to employ phenomelogical laws for both shear stress evolution and hardening, I changed to the CrystalPlasticityKalidindiUpdate even for hcp phase. I am able to run simulations after making some modifications to the hardening matrix in relation to the original CrystalPlasticityKalidindiUpdate. However, I am still passing two separate crystal plasticity models—one for prismatic planes and one for basal planes—so I am not taking into account the latent hardening between different modes.

  [stress_alpha]
    type = ComputeMultipleCrystalPlasticityStress
    crystal_plasticity_models = 'alpha_xtalpl_basal alpha_xtalpl_prism'
    tan_mod_type = exact
    maximum_substep_iteration = 10
    block = '2 4 6 8 10 12 14 16 18 20 22 24 26'
  []

Do you believe it would be better to work on CrystalPlasticityHCPDislocationSlipBeyerleinUpdate to make it a phenomenological model or to adapt CrystalPlasticityKalidindiUpdate to get around this problem and maintain a phenomenological model?
Once again, thank you.

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Basal-prismatic slip systems Crystal Plasticity #27763

{{title}}

Replies: 4 comments 15 replies

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

Select a reply

Basal-prismatic slip systems Crystal Plasticity #27763

Replies: 4 comments · 15 replies

GiudGiud May 31, 2024 Collaborator

p-antonioni Jun 10, 2024 Author

sapitts Jun 24, 2024 Collaborator

sapitts Jun 24, 2024 Collaborator

p-antonioni Jun 26, 2024 Author

Replies: 4 comments 15 replies

GiudGiud
May 31, 2024
Collaborator

p-antonioni Jun 10, 2024
Author

sapitts
Jun 24, 2024
Collaborator

sapitts Jun 24, 2024
Collaborator

p-antonioni Jun 26, 2024
Author