I am trying to add the plastic_strain auxvariable to my calculations and am getting the following error:
Does the FiniteStrainCrystalPlasticityThermal material not include plastic_strain?
[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'
[../]
[]