Solid Mechanics Failure Analysis Report

[./1d_finite_rz_coord_axis_error]
  type = 'RunException'
  prereq = axisymmetric_gps_finite
  input = 'axisymm_gps_finite.i'
  cli_args = 'Problem/rz_coord_axis=X'
  expect_err = 'rz_coord_axis=Y is the only supported option for StressDivergenceRZTensors'
  design = 'source/materials/ComputeAxisymmetric1DFiniteStrain.md'
  requirement = 'The StressDivergenceRZTensors class shall generate an error if used with Problem/rz_coord_axis set to anything other than Y'
[../]

[two_block.error_unrestricted]
  type = RunException
  input = "two_block.i"
  expect_err = "No SolidMechanics/QuasiStatic action can be block unrestricted if more than one SolidMechanics/QuasiStatic action is specified."
  design = "syntax/Physics/SolidMechanics/QuasiStatic/index.md"
  requirement = "The QuasiStatic SolidMechanics Physics shall error if an input file does not specify block restrictions for the MasterAction in input files with more than one instance of the MasterAction block."
[]

[two_block.error_overlapping]
  type = RunException
  input = "two_block.i"
  cli_args = "Physics/SolidMechanics/QuasiStatic/block1/block=2"
  expect_err = "The block restrictions in the SolidMechanics/QuasiStatic actions must be non-overlapping."
  design = "syntax/Physics/SolidMechanics/QuasiStatic/index.md"
  requirement = "The QuasiStatic SolidMechanics Physics shall error if an input file specifies overlapping block restrictions for the MasterAction in input files with more than one instance of the MasterAction block."
[]

[no_block]
  type = RunException
  input = "no_block.i"
  cli_args = "--error"
  expect_err = "Common parameters are supplied, but not used in"
  design = "syntax/Physics/SolidMechanics/QuasiStatic/index.md"
  requirement = "The QuasiStatic SolidMechanics Physics shall warn if global Master action parameters are supplied but no Master action subblock have been added."
[]

[two_coord.error_different_coords]
  type = RunException
  input = "two_coord.i"
  cli_args = "Physics/SolidMechanics/QuasiStatic/active=error"
  expect_err = "The SolidMechanics action requires all subdomains to have the same coordinate system."
  design = "syntax/Physics/SolidMechanics/QuasiStatic/index.md"
  requirement = "The QuasiStatic SolidMechanics Physics shall error if an input file assigns the same QuasiStatic SolidMechanics Physics block to mesh blocks with different coordinate systems."
[]

[material_output_order.error1]
  type = RunException
  input = "material_output_order.i"
  cli_args = 'Physics/SolidMechanics/QuasiStatic/material_output_family="LAGRANGE MONOMIAL"'
  expect_err = "The number of families assigned to material outputs must be: 0 to be assigned MONOMIAL; 1 to assign all outputs the same value, or the same size as the number of generate outputs listed."
  design = "syntax/Physics/SolidMechanics/QuasiStatic/index.md"
  requirement = "The QuasiStatic SolidMechanics Physics shall error if an the material outputs and or families do not match the number of material outputs."
[]

[./error_unrestricted]
  type = RunException
  input = 'two_block.i'
  expect_err = 'No SolidMechanics/QuasiStatic action can be block unrestricted if more than one SolidMechanics/QuasiStatic action is specified.'
  prereq = 'two_block'
  requirement = 'The QuasiStatic SolidMechanics Physics shall error if an input file does not specify block restrictions for the MasterAction in input files with more than one instance of the MasterAction block.'
[../]

[./error_overlapping]
  type = RunException
  input = 'two_block.i'
  cli_args = 'Physics/SolidMechanics/QuasiStatic/block1/block=2'
  expect_err = 'The block restrictions in the SolidMechanics/QuasiStatic actions must be non-overlapping.'
  prereq = 'error_unrestricted'
  requirement = 'The QuasiStatic SolidMechanics Physics shall error if an input file specifies overlapping block restrictions for the MasterAction in input files with more than one instance of the MasterAction block.'
[../]

[./error_coord]
  type = RunException
  input = 'two_coord.i'
  cli_args = 'Physics/SolidMechanics/QuasiStatic/active=error'
  expect_err = 'The SolidMechanics action requires all subdomains to have the same coordinate system.'
  prereq = 'two_coord'
  requirement = 'The QuasiStatic SolidMechanics Physics shall error if an input file assigns the same QuasiStatic SolidMechanics Physics block to mesh blocks with different coordinate systems.'
[../]

[bending_exception]
  type = 'RunException'
  input = bending_jacobian.i
  cli_args = 'BCs/move_y/function="-100*t"'
  expect_err = 'Cannot take square root of a number less than or equal to zero in the calculation of C3_test for the Rashid approximation for the rotation tensor. This zero or negative number may occur when elements become heavily distorted.'
  issues = '#19067'
  requirement = 'Finite strain methods in Tensor Mechanics, using the auto differentiation capabilities, shall cut the timestep through  a mooseException when the loading conditions deform the elements so much as to produce a negative number under the square root term in the Rashid approximation for the rotation tensor calcuation'
[]

[cracking_err_diagonal]
  prereq = cracking
  type = RunException
  input = cracking.i
  expect_err = "Defaulting to the legacy option of 'DIAGONAL', but the 'FULL' option is preferred"
  requirement = "The ComputeSmearedCrackingStress model shall generate a warning if 'cracked_elasticity_type' is not specified"
[]

[cracking_err_softening]
  prereq = cracking_err_diagonal
  allow_warnings = true #Ignore warning about missing cracked_elasticity_type
  type = RunException
  input = cracking.i
  cli_args = "Materials/elastic_stress/softening_models='abrupt_softening abrupt_softening'"
  expect_err = "Either 1 or 3 softening models must be specified"
  requirement = "The ComputeSmearedCrackingStress model shall generate an error if an invalid number of softening models is specified"
[]

[zero]
  type = CSVDiff
  input = 'negative_porosity.i'
  csvdiff = 'negative_porosity_zero_out.csv'
  cli_args = 'Outputs/file_base=negative_porosity_zero_out Materials/lps/negative_behavior=ZERO'
  detail = 'by setting the porosity to zero.'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[initial]
  type = CSVDiff
  input = 'negative_porosity.i'
  csvdiff = 'negative_porosity_initial_out.csv'
  cli_args = 'Outputs/file_base=negative_porosity_initial_out Materials/lps/negative_behavior=INITIAL_CONDITION'
  detail = 'by setting the porosity to the initial condition.'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[exception]
  type = RunException
  input = 'negative_porosity.i'
  expect_err = 'porosity is negative'
  cli_args = 'Materials/lps/negative_behavior=EXCEPTION'
  detail = 'by throwing an exception.'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[./beam_action_test1]
  type = 'RunException'
  input = 'beam_action_chk.i'
  expect_err = 'LineElementAction: A vector of displacement variable names should be provided as'

  requirement = "The LineElementAction class shall produce an error when the displacement variables are not provided by the user."
[../]

[./beam_action_test2]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x"
    Physics/SolidMechanics/LineElement/QuasiStatic/block_1/use_displaced_mesh=true Physics/SolidMechanics/LineElement/QuasiStatic/block_1/strain_type=small Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotation_type=small'
  expect_err = 'LineElementAction: Wrong combination of `use_displaced_mesh`, `strain_type` and `rotation_type`.'

  requirement = "The LineElementAction class shall produce an error if the user provided inputs for `strain_type`, `rotation_type` and `use_displaced_mesh` parameters are not compatible."
[../]

[./beam_action_test3]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/save_in="res_x res_y"'
  expect_err = 'LineElementAction: Number of save_in variables should equal to the number of displacement variables'

  requirement = "The LineElementAction class shall produce an error if the number of variables listed in the `save_in` parameter differs from the number of displacement variables."
[../]

[./beam_action_test4]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/diag_save_in="res_x res_y"'
  expect_err = 'LineElementAction: Number of diag_save_in variables should equal to the number of displacement variables'

  requirement = "The LineElementAction class shall produce an error if the number of variables listed in the `diag_save_in` parameter differs from the number of displacement variables."
[../]

[./beam_action_test5]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements=""'
  expect_err = 'LineElementAction: Rotational variable names should be provided for beam elements using `rotations` parameter'

  requirement = "The LineElementAction class shall produce an error if the names for the rotational degrees of freedom are not provided by the user."
[../]

[./beam_action_test6]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations=""'
  expect_err = 'LineElementAction: Number of rotational and displacement variable names provided as input for beam should be same.'

  requirement = "The LineElementAction class shall produce an error if the number of rotational variables provided as input differs from the number of displacement variables."
[../]

[./beam_action_test7]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x"'
  expect_err = 'LineElementAction: `y_orientation`, `area`, `Iy` and `Iz` should be provided for'

  requirement = "The LineElementAction class shall produce an error if the moment of inertia, area and orientation of the beam are not provided as input."
[../]

[./beam_action_test8]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/area=1.0'
  expect_err = 'LineElementAction: Variable names for translational and rotational velocities and accelerations should be provided as input to perform dynamic simulation using beam elements using `velocities`, `accelerations`, `rotational_velocities`'

  requirement = "The LineElementAction class shall produce an error if translational and rotational velocities and accelerations are not provided as input for dynamic simulations using beam elements."
[../]

[./beam_action_test9]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/area=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotational_velocities="rot_vel_x rot_vel_y"'
  expect_err = 'LineElementAction: Number of translational and rotational velocity and acceleration variable names provided as input for the beam should be same as number of displacement variables'

  requirement = "The LineElementAction class shall produce an error if the number of translational and rotational velocities and accelerations differs from the number of displacement variables."
[../]

[./beam_action_test10]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/area=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotational_velocities="rot_vel_x"'
  expect_err = 'LineElementAction: Newmark time integration parameters `beta` and `gamma` should be provided as input to perform dynamic simulations using beams.'

  requirement = "The LineElementAction class shall produce an error if Newmark time integration parameters (`beta` and `gamma`) are not provided as input for dynamic simulations using beam elements."
[../]

[./beam_action_test11]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/area=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotational_velocities="rot_vel_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/beta=0.25'
  expect_err = 'LineElementAction: Either name of the density material property or a constant density value should be provided as input using `density` for creating the consistent mass/inertia matrix required for dynamic beam simulation.'

  requirement = "The LineElementAction class shall produce an error if density is not provided as input for dynamic beam simulations using beams elements with consistent mass/inertia matrix."
[../]

[./beam_action_test12]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/area=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotational_velocities="rot_vel_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/beta=0.25 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/density=1.0'
  expect_err = 'LineElementAction: `nodal_mass` or `nodal_mass_file` should be provided as input to calculate inertial forces on beam due to nodal mass.'

  requirement = "The LineElementAction class shall produce an error if nodal mass is not provided as input for dynamic beam simulations using beam elements with nodal mass matrix."
[../]

[./beam_action_test13]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/area=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotational_velocities="rot_vel_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/beta=0.25 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/density=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/nodal_mass=0.1'
  expect_err = 'LineElementAction: `nodal_Ixx`, `nodal_Iyy`, `nodal_Izz` should be provided as input to calculate inertial torque on beam due to nodal inertia.'

  requirement = "The LineElementAction class shall produce an error if nodal inertia is not provided as input for dynamic beam simulations using beam elements with nodal inertia matrix."
[../]

[./beam_action_test14]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/area=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotational_velocities="rot_vel_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/beta=0.25 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/density=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/nodal_mass=0.1 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/nodal_Ixx=0.1 Physics/SolidMechanics/LineElement/QuasiStatic/block_all/block="1 2"'
  expect_err = 'LineElementAction: The block restrictions in the LineElement actions must be non-overlapping.'

  requirement = "The LineElementAction class shall produce an error if multiple subblocks specify properties for the same mesh block."
[../]

[./beam_action_test15]
  type = 'RunException'
  input = 'beam_action_chk.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/block_1/displacements="disp_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotations="rot_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/area=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/rotational_velocities="rot_vel_x" Physics/SolidMechanics/LineElement/QuasiStatic/block_1/beta=0.25 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/density=1.0 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/nodal_mass=0.1 Physics/SolidMechanics/LineElement/QuasiStatic/block_1/nodal_Ixx=0.1'
  expect_err = 'LineElementAction: No LineElement action can be block unrestricted if more than one LineElement action is specified.'

  requirement = "The LineElementAction class shall produce an error if an action subblock is mesh block restricted while another is not."
[../]

[./beam_action_test16]
  type = 'RunException'
  input = '2_block_common.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/dynamic_nodal_translational_inertia=true'
  expect_err = 'LineElementAction: Variable names for translational and rotational velocities and accelerations should be provided as input to perform dynamic simulation using beam elements'

  requirement = "The LineElementAction class shall produce an error if `dynamic_nodal_translational_inertia` is set to true in the common action block but the subblocks do not have the parameters required for a dynamic beam simulation using beam elements."
[../]

[./error_1]
  type = 'RunException'
  input = 'dyn_euler_small.i'
  cli_args = "Kernels/inertial_force_x/rotations='rot_x rot_y' Outputs/csv=false Outputs/exodus=false"
  expect_err = "InertialForceBeam: The number of variables supplied in 'displacements' and 'rotations' must match."

  requirement = "The InertialForceBeam class shall produce an error if the number of variables provided for rotations differs from that provided for displacements."
[../]

[./error_2]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_inertia_damping.i'
  cli_args = "NodalKernels/rot_x_inertial/rotational_velocities='rot_vel_x rot_vel_y' Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalRotationalInertia: rotational_velocities and rotational_accelerations should be same size as rotations."
  superlu = true
  requirement = "The NodalRotatioanlInertia class shall produce an error if the number of rotational velocities and accelerations provided as input differ from the number of rotations."
[../]

[./error_3]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_inertia_damping.i'
  cli_args = "NodalKernels/rot_x_inertial/Ixx=1 NodalKernels/rot_x_inertial/Iyy=1 NodalKernels/rot_x_inertial/Izz=1 NodalKernels/rot_x_inertial/Ixy=1 Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalRotationalInertia: The moment of inertia tensor should be positive definite."
  superlu = true
  requirement = "The NodalRotationalInertia class shall produce an error if the user provided nodal inertia is not positive definite."
[../]

[./error_4]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_inertia_damping.i'
  cli_args = "NodalKernels/rot_x_inertial/x_orientation='2 0 0' NodalKernels/rot_x_inertial/y_orientation='0 1 0' Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalRotationalInertia: x_orientation and y_orientation must be unit vectors."
  superlu = true
  requirement = "The NodalRotatioanlInertia class shall produce an error if the user provided x and y orientations are not unit vectors."
[../]

[./error_5]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_inertia_damping.i'
  cli_args = "NodalKernels/rot_x_inertial/x_orientation='1 0 0' NodalKernels/rot_x_inertial/y_orientation='1 0 0' Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalRotationalInertia: x_orientation and y_orientation should be perpendicular to each other."
  superlu = true
  requirement = "The NodalRotatioanlInertia class shall produce an error if the user provided x and y orientations are not perpendicular to each other."
[../]

[./error_6]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_inertia_damping.i'
  cli_args = "NodalKernels/rot_x_inertial/x_orientation='1 0 0' Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalRotationalInertia: Both x_orientation and y_orientation should be provided if x_orientation or y_orientation is different from global x or y direction, respectively."
  superlu = trues
  requirement = "The NodalRotatioanlInertia class shall produce an error if only x or y orientation is provided as input by the user."
[../]

[./error_7]
  type = 'RunException'
  input = 'dyn_euler_small.i'
  cli_args = "Kernels/inertial_force_x/accelerations='accel_x accel_y' Outputs/csv=false Outputs/exodus=false"
  expect_err = "InertialForceBeam: The number of variables supplied in 'velocities', 'accelerations', 'rotational_velocities' and 'rotational_accelerations' must match the number of displacement variables."

  requirement = "The InertialForceBeam class shall produce an error if the number of translational and rotational velocities and accelerations provided as input differ from the number of displacement variables."
[../]

[./error_8]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_file.i'
  cli_args = "NodalKernels/x_inertial/nodal_mass_file=nodal_mass.csv NodalKernels/x_inertial/mass=2.0 Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalTranslationalInertia: Please provide either mass or nodal_mass_file as input, not both."
  superlu = true
  requirement = "The NodalTranslationalInertia class shall produce an error if nodal mass is provided as input both as a constant value and also using a csv file."
[../]

[./error_9]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_file.i'
  cli_args = "Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalTranslationalInertia: Please provide either mass or nodal_mass_file as input."
  superlu = true
  requirement = "The NodalTranslationalInertia class shall produce an error if nodal mass is not provided as input either as a constant value or using a csv file."
[../]

[./error_10]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_file.i'
  cli_args = "NodalKernels/x_inertial/nodal_mass_file=nodal_mass_col_test.csv Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalTranslationalInertia: The number of columns in"
  superlu = true
  requirement = "The NodalTranslationalInertia class shall produce an error if the number of columns in the nodal mass file is not 4."
[../]

[./error_11]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_file.i'
  cli_args = "NodalKernels/x_inertial/nodal_mass_file=nodal_mass.csv NodalKernels/x_inertial/boundary='right' Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalTranslationalInertia: Out of 2 nodal positions in"
  superlu = true
  requirement = "The NodalTranslationalInertia class shall produce an error if all the nodal positions provided in the nodal mass file cannot be found in the given boundary or mesh block."
[../]

[./error_12]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_gravity.i'
  cli_args = "NodalKernels/force_y2/nodal_mass_file=nodal_mass.csv NodalKernels/force_y2/mass=2.0 Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalGravity: Please provide either mass or nodal_mass_file as input, not both."
  superlu = true
  requirement = "The NodalGravity class shall produce an error if nodal mass is provided as input both as a constant value and also using a csv file."
[../]

[./error_13]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_gravity.i'
  cli_args = "Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalGravity: Please provide either mass or nodal_mass_file as input."
  superlu = true
  requirement = "The NodalGravity class shall produce an error if nodal mass is not provided as input either as a constant value or using a csv file."
[../]

[./error_14]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_gravity.i'
  cli_args = "NodalKernels/force_y2/nodal_mass_file=nodal_mass_col_test.csv Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalGravity: The number of columns in"
  superlu = true
  requirement = "The NodalGravity class shall produce an error if the number of columns in the nodal mass file is not 4."
[../]

[./error_15]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_gravity.i'
  cli_args = "NodalKernels/force_y2/nodal_mass_file=nodal_mass.csv NodalKernels/force_y2/boundary='right' Outputs/csv=false Outputs/exodus=false"
  expect_err = "NodalGravity: Out of 2 nodal positions in"
  superlu = true
  requirement = "The NodalGravity class shall produce an error if all the nodal positions provided in the nodal mass file cannot be found in the given boundary or mesh block."
[../]

[./error_16]
  type = 'RunException'
  input = 'dyn_euler_small_rayleigh_hht_action.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/all/add_dynamic_variables=false Outputs/csv=false Outputs/exodus=false'
  expect_err = "LineElementAction: When using 'dynamic_consistent_inertia', 'dynamic_nodal_rotational_inertia' or '_dynamic_nodal_translational_inertia', the velocity and acceleration AuxVariables and the corresponding AuxKernels are automatically set by"

  requirement = "The LineElementAction class shall produce an error if `add_dynamic_variables` option is set to false while `dynamic_consistent_inertia`, `dynamic_nodal_rotational_inertia` or `dynamic_nodal_translational_inertia` options are set to true."
  design = "C0TimoshenkoBeam.md LineElement/QuasiStatic/index.md"
[../]

[./error_17]
  type = 'RunException'
  input = 'dyn_euler_small_added_mass_inertia_damping_action.i'
  cli_args = 'Physics/SolidMechanics/LineElement/QuasiStatic/all/nodal_mass_file=nodal_mass.csv Outputs/csv=false Outputs/exodus=false'
  expect_err = "NodalTranslationalInertia: Please provide either mass or nodal_mass_file as input, not both"
  superlu = true
  requirement = "The NodalTranslationalInertia class shall produce an error if nodal mass is provided as input both as constant value and also using a csv file."
[../]

[./eigenstrain_from_var_test1]
  type = 'RunException'
  input = 'eigenstrain_from_var.i'
  cli_args = "Materials/thermal/displacement_eigenstrain_variables='thermal_eig zero1' Outputs/csv=false Outputs/exodus=false"
  expect_err = "ComputeEigenstrainBeamFromVariable: If the displacement or rotational eigenstrains are provided, it should contain 3 variables corresponding to the three components in the global coordinate system."

  requirement = "The ComputeEigenstrainBeamFromVariable class shall report an error if less than 3 displacement or rotational eigenstrains are provided by the user."
  design = "ComputeEigenstrainBeamFromVariable.md"
  issues = "#10313"
[../]

[./error_displacements1]
  type = 'RunException'
  input = 'euler_small_strain_y.i'
  cli_args = 'Materials/strain/rotations="rot_x" Outputs/exodus=false Outputs/csv=false'
  expect_err = "ComputeIncrementalBeamStrain: The number of variables supplied in 'displacements' and 'rotations' must match."

  requirement = "The ComputeIncrementalBeamStrain class shall produce an error if the number of supplied displacements and rotations do not match."
[../]

[./error_displacements2]
  type = 'RunException'
  input = 'euler_small_strain_y.i'
  cli_args = 'Kernels/solid_disp_x/rotations="rot_x" Outputs/exodus=false Outputs/csv=false'
  expect_err = "StressDivergenceBeam: The number of displacement and rotation variables should be same."

  requirement = "The StressDivergenceBeam class shall produce an error if the number of supplied displacements and rotations do not match."
[../]

[./error_large_strain]
  type = 'RunException'
  input = 'euler_small_strain_y.i'
  cli_args = 'Materials/strain/large_strain=true Materials/strain/Ay=1.0 Outputs/exodus=false Outputs/csv=false'
  expect_err = 'ComputeIncrementalBeamStrain: Large strain calculation does not currently support asymmetric beam configurations with non-zero first or third moments of area.'

  requirement = "The ComputeIncrementalBeamStrain class shall produce an error if large strain calculation is requested for asymmetric beam configurations with non-zero first or third moments of area."
[../]

[./error_y_orientation]
  type = 'RunException'
  input = 'euler_small_strain_y.i'
  cli_args = 'Materials/strain/y_orientation="0.5 0.5 0.0" Outputs/exodus=false Outputs/csv=false'
  expect_err = 'ComputeIncrementalBeamStrain: y_orientation should be perpendicular to the axis of the beam.'

  requirement = "The ComputeIncrementalBeamStrain class shall produce an error if the y orientation provided is not perpendicular to the beam axis."
[../]

[./except1]
  type = 'RunException'
  input = 'except1.i'
  expect_err = 'CappedWeakPlaneStressUpdate: Weak-plane friction and dilation angles must lie in \[0, Pi/2\]'
  requirement = "The CappedWeakPlaneStressUpdate model shall generate an error if the friction angle is negative"
  design = "CappedWeakPlaneStressUpdate.md"
  issues = "#7784"
[../]

[./except2]
  type = 'RunException'
  input = 'except2.i'
  expect_err = 'CappedWeakPlaneStressUpdate: Weak-plane friction and dilation angles must lie in \[0, Pi/2\]'
  requirement = "The CappedWeakPlaneStressUpdate model shall generate an error if the dilation angle is negative"
  design = "CappedWeakPlaneStressUpdate.md"
  issues = "#7784"
[../]

[./except3]
  type = 'RunException'
  input = 'except3.i'
  expect_err = 'CappedWeakPlaneStressUpdate: Weak-plane friction angle must not be less than dilation angle'
  requirement = "The CappedWeakPlaneStressUpdate model shall generate an error if the friction angle is less than the dilation angle"
  design = "CappedWeakPlaneStressUpdate.md"
  issues = "#7784"
