Displacement About Axis

Implements a boundary condition that enforces rotationaldisplacement around an axis on a boundary

Description

The boundary condition class DisplacementAboutAxis applies a rotating displacement to the specified mesh surface according to the user defined rotation function. The boundary condition is always applied to the displaced mesh. The rotation function can be given in either radians or in angles, and an axis of rotation can be specified with the axis_origin and axis_direction parameters.

The rotating displacement value at the current node is calculated according to Eq. (1): $(1)var element = document.getElementById("moose-equation-5a4af6d0-2187-4662-b1be-2bca51902867");katex.render(" u_{rotation} = \\boldsymbol{T}^{-1} \\cdot \\boldsymbol{R}_x^{-1} \\cdot \\boldsymbol{R}_y^{-1} \\cdot \\boldsymbol{R}_z \\cdot \\boldsymbol{R}_y \\cdot \\boldsymbol{R}_x \\cdot \\boldsymbol{T}", element, {displayMode:true,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ where $var element = document.getElementById("moose-equation-915370e8-d865-42d0-8392-cf29d75f139e");katex.render("\\boldsymbol{T}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ is the translation matrix for axes of rotation not centered at the coordinate system origin, and $var element = document.getElementById("moose-equation-3d855cdf-04ec-43dd-8e02-d6a62d9880b5");katex.render("\\boldsymbol{R}_x", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ , $var element = document.getElementById("moose-equation-d8c2ee5c-2df2-4444-b37a-614aaa0dc1e9");katex.render("\\boldsymbol{R}_y", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ , and $var element = document.getElementById("moose-equation-01e497a2-7619-4e26-8d58-a8adeeab4c8e");katex.render("\\boldsymbol{R}_z", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ are rotation matrices about the $var element = document.getElementById("moose-equation-2ff4fc87-8d59-4900-9d74-454e07af9dcd");katex.render("\\hat{x}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ , $var element = document.getElementById("moose-equation-c809f8c4-7518-4f76-825f-7f6894cb2207");katex.render("\\hat{y}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ , and $var element = document.getElementById("moose-equation-928c302c-05f5-4000-80ba-82ac4723ead4");katex.render("\\hat{z}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ coordinate system axes, respectively.

Example Input File Syntax

[./top_x]
  type = DisplacementAboutAxis
  boundary = 2
  function = rampConstantAngle
  angle_units = degrees
  axis_origin = '10. 10. 10.'
  axis_direction = '0 -1.0 1.0'
  component = 0
  variable = disp_x
[../]

/opt/civet/build_0/moose/modules/tensor_mechanics/test/tests/torque_reaction/torque_reaction_cylinder.i

# This test uses the DisplacementAboutAxis boundary condition to twist the top
# of a cylinder while the bottom face of the cylinder remains fixed.  The
# TorqueReaction postprocessor is used to calculate the applied torque acting
# on the cylinder at the top face.  This test can be extended, with a new mesh,
# to model a crack in the center of the cylinder face under type III loading.

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  file = cylinder.e
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
  [./disp_z]
  [../]
[]

[GlobalParams]
  volumetric_locking_correction=true
[]

[AuxVariables]
  [./stress_xx]      # stress aux variables are defined for output; this is a way to get integration point variables to the output file
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./saved_x]
  [../]
  [./saved_y]
  [../]
  [./saved_z]
  [../]
[]

[Functions]
  [./rampConstantAngle]
    type = PiecewiseLinear
    x = '0. 1.'
    y = '0. 1.'
    scale_factor = 0.1
  [../]
[]

[Kernels]
  [./TensorMechanics]
    save_in = 'saved_x saved_y saved_z'
    use_displaced_mesh = true
  [../]
[]

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

[BCs]
  [./bottom_x]
    type = DirichletBC
    variable = disp_x
    boundary = 1
    value = 0.0
  [../]
  [./bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1
    value = 0.0
  [../]
  [./bottom_z]
    type = DirichletBC
    variable = disp_z
    boundary = 1
    value = 0.0
  [../]
  [./top_x]
    type = DisplacementAboutAxis
    boundary = 2
    function = rampConstantAngle
    angle_units = degrees
    axis_origin = '10. 10. 10.'
    axis_direction = '0 -1.0 1.0'
    component = 0
    variable = disp_x
  [../]
  [./top_y]
    type = DisplacementAboutAxis
    boundary = 2
    function = rampConstantAngle
    angle_units = degrees
    axis_origin = '10. 10. 10.'
    axis_direction = '0 -1.0 1.0'
    component = 1
    variable = disp_y
  [../]
  [./top_z]
    type = DisplacementAboutAxis
    boundary = 2
    function = rampConstantAngle
    angle_units = degrees
    axis_origin = '10. 10. 10.'
    axis_direction = '0 -1.0 1.0'
    component = 2
    variable = disp_z
  [../]
[] # BCs

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 1
    youngs_modulus = 207000
    poissons_ratio = 0.3
  [../]
  [./strain]
    type = ComputeFiniteStrain
    block = 1
  [../]
  [./elastic_stress]
    type = ComputeFiniteStrainElasticStress
    block = 1
  [../]
[]

[Executioner]

  type = Transient
  #Preconditioned JFNK (default)
  solve_type = 'PJFNK'

  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-ksp_gmres_restart'
  petsc_options_value = '101'

  line_search = 'none'

  l_max_its = 50
  nl_max_its = 20
  nl_abs_tol = 1e-12
  nl_rel_tol = 1e-11
  l_tol = 1e-10

  start_time = 0.0
  dt = 0.25

  end_time = 0.5
[]

[Postprocessors]
  [./torque]
    type = TorqueReaction
    boundary = 2
    reaction_force_variables = 'saved_x saved_y saved_z'
    axis_origin = '10. 10. 10.'
    direction_vector = '0 -1.0 1.0'
  [../]
[]

[Outputs]
  file_base = torque_reaction_cylinder_out
  exodus = true
[]

A Function is required to prescribe the rate of the DisplacementAboutAxis boundary condition application to the mesh. Note that the name of the Function is used as the argument for the function input parameter in the DisplacementAboutAxis block.

[./rampConstantAngle]
  type = PiecewiseLinear
  x = '0. 1.'
  y = '0. 1.'
  scale_factor = 0.1
[../]

/opt/civet/build_0/moose/modules/tensor_mechanics/test/tests/torque_reaction/torque_reaction_cylinder.i

# This test uses the DisplacementAboutAxis boundary condition to twist the top
# of a cylinder while the bottom face of the cylinder remains fixed.  The
# TorqueReaction postprocessor is used to calculate the applied torque acting
# on the cylinder at the top face.  This test can be extended, with a new mesh,
# to model a crack in the center of the cylinder face under type III loading.

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  file = cylinder.e
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
  [./disp_z]
  [../]
[]

[GlobalParams]
  volumetric_locking_correction=true
[]

[AuxVariables]
  [./stress_xx]      # stress aux variables are defined for output; this is a way to get integration point variables to the output file
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./saved_x]
  [../]
  [./saved_y]
  [../]
  [./saved_z]
  [../]
[]

[Functions]
  [./rampConstantAngle]
    type = PiecewiseLinear
    x = '0. 1.'
    y = '0. 1.'
    scale_factor = 0.1
  [../]
[]

[Kernels]
  [./TensorMechanics]
    save_in = 'saved_x saved_y saved_z'
    use_displaced_mesh = true
  [../]
[]

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

[BCs]
  [./bottom_x]
    type = DirichletBC
    variable = disp_x
    boundary = 1
    value = 0.0
  [../]
  [./bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1
    value = 0.0
  [../]
  [./bottom_z]
    type = DirichletBC
    variable = disp_z
    boundary = 1
    value = 0.0
  [../]
  [./top_x]
    type = DisplacementAboutAxis
    boundary = 2
    function = rampConstantAngle
    angle_units = degrees
    axis_origin = '10. 10. 10.'
    axis_direction = '0 -1.0 1.0'
    component = 0
    variable = disp_x
  [../]
  [./top_y]
    type = DisplacementAboutAxis
    boundary = 2
    function = rampConstantAngle
    angle_units = degrees
    axis_origin = '10. 10. 10.'
    axis_direction = '0 -1.0 1.0'
    component = 1
    variable = disp_y
  [../]
  [./top_z]
    type = DisplacementAboutAxis
    boundary = 2
    function = rampConstantAngle
    angle_units = degrees
    axis_origin = '10. 10. 10.'
    axis_direction = '0 -1.0 1.0'
    component = 2
    variable = disp_z
  [../]
[] # BCs

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 1
    youngs_modulus = 207000
    poissons_ratio = 0.3
  [../]
  [./strain]
    type = ComputeFiniteStrain
    block = 1
  [../]
  [./elastic_stress]
    type = ComputeFiniteStrainElasticStress
    block = 1
  [../]
[]

[Executioner]

  type = Transient
  #Preconditioned JFNK (default)
  solve_type = 'PJFNK'

  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-ksp_gmres_restart'
  petsc_options_value = '101'

  line_search = 'none'

  l_max_its = 50
  nl_max_its = 20
  nl_abs_tol = 1e-12
  nl_rel_tol = 1e-11
  l_tol = 1e-10

  start_time = 0.0
  dt = 0.25

  end_time = 0.5
[]

[Postprocessors]
  [./torque]
    type = TorqueReaction
    boundary = 2
    reaction_force_variables = 'saved_x saved_y saved_z'
    axis_origin = '10. 10. 10.'
    direction_vector = '0 -1.0 1.0'
  [../]
[]

[Outputs]
  file_base = torque_reaction_cylinder_out
  exodus = true
[]

Input Parameters

angle_unitsThe units of the angle of rotation. Choices are:degrees radians
C++ Type:MooseEnum
Options:degrees radians
Description:The units of the angle of rotation. Choices are:degrees radians
axis_directionDirection of the axis of rotation
C++ Type:libMesh::VectorValue
Options:
Description:Direction of the axis of rotation
axis_originOrigin of the axis of rotation
C++ Type:libMesh::VectorValue
Options:
Description:Origin of the axis of rotation
boundaryThe list of boundary IDs from the mesh where this boundary condition applies
C++ Type:std::vector
Options:
Description:The list of boundary IDs from the mesh where this boundary condition applies
componentThe component for the rotational displacement
C++ Type:int
Options:
Description:The component for the rotational displacement
functionThe function providing the angle of rotation.
C++ Type:FunctionName
Options:
Description:The function providing the angle of rotation.
variableThe name of the variable that this boundary condition applies to
C++ Type:NonlinearVariableName
Options:
Description:The name of the variable that this boundary condition applies to

Required Parameters

diag_save_inThe name of auxiliary variables to save this BC's diagonal jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector
Options:
Description:The name of auxiliary variables to save this BC's diagonal jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
displacementsThe displacements
C++ Type:std::vector
Options:
Description:The displacements
presetTrueWhether or not to preset the BC (apply the value before the solve begins). Note that the default value of this parameter is handled by the use_legacy_dirichlet_bc parameter on the MooseApp.
Default:True
C++ Type:bool
Options:
Description:Whether or not to preset the BC (apply the value before the solve begins). Note that the default value of this parameter is handled by the use_legacy_dirichlet_bc parameter on the MooseApp.
save_inThe name of auxiliary variables to save this BC's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector
Options:
Description:The name of auxiliary variables to save this BC's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Options:
Description:Determines whether this object is calculated using an implicit or explicit form
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Options:
Description:The seed for the master random number generator
use_displaced_meshTrueWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:True
C++ Type:bool
Options:
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

extra_matrix_tagsThe extra tags for the matrices this Kernel should fill
C++ Type:std::vector
Options:
Description:The extra tags for the matrices this Kernel should fill
extra_vector_tagsThe extra tags for the vectors this Kernel should fill
C++ Type:std::vector
Options:
Description:The extra tags for the vectors this Kernel should fill
matrix_tagssystem timeThe tag for the matrices this Kernel should fill
Default:system time
C++ Type:MultiMooseEnum
Options:nontime system time
Description:The tag for the matrices this Kernel should fill
vector_tagsresidualThe tag for the vectors this Kernel should fill
Default:residual
C++ Type:MultiMooseEnum
Options:nontime time residual
Description:The tag for the vectors this Kernel should fill

Tagging Parameters

Input Files

modules/tensor_mechanics/test/tests/torque_reaction/disp_about_axis_axial_motion.i
modules/tensor_mechanics/test/tests/torque_reaction/torque_reaction_cylinder.i
modules/tensor_mechanics/test/tests/torque_reaction/disp_about_axis_errors.i

modules/tensor_mechanics/test/tests/torque_reaction/disp_about_axis_axial_motion.i

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  type = GeneratedMesh
  dim = 3
  nx = 1
  ny = 1
  nz = 1
[]


[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
  [./disp_z]
  [../]
[]

[GlobalParams]
  volumetric_locking_correction=true
[]

[AuxVariables]
  [./stress_xx]      # stress aux variables are defined for output; this is a way to get integration point variables to the output file
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_yy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_zz]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./vonmises]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]


[Functions]
  [./rampConstant]
    type = PiecewiseLinear
    x = '0. 1.'
    y = '0. 1.'
    scale_factor = 1.
  [../]
[]

[Kernels]
  [./TensorMechanics]
    use_displaced_mesh = true
  [../]
[]

[AuxKernels]
  [./stress_xx]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_i = 0
    index_j = 0
    execute_on = timestep_end
  [../]
  [./stress_yy]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_yy
    index_i = 1
    index_j = 1
    execute_on = timestep_end
  [../]
  [./stress_zz]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_zz
    index_i = 2
    index_j = 2
    execute_on = timestep_end
  [../]
  [./vonmises]
    type =  RankTwoScalarAux
    rank_two_tensor = stress
    variable = vonmises
    scalar_type = VonMisesStress
    execute_on = timestep_end
  [../]
[]

[BCs]

  [./bottom_x]
    type = DirichletBC
    variable = disp_x
    boundary = bottom
    value = 0.0
  [../]

  [./bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0.0
  [../]

  [./bottom_z]
    type = DirichletBC
    variable = disp_z
    boundary = bottom
    value = 0.0
  [../]

  [./top_x]
    type = DisplacementAboutAxis
    boundary = top
    function = rampConstant
    angle_units = degrees
    axis_origin = '0. 0. 0.'
    axis_direction = '0. 0. 1.'
    component = 0
    variable = disp_x
  [../]

  [./top_y]
    type = DisplacementAboutAxis
    boundary = top
    function = rampConstant
    angle_units = degrees
    axis_origin = '0. 0. 0.'
    axis_direction = '0. 0. 1.'
    component = 1
    variable = disp_y
  [../]

# Because want to keep the rotation fixed about the z axis,
# do not apply a DisplacementAboutAxis BC on the disp_z
[] # BCs

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 0
    youngs_modulus = 207000
    poissons_ratio = 0.3
  [../]
  [./strain]
    type = ComputeFiniteStrain
    block = 0
  [../]
  [./elastic_stress]
    type = ComputeFiniteStrainElasticStress
    block = 0
  [../]
[]


[Executioner]
  type = Transient

  solve_type = 'PJFNK'

  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-ksp_gmres_restart'
  petsc_options_value = '101'

  line_search = 'none'

  l_max_its = 30
  nl_max_its = 20
  nl_rel_tol = 1e-10
  nl_abs_tol = 1e-12
  l_tol = 1e-8

  start_time = 0.0
  dt = 0.1
  dtmin = 0.1 # die instead of cutting the timestep

  end_time = 0.5
[]

[Postprocessors]
  [./_dt]
    type = TimestepSize
  [../]
[]


[Outputs]
  file_base = disp_about_axis_axial_motion_out
  exodus = true
[]

modules/tensor_mechanics/test/tests/torque_reaction/torque_reaction_cylinder.i

# This test uses the DisplacementAboutAxis boundary condition to twist the top
# of a cylinder while the bottom face of the cylinder remains fixed.  The
# TorqueReaction postprocessor is used to calculate the applied torque acting
# on the cylinder at the top face.  This test can be extended, with a new mesh,
# to model a crack in the center of the cylinder face under type III loading.

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  file = cylinder.e
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
  [./disp_z]
  [../]
[]

[GlobalParams]
  volumetric_locking_correction=true
[]

[AuxVariables]
  [./stress_xx]      # stress aux variables are defined for output; this is a way to get integration point variables to the output file
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./saved_x]
  [../]
  [./saved_y]
  [../]
  [./saved_z]
  [../]
[]

[Functions]
  [./rampConstantAngle]
    type = PiecewiseLinear
    x = '0. 1.'
    y = '0. 1.'
    scale_factor = 0.1
  [../]
[]

[Kernels]
  [./TensorMechanics]
    save_in = 'saved_x saved_y saved_z'
    use_displaced_mesh = true
  [../]
[]

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

[BCs]
  [./bottom_x]
    type = DirichletBC
    variable = disp_x
    boundary = 1
    value = 0.0
  [../]
  [./bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1
    value = 0.0
  [../]
  [./bottom_z]
    type = DirichletBC
    variable = disp_z
    boundary = 1
    value = 0.0
  [../]
  [./top_x]
    type = DisplacementAboutAxis
    boundary = 2
    function = rampConstantAngle
    angle_units = degrees
    axis_origin = '10. 10. 10.'
    axis_direction = '0 -1.0 1.0'
    component = 0
    variable = disp_x
  [../]
  [./top_y]
    type = DisplacementAboutAxis
    boundary = 2
    function = rampConstantAngle
    angle_units = degrees
    axis_origin = '10. 10. 10.'
    axis_direction = '0 -1.0 1.0'
    component = 1
    variable = disp_y
  [../]
  [./top_z]
    type = DisplacementAboutAxis
    boundary = 2
    function = rampConstantAngle
    angle_units = degrees
    axis_origin = '10. 10. 10.'
    axis_direction = '0 -1.0 1.0'
    component = 2
    variable = disp_z
  [../]
[] # BCs

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 1
    youngs_modulus = 207000
    poissons_ratio = 0.3
  [../]
  [./strain]
    type = ComputeFiniteStrain
    block = 1
  [../]
  [./elastic_stress]
    type = ComputeFiniteStrainElasticStress
    block = 1
  [../]
[]

[Executioner]

  type = Transient
  #Preconditioned JFNK (default)
  solve_type = 'PJFNK'

  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-ksp_gmres_restart'
  petsc_options_value = '101'

  line_search = 'none'

  l_max_its = 50
  nl_max_its = 20
  nl_abs_tol = 1e-12
  nl_rel_tol = 1e-11
  l_tol = 1e-10

  start_time = 0.0
  dt = 0.25

  end_time = 0.5
[]

[Postprocessors]
  [./torque]
    type = TorqueReaction
    boundary = 2
    reaction_force_variables = 'saved_x saved_y saved_z'
    axis_origin = '10. 10. 10.'
    direction_vector = '0 -1.0 1.0'
  [../]
[]

[Outputs]
  file_base = torque_reaction_cylinder_out
  exodus = true
[]

modules/tensor_mechanics/test/tests/torque_reaction/disp_about_axis_errors.i

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  type = GeneratedMesh
  dim = 3
  nx = 1
  ny = 1
  nz = 1
[]


[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
  [./disp_z]
  [../]
[]


[AuxVariables]
  [./stress_xx]      # stress aux variables are defined for output; this is a way to get integration point variables to the output file
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_yy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_zz]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./vonmises]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]


[Functions]
  [./rampConstant]
    type = PiecewiseLinear
    x = '0. 1.'
    y = '0. 1.'
    scale_factor = 1.
  [../]
[]

[Kernels]
  [./TensorMechanics]
    use_displaced_mesh = true
  [../]
[]

[AuxKernels]
  [./stress_xx]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_i = 0
    index_j = 0
    execute_on = timestep_end
  [../]
  [./stress_yy]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_yy
    index_i = 1
    index_j = 1
    execute_on = timestep_end
  [../]
  [./stress_zz]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_zz
    index_i = 2
    index_j = 2
    execute_on = timestep_end
  [../]
  [./vonmises]
    type =  RankTwoScalarAux
    rank_two_tensor = stress
    variable = vonmises
    scalar_type = VonMisesStress
    execute_on = timestep_end
  [../]
[]

[BCs]

  [./bottom_x]
    type = DirichletBC
    variable = disp_x
    boundary = bottom
    value = 0.0
  [../]

  [./bottom_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0.0
  [../]

  [./bottom_z]
    type = DirichletBC
    variable = disp_z
    boundary = bottom
    value = 0.0
  [../]

  [./top_x]
    type = DisplacementAboutAxis
    boundary = top
    function = rampConstant
    angle_units = degrees
    axis_origin = '0. 0. 0.'
    axis_direction = '0. 0. 1.'
    component = 0
    variable = disp_x
  [../]

  [./top_y]
    type = DisplacementAboutAxis
    boundary = top
    function = rampConstant
    angle_units = degrees
    axis_origin = '0. 0. 0.'
    variable = disp_y
  [../]

[] # BCs

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 0
    youngs_modulus = 207000
    poissons_ratio = 0.3
  [../]
  [./strain]
    type = ComputeFiniteStrain
    block = 0
  [../]
  [./elastic_stress]
    type = ComputeFiniteStrainElasticStress
    block = 0
  [../]
[]


[Executioner]
  type = Transient

  solve_type = 'PJFNK'

  petsc_options = '-snes_ksp_ew'
  petsc_options_iname = '-ksp_gmres_restart'
  petsc_options_value = '101'

  line_search = 'none'

  l_max_its = 30
  nl_max_its = 20
  nl_rel_tol = 1e-10
  nl_abs_tol = 1e-12
  l_tol = 1e-8

  start_time = 0.0
  dt = 0.1
  dtmin = 0.1 # die instead of cutting the timestep

  end_time = 0.5
[]


[Outputs]
  file_base = disp_about_axis_errors_out
  exodus = true
[]

Description
Example Input File Syntax
Input Parameters
Input Files

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools