Generalized Plane Strain Action System

Set up the GeneralizedPlaneStrain environment

Description

This action sets up a generalized plane strain model. A detailed description of the formulation is available on the generalized plane strain page.

warning:For 2D and 1D Simulations

GeneralizedPlaneStrainAction supports 1D axisymmetric or 2D generalized plane strain cases. For 1D axisymmetric and 2D Cartesian cases in the x-y plane, the number of displacement variables must be one or two, respectively.

For 2D generalized plane strain cases in the x-z or y-z planes, the number of displacement variables must be three.

Constructed MooseObjects

The GeneralizedPlaneStrain Action constructs the objects and scalar out-of-plane strain variable required for a generalized plane strain simulation. When the scalar variable named by "scalar_out_of_plane_strain" does not exist, the action creates a FIRST order nonlinear scalar variable in nonlinear system nl0. If that scalar variable already exists in nl0, the action reuses it.

Table 1: Legacy objects created by GeneralizedPlaneStrainAction

Object	Purpose
Generalized Plane Strain Off Diagonal	Couples in-plane displacement variables and the scalar out-of-plane strain variable in the off-diagonal Jacobian.
Generalized Plane Strain Scalar Kernel	Assembles the scalar out-of-plane equilibrium residual.
Generalized Plane Strain User Object	Computes residual and diagonal Jacobian data for the scalar kernel.

When "use_automatic_differentiation" is true, the action creates exactly one AD object and skips the legacy UserObject, ScalarKernel, and off-diagonal Kernel objects. ADKernelScalarBase owns the scalar residual assembly in this mode; when "use_displaced_mesh" is true, the AD kernel uses displaced-mesh quadrature weights for the scalar equation.

Table 2: AD objects created by GeneralizedPlaneStrainAction

Object	Purpose
ADGeneralizedPlaneStrain	Assembles the scalar out-of-plane equilibrium residual and AD couplings.

The scalar residual is $var element = document.getElementById("moose-equation-f729a635-9a26-406e-a8e3-f8635b5adff1");katex.render(" R_s = \\int_A \\left( \\sigma_{oo} + p_o \\right) dA = 0,", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ where $var element = document.getElementById("moose-equation-0b00ba33-bd77-4e2d-b3d9-d38809512d1f");katex.render("p_o", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is positive when applied toward the body. In Cartesian coordinates, the out-of-plane x, y, and z directions use $var element = document.getElementById("moose-equation-0c77f851-8c26-4eb6-8ac2-e3675f77e8f5");katex.render("\\sigma_{xx}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , $var element = document.getElementById("moose-equation-bfdd317e-f1bb-409c-b119-29051e72ff85");katex.render("\\sigma_{yy}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , and $var element = document.getElementById("moose-equation-72bb79e1-7600-4bcf-82f4-413cd7144693");katex.render("\\sigma_{zz}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , respectively. In RZ axisymmetry, the out-of-plane direction is the axial y direction and the residual uses axisymmetric coordinate weighting.

Example Input Syntax

Subblocks

The subblocks of the GeneralizedPlaneStrain action trigger object construction. If a generalized plane strain model is applied to the whole simulation domain, use a single subblock.

[./GeneralizedPlaneStrain]
  [./gps]
    use_displaced_mesh = true
    displacements = 'disp_x disp_y'
    scalar_out_of_plane_strain = scalar_strain_zz
    out_of_plane_pressure_function = traction_function
    pressure_factor = 1e5
  [../]
[../]

If different mesh subdomains use different generalized plane strain models, use multiple subblocks with subdomain restrictions.

[./GeneralizedPlaneStrain]
  [./gps1]
    use_displaced_mesh = true
    displacements = 'disp_x disp_y'
    scalar_out_of_plane_strain = scalar_strain_zz1
    block = '1'
  [../]
  [./gps2]
    use_displaced_mesh = true
    displacements = 'disp_x disp_y'
    scalar_out_of_plane_strain = scalar_strain_zz2
    block = '2'
  [../]
[../]

An example of using generalized plane strain through the Solid Mechanics QuasiStatic physics block with an "out_of_plane_direction" different from the default $var element = document.getElementById("moose-equation-ed5a1191-5385-42db-b5b2-2162a8432268");katex.render("z", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ direction is given by:

[./generalized_plane_strain]
  block = 1
  strain = SMALL
  scalar_out_of_plane_strain = scalar_strain_yy
  out_of_plane_direction = y
  planar_formulation = GENERALIZED_PLANE_STRAIN
  eigenstrain_names = 'eigenstrain'
  generate_output = 'stress_xx stress_xz stress_yy stress_zz strain_xx strain_xz strain_yy strain_zz'
[../]

Parameters supplied at the [Physics/SolidMechanics/GeneralizedPlaneStrain] level act as defaults for the QuasiStatic Physics subblocks.

The following test-spec block shows the standard small-strain generalized plane strain input run with automatic differentiation.

  [generalized_plane_strain_small_ad]
    type = 'Exodiff'
    input = 'generalized_plane_strain_small.i'
    exodiff = 'generalized_plane_strain_small_out.e'
    custom_cmp = 'generalized.exodiff'
    prereq = generalized_plane_strain_small
    cli_args = 'Physics/SolidMechanics/QuasiStatic/all/use_automatic_differentiation=true '
               'Postprocessors/react_z/type=ADMaterialTensorIntegral '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    requirement = 'The system shall support an automatic-differentiation generalized plane strain mechanics solution'
    design = 'ADComputePlaneSmallStrain.md'
  []

The following standalone action input lets the action create the missing scalar out-of-plane strain variable in AD mode.

[Physics<<<{"href": "../../index.html"}>>>]
  [SolidMechanics<<<{"href": "../index.html"}>>>]
    [GeneralizedPlaneStrain<<<{"href": "index.html"}>>>]
      [gps]
        use_automatic_differentiation<<<{"description": "Use automatic differentiation to assemble the generalized plane strain equation and its coupling terms"}>>> = true
        use_displaced_mesh<<<{"description": "Whether to use displaced mesh"}>>> = true
        displacements<<<{"description": "The displacement variables"}>>> = 'disp_x disp_y'
        scalar_out_of_plane_strain<<<{"description": "Scalar variable for the out-of-plane strain (in y direction for 1D Axisymmetric or in z direction for 2D Cartesian problems)"}>>> = scalar_strain_zz
        out_of_plane_pressure_function<<<{"description": "Function used to prescribe pressure (applied toward the body) in the out-of-plane direction (y for 1D Axisymmetric or z for 2D Cartesian problems)"}>>> = traction_function
        pressure_factor<<<{"description": "Scale factor applied to prescribed out-of-plane pressure (both material and function)"}>>> = 1e5
      []
    []
  []
[]

Input Parameters

displacementsThe displacement variables
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:The displacement variables
scalar_out_of_plane_strainScalar variable for the out-of-plane strain (in y direction for 1D Axisymmetric or in z direction for 2D Cartesian problems)
C++ Type:VariableName
Unit:(no unit assumed)
Controllable:No
Description:Scalar variable for the out-of-plane strain (in y direction for 1D Axisymmetric or in z direction for 2D Cartesian problems)

Required Parameters

absolute_value_vector_tagsThe tag names for extra vectors that the absolute value of the residual should be accumulated into
C++ Type:std::vector<TagName>
Controllable:No
Description:The tag names for extra vectors that the absolute value of the residual should be accumulated into
active__all__ If specified only the blocks named will be visited and made active
Default:__all__
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified only the blocks named will be visited and made active
base_nameMaterial property base name
C++ Type:std::string
Controllable:No
Description:Material property base name
blockThe list of ids of the blocks (subdomain) that the GeneralizedPlaneStrain kernels will be applied to
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of ids of the blocks (subdomain) that the GeneralizedPlaneStrain kernels will be applied to
extra_vector_tagsThe tag names for extra vectors that residual data should be saved into
C++ Type:std::vector<TagName>
Controllable:No
Description:The tag names for extra vectors that residual data should be saved into
inactiveIf specified blocks matching these identifiers will be skipped.
C++ Type:std::vector<std::string>
Controllable:No
Description:If specified blocks matching these identifiers will be skipped.
out_of_plane_directionzThe direction of the out-of-plane strain.
Default:z
C++ Type:MooseEnum
Options:x, y, z
Controllable:No
Description:The direction of the out-of-plane strain.
out_of_plane_pressure_functionFunction used to prescribe pressure (applied toward the body) in the out-of-plane direction (y for 1D Axisymmetric or z for 2D Cartesian problems)
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Function used to prescribe pressure (applied toward the body) in the out-of-plane direction (y for 1D Axisymmetric or z for 2D Cartesian problems)
out_of_plane_pressure_material0Material used to prescribe pressure (applied toward the body) in the out-of-plane direction
Default:0
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:Material used to prescribe pressure (applied toward the body) in the out-of-plane direction
pressure_factorScale factor applied to prescribed out-of-plane pressure (both material and function)
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scale factor applied to prescribed out-of-plane pressure (both material and function)
temperatureThe temperature variable
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:The temperature variable
use_automatic_differentiationFalseUse automatic differentiation to assemble the generalized plane strain equation and its coupling terms
Default:False
C++ Type:bool
Controllable:No
Description:Use automatic differentiation to assemble the generalized plane strain equation and its coupling terms
use_displaced_meshFalseWhether to use displaced mesh
Default:False
C++ Type:bool
Controllable:No
Description:Whether to use displaced mesh

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.

Advanced Parameters

Associated Actions

Available Actions

Solid Mechanics App
GeneralizedPlaneStrainActionSet up the GeneralizedPlaneStrain environment

(modules/solid_mechanics/test/tests/generalized_plane_strain/out_of_plane_pressure.i)

# Tests for application of out-of-plane pressure in generalized plane strain.

[Mesh]
  [./square]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 2
    ny = 2
  [../]
  displacements = 'disp_x disp_y'
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
  [./scalar_strain_zz]
    order = FIRST
    family = SCALAR
  [../]
[]

[AuxVariables]
  [./saved_x]
    order = FIRST
    family = LAGRANGE
  [../]
  [./saved_y]
    order = FIRST
    family = LAGRANGE
  [../]

  [./stress_xx]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_xy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_yy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_zz]
    order = CONSTANT
    family = MONOMIAL
  [../]

  [./strain_xx]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./strain_xy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./strain_yy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./strain_zz]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[Postprocessors]
  [./react_z]
    type = MaterialTensorIntegral
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
  [../]
[]

[Physics]
  [SolidMechanics]
    [./GeneralizedPlaneStrain]
      [./gps]
        use_displaced_mesh = true
        displacements = 'disp_x disp_y'
        scalar_out_of_plane_strain = scalar_strain_zz
        out_of_plane_pressure_function = traction_function
        pressure_factor = 1e5
      [../]
    [../]
  [../]
[]

[Kernels]
  [SolidMechanics]
    use_displaced_mesh = false
    displacements = 'disp_x disp_y'
    save_in = 'saved_x saved_y'
  [../]
[]

[AuxKernels]
  [./stress_xx]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_i = 0
    index_j = 0
  [../]
  [./stress_xy]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xy
    index_i = 0
    index_j = 1
  [../]
  [./stress_yy]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_yy
    index_i = 1
    index_j = 1
  [../]
  [./stress_zz]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_zz
    index_i = 2
    index_j = 2
  [../]

  [./strain_xx]
    type = RankTwoAux
    rank_two_tensor = total_strain
    variable = strain_xx
    index_i = 0
    index_j = 0
  [../]
  [./strain_xy]
    type = RankTwoAux
    rank_two_tensor = total_strain
    variable = strain_xy
    index_i = 0
    index_j = 1
  [../]
  [./strain_yy]
    type = RankTwoAux
    rank_two_tensor = total_strain
    variable = strain_yy
    index_i = 1
    index_j = 1
  [../]
  [./strain_zz]
    type = RankTwoAux
    rank_two_tensor = total_strain
    variable = strain_zz
    index_i = 2
    index_j = 2
  [../]
[]

[Functions]
  [./traction_function]
    type = PiecewiseLinear
    x = '0  2'
    y = '0  1'
  [../]
[]

[BCs]
  [./leftx]
    type = DirichletBC
    boundary = 3
    variable = disp_x
    value = 0.0
  [../]
  [./bottomy]
    type = DirichletBC
    boundary = 0
    variable = disp_y
    value = 0.0
  [../]
[]

[Materials]
  [./elastic_tensor]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = 0.3
    youngs_modulus = 1e6
  [../]
  [./strain]
    type = ComputePlaneSmallStrain
    displacements = 'disp_x disp_y'
    scalar_out_of_plane_strain = scalar_strain_zz
  [../]
  [./stress]
    type = ComputeLinearElasticStress
  [../]
  # This material is not used for anything in the base verison of this test,
  # but is used in a variant of the test with cli_args
  [./traction_material]
    type = GenericFunctionMaterial
    prop_names = traction_material
    prop_values = traction_function
  [../]
[]

[Executioner]
  type = Transient

  solve_type = PJFNK
  line_search = none

# controls for linear iterations
  l_max_its = 100
  l_tol = 1e-4

# controls for nonlinear iterations
  nl_max_its = 15
  nl_rel_tol = 1e-14
  nl_abs_tol = 1e-11

# time control
  start_time = 0.0
  dt = 1.0
  dtmin = 1.0
  end_time = 2.0
  num_steps = 5000
[]

[Outputs]
  exodus = true
[]

(modules/solid_mechanics/test/tests/generalized_plane_strain/generalized_plane_strain_squares.i)

[Mesh]
  file = 2squares.e
  displacements = 'disp_x disp_y'
[]

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

  [./scalar_strain_zz1]
    order = FIRST
    family = SCALAR
  [../]
  [./scalar_strain_zz2]
    order = FIRST
    family = SCALAR
  [../]
[]

[AuxVariables]
  [./temp]
    order = FIRST
    family = LAGRANGE
  [../]
  [./saved_x]
    order = FIRST
    family = LAGRANGE
  [../]
  [./saved_y]
    order = FIRST
    family = LAGRANGE
  [../]

  [./stress_xx]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_xy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_yy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_zz]
    order = CONSTANT
    family = MONOMIAL
  [../]

  [./strain_xx]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./strain_xy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./strain_yy]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./aux_strain_zz]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[Postprocessors]
  [./react_z1]
    type = MaterialTensorIntegral
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
    block = 1
  [../]
  [./react_z2]
    type = MaterialTensorIntegral
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
    block = 2
  [../]
[]

[Physics]
  [SolidMechanics]
    [./GeneralizedPlaneStrain]
      [./gps1]
        use_displaced_mesh = true
        displacements = 'disp_x disp_y'
        scalar_out_of_plane_strain = scalar_strain_zz1
        block = '1'
      [../]
      [./gps2]
        use_displaced_mesh = true
        displacements = 'disp_x disp_y'
        scalar_out_of_plane_strain = scalar_strain_zz2
        block = '2'
      [../]
    [../]
  [../]
[]

[Kernels]
  [SolidMechanics]
    use_displaced_mesh = false
    displacements = 'disp_x disp_y'
    temperature = temp
    save_in = 'saved_x saved_y'
    block = '1 2'
  [../]
[]

[AuxKernels]
  [./tempfuncaux]
    type = FunctionAux
    variable = temp
    function = tempfunc
    use_displaced_mesh = false
  [../]
  [./stress_xx]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_i = 0
    index_j = 0
  [../]
  [./stress_xy]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xy
    index_i = 0
    index_j = 1
  [../]
  [./stress_yy]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_yy
    index_i = 1
    index_j = 1
  [../]
  [./stress_zz]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_zz
    index_i = 2
    index_j = 2
  [../]

  [./strain_xx]
    type = RankTwoAux
    rank_two_tensor = total_strain
    variable = strain_xx
    index_i = 0
    index_j = 0
  [../]
  [./strain_xy]
    type = RankTwoAux
    rank_two_tensor = total_strain
    variable = strain_xy
    index_i = 0
    index_j = 1
  [../]
  [./strain_yy]
    type = RankTwoAux
    rank_two_tensor = total_strain
    variable = strain_yy
    index_i = 1
    index_j = 1
  [../]
  [./aux_strain_zz]
    type = RankTwoAux
    rank_two_tensor = total_strain
    variable = aux_strain_zz
    index_i = 2
    index_j = 2
  [../]
[]

[Functions]
  [./tempfunc]
    type = ParsedFunction
    expression = '(1-x)*t'
  [../]
[]

[BCs]
  [./bottom1x]
    type = DirichletBC
    boundary = 1
    variable = disp_x
    value = 0.0
  [../]
  [./bottom1y]
    type = DirichletBC
    boundary = 1
    variable = disp_y
    value = 0.0
  [../]
  [./bottom2x]
    type = DirichletBC
    boundary = 2
    variable = disp_x
    value = 0.0
  [../]
  [./bottom2y]
    type = DirichletBC
    boundary = 2
    variable = disp_y
    value = 0.0
  [../]
[]

[Materials]
  [./elastic_tensor]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = 0.3
    youngs_modulus = 1e6
    block = '1 2'
  [../]
  [./strain1]
    type = ComputePlaneSmallStrain
    displacements = 'disp_x disp_y'
    scalar_out_of_plane_strain = scalar_strain_zz1
    block = 1
    eigenstrain_names = eigenstrain
  [../]
  [./strain2]
    type = ComputePlaneSmallStrain
    displacements = 'disp_x disp_y'
    scalar_out_of_plane_strain = scalar_strain_zz2
    block = 2
    eigenstrain_names = eigenstrain
  [../]
  [./thermal_strain]
    type = ComputeThermalExpansionEigenstrain
    temperature = temp
    thermal_expansion_coeff = 0.02
    stress_free_temperature = 0.5
    block = '1 2'
    eigenstrain_name = eigenstrain
  [../]
  [./stress]
    type = ComputeLinearElasticStress
    block = '1 2'
  [../]
[]

[Executioner]
  type = Transient

  solve_type = PJFNK
  line_search = none

# controls for linear iterations
  l_max_its = 100
  l_tol = 1e-10

# controls for nonlinear iterations
  nl_max_its = 15
  nl_rel_tol = 1e-12
  nl_abs_tol = 1e-10

# time control
  start_time = 0.0
  dt = 1.0
  dtmin = 1.0
  end_time = 2.0
  num_steps = 5000
[]

[Outputs]
  exodus = true
[]

(modules/solid_mechanics/test/tests/2D_different_planes/gps_xz.i)

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  file = square_xz_plane.e
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_z]
  [../]
  [./scalar_strain_yy]
    order = FIRST
    family = SCALAR
  [../]
[]

[AuxVariables]
  [./temp]
  [../]
  [./disp_y]
  [../]
[]

[Physics/SolidMechanics/QuasiStatic]
  [./generalized_plane_strain]
    block = 1
    strain = SMALL
    scalar_out_of_plane_strain = scalar_strain_yy
    out_of_plane_direction = y
    planar_formulation = GENERALIZED_PLANE_STRAIN
    eigenstrain_names = 'eigenstrain'
    generate_output = 'stress_xx stress_xz stress_yy stress_zz strain_xx strain_xz strain_yy strain_zz'
  [../]
[]

[AuxKernels]
  [./tempfuncaux]
    type = FunctionAux
    variable = temp
    function = tempfunc
  [../]
[]

[Functions]
  [./tempfunc]
    type = ParsedFunction
    expression = '(1-x)*t'
  [../]
[]

[BCs]
  [./bottomx]
    type = DirichletBC
    boundary = 3
    variable = disp_x
    value = 0.0
  [../]
  [./bottomy]
    type = DirichletBC
    boundary = 3
    variable = disp_z
    value = 0.0
  [../]
[]

[Materials]
  [./elastic_stress]
    type = ComputeLinearElasticStress
    block = 1
  [../]
  [./thermal_strain]
    type = ComputeThermalExpansionEigenstrain
    temperature = temp
    thermal_expansion_coeff = 0.02
    stress_free_temperature = 0.5
    eigenstrain_name = eigenstrain
  [../]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = 1
    poissons_ratio = 0.3
    youngs_modulus = 1e6
  [../]
[]

[Postprocessors]
  [./react_y]
    type = MaterialTensorIntegral
    use_displaced_mesh = false
    rank_two_tensor = stress
    index_i = 1
    index_j = 1
  [../]
[]

[Executioner]
  type = Transient

  solve_type = PJFNK
  line_search = none

# controls for linear iterations
  l_max_its = 100
  l_tol = 1e-10

# controls for nonlinear iterations
  nl_max_its = 10
  nl_rel_tol = 1e-12

# time control
  start_time = 0.0
  dt = 1.0
  dtmin = 1.0
  end_time = 2.0
[]

[Outputs]
  file_base = gps_xz_small_out
  exodus = true
[]

(modules/solid_mechanics/test/tests/generalized_plane_strain/tests)

[Tests]
  issues = '#5042'
  parallel_scheduling = true
  [plane_strain]
    type = 'Exodiff'
    input = 'plane_strain.i'
    exodiff = 'plane_strain_out.e'
    custom_cmp = 'plane_strain.exodiff'
    requirement = 'The system shall support a traditional plane strain mechanics solution'
    design = 'ComputePlaneSmallStrain.md'
  []
  [plane_strain_prescribed]
    type = 'Exodiff'
    input = 'plane_strain_prescribed.i'
    exodiff = 'plane_strain_prescribed_out.e'
    custom_cmp = 'plane_strain.exodiff'
    requirement = 'The system shall support a traditional plane strain mechanics solution where the out-of-plane strain is prescribed'
    design = 'ComputePlaneSmallStrain.md'
  []
  [generalized_plane_strain_small]
    type = 'Exodiff'
    input = 'generalized_plane_strain_small.i'
    exodiff = 'generalized_plane_strain_small_out.e'
    custom_cmp = 'generalized.exodiff'
    requirement = 'The system shall support a generalized plane strain mechanics solution'
    design = 'ComputePlaneSmallStrain.md'
  []
  [generalized_plane_strain_small_ad]
    type = 'Exodiff'
    input = 'generalized_plane_strain_small.i'
    exodiff = 'generalized_plane_strain_small_out.e'
    custom_cmp = 'generalized.exodiff'
    prereq = generalized_plane_strain_small
    cli_args = 'Physics/SolidMechanics/QuasiStatic/all/use_automatic_differentiation=true '
               'Postprocessors/react_z/type=ADMaterialTensorIntegral '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    requirement = 'The system shall support an automatic-differentiation generalized plane strain mechanics solution'
    design = 'ADComputePlaneSmallStrain.md'
  []
  [generalized_plane_strain_small_ad_jacobian]
    type = 'PetscJacobianTester'
    input = 'generalized_plane_strain_small.i'
    run_sim = 'True'
    ratio_tol = 1e-7
    difference_tol = 1e-2
    cli_args = 'Physics/SolidMechanics/QuasiStatic/all/use_automatic_differentiation=true '
               'Postprocessors/react_z/type=ADMaterialTensorIntegral '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    requirement = 'The system shall compute the Jacobian for an automatic-differentiation generalized plane strain mechanics solution'
    design = 'ADComputePlaneSmallStrain.md'
  []
  [generalized_plane_strain_ref_resid_absolute_value_vector_tags]
    type = 'Exodiff'
    input = 'generalized_plane_strain_ref_resid.i'
    exodiff = 'generalized_plane_strain_ref_resid_out.e'
    custom_cmp = 'generalized.exodiff'
    requirement = 'The system shall support a generalized plane strain mechanics solution using the reference residual approach to check solution convergence of the field and scalar variables when the reference vector is constructed using absolute value vector tags.'
    design = 'ComputePlaneSmallStrain.md'
  []
  [generalized_plane_strain_ref_resid_absolute_value_vector_tags_ad]
    type = 'Exodiff'
    input = 'generalized_plane_strain_ref_resid.i'
    exodiff = 'generalized_plane_strain_ref_resid_out.e'
    custom_cmp = 'generalized.exodiff'
    prereq = generalized_plane_strain_ref_resid_absolute_value_vector_tags
    cli_args = 'Physics/SolidMechanics/QuasiStatic/all/use_automatic_differentiation=true '
               'Postprocessors/react_z/type=ADMaterialTensorIntegral '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    requirement = 'The system shall support an automatic-differentiation generalized plane strain mechanics solution using the reference residual approach to check solution convergence of the field and scalar variables when the reference vector is constructed using absolute value vector tags.'
    design = 'ADComputePlaneSmallStrain.md'
  []
  [generalized_plane_strain_ref_resid_extra_vector_tags]
    prereq = generalized_plane_strain_ref_resid_absolute_value_vector_tags
    type = 'Exodiff'
    input = 'generalized_plane_strain_ref_resid.i'
    cli_args = "Physics/SolidMechanics/QuasiStatic/all/absolute_value_vector_tags='' "
               "Physics/SolidMechanics/QuasiStatic/all/extra_vector_tags='ref' "
               "Outputs/file_base=generalized_plane_strain_ref_resid_extra_vector_tags_out"
    exodiff = 'generalized_plane_strain_ref_resid_extra_vector_tags_out.e'
    custom_cmp = 'generalized.exodiff'
    requirement = 'The system shall support a generalized plane strain mechanics solution using the reference residual approach to check solution convergence of the field and scalar variables when the reference vector is constructed using unbalanced kernel contributions.'
    design = 'ComputePlaneSmallStrain.md'
  []
  [generalized_plane_strain_ref_resid_extra_vector_tags_ad]
    prereq = 'generalized_plane_strain_ref_resid_extra_vector_tags generalized_plane_strain_ref_resid_absolute_value_vector_tags_ad'
    type = 'Exodiff'
    input = 'generalized_plane_strain_ref_resid.i'
    cli_args = "Physics/SolidMechanics/QuasiStatic/all/use_automatic_differentiation=true "
               "Physics/SolidMechanics/QuasiStatic/all/absolute_value_vector_tags='' "
               "Physics/SolidMechanics/QuasiStatic/all/extra_vector_tags='ref' "
               "Postprocessors/react_z/type=ADMaterialTensorIntegral "
               "Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor "
               "Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain "
               "Materials/stress/type=ADComputeLinearElasticStress "
               "Outputs/file_base=generalized_plane_strain_ref_resid_extra_vector_tags_out"
    exodiff = 'generalized_plane_strain_ref_resid_extra_vector_tags_out.e'
    custom_cmp = 'generalized.exodiff'
    requirement = 'The system shall support an automatic-differentiation generalized plane strain mechanics solution using the reference residual approach to check solution convergence of the field and scalar variables when the reference vector is constructed using unbalanced kernel contributions.'
    design = 'ADComputePlaneSmallStrain.md'
  []
  [generalized_plane_strain_increment]
    type = 'Exodiff'
    input = 'generalized_plane_strain_increment.i'
    exodiff = 'generalized_plane_strain_increment_out.e'
    custom_cmp = 'generalized.exodiff'
    requirement = 'The system shall support a generalized plane strain mechanics solution with incremental strain'
    design = 'ComputePlaneIncrementalStrain.md'
  []
  [generalized_plane_strain_increment_ad]
    type = 'Exodiff'
    input = 'generalized_plane_strain_increment.i'
    exodiff = 'generalized_plane_strain_increment_out.e'
    custom_cmp = 'generalized.exodiff'
    prereq = generalized_plane_strain_increment
    cli_args = 'Physics/SolidMechanics/QuasiStatic/all/use_automatic_differentiation=true '
               'Postprocessors/react_z/type=ADMaterialTensorIntegral '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain '
               'Materials/stress/type=ADComputeStrainIncrementBasedStress'
    requirement = 'The system shall support an automatic-differentiation generalized plane strain mechanics solution with incremental strain'
    design = 'ADComputePlaneIncrementalStrain.md'
  []
  [generalized_plane_strain_finite]
    type = 'Exodiff'
    input = 'generalized_plane_strain_finite.i'
    exodiff = 'generalized_plane_strain_finite_out.e'
    custom_cmp = 'generalized.exodiff'
    requirement = 'The system shall support a generalized plane strain mechanics solution with finite strain'
    design = 'ComputePlaneFiniteStrain.md'
  []
  [generalized_plane_strain_finite_ad]
    type = 'Exodiff'
    input = 'generalized_plane_strain_finite.i'
    exodiff = 'generalized_plane_strain_finite_out.e'
    custom_cmp = 'generalized.exodiff'
    prereq = generalized_plane_strain_finite
    cli_args = 'Physics/SolidMechanics/QuasiStatic/all/use_automatic_differentiation=true '
               'Postprocessors/react_z/type=ADMaterialTensorIntegral '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain '
               'Materials/stress/type=ADComputeFiniteStrainElasticStress'
    requirement = 'The system shall support an automatic-differentiation generalized plane strain mechanics solution with finite strain'
    design = 'ADComputePlaneFiniteStrain.md'
  []
  [generalized_plane_strain_squares]
    type = 'Exodiff'
    input = 'generalized_plane_strain_squares.i'
    exodiff = 'generalized_plane_strain_squares_out.e'
    abs_zero = 1e-7
    requirement = 'The system shall support setting up a generalized plane strain problem through an action'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#7840'
  []
  [generalized_plane_strain_squares_ad]
    type = 'Exodiff'
    input = 'generalized_plane_strain_squares.i'
    exodiff = 'generalized_plane_strain_squares_out.e'
    prereq = generalized_plane_strain_squares
    abs_zero = 1e-7
    cli_args = 'Physics/SolidMechanics/GeneralizedPlaneStrain/gps1/use_automatic_differentiation=true '
               'Physics/SolidMechanics/GeneralizedPlaneStrain/gps2/use_automatic_differentiation=true '
               'Kernels/SolidMechanics/use_automatic_differentiation=true '
               'Postprocessors/react_z1/type=ADMaterialTensorIntegral '
               'Postprocessors/react_z2/type=ADMaterialTensorIntegral '
               'AuxKernels/stress_xx/type=ADRankTwoAux AuxKernels/stress_xy/type=ADRankTwoAux '
               'AuxKernels/stress_yy/type=ADRankTwoAux AuxKernels/stress_zz/type=ADRankTwoAux '
               'AuxKernels/strain_xx/type=ADRankTwoAux AuxKernels/strain_xy/type=ADRankTwoAux '
               'AuxKernels/strain_yy/type=ADRankTwoAux AuxKernels/aux_strain_zz/type=ADRankTwoAux '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/strain1/type=ADComputePlaneSmallStrain '
               'Materials/strain2/type=ADComputePlaneSmallStrain '
               'Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    requirement = 'The system shall support setting up an automatic-differentiation generalized plane strain problem through an action'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#7840'
  []
  [generalized_plane_strain_squares_ad_jacobian]
    type = 'PetscJacobianTester'
    input = 'generalized_plane_strain_squares.i'
    run_sim = 'True'
    ratio_tol = 1e-7
    difference_tol = 1e0
    cli_args = 'Physics/SolidMechanics/GeneralizedPlaneStrain/gps1/use_automatic_differentiation=true '
               'Physics/SolidMechanics/GeneralizedPlaneStrain/gps2/use_automatic_differentiation=true '
               'Kernels/SolidMechanics/use_automatic_differentiation=true '
               'Postprocessors/react_z1/type=ADMaterialTensorIntegral '
               'Postprocessors/react_z2/type=ADMaterialTensorIntegral '
               'AuxKernels/stress_xx/type=ADRankTwoAux AuxKernels/stress_xy/type=ADRankTwoAux '
               'AuxKernels/stress_yy/type=ADRankTwoAux AuxKernels/stress_zz/type=ADRankTwoAux '
               'AuxKernels/strain_xx/type=ADRankTwoAux AuxKernels/strain_xy/type=ADRankTwoAux '
               'AuxKernels/strain_yy/type=ADRankTwoAux AuxKernels/aux_strain_zz/type=ADRankTwoAux '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/strain1/type=ADComputePlaneSmallStrain '
               'Materials/strain2/type=ADComputePlaneSmallStrain '
               'Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    requirement = 'The system shall compute the Jacobian for an automatic-differentiation generalized plane strain action with multiple action blocks'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#7840'
  []
  [out_of_plane_pressure_function]
    type = 'Exodiff'
    input = 'out_of_plane_pressure.i'
    exodiff = 'out_of_plane_pressure_out.e'
    abs_zero = 1e-9
    requirement = 'The system shall support setting the out-of-plane pressure using a function for generalized plane strain problems'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#7840'
  []
  [out_of_plane_pressure_function_ad]
    type = 'Exodiff'
    input = 'out_of_plane_pressure.i'
    exodiff = 'out_of_plane_pressure_out.e'
    prereq = out_of_plane_pressure_material
    abs_zero = 1e-9
    cli_args = 'Physics/SolidMechanics/GeneralizedPlaneStrain/gps/use_automatic_differentiation=true '
               'Kernels/SolidMechanics/use_automatic_differentiation=true '
               'Postprocessors/react_z/type=ADMaterialTensorIntegral '
               'AuxKernels/stress_xx/type=ADRankTwoAux AuxKernels/stress_xy/type=ADRankTwoAux '
               'AuxKernels/stress_yy/type=ADRankTwoAux AuxKernels/stress_zz/type=ADRankTwoAux '
               'AuxKernels/strain_xx/type=ADRankTwoAux AuxKernels/strain_xy/type=ADRankTwoAux '
               'AuxKernels/strain_yy/type=ADRankTwoAux AuxKernels/strain_zz/type=ADRankTwoAux '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/strain/type=ADComputePlaneSmallStrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    requirement = 'The system shall support setting the out-of-plane pressure using a function for automatic-differentiation generalized plane strain problems'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#7840'
  []
  [out_of_plane_pressure_material]
    type = 'Exodiff'
    input = 'out_of_plane_pressure.i'
    exodiff = 'out_of_plane_pressure_out.e'
    cli_args = 'Physics/SolidMechanics/GeneralizedPlaneStrain/gps/out_of_plane_pressure_function=0 Physics/SolidMechanics/GeneralizedPlaneStrain/gps/out_of_plane_pressure_material=traction_material'
    abs_zero = 1e-9
    requirement = 'The system shall support setting the out-of-plane pressure using a material property for generalized plane strain problems'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#17250'
    prereq = out_of_plane_pressure_function
  []
  [out_of_plane_pressure_material_ad]
    type = 'Exodiff'
    input = 'out_of_plane_pressure.i'
    exodiff = 'out_of_plane_pressure_out.e'
    cli_args = 'Physics/SolidMechanics/GeneralizedPlaneStrain/gps/use_automatic_differentiation=true '
               'Physics/SolidMechanics/GeneralizedPlaneStrain/gps/out_of_plane_pressure_function=0 '
               'Physics/SolidMechanics/GeneralizedPlaneStrain/gps/out_of_plane_pressure_material=traction_material '
               'Kernels/SolidMechanics/use_automatic_differentiation=true '
               'Postprocessors/react_z/type=ADMaterialTensorIntegral '
               'AuxKernels/stress_xx/type=ADRankTwoAux AuxKernels/stress_xy/type=ADRankTwoAux '
               'AuxKernels/stress_yy/type=ADRankTwoAux AuxKernels/stress_zz/type=ADRankTwoAux '
               'AuxKernels/strain_xx/type=ADRankTwoAux AuxKernels/strain_xy/type=ADRankTwoAux '
               'AuxKernels/strain_yy/type=ADRankTwoAux AuxKernels/strain_zz/type=ADRankTwoAux '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/strain/type=ADComputePlaneSmallStrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    abs_zero = 1e-9
    requirement = 'The system shall support setting the out-of-plane pressure using a material property for automatic-differentiation generalized plane strain problems'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#17250'
    prereq = out_of_plane_pressure_function_ad
  []
  [out_of_plane_pressure_material_ad_jacobian]
    type = 'PetscJacobianTester'
    input = 'out_of_plane_pressure.i'
    run_sim = 'True'
    ratio_tol = 1e-7
    difference_tol = 1e-1
    cli_args = 'Physics/SolidMechanics/GeneralizedPlaneStrain/gps/use_automatic_differentiation=true '
               'Physics/SolidMechanics/GeneralizedPlaneStrain/gps/out_of_plane_pressure_function=0 '
               'Physics/SolidMechanics/GeneralizedPlaneStrain/gps/out_of_plane_pressure_material=traction_material '
               'Kernels/SolidMechanics/use_automatic_differentiation=true '
               'Postprocessors/react_z/type=ADMaterialTensorIntegral '
               'AuxKernels/stress_xx/type=ADRankTwoAux AuxKernels/stress_xy/type=ADRankTwoAux '
               'AuxKernels/stress_yy/type=ADRankTwoAux AuxKernels/stress_zz/type=ADRankTwoAux '
               'AuxKernels/strain_xx/type=ADRankTwoAux AuxKernels/strain_xy/type=ADRankTwoAux '
               'AuxKernels/strain_yy/type=ADRankTwoAux AuxKernels/strain_zz/type=ADRankTwoAux '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/strain/type=ADComputePlaneSmallStrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    requirement = 'The system shall compute the Jacobian for automatic-differentiation generalized plane strain with material-property out-of-plane pressure'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#17250'
  []
  [generalized_plane_strain_scalar_vector]
    type = 'Exodiff'
    input = 'generalized_plane_strain_scalar_vector.i'
    exodiff = 'generalized_plane_strain_scalar_vector_out.e'
    prereq = 'generalized_plane_strain_squares'
    abs_zero = 1e-7
    allow_test_objects = true
    requirement = 'The system shall support listing all of the out-of-plane strain variables in the strain calculator'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#7840'
  []
  [generalized_plane_strain_scalar_vector_ad]
    type = 'Exodiff'
    input = 'generalized_plane_strain_scalar_vector.i'
    exodiff = 'generalized_plane_strain_scalar_vector_out.e'
    prereq = 'generalized_plane_strain_scalar_vector'
    abs_zero = 1e-7
    allow_test_objects = true
    cli_args = 'Physics/SolidMechanics/GeneralizedPlaneStrain/gps1/use_automatic_differentiation=true '
               'Physics/SolidMechanics/GeneralizedPlaneStrain/gps2/use_automatic_differentiation=true '
               'Kernels/SolidMechanics/use_automatic_differentiation=true '
               'Postprocessors/react_z1/type=ADMaterialTensorIntegral '
               'Postprocessors/react_z2/type=ADMaterialTensorIntegral '
               'AuxKernels/stress_xx/type=ADRankTwoAux AuxKernels/stress_xy/type=ADRankTwoAux '
               'AuxKernels/stress_yy/type=ADRankTwoAux AuxKernels/stress_zz/type=ADRankTwoAux '
               'AuxKernels/strain_xx/type=ADRankTwoAux AuxKernels/strain_xy/type=ADRankTwoAux '
               'AuxKernels/strain_yy/type=ADRankTwoAux AuxKernels/aux_strain_zz/type=ADRankTwoAux '
               'Materials/elastic_tensor/type=ADComputeIsotropicElasticityTensor '
               'Materials/strain1/type=ADComputePlaneSmallStrain '
               'Materials/thermal_strain/type=ADComputeThermalExpansionEigenstrain '
               'Materials/stress/type=ADComputeLinearElasticStress'
    requirement = 'The system shall support listing all of the out-of-plane strain variables in the automatic-differentiation strain calculator'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#7840'
  []
  [generalized_plane_strain_auto_scalar]
    type = 'Exodiff'
    input = 'out_of_plane_pressure.i'
    cli_args = "Variables/inactive=scalar_strain_zz
                Problem/nl_sys_names=mechanics
                Variables/disp_x/solver_sys=mechanics Variables/disp_y/solver_sys=mechanics"
    exodiff = 'out_of_plane_pressure_out.e'
    prereq = 'out_of_plane_pressure_material_ad'
    abs_zero = 1e-9
    requirement = 'The system shall support automatic creation of the scalar out-of-plane strain variable for legacy generalized plane strain problems'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#7840'
  []
  [generalized_plane_strain_auto_scalar_ad]
    type = 'Exodiff'
    input = 'generalized_plane_strain_auto_scalar.i'
    exodiff = 'out_of_plane_pressure_out.e'
    prereq = 'generalized_plane_strain_auto_scalar'
    abs_zero = 1e-9
    requirement = 'The system shall support automatic creation of the scalar out-of-plane strain variable for automatic-differentiation generalized plane strain problems'
    design = 'GeneralizedPlaneStrainAction.md'
    issues = '#7840'
  []
[]

(modules/solid_mechanics/test/tests/generalized_plane_strain/generalized_plane_strain_auto_scalar.i)

[Mesh]
  [square]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 2
    ny = 2
  []
  displacements = 'disp_x disp_y'
[]

[Problem]
  nl_sys_names = 'unused mechanics'
[]

[Variables]
  [unused]
    solver_sys = unused
    outputs = none
  []
  [disp_x]
    solver_sys = mechanics
  []
  [disp_y]
    solver_sys = mechanics
  []
[]

[AuxVariables]
  [saved_x]
    order = FIRST
    family = LAGRANGE
  []
  [saved_y]
    order = FIRST
    family = LAGRANGE
  []

  [stress_xx]
    order = CONSTANT
    family = MONOMIAL
  []
  [stress_xy]
    order = CONSTANT
    family = MONOMIAL
  []
  [stress_yy]
    order = CONSTANT
    family = MONOMIAL
  []
  [stress_zz]
    order = CONSTANT
    family = MONOMIAL
  []

  [strain_xx]
    order = CONSTANT
    family = MONOMIAL
  []
  [strain_xy]
    order = CONSTANT
    family = MONOMIAL
  []
  [strain_yy]
    order = CONSTANT
    family = MONOMIAL
  []
  [strain_zz]
    order = CONSTANT
    family = MONOMIAL
  []
[]

[Postprocessors]
  [react_z]
    type = ADMaterialTensorIntegral
    rank_two_tensor = stress
    index_i = 2
    index_j = 2
  []
[]

[Physics]
  [SolidMechanics]
    [GeneralizedPlaneStrain]
      [gps]
        use_automatic_differentiation = true
        use_displaced_mesh = true
        displacements = 'disp_x disp_y'
        scalar_out_of_plane_strain = scalar_strain_zz
        out_of_plane_pressure_function = traction_function
        pressure_factor = 1e5
      []
    []
  []
[]

[Kernels]
  [unused]
    type = ADReaction
    variable = unused
  []
  [SolidMechanics]
    use_automatic_differentiation = true
    use_displaced_mesh = false
    displacements = 'disp_x disp_y'
    save_in = 'saved_x saved_y'
  []
[]

[AuxKernels]
  [stress_xx]
    type = ADRankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_i = 0
    index_j = 0
  []
  [stress_xy]
    type = ADRankTwoAux
    rank_two_tensor = stress
    variable = stress_xy
    index_i = 0
    index_j = 1
  []
  [stress_yy]
    type = ADRankTwoAux
    rank_two_tensor = stress
    variable = stress_yy
    index_i = 1
    index_j = 1
  []
  [stress_zz]
    type = ADRankTwoAux
    rank_two_tensor = stress
    variable = stress_zz
    index_i = 2
    index_j = 2
  []

  [strain_xx]
    type = ADRankTwoAux
    rank_two_tensor = total_strain
    variable = strain_xx
    index_i = 0
    index_j = 0
  []
  [strain_xy]
    type = ADRankTwoAux
    rank_two_tensor = total_strain
    variable = strain_xy
    index_i = 0
    index_j = 1
  []
  [strain_yy]
    type = ADRankTwoAux
    rank_two_tensor = total_strain
    variable = strain_yy
    index_i = 1
    index_j = 1
  []
  [strain_zz]
    type = ADRankTwoAux
    rank_two_tensor = total_strain
    variable = strain_zz
    index_i = 2
    index_j = 2
  []