[../]

[./except4]
  type = 'RunException'
  input = 'except4.i'
  expect_err = 'CappedWeakPlaneStressUpdate: Weak-plane cohesion must not be negative'
  requirement = "The CappedWeakPlaneStressUpdate model shall generate an error if the cohesion is negative"
  design = "CappedWeakPlaneStressUpdate.md"
  issues = "#7784"
[../]

[./except5]
  type = 'RunException'
  input = 'except5.i'
  expect_err = 'CappedWeakPlaneStressUpdate: Weak plane tensile strength plus compressive strength must be greater than smoothing_tol'
  requirement = "The CappedWeakPlaneStressUpdate model shall generate an error if the sum of the tensile and compressive strength is less than smoothing_tol"
  design = "CappedWeakPlaneStressUpdate.md"
  issues = "#7784"
[../]

[./except6]
  type = 'RunException'
  input = 'except6.i'
  expect_err = 'CappedWeakInclinedPlaneStressUpdate: normal_vector must not have zero length'
  requirement = "The CappedWeakPlaneStressUpdate model shall generate an error if the normal vector has zero length"
  design = "CappedWeakPlaneStressUpdate.md"
  issues = "#7784"
[../]

[./throw_test]
  type = 'RunException'
  input = 'throw_test.i'
  expect_err = "Solve failed and timestep already at or below dtmin"
  recover = false
  requirement = "The system shall permit exceptions to be thrown from material models with stateful properties without reading/writing to/from uninitialized memory"
  design = "CappedWeakPlaneStressUpdate.md"
  issues = "#7960"
[../]

[./num_constants]
  type = 'RunException'
  input = 'num_constants.i'
  expect_err = "Exactly two elastic constants must be defined for material 'elasticity_tensor'"
  design = "ComputeIsotropicElasticityTensor.md"
  issues = "#9438"
  requirement = "The system shall generate an error if a number of elastic constants other than two is supplied for an isotropic elasticity tensor"
[../]

[./youngs_modulus]
  type = 'RunException'
  input = 'youngs_modulus.i'
  expect_err = "Youngs modulus must be positive in material 'elasticity_tensor'"
  design = "ComputeIsotropicElasticityTensor.md"
  issues = "#9438"
  requirement = "The system shall generate an error if a non-positive Youngs modulus is supplied for an isotropic elasticity tensor"
[../]

[./bulk_modulus]
  type = 'RunException'
  input = 'bulk_modulus.i'
  expect_err = "Bulk modulus must be positive in material 'elasticity_tensor'"
  design = "ComputeIsotropicElasticityTensor.md"
  issues = "#9438"
  requirement = "The system shall generate an error if a non-positive bulk modulus is supplied for an isotropic elasticity tensor"
[../]

[./poissons_ratio]
  type = 'RunException'
  input = 'poissons_ratio.i'
  expect_err = "Poissons ratio must be greater than -1 and less than 0.5 in material 'elasticity_tensor'"
  design = "ComputeIsotropicElasticityTensor.md"
  issues = "#9438"
  requirement = "The system shall generate an error if a Poissons ratio outside the range from -1 to 0.5 is supplied for an isotropic elasticity tensor"
[../]

[./shear_modulus]
  type = 'RunException'
  input = 'shear_modulus.i'
  expect_err = "Shear modulus must not be negative in material 'elasticity_tensor'"
  design = "ComputeIsotropicElasticityTensor.md"
  issues = "#9438"
  requirement = "The system shall generate an error if a non-positive shear modulus is supplied for an isotropic elasticity tensor"
[../]

[jacobian_small_strain]
  type = CSVDiff
  input = 'czm_patch_test_base.i'
  csvdiff = 'czm_jacobian_small_strain_out.csv'
  cli_args = "Materials/czm_mat/type=PureElasticTractionSeparation Materials/czm_mat/normal_stiffness=100e2 Materials/czm_mat/tangent_stiffness=70e2   Physics/SolidMechanics/CohesiveZone/czm_ik/strain=SMALL Physics/SolidMechanics/CohesiveZone/czm_ik/generate_output='traction_x traction_y traction_z normal_traction tangent_traction jump_x jump_y jump_z normal_jump tangent_jump' Physics/SolidMechanics/QuasiStatic/all/strain=SMALL  Physics/SolidMechanics/QuasiStatic/all/use_finite_deform_jacobian=false Materials/stress/type=ADComputeLinearElasticStress Outputs/file_base=czm_jacobian_small_strain_out Outputs/exodus=true"
  mesh_mode = 'REPLICATED'
  recover = false
  detail = "The small strain cohesive zone model shall converge quadratically when using the hand-coded Jacobian"
[]

[jacobian_total_lagrangian]
  type = CSVDiff
  input = 'czm_patch_test_base.i'
  csvdiff = 'czm_jacobian_total_lagrangian_out.csv'
  cli_args = "Materials/czm_mat/type=PureElasticTractionSeparation Materials/czm_mat/normal_stiffness=100e2 Materials/czm_mat/tangent_stiffness=70e2 Physics/SolidMechanics/CohesiveZone/czm_ik/strain=FINITE Physics/SolidMechanics/CohesiveZone/czm_ik/generate_output='traction_x traction_y traction_z normal_traction tangent_traction jump_x jump_y jump_z normal_jump tangent_jump' Outputs/file_base=czm_jacobian_total_lagrangian_out Outputs/exodus=true"
  mesh_mode = 'REPLICATED'
  recover = false
  detail = "The finite strain cohesive zone model shall converge quadratically when using the hand-coded Jacobian"
[]

[small_strain_incremental]
  type = Exodiff
  input = 'czm_patch_test_base.i'
  exodiff = czm_small_strain_incremental_out.e
  cli_args = "Materials/czm_mat/type=PureElasticTractionSeparationIncremental Materials/czm_mat/normal_stiffness=100e2 Materials/czm_mat/tangent_stiffness=70e2   Physics/SolidMechanics/CohesiveZone/czm_ik/strain=SMALL Physics/SolidMechanics/CohesiveZone/czm_ik/generate_output='traction_x traction_y traction_z normal_traction tangent_traction jump_x jump_y jump_z normal_jump tangent_jump' Physics/SolidMechanics/QuasiStatic/all/strain=SMALL  Physics/SolidMechanics/QuasiStatic/all/use_finite_deform_jacobian=false Materials/stress/type=ADComputeLinearElasticStress Outputs/file_base=czm_small_strain_incremental_out Outputs/exodus=true"
  mesh_mode = 'REPLICATED'
  recover = false
  detail = "The small strain cohesive zone model shall return the same results of the total formulation when using an incremental material"
  allow_test_objects = true
[]

[total_lagrangian_incremental]
  type = Exodiff
  input = 'czm_patch_test_base.i'
  exodiff = czm_total_lagrangian_incremental_out.e
  cli_args = "Materials/czm_mat/type=PureElasticTractionSeparationIncremental Materials/czm_mat/normal_stiffness=100e2 Materials/czm_mat/tangent_stiffness=70e2 Physics/SolidMechanics/CohesiveZone/czm_ik/strain=FINITE Physics/SolidMechanics/CohesiveZone/czm_ik/generate_output='traction_x traction_y traction_z normal_traction tangent_traction jump_x jump_y jump_z normal_jump tangent_jump' Outputs/file_base=czm_total_lagrangian_incremental_out Outputs/exodus=true"
  mesh_mode = 'REPLICATED'
  recover = false
  detail = "The finite strain cohesive zone model shall return the same results of the total formulation when using an incremental material"
  allow_test_objects = true
[]

[multiple_action_and_materials]
  type = Exodiff
  input = 'czm_multiple_action_and_materials.i'
  exodiff = 'czm_multiple_action_and_materials_out.e'
  detail = "The system shall allow using multiple CohesiveZoneAction and different cohesive materials for different boundaries"
  cli_args = "Physics/SolidMechanics/CohesiveZone/strain=FINITE"
  mesh_mode = 'REPLICATED'
  recover = false
  allow_test_objects = true
[]

[multiple_action_different_strain]
  type = RunException
  input = 'czm_multiple_action_and_materials.i'
  expect_err = 'All blocks of the CohesiveZoneAction should have the same strain formulation'
  detail = "The system shall report an error when different strains are provided in different blocks of the CohesiveZoneAction"
  cli_args = "Physics/SolidMechanics/CohesiveZone/czm_ik_012/strain=FINITE Physics/SolidMechanics/CohesiveZone/czm_ik_23/strain=SMALL"
  mesh_mode = 'REPLICATED'
  recover = false
  allow_test_objects = true
[]

[czm_multiple_action_overlapping_boundary]
  type = RunException
  input = 'czm_multiple_action_and_materials.i'
  expect_err = 'The boundary restrictions in the CohesiveZoneAction actions must be non-overlapping.'
  detail = "The system shall report an error when overlapping boundaries are provided in different blocks of the CohesiveZoneAction"
  cli_args = "Physics/SolidMechanics/CohesiveZone/czm_ik_012/boundary='Block0_Block1' Physics/SolidMechanics/CohesiveZone/czm_ik_23/boundary='Block0_Block1'"
  mesh_mode = 'REPLICATED'
  recover = false
  allow_test_objects = true
[]

[multiple_action_common_outputs]
  type = Exodiff
  input = 'czm_multiple_action_and_materials.i'
  exodiff = czm_multiple_action_common_output_out.e
  detail = "The system shall use the common block generate output and additional outputs for all CZM blocks"
  cli_args = "Physics/SolidMechanics/CohesiveZone/generate_output='traction_x traction_y traction_z jump_x jump_y jump_z' Physics/SolidMechanics/CohesiveZone/additional_generate_output='normal_traction tangent_traction normal_jump tangent_jump'  Physics/SolidMechanics/CohesiveZone/material_output_order=FIRST Outputs/file_base=czm_multiple_action_common_output_out"
  mesh_mode = 'REPLICATED'
  recover = false
  allow_test_objects = true
[]

[czm_multiple_action_additional_outputs_on_different_block]
  type = Exodiff
  input = 'czm_multiple_action_and_materials.i'
  exodiff = czm_multiple_action_additional_outputs_on_different_block_out.e
  detail = "The system shall add additional outputs requested within each block"
  cli_args = "Physics/SolidMechanics/CohesiveZone/generate_output='traction_x traction_y traction_z jump_x jump_y jump_z' Physics/SolidMechanics/CohesiveZone/czm_ik_012/additional_generate_output='normal_traction tangent_traction' Physics/SolidMechanics/CohesiveZone/czm_ik_23/additional_generate_output='normal_jump tangent_jump' Outputs/file_base=czm_multiple_action_additional_outputs_on_different_block_out"
  mesh_mode = 'REPLICATED'
  recover = false
  allow_test_objects = true
[]

[multiple_action_differnt_outputs_orders_on_different_block]
  type = Exodiff
  input = 'czm_multiple_action_and_materials.i'
  exodiff = czm_multiple_action_differnt_outputs_orders_on_different_block_out.e
  detail = "The system shall be able to use different output orders on different blocks"
  cli_args = "Physics/SolidMechanics/CohesiveZone/czm_ik_012/material_output_order=CONSTANT Physics/SolidMechanics/CohesiveZone/czm_ik_23/material_output_order=FIRST Outputs/file_base=czm_multiple_action_differnt_outputs_orders_on_different_block_out"
  mesh_mode = 'REPLICATED'
  recover = false
  allow_test_objects = true
[]

[action_3D_small_strain_normal]
  type = Exodiff
  input = czm_multiple_dimension_base.i
  exodiff = czm_action_3D_small_strain_normal_out.e
  cli_args = 'BCs/right_x/function=1*t BCs/right_y/function=0*t BCs/right_z/function=0*t Outputs/out/file_base=czm_action_3D_small_strain_normal_out Mesh/msh/nx=2 Mesh/msh/ny=2 Mesh/msh/nz=2 Mesh/msh/dim=3'
  detail = 'The system shall allow for cohesive zone laws to represent the traction-separation behavior at an interface between two bodies represented by continuum elements in 3D using the PureElasticTractionSeparation model, and only compute a normal gap under purely normal loading.'
  mesh_mode = 'REPLICATED'
  recover = false
  abs_zero = 1e-6
