[Mesh]
  [Cube]
  type = FileMeshGenerator
  file = DRAGen_RVE_16G.e
 # construct_side_list_from_node_list = true
  []
  [Normal_SS]
  type = AllSideSetsByNormalsGenerator
  input = Cube
  show_info = True
  []
[]


[GlobalParams]
  displacements = 'u_x u_y u_z'
[]

[Variables]
  [./u_x]
    order = FIRST
    family = LAGRANGE
  [../]
  [./u_y]
    order = FIRST
    family = LAGRANGE
  [../]
  [./u_z]
    order = FIRST
    family = LAGRANGE
  [../]
[]


[Modules/TensorMechanics/Master]
  [./block1]
    strain = FINITE
    add_variables = true
    generate_output = 'stress_zz stress_xx stress_xy stress_xz stress_yy stress_yz strain_xx strain_xy strain_xz strain_yy strain_yz strain_zz vonmises_stress'
  [../]
[]



[Materials]
  [./elasticity_tensor]
    type = ComputeElasticityTensorCP
    C_ijkl = '2.810e5 1.324e5 1.324e5 2.810e5 1.324e5 2.810e5 0.754e5 0.754e5 0.754e5'
    fill_method = symmetric9
    euler_angle_1 = 173.0
    euler_angle_2 = 67.0
    euler_angle_3 = 266.0
  [../]
[]

[BCs]
  [./Periodic]
    [./all]
      variable = 'u_x u_y u_z'
      auto_direction = 'x y'
    [../]
  [../]

  [./top]
    type = FunctionDirichletBC
    variable = 'u_z'
    boundary = 'front'
    function = 0.001*t
  [../]
  [./bot]
    type = FunctionDirichletBC
    variable = 'u_z'
    boundary = 'back'
    function = -0.001*t
  [../]
[]


[Executioner]

  type = Transient
  solve_type = 'PJFNK'
  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  petsc_options_value = 'hypre    boomeramg          31'
  line_search = 'none'
  l_max_its = 30
  nl_max_its = 50
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-6
  l_tol = 1e-8

  start_time = 0.0
  end_time = 0.125
  dt = 0.015
[]

[Outputs]
  execute_on = 'timestep_end'
  file_base = out_TM
  print_linear_residuals = false
  exodus = true
[]