ngrilli / c_pfor_am

Crystal plasticity for additive manufacturing
GNU Lesser General Public License v2.1
50 stars 26 forks source link

Plastic Strain Variable #3

Closed lukuna closed 2 years ago

lukuna commented 2 years ago

I am trying to add the plastic_strain auxvariable to my calculations and am getting the following error:

Material property 'plastic_strain', requested by 'plastic_xx' is not defined on block 1

Does the FiniteStrainCrystalPlasticityThermal material not include plastic_strain?

A snippet of the input and relevant blocks is included below:


[AuxKernels]
  [./plastic_xx]
    type = RankTwoAux
    rank_two_tensor = plastic_strain
    variable = plastic_xx
    index_i = 0
    index_j = 0
    block = '1'
  [../]
  [./plastic_yy]
    type = RankTwoAux
    rank_two_tensor = plastic_strain
    variable = plastic_yy
    index_i = 1
    index_j = 1
    block = '1'
  [../]
  [./plastic_zz]
    type = RankTwoAux
    rank_two_tensor = plastic_strain
    variable = plastic_zz
    index_i = 2
    index_j = 2
    block = '1'
  [../]
[]

[AuxVariables]
  [./plastic_yy]
    order = CONSTANT
    family = MONOMIAL
     block = '1'
  [../]
  [./plastic_xx]
    order = CONSTANT
    family = MONOMIAL
     block = '1'
  [../]
  [./plastic_zz]
    order = CONSTANT
    family = MONOMIAL
     block = '1'
  [../]
[]

[Materials]
  [./crysp]
    type = FiniteStrainCrystalPlasticityThermal
    block = '1'
    slip_sys_file_name = input_slip_sys.txt # no need to normalize vectors
    nss = 12 #Number of slip systems
    num_slip_sys_flowrate_props = 2 #Number of flow rate properties in a slip system
    flowprops = '1 12 0.0005 0.1' # slip rate equations parameters
    hprops = '1.0 3839.0 180.0 302.0 2.5' # hardening properties
    gprops = '1 12 180.0' # initial values of slip system resistances (start_slip_sys, end_slip_sys, value)
    tan_mod_type = none
    thermal_expansion = '44.73e-6'
    reference_temperature = '403.0'
    temp = temp
    slip_incr_tol = 0.001
    maxiter = 200
    maxitergss = 200
    maximum_substep_iteration = 5
    gen_random_stress_flag = true
    rtol = '1.0e-4'
    abs_tol = '1.0e-4'
    dCRSS_dT_A = 0.53
    dCRSS_dT_B = 0.47
    dCRSS_dT_C = 0.008
    dCTE_dT = '0.01011e-6'
  [../]

  [./elasticity_tensor]
    type = ComputeElasticityTensorMelting
    C_ijkl = '2.046e5 1.377e5 1.377e5 2.046e5 1.377e5 2.046e5 1.262e5 1.262e5 1.262e5'
    fill_method = symmetric9
    read_prop_user_object = prop_read
    temp = temp
    dC11_dT = 0.0004415
    dC12_dT = 0.0003275
    dC44_dT = 0.0004103
    residual_stiffness = 0.001
    temperature_read_user_object = temperature_read
    temperature_time_step = 95.0
    reference_temperature = '403.0'
    activate_elems = true
    block = '1'
  [../]
  [./strain]
    type = ComputeFiniteStrain
    displacements = 'disp_x disp_y disp_z'
    block = '1'
  [../]
  [./dummy_mat_inactive]
    type = GenericConstantMaterial
    prop_names = 'dummy_mat'
    prop_values = '0.0'
    block = '2'
  [../]
[]
ngrilli commented 2 years ago

@lukuna there is no material property called "plastic_strain" in FiniteStrainCrystalPlasticityThermal, to see the list of all material properties related to plasticity, you need to check in the parent class, FiniteStrainCrystalPlasticity.

You can try with: rank_two_tensor = fp

which is the plastic deformation gradient.

Trust this helps, Nicolò