[]

[action_2D_small_strain_normal]
  type = Exodiff
  input = czm_multiple_dimension_base.i
  exodiff = czm_action_2D_small_strain_normal_out.e
  cli_args = "BCs/right_x/function=1*t BCs/right_y/function=0*t BCs/inactive='left_z right_z' Outputs/out/file_base=czm_action_2D_small_strain_normal_out Mesh/msh/nx=2 Mesh/msh/ny=2 Mesh/msh/dim=2 Mesh/subdomain_1/top_right='0.5 1 0' Mesh/subdomain_2/top_right='1 1 0' GlobalParams/displacements='disp_x disp_y' Physics/SolidMechanics/CohesiveZone/czm1/displacements='disp_x disp_y' Materials/czm_mat/displacements='disp_x disp_y'"
  detail = 'The system shall allow for cohesive zone laws to represent the traction-separation behavior at an interface between two bodies represented by continuum elements in 2D using an elastic traction sepration law, and only compute a normal gap under purely normal loading when using the small strain formulation.'
  mesh_mode = 'REPLICATED'
  recover = false
  abs_zero = 1e-6
[]

[action_1D_small_strain_normal]
  type = Exodiff
  input = czm_multiple_dimension_base.i
  exodiff = czm_action_1D_small_strain_normal_out.e
  cli_args = "BCs/right_x/function=1*t BCs/inactive='left_y right_y left_z right_z' Outputs/out/file_base=czm_action_1D_small_strain_normal_out Mesh/msh/nx=2  Mesh/msh/dim=1 Mesh/subdomain_1/top_right='0.5 0 0' Mesh/subdomain_2/top_right='1 0 0' GlobalParams/displacements='disp_x' Physics/SolidMechanics/CohesiveZone/czm1/displacements='disp_x' Physics/SolidMechanics/CohesiveZone/czm1/generate_output='traction_x normal_traction tangent_traction jump_x normal_jump tangent_jump' Materials/czm_mat/displacements='disp_x'"
  detail = 'The system shall allow for cohesive zone laws to represent the traction-separation behavior at an interface between two bodies represented by continuum elements in 1D using an elastic traction sepration law, and only compute a normal gap under purely normal loading when using the small strain formulation.'
  mesh_mode = 'REPLICATED'
  recover = false
[]

[action_3D_total_lagrangian_normal]
  type = Exodiff
  input = czm_multiple_dimension_base.i
  exodiff = czm_action_3D_total_lagrangian_normal_out.e
  cli_args = 'BCs/right_x/function=1*t BCs/right_y/function=0*t BCs/right_z/function=0*t Outputs/out/file_base=czm_action_3D_total_lagrangian_normal_out Mesh/msh/nx=2 Mesh/msh/ny=2 Mesh/msh/nz=2 Mesh/msh/dim=3 Physics/SolidMechanics/CohesiveZone/czm1/strain=FINITE'
  detail = 'The system shall allow for cohesive zone laws to represent the traction-separation behavior at an interface between two bodies represented by continuum elements in 3D using an elastic traction sepration law, and only compute a normal gap under purely normal loading when using the total Lagrangian formulation.'
  mesh_mode = 'REPLICATED'
  recover = false
  abs_zero = 1e-6
[]

[action_2D_total_lagrangian_normal]
  type = Exodiff
  input = czm_multiple_dimension_base.i
  exodiff = czm_action_2D_total_lagrangian_normal_out.e
  cli_args = "BCs/right_x/function=1*t BCs/right_y/function=0*t BCs/inactive='left_z right_z' Outputs/out/file_base=czm_action_2D_total_lagrangian_normal_out Mesh/msh/nx=2 Mesh/msh/ny=2 Mesh/msh/dim=2 Mesh/subdomain_1/top_right='0.5 1 0' Mesh/subdomain_2/top_right='1 1 0' GlobalParams/displacements='disp_x disp_y' Physics/SolidMechanics/CohesiveZone/czm1/displacements='disp_x disp_y' Materials/czm_mat/displacements='disp_x disp_y'  Physics/SolidMechanics/CohesiveZone/czm1/strain=FINITE"
  detail = 'The system shall allow for cohesive zone laws to represent the traction-separation behavior at an interface between two bodies represented by continuum elements in 2D using an elastic traction sepration law, and only compute a normal gap under purely normal loading when using the total Lagrangian formulation.'
  mesh_mode = 'REPLICATED'
  recover = false
  abs_zero = 1e-6
[]

[action_1D_total_lagrangian_normal]
  type = Exodiff
  input = czm_multiple_dimension_base.i
  exodiff = czm_action_1D_total_lagrangian_normal_out.e
  cli_args = "BCs/right_x/function=1*t BCs/inactive='left_y right_y left_z right_z' Outputs/out/file_base=czm_action_1D_total_lagrangian_normal_out Mesh/msh/nx=2  Mesh/msh/dim=1 Mesh/subdomain_1/top_right='0.5 0 0' Mesh/subdomain_2/top_right='1 0 0' GlobalParams/displacements='disp_x' Physics/SolidMechanics/CohesiveZone/czm1/displacements='disp_x' Physics/SolidMechanics/CohesiveZone/czm1/generate_output='traction_x normal_traction tangent_traction jump_x normal_jump tangent_jump' Materials/czm_mat/displacements='disp_x' Physics/SolidMechanics/CohesiveZone/czm1/strain=FINITE"
  detail = 'The system shall allow for cohesive zone laws to represent the traction-separation behavior at an interface between two bodies represented by continuum elements in 1D using an elastic traction sepration law, and only compute a normal gap under purely normal loading when using the total Lagrangian formulation.'
  mesh_mode = 'REPLICATED'
  recover = false
[]

[czm_incremental_stretch_rotate_total_lagrangian]
  type = Exodiff
  input = 'stretch_rotate_large_deformation.i'
  exodiff = 'czm_incremental_stretch_rotate_total_lagrangian_out.e'
  cli_args = "Materials/czm_mat/type=PureElasticTractionSeparationIncremental  Outputs/file_base=czm_incremental_stretch_rotate_total_lagrangian_out"
  detail = "The CZMInterfaceKernelTotalLagrangian should properly account for large area changes and interface rotations when using an incremental material."
  mesh_mode = 'REPLICATED'
  recover = false
  allow_test_objects = true
[]

[czm_total_stretch_rotate_total_lagrangian]
  type = Exodiff
  input = 'stretch_rotate_large_deformation.i'
  exodiff = 'czm_total_stretch_rotate_total_lagrangian_out.e'
  cli_args = "Outputs/file_base=czm_total_stretch_rotate_total_lagrangian_out"
  detail = "The CZMInterfaceKernelTotalLagrangian should properly account for large area changes and interface rotations when using a total strain material."
  mesh_mode = 'REPLICATED'
  recover = false
[]

[./except1]
  type = 'RunException'
  input = 'non-isotropic_error_test.i'
  expect_err = 'CriticalTimeStep can only be used with elasticity tensor materials that guarantee isotropic tensors.'
  requirement = "The system shall produce an error if the input elasticity tensor is non-isotropic."
[../]

[test_adapt_err]
  type = 'RunException'
  input = 'test_adapt.i'
  cli_args = 'VectorPostprocessors/section_output/require_equal_node_counts=true Outputs/file_base=test_adapt_err_out'
  expect_err = "Node counts do not match for all axial positions. If this behavior is desired to accommodate nonuniform meshes, set 'require_equal_node_counts=false'"
  requirement = 'The system shall compute the average of the nodal displacements of a cross section at automatically located points that are updated with mesh adaptivity.'
[]

[test_one_step_empty_positions_err]
  type = 'RunException'
  input = 'test_one_step.i'
  cli_args = "VectorPostprocessors/section_output/positions='' Outputs/file_base=err_pos_out"
  expect_err = "If 'positions' is specified, at least one value must be provided"
  requirement = 'The system shall generate an error if the user defines an empty set of positions in AverageSectionValueSampler.'
[]

[test_one_step_nonexistent_variable]
  type = 'RunException'
  input = 'test_one_step.i'
  cli_args = "VectorPostprocessors/section_output/variables='dummy' Outputs/file_base=err_var_out"
  expect_err = "One or more of the specified variables do not exist in this system."
  requirement = 'The system shall generate an error if the user specifies a nonexistent variable in AverageSectionValueSampler.'
[]

[volumetric_eigenstrain_density_check]
  type = 'RunException'
  input = 'volumetric_eigenstrain_increase.i'
  cli_args = 'Materials/void_density/expression="-0.01"'
  expect_err = 'A negative, non-physical spherical void number density has been provided: -0.01. This value is non-physical and violates the assumptions of this eigenstrain class'
  requirement = 'The system shall not attempt to calculate the volumetric eigenstrain from a negative number density of voids.'
  issues = '#25932'
  design = '/ComputeCrystalPlasticityVolumetricEigenstrain.md'
[]

[volumetric_eigenstrain_radius_check]
  type = 'RunException'
  input = 'volumetric_eigenstrain_increase.i'
  cli_args = 'Materials/void_radius/prop_values="-1.0e-6"'
  expect_err = 'A negative mean spherical void radius value, -1e-06, has been provided; this value is non-physical and violates the assumptions of this eigenstrain class'
  requirement = 'The system shall not attempt to calculate the volumetric eigenstrain from a negative average radius value.'
  issues = '#25932'
  design = '/ComputeCrystalPlasticityVolumetricEigenstrain.md'
[]

[euler_angle_auxvar_conflict]
  type = 'RunException'
  input = 'euler_angle_conflict.i'
  expect_err = 'Euler angles cannot be supplied from both coupled variables and auxiliary input file in the option `read_prop_user_object`.'
  requirement = 'An error will be thrown when trying to set Euler angles information from a text file and AuxVariables at the same time.'
[]

[euler_angle_components]
  type = 'RunException'
  input = 'euler_angle_auxvar.i'
  expect_err = 'The Euler angles should have three components.'
  cli_args = 'Materials/elasticity_tensor/euler_angle_variables="euler_angle_1"'
  requirement = 'An error will be thrown when trying to set Euler angles information with less than three components.'
[]

[update_method_hcp_aprismatic_capyramidal]
  type = 'CSVDiff'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  csvdiff = 'update_method_hcp_aprismatic_capyramidal_out.csv'
  cli_args = 'Outputs/out/type=Checkpoint Outputs/out/execute_on=Final'
  allow_warnings = true
  detail = 'compute the slip system strength and resolved applied shear stress for the prismatic<a> and pyramidal<c+a> slip systems when loaded along the c-axis;'
[]

[update_method_hcp_aprismatic_active]
  type = 'CSVDiff'
  input = 'update_method_hcp_aprismatic_active.i'
  csvdiff = 'update_method_hcp_aprismatic_active_out.csv'
  allow_warnings = true
  detail = 'evolve the forest dislocation and substructure densities on prismatic<a> slip planes when loaded in an appropriate orientation;'
[]

[update_method_hcp_capyramidal_active]
  type = 'CSVDiff'
  input = 'update_method_hcp_capyramidal_active.i'
  csvdiff = 'update_method_hcp_capyramidal_active_out.csv'
  allow_warnings = true
  detail = 'evolve the forest dislocation densities on pyramidal <c+a> slip planes and the substructure density when the crystal is loaded in an appropriate orientation.'
[]

[update_method_hcp_basal_active]
  type = 'CSVDiff'
  input = 'update_method_hcp_basal_active.i'
  csvdiff = 'update_method_hcp_basal_active_out.csv'
  cli_args = 'Outputs/out/type=Checkpoint Outputs/out/execute_on=Final'
  allow_warnings = true
  detail = 'evolve the forest dislocation densities on basal <a> slip planes when the crystal is loaded in an appropriate orientation.'
[]