[]

[Functions]
  [traction_function]
    type = PiecewiseLinear
    x = '0  2'
    y = '0  1'
  []
[]

[BCs]
  [leftx]
    type = DirichletBC
    boundary = 3
    variable = disp_x
    value = 0.0
  []
  [bottomy]
    type = DirichletBC
    boundary = 0
    variable = disp_y
    value = 0.0
  []
[]

[Materials]
  [elastic_tensor]
    type = ADComputeIsotropicElasticityTensor
    poissons_ratio = 0.3
    youngs_modulus = 1e6
  []
  [strain]
    type = ADComputePlaneSmallStrain
    displacements = 'disp_x disp_y'
    scalar_out_of_plane_strain = scalar_strain_zz
  []
  [stress]
    type = ADComputeLinearElasticStress
  []
  [traction_material]
    type = ADGenericFunctionMaterial
    prop_names = traction_material
    prop_values = traction_function
  []
[]

[Executioner]
  type = Transient

  solve_type = PJFNK
  line_search = none

  l_max_its = 100
  l_tol = 1e-4

  nl_max_its = 15
  nl_rel_tol = 1e-14
  nl_abs_tol = 1e-11

  start_time = 0.0
  dt = 1.0
  dtmin = 1.0
  end_time = 2.0
  num_steps = 2
[]

[Outputs]
  exodus = true
  file_base = out_of_plane_pressure_out
[]

Description
Constructed MooseObjects
Example Input Syntax
Input Parameters
Associated Actions
Available Actions