ngrilli / c_pfor_am Goto Github PK

[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'
  [../]
[]