[update_method_hcp_no_substructure]
  type = 'CSVDiff'
  input = 'update_method_hcp_no_substructure.i'
  csvdiff = 'update_method_hcp_no_substructure_out.csv'
  allow_warnings = true
  detail = 'prevent the calculation of a negative substructure density increment.'
[]

[update_method_hcp_no_negative_aprismatic]
  type = 'CSVDiff'
  input = 'update_method_hcp_no_negative_aprismatic.i'
  csvdiff = 'update_method_hcp_no_negative_aprismatic_out.csv'
  allow_warnings = true
  detail = 'prevent the calculation of a negative total forest density.'
[]

[update_method_hcp_size_slip_modes]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/number_slip_systems_per_mode="3 3 9"'
  expect_err = 'The size the number of slip systems per mode is not equal to the number of slip system types'
  detail = 'ensure the user specifies a number of the slip systems per mode which is equal to the number of modes;'
[]

[update_method_hcp_sum_slip_modes_vs_systems]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/number_slip_systems_per_mode="3 15"'
  expect_err = 'The number of slip systems and the sum of the slip systems in each of the slip system modes are not equal'
  detail = 'ensure the user specifies an equal total number of slip systems  and the sum of the different slip mode systems;'
[]

[update_method_hcp_size_burgers_vectors]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/burgers_vector_per_mode="2"'
  expect_err = 'Please ensure that the size of burgers_vector_per_mode equals the value supplied for slip_system_modes'
  detail = 'ensure the user specifies a number of burgers vectors per mode that is equal to the number of modes;'
[]

[update_method_hcp_size_slip_generation]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/slip_generation_coefficient_per_mode="1 2 3"'
  expect_err = 'Please ensure that the size of slip_generation_coefficient_per_mode equals the value supplied for slip_system_modes'
  detail = 'ensure the user specifies a number of slip generation coefficient per mode that is equal to the number of modes;'
[]

[update_method_hcp_size_activation_energy]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/normalized_slip_activiation_energy_per_mode="1 2 3 4"'
  expect_err = 'Please ensure that the size of normalized_slip_activiation_energy_per_mode equals the value supplied for slip_system_modes'
  detail = 'ensure the user specifies a size of slip activation energy per mode that is equal to the specified number of modes;'
[]

[update_method_hcp_size_proportionality_factor]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/slip_energy_proportionality_factor_per_mode="1"'
  expect_err = 'Please ensure that the size of slip_energy_proportionality_factor_per_mode equals the value supplied for slip_system_modes'
  detail = 'ensure the user specifies a size of slip energy proportionality factor per mode that is equal to the specified number of modes;'
[]

[update_method_hcp_size_shear_modulus]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/effective_shear_modulus_per_mode="1"'
  expect_err = 'Please ensure that the size of effective_shear_modulus_per_mode equals the value supplied for slip_system_modes'
  detail = 'ensure the user specifies a size of effective shear modulus per mode that is equal to the specified number of modes;'
[]

[update_method_hcp_size_substructure_rate]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/substructure_rate_coefficient_per_mode="3 2 1"'
  expect_err = 'Please ensure that the size of substructure_rate_coefficient_per_mode equals the value supplied for slip_system_modes'
  detail = 'ensure the user specifies a size of substructure generation rate coefficient per mode that is equal to the specified number of modes;'
[]

[update_method_hcp_size_hall_petch]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/Hall_Petch_like_constant_per_mode="1"'
  expect_err = 'Please ensure that the size of Hall_Petch_like_constant_per_mode equals the value supplied for slip_system_modes'
  detail = 'ensure the user specifies a size of Hall-Petch like coefficient per mode that is equal to the specified number of modes;'
[]

[update_method_hcp_size_lattice_friction]
  type = 'RunException'
  input = 'update_method_hcp_aprismatic_capyramidal.i'
  cli_args = 'Materials/trial_xtalpl/lattice_friction_per_mode="3"'
  expect_err = 'Please ensure that the size of lattice_friction_per_mode equals the value supplied for slip_system_modes'
  detail = 'ensure the user specifies a size of lattice friction values per mode that is equal to the specified number of modes;'
[]

[update_method_hcp_representative_slip_systems]
  type = 'CSVDiff'
  input = 'update_method_hcp_representative_slip_systems.i'
  csvdiff = 'update_method_hcp_representative_slip_systems_out.csv'
  allow_warnings = true
  detail = ' compute the Schmid tensors and resulting resolved applied shear stress for a non-physical representative set of basal, prismatic, and pyramidal slip systems;'
[]

[update_method_hcp_representative_equal_lattice_parameters]
  type = 'RunException'
  input = 'update_method_hcp_representative_slip_systems.i'
  cli_args = 'Materials/trial_xtalpl/unit_cell_dimension="3e-7 3e-7 3e-7"'
  expect_err = 'The specified unit cell dimensions are not consistent with expectations for HCP crystal hexagonal lattices'
  detail = 'ensure that the provided lattice parameters match the expectation a = b != c for HCP crystals;'
[]

[update_method_hcp_representative_incorrect_index_numbers]
  type = 'RunException'
  input = 'update_method_hcp_representative_slip_systems.i'
  cli_args = 'Materials/trial_xtalpl/slip_sys_file_name=input_slip_sys.txt Materials/trial_xtalpl/number_slip_systems=12'
  expect_err = 'CrystalPlasticityStressUpdateBase Error: The number of entries in the first row of the slip system file is not consistent with the expectations for the 4-index Miller-Bravais assumption for HCP crystals. This file should represent both the slip plane normal and the slip direction with 4-indices each.'
  detail = 'ensure that eight indices are given for the Miller-Bravais notation for HCP slip systems;'
[]

[update_method_hcp_representative_incorrect_basal_plane]
  type = 'RunException'
  input = 'update_method_hcp_representative_slip_systems.i'
  cli_args = 'Materials/trial_xtalpl/slip_sys_file_name=error_slip_sys_hcp.txt'
  expect_err = 'CrystalPlasticityStressUpdateBase Error: The specified HCP basal plane Miller-Bravais indices do not sum to zero. Check the values supplied in the associated text file.'
  detail = ' ensure that the Miller-Bravais basal plane indices sum to zero.'
[]

[update_method_hcp_representative_incorrect_basal_slip_direction]
  type = 'RunException'
  input = 'update_method_hcp_representative_slip_systems.i'
  cli_args = 'Materials/trial_xtalpl/slip_sys_file_name=error_slip_dir_sys_hcp.txt Materials/trial_xtalpl/number_slip_systems=1'
  expect_err = 'CrystalPlasticityStressUpdateBase Error: The specified HCP slip direction Miller-Bravais indices in the basal plane '
  detail = ' ensure that the Miller-Bravais indices for the slip direction  in the basal plane sum to zero.'
[]

[update_method_hcp_convergence_issue_flag]
  type = 'RunException'
  input = 'update_method_hcp_convergence_issue_flag.i'
  expect_err = 'Maximum allowable slip increment exceeded'
  detail = ' provide warning messages about constitutive model and slip increment convergence issues when the corresponding input file parameter is set to true.'
[]

[update_method_hcp_stress_convergence_issue_flag]
  type = 'RunException'
  input = 'update_method_hcp_convergence_issue_flag.i'
  cli_args = 'Materials/trial_xtalpl/slip_increment_tolerance=2.0 Materials/stress/maxiter=1'
  expect_err = 'ComputeMultipleCrystalPlasticityStress: Stress Integration error rmax ='
  detail = 'provide warning messages about constitutive model stress convergence issues with the appropriate input file parameter setting.'
[]

[rotation_matrix_and_euler_angles_error]
  type = RunException
  input = 'rotation_matrix_update_euler_angle_111_orientation.i'
  cli_args = 'Materials/elasticity_tensor/euler_angle_2=50'
  expect_err = 'Bunge Euler angle information and the rotation_matrix cannot both be specified. Provide only one type of orientation information in the input file.'
  requirement = 'The system shall not allow specification of both the rotation matrix and the Bunge Euler angles within the same input file.'
  issues = '#23451'
  design = '/ComputeElasticityTensorCP.md'
[]

[only_twinning_fcc]
  type = 'CSVDiff'
  input = 'only_twinning_fcc.i'
  csvdiff = 'only_twinning_fcc_out.csv'
  allow_warnings = true
  cli_args = 'Outputs/out/type=Checkpoint Outputs/out/execute_on=Final'
  detail = 'compute the plastic shear due to twinning as a function of the applied shear stress and resistance to twin propagation, the corresponding twin volume fraction per twin system, and the total value of the twin volume fraction across all twin systems.'
[]

[check_direction_twin_propagation]
  type = 'CSVDiff'
  input = 'check_direction_twin_propagation.i'
  csvdiff = 'check_direction_twin_propagation_out.csv'
  allow_warnings = true
  detail = 'compute a non-zero twin volume fraction, and associate plastic shear increment, for only postive applied shear stress values.'
[]

[coplanar_twin_hardening]
  type = 'CSVDiff'
  input = 'coplanar_twin_hardening.i'
  csvdiff = 'coplanar_twin_hardening_out.csv'
  detail = 'calculate the coplanar twin propogation resistance contribution as a function of the twin volume fractions on coplanar systems alone.'
[]

[non_coplanar_twin_hardening]
  type = 'CSVDiff'
  input = 'non_coplanar_twin_hardening.i'
  csvdiff = 'non_coplanar_twin_hardening_out.csv'
  detail = 'calculate the contribution to twin propagation resistance from non-coplanar systems as a function of only the twin volume fraction on non-coplanar twin systems.'
[]

[upper_twin_fraction_limit]
  type = CSVDiff
  input = 'upper_twin_fraction_limit.i'
  csvdiff = 'upper_twin_fraction_limit_out.csv'
  allow_warnings = true
  detail = 'set the plastic shear increments due to twinning to zero once the total twin volume fraction limit is reached.'
[]

[upper_twin_fraction_limit_error]
  type = RunException
  input = 'upper_twin_fraction_limit.i'
  cli_args = 'Materials/twin_only_xtalpl/print_state_variable_convergence_error_messages=true Materials/twin_only_xtalpl/initial_twin_lattice_friction=3.0 Executioner/dtmin=0.025 Outputs/csv=false'
  expect_err = 'Maximum allowable twin volume fraction limit exceeded with a value of'
  detail = 'prevent the accumulation of a total twin volume fraction above the user-defined limit.'
[]

[strain_increment_upper_limit]
  type = RunException
  input = 'upper_twin_fraction_limit.i'
  cli_args = 'BCs/tdisp/function=5.0e-3*t'
  expect_err = 'Maximum allowable plastic slip increment due to twinning exceeded'
  detail = 'prevent the calculation of a plastic shear increment in excess of the user-defined tolerance value.'
[]

[./test_ii_arb_eig_grad_err_jint]
  prereq = test_ii_arb_eig_grad
  type = 'RunException'
  input = 'interaction_integral_2d_eig_grad.i'
  cli_args = "DomainIntegral/integrals='JIntegral InteractionIntegralKI'"
  expect_err = "'eigenstrain_gradient' cannot be specified when the computed integrals include JIntegral, CIntegral, or KFromJIntegral"
  requirement = 'The DomainIntegralAction shall generate an error if a user specifies the eigenstrain gradient and also request the J integral'
  allow_test_objects = true
  issues = '#18804'
[../]

[./test_ii_arb_eig_grad_err_temp]
  prereq = test_ii_arb_eig_grad_err_jint
  type = 'RunException'
  input = 'interaction_integral_2d_eig_grad.i'
  cli_args = "DomainIntegral/temperature=temp DomainIntegral/eigenstrain_names=thermal_expansion"
  expect_err = "'eigenstrain_gradient' cannot be specified together with 'temperature' or 'eigenstrain_names'. These are for separate, mutually exclusive systems for including the effect of eigenstrains"
  requirement = 'The DomainIntegralAction shall generate an error if a user specifies the temperature and also provides the eigenstrain_gradient'
  allow_test_objects = true
  issues = '#18804'
[../]

[./test_ii_bf_err_jint]
  prereq = test_ii_bf
  type = 'RunException'
  input = 'interaction_integral_2d_bf.i'
  cli_args = "DomainIntegral/integrals='JIntegral InteractionIntegralKI'"
  expect_err = "'body_force' cannot be specified when the computed integrals include JIntegral, CIntegral, or KFromJIntegral"
  requirement = 'The DomainIntegralAction shall generate an error if a user specifies the body force and also request the J integral'
  issues = '#18804'
[../]

[stress_errorcheck]
  type = 'RunException'
  input = 'finite_strain_stress_errorcheck.i'
  expect_err = 'This linear elastic stress calculation only works for small strains; use ComputeFiniteStrainElasticStress for simulations using incremental and finite strains.'
  issues = '#4716 #7555'
  design = 'source/materials/ComputeFiniteStrainElasticStress.md source/materials/ComputeLinearElasticStress.md syntax/Physics/SolidMechanics/QuasiStatic/index.md'
  requirement = 'The ComputeLinearElasticStress class shall generate an error if a user attempts to run a problem using ComputeLinearElasticStress with a finite strain formulation.'
[]

[./except01]
  type = 'RunException'
  input = 'except01.i'
  expect_err = 'ComputeEigenstrainFromInitialStress: 9 initial stress functions must be provided.  You supplied 3'
  design = 'ComputeEigenstrainFromInitialStress.md'
  issues = "#9749"
  requirement = "SolidMechanics shall allow users to specify initial stresses, but shall error-out with appropriate message if the user does not supply the correct number of functions to define the initial stress tensor"
[../]

[./except02]
  type = 'RunException'
  input = 'except02.i'
  expect_err = 'ComputeEigenstrainFromInitialStress: If you supply initial_stress_aux, 9 values must be given.  You supplied 3'
  design = 'ComputeEigenstrainFromInitialStress.md'
  issues = "#13087"
  requirement = "SolidMechanics shall allow users to specify initial stresses, but shall error-out with appropriate message if the user does not supply the correct number of AuxVariables to define the initial stress tensor"
[../]

[j_2d_block_restrict_error]
  type = 'RunException'
  issues = '#24795'
  input = 'j_integral_2d_block_restrict.i'
  cli_args = "DomainIntegral/block='1 2'"
  expect_err = "'EshelbyTensor' is not defined on block lowerd_dummy."
  requirement = 'The domain integral action shall error out if the blocks selected do not have the adequate material objects.'
[]

[j_int_fgm_sif_error]
  type = 'RunException'
  input = 'j_int_fgm_sif.i'
  expect_err = 'You have selected to compute the interaction integral for a crack in FGM'
  requirement = "The Domain Integral shall error out if the user input does not consistently define the required material properties to consider extra interaction integral terms for functionally graded materials."
[]

[./with_correction]
  issues = '#20281'
  type = RunException
  input = 'simple_test.i'
  design = 'source/actions/TensorMechanicsAction.md'
  requirement = "Cannot run second order problem with stabilization"
  cli_args = "Physics/SolidMechanics/QuasiStatic/all/volumetric_locking_correction=true"
  expect_err = "Volumetric locking correction should not be used for higher-order elements."
[../]

[./stvenantkirchhoff]
  type = RunException
  input = stvenantkirchhoff.i
  expect_err = 'ComputeStVenantKirchhoffStress requires an isotropic elasticity'
  requirement = 'Throw an error if you do not provide an isotropic tensor to ComputeStVenantKirchhoffStress'
  issues = '#17472'
  design = 'source/materials/lagrangian/ComputeStVenantKirchhoffStress.md'
[../]

[./elements_damage_limit_unknown_prop]
  type = 'RunException'
  prereq = elements_damage_limit
  input = 'elements_changed_timestep_limit.i'
  cli_args = 'Postprocessors/time_step_limit/elements_changed_property=unknown_property' #Point it to an unknown property
  expect_err = "Material property 'unknown_property', requested by 'time_step_limit' is not defined on block 0"
  requirement = "The MaterialTimeStepPostprocessor shall generate an error if an unknown property is requested with the 'time_step_limit' parameter"
  issues = '#11662'
  design = 'MaterialTimeStepPostprocessor.md'
[../]

[./e]
  type = 'RunException'
  prereq = scalar_damage_material_limit
  input = 'scalar_material_damage_timestep_limit.i'
  cli_args = 'Postprocessors/time_step_limit/use_material_timestep_limit=false' #disable botg matl_timestep and elements count
  expect_err = "either use_material_timestep_limit needs to be true or elements_changed_property defined"
  requirement = "The MaterialTimeStepPostprocessor shall generate an error if neither the material time step limit nor the elements changed limit is specified."
  issues = '#11662'
  design = 'MaterialTimeStepPostprocessor.md'
[../]

[./except1]
  type = 'RunException'
  input = 'except1.i'
  expect_err = 'Mohr-Coulomb friction angle must not be less than Mohr-Coulomb dilation angle'
[../]

[./except2]
  type = 'RunException'
  input = 'except2.i'
  expect_err = 'mc_lode_cutoff must not be negative'
[../]

[./except3]
  type = 'RunException'
  input = 'except3.i'
  expect_err = 'Mohr-Coulomb edge smoothing angle is too small and a non-convex yield surface will result.  Please choose a larger value'
[../]

[./except4]
  type = 'RunException'
  input = 'except4.i'
  expect_err = 'Mohr-Coulomb edge smoothing angle is too small and a non-convex yield surface will result.  Please choose a larger value'
[../]

[multipowerlaw_user_input_error]
  type = 'RunException'
  input = 'power_law_creep.i'
  expect_err = 'Stress thresholds input must be provided in increasing ordered'
  cli_args = "Materials/creep_one/type=ADMultiplePowerLawCreepStressUpdate Materials/creep_one/stress_thresholds='0 53000 52000' Materials/creep_one/coefficient='1e-24 1.92307692e-29' Materials/creep_one/n_exponent='4 9.0' Materials/creep_one/m_exponent='0 0' Materials/creep_one/activation_energy='0 0'
              Materials/creep_nine/type=ADMultiplePowerLawCreepStressUpdate Materials/creep_nine/stress_thresholds='0 53000 300000' Materials/creep_nine/coefficient='9e-24 3.32840237e-33' Materials/creep_nine/n_exponent='4 9.0' Materials/creep_nine/m_exponent='0 0' Materials/creep_nine/activation_energy='0 0'
              Outputs/file_base=multi_power_law_creep_two_intervals_different_no_continuity_out"

  requirement = "Stress thresholds must be provided in increasing order, otherwise an error will inform the user of such a requirement. "
[]

[moose_to_neml2_parameter_error]
  type = RunException
  input = 'moose_to_neml2_parameters.i'
  cli_args = 'UserObjects/gather_E/neml2_parameter=p'
  requirement = 'The system shall error out while trying to set a parameter that does not exist.'
  expect_err = 'there is not such parameter in the NEML2 material model.'
  required_submodule = 'contrib/neml2'
  libtorch = true
  recover = false
[]

[zero]
  type = CSVDiff
  input = 'negative_porosity.i'
  csvdiff = 'negative_porosity_zero_out.csv'
  cli_args = 'Outputs/file_base=negative_porosity_zero_out Materials/porosity/negative_behavior=ZERO'
  detail = 'by setting the porosity to zero.'
[]

[initial]
  type = CSVDiff
  input = 'negative_porosity.i'
  csvdiff = 'negative_porosity_initial_out.csv'
  cli_args = 'Outputs/file_base=negative_porosity_initial_out Materials/porosity/negative_behavior=INITIAL_CONDITION'
  detail = 'by setting the porosity to the initial condition.'
[]

[exception]
  type = RunException
  input = 'negative_porosity.i'
  expect_err = 'porosity is negative'
  cli_args = 'Materials/porosity/negative_behavior=EXCEPTION'
  detail = 'by throwing an exception.'
[]

[exception]
  type = RunException
  input = 'nan.i'
  expect_err = 'porosity is nan'
  detail = 'by throwing an exception.'
[]

[./sphere_1d_spherical_except_1]
  type = 'RunException'
  input = 'sphere_1d_spherical.i'
  cli_args = "AuxKernels/raddispaux/origin='0 0 0'"
  expect_err = "The 'origin' parameter is only valid for models with Cartesian or axisymmetric coordinate systems."
  prereq = sphere_1d_spherical
  design = 'RadialDisplacementSphereAux.md'
  requirement = 'The system shall report an error if "origin" is supplied to RadialDisplacementSphereAux when the coordinate system is not Cartesian or axisymmetric'
[../]

[./cylinder_2d_axisymmetric_except_1]
  type = 'RunException'
  input = 'cylinder_2d_axisymmetric.i'
  cli_args = "AuxKernels/raddispaux/axis_vector='0 0 1'"
  expect_err = "The 'axis_vector' parameter is only valid for 3D Cartesian models."
  prereq = cylinder_2d_axisymmetric
  design = 'RadialDisplacementCylinderAux.md'
  requirement = 'The system shall report an error if "axis_vector" is supplied to RadialDisplacementCylinderAux and the model is not 3D Cartesian'
[../]

[./cylinder_3d_cartesian_except_1]
  type = 'RunException'
  input = 'cylinder_3d_cartesian.i'
  cli_args = "AuxKernels/raddispaux/axis_vector='0 0 0'"
  expect_err = "axis_vector must have nonzero length"
  prereq = cylinder_3d_cartesian
  design = 'RadialDisplacementCylinderAux.md'
  requirement = 'The system shall report an error if "axis_vector" has length of zero'
[../]

[isotropic_plasticity_error1]
  type = 'RunException'
  input = 'isotropic_plasticity_errors.i'
  cli_args = 'Materials/isotropic_plasticity/hardening_function=harden_func
                Materials/isotropic_plasticity/hardening_constant=2
                Materials/isotropic_plasticity/yield_stress=50
                Materials/elasticity_tensor/type=ComputeIsotropicElasticityTensor Materials/elasticity_tensor/youngs_modulus=2.1e5
                Materials/elasticity_tensor/poissons_ratio=0.0'
  expect_err = 'Only the hardening_constant or only the hardening_function can be defined but not both'
  design = 'IsotropicPlasticityStressUpdate.md'
  requirement = 'The system shall only allow the calculation of a J2 isotropic plasticity stress response with only one but not both of a hardening function or hardening constant specified to define the evolving yield surface.'
[]

[isotropic_plasticity_error2]
  type = 'RunException'
  input = 'isotropic_plasticity_errors.i'
  cli_args = 'Materials/isotropic_plasticity/yield_stress=50
                Materials/elasticity_tensor/type=ComputeIsotropicElasticityTensor Materials/elasticity_tensor/youngs_modulus=2.1e5
                Materials/elasticity_tensor/poissons_ratio=0.0'
  expect_err = 'Either hardening_constant or hardening_function must be defined'
  prereq = 'isotropic_plasticity_error1'
  design = 'IsotropicPlasticityStressUpdate.md'
  requirement = 'The system shall only calculate the J2 isotropic plasticity stress response when either a hardening function or a hardening constant is specified to define the evolving yield surface.'
[]

[isotropic_plasticity_error3]
  type = 'RunException'
  input = 'isotropic_plasticity_errors.i'
  cli_args = 'Materials/isotropic_plasticity/hardening_function=harden_func
                Materials/elasticity_tensor/type=ComputeIsotropicElasticityTensor Materials/elasticity_tensor/youngs_modulus=2.1e5
                Materials/elasticity_tensor/poissons_ratio=0.0'
  expect_err = 'Either yield_stress or yield_stress_function must be given'
  prereq = 'isotropic_plasticity_error2'
  design = 'IsotropicPlasticityStressUpdate.md'
  requirement = 'The system shall only calculate the J2 isotropic plasticity stress response when either a constant yield stress or a yield stress function is specified to define the initial yield surface.'
[]

[isotropic_plasticity_error4]
  type = 'RunException'
  input = 'isotropic_plasticity_errors.i'
  cli_args = 'Materials/isotropic_plasticity/hardening_function=harden_func
                Materials/elasticity_tensor/type=ComputeIsotropicElasticityTensor Materials/elasticity_tensor/youngs_modulus=2.1e5
                Materials/elasticity_tensor/poissons_ratio=0.0
                Materials/isotropic_plasticity/yield_stress=-5.0'
  expect_err = 'Yield stress must be greater than zero'
  prereq = 'isotropic_plasticity_error3'
  design = 'IsotropicPlasticityStressUpdate.md'
  requirement = 'The system shall return an error if a negative yield stress is supplied when calculating the J2 isotropic plasticity stress response.'
[]

[isotropic_plasticity_error5]
  type = 'RunException'
  input = 'isotropic_plasticity_errors.i'
  cli_args = 'Materials/isotropic_plasticity/hardening_function=harden_func
                Materials/isotropic_plasticity/yield_stress=50
                Materials/elasticity_tensor/type=ComputeElasticityTensor Materials/elasticity_tensor/fill_method=symmetric21
                Materials/elasticity_tensor/C_ijkl="1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21" '
  expect_err = 'Model isotropic_plasticity requires an isotropic elasticity tensor, but the one supplied is not guaranteed isotropic'
  prereq = 'isotropic_plasticity_error4'
  design = 'IsotropicPlasticityStressUpdate.md'
  requirement = 'The system shall return an error if anisotropic elasticity tensor is supplied when the J2 isotropic plasticity stress response calculation is requested.'
[]

[isotropic_plasticity_substepping_error]
  type = 'RunException'
  input = 'isotropic_plasticity_finite_strain.i'
  cli_args = 'Materials/isotropic_plasticity/use_substep=true'
  expect_err = 'Usage of substepping has been requested, but the inelastic model'
  design = 'RadialReturnStressUpdate.md'
  issues = '#14757'
  requirement = 'The system shall return an error if the user requests substepping for a radial return inelastic model that does not implement it yet.'
[]

[3tile_strain]
  type = 'CSVDiff'
  input = '3tile_strain.i'
  csvdiff = '3tile_strain_out.csv'
  cli_args = 'Materials/rom_stress_prediction/finite_difference_width=1e-2'
  detail = 'and evolve over time with an good finite difference width.'
  petsc_version = '>=3.9.0'
  allow_test_objects = true
[]

[3tile_strain_small_fd]
  type = 'CSVDiff'
  input = '3tile_strain.i'
  csvdiff = '3tile_strain_small_fd_out.csv'
  cli_args = 'Outputs/file_base=3tile_strain_small_fd_out Materials/rom_stress_prediction/finite_difference_width=1e-20'
  detail = 'and evolve over time with an poor finite difference width with more nonlinear iterations than with a good finite difference width.'
  petsc_version = '>=3.9.0'
  allow_test_objects = true
[]

[3d]
  type = 'CSVDiff'
  input = '3d.i'
  csvdiff = '3d_out.csv'
  detail = 'in 3D.'
  petsc_version = '>=3.9.0'
  allow_test_objects = true
[]

[3d_MPA]
  type = 'CSVDiff'
  prereq = 'rom/3d'
  input = '3d.i'
  cli_args = 'stress_unit=MPa E=3.3e5 p=1e-1'
  csvdiff = '3d_out.csv'
  detail = 'with units of stress other than Pascal.'
  allow_test_objects = true
[]

[2drz]
  type = 'CSVDiff'
  input = '2drz.i'
  csvdiff = '2drz_out.csv'
  detail = 'in 2DRz.'
  petsc_version = '>=3.9.0'
  allow_test_objects = true
[]

[2drz_json]
  type = 'CSVDiff'
  input = '2drz_json.i'
  csvdiff = '2drz_json_out.csv'
  detail = 'in 2DRz with a model loaded from a JSON data file.'
  petsc_version = '>=3.9.0'
[]

[nonad_verification]
  type = 'CSVDiff'
  input = 'nonad_verification.i'
  csvdiff = 'nonad_verification_out.csv'
  petsc_version = '>=3.9.0'
  max_parallel = 1
  allow_test_objects = true
  detail = 'in isolation (i.e. without a full displacement solve), and match with code-to-code comparison with a small set of input parameters not using AD.'
[]

[nonad_verification_heavy]
  type = 'CSVDiff'
  input = 'nonad_verification.i'
  csvdiff = 'nonad_verification_heavy_out.csv'
  heavy = true
  max_time = 500
  petsc_version = '>=3.9.0'
  cli_args = 'Executioner/dt=1e-5 Outputs/file_base=nonad_verification_heavy_out Outputs/interval=100'
  max_parallel = 1
  allow_test_objects = true
  detail = 'in isolation (i.e. without a full displacement solve), and match with code-to-code comparison with a large set of input parameters not using AD.'
[]

[lower_limit]
  type = 'CSVDiff'
  input = lower_limit.i
  csvdiff = lower_limit_out.csv
  allow_test_objects = true
  detail = 'at the lower temperature limit.'
[]

[lower_limit_extrapolated]
  type = 'CSVDiff'
  input = lower_limit.i
  csvdiff = extrapolated_out.csv
  cli_args = 'temp=500 disp=0.5e6 Outputs/file_base=extrapolated_out'
  allow_test_objects = true
  detail = 'below the lower temperature limit and extrapolate correctly.'
[]

[lower_limit_error]
  type = 'RunException'
  input = lower_limit.i
  expect_err = "In rom_stress_prediction: temperature input parameter with value \(500\) is out of lower global range"
  cli_args = 'temp=500 disp=0.5e6 Materials/rom_stress_prediction/temperature_input_window_low_failure=ERROR'
  allow_test_objects = true
  detail = 'but should error out below the lower temperature limit if requested by the user.'
[]

[upper_limit_error]
  type = 'RunException'
  input = lower_limit.i
  expect_err = "In rom_stress_prediction: temperature input parameter with value \(2500\) is out of upper global range"
  cli_args = 'temp=2500 disp=0.5e6 Materials/rom_stress_prediction/temperature_input_window_low_failure=ERROR'
  allow_test_objects = true
  detail = 'but should error out above the upper temperature limit if requested by the user.'
[]

[creep_ramp_sub_false]
  type = 'CSVDiff'
  input = creep_ramp_sub_false.i
  csvdiff = creep_ramp_sub_false_out.csv
  allow_test_objects = true
  detail = 'when substepping is off and the time step size does not guarantee the most accurate results.'
[]

[creep_ramp_sub_false_more_steps]
  type = 'CSVDiff'
  input = creep_ramp_sub_false_more_steps.i
  csvdiff = creep_ramp_sub_false_more_steps_out.csv
  allow_test_objects = true
  detail = 'when substepping is off and the time step size guarantees accurate results.'
[]

[creep_ramp_sub_true]
  type = 'CSVDiff'
  input = creep_ramp_sub_true.i
  csvdiff = creep_ramp_sub_true_out.csv
  allow_test_objects = true
  allow_warnings = true
  detail = 'when substepping is on and the numerical integration error selected by the user guarantees accurate results.'
[]

[creep_ramp_sub_true_exception]
  type = 'RunApp'
  input = creep_ramp_sub_true.i
  expect_out = "The number of substeps computed exceeds 'maximum_number_substeps'"
  cli_args = 'Materials/rom_stress_prediction/substep_strain_tolerance=5.0e-6'
  allow_test_objects = true
  allow_warnings = true
  detail = 'when substepping is on and the numerical integration error selected by the user guarantees accurate results but substepping settings require a system-wise time step cut because the computed number of substeps exceeds the maximum_number_substeps parameter.'
[]

[3tile]
  type = 'CSVDiff'
  input = '3tile.i'
  csvdiff = '3tile_out.csv'
  detail = 'for a variety of overlapping tiles and two overlapping partitions.'
  petsc_version = '>=3.9.0'
  allow_test_objects = true
[]

[3tile_json]
  type = 'CSVDiff'
  input = '3tile_json.i'
  csvdiff = '3tile_json_out.csv'
  detail = 'for a variety of overlapping tiles and two overlapping partitions with a model loaded from a JSON data file.'
  petsc_version = '>=3.9.0'
[]

[3tile_json_relative]
  type = 'CSVDiff'
  input = '3tile_json.i'
  csvdiff = '3tile_json_out.csv'
  cli_args = 'Materials/rom_stress_prediction/model=gold/3tile_json_export.json'
  detail = 'for a variety of overlapping tiles and two overlapping partitions with a model loaded from a JSON data file relative to the input.'
  petsc_version = '>=3.9.0'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[3tile_json_export]
  type = 'JSONDiff'
  input = '3tile.i'
  jsondiff = '3tile_json_export.json'
  detail = 'for a variety of overlapping tiles and two overlapping partitions, and be able to export the model data to a JSON file.'
  petsc_version = '>=3.9.0'
  cli_args = 'Materials/rom_stress_prediction/export_model=3tile_json_export.json'
  allow_test_objects = true
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[scalar_damage_material_out_of_bounds]
  type = 'RunException'
  prereq = scalar_damage_material
  input = 'scalar_material_damage.i'
  cli_args = 'Executioner/end_time=1.2'
  expect_err = 'damage_index must be between 0 and 1. Current value is:'
  requirement = "The system shall report an error when damage index is greater than 1."
[]

[ad_scalar_damage_material_out_of_bounds]
  type = 'RunException'
  prereq = ad_scalar_damage_material
  input = 'ad_scalar_material_damage.i'
  cli_args = 'Executioner/end_time=1.2'
  expect_err = 'damage_index must be between 0 and 1. Current value is:'
  requirement = "The system shall report an error when damage index is greater than 1 also when using automatic differentiation."
[]

[scalar_damage_incompatible_model]
  type = 'RunException'
  prereq = scalar_damage_material
  input = 'scalar_material_damage.i'
  cli_args = 'Materials/stress/damage_model=damage_index' #Point it to an incompatible model
  expect_err = 'Damage Model damage_index is not compatible with ComputeDamageStress'
  allow_warnings = true #Avoid failing on another warning before hitting the error we're testing for
  requirement = "The system shall make sure that the damage model is derived from DamageBase and error out when incompatible damage model is used in conjunction with ComputeDamageStress"
[]

[ad_scalar_damage_incompatible_model]
  type = 'RunException'
  prereq = ad_scalar_damage_material
  input = 'ad_scalar_material_damage.i'
  cli_args = 'Materials/stress/damage_model=damage_index' #Point it to an incompatible model
  expect_err = 'Damage Model damage_index is not compatible with ComputeDamageStress'
  allow_warnings = true #Avoid failing on another warning before hitting the error we're testing for
  requirement = "The system shall make sure that the damage model is derived from DamageBase and error out when incompatible damage model is used in conjunction with ComputeDamageStress when using automatic differentiation"
[]

[multiple_inelastic_bad_damage_model]
  type = RunException
  input = ad_scalar_material_damage_creep_power.i
  expect_err = "Damage Model.*damage_index.*is not compatible with ADComputeMultipleInelasticStress"
  requirement = "The system shall provide a meaningful error if the user selects a material model that is not recognized by the automatic differentiation multiple stress object."
  issues = '#20722'
  cli_args = 'Materials/stress/damage_model=damage_index'
  allow_warnings = True
[]

[qp_error]
  type = 'RunException'
  input = pressure_error.i
  expect_err = 'Shell element needs to have exactly four quadrature points.'

  requirement = "The system shall generate an error if the shell strain material is evaluated when the number of quadrature points does not equal 4."
[]

[cracking_err_diagonal]
  prereq = cracking
  type = RunException
  input = cracking.i
  expect_err = "Defaulting to the legacy option of 'DIAGONAL', but the 'FULL' option is preferred"
  requirement = "The ComputeSmearedCrackingStress model shall generate a warning if 'cracked_elasticity_type' is not specified"
[]

[cracking_err_softening]
  prereq = cracking_err_diagonal
  allow_warnings = true #Ignore warning about missing cracked_elasticity_type
  type = RunException
  input = cracking.i
  cli_args = "Materials/elastic_stress/softening_models='abrupt_softening abrupt_softening'"
  expect_err = "Either 1 or 3 softening models must be specified"
  requirement = "The ComputeSmearedCrackingStress model shall generate an error if an invalid number of softening models is specified"
[]

[./except01]
  type = RunException
  input = cosserat_tension.i
  cli_args = 'Kernels/x_moment/component=1 Outputs/file_base=except01'
  expect_err = 'MomentBalancing: The variable for this Kernel must be equal to the Cosserat rotation variable defined by the "component" and the "Cosserat_rotations" parameters'
  prereq = tension
[../]

[./except1]
  type = RunException
  input = except1.i
  expect_err = 'StickyBC: min_value must not be greater than max_value'
[../]

[incr_chk1]
  type = 'RunException'
  input = 'incr_model.i'
  cli_args = 'Materials/strain_energy_density/incremental=false'
  design = 'StrainEnergyDensity.md'
  expect_err = 'StrainEnergyDensity: Specified incremental = false, but material model is incremental.'
  requirement = 'The system shall be capable of informing a user when they incorrectly choose not to use the incremental strain energy density formulation with an incremental material model.'
[]

[tot_chk1]
  type = 'RunException'
  input = 'tot_model.i'
  cli_args = 'Materials/strain_energy_density/incremental=true'
  design = 'StrainEnergyDensity.md'
  expect_err = "StrainEnergyDensity: Specified incremental = true, but material model is not incremental."
  requirement = 'The system shall be capable of informing a user when they incorrectly choose to use the incremental strain energy density formulation in a material utilizing small strain.'
[]

[ad_power_law]
  type = 'RunException'
  input = 'ad_power_law_creep.i'
  expect_err = 'Exceeded maximum iterations in material return mapping iterations.'
  cli_args = 'Materials/power_law_creep/use_substepping=NONE Executioner/dtmin=0.1 Outputs/file_base=ad_power_law_out'
  requirement = 'The system shall report diverged return mapping solve due to large deformation imposed in one step using automatic differentiation'
[]

[ad_power_law_adaptive_substepping_increment]
  # We only check that the simulation converges with adaptive substepping
  type = 'RunException'
  input = 'ad_power_law_creep.i'
  cli_args = "Materials/power_law_creep/use_substepping=INCREMENT_BASED Materials/power_law_creep/adaptive_substepping=true Materials/power_law_creep/max_inelastic_increment=100 Materials/power_law_creep/maximum_number_substeps=2 Executioner/dt=100 Executioner/dtmin=100 Executioner/end_time=100 Outputs/file_base=ad_power_law_adaptive_substepping_out"
  expect_err = 'Adaptive substepping failed. Maximum number of substeps exceeded.'
  requirement = "The system shall converge under large deformation owing to the help of adaptive substepping using automatic differentiation with an increment based substepping criterion"
[]

[ad_power_law_adaptive_substepping_error]
  # We only check that the simulation converges with adaptive substepping
  type = 'RunException'
  input = 'ad_power_law_creep.i'
  cli_args = "Materials/power_law_creep/use_substepping=ERROR_BASED Materials/power_law_creep/adaptive_substepping=true Materials/power_law_creep/substep_strain_tolerance=100 Materials/power_law_creep/maximum_number_substeps=2 Executioner/dt=100 Executioner/dtmin=100 Executioner/end_time=100 Outputs/file_base=ad_power_law_adaptive_substepping_out"
  expect_err = 'Adaptive substepping failed. Maximum number of substeps exceeded.'
  requirement = "The system shall converge under large deformation owing to the help of adaptive substepping using automatic differentiation with an error based substepping criterion"
[]

[power_law]
  type = 'RunException'
  input = 'power_law_creep.i'
  expect_err = 'Exceeded maximum iterations in material return mapping iterations.'
  cli_args = 'Executioner/dtmin=0.1 Outputs/file_base=power_law_out'
  requirement = 'The system shall report diverged return mapping solve due to large deformation imposed in one step'
[]

[power_law_adaptive_substepping_increment]
  # We only check that the simulation converges with adaptive substepping
  type = 'RunException'
  input = 'power_law_creep.i'
  cli_args = "Materials/power_law_creep/use_substepping=INCREMENT_BASED Materials/power_law_creep/adaptive_substepping=true Materials/power_law_creep/max_inelastic_increment=100 Materials/power_law_creep/maximum_number_substeps=2 Executioner/dt=100 Executioner/dtmin=100 Executioner/end_time=100 Outputs/file_base=power_law_adaptive_substepping_out"
  expect_err = 'Adaptive substepping failed. Maximum number of substeps exceeded.'
  requirement = "The system shall converge under large deformation owing to the help of adaptive substepping with an increment based substepping criterion"
[]

[power_law_adaptive_substepping_error]
  # We only check that the simulation converges with adaptive substepping
  type = 'RunException'
  input = 'power_law_creep.i'
  cli_args = "Materials/power_law_creep/use_substepping=INCREMENT_BASED Materials/power_law_creep/adaptive_substepping=true Materials/power_law_creep/substep_strain_tolerance=100 Materials/power_law_creep/maximum_number_substeps=2 Executioner/dt=100 Executioner/dtmin=100 Executioner/end_time=100 Outputs/file_base=power_law_adaptive_substepping_out"
  expect_err = 'Adaptive substepping failed. Maximum number of substeps exceeded.'
  requirement = "The system shall converge under large deformation owing to the help of adaptive substepping with an error based substepping criterion"
[]

[./disp_about_axis_error1]
  design = DisplacementAboutAxis.md
  issues = '#15092'
  requirement = 'The system shall generate an error if the DisplacementAboutAxis boundary condition is applied to a component outside the range applicable for the dimensionality of the model'
  type = 'RunException'
  input = 'disp_about_axis_errors.i'
  cli_args = 'BCs/top_y/component=3 BCs/top_y/axis_direction="0. 0. 1."'
  expect_err = 'Invalid component given for top_y: 3.'
[../]

[./disp_about_axis_error2]
  design = DisplacementAboutAxis.md
  issues = '#15092'
  requirement = 'The system shall generate an error if the vector prescribed using the axis_direction parameter in the DisplacementAboutAxis boundary condition has a length of zero'
  type = 'RunException'
  input = 'disp_about_axis_errors.i'
  cli_args = 'BCs/top_y/component=1 BCs/top_y/axis_direction="0. 0. 0."'
  expect_err = 'Please specify a non-zero direction vector for the axis_direction in top_y'
[../]

[./action_error_1]
  type = 'RunException'
  input = 'truss_2d_action.i'
  cli_args = 'Outputs/exodus=false Outputs/csv=false'
  expect_err = 'LineElementAction: `area` should be provided as input for'
  requirement = "The LineElementAction shall produce an error if `area` is not provided as input for truss elements."
[../]

[./action_error_2]
  type = 'RunException'
  input = 'truss_2d_action.i'
  cli_args = 'Outputs/exodus=false Outputs/csv=false Physics/SolidMechanics/LineElement/QuasiStatic/block/area=area Physics/SolidMechanics/LineElement/QuasiStatic/block/rotations="rot_x"'
  expect_err = 'LineElementAction: Rotational variables cannot be set for truss elements.'
  requirement = "The LineElementAction shall produce an error if rotational variables are provided as input for truss elements."
[../]

[./error_nohardening]
  type = 'RunException'
  input = 'truss_plastic.i'
  cli_args = 'Outputs/exodus=false Outputs/csv=false'
  expect_err = 'PlasticTruss: Either hardening_constant or hardening_function must be defined'
  design = "source/materials/PlasticTruss.md"
  requirement = "PlasticTruss material should produce error if neither the hardening constant nor a hardening function is provided."
[../]

[./error_hardening]
  type = 'RunException'
  input = 'truss_plastic.i'
  cli_args = 'Outputs/exodus=false Outputs/csv=false Materials/truss/hardening_function=hf Materials/truss/hardening_constant=0.0'
  expect_err = 'PlasticTruss: Only the hardening_constant or only the hardening_function can be defined but not both'
  design = "source/materials/PlasticTruss.md"
  requirement = "PlasticTruss material should produce error if both hardening constant and hardening function are provided."
[../]

[info]
  type = 'RunApp'
  input = 'utility_functions.i'
  library_mode = 'DYNAMIC'
  valgrind = 'NONE'
  cli_args = 'UserObjects/uexternaldb/execute_on=INITIAL'
  expect_out = '\*\* Abaqus Info \*\*\nInline message 1 2 0.100000 0.200000 0.300000 ONE      TWO     !'
  detail = 'info level.'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[warning]
  type = 'RunApp'
  input = 'utility_functions.i'
  library_mode = 'DYNAMIC'
  valgrind = 'NONE'
  cli_args = "UserObjects/step_uo/step_start_times='0.0 0.1 0.2' UserObjects/uexternaldb/execute_on=TIMESTEP_END"
  expect_out = '\*\* Abaqus Warning \*\*\nInline message 1 2 0.100000 0.200000 0.300000 ONE      TWO     !'
  detail = 'warning level.'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[non_fatal]
  type = 'RunApp'
  input = 'utility_functions.i'
  library_mode = 'DYNAMIC'
  valgrind = 'NONE'
  cli_args = "UserObjects/step_uo/step_start_times='0.0 0.1 0.2 0.3' UserObjects/uexternaldb/execute_on=TIMESTEP_END"
  expect_out = '\*\* Abaqus Non-fatal Error \*\*\nInline message 1 2 0.100000 0.200000 0.300000 ONE      TWO     !'
  detail = 'non-fatal error level.'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[error]
  type = 'RunException'
  input = 'utility_functions.i'
  library_mode = 'DYNAMIC'
  valgrind = 'NONE'
  cli_args = "UserObjects/step_uo/step_start_times='0.0 0.1 0.2 0.3 0.4' UserObjects/uexternaldb/execute_on=TIMESTEP_END"
  expect_err = '\*\*\* ERROR \*\*\*\nInline message 1 2 0.100000 0.200000 0.300000 ONE      TWO     !'
  detail = 'fatal error level.'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[total_strain_error]
  type = 'RunException'
  input = 'linear_strain_hardening.i'
  library_mode = 'DYNAMIC'
  expect_err = 'AbaqusUMATStress \'constant\': Incremental strain quantities are not available'
  cli_args = 'Physics/SolidMechanics/QuasiStatic/all/strain=SMALL'
  valgrind = 'NONE'
  requirement = 'The system shall raise an error if a UMAT that requires incremental quantities is being used with a total strain formulation'
  # skip test if test is being run out-of-tree. Issue Ref: #25279
  installation_type = in_tree
[]

[except1]
  type = 'RunException'
  input = 'except1.i'
  expect_err = 'Weak-Plane-Shear friction angle must not be less than Weak-Plane-Shear dilation angle'
  requirement = 'The system shall not allow a weak plane shear friction angle less than the dilation angle.'
[]

[except2]
  type = 'RunException'
  input = 'except2.i'
  expect_err = 'ComputeMultiPlasticityStress: transverse_direction vector must not have zero length'
  requirement = 'The system shall not allow a zero length transverse direction vector'
[]

[except3]
  type = 'RunException'
  input = 'large_deform1.i'
  cli_args = 'Physics/SolidMechanics/QuasiStatic/all/decomposition_method=TaylorExpansion'
  expect_err = 'Cannot take square root of a number less than or equal to zero in the calculation of C1 for the Rashid approximation for the rotation tensor. This zero or negative number may occur when elements become heavily distorted.'
  issues = '#19067'
  design = '/ComputeFiniteStrain.md'
  requirement = 'The system shall cut the timestep through a mooseException when the loading conditions deform the element enough to produce a negative number under a square root term in the Rashid approximation for the rotation tensor'
[]

[except1]
  type = 'RunException'
  input = 'except1.i'
  expect_err = 'Weak plane tensile strength must not be negative'
  requirement = 'The system shall not allow a negative tensile strength value.'
[]