Solid Mechanics System Requirements Specification

This template follows INL template TEM-135, "IT System Requirements Specification".

note

This document serves as an addendum to Framework System Requirements Specification and captures information for Software Requirement Specification (SRS) specific to the Solid Mechanics module.

Framework System Requirements Specification

Introduction

System Purpose

The purpose of the MOOSE Solid Mechanics module is to provide the foundational capabilities for computing mechanical deformation of solids and structures. It computes strains and stresses, and solves for the displacement field that satisfies equilibrium. It is intended to both provide a basic set of capabilities and also be readily extensible by applications based on it to represent specialized material behavior.

System Scope

The MOOSE Solid Mechanics module provides an extensible set of capabilities for solving for mechanical deformation of solids and structures. It provides a set of C++ base classes that define interfaces for MOOSE Material objects that compute various mechanical behavior of materials at quadrature points, which include elastic properties, strains, eigenstrains, stresses, inelastic behavior, and damage. It also provides the needed Kernel classes to account for the contributions of the stress, inertia and damping in the solution for the displacement field that satisfies equilibrium. These models support one-, two-, and three-dimensional models of continuous materials, with variety of options for the treatment of the lower-dimensional models, including plane stress, plane strain, axisymmetry, and generalized plane strain. These models support both small- and finite-strain assumptions.

For modeling discrete interfaces between solid elements, this module provides a system similar to that for continuous materials that permits the definition of traction-separation laws for cohesive behavior. It also has similar support for lower-dimensional elements to represent structural elements such as beams and shells. As for the continuum models, sets of Material and Kernel classes are defined for these cases.

This module also provides a comprehensive set of boundary conditions relevant for mechanics modeling, such as pressure and traction boundary conditions. It also provides extensive postprocessing capabilities for computing quantities relevant to mechanics, such as fracture integrals.

In addition to defining the base classes that enable modeling arbitrary materials, this module also provides a set of specializations of those models for widely-used assumptions of material behavior. These include elasticity tensors defined in a variety of ways for isotropic and anisotropic materials, basic creep, plasticity and damage models, and models for eigenstrains due to thermal expansion.

System Overview

System Context

The Solid Mechanics module is command-line driven. Like MOOSE, this is typical for a high-performance software that is designed to run across several nodes of a cluster system. As such, all usage of the software is through any standard terminal program generally available on all supported operating systems. Similarly, for the purpose of interacting through the software, there is only a single user, "the user", which interacts with the software through the command-line. The Solid Mechanics module does not maintain any back-end database or interact with any system daemons. It is an executable, which may be launched from the command line and writes out various result files as it runs.

Figure 1: Usage of the Solid Mechanics module and other MOOSE-based applications.

System Functions

Since the Solid Mechanics module is a command-line driven application, all functionality provided in the software is operated through the use of standard UNIX command line flags and the extendable MOOSE input file. The Solid Mechanics module is completely extendable so individual design pages should be consulted for specific behaviors of each user-defined object.

User Characteristics

Like MOOSE, there are three kinds of users working on the Solid Mechanics module:

Solid Mechanics module Developers: These are the core developers of the Solid Mechanics module. They are responsible for following and enforcing the software development standards of the module, as well as designing, implementing, and maintaining the software.
Developers: A scientist or engineer that uses the Solid Mechanics module alongside MOOSE to build their own application. This user will typically have a background in modeling or simulation techniques (and perhaps numerical analysis) but may only have a limited skillset when it comes to code development using the C++ language. This is the primary focus group of the module. In many cases, these developers will be encouraged to contribute module-appropriate code back to the Solid Mechanics module, or to MOOSE itself.
Analysts: These are users that will run the code and perform analysis on the simulations they perform. These users may interact with developers of the system requesting new features and reporting bugs found and will typically make heavy use of the input file format.

Assumptions and Dependencies

The Solid Mechanics module is developed using MOOSE and can itself be based on various MOOSE modules, as such the SRS for the Solid Mechanics module is dependent upon the files listed at the beginning of this document. Any further assumptions or dependencies are outlined in the remainder of this section.

The Solid Mechanics module has no constraints on hardware and software beyond those of the MOOSE framework. The Solid Mechanics module provides access to a number of code objects that perform computations such as material behavior and boundary conditions. These objects each make their own physics-based assumptions, such as the units of the inputs and outputs. Those assumptions are described in the documentation for those individual objects.

References

ISO/IEC/IEEE 24765:2010(E). Systems and software engineering—Vocabulary. first edition, December 15 2010.

@manual{ISO-systems-software,
    author = "24765:2010(E), ISO/IEC/IEEE",
    title = "Systems and software engineering---Vocabulary",
    edition = "First",
    year = "2010",
    month = "December 15"
}

ASME NQA-1. ASME NQA-1-2008 with the NQA-1a-2009 addenda: Quality Assurance Requirements for Nuclear Facility Applications. first edition, August 31 2009.

@manual{ASME-NQA-1-2008,
    author = "NQA-1, ASME",
    title = "ASME NQA-1-2008 with the NQA-1a-2009 addenda: Quality Assurance Requirements for Nuclear Facility Applications",
    orginization = "American Socieity of Mechanical Engineers",
    year = "2009",
    month = "August 31",
    edition = "First"
}

Definitions and Acronyms

This section defines, or provides the definition of, all terms and acronyms required to properly understand this specification.

Definitions

Verification: (1) The process of: evaluating a system or component to determine whether the products of a given development phase satisfy the conditions imposed at the start of that phase. (2) Formal proof of program correctness (e.g., requirements, design, implementation reviews, system tests) (24765:2010(E), 2010).

Acronyms

Acronym	Description
INL	Idaho National Laboratory
LGPL	GNU Lesser General Public License
MOOSE	Multiphysics Object Oriented Simulation Environment
NQA-1	Nuclear Quality Assurance Level 1
POSIX	Portable Operating System Interface
SRS	Software Requirement Specification

System Requirements

In general, the following is required for MOOSE-based development:

A Portable Operating System Interface (POSIX) compliant Unix-like operating system. This includes any modern Linux-based operating system (e.g., Ubuntu, Fedora, Rocky, etc.), or a Macintosh machine running either of the last two MacOS releases.

Hardware	Information
CPU Architecture	x86_64, ARM (Apple Silicon)
Memory	8 GB (16 GBs for debug compilation)
Disk Space	30GB

Libraries	Version / Information
GCC	8.5.0 - 12.2.1
LLVM/Clang	10.0.1 - 16.0.6
Intel (ICC/ICX)	Not supported at this time
Python	3.7 - 3.11
Python Packages	packaging pyaml jinja2

Functional Requirements

solid_mechanics: 1D Axisymmetric
8.1.1The system shall support generalized plane strain with incremental strain for 1D meshes using the SolidMechanics/QuasiStatic Physics.
8.1.2The system shall support generalized plane strain with small strain for 1D meshes using the SolidMechanics/QuasiStatic Physics.
8.1.3The system shall support generalized plane strain with finite strain for 1D meshes using the SolidMechanics/QuasiStatic Physics.
8.1.4The ComputeAxisymmetric1DIncrementalStrain class shall compute the elastic stress for a 1D axisymmetric small incremental strain formulation under a combination of applied tensile displacement and thermal expansion loading using the SolidMechanics/QuasiStatic Physics.
8.1.5The ComputeAxisymmetric1DSmallStrain class shall compute the elastic stress for a 1D axisymmetric small total strain formulation under a combination of applied tensile displacement and thermal expansion loading using the SolidMechanics/QuasiStatic Physics.
8.1.6The ComputeAxisymmetric1DFiniteStrain class shall compute the elastic stress for a 1D axisymmetric incremental finite strain formulation under a combination of applied tensile displacement and thermal expansion loading using the SolidMechanics/QuasiStatic Physics.
8.1.7The ComputeAxisymmetric1DIncrementalStrain class shall, under generalized plane strain conditions, compute the elastic stress for a 1D axisymmetric small incremental strain formulation under a combination of applied tensile displacement and thermal expansion loading.
8.1.8The ComputeAxisymmetric1DSmallStrain class shall, under generalized plane strain conditions, compute the elastic stress for a 1D axisymmetric small total strain formulation under a combination of applied tensile displacement and thermal expansion loading.
8.1.9The ComputeAxisymmetric1DFiniteStrain class shall, under generalized plane strain conditions, compute the elastic stress for a 1D axisymmetric incremental finite strain formulation under a combination of applied tensile displacement and thermal expansion loading.
8.1.10The StressDivergenceRZTensors class shall generate an error if used with Problem/rz_coord_axis set to anything other than Y

solid_mechanics: 1D Spherical
8.2.1The ComputeRSphericalSmallStrain class, called through the QuasiStatic SolidMechanics Physics action, shall compute the total linearized solution for the displacement of a solid isotropic elastic sphere with a pressure applied to the outer surface using a 1D spherical symmetric formulation with total small strain assumptions. The RankTwoSphericalComponent class, called through the QuasiStatic SolidMechanics Physics, shall output the spherical radial and hoop stress.
8.2.2The ComputeRSphericalIncrementalStrain class, called through the QuasiStatic SolidMechanics Physics, shall find the linearized incremental strain displacement of a solid isotropic elastic sphere with a pressure applied to the outer surface using a 1D spherical symmetric formulation with incremental small strain assumptions. The RankTwoSphericalComponent class, called through the QuasiStatic SolidMechanics Physics, shall output the spherical radial and hoop stress.
8.2.3The ComputeRSphericalFiniteStrain class, called through the QuasiStatic SolidMechanics Physics action, shall find the finite incremental strain displacement of a thick walled hollow isotropic elastic sphere under an applied load using a 1D spherical symmetric fomulation with incremental finite strain assumptions. The RankTwoSphericalComponent class, called through the QuasiStatic SolidMechanics Physics, shall output the spherical radial stress.

solid_mechanics: 2D Different Planes
8.3.1The tensor mechanics strain calculators shall solve plane strain in the x-y plane for small strain
8.3.2The tensor mechanics strain calculators shall solve plane strain in the x-y plane for incremental strain
8.3.3The tensor mechanics strain calculators shall solve plane strain in the x-y plane for finite strain
8.3.4The tensor mechanics strain calculators shall solve plane strain in the x-z plane for small strain
8.3.5The tensor mechanics strain calculators shall solve plane strain in the x-z plane for incremental strain
8.3.6The tensor mechanics strain calculators shall solve plane strain in the x-z plane for finite strain
8.3.7The tensor mechanics strain calculators shall solve plane strain in the y-z plane for small strain
8.3.8The tensor mechanics strain calculators shall solve plane strain in the y-z plane for incremental strain
8.3.9The tensor mechanics strain calculators shall solve plane strain in the y-z plane for finite strain
8.3.10The tensor mechanics strain calculators shall solve generalized plane strain in the x-y plane for small strain
8.3.11The tensor mechanics strain calculators shall solve generalized plane strain in the x-y plane for incremental strain
8.3.12The tensor mechanics strain calculators shall solve generalized plane strain in the x-y plane for finite strain
8.3.13The tensor mechanics strain calculators shall solve generalized plane strain in the x-z plane for small strain
8.3.14The tensor mechanics strain calculators shall solve generalized plane strain in the x-z plane for incremental strain
8.3.15The tensor mechanics strain calculators shall solve generalized plane strain in the x-z plane for finite strain
8.3.16The tensor mechanics strain calculators shall solve generalized plane strain in the y-z plane for small strain
8.3.17The tensor mechanics strain calculators shall solve generalized plane strain in the y-z plane for incremental strain
8.3.18The tensor mechanics strain calculators shall solve generalized plane strain in the y-z plane for finite strain
8.3.19The Jacobian for plane strain in the x-y plane shall be correct
8.3.20The Jacobian for plane strain in the x-z plane shall be correct
8.3.21The Jacobian for plane strain in the y-z plane shall be correct
8.3.22The Jacobian for generalized plane strain in the x-y plane shall be correct
8.3.23The Jacobian for generalized plane strain in the x-z plane shall be correct
8.3.24The Jacobian for generalized plane strain in the y-z plane shall be correct

solid_mechanics: 2D Geometries
8.4.1The ComputePlaneSmallStrain class shall compute the elastic stress and strain for a planar square geometry under tension using a total small plane strain formulation.
8.4.2The ComputePlaneSmallStrain class shall compute the same result for elastic strain and stress when using the B-bar volumentric locking correction as computed without the volumetric locking correction for a planar geometry using a total small plane strain formulation.
8.4.3The ComputePlaneFiniteStrain class shall compute the elastic stress and strain for a planar square geometry under tension using a finite incremental plane strain formulation.
8.4.4The ComputePlaneFiniteStrain class shall compute the same result for elastic strain and stress when using the B-bar volumentric locking correction as computed without the volumetric locking correction for a planar geometry using a finite incremental plane strain formulation.
8.4.5The ComputeAxisymmetricRZSmallStrain class shall compute the mechanical response for a pressurized hollow cylinder with a small total axisymmetric strain formulation.
8.4.6The ComputeAxisymmetricRZIncrementalStrain class shall compute the mechanical response for a pressurized hollow cylinder with a small incremental axisymmetric strain formulation.
8.4.7The ComputeAxisymmetricRZFiniteStrain class shall compute the mechanical response for a pressurized hollow cylinder with a small incremental axisymmetric strain formulation.
8.4.8The QuasiStatic SolidMechanics Physics shall calculate the elastic stress and strain response for a 3D pressurized hollow cylinder with a large strain incremental strain formulation.
8.4.9The ComputeAxisymmetricRZFiniteStrain class shall compute the reaction forces on the top surface of a cylinder which is loaded axially in tension.
8.4.10The ComputeAxisymmetricRZFiniteStrain class shall compute the reaction forces on the top surface of a cylinder which is loaded axially in tension when using the B-bar volumetric locking correction.
8.4.11The volumetric locking correction option in ComputeAxisymmetricRZFiniteStrain shall reinit material properties without inverting a zero tensor when called from a side postprocessor applied to the axis of rotation in an axisymmetric simulation.

solid_mechanics: Cylindricalranktwoaux
8.5.1The Tensor Mechanics system shall support transformations of a rank two tensor into cylindircal coordinates.
8.5.2The Tensor Mechanics system including volumetric locking correction shall support transformations of a rank two tensor into cylindircal coordinates.

solid_mechanics: Accumulate Aux
8.6.1The system shall provide an aux kernel that accumulates the values of a given variable.

solid_mechanics: Action
8.7.1The QuasiStatic SolidMechanics Physics shall support changing the base name when creating a consistent strain calculator material and stress divergence kernel and shall generate different sets of outputs for different mesh subblocks with the appropriate base name.
8.7.2The QuasiStatic SolidMechanics Physics shall create a consistent strain calculator material and stress divergence kernel and shall generate different sets of outputs for different mesh subblocks.
8.7.3The QuasiStatic SolidMechanics Physics shall create different sets of consistent strain calculator material and stress divergence kernel pairs for different mesh subblocks requesting different strain formulations.
8.7.4The QuasiStatic SolidMechanics Physics shall create different sets of consistent strain calculator material and stress divergence kernel pairs for different mesh subblocks requesting different strain formulations with the new kernel system.
8.7.5The 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.
8.7.6The 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.
8.7.7The QuasiStatic SolidMechanics Physics shall warn if global Master action parameters are supplied but no Master action subblock have been added.
8.7.8The QuasiStatic SolidMechanics Physics shall create different sets of consistent strain calculator material and stress divergence kernel pairs for different mesh subblocks using different coordinate systems.
8.7.9The QuasiStatic SolidMechanics Physics shall error if an input file assigns the same QuasiStatic SolidMechanics Physics block to mesh blocks with different coordinate systems.
8.7.10The QuasiStatic SolidMechanics Physics shall extract eigenstrain names from material classes and correctly output these names to the console.
8.7.11The QuasiStatic SolidMechanics Physics shall extract eigenstrain names from material classes and correctly output these names to the console.
8.7.12The QuasiStatic SolidMechanics Physics shall extract eigenstrain names from material classes and correctly output these names to the console.
8.7.13The QuasiStatic SolidMechanics Physics shall extract eigenstrain names from material classes and correctly output these names to the console.
8.7.14The QuasiStatic SolidMechanics Physics shall extract eigenstrain names from material classes and correctly output these names to the console.
8.7.15The QuasiStatic SolidMechanics Physics shall extract eigenstrain names from material classes and correctly output these names to the console.
8.7.16The QuasiStatic SolidMechanics Physics shall determine the necessary orders and families to apply to material outputs.
8.7.17The QuasiStatic SolidMechanics Physics shall apply a single given output order to all output variables.
8.7.18The QuasiStatic SolidMechanics Physics shall apply a single given output family to all output variables.
8.7.19The QuasiStatic SolidMechanics Physics shall error if an the material outputs and or families do not match the number of material outputs.
8.7.20The QuasiStatic SolidMechanics Physics shall determine the necessary orders and families to apply to material outputs.
8.7.21The QuasiStatic SolidMechanics Physics shall determine the necessary orders and families to apply to material outputs.
8.7.22The QuasiStatic SolidMechanics Physics shall permit scalar quantity output from arbitrary tensors.

solid_mechanics: Ad 1D Spherical
8.8.1The ComputeRSphericalSmallStrain class, called through the QuasiStatic SolidMechanics Physics action, shall compute the total linearized solution for the displacement of a solid isotropic elastic sphere with a pressure applied to the outer surface using a 1D spherical symmetric formulation with total small strain assumptions. The RankTwoSphericalComponent class, called through the QuasiStatic SolidMechanics Physics, shall output the spherical radial and hoop stress.
8.8.2The ComputeRSphericalIncrementalStrain class, called through the QuasiStatic SolidMechanics Physics, shall find the linearized incremental strain displacement of a solid isotropic elastic sphere with a pressure applied to the outer surface using a 1D spherical symmetric formulation with incremental small strain assumptions. The RankTwoSphericalComponent class, called through the QuasiStatic SolidMechanics Physics, shall output the spherical radial and hoop stress.
8.8.3The ComputeRSphericalFiniteStrain class, called through the QuasiStatic SolidMechanics Physics action, shall find the finite incremental strain displacement of a thick walled hollow isotropic elastic sphere under an applied load using a 1D spherical symmetric fomulation with incremental finite strain assumptions. The RankTwoSphericalComponent class, called through the QuasiStatic SolidMechanics Physics, shall output the spherical radial stress.
8.8.4The Jacobian for the AD small strain elasticity problem with Pressure BC in spherical coordinates shall be perfect
8.8.5The Jacobian for the AD small incremental strain elasticity problem with Pressure BC in spherical coordinates shall be perfect
8.8.6The Jacobian for the AD small incremental strain elasticity problem with Pressure BC in spherical coordinates shall be perfect

solid_mechanics: Ad 2D Geometries
8.9.1The ADComputeAxisymmetricRZSmallStrain class shall compute the mechanical response for a pressurized hollow cylinder with a small total axisymmetric strain formulation.
8.9.2The ADComputeAxisymmetricRZIncrementalStrain class shall compute the mechanical response for a pressurized hollow cylinder with a small incremental axisymmetric strain formulation.
8.9.3The ADComputeAxisymmetricRZFiniteStrain class shall compute the mechanical response for a pressurized hollow cylinder with a small incremental axisymmetric strain formulation.
8.9.4The ADComputeAxisymmetricRZFiniteStrain class shall compute the reaction forces on the top surface of a cylinder which is loaded axially in tension.
8.9.5The ADComputeAxisymmetricRZFiniteStrain class shall compute the reaction forces on the top surface of a cylinder which is loaded axially in tension when using the B-bar volumetric locking correction.
8.9.6The volumetric locking correction option in ADComputeAxisymmetricRZFiniteStrain shall reinit material properties without inverting a zero tensor when called from a side postprocessor applied to the axis of rotation in an axisymmetric simulation.
8.9.7The ADComputeAxisymmetricRZSmallStrain class shall compute the mechanical response for a pressurized hollow cylinder with a small total axisymmetric strain formulation and shall produce perfect jacobians.
8.9.8The ADComputeAxisymmetricRZIncrementalStrain class shall compute the mechanical response for a pressurized hollow cylinder with a small incremental axisymmetric strain formulation and shall produce perfect jacobians.
8.9.9The ADComputeAxisymmetricRZFiniteStrain class shall compute the mechanical response for a pressurized hollow cylinder with a small incremental axisymmetric strain formulation and shall produce perfect jacobians.
8.9.10The ADComputeAxisymmetricRZFiniteStrain class shall compute the reaction forces on the top surface of a cylinder which is loaded axially in tension and shall produce perfect jacobians.
8.9.11The ADComputeAxisymmetricRZFiniteStrain class shall compute the reaction forces on the top surface of a cylinder which is loaded axially in tension when using the B-bar volumetric locking correction and shall produce perfect jacobians.
8.9.12The volumetric locking correction option in ADComputeAxisymmetricRZFiniteStrain shall reinit material properties without inverting a zero tensor when called from a side postprocessor applied to the axis of rotation in an axisymmetric simulation and shall produce perfect jacobians.

solid_mechanics: Ad Action
8.10.1The QuasiStatic SolidMechanics Physics shall create a consistent strain calculator material and stress divergence kernel and shall generate different sets of outputs for different mesh subblocks.
8.10.2The QuasiStatic SolidMechanics Physics shall create different sets of consistent strain calculator material and stress divergence kernel pairs for different mesh subblocks requesting different strain formulations.
8.10.3The 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.
8.10.4The 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.
8.10.5The QuasiStatic SolidMechanics Physics shall create different sets of consistent strain calculator material and stress divergence kernel pairs for different mesh subblocks using different coordinate systems.
8.10.6The QuasiStatic SolidMechanics Physics shall error if an input file assigns the same QuasiStatic SolidMechanics Physics block to mesh blocks with different coordinate systems.
8.10.7The Jacobian for the automatic differentiation in the two_block testproblem shall be perfect
8.10.8The Jacobian for the automatic differentiation in the two_block testproblem shall be perfect (non action test case)
8.10.9The Jacobian for the automatic differentiation in the two_block_new problem shall be perfect
8.10.10The Jacobian for the automatic differentiation two_coord problem shall be perfect

solid_mechanics: Ad Anisotropic Creep
8.11.1The system shall avoid regression on material time step and combined anisotropic creep computations
8.11.2The system shall be capable of reproducing verification results of creep strain ratios along the X direction.
8.11.3The system shall be capable of reproducing verification results of creep strain ratios along the X direction using hand-coded Jacobians.
8.11.4The system shall reproduce the creep strain as those evaluated with residual and derivative computations used for isotropic elasticity.
8.11.5The system shall compute shear creep strain in an element subject to shear load in one direction.
8.11.6The system shall compute shear creep strain in an element subject to shear load in two directions.
8.11.7The system shall compute creep strains in an element with anisotropic elasticity and subject to biaxial tensile load.
8.11.8The system shall be capable of reproducing verification results of creep strain ratios along the X direction using hand-coded Jacobians when no transformation is used to rotate the Hill tensor for small deformation.
8.11.9The system shall be capable of reproducing verification results of creep strain ratios along the Y direction.
8.11.10The system shall be capable of reproducing verification results of creep strain ratios along the Z direction.
8.11.11The system shall be capable of reproducing isotropic creep with the right anisotropic creep parameters: Baseline
8.11.12The system shall provide a perfect Jacobian while calculating large deformation creep.
8.11.13The system shall provide a reasonable Jacobian while calculating large deformation creep using the elasticity tensor as the tangent operator.
8.11.14The system shall be capable of reproducing isotropic creep with the right anisotropic creep parameters.
8.11.15The system shall be capable of enforcing a time step such that the creep rate integration error is controled by the user with the aid of a soft terminator when regular (non-AD) objects are used.
8.11.16The system shall be capable of enforcing a time step such that the creep rate integration error is controled by the user with the aid of a soft terminator.
8.11.17The system shall be capable of running Hill materials with constant coefficients and coefficients that are function of temperature to capture material texture. This test performs a reference run with constant coefficients
8.11.18The system shall be capable of running Hill materials with constant coefficients and coefficients that are function of temperature to capture material texture. This test performs a run using a constant temperature function.
8.11.19The system shall be capable of running anisotropic Hill creep with temperature-dependent Hill coefficients and generating correct output.
8.11.20The system shall be capable of using large deformation rotation to update the orientation of Hill parameters even when there is no rigid body rotation. This test also serves as a reference to verify the updates to the Hill parameters when there is significant rigid body motion
8.11.21The system shall be capable of using large deformation rotation to update the orientation of Hill parameters when there is 90-degree rigid body orientation change about the Z axis. The creep strains have to be very similar to the output when there is no rigid body motion (Note: The mere constraints enforcing rigid body rotation induce numerical errors). In this case using Hill constants.
8.11.22The system shall be capable of using large deformation rotation to update the orientation of Hill parameters even when there is no rigid body rotation. This test also serves as a reference to verify the updates to the Hill parameters when there is significant rigid body motion
8.11.23The system shall be capable of not using large deformation rotation to define Hill anisotropic parameters.
8.11.24The system shall be capable of using large deformation rotation to update the orientation of Hill parameters when there is 90-degree rigid body orientation change about the X axis. The creep strains have to be very similar to the output when there is no rigid body motion (Note: The mere constraints enforcing rigid body rotation induce numerical errors). In this case using Hill constants.

solid_mechanics: Ad Elastic
8.12.1We shall be able to reproduce finite strain elasticity results of the hand-coded simulation using automatic differentiation. (non-AD reference)
8.12.2We shall be able to reproduce finite strain elasticity results of the hand-coded simulation using automatic differentiation.
8.12.3The Jacobian for the AD finite strain elasticity problem shall be perfect
8.12.4We shall be able to reproduce incremental small strain elasticity results of the hand-coded simulation using automatic differentiation. (non-AD reference)
8.12.5We shall be able to reproduce incremental small strain elasticity results of the hand-coded simulation using automatic differentiation.
8.12.6The Jacobian for the AD incremental small strain elasticity problem shall be perfect
8.12.7The system shall provide an AD enabled Green-Lagrange strain calculator
8.12.8The Jacobian for the Green-Lagrange strain calculator shall be perfect
8.12.9We shall be able to reproduce finite strain elasticity results of the hand-coded simulation in cylindrical coordinates using automatic differentiation. (non-AD reference)
8.12.10We shall be able to reproduce finite strain elasticity results of the hand-coded simulation in cylindrical coordinates using automatic differentiation.
8.12.11The Jacobian for the AD finite strain elasticity problem in cylindrical coordinates shall be perfect
8.12.12We shall be able to reproduce incremental small strain elasticity results of the hand-coded simulation in cylindrical coordinates using automatic differentiation. (non-AD reference)
8.12.13We shall be able to reproduce incremental small strain elasticity results of the hand-coded simulation in cylindrical coordinates using automatic differentiation.
8.12.14The Jacobian for the AD incremental small strain elasticity problem in cylindrical coordinates shall be perfect
8.12.15We shall be able to reproduce small strain elasticity results of the hand-coded simulation in cylindrical coordinates using automatic differentiation. (non-AD reference)
8.12.16We shall be able to reproduce small strain elasticity results of the hand-coded simulation in cylindrical coordinates using automatic differentiation.
8.12.17The Jacobian for the AD small strain elasticity problem in cylindrical coordinates shall be perfect
8.12.18We shall be able to reproduce finite strain elasticity results of the hand-coded simulation in spherical coordinates using automatic differentiation. (non-AD reference)
8.12.19We shall be able to reproduce finite strain elasticity results of the hand-coded simulation in spherical coordinates using automatic differentiation.
8.12.20The Jacobian for the AD finite strain elasticity problem in spherical coordinates shall be perfect
8.12.21We shall be able to reproduce incremental small strain elasticity results of the hand-coded simulation in spherical coordinates using automatic differentiation. (non-AD reference)
8.12.22We shall be able to reproduce incremental small strain elasticity results of the hand-coded simulation in spherical coordinates using automatic differentiation.
8.12.23The Jacobian for the AD incremental small strain elasticity in spherical coordinates problem shall be perfect
8.12.24We shall be able to reproduce small strain elasticity results of the hand-coded simulation in spherical coordinates using automatic differentiation. (non-AD reference)
8.12.25We shall be able to reproduce small strain elasticity results of the hand-coded simulation in spherical coordinates using automatic differentiation.
8.12.26The Jacobian for the AD small strain elasticity problem in spherical coordinates shall be perfect

solid_mechanics: Ad Finite Strain Jacobian
8.13.1Finite strain methods in Tensor Mechanics should be able to adequately simulate a bar bending simulation in 2D using AD and match non-AD methods
8.13.2Finite strain methods in Tensor Mechanics should be able to adequately simulate a bar bending simulation in 2D using a volumetric locking correction using AD and match non-AD methods
8.13.3Finite strain methods in Tensor Mechanics should be able to adequately simulate a tensile test simulation in 3D using AD and match non-AD methods
8.13.4Finite strain methods in Tensor Mechanics should be able to adequately simulate a tensile test simulation in 3D using a volumetric locking correction using AD and match non-AD methods
8.13.5Finite strain methods in Tensor Mechanics should be able to adequately simulate a bar bending simulation in 2D using AD and calculate perfect Jacobians
8.13.6Finite strain methods in Tensor Mechanics should be able to adequately simulate a bar bending simulation in 2D using a volumetric locking correction using AD and calculate perfect Jacobians
8.13.7Finite strain methods in Tensor Mechanics should be able to adequately simulate a tensile test simulation in 3D using AD and calculate perfect Jacobians
8.13.8Finite strain methods in Tensor Mechanics should be able to adequately simulate a tensile test simulation in 3D using a volumetric locking correction using AD and calculate perfect Jacobians
8.13.9Finite 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

solid_mechanics: Ad Isotropic Elasticity Tensor
8.14.1The ComputeIsotropicElasticityTensor class shall correctly compute the elasticity tensor from the lambda and shear modulus for an isotropic material using AD formulations.
8.14.2The ComputeIsotropicElasticityTensor class shall correctly compute the elasticity tensor from the Young's modulus and Poisson's ratio for an isotropic material using AD formulations.
8.14.3The ComputeIsotropicElasticityTensor class shall correctly compute the elasticity tensor from their bulk modulus and shear modulus for an isotropic material using AD formulations.
8.14.4The ComputeElasticityTensor class shall correctly compute the elasticity tensor for an isotropic axisymmetric problem.
8.14.5The ComputeIsotropicElasticityTensor class shall correctly compute the elasticity tensor from the lambda and shear modulus for an isotropic material using AD formulations and produce a perfect Jacobian.
8.14.6The ComputeIsotropicElasticityTensor class shall correctly compute the elasticity tensor from the Young's modulus and Poisson's ratio for an isotropic material using AD formulations and produce a perfect Jacobian.
8.14.7The ComputeIsotropicElasticityTensor class shall correctly compute the elasticity tensor from their bulk modulus and shear modulus for an isotropic material using AD formulations and produce a perfect Jacobian.
8.14.8The ComputeElasticityTensor class shall correctly compute the elasticity tensor for an isotropic axisymmetric problem and produce a perfect Jacobian.

solid_mechanics: Ad Linear Elasticity
8.15.1We shall be able to reproduce linear elastic stress results of the hand-coded simulation using automatic differentiation.
8.15.2The Jacobian for the AD linear elastic stress problem shall be perfect
8.15.3We shall be able to introduce extra stresses into the stress calculators using automatic differentiation.
8.15.4The Jacobian for the AD linear elastic stress problem shall be perfect
8.15.5We shall be able to reproduce eigenstrain results of the hand-coded simulation using automatic differentiation.
8.15.6The Jacobian for the AD eigenstrain problem shall be perfect
8.15.7We shall be able to reproduce small strain with specified tensors results of the hand-coded simulation using automatic differentiation.
8.15.8The Jacobian for the AD small strain with specified tensors problem shall be perfect
8.15.9We shall be able to reproduce thermal eigenstrain results of the hand-coded simulation using automatic differentiation.
8.15.10The Jacobian for the AD thermal eigenstrain problem shall be perfect

solid_mechanics: Ad Plastic
8.16.1The AD multiple inelastic stress calculator shall provide a correct stress for a single power law creep model (reference computation)
8.16.2The AD multiple inelastic stress calculator shall provide a correct stress for a single power law creep model and an additional zero creep power law model
8.16.3The AD multiple inelastic stress calculator shall provide a correct stress for the linear combination of two power law creep models
8.16.4The AD multiple inelastic stress calculator shall provide a correct stress when cycling through two identical power law creep models
8.16.5The AD multiple inelastic stress calculator shall provide a correct jacobian for a single power law creep model
8.16.6The AD multiple inelastic stress calculator shall provide a correct jacobian for a single power law creep model and an additional zero creep power law model
8.16.7The AD multiple inelastic stress calculator shall provide a correct jacobian for the linear combination of two power law creep models
8.16.8The AD multiple inelastic stress calculator shall provide a correct jacobian when cycling through two identical power law creep models

solid_mechanics: Ad Pressure
8.17.1The Pressure boundary condition action shall create the objects needed to apply automatic differentiation pressure boundary conditions on a 3D model as demonstrated by correctly computing the response of an elastic small-strain isotropic unit cube with pressure applied on three faces to create a hydrostatic pressure and match non-AD methods.
8.17.2The Pressure boundary condition action shall create the objects needed to apply automatic differentiation pressure boundary conditions on a 3D model as demonstrated by correctly computing the response of an elastic small-strain isotropic unit cube with pressure applied on three faces to create a hydrostatic pressure using the volumetric locking correction b-bar formulation and match non-AD methods.
8.17.3The Pressure boundary condition action shall create the objects needed to apply automatic differentiation pressure boundary conditions on a 3D model as demonstrated by correctly computing the response of an elastic small-strain isotropic unit cube with pressure applied on three faces to create a hydrostatic pressure and calculate a perfect Jacobian.
8.17.4The Pressure boundary condition action shall create the objects needed to apply automatic differentiation pressure boundary conditions on a 3D model as demonstrated by correctly computing the response of an elastic small-strain isotropic unit cube with pressure applied on three faces to create a hydrostatic pressure using the volumetric locking correction b-bar formulation and calculate a perfect Jacobian.

solid_mechanics: Ad Return Mapping
8.18.1The return mapping algorithm shall use automatic differentiation to compute the derivative of the yield function with respect to the internal variable, and the solution should be the same as existing hand coded derivative.

solid_mechanics: Ad Simple Linear
8.19.1We shall be able to run a simple linear small-strain problem using a hand-coded Jacobian
8.19.2We shall be able to reproduce the results of the hand-coded simulation using automatic differentiation in the production stress divergence kernel
8.19.3We shall be able to reproduce the results of the hand-coded simulation using automatic differentiation with reversed stress and strain materials
8.19.4The Jacobian for the hand-coded problem shall be perfect
8.19.5The Jacobian for the automatic differentiation problem shall be perfect
8.19.6The Jacobian for the automatic differentiation problem with reversed stress and strain materials shall be perfect

solid_mechanics: Ad Smeared Cracking
8.20.1The MOOSE SolidMechanics module shall simulate cracking on a specimen under tension in cartesian coordinates using AD and match non-AD methods.
8.20.2The MOOSE SolidMechanics module shall simulate cracking on a specimen under tension in rz coordinates using AD and match non-AD methods.
8.20.3The MOOSE SolidMechanics module shall simulate cracking while the cracking strength is prescribed by an elemental AuxVariable using AD and match non-AD methods.
8.20.4The MOOSE SolidMechanics module shall simulate exponential stress release using AD and match non-AD methods.
8.20.5The MOOSE SolidMechanics module shall simulate exponential stress relase while using the rz coordinate system using AD and match non-AD methods.
8.20.6The MOOSE SolidMechanics module shall demonstrate softening using the power law for smeared cracking using AD and match non-AD methods.
8.20.7The MOOSE SolidMechanics module shall demonstrate the prescribed softening laws in three directions, power law (x), exponential (y), and abrupt (z) using AD and match non-AD methods.
8.20.8The MOOSE SolidMechanics module shall simulate smeared cracking in the x y and z directions using AD and match non-AD methods.
8.20.9The MOOSE SolidMechanics module shall simulate smeared cracking under plane stress conditions using AD and match non-AD methods.
8.20.10The MOOSE SolidMechanics module shall demonstrate that the smeared cracking model correctly handles finite rotation of cracked elements using AD and match non-AD methods.
8.20.11The MOOSE SolidMechanics module shall demonstrate the finite rotation of cracked elements where the crack is prescribed in x using AD and match non-AD methods.
8.20.12The MOOSE SolidMechanics module shall demonstrate the finite rotation of cracked elements where the crack is prescribed in z using AD and match non-AD methods.
8.20.13The MOOSE SolidMechanics module shall demonstrate the finite rotation of cracked elements where two cracks are prescribed in x and z using AD and match non-AD methods.
8.20.14The MOOSE SolidMechanics module shall compute accurate AD Jacobian of system with multiple softening laws.
8.20.15The ComputeSmearedCrackingStress model shall generate a warning if 'cracked_elasticity_type' is not specified
8.20.16The ComputeSmearedCrackingStress model shall generate an error if an invalid number of softening models is specified

solid_mechanics: Ad Thermal Expansion Function
8.21.1
The system shall compute an eigenstrain due to thermal expansion using a function that describes a constant mean and instantaneous thermal expansion using the AD formulation
1. and the finite strain formulation
2. and the small strain formulation
8.21.2
The system shall compute an eigenstrain due to thermal expansion using a function that describes a mean and instantaneous thermal expansion with a linear relationship to temperature using the AD formulation
1. and the finite strain formulation
2. and the small strain formulation
8.21.3
The system shall compute an eigenstrain due and allow a smooth transition from negative to positive strain across the reference temperature and compare favorably to hand calculations
1. using a mean thermal expansion coefficient
2. using a complex mean thermal expansion coefficient
3. using a complex mean thermal expansion coefficient and compute a perfect jacobian
4. using a instantaneous thermal expansion coefficient
5. using a complex instantaneous thermal expansion coefficient
6. using a complex instantaneous thermal expansion coefficient and compute a perfect jacobian
7. using a dilatation thermal expansion coefficient

solid_mechanics: Ad Viscoplasticity Stress Update
8.22.1The PowerLawCreepStressUpdate, called through the ADComputeMultipleInelasticStress, shall compute a creep strain based on an extrenal loading.
8.22.2The Jacobian for the AD regular creep problem shall be perfect
8.22.3The ADViscoplasticityStressUpdate class shall compute a ratio between the gauge stress, equilvalent stress, and hydrostatic stress across a wide swath of exponents and stress states using spherical pore geometry.
8.22.4The Jacobian for the AD exact spherical problem shall be perfect
8.22.5The ADViscoplasticityStressUpdate class shall compute a ratio between the gauge stress, equilvalent stress, and hydrostatic stress across a wide swath of exponents and stress states using spherical pore geometry.
8.22.6The Jacobian for the AD exact cylindrical problem shall be perfect
8.22.7The ADViscoplasticityStressUpdate class shall compute the viscoplastic response using a single model with LPS spherical formulation that increases the porosity due to an external strain.
8.22.8The Jacobian for the AD lps single problem shall be perfect
8.22.9The ADViscoplasticityStressUpdate class shall compute the viscoplastic response using two LPS models with spherical formulations and the same stress exponential that is close to combining the models into a single ADViscoplasticityStressUpdate instance.
8.22.10The Jacobian for the AD lps single split problem shall be perfect
8.22.11The ADViscoplasticityStressUpdate class shall compute the viscoplastic response using two LPS models with spherical formulations and two different stress exponents that increases the porosity due to an external strain.
8.22.12The Jacobian for the AD lps dual problem shall be perfect
8.22.13The ADViscoplasticityStressUpdate class shall compute the viscoplastic response using a single model with GTN formulation that increases the porosity due to an external strain.
8.22.14The Jacobian for the AD gtn single problem shall be perfect
8.22.15
The system shall handle a negative calculated porosity computed from a strain tensor
1. by setting the porosity to zero.
2. by setting the porosity to the initial condition.
3. by throwing an exception.

solid_mechanics: Anisotropic Elastoplasticity
8.23.1Anisotropic plasticity and anisotropic elastoplasticity must give almost identical results if elastic behavior is isotropic – plastic anisotropy
8.23.2Anisotropic plasticity and anisotropic elastoplasticity must give almost identical results if elastic behavior is isotropic – elastoplastic anisotropy
8.23.3The system shall be able to solve anisotropic plasticity and anisotropic elastoplasticity problems with elastoplastic anisotropy
8.23.4Uniaxial tensile test run in x direction with anisotropic elasticity and linear strain hardening with automatic differentiation
8.23.5Uniaxial tensile test run in y direction with anisotropic elasticity and linear strain hardening with automatic differentiation
8.23.6Uniaxial tensile test run in x direction with anisotropic elasticity and non-linear strain hardening with automatic differentiation
8.23.7Compute hoop plastic strain in a cylinder and compare it for two elements located separately in hoop direction
8.23.8Compute hoop plastic strain in a cylinder oriented along the x axis and compare it for two elements located separately in hoop direction for coarsely meshed model
8.23.9Compute hoop plastic strain in a cylinder oriented along the y axis and compare it for two elements located separately in hoop direction for coarsely meshed model
8.23.10Compute hoop plastic strain in a cylinder oriented along the z axis and compare it for two elements located separately in hoop direction for coarsely meshed model

solid_mechanics: Anisotropic Patch
8.24.1The mechanics system shall be capable of accurately computing the elastic response of an anisotropic elastic material where 6 components of a symmetric elasticity tensor are output on an irregular patch of elements with total small strain assumptions

solid_mechanics: Anisotropic Plasticity
8.25.1Anisotropic plasticity must yield same results as finite strain elasticity if yield condition is never met
8.25.2Anisotropic plasticity must yield same results as finite strain elasticity if yield condition is never met –elasticity
8.25.3Anisotropic plasticity with isotropic elasticity must reproduce simplified verification problem along x axis
8.25.4Anisotropic plasticity with hill tensor rotation must produce simplified verification problem along x axis
8.25.5Anisotropic plasticity with isotropic elasticity must reproduce simplified verification problem along y axis

solid_mechanics: Auxkernels
8.26.1The system shall compute the VonMises value of a RankTwoTensor
8.26.2The system shall allow RankTwoScalarAux to output principal stresses
8.26.3The system shall compute the local elastic energy

solid_mechanics: Beam
8.27.1The LineElementAction class shall correctly create the objects required for a mechanics simulation using beam or truss elements.
8.27.2The LineElementAction class shall correctly set the common parameters in the action subblocks.
8.27.3The LineElementAction class shall produce an error when the displacement variables are not provided by the user.
8.27.4The LineElementAction class shall produce an error if the user provided inputs for strain_type, rotation_type and use_displaced_mesh parameters are not compatible.
8.27.5The LineElementAction class shall produce an error if the number of variables listed in the save_in parameter differs from the number of displacement variables.
8.27.6The 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.
8.27.7The LineElementAction class shall produce an error if the names for the rotational degrees of freedom are not provided by the user.
8.27.8The LineElementAction class shall produce an error if the number of rotational variables provided as input differs from the number of displacement variables.
8.27.9The LineElementAction class shall produce an error if the moment of inertia, area and orientation of the beam are not provided as input.
8.27.10The LineElementAction class shall produce an error if translational and rotational velocities and accelerations are not provided as input for dynamic simulations using beam elements.
8.27.11The LineElementAction class shall produce an error if the number of translational and rotational velocities and accelerations differs from the number of displacement variables.
8.27.12The 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.
8.27.13The 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.
8.27.14The 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.
8.27.15The 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.
8.27.16The LineElementAction class shall produce an error if multiple subblocks specify properties for the same mesh block.
8.27.17The LineElementAction class shall produce an error if an action subblock is mesh block restricted while another is not.
8.27.18The 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.
8.27.19The mechanics system shall accurately predict the displacement of a beam element with a frictionless contact constraint.
8.27.20The mechanics system shall accurately predict the displacement of a beam element with a glued contact constraint.
8.27.21The mechanics system shall accurately predict the displacement of a beam element with a frictional contact constraint.
8.27.22The mechanics system shall correctly predict the natural frequencies of an Euler-Bernoulli beam modeled using beam elements with consistent mass/inertia.
8.27.23The mechanics system shall correctly predict the natural frequencies of a Timoshenko beam modeled using beam elements with consistent mass/inertia.
8.27.24The mechanics system shall correctly predict the natural frequencies of an Euler-Bernoulli beam modeled using beam elements in the presence of Rayleigh damping and numerical damping introduced by Hilber-Hughes-Taylor (HHT) time integration.
8.27.25The mechanics system shall correctly predict the natural frequencies of an Euler-Bernoulli beam modeled using beam elements in the presence of Rayleigh damping and numerical damping introduced by Hilber-Hughes-Taylor (HHT) time integration when using the velocity and acceleration computed using the Newmark-Beta time integrator.
8.27.26The mechanics system shall correctly predict the natural frequencies of a massless Euler-Bernoulli beam modeled using beam elements with a nodal masses placed at the ends.
8.27.27The mechanics system shall correctly predict the natural frequencies of a massless Euler-Bernoulli beam modeled using beam elements with added nodal masses when the location and values of the masses are provided using a csv file.
8.27.28The mechanics system shall correctly model the response of a beam modeled using beam elements when gravitational force (proportional to nodal mass) is applied to the beam.
8.27.29The mechanics system shall correctly model the response of a beam modeled using beam elements under gravitational force when the nodal mass distribution is provided using a csv file.
8.27.30The LineElementAction shall create the translational and rotational velocities and accelerations required for a dynamic simulation using beam elements.
8.27.31The mechanics system shall correctly model the response of a beam modeled using beam elements in the presence of nodal mass, nodal inertia and Rayleigh damping.
8.27.32The mechanics system shall correctly model the response of a beam modeled using beam elements in the presence of nodal mass, nodal inertia and Rayleigh damping when using the velocity and accelerations computed by the Newmark-Beta time integrator.
8.27.33The LineElementAction shall correctly create the input blocks required for a dynamic beam simulation using beam elements and a consistent mass/inertia matrix in the presence of Rayleigh damping and numerical damping in the form of Hilber-Hughes-Taylor (HHT) time integration.
8.27.34The LineElmentAction shall correctly create the input blocks required for a dynamic beam simulation using beam elements and nodal mass/inertia matrix in the presence of Rayleigh damping and numerical damping in the form of Hilber-Hughes-Taylor (HHT) time integration.
8.27.35The mechanics system shall correctly predict the natural frequency of a cantilever beam modeled using beam elements with a mass at the free end.
8.27.36The InertialForceBeam class shall produce an error if the number of variables provided for rotations differs from that provided for displacements.
8.27.37The NodalRotatioanlInertia class shall produce an error if the number of rotational velocities and accelerations provided as input differ from the number of rotations.
8.27.38The NodalRotationalInertia class shall produce an error if the user provided nodal inertia is not positive definite.
8.27.39The NodalRotatioanlInertia class shall produce an error if the user provided x and y orientations are not unit vectors.
8.27.40The NodalRotatioanlInertia class shall produce an error if the user provided x and y orientations are not perpendicular to each other.
8.27.41The NodalRotatioanlInertia class shall produce an error if only x or y orientation is provided as input by the user.
8.27.42The 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.
8.27.43The NodalTranslationalInertia class shall produce an error if nodal mass is provided as input both as a constant value and also using a csv file.
8.27.44The NodalTranslationalInertia class shall produce an error if nodal mass is not provided as input either as a constant value or using a csv file.
8.27.45The NodalTranslationalInertia class shall produce an error if the number of columns in the nodal mass file is not 4.
8.27.46The 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.
8.27.47The NodalGravity class shall produce an error if nodal mass is provided as input both as a constant value and also using a csv file.
8.27.48The NodalGravity class shall produce an error if nodal mass is not provided as input either as a constant value or using a csv file.
8.27.49The NodalGravity class shall produce an error if the number of columns in the nodal mass file is not 4.
8.27.50The 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.
8.27.51The 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.
8.27.52The NodalTranslationalInertia class shall produce an error if nodal mass is provided as input both as constant value and also using a csv file.
8.27.53The ComputeThermalExpansionEigenstrainBeam class shall correctly calculate eigenstrains due to changes in temperature.
8.27.54The ComputeEigenstrainBeamFromVariable class shall correctly transfer eigenstrains from auxvariables into eigenstrain material property.
8.27.55The ComputeEigenstrainBeamFromVariable class shall report an error if less than 3 displacement or rotational eigenstrains are provided by the user.
8.27.56The mechanics system shall accurately predict the static bending response of a Timoshenko beam modeled using beam elements under small deformations in the y direction.
8.27.57The mechanics system shall accurately predict the static bending response of a Timoshenko beam modeled using beam elements under small deformations in the z direction.
8.27.58The mechanics system shall accurately predict the static bending response of a Euler-Bernoulli beam modeled using beam elements under small deformations in the y direction.
8.27.59The mechanics system shall accurately predict the static bending response of a Euler-Bernoulli beam modeled using beam elements under small deformations in the z direction.
8.27.60The mechanics system shall accurately predict the static bending response of a Euler-Bernoulli beam modeled using beam elements under finite deformations in the y direction.
8.27.61The mechanics system shall accurately predict the static bending response of a Euler-Bernoulli beam modeled using beam elements under finite deformations in the z direction.
8.27.62The LineElementAction class shall accurately create the objects required to model the static bending response of an Euler-Bernoulli beam modeled using beam elements under small deformations.
8.27.63The LineElementAction class shall accurately create the objects required to model the static bending response of an Euler-Bernoulli beam modeled using beam elements under finite deformations.
8.27.64The mechanics system shall accurately predict the axial displacement of an Euler-Bernoulli pipe modeled using beam elements.
8.27.65The mechanics system shall accurately predict the axial forces on an Euler-Bernoulli pipe modeled using beam elements.
8.27.66The mechanics system shall accurately predict the bending response of an Euler-Bernoulli pipe modeled using beam elements.
8.27.67The ComputeIncrementalBeamStrain class shall produce an error if the number of supplied displacements and rotations do not match.
8.27.68The StressDivergenceBeam class shall produce an error if the number of supplied displacements and rotations do not match.
8.27.69The 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.
8.27.70The ComputeIncrementalBeamStrain class shall produce an error if the y orientation provided is not perpendicular to the beam axis.
8.27.71The mechanics system shall accurately predict the torsional response of a beam modeled using beam elements with auto-calculated polar moment of inertia.
8.27.72The mechanics system shall accurately predict the torsional response of a beam modeled using beam elements with user provided polar moment of inertia.
8.27.73
The mechanics system shall accurately predict the static bending response of an Euler beam modeled using beam elements under small deformation when the beam is
1. oriented along the global Z axis.
2. oriented on the YZ plane at a 45 deg. angle.
3. oriented on the YZ plane at a 45 deg. angle and has in-plane loading.
4. oriented on the YZ plane at a 45 deg. angle and has in-plane loading with non-symmetric cross section geometry.
5. oriented on the YZ plane at a 45 deg. angle and has in-plane loading and the cross section geometry is non-symmetric.
6. oriented along the global Y axis.
7. oriented on the XZ plane at a 45 deg. angle.
8. oriented on the XZ plane at a 45 deg. angle, and the external loading takes place on the same plane.
9. oriented on the XY plane at a 45 deg. angle.
10. oriented on the XY plane at a 45 deg. angle, and the external loading takes place on the same plane.
8.27.74
The mechanics system shall accurately predict the static bending response of an Euler beam modeled using beam elements under small deformations when the beam is
1. subjected to simply supported BCs and distributed loading.
2. subjected to combined bending and torsion loading.

solid_mechanics: Capped Drucker Prager
8.28.1
8.28.2
8.28.3
8.28.4
8.28.5
8.28.6
8.28.7
8.28.8
8.28.9
8.28.10
8.28.11
8.28.12

solid_mechanics: Capped Mohr Coulomb
8.29.1
8.29.2
8.29.3
8.29.4
8.29.5
8.29.6
8.29.7
8.29.8
8.29.9
8.29.10
8.29.11
8.29.12
8.29.13
8.29.14
8.29.15
8.29.16
8.29.17
8.29.18
8.29.19
8.29.20
8.29.21
8.29.22
8.29.23
8.29.24
8.29.25
8.29.26
8.29.27
8.29.28
8.29.29
8.29.30
8.29.31
8.29.32
8.29.33
8.29.34
8.29.35

solid_mechanics: Capped Weak Plane
8.30.1The CappedWeakPlaneStressUpdate model shall generate an error if the friction angle is negative
8.30.2The CappedWeakPlaneStressUpdate model shall generate an error if the dilation angle is negative
8.30.3The CappedWeakPlaneStressUpdate model shall generate an error if the friction angle is less than the dilation angle
8.30.4The CappedWeakPlaneStressUpdate model shall generate an error if the cohesion is negative
8.30.5The CappedWeakPlaneStressUpdate model shall generate an error if the sum of the tensile and compressive strength is less than smoothing_tol
8.30.6The CappedWeakPlaneStressUpdate model shall generate an error if the normal vector has zero length
8.30.7The CappedWeakPlaneStressUpdate model shall correctly compute stresses in the elastic regime
8.30.8The CappedWeakPlaneStressUpdate model shall correctly represent tensile failure with the Lame coefficient lambda=0
8.30.9The CappedWeakPlaneStressUpdate model shall correctly represent tensile failure with the Lame coefficient lambda=4
8.30.10The CappedWeakPlaneStressUpdate model shall correctly represent compression failure
8.30.11The CappedWeakPlaneStressUpdate model shall correctly represent shear failure
8.30.12The CappedWeakPlaneStressUpdate model shall correctly represent both tensile and shear failure
8.30.13The CappedWeakPlaneStressUpdate model shall correctly represent tensile behavior with hardening
8.30.14The CappedWeakPlaneStressUpdate model shall correctly represent compression behavior with hardening
8.30.15The CappedWeakPlaneStressUpdate model shall correctly represent shear behavior with hardening
8.30.16The CappedWeakPlaneStressUpdate model shall correctly represent hardening under combined tension and shear
8.30.17The CappedWeakPlaneStressUpdate model shall correctly represent hardening under combined tension and shear with an initial stress
8.30.18The CappedWeakPlaneStressUpdate model shall correctly represent the behavior of a column of elements that is pulled, then pushed
8.30.19The CappedWeakPlaneStressUpdate model shall correctly represent the behavior of a column of elements that is pulled, then pushed, with tensile hardening
8.30.20The CappedWeakPlaneStressUpdate model shall correctly represent the behavior of a beam with its ends fully clamped
8.30.21The CappedWeakPlaneStressUpdate model shall correctly represent the tensile failure of a single layer of elements in 1 nonlinear step
8.30.22The CappedWeakPlaneStressUpdate model shall correctly represent a dynamic problem with plasticity in which a column of material is pulled in tension
8.30.23The CappedWeakPlaneStressUpdate model shall correctly represent a dynamic problem with plasticity in which a column of material is pushed in compression
8.30.24The system shall permit exceptions to be thrown from material models with stateful properties without reading/writing to/from uninitialized memory
8.30.25The CappedWeakInclinedPlaneStressUpdate model shall correctly represent tensile failure with a specified normal=(1,0,0)
8.30.26The CappedWeakPlaneStressUpdate model shall correctly represent tensile failure with a specified normal=(0,1,0)
8.30.27The CappedWeakPlaneStressUpdate model shall correctly represent shear failure with a specified normal=(1,0,0)
8.30.28The CappedWeakPlaneCosseratStressUpdate model shall correctly represent plastic behavior under a first set of loading conditions
8.30.29The CappedWeakPlaneCosseratStressUpdate model shall correctly represent plastic behavior under a second set of loading conditions
8.30.30The CappedWeakPlaneCosseratStressUpdate model shall correctly represent plastic behavior under a third set of loading conditions
8.30.31The CappedWeakPlaneCosseratStressUpdate model shall correctly represent plastic behavior under a fourth set of loading conditions

solid_mechanics: Central Difference
8.31.1The NewmarkBeta timeintegrator shall correctly calculate the response of a 1D mesh.
8.31.2The CentralDifference timeintegrator shall correctly calculate the response of a 1D mesh when the the consistent mass matrix option is used.
8.31.3The NewmarkBeta timeintegrator shall correctly calculate the response of a 2D mesh.
8.31.4The CentralDifference timeintegrator shall correctly calculate the response of a 2D mesh when the the consistent mass matrix option is used.
8.31.5The NewmarkBeta timeintegrator shall correctly calculate the response of a 3D mesh.
8.31.6The CentralDifference timeintegrator shall correctly calculate the response of a 3D mesh when the the consistent mass matrix option is used.
8.31.7The central difference time integrator time step can be increased by selecting density scaling.
8.31.8The system shall include central difference time integration that correctly calculates the response of a 1D mesh with nodal masses equal to those of a corresponding lumped mass system.
8.31.9The system shall include Newmark-beta time integration that correctly calculate the response of a 1D mesh with nodal masses equal to those of a corresponding lumped mass system.
8.31.10The system shall include a central difference time integration that when used with the lumped mass shall correctly calculate the response of a 1D mesh and produce results that are identical to those calculated using equivalent nodal masses.
8.31.11The system shall include central difference time integration that when used with the constant mass option shall correctly calculate the response of a 1D mesh and produce results that are identical to those calculated using equivalent nodal masses.
8.31.12The system shall include central difference time integration that correctly calculates the response of a 2D mesh with nodal masses equal to those of a corresponding lumped mass system.
8.31.13The system shall include Newmark-beta time integration that correctly calculates the response of a 2D mesh with nodal masses equal to those of a corresponding lumped mass system.
8.31.14The system syall include central difference time integration that when used with the lumped mass option shall correctly calculate the response of a 2D mesh and produce results that are identical to those calculated using equivalent nodal masses.
8.31.15The system syall include central difference time integration that when used with the constant mass option shall correctly calculate the response of a 2D mesh and produce results that are identical to those calculated using equivalent nodal masses.
8.31.16The system shall include central difference time integration that correctly calculates the response of a 3D mesh with nodal masses equal to those of a corresponding lumped mass system.
8.31.17The system shall include Newmar-beta time integration that correctly calculates the response of a 3D mesh with nodal masses equal to those of a corresponding lumped mass system.
8.31.18The system shall include central difference time integration that when used with the lumped mass option shall correctly calculate the response of a 3D mesh and produce results that are identical to those calculated using equivalent nodal masses.
8.31.19The system shall include central difference time integration that when used with the constant mass option shall correctly calculate the response of a 3D mesh and produce results that are identical to those calculated using equivalent nodal masses.

solid_mechanics: Check Error
8.32.1The system shall generate an error if a number of elastic constants other than two is supplied for an isotropic elasticity tensor
8.32.2The system shall generate an error if a non-positive Youngs modulus is supplied for an isotropic elasticity tensor
8.32.3The system shall generate an error if a non-positive bulk modulus is supplied for an isotropic elasticity tensor
8.32.4The system shall generate an error if a Poissons ratio outside the range from -1 to 0.5 is supplied for an isotropic elasticity tensor
8.32.5The system shall generate an error if a non-positive shear modulus is supplied for an isotropic elasticity tensor

solid_mechanics: Cohesive Zone Model
8.33.1
The system shall be able to use the small strain and total lagrangian cohesive zone model formulation, demonstrate an exact jacobian, and use the parent action to simplify input an doutput and provide the correct results
1. The small strain cohesive zone model shall converge quadratically when using the hand-coded Jacobian
2. The finite strain cohesive zone model shall converge quadratically when using the hand-coded Jacobian
3. The small strain cohesive zone model shall return the same results of the total formulation when using an incremental material
4. The finite strain cohesive zone model shall return the same results of the total formulation when using an incremental material
5. The system shall allow using multiple CohesiveZoneAction and different cohesive materials for different boundaries
6. The system shall report an error when different strains are provided in different blocks of the CohesiveZoneAction
7. The system shall report an error when overlapping boundaries are provided in different blocks of the CohesiveZoneAction
8. The system shall use the common block generate output and additional outputs for all CZM blocks
9. The system shall add additional outputs requested within each block
10. The system shall be able to use different output orders on different blocks
11. 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.
12. 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.
13. 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.
14. 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.
15. 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.
16. 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.
17. The CZMInterfaceKernelTotalLagrangian should properly account for large area changes and interface rotations when using an incremental material.
18. The CZMInterfaceKernelTotalLagrangian should properly account for large area changes and interface rotations when using a total strain material.
8.33.2
The system shall be able to model decohesion of composites using cohesive zone model
1. under mode I loading.
2. under mode I loading with scaled normal strength.
3. under mode II loading.
4. under mode I loading by lagging the seperation state.
5. under mode I loading with power law propagation criterion.
8.33.3
The system shall be able to use the small strain and total lagrangian cohesive zone model formulation, demonstrate an exact jacobian, and use the parent action to simplify input an doutput and provide the correct results
1. The system shall converge quadratically when using the handcoded Jacobian of the 3DC traction separation model
2. The 3DC traction separation shall produce only a normal traction when a normal displacement jump is applied.
3. The 3DC traction separation shall produce only a shear traction y when a only a displacement jump along y is applied.
4. The 3DC traction separation shall produce only a shear traction z when a only a displacement jump along z is applied.
8.33.4
The system shall be able to use the small strain cohesive zone model formulation using automatic differentiation
1. The solution obtained using AD CZM shall match
2. the non-AD solution
8.33.5
The system shall be able to use the total lagrangian cohesive zone model formulation using automatic differentiation
1. The solution obtained using AD CZM shall match
2. the non-AD solution

solid_mechanics: Combined Creep Plasticity
8.34.1MOOSE tensor mechanics module shall solve a combined creep and plasticity 1-d bar problem.
8.34.2MOOSE tensor mechanics module shall solve a combined creep and plasticity 1-d bar problem with a non-zero start time.
8.34.3MOOSE tensor mechanics module shall solve a combined creep and plasticity 3D cube problem with a time-varying pressure BC.
8.34.4MOOSE tensor mechanics module shall solve a combined creep and plasticity 3D cube problem with a constant displacement BC and stress relaxation.
8.34.5MOOSE tensor mechanics module shall solve a combined creep and plasticity 1-d bar problem with the creep and plasticity models tightly coupled.

solid_mechanics: Coupled Pressure
8.35.1The system shall allow to apply a pressure boundary condition from a variable

solid_mechanics: Crack Front Stress
8.36.1The system shall compute the average stress at each crack point in a 3D domain.

solid_mechanics: Crack Loop
8.37.1The system shall report crack front data for a curved crack front to the console.

solid_mechanics: Creep Tangent Operator
8.38.1The system shall compute the proper stress update using the radial return isotropic power law creep model.
8.38.2The system shall compute the proper stress update using multiple radial return isotropic power law creep models where one of the models returns zero.
8.38.3The system shall compute the sum of multiple stress updates using multiple radial return isotropic power law creep models.
8.38.4The system shall support the cycling of multiple creep models when computing stress updates.
8.38.5The system shall produce the correct Jacobians for radial return isotropic power law creep models.
8.38.6The system shall produce the correct Jacobians for radial return isotropic power law creep models when one of the models returns zero.
8.38.7The system shall produce the correct Jacobians for summed radial return isotropic power law creep models.
8.38.8The systam shall produce the correct Jacobians for cyclic creep model evaluation.

solid_mechanics: Critical Time Step
8.39.1
The system shall correctly compute the critical time step for a solid with:
1. uniform properties and
2. variable properties.
8.39.2The system shall correctly compute the critical time step for a beam.
8.39.3The system shall produce an error if the input elasticity tensor is non-isotropic.

solid_mechanics: Cross Section Deflection
8.40.1The system shall compute the average of the nodal displacements of a cross section defined by the user via a nodal vector.
8.40.2The system shall compute the average of the nodal displacements of a cross section defined by the user via a nodal vector at automatically determined positions.
8.40.3The system shall compute the average of the nodal displacements of a cross section at automatically located points that are updated with mesh adaptivity.
8.40.4The system shall compute the average of the nodal displacements of a cross section at automatically located points that are updated with mesh adaptivity.
8.40.5The system shall generate an error if the user defines an empty set of positions in AverageSectionValueSampler.
8.40.6The system shall generate an error if the user specifies a nonexistent variable in AverageSectionValueSampler.
8.40.7The system shall compute the average of the nodal displacements of a cross section defined by the user via a nodal vector at multiple steps.
8.40.8The system shall compute the average of the nodal displacements of a cross section defined by the user via two nodal vector and two reference points for two respective ducts.
8.40.9The system shall compute the average of the nodal displacements of a cross section with mechanically and thermally induced deformation.
8.40.10The system shall compute the average of the nodal displacements of a cross section with mechanically and thermally induced deformation with a symmetry plane not aligned with a Cartesian coordinate, and match a reference solution without the symmetry plane.

solid_mechanics: Crystal Plasticity
8.41.1The system shall produce correct thermal deformation gradient values with changing temperature
8.41.2The system shall produce correct thermal deformation gradient values when the simulation is restarted
8.41.3The system shall produce correct deformation gradient values for the total and individual eigenstrains
8.41.4The system shall produce correct thermal deformation gradient values with changing temperature under a 011 loading orientation
8.41.5The system shall produce correct thermal deformation gradient values with changing temperature for an anisotropic HCP crystal lattice
8.41.6The system shall produce a positive isotropic eigenstrain in response to an increasing number density and shall demonstrate an applied shear stress on the slip systems with a z-component in response to the BCs on displacement in the z-direction.
8.41.7The system shall continue the isotropic eigenstrain response calculation upon simulation restart.
8.41.8The system shall produce a negative isotropic eigenstrain in response to an decreasing average radius and shall demonstrate an applied stress on the slip systems with a z-component in response to the BC constraints on displacement of the crystal in the z-direction.
8.41.9The system shall produce an isotropic eigenstrain in response to parabolic evolution of both the void number density and average radius, and shall induce plastic deformation in the crystal plasticity reponse through the state dependance of the deformation gradient that final strain state of the FCC crystal is not zero.
8.41.10The system shall not attempt to calculate the volumetric eigenstrain from a negative number density of voids.
8.41.11The system shall not attempt to calculate the volumetric eigenstrain from a negative average radius value.
8.41.12The crystal plasticity models should be able to utilize Euler angle values supplied through AuxVariables.
8.41.13An error will be thrown when trying to set Euler angles information from a text file and AuxVariables at the same time.
8.41.14An error will be thrown when trying to set Euler angles information with less than three components.
8.41.15
The model for HCP crystal plasticity, using the Beyerlein and Tome (2008) constitutive model, shall:
1. 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;
2. evolve the forest dislocation and substructure densities on prismatic<a> slip planes when loaded in an appropriate orientation;
3. evolve the forest dislocation densities on pyramidal <c+a> slip planes and the substructure density when the crystal is loaded in an appropriate orientation.
4. evolve the forest dislocation densities on basal <a> slip planes when the crystal is loaded in an appropriate orientation.
5. prevent the calculation of a negative substructure density increment.
6. prevent the calculation of a negative total forest density.
7. ensure the user specifies a number of the slip systems per mode which is equal to the number of modes;
8. ensure the user specifies an equal total number of slip systems and the sum of the different slip mode systems;
9. ensure the user specifies a number of burgers vectors per mode that is equal to the number of modes;
10. ensure the user specifies a number of slip generation coefficient per mode that is equal to the number of modes;
11. ensure the user specifies a size of slip activation energy per mode that is equal to the specified number of modes;
12. ensure the user specifies a size of slip energy proportionality factor per mode that is equal to the specified number of modes;
13. ensure the user specifies a size of effective shear modulus per mode that is equal to the specified number of modes;
14. ensure the user specifies a size of substructure generation rate coefficient per mode that is equal to the specified number of modes;
15. ensure the user specifies a size of Hall-Petch like coefficient per mode that is equal to the specified number of modes;
16. ensure the user specifies a size of lattice friction values per mode that is equal to the specified number of modes;
8.41.16
Upon simulation restart, the HCP crystal plasticity capability shall:
1. continue to evolve the forest dislocation densities on basal <a> slip planes.
2. continue to compute the slip system strength and resolved applied shear stress for the prismatic<a> and pyramidal<c+a> slip systems.
8.41.17
The HCP crystal plasticity capability shall:
1. compute the Schmid tensors and resulting resolved applied shear stress for a non-physical representative set of basal, prismatic, and pyramidal slip systems;
2. ensure that the provided lattice parameters match the expectation a = b != c for HCP crystals;
3. ensure that eight indices are given for the Miller-Bravais notation for HCP slip systems;
4. ensure that the Miller-Bravais basal plane indices sum to zero.
5. ensure that the Miller-Bravais indices for the slip direction in the basal plane sum to zero.
6. provide warning messages about constitutive model and slip increment convergence issues when the corresponding input file parameter is set to true.
7. provide warning messages about constitutive model stress convergence issues with the appropriate input file parameter setting.
8.41.18
The crystal plasticity HCP system shall:
1. calculate the Schmid tensor for a set of HCP twin systems given in Miller-Bravais indices and shall calculate a twin volume fraction when using the FCC Kalidindi twinning model with a modified input file.
2. compute a plastic velocity gradient that accounts for deformation contributions from both slip and twinning in an HCP crystal.
8.41.19The crystal plasticity system shall compute the plastic velocity gradient from both slip and twinning in anHCP crystal upon simulation restart.
8.41.20The system shall compute the stress and strain response of a single crystal using the Kalidindi constitutive equation for hardening as a function of slip rate on each slip system.
8.41.21The system shall compute the stress and strain response of a single crystal using the Kalidindi constitutive equation for hardening as a function of slip rate on each slip system with the substepping capability that reduces the intermediate time step size to aid with convergence within the crystal plasticity hardening model.
8.41.22The system shall compute the stress and strain response of a single crystal using the Kalidindi constitutive equation for hardening as a function of slip rate on each slip system with the bisection line search method within the crystal plasticity hardening model to aid with convergence.
8.41.23This is a deprecated system that has been replaced by the stress update material-based crystal plasticity and should be removed.
8.41.24This is a deprecated system that has been replaced by the stress update material-based crystal plasticity and should be removed.
8.41.25This is a deprecated system that has been replaced by the stress update material-based crystal plasticity and should be removed.
8.41.26This is a deprecated system that has been replaced by the stress update material-based crystal plasticity and should be removed.
8.41.27This is a deprecated system that has been replaced by the stress update material-based crystal plasticity and should be removed.
8.41.28This is a deprecated system that has been replaced by the stress update material-based crystal plasticity and should be removed.
8.41.29This is a deprecated system that has been replaced by the stress update material-based crystal plasticity and should be removed.
8.41.30This is a deprecated system that has been replaced by the stress update material-based crystal plasticity and should be removed.
8.41.31This is a deprecated system that has been replaced by the stress update material-based crystal plasticity and should be removed.
8.41.32The system shall provide a StressUpdate material-based extensible crystal plasticity system
8.41.33The system shall provide a StressUpdate material-based extensible crystal plasticity framework that shall continue the state-dependent calculation upon restart
8.41.34The material-based crystal plasticity model shall compute the lagrangian stress, PK2 stress, and slip for an orientation not aligned with the Cartesian coordinate system.
8.41.35A constitutive failure shall trigger an exception leading to a cut time step rather than a failed solve using the stress-update based crystal plasticity implementation
8.41.36The crystal plasticity system shall provide a substepping capability for improved convergence using the stress-update based crystal plasticity implementation
8.41.37The crystal plasticity system shall implement a line search capability for accelerated internal Newton solves using the stress-update based crystal plasticity implementation
8.41.38The crystal plasticity rotations shall correctly rotate the elasticity tensors at each material point using the stress-update based crystal plasticity implementation
8.41.39The Zienkiewicz-Zhu patch shall calculate the stress components at the nodes, with equivalent results in both serial and parallel simulations, in the stress-update based crystal plasticity implementation
8.41.40The system shall apply correct boundary condition after timestep being cut during substepping using the stress-update based crystal plasticity implementation
8.41.41The system shall produce correct Euler angles during the crystal plastic simulation for single crystal oriented along (111).
8.41.42The system shall produce correct Euler angles during the crystal plastic simulation for single crystal oriented along (112).
8.41.43The system shall produce correct Euler angles during the crystal plastic simulation for single crystal oriented along (123).
8.41.44The system shall allow multiple crystal plasticity models to be used in a polycrystalline simulation through the use of the base_name input parameter.
8.41.45The system shall allow multiple crystal plasticity models to be used in a polycrystalline simulation with the restart capability.
8.41.46The system shall produce the initial Euler angles phi1=0, psi=54.7356103, and phi2=45 when supplied with the associated rotation matrix and shall calculate a relatively higher 2nd Piola Kirchoff stress in agreement with the low Schmid factor for this loading orientation.
8.41.47The system shall produce the initial Euler angles of phi1=45, psi=90, and phi2=45, which shall quickly update as the crystal rotates, when supplied with the associated rotation matrix. The system shall calculate a relatively low 2nd Piola Kirchoff stress in the z-direction in agreement with the high Schmid factor for this loading orientation
8.41.48The system shall produce the initial Euler angles of phi1=90, psi=30, and phi2=35, when supplied with the associated rotation matrix, and the updated Euler angles values as the crystal under goes rotation during deformation.
8.41.49The system shall not allow specification of both the rotation matrix and the Bunge Euler angles within the same input file.
8.41.50
The crystal plasticity model for twin propagation in FCC materials, using the Kalidindi (2001) constiutive model, shall:
1. 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.
2. compute a non-zero twin volume fraction, and associate plastic shear increment, for only postive applied shear stress values.
3. calculate the coplanar twin propogation resistance contribution as a function of the twin volume fractions on coplanar systems alone.
4. 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.
5. set the plastic shear increments due to twinning to zero once the total twin volume fraction limit is reached.
6. prevent the accumulation of a total twin volume fraction above the user-defined limit.
7. prevent the calculation of a plastic shear increment in excess of the user-defined tolerance value.
8.41.51
When combined with a dislocation slip constitutive model in a simulation of an FCC material, the Kalidindi (2001) twin propagation model shall:
1. compute plastic shear due to both slip and twin propagation when loaded under tension in a 111 orientation.
2. compute a near-zero twin volume fraction when loaded under compression in a 111 orientation.
3. compute a near-zero volume fraction of twins for all twin system when a 100 oriented crystal is loaded under tension.
4. compute both dislocation glide and twin propagation in a 100 oriented crystal loaded under compression.
8.41.52
Upon simulation restart, the crystal plasticity capability shall:
1. continue to compute the plastic shear due to twinning after a simulation restart.
2. continue to calculate plastic shear due to both dislocation slip and twinning when loaded under tension in a 111 orientation in a simulation restart.
8.41.53The system shall provide a plug-in based extensible crystal plasticity system
8.41.54The material-based crystal plasticity model shall compute the lagrangian stress, PK2 stress, and slip for an orientation not aligned with the Cartesian coordinate system.
8.41.55A constitutive failure shall trigger an exception leading to a cut time step rather than a failed solve
8.41.56The crystal plasticity system shall provide a function to read slip system parameters from files
8.41.57The crystal plasticity system shall use pluggable user objects to dtermine the plasticity state variable evolution
8.41.58The crystal plasticity system shall make local Euler angles at material points available for output
8.41.59The crystal plasticity system shall provide a substepping capability for improved convergence
8.41.60The crystal plasticity system shall implement a line search capability for accellerated internal Newton solves
8.41.61The crystal plasticity rotations shall correctly rotate the elasticity tensors at each material point
8.41.62The crystal plasticity system shall provide a plugin user object implementing the Voce hardening law
8.41.63The Zienkiewicz-Zhu patch shall calculate the stress components at the nodes, with equivalent results in both serial and parallel simulations, in a crystal plasticity finite strain application.
8.41.64The system shall provide an object to read values from a file and map them onto a mesh besed on mesh block IDs
8.41.65The system shall apply correct boundary condition after timestep being cut during substepping

solid_mechanics: Domain Integral Thermal
8.42.1The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D.
8.42.2The domain integral action shall compute all of the fracture domain integrals including the C integral for problems in 2D.
8.42.3The domain integral action shall compute all of the fracture domain integrals including the interaction integral for problems in 2D.
8.42.4The domain integral action shall compute the stress intensity factor using the interaction integral approach for a problem dominated by thermal strains using an approach that integrates eigenstrain gradients in a general manner.
8.42.5The interaction integral shall account for the contributions of the gradients of arbitrary eigenstrains in 2D
8.42.6The DomainIntegralAction shall generate an error if a user specifies the eigenstrain gradient and also request the J integral
8.42.7The DomainIntegralAction shall generate an error if a user specifies the temperature and also provides the eigenstrain_gradient
8.42.8The interaction integral shall account for the contributions of body forces in 2D
8.42.9The DomainIntegralAction shall generate an error if a user specifies the body force and also request the J integral
8.42.10The domain integral action shall compute all of the fracture domain integrals including the interaction integral for problems in any plane for 2D.
8.42.11The domain integral action shall compute all of the fracture domain integrals including the C(t) integral for problems in 2D.
8.42.12The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D using the instantaneous thermal expansion function eigenstrain.
8.42.13The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D using the mean thermal expansion function eigenstrain.
8.42.14The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D using the instantaneous thermal expansion function eigenstrain.

solid_mechanics: Drucker Prager
8.43.1
8.43.2
8.43.3
8.43.4
8.43.5
8.43.6
8.43.7
8.43.8
8.43.9
8.43.10
8.43.11
8.43.12

solid_mechanics: Dynamics
8.44.1The PresetAcceleration class shall accurately prescribe the acceleration at the given boundary.
8.44.2The PresetAcceleration class shall accurately prescribe the acceleration at the given boundary when the Newmark-Beta time integrator is used to calculate the velocity and acceleration.
8.44.3The LinearNodalConstraint class shall constrain the secondary nodes to move as a linear combination of the primary nodes.
8.44.4The PresetDisplacement class shall accurately prescribe the displacement at the given boundary.
8.44.5The PresetDisplacement class shall accurately prescribe the displacement at the given boundary using the velocity and and acceleration computed using the Newmark-Beta time integrator.
8.44.6The mechanics system shall accurately conduct a static analysis in a small number of time steps to equilibrate the system under gravity before starting the dynamic analysis.
8.44.7The mechanics system shall accurately predict the dynamic response of a linear elastic system with both Rayleigh damping and numerical damping resulting from Hilber-Hughes-Taylor (HHT) time integration.
8.44.8The mechanics system shall accurately predict the dynamic response of a linear elastic system with both Rayleigh damping and numerical damping resulting from Hilber-Hughes-Taylor (HHT) time integration when using the velocity and acceleration computed using the Newmark-Beta time integrator.
8.44.9The mechanics system shall accurately predict the dynamic response of a linear elastic system with a constant Rayleigh damping.
8.44.10The mechanics system shall accurately predict the dynamic response of a linear elastic system with Rayleigh damping provided as a material property.
8.44.11The mechanics system shall accurately predict the dynamic response of a linear elastic system using Hilber-Hughes-Taylor (HHT) time integration.
8.44.12The mechanics system shall accurately predict the dynamic response of a linear elastic system using Hilber-Hughes-Taylor (HHT) time integration when using the dynamic tensor mechanics action.
8.44.13The mechanics system shall accurately predict the dynamic response of a linear elastic system using Newmark time integration.
8.44.14The mechanics system shall accurately predict the dynamic response of a linear elastic system using Newmark time integration and the dynamic tensor mechanics action.
8.44.15The mechanics system shall accurately predict the dynamic response of a linear elastic system using Hilber-Hughes-Taylor (HHT) time integration when velocity and acceleration of the system are calculated using the Newmark-Beta time integrator.
8.44.16The dynamic tensor mechanics action shall support automatic differentiation using the Newmark-Beta time integrator.
8.44.17The solid mechanics app shall support a form of central difference time integration that uses a direct calculation of the acceleration from the residual forces to perform a solution update.
8.44.18The solid mechanics app shall support a form of central difference time integration that uses a direct calculation of the acceleration from the residual forces to perform a solution update. The direct central difference method should properly use a time step averaging method to handle non-constant time steps.
8.44.19The solid mechanics app shall support dirichlet boundary conditions for direct central difference time integration that correctly enforce constant and function boundary conditions
8.44.20The mechanics system shall correctly predict 1D wave propagation in a linear elastic material with numerical damping resulting from Hilber-Hughes-Taylor (HHT) time integration.
8.44.21The mechanics system shall correctly predict 1D wave propagation in a linear elastic material with no numerical or structural damping.
8.44.22The mechanics system shall correctly predict 1D wave propagation in a linear elastic material with both Rayleigh damping and numerical damping resulting from Hilber-Hughes-Taylor (HHT) time integration.
8.44.23The mechanics system shall correctly predict 1D wave propagation in a linear elastic material with both Rayleigh damping and numerical damping resulting from Hilber-Hughes-Taylor (HHT) time integration when automatic differentiation is used.
8.44.24The mechanics system shall correctly compute the Jacobian for 1D wave propagation in a linear elastic material with both Rayleigh damping and numerical damping resulting from Hilber-Hughes-Taylor (HHT) time integration when automatic differentiation is used.
8.44.25The mechanics system shall correctly predict 1D wave propagation in a linear elastic material with both Rayleigh damping and numerical damping resulting from Hilber-Hughes-Taylor (HHT) time integration when using the velocity and acceleration computed using the Newmark-Beta time integrator.
8.44.26The mechanics system shall correctly predict 1D wave propagation in a linear elastic material with Rayleigh damping.
8.44.27The mechanics system shall correctly predict 1D wave propagation in a linear elastic material with Rayleigh damping when using the dynamic tensor mechanics parent action.

solid_mechanics: Eigenstrain
8.45.1
8.45.2
8.45.3

solid_mechanics: Elastic Patch
8.46.1The tensor mechanics module shall have the ability to compute spatially uniform stresses under prescribed linearly varying displacements on a set of irregular hexes.
8.46.2The tensor mechanics module shall have the ability to compute spatially uniform stresses under prescribed linearly varying displacements on a set of irregular hexes using an incremental small-strain calculation.
8.46.3The tensor mechanics module shall have the ability to compute spatially uniform stresses under prescribed linearly varying displacements on a set of irregular hexes using an total small-strain calculation.
8.46.4The tensor mechanics module shall have the ability to compute spatially uniform stresses under prescribed linearly varying displacements on a set of irregular hexes using an incremental small-strain calculation with no displaced mesh created.
8.46.5The tensor mechanics module shall have the ability to compute spatially uniform stresses under prescribed linearly varying displacements on a set of irregular hexes when using volumetric locking correction.
8.46.6The tensor mechanics module shall have the ability to compute spatially uniform stresses under prescribed linearly varying displacements on a set of irregular hexes when running on 2 processors in parallel.
8.46.7The tensor mechanics module shall have the ability to compute spatially uniform stresses under prescribed linearly varying displacements on a set of irregular hexes when employing volumetric locking correction and running on 2 processors in parallel.
8.46.8The tensor mechanics module shall have the ability to compute spatially uniform stresses under prescribed linearly varying displacements on a set of irregular 20-noded quadratic hexes.

solid_mechanics: Elasticitytensor
8.47.1The system shall provide a method for assembling an elasticity tensor from multiple tensor contributions weighted by material properties.
8.47.2The system shall provide a method for rotating a tensor using a rotation matrix.
8.47.3The system shall provide a method for rotating a tensor using a rotation matrix requiring more than one rotation.

solid_mechanics: Elem Prop Read User Object
8.48.1The system shall provide an object to read values from a file and map them onto a mesh based on mesh element IDs
8.48.2The system shall provide an object to read values from a file and map them onto a mesh based on grain IDs determined by a random Voronoi tessellation

solid_mechanics: Finite Strain Elastic
8.49.1The ComputeFiniteStrainElasticStress class shall compute the elastic stress for a finite strain formulation found with the Taylor expansion from Rashid(1993) on a unit 3D cube in a Cartesian system using the QuasiStatic SolidMechanics Physics.
8.49.2The ComputeFiniteStrainElasticStress class shall compute the elastic stress for a finite strain formulation found with the Taylor expansion from Rashid(1993) on a unit 3D cube in a Cartesian system using the volumetric locking correction b-bar formulation.
8.49.3The ComputeMultiPlasticityStress class shall, when supplied with no plastic models, reduce to and produce the solely elastic stress solution for a finite strain fomulation, using the QuasiStatic SolidMechanics Physics.
8.49.4The ComputeMultiPlasticityStress class shall, when supplied with no plastic models, reduce to and produce the solely elastic stress solution for a finite strain fomulation, using the volumetric locking correction b-bar formulation.
8.49.5The ComputeFiniteStrainElasticStress class shall compute the elastic stress based on a finite strain fomulation and then follow the stress as the mesh is rotated 90 degrees in accordance with Kamojjala et al.(2015) using the QuasiStatic SolidMechanics Physics.
8.49.6The ComputeFiniteStrainElasticStress class shall compute the elastic stress based on a finite strain fomulation and then follow the stress as the mesh is rotated 90 degrees in accordance with Kamojjala et al.(2015) using the volumetric locking correction b-bar formulation.
8.49.7The ComputeFiniteStrainElasticStress class shall compute the elastic stress for a finite strain formulation using a direct eigensolution to perform the polar decomposition of stretch and rotation on a unit 3D cube in a Cartesian system using the QuasiStatic SolidMechanics Physics.
8.49.8The ComputeFiniteStrainElasticStress class shall compute the elastic stress for a finite strain formulation using a direct eigensolution to perform the polar decomposition of stretch and rotation on a unit 3D cube in a Cartesian system using the volumetric locking correction b-bar formulation.
8.49.9The ComputeLinearElasticStress class shall generate an error if a user attempts to run a problem using ComputeLinearElasticStress with a finite strain formulation.

solid_mechanics: Finite Strain Elastic Anisotropy
8.50.1Finite strain methods in Tensor Mechanics should be able to adequately simulate a complex strain state simulation in 3D using an orhotropic filling with isotropic properties.
8.50.2Finite strain methods in Tensor Mechanics should be able to adequately simulate a complex strain state simulation in 3D using an orhotropic filling with isotropic properties.
8.50.3Finite strain methods in Tensor Mechanics should be able to rotate an orthotropic beam-like element 90 degrees and retrieve the proper displacement after being solicited by a pressure boundary condition.
8.50.4Finite strain methods in Tensor Mechanics should be able to rotate an orthotropic beam-like element 90 degrees and retrieve the proper displacement after being solicited by a pressure boundary condition when automatic differentiation is used.
8.50.5Finite strain methods in Tensor Mechanics should be able to rotate an orthotropic beam-like element 360 degrees and retrieve the proper displacement after being solicited by a pressure boundary condition.
8.50.6Finite strain methods in Tensor Mechanics should be able to rotate an orthotropic beam-like element 360 degrees and retrieve the proper displacement after being solicited by a pressure boundary condition when automatic differentiation is used.

solid_mechanics: Finite Strain Jacobian
8.51.1Finite strain methods in Tensor Mechanics should be able to adequately simulate a bar bending simulation in 2D
8.51.2Finite strain methods in Tensor Mechanics should be able to adequately simulate a bar bending simulation in 2D using a volumetric locking correction
8.51.3Finite strain methods in Tensor Mechanics should be able to adequately simulate a tensile test simulation in 3D
8.51.4Finite strain methods in Tensor Mechanics should be able to adequately simulate a tensile test simulation in 3D using a volumetric locking correction

solid_mechanics: Finite Strain Tensor Mechanics Tests
8.52.1The system shall track a changing global stress state when a model undergoes rigid body rotation
8.52.2The system shall compute a uniform stress state given a uniform strain state with finite strains

solid_mechanics: Free Expansion Abs Ref
8.53.1The tensor mechanics action shall set up the mechanics kernels with the appropriate options to accumulate absolute values of the residual to a tag vector usable for relative convergence checks.

solid_mechanics: Generalized Plane Strain
8.54.1The system shall support a traditional plane strain mechanics solution
8.54.2The system shall support a traditional plane strain mechanics solution where the out-of-plane strain is prescribed
8.54.3The system shall support a generalized plane strain mechanics solution
8.54.4The system shall support a generalized plane strain mechanics solution using the reference residual approach to check solution convergence of the field and scalar variables when the reference vector is constructed using absolute value vector tags.
8.54.5The system shall support a generalized plane strain mechanics solution using the reference residual approach to check solution convergence of the field and scalar variables when the reference vector is constructed using unbalanced kernel contributions.
8.54.6The system shall support a generalized plane strain mechanics solution with incremental strain
8.54.7The system shall support a generalized plane strain mechanics solution with finite strain
8.54.8The system shall support setting up a generalized plane strain problem through an action
8.54.9The system shall support setting the out-of-plane pressure using a function for generalized plane strain problems
8.54.10The system shall support setting the out-of-plane pressure using a material property for generalized plane strain problems
8.54.11The system shall support listing all of the out-of-plane strain variables in the strain calculator

solid_mechanics: Global Strain
8.55.1The globalstrain system shall correctly compute the volume change due to applied stress while still maintaining periodicity in 2D.
8.55.2The globalstrain system shall correctly compute the volume change under uniaxial stress while still maintaining periodicity in all the directions in 3D.
8.55.3The globalstrain system shall correctly compute the volume change under hydrostratic stress while still maintaining periodicity in all the directions in 3D.
8.55.4The globalstrain system shall correctly compute the shear deformation due to applied stress while still maintaining periodicity in all the directions in 3D.
8.55.5The globalstrain system shall correctly compute the deformation behavior in 2D with applied displacement boundary condition in one direction while still maintaining periodicity in the other.
8.55.6The globalstrain system shall correctly compute the deformation behavior in 3D with applied displacement boundary condition in one direction while still maintaining periodicity in the others.
8.55.7The globalstrain system shall correctly compute the deformation behavior in 3D with pressure boundary condition in one direction while still maintaining periodicity in the others.
8.55.8The 'GlobalStrainAction' should set all the objects reqiured for the globalstrain system to correctly compute the deformation behavior maintaining strain periodicity.

solid_mechanics: Gravity
8.56.1The tensor mechanics module shall have the capability of applying a body force term in the stress divergence equilibrium equation that accounts for the force of gravity on a solid object due to its own weight.
8.56.2The tensor mechanics module shall be able to reproduce gravity test results of the hand-coded jacobian using automatic differentiation.
8.56.3The Jacobian for the AD gravity problem shall be perfect
8.56.4The system shall provide a way to apply a vector valued body force density to a volume
8.56.5The system shall provide a way to compute the system total kinetic energy through the use of auxiliary kernels and postprocessors. This test verifies that a body falls at the right acceleration under the action of gravity and that the computed kinetic energy matches the analytical expression.

solid_mechanics: Homogenization
8.57.1The system shall compute homogenized elastic constants using the asymptotic expansion homogenization approach and match values for the so-called long fiber problem
8.57.2The system shall compute homogenized elastic constants using the asymptotic expansion homogenization approach and match values for the so-called short fiber problem

solid_mechanics: Hyperelastic Viscoplastic
8.58.1
8.58.2
8.58.3
8.58.4

solid_mechanics: Ics
8.59.1VolumeWeightedWeibull shall generate a randomly distributed field that approximates the analytic expression for the Weibull distribution when the mesh is uniform and the reference volume is set equal to the element size
8.59.2VolumeWeightedWeibull shall generate a randomly distributed field that approaches the analytic expression for the Weibull distribution when the mesh is uniform and the reference volume is set equal to the element size as the mesh density is increased
8.59.3VolumeWeightedWeibull shall generate a randomly distributed field that approximates the analytic expression for the Weibull distribution when the mesh is uniform, the reference volume is set to a value different from the element size, and the median is adjusted to account for the different reference volume

solid_mechanics: Inclined Bc
8.60.1The SolidMechanics module shall have the capabilty to enforce inclined boundary conditions on a 2D model using a penalty method.
8.60.2The SolidMechanics module shall have the capabilty to enforce inclined boundary conditions on a 3D model using a penalty method.
8.60.3The SolidMechanics module shall have the capabilty to enforce inclined boundary conditions while using automatic differentiation.

solid_mechanics: Inertial Torque
8.61.1The tensor mechanics module computes residual for a simplesituation correctly
8.61.2The tensor mechanics module computes the ith component ofinertial torque where the only degree of freedom in y

solid_mechanics: Initial Stress
8.62.1SolidMechanics 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
8.62.2SolidMechanics 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
8.62.3SolidMechanics shall allow users to specify initial stresses using Functions
8.62.4SolidMechanics shall allow users to specify initial stresses using AuxVariables
8.62.5SolidMechanics shall allow users to specify initial stresses for problems with Cosserat mechanics
8.62.6SolidMechanics shall allow users to specify initial stresses for problems with plasticity, and if the initial stresses are inadmissible, the return-map algorithm will be applied, perhaps incrementally, to bring the initial stresses back to the admissible region

solid_mechanics: Interaction Integral
8.63.1The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in 2D.
8.63.2The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems for all planes in 2D.
8.63.3The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in 3d evaluated as 2D.
8.63.4The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in 3D.
8.63.5The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in 3D while supressing the output of q function values.
8.63.6The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in 3D at specified points.
8.63.7The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in any plane in 3D.
8.63.8The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in 2D while outputting q vlaues.
8.63.9The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in any plane 2D while outputting q values.
8.63.10The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in 3D evaluated as 2D.
8.63.11The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in 3D while outputting q values.
8.63.12The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in 3D for specified points, while outputting q values.
8.63.13The Domain Integral Action shall compute all of the fracture domain integrals including the interaction integral for problems in any plane in 3D while outputting q values.

solid_mechanics: Interaction Integral Benchmark
8.64.1The system shall compute mixed-mode fracture integrals that match values of K fields applied at the boundary of a 3D disk cut on one side by a slit

solid_mechanics: Interface Stress
8.65.1
8.65.2

solid_mechanics: Isotropicsd Plasticity
8.66.1
8.66.2

solid_mechanics: Isotropic Elasticity Tensor
8.67.1The ComputeIsotropicElasticityTensor class shall correctly compute the elasticity tensor from the lambda and shear modulus for an isotropic material.
8.67.2The ComputeIsotropicElasticityTensor class shall correctly compute the elasticity tensor from the Young's modulus and Poisson's ratio for an isotropic material.
8.67.3The ComputeIsotropicElasticityTensor class shall correctly compute the elasticity tensor from their bulk modulus and shear modulus for an isotropic material.
8.67.4The ComputeElasticityTensor class shall correctly compute the elasticity tensor for an isotropic axisymmetric problem.

solid_mechanics: J2 Plasticity
8.68.1
8.68.2
8.68.3
8.68.4
8.68.5
8.68.6

solid_mechanics: J Integral
8.69.1The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D.
8.69.2The domain integral action shall compute the fracture domain integrals including the J integral for problems in which the fracture domains of interest represent a subset of all the domains in the system.
8.69.3The domain integral action shall error out if the blocks selected do not have the adequate material objects.
8.69.4The domain integral action shall suppress the vector postprocessor output and stop generating a csv file for each integration ring at each time step.
8.69.5The domain integral action shall suppress the vector postprocessor output while the postprocessor outputs are not affected.
8.69.6The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D using small strain.
8.69.7The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D at specified points.
8.69.8The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D given a mouth direction.
8.69.9The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D using the topology type q function.
8.69.10The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D evaluated as a 2D problem.
8.69.11The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D evaluated as a 2D problem using the topology type q function.
8.69.12The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D.
8.69.13The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D with the q function turned off.
8.69.14The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D with specified points.
8.69.15The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D given a crack mouth direction.
8.69.16The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D given a crack mouth direction and end direction vector.
8.69.17The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D with a topology type q function.
8.69.18The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D evaluated as a 2D problem using the topology type q function.
8.69.19The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D given a crack mouth direction.
8.69.20The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D while supressing the output of the q function values.
8.69.21The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D while outputting the q function values.
8.69.22The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 2D with the topology type q function and outputting the values.
8.69.23The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D while supressing the output of the q values.
8.69.24The domain integral action shall compute all of the fracture domain integrals including the J integral for problems in 3D with the topology type q function and outputting the values.

solid_mechanics: J Integral Vtest
8.70.1The Domain Integral Action shall compute all of the fracture domain integrals including the J integral for surface breaking elliptical cracks.
8.70.2The Domain Integral Action shall compute all of the fracture domain integrals including the J integral for surface breaking elliptical cracks using the crack mouth specification.
8.70.3The Domain Integral Action shall compute all of the fracture domain integrals including the J integral for surface breaking elliptical cracks using the crack mouth specification computing the system Jacobian via automatic differentiation.
8.70.4The Domain Integral Action shall compute all of the fracture domain integrals including the J integral for surface breaking elliptical cracks with crack face pressure.
8.70.5The Domain Integral Action shall compute all of the fracture domain integrals including the J integral for surface breaking elliptical cracks with crack face pressure and crack mouth boundary specified.
8.70.6The Domain Integral Action shall compute all of the fracture domain integrals including the C integral for surface breaking elliptical cracks.
8.70.7The Domain Integral Action shall compute all of the fracture domain integrals including the C integral for surface breaking elliptical cracks using automatic differentiation.
8.70.8The Domain Integral Action shall compute stress intensity factors from the interaction integral for a crack perpendicular to a bi-material interface, treating the interface as a functionally-graded material.
8.70.9The 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.
8.70.10The Domain Integral Action shall compute stress intensity factors from the interaction integral, with an added term stemming for spatially varying properties, that accounts for the crack perpendicularity to the property grading and yield verified results.
8.70.11The Domain Integral Action shall compute verified stress intensity factors from the interaction integral for axisymmetric geometries and circumferential cracks.

solid_mechanics: Jacobian
8.71.1
8.71.2
8.71.3
8.71.4
8.71.5
8.71.6
8.71.7
8.71.8
8.71.9
8.71.10
8.71.11
8.71.12
8.71.13
8.71.14
8.71.15
8.71.16
8.71.17
8.71.18
8.71.19
8.71.20
8.71.21
8.71.22
8.71.23
8.71.24
8.71.25
8.71.26
8.71.27
8.71.28
8.71.29
8.71.30
8.71.31
8.71.32
8.71.33
8.71.34
8.71.35
8.71.36
8.71.37
8.71.38
8.71.39
8.71.40
8.71.41
8.71.42
8.71.43
8.71.44
8.71.45
8.71.46
8.71.47
8.71.48
8.71.49
8.71.50
8.71.51
8.71.52
8.71.53
8.71.54
8.71.55
8.71.56
8.71.57
8.71.58
8.71.59
8.71.60
8.71.61
8.71.62
8.71.63
8.71.64
8.71.65
8.71.66
8.71.67
8.71.68
8.71.69
8.71.70
8.71.71
8.71.72
8.71.73
8.71.74
8.71.75
8.71.76
8.71.77
8.71.78
8.71.79
8.71.80
8.71.81
8.71.82
8.71.83
8.71.84
8.71.85
8.71.86
8.71.87
8.71.88
8.71.89
8.71.90
8.71.91
8.71.92
8.71.93
8.71.94
8.71.95
8.71.96

solid_mechanics: Jacobian Damper
8.72.1The module shall be able to apply a damper on the displaced mesh to limit the element Jacobian change over nonlinear iterations.
8.72.2The element Jacobian damper shall be able to apply a partial nonlinear update.
8.72.3The module shall be able to apply a damper on the undisplaced mesh to limit the element Jacobian change over nonlinear iterations.
8.72.4The module shall be able to apply a damper on certain blocks of the undisplaced mesh to limit the element Jacobian change over nonlinear iterations.

solid_mechanics: Lagrangian
8.73.1Can run second order problem without stabilization
8.73.2Cannot run second order problem with stabilization
8.73.3Homogenization system can be block restricted
8.73.4Action produces identical results for small deformations/without stabilization in axisymmetric RZ coordinates
8.73.5Action produces identical results for small deformations/with stabilization in axisymmetric RZ coordinates
8.73.6Action produces identical results for large deformations/without stabilization in axisymmetric RZ coordinates
8.73.7Action produces identical results for large deformations/with stabilization in axisymmetric RZ coordinates
8.73.8Small strain solution in 2D axisymmetric RZ coordinates matches analytical solution
8.73.9Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Dirichlet BCs in 2D axisymmetric RZ coordinates
8.73.10Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Dirichlet BCs in 2D axisymmetric RZ coordinates
8.73.11Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Dirichlet BCs in 2D axisymmetric RZ coordinates
8.73.12Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Dirichlet BCs in 2D axisymmetric RZ coordinates
8.73.13Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Neumann BCs in 2D axisymmetric RZ coordinates
8.73.14Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Neumann BCs in 2D axisymmetric RZ coordinates
8.73.15Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Neumann BCs in 2D axisymmetric RZ coordinates
8.73.16Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Neumann BCs in 2D axisymmetric RZ coordinates
8.73.17Small strain uniaxial patch test in 2D axisymmetric RZ coordinates passes
8.73.18Large strain uniaxial patch test in 2D axisymmetric RZ coordinates passes
8.73.19Jacobian is exact with relative tolerance of 5e-6 for small displacements in 2D axisymmetric cylindrical coordinates, no stabilization, including cross-terms
8.73.20Jacobian is exact with relative tolerance of 5e-6 for large displacements in 2D axisymmetric cylindrical coordinates, no stabilization, including cross-terms
8.73.21Jacobian is exact with relative tolerance of 5e-6 for small displacements in 2D axisymmetric cylindrical coordinates, with stabilization, including cross-terms
8.73.22Jacobian is exact with relative tolerance of 5e-6 for large displacements in 2D axisymmetric cylindrical coordinates, with stabilization, including cross-terms
8.73.23Action produces identical results for large deformations/with stabilization
8.73.24Action produces identical results for large deformations/no stabilization
8.73.25Action produces identical results for small deformations/with stabilization
8.73.26Action produces identical results for small deformations/no stabilization
8.73.27Action produces identical results for 1D Cartesian problem
8.73.28Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Dirichlet BCs
8.73.29Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Neumann BCs
8.73.30Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Dirichlet BCs
8.73.31Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Neumann BCs
8.73.32Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Dirichlet BCs
8.73.33Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Neumann BCs
8.73.34Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Dirichlet BCs
8.73.35Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Neurmann BCs
8.73.36Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Dirichlet BCs
8.73.37Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Neumann BCs
8.73.38Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Dirichlet BCs
8.73.39Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Neumann BCs
8.73.40Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Dirichlet BCs
8.73.41Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Neumann BCs
8.73.42Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Dirichlet BCs
8.73.43Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Neurmann BCs
8.73.44Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Dirichlet BCs
8.73.45Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Neumann BCs
8.73.46Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Dirichlet BCs
8.73.47Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Neumann BCs
8.73.48Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Dirichlet BCs
8.73.49Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Neumann BCs
8.73.50Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Dirichlet BCs
8.73.51Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Neurmann BCs
8.73.52Quadratic convergence in 1D for DirichletBCs, SMP does better than FDP for this problem.
8.73.53Quadratic convergence in 1D for NeumannBCs, SMP does better than FDP for this problem.
8.73.54Quadratic convergence in 2D for DirichletBCs, SMP matches FDP.
8.73.55Quadratic convergence in 2D for NeumannBCs, SMP matches FDP.
8.73.56Quadratic convergence in 3D for DirichletBCs, SMP matches FDP.
8.73.57Quadratic convergence in 3D for NeumannBCs, SMP matches FDP.
8.73.58Quadratic convergence in 3D for PresetBCs on a complicated geometry, SMP matches FDP.
8.73.59Quadratic convergence in 3D for PresetBCs on a complicated geometry, SMP matches FDP.
8.73.60Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and an elastic material
8.73.61Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and an elastic material
8.73.62Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and an elastic material
8.73.63Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and an elastic material
8.73.64Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and an plastic_j2 material
8.73.65Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and an plastic_j2 material
8.73.66Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and an plastic_j2 material
8.73.67Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and an plastic_j2 material
8.73.68Resulting flow curve matches analytic calculation
8.73.69Action produces identical results for a 3D, large deformation problem with mixed constraints compared to a case where the components are specified manually.
8.73.70Action produces identical results for a 2D, small deformation problem with mixed constraints compared to a case where the components are specific manually.
8.73.71Jacobian is exact with rtol 1e-7 for small deformations, imposed strains, and no stabilization
8.73.72Jacobian is exact with rtol 1e-7 for small deformations, imposed stresses, and no stabilization
8.73.73Jacobian is exact with rtol 1e-7 for large deformations, imposed strains, and no stabilization
8.73.74Jacobian is exact with rtol 1e-7 for large deformations, imposed stresses, and no stabilization
8.73.75Jacobian is exact with rtol 1e-7 for small deformations, imposed strains, and stabilization on
8.73.76Jacobian is exact with rtol 1e-7 for small deformations, imposed stresses, and stabilization on
8.73.77Jacobian is exact with rtol 1e-6 for large deformations, imposed strains, and stabilization on
8.73.78Jacobian is exact with rtol 1e-7 for large deformations, imposed stresses, and stabilization on
8.73.79Homogenized strain field hits the targets in 1D using large kinematics
8.73.80Homogenized stress field hits the targets in 1D using large kinematics
8.73.81Homogenization with strain constraints hits the targets in 2D using large kinematics
8.73.82Homogenization with stress constraints hits the targets in 2D using large kinematics
8.73.83Homogenization with strain constraints hits the targets in 3D using large kinematics
8.73.84Homogenization with stress constraints hits the targets in 3D using large kinematics
8.73.85Homogenization from scalar wrapper class with strain constraints hits the targets in 2D
8.73.86Homogenization from scalar wrapper class with stress constraints hits the targets in 2D
8.73.87Homogenization from scalar wrapper class with mixed stress and strain constraints hits the targets in 2D
8.73.88Framework scalar kernel wrapper correctly assembles scalar-to-scalar coupling Jacobian
8.73.89Framework scalar kernel wrapper correctly assembles Jacobian by rows
8.73.90Homogenized strain field hits the targets in 1D
8.73.91Homogenized stress field hits the targets in 1D
8.73.92Homogenization with strain constraints hits the targets in 2D
8.73.93Homogenization with stress constraints hits the targets in 2D
8.73.94Homogenization with mixed stress and strain constraints hits the targets in 2D
8.73.95Homogenization with strain constraints hits the targets in 3D
8.73.96Homogenization with stress constraints hits the targets in 3D
8.73.97Homogenization with mixed stress and strain constraints hits the targets in 3D
8.73.98Small strain uniaxial patch test passes
8.73.99Large strain uniaxial patch test passes
8.73.100The new system shall produce the correct results if residual and jacobian are computed together.
8.73.1012D generalized plane strain solution shall match the 3D results in the limit of infinite thickness.
8.73.102Weak plane stress Jacobian is exact with relative tolerance of 1e-6 for small displacements without stabilization.
8.73.103Weak plane stress Jacobian is exact with relative tolerance of 1e-6 for small displacements with stabilization.
8.73.104Weak plane stress Jacobian is exact with relative tolerance of 1e-6 for large displacements without stabilization.
8.73.105Weak plane stress Jacobian is exact with relative tolerance of 1e-6 for large displacements with stabilization.
8.73.106
Using weak plane stress,
1. the 2D weak plane stress result should match
2. the 3D result where the out-of-plane stress is zero.
8.73.107Classical Truesdell response to pure shear matches known solution
8.73.108Classical Jaumann response to pure shear matches known solution
8.73.109Classical Green-Naghdi response to pure shear matches known solution
8.73.110Classical Truesdell response to pure rotation matches known solution
8.73.111Classical Jaumann response to pure rotation matches known solution
8.73.112Classical Green-Naghdi response to pure rotation matches known solution
8.73.113Treusdell tangent is exact within rtol of 1e-7
8.73.114Jaumann tangent is exact within rtol of 1e-7
8.73.115Green-Naghdi tangent is exact within rtol of 1e-7
8.73.116Stresses are being computed on the deformed volume
8.73.117Stress state correctly rotates from step 10 to the final step
8.73.118Classical demonstration that the shear stress/strain response remains linear for a Truesdell rate.
8.73.119System still functions if you supply a base_name
8.73.120Stabilized solution matches Cook's solution for stabilized Q4 elements.
8.73.121Stabilized solution matches Cook's solution for stabilized Q4 elements.
8.73.122Strains are as-expected (0.1 total, 0 mechanical) for a free expansion
8.73.123Stresses and strain as expected for thermal expansion constrained in one direction
8.73.124Jacobian is exact with relative tolerance of 1e-7 for small displacements, no stabilization, including cross-terms
8.73.125Jacobian is exact with relative tolerance of 1e-7 for large displacements, no stabilization, including cross-terms
8.73.126Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, including cross-terms
8.73.127Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, including cross-terms
8.73.128Action produces identical results for large deformations/with stabilization
8.73.129Action produces identical results for large deformations/no stabilization
8.73.130Action produces identical results for small deformations/with stabilization
8.73.131Action produces identical results for small deformations/no stabilization
8.73.132Action produces identical results for 1D Cartesian problem
8.73.133Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Dirichlet BCs
8.73.134Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Neumann BCs
8.73.135Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Dirichlet BCs
8.73.136Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Neumann BCs
8.73.137Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Dirichlet BCs
8.73.138Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Neumann BCs
8.73.139Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Dirichlet BCs
8.73.140Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Neurmann BCs
8.73.141Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Dirichlet BCs
8.73.142Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Neumann BCs
8.73.143Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Dirichlet BCs
8.73.144Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Neumann BCs
8.73.145Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Dirichlet BCs
8.73.146Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Neumann BCs
8.73.147Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Dirichlet BCs
8.73.148Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Neurmann BCs
8.73.149Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Dirichlet BCs
8.73.150Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Neumann BCs
8.73.151Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Dirichlet BCs
8.73.152Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Neumann BCs
8.73.153Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Dirichlet BCs
8.73.154Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Neumann BCs
8.73.155Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Dirichlet BCs
8.73.156Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Neurmann BCs
8.73.157Quadratic convergence in 1D for DirichletBCs, SMP does better than FDP for this problem.
8.73.158Quadratic convergence in 1D for NeumannBCs, SMP does better than FDP for this problem.
8.73.159Quadratic convergence in 2D for DirichletBCs, SMP matches FDP.
8.73.160Quadratic convergence in 2D for NeumannBCs, SMP matches FDP.
8.73.161Quadratic convergence in 3D for DirichletBCs, SMP matches FDP.
8.73.162Quadratic convergence in 3D for NeumannBCs, SMP matches FDP.
8.73.163The system shall achieve quadratic convergence in 3D with preset Dirichlet displacement boundary conditions on a complicated geometry for small deformation, and the preconditioning for a single matrix shall match finite-difference preconditioning results.
8.73.164The system shall achieve quadratic convergence in 3D with preset Dirichlet displacement boundary conditions for large deformation, and the preconditioning for a single matrix shall match finite-difference preconditioning results.
8.73.165Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and an elastic material
8.73.166Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and an elastic material
8.73.167Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and an elastic material
8.73.168Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and an elastic material
8.73.169Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and an plastic_j2 material
8.73.170Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and an plastic_j2 material
8.73.171Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and an plastic_j2 material
8.73.172Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and an plastic_j2 material
8.73.173Resulting flow curve matches analytic calculation
8.73.174Small strain uniaxial patch test passes
8.73.175Large strain uniaxial patch test passes
8.73.176Stresses are being computed on the deformed volume
8.73.177Stress state correctly rotates from step 10 to the final step
8.73.178Classical demonstration that the shear stress/strain response remains linear for a Truesdell rate.
8.73.179System still functions correctly if you supply a base_name
8.73.180Stabilized solution matches Cook's solution for stabilized Q4 elements for small deformations.
8.73.181Stabilized solution matches Cook's solution for stabilized Q4 elements for large deformations.
8.73.182Strains are as-expected (0.1 total, 0 mechanical) for a free expansion
8.73.183Stresses and strain as expected for thermal expansion constrained in one direction
8.73.184Jacobian is exact with relative tolerance of 1e-7 for small displacements, no stabilization, including cross-terms
8.73.185Jacobian is exact with relative tolerance of 1e-7 for large displacements, no stabilization, including cross-terms
8.73.186Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, including cross-terms
8.73.187Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, including cross-terms
8.73.188Action produces identical results for small deformations/without stabilization in axisymmetric RZ coordinates
8.73.189Action produces identical results for large deformations/without stabilization in axisymmetric RZ coordinates
8.73.190Small strain solution in 1D centrosymmetric spherical coordinates matches analytical solution
8.73.191Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Dirichlet BCs in 1D centrosymmetric spherical coordinates
8.73.192Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Dirichlet BCs in 1D centrosymmetric spherical coordinates
8.73.193Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Dirichlet BCs in 1D centrosymmetric spherical coordinates
8.73.194Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Dirichlet BCs in 1D centrosymmetric spherical coordinates
8.73.195Jacobian is exact with relative tolerance of 1e-7 for small displacements, without stabilization, and Neumann BCs in 1D centrosymmetric spherical coordinates
8.73.196Jacobian is exact with relative tolerance of 1e-7 for small displacements, with stabilization, and Neumann BCs in 1D centrosymmetric spherical coordinates
8.73.197Jacobian is exact with relative tolerance of 1e-7 for large displacements, without stabilization, and Neumann BCs in 1D centrosymmetric spherical coordinates
8.73.198Jacobian is exact with relative tolerance of 1e-7 for large displacements, with stabilization, and Neumann BCs in 1D centrosymmetric spherical coordinates
8.73.199Jacobian is exact with relative tolerance of 5e-6 for small displacements in 1D centrosymmetric spherical coordinates, no stabilization, including cross-terms
8.73.200Jacobian is exact with relative tolerance of 5e-6 for large displacements in 1D centrosymmetric spherical coordinates, no stabilization, including cross-terms
8.73.201Jacobian is exact with relative tolerance of 5e-6 for small displacements in 1D centrosymmetric spherical coordinates, with stabilization, including cross-terms
8.73.202Jacobian is exact with relative tolerance of 5e-6 for large displacements in 1D centrosymmetric spherical coordinates, with stabilization, including cross-terms
8.73.203Throw an error if you do not provide an isotropic tensor to ComputeStVenantKirchhoffStress
8.73.204Correct Jacobian coming from the material model for nonlinear kinematics with relative error less than 1e-6
8.73.205Correct Jacobian coming from the material model for linear kinematics with relative error less than 1e-6
8.73.206Correct Jacobian coming from the material model for nonlinear kinematics with relative error less than 1e-6
8.73.207Correct Jacobian coming from the material model for linear kinematics with relative error less than 1e-6
8.73.208Correct Jacobian coming from the material model for nonlinear kinematics with relative error less than 1e-6
8.73.209Correct Jacobian coming from the material model for linear kinematics with relative error less than 1e-6
8.73.210Correct Jacobian coming from the material model for nonlinear kinematics with relative error less than 1e-6
8.73.211Correct Jacobian coming from the material model for linear kinematics with relative error less than 1e-6
8.73.212Correct Jacobian for elastoplasticity with linear hardening with relative error less than 1e-6
8.73.213Correct Jacobian for elastoplasticity with power-law hardening with relative error less than 1e-6
8.73.214Stress matches elastic solution (xx=1e4, rest zero)
8.73.215Stress matches exact solution (xx=1.6144e3, rest zero)
8.73.216Stress matches exact analytic solution, rest zero
8.73.217Stress matches exact solution
8.73.218Stress matches exact analytic solution, rest zero
8.73.219The hypoelastic equivalent of the St. Venant-Kirchhoff model gives the same result as its hyperelastic counterpart.
8.73.220Stress matches exact analytic solution for elastoplasticity with linear hardening
8.73.221Stress matches exact analytic solution for elastoplasticity with power-law hardening
8.73.222The Lagrangian kernels can use a small strain elastic update with a large kinematics formulation

solid_mechanics: Line Material Rank Two Sampler
8.74.1The system shall allow sampling of material properties derived from rank two tensors along a line and output those quantities via a vectorpostprocessor.
8.74.2The system shall allow sampling of scalar material properties along a line and output those quantities via a vectorpostprocessor.

solid_mechanics: Mandel Notation
8.75.1
The system shall solve the continuum mechanics equations in the small strain formulation
1. with full tensors
2. with symmetric Mandel notation tensors
8.75.2
The system shall solve the continuum mechanics equations in the finite strain formulation
1. with full tensors
2. with symmetric Mandel notation tensors

solid_mechanics: Material Limit Time Step
8.76.1The system shall compute a timestep size that is limited by a power law creep model and a specified maximum inelastic strain increment.
8.76.2The system shall not impose a limit on the time step during the initial evaluation when the creep limiting time step option is used.
8.76.3The ScalarMaterialDamage model shall be capable of informing the material-based time step calculation based on the damage evolution
8.76.4The ScalarMaterialDamage model shall be capable of informing the material-based time step calculation based on the number of elements changing damage state in a step
8.76.5The MaterialTimeStepPostprocessor model shall be capable of computing a time step based on both the material time step limited by damage evoluation and the number of elements changing damage state in a step
8.76.6The MaterialTimeStepPostprocessor shall generate an error if an unknown property is requested with the 'time_step_limit' parameter
8.76.7The MaterialTimeStepPostprocessor shall generate an error if neither the material time step limit nor the elements changed limit is specified.
8.76.8The system to limit the analysis time step based on material behavior shall correctly function when used with an isotropic plasticity model.
8.76.9The system to limit the analysis time step based on material behavior shall correctly function when used with an isotropic plasticity model when the MaterialTimeStepPostprocessor is run at initialization.
8.76.10The ComputeMultipleInelasticStress model shall compute a time step equal to the maximum real number if no inelastic model is provided
8.76.11The system shall compute a time step equal to a user input maximum real number if no inelastic model is provided

solid_mechanics: Mean Cap
8.77.1
8.77.2
8.77.3
8.77.4

solid_mechanics: Mean Cap Tc
8.78.1
8.78.2
8.78.3
8.78.4
8.78.5
8.78.6
8.78.7
8.78.8
8.78.9
8.78.10
8.78.11
8.78.12

solid_mechanics: Modal Analysis
8.79.1The module shall be able to perform modal analysis with the Eigenvalue Executioner.

solid_mechanics: Mohr Coulomb
8.80.1
8.80.2
8.80.3
8.80.4
8.80.5
8.80.6
8.80.7
8.80.8
8.80.9
8.80.10
8.80.11
8.80.12
8.80.13
8.80.14
8.80.15
8.80.16
8.80.17
8.80.18
8.80.19
8.80.20
8.80.21
8.80.22
8.80.23
8.80.24
8.80.25
8.80.26
8.80.27
8.80.28
8.80.29
8.80.30
8.80.31
8.80.32
8.80.33

solid_mechanics: Multi
8.81.1
8.81.2
8.81.3
8.81.4
8.81.5
8.81.6
8.81.7
8.81.8
8.81.9
8.81.10
8.81.11
8.81.12
8.81.13
8.81.14
8.81.15
8.81.16
8.81.17
8.81.18
8.81.19
8.81.20
8.81.21
8.81.22
8.81.23
8.81.24
8.81.25
8.81.26
8.81.27
8.81.28
8.81.29
8.81.30
8.81.31
8.81.32
8.81.33
8.81.34
8.81.35

solid_mechanics: Multi Power Law
8.82.1The AD multiple power law creep must reproduce the results of the AD single power law creep: Reference results.
8.82.2The AD multiple power law creep must reproduce the results of the AD single power law creep when a single stress increment is provided by the user.
8.82.3The AD multiple power law creep must reproduce the results of the AD single power law creep with multiple stress intervals with same power law material parameters.
8.82.4The AD multiple power law creep must avoid regression for a simple few-element problem whose solution has stress levels encompassing two different sets of power law parameters
8.82.5Stress thresholds must be provided in increasing order, otherwise an error will inform the user of such a requirement.

solid_mechanics: Multiple Two Parameter Plasticity
8.83.1
8.83.2
8.83.3

solid_mechanics: Neml2
8.84.4The system shall provide a modular way to transfer variables and/or real material property to NEML2 for use as model parameters.
8.84.5The system shall error out while trying to set a parameter that does not exist.
8.84.6The system shall provide for a modular way to compute NEML2 output with respect to specific model parameters and save that it as a MOOSE material property.
8.84.7The system shall provide for a modular way to transfer NEML2 results to material properties..

solid_mechanics: Nodal Patch Recovery
8.85.1The SolidMechanics module shall be able to construct nodal variables from material properties that are defined at quadrature points.

solid_mechanics: Notched Plastic Block
8.86.1
SolidMechanics shall be able to simulate the confined, uniaxial extension of a notched block which has constitutive law described by:
1. unsmoothed capped-Mohr-Coulomb plasticity
2. smoothed capped-Mohr-Coulomb plasticity, with smoothing performed by the novel MOOSE smoothing method described in Wilkins et al
3. unsmoothed Mohr-Coulomb plasticity
4. smoothed Mohr-Coulomb plasticity, with smoothing performed by the Abbo et al method
5. smoothed Mohr-Coulomb plasticity, with smoothing performed by the novel MOOSE smoothing method described in Wilkins et al, using the cosine smoother
6. smoothed Mohr-Coulomb plasticity, with smoothing performed by the novel MOOSE smoothing method described in Wilkins et al, using the poly1 smoother
7. smoothed Mohr-Coulomb plasticity, with smoothing performed by the novel MOOSE smoothing method described in Wilkins et al, using the poly2 smoother
8. smoothed Mohr-Coulomb plasticity, with smoothing performed by the novel MOOSE smoothing method described in Wilkins et al, using the poly3 smoother

solid_mechanics: Orthotropic Plasticity
8.87.1The system shall be capable of simulating materials that exhibit orthotropic plasticity with constant hardening and linear strain applied in the x and y directions.
8.87.2The system shall be capable of simulating materials that exhibit orthotropic plasticity with power rule hardening and linear strain applied in the x direction.

solid_mechanics: Plane Stress
8.88.1The system shall compute the correct Jacobian for plane stress conditions
8.88.2The system shall compute the correct Jacobian for plane stress conditions
8.88.3The system shall compute plane stress conditions with small strains with input provided using the Master action
8.88.4The system shall compute plane stress conditions with small strains
8.88.5The system shall compute plane stress conditions with incremental strains with input provided using the Master action
8.88.6The system shall compute plane stress conditions with incremental strains with input provided using the Master action
8.88.7The system shall compute plane stress conditions with finite strains with input provided using the Master action
8.88.8The system shall compute plane stress conditions with finite strains with input provided using the Master action
8.88.9The system shall compute the response of a 3D cube in uniaxial tension with finite strain to provide a benchmark for a 2D plane stress, finite strain model
8.88.10The system shall compute the response of a cube in uniaxial tension using a 2D plane stress, finite strain model, and produce the same result as a 3D model

solid_mechanics: Poro
8.89.1
8.89.2

solid_mechanics: Porosity
8.90.1
The system shall compute a Porosity material property from a strain tensor
1. using traditional materials.
2. using automatic differentiation materials.
8.90.2
The system shall handle a negative calculated porosity computed from a strain tensor
1. by setting the porosity to zero.
2. by setting the porosity to the initial condition.
3. by throwing an exception.
8.90.3
The system shall handle a nan calcualted porosity computed from a strain tensor
1. by throwing an exception.

solid_mechanics: Postprocessors
8.91.1The system shall be able to compute the force on a sideset.
8.91.2The system shall be able to compute the volume average of rank two tensor component.
8.91.3The system shall be able to compute the volume average of rank two tensor component, which was computed using automatic differentiation.
8.91.4The system shall compute normal displacements on a surface for checking correct application of penalty boundary conditions.

solid_mechanics: Power Law Creep
8.92.1
The system shall gracefully allow exceptions from the radial return mapping algorithm
1. not using automatic differentiation methods.
2. using automatic differentiation methods.
8.92.2
The system shall gracefully correct for out of bounds guess from the radial return mapping algorithm
1. from the high bound.
2. from the low bound.
8.92.3
The system shall provide a power law creep mechanical constitutive model not using AD
1. that works with a finite strain formulation.
2. that works with a coupled plasticty model.
3. that shall be restartable.
4. that shall be restartable with a response that matches the original run.
5. that works with a small strain formulation.
8.92.4
The system shall provide a power law creep mechanical constitutive model using AD
1. that works with a finite strain formulation.
2. that shall be restartable.
3. that shall be restartable with a response that matches the original run.
4. that works with a small strain formulation.
5. calculate a perfect Jacobian
6. calculate a perfect Jacobian with the small strain formulation

solid_mechanics: Preconditioner Reuse
8.93.1Convergence matches previous version of MOOSE without the preconditioner reuse system
8.93.2Preconditioner is reused until the linear iterations exceed the value of reuse_preconditioner_max_its upon which the system recalculates the preconditioner

solid_mechanics: Pressure
8.94.1The Pressure boundary condition action shall create the objects needed to apply pressure boundary conditions on a 3D model as demonstrated by correctly computing the response of an elastic small-strain isotropic unit cube with pressure applied on three faces to create a hydrostatic pressure.
8.94.2The Pressure boundary condition action shall create the objects needed to apply pressure boundary conditions on a 3D model as demonstrated by correctly computing the response of an elastic small-strain isotropic unit cube with pressure applied on three faces to create a hydrostatic pressure using the volumetric locking correction b-bar formulation.
8.94.3The Pressure boundary condition shall compute the correct Jacobian for a problem in 3D space using hex8 elements.
8.94.4The Pressure boundary condition shall compute the correct Jacobian for a problem in 3D space using hex20 elements.
8.94.5The Pressure boundary condition shall compute the correct Jacobian for a problem in RZ coordinates with quad4 elements.
8.94.6The Pressure boundary condition shall compute the correct Jacobian for a problem in RZ coordinates with quad8 elements.

solid_mechanics: Radial Disp Aux
8.95.1The system shall compute the radial component of displacement for axisymmetric cylindrical models
8.95.2The system shall compute the radial component of displacement for 2D Cartesian cylindrical models
8.95.3The system shall compute the radial component of displacement for 3D Cartesian cylindrical models
8.95.4The system shall compute the radial component of displacement for 1D spherical models
8.95.5The system shall compute the radial component of displacement for axisymmetric spherical models
8.95.6The system shall compute the radial component of displacement for 3D Cartesian spherical models
8.95.7The system shall report an error if "origin" is supplied to RadialDisplacementSphereAux when the coordinate system is not Cartesian or axisymmetric
8.95.8The system shall report an error if "axis_vector" is supplied to RadialDisplacementCylinderAux and the model is not 3D Cartesian
8.95.9The system shall report an error if "axis_vector" has length of zero

solid_mechanics: Rate Independent Cyclic Hardening
8.96.1The system shall compute the J2 linear kinematic plasticity stress based on the imposed cyclic nonsymmetric strain
8.96.2The system shall compute the J2 linear kinematic plasticity stress based on the imposedcyclic symmetric strain
8.96.3The system shall compute the J2 nonlinear isotropic plasticity stress based on the imposedcyclic symmetric strain
8.96.4The system shall compute the J2 nonlinear combined plasticity stress based on the imposed cyclic symmetric strain
8.96.5The system shall compute the J2 nonlinear kinematic plasticity stress based on the imposed cyclic nonsymmetric strain
8.96.6The system shall compute the J2 nonlinear kinematic plasticity stress based on the imposed cyclic nonsymmetric strain
8.96.7The system shall compute the J2 nonlinear kinematic plasticity strain based on the imposed cyclic nonsymmetric loading
8.96.8The system shall compute the J2 linear kinematic plasticity strain based on the imposed cyclic nonsymmetric loading

solid_mechanics: Recompute Radial Return
8.97.1The system shall compute the J2 isotropic plasticity stress and plastic strain response under tensile loading within the small incremental strain formulation.
8.97.2The system shall compute the J2 isotropic plasticity stress and plastic strain response under tensile loading within the small incremental strain formulation while prescribing a base name for the isotropic plasticity material properties.
8.97.3The system shall compute the J2 isotropic plasticity stress and plastic strain response under tensile loading within the small incremental strain formulation and using the b-bar element volume correction.
8.97.4The system shall compute the J2 isotropic plasticity stress and plastic strain response under tensile loading within the finite incremental strain formulation.
8.97.5The system shall compute the J2 isotropic plasticity stress and plastic strain response under tensile loading within the finite incremental strain formulation and using the b-bar element volume correction.
8.97.6The system shall compute the hyperbolic visoplastic stress response for a time-dependent linear strain hardening plasticity model in a small incremental strain formulation in a manner equivalent to the ABAQUS result.
8.97.7The 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.
8.97.8The 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.
8.97.9The 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.
8.97.10The system shall return an error if a negative yield stress is supplied when calculating the J2 isotropic plasticity stress response.
8.97.11The system shall return an error if anisotropic elasticity tensor is supplied when the J2 isotropic plasticity stress response calculation is requested.
8.97.12The system shall calculate, with J2 isotropic plasticity, the transient stress eigenvalues with stationary eigenvectors verification test from K. Jamojjala, R. Brannon, A. Sadeghirad, J. Guilkey, Verification tests in solid mechanics, Engineering with Computers, Vol 31., p. 193-213.
8.97.13The system shall calculate, with J2 isotropic plasticity, the transient stress eigenvalues with stationary eigenvectors verification test from K. Jamojjala, R. Brannon, A. Sadeghirad, J. Guilkey, Verification tests in solid mechanics, Engineering with Computers, Vol 31., p. 193-213. This shall be done with a coupled creep and plasticity model where the creep model is inactive.
8.97.14The system shall return an error if the user requests substepping for a radial return inelastic model that does not implement it yet.

solid_mechanics: Rom Stress Update
8.98.1
The system shall compute a creep rate based on a reduced order model
1. and evolve over time with an good finite difference width.
2. and evolve over time with an poor finite difference width with more nonlinear iterations than with a good finite difference width.
3. in 3D.
4. with units of stress other than Pascal.
5. in 2DRz.
6. in 2DRz with a model loaded from a JSON data file.
7. 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.
8. 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.
9. at the lower temperature limit.
10. below the lower temperature limit and extrapolate correctly.
11. but should error out below the lower temperature limit if requested by the user.
12. but should error out above the upper temperature limit if requested by the user.
13. when substepping is off and the time step size does not guarantee the most accurate results.
14. when substepping is off and the time step size guarantees accurate results.
15. when substepping is on and the numerical integration error selected by the user guarantees accurate results.
16. 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.
17. for a variety of overlapping tiles and two overlapping partitions.
18. for a variety of overlapping tiles and two overlapping partitions with a model loaded from a JSON data file.
19. for a variety of overlapping tiles and two overlapping partitions with a model loaded from a JSON data file relative to the input.
20. for a variety of overlapping tiles and two overlapping partitions, and be able to export the model data to a JSON file.
8.98.2
The system shall compute a creep rate based on a reduced order model using AD
1. in 3D.
2. with units of stress other than Pascal.
3. in 3D and compute a perfect Jacobian.
4. in 2DRz.
5. in 2DRz and compute a perfect Jacobian.
6. in isolation (i.e. without a full displacement solve), and match with code-to-code comparison with a small set of input parameters using AD.
7. in isolation (i.e. without a full displacement solve), and match with code-to-code comparison with a large set of input parameters using AD.
8. at the lower temperature limit.
9. below the lower temperature limit and extrapolate correctly.

solid_mechanics: Scalar Material Damage
8.99.1The system shall calculate the effect of damage on the stress of a elastic material.
8.99.2The system shall calculate the effect of damage on the stress of a elastic material when using automatic differentiation.
8.99.3The system shall calculate damaged stress based on old damage index.
8.99.4The system shall calculate damaged stress based on old damage index when using automatic differentiation.
8.99.5The system shall report an error when damage index is greater than 1.
8.99.6The system shall report an error when damage index is greater than 1 also when using automatic differentiation.
8.99.7The 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
8.99.8The 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
8.99.9The system shall calculate the maximum value of the damage index comparing different damage models.
8.99.10The system shall calculate the effective damage index from different damage models.
8.99.11The system shall allow for damage to be limited to a user-defined value for an individual damage model that is combined with another damage model
8.99.12The system shall allow for damage to be limited to a user-defined value for the combination of multiple damage models.
8.99.13The system shall calculate the maximum value of the damage index comparing different damage models when using automatic differentiation.
8.99.14The system shall calculate the effective damage index from different damage models when using automatic differentiation.
8.99.15The system shall allow for damage to be limited to a user-defined value for an individual damage model that is combined with another damage model when using automatic differentiation.
8.99.16The system shall allow for damage to be limited to a user-defined value for the combination of multiple damage models when using automatic differentiation.
8.99.17The system shall calculate the effect of damage on the stress of a inelastic material in conjunction with the creep or plastic deformation.
8.99.18The system shall calculate the effect of damage on the stress of a inelastic material in conjunction with the power law creep deformation.
8.99.19The system shall calculate the effect of damage on the stress of a inelastic material in conjunction with the power law creep deformation with automatic differentiation.
8.99.20The system shall provide a meaningful error if the user selects a material model that is not recognized by the automatic differentiation multiple stress object.
8.99.21The system shall calculate the effect of damage on the stress of a nonlocal damage model.
8.99.22The system shall calculate the effect of damage on the stress of a nonlocal damage model. When using automatic differentiation.

solid_mechanics: Shell
8.100.1
The mechanics system shall accurately predict the dynamic behavior of a shell element with clamped-free-free-free conditions when
1. the surface is on the XY plane
2. the surface is at a 45 deg. angle w.r.t. the Y axis
3. the surface is at a 45 deg. angle w.r.t. the Y axis and HHT input is provided
8.100.2The mechanics system shall accurately predict the dynamic behavior of a shell element.
8.100.3The mechanics system shall accurately compute the deflection of a cantilever beam when it is modeled using shell elements.
8.100.4The mechanics system shall accurately compute the deflection of a rotated cantilever beam when it is modeled using shell elements.
8.100.5The mechanics system shall accurately compute the deflection of a cantilever beam when it is modeled using shell elements under large strain and rotations are included.
8.100.6The mechanics system shall accurately compute the Jacobian for a small strain quasi-static shell element.
8.100.7The mechanics system shall accurately compute the Jacobian for a large strain quasi-static shell element.
8.100.8The mechanics system shall accurately model the deflection of a simply supported under uniform loading.
8.100.9The mechanics system shall accurately model deflection of a plate with multiple force and moment boundary conditions.
8.100.10The mechanics system shall converge to the analytic displacement of a cylinder modeled with shell elements, with point loading applied in the x direction.
8.100.11The mechanics system shall converge to the analytic displacement of a cylinder modeled with shell elements, with point loading applied in the y direction.
8.100.12The mechanics system shall accurately model the deflection of a inclined shell elements.
8.100.13The mechanics system shall accurately compute the axial stress and strain of the shell elements.
8.100.14The mechanics system shall accurately compute the axial stress and strain of the shell elements using the finite strain formulation.
8.100.15The mechanics system shall accurately compute the stress and strain of the shell elements subjected to shear deformation.
8.100.16The mechanics system shall accurately compute the stress and strain of the tapered shell elements.
8.100.17The system shall generate an error if the shell strain material is evaluated when the number of quadrature points does not equal 4.
8.100.18Shell elements should have compatible stresses even for unstructured meshes.

solid_mechanics: Smeared Cracking
8.101.1The MOOSE SolidMechanics module shall simulate cracking on a specimen under tension in cartesian coordinates.
8.101.2The MOOSE SolidMechanics module shall simulate cracking on a specimen under tension in rz coordinates.
8.101.3The MOOSE SolidMechanics module shall simulate cracking while the cracking strength is prescribed by an elemental AuxVariable.
8.101.4The MOOSE SolidMechanics module shall simulate exponential stress release.
8.101.5The MOOSE SolidMechanics module shall simulate exponential stress relase while using the rz coordinate system.
8.101.6The MOOSE SolidMechanics module shall demonstrate softening using the power law for smeared cracking.
8.101.7The MOOSE SolidMechanics module shall demonstrate the prescribed softening laws in three directions, power law (x), exponential (y), and abrupt (z).
8.101.8The MOOSE SolidMechanics module shall simulate smeared cracking in the x y and z directions.
8.101.9The MOOSE SolidMechanics module shall simulate smeared cracking under plane stress conditions.
8.101.10The MOOSE SolidMechanics module shall demonstrate that the smeared cracking model correctly handles finite rotation of cracked elements.
8.101.11The MOOSE SolidMechanics module shall demonstrate the finite rotation of cracked elements where the crack is prescribed in x.
8.101.12The MOOSE SolidMechanics module shall demonstrate the finite rotation of cracked elements where the crack is prescribed in z.
8.101.13The MOOSE SolidMechanics module shall demonstrate the finite rotation of cracked elements where two cracks are prescribed in x and z.
8.101.14The ComputeSmearedCrackingStress model shall generate a warning if 'cracked_elasticity_type' is not specified
8.101.15The ComputeSmearedCrackingStress model shall generate an error if an invalid number of softening models is specified

solid_mechanics: Static Deformations
8.102.1
8.102.2
8.102.3
8.102.4
8.102.5
8.102.6
8.102.7
8.102.8
8.102.9
8.102.10
8.102.11
8.102.12

solid_mechanics: Stickybc
8.103.1
8.103.2
8.103.3

solid_mechanics: Strain Energy Density
8.104.1The system shall be capable of calculating strain energy density incrementally in materials with elastic stress and finite strain.
8.104.2The 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.
8.104.3The system shall be able to create a compliance sensitivity material.
8.104.4The system shall be capable of calculating strain energy density for materials with elastic stress and small strain.
8.104.5The 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.
8.104.6The system shall be capable of calculating strain energy density incrementally in materials with inelastic stress and isotropic plasticity.
8.104.7The system shall be capable of calculating strain energy rate density with elastic stress and finite strain.
8.104.8The system shall be capable of calculating strain energy rate density when using small strain assumptions.
8.104.9The system shall be capable of calculating strain energy rate density when using automatic differentiation and weak plane stress.
8.104.10The system shall be capable of calculating strain energy rate density when using hand-coded Jacobian and weak plane stress.
8.104.11The system shall be capable of calculating strain energy density for materials with a user-specified stress by name and small strain.

solid_mechanics: Stress Recovery
8.105.1The Zienkiewicz-Zhu patch shall calculate the stress components at the nodes, with equivalent results in both serial and parallel simulations, in a small strain application.
8.105.2The Zienkiewicz-Zhu patch shall calculate the stress components at the nodes, with equivalent results in both serial and parallel simulations, in a finite strain application.
8.105.3In areas of high concentration gradients, the Zienkiewicz-Zhu implementation shall recover the specified material property.

solid_mechanics: Substepping
8.106.1The system shall report diverged return mapping solve due to large deformation imposed in one step using automatic differentiation
8.106.2The system shall converge under large deformation owing to the help of substepping using automatic differentiation
8.106.3The system shall converge under large deformation owing to the help of adaptive substepping using automatic differentiation with an increment based substepping criterion
8.106.4The system shall converge under large deformation owing to the help of adaptive substepping using automatic differentiation with an error based substepping criterion
8.106.5The system shall report diverged return mapping solve due to large deformation imposed in one step
8.106.6The system shall converge under large deformation owing to the help of substepping
8.106.7The system shall converge under large deformation owing to the help of adaptive substepping with an increment based substepping criterion
8.106.8The system shall converge under large deformation owing to the help of adaptive substepping with an error based substepping criterion
8.106.9The system shall output verbose information about radial return iterations.

solid_mechanics: T Stress
8.107.1The Domain Integral Action shall compute all of the fracture domain integrals including the T stress for cracks in an infinite plate.
8.107.2The Domain Integral Action shall compute all of the fracture domain integrals including the T stress for an elliptical crack in 3D.
8.107.3The Domain Integral Action shall compute all of the fracture domain integrals including the T stress for an elliptical crack in 3D using automatic differentiation objects.

solid_mechanics: Temperature Dependent Hardening
8.108.1The system shall compute the stress as a function of temperature and plastic strain from user-supplied hardening functions not using automatic differentiation.
8.108.2The system shall compute the stress as a function of temperature and plastic strain from user-supplied hardening functions using automatic differentiation.

solid_mechanics: Tensile
8.109.1
8.109.2
8.109.3
8.109.4
8.109.5
8.109.6
8.109.7
8.109.8
8.109.9
8.109.10
8.109.11
8.109.12
8.109.13
8.109.14
8.109.15
8.109.16
8.109.17
8.109.18
8.109.19
8.109.20
8.109.21
8.109.22
8.109.23
8.109.24
8.109.25
8.109.26
8.109.27
8.109.28
8.109.29
8.109.30
8.109.31
8.109.32
8.109.33
8.109.34
8.109.35
8.109.36
8.109.37
8.109.38
8.109.39

solid_mechanics: Test Jacobian
8.110.1The mechanics system shall correctly compute the jacobian for 3D problems using small strain.
8.110.2The mechanics system shall correctly compute the jacobian for 3D problems using incremental small strain.
8.110.3The mechanics system shall correctly compute the jacobian for 3D problems using finite strain.
8.110.4The mechanics system shall correctly compute the jacobian for 3D problems using small strain and volumetric locking correction.
8.110.5The mechanics system shall correctly compute the jacobian for 3D problems using incremental small strain and volumetric locking correction.
8.110.6The mechanics system shall correctly compute the jacobian for 3D problems using finite strain and volumetric locking correction.
8.110.7The mechanics system shall correctly compute the jacobian for RZ problems using small strain.
8.110.8The mechanics system shall correctly compute the jacobian for RZ problems using incremental small strain.
8.110.9The mechanics system shall correctly compute the jacobian for RZ problems using finite strain.
8.110.10The mechanics system shall correctly compute the jacobian for RZ problems using small strain and volumetric locking correction.
8.110.11The mechanics system shall correctly compute the jacobian for RZ problems using incremental small strain and volumetric locking correction.
8.110.12The mechanics system shall correctly compute the jacobian for RZ problems using finite strain and volumetric locking correction.
8.110.13The mechanics system shall correctly compute the jacobian for planestrain problems using small strain.
8.110.14The mechanics system shall correctly compute the jacobian for planestrain problems using incremental small strain.
8.110.15The mechanics system shall correctly compute the jacobian for planestrain problems using finite strain.
8.110.16The mechanics system shall correctly compute the jacobian for planestrain problems using small strain and volumetric locking correction.
8.110.17The mechanics system shall correctly compute the jacobian for planestrain problems using incremental small strain and volumetric locking correction.
8.110.18The mechanics system shall correctly compute the jacobian for planestrain problems using finite strain and volumetric locking correction.
8.110.19The mechanics system shall correctly compute the jacobian for spherical problems using small strain.
8.110.20The mechanics system shall correctly compute the jacobian for spherical problems using small, incremental strain.
8.110.21The mechanics system shall correctly compute the jacobian for spherical problems using small, incremental strain.

solid_mechanics: Thermal Expansion
8.111.1The tensor mechanics module shall have the capability to calculate the eigenstrain tensor resulting from isotropic thermal expansion.
8.111.2The tensor mechanics module shall have the capability to calculate the eigenstrain tensor resulting from isotropic thermal expansion when restarting the simulation.
8.111.3The tensor mechanics module shall have the capability to calculate the eigenstrain tensor resulting from isotropic thermal expansion with an initial strain due to the difference between the stress free temperature and initial temperature of the material.
8.111.4The tensor mechanics module shall have the capability to combine multiple eigenstrains to correctly calculate an eigenstrain tensor resulting from isotropic thermal expansion.
8.111.5The tensor mechanics module shall have the capability to calculate the eigenstrain tensor resulting from isotropic thermal expansion.
8.111.6The tensor mechanics module shall have the capability to calculate the eigenstrain tensor resulting from isotropic thermal expansion using the previous timestep.
8.111.7The Jacobian for the AD eigenstrain tensor resulting from isotropic thermal expansion shall be perfect
8.111.8The tensor mechanics module shall have the capability to calculate the eigenstrain tensor resulting from isotropic thermal expansion with an initial strain due to the difference between the stress free temperature and initial temperature of the material.
8.111.9The Jacobian for the AD eigenstrain tensor resulting from isotropic thermal expansion with an initial strain shall be perfect

solid_mechanics: Thermal Expansion Function
8.112.1
The system shall compute an eigenstrain due to thermal expansion using a function that describes a constant mean and instantaneous thermal expansion
1. using finite strain formulation
2. using small strain formulation
8.112.2
The system shall compute an eigenstrain due to thermal expansion using a function that describes a mean and instantaneous thermal expansion with a linear relationship to temperature
1. using finite strain formulation
2. using small strain formulation
8.112.3
The system shall compute an eigenstrain due and allow a smooth transition from negative to positive strain across the reference temperature and compare favorably to hand calculations
1. using a mean thermal expansion coefficient
2. using a mean thermal expansion coefficient with a scaling factor for sensitivity studies
3. using a instantaneous thermal expansion coefficient
4. using a dilatation thermal expansion coefficient

solid_mechanics: Torque
8.113.1The mechanics system shall provide a way to apply a torque to a boundary for small strain simulations.
8.113.2The mechanics system shall provide a way to apply a torque to a boundary for small strain simulations with automatic differentiation.

solid_mechanics: Torque Reaction
8.114.1The system shall compute the torque for an applied deformation at the end of a beam in 2D
8.114.2The system shall compute the torque for an applied deformation at the end of a beam in 3D
8.114.3The system shall compute the torque resulting from twisting applied using DisplacementAboutAxis to the end of an elastic cylinder
8.114.4The system shall provide a boundary condition that applies pure finite rigid body rotation to a surface of an elastic body
8.114.5The system shall provide a boundary condition that applies pure finite rigid body rotation to a surface of an elastic body incrementally
8.114.6The system shall be able to use disaplcement about axis after applying an arbitraty deformation
8.114.7The 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
8.114.8The system shall generate an error if the vector prescribed using the axis_direction parameter in the DisplacementAboutAxis boundary condition has a length of zero

solid_mechanics: Transfer From Displaced
8.115.1The system shall be able to transfer a field from a displaced source domain to a target domain, using shape function evaluations in the source domain.

solid_mechanics: Truss
8.116.1The mechanics system shall accurately model the axial response of 3D truss elements.
8.116.2The truss element shall work correctly when hex elements are also included in the same input or mesh file.
8.116.3The mechanics system shall accurately model the static response of a 2D frame modeled using truss elements.
8.116.4The LineElementAction shall correctly create the objects required for modeling the response of a mechanics system using 3D truss elements.
8.116.5The LineElementAction shall correctly create the objects required for modeling the response of a mechanics system using truss and hex elements.
8.116.6The LineElementAction shall correctly create the objects required for modeling the response of a mechanics system using 2D truss elements.
8.116.7The LineElementAction shall produce an error if area is not provided as input for truss elements.
8.116.8The LineElementAction shall produce an error if rotational variables are provided as input for truss elements.
8.116.9The system shall correctly model the plastic response of truss elements with a linear hardening model under tension.
8.116.10The system shall correctly model the plastic response of truss elements with perfect plasticity under tension.
8.116.11The system shall correctly model the plastic response of truss elements with a user-defined hardening function model under tension.
8.116.12The system shall correctly model the plastic response of truss elements with a linear hardening model under compression.
8.116.13The system shall correctly model the plastic response of truss elements with perfect plasticity under compression.
8.116.14The system shall correctly model the plastic response of truss elements with a user-defined hardening function model under compression.
8.116.15PlasticTruss material should produce error if neither the hardening constant nor a hardening function is provided.
8.116.16PlasticTruss material should produce error if both hardening constant and hardening function are provided.

solid_mechanics: Uel
8.117.1The UEL tri tests example shall be built.
8.117.2The UEL tri states tests example shall be built.
8.117.3The UEL build tests example shall be built.
8.117.4The UEL tri tests example shall be built.
8.117.5The UEL tri states tests example shall be built.
8.117.6The UEL build tests example shall be built.
8.117.7The system shall be able to solve a simple mechanics problem with small elastic deformation with a triangular mesh using MOOSE capabilities.
8.117.8The system shall be able to solve a simple mechanics problem with small elastic deformation with a triangular mesh using a UEL plugin.
8.117.9The system shall be able to solve a simple mechanics problem with small elastic deformation with a triangular mesh using a UEL plugin calling a UMAT routine with its standard interface.
8.117.10The Abaqus UEL interface shall be able to pass initial correct values of coordinates into a user-defined UEL routine for a hexahedral element and a non-trivial setup.
8.117.11The Abaqus UEL interface shall be able to interface with a UEL that calls a UMAT interface and must yield the same results for a triangular element as a UEL routine that computes internal forces without calling a UMAT routine.
8.117.12The Abaqus UMAT interface shall be able to yield the same results for a triangular element as a UEL routine that computes internal forces without calling a UMAT routine or an analogous UEL interface that calls a UMAT routine for computing the stress vector and the Jacobian.
8.117.13The Abaqus UEL interface, calling a UMAT routine, shall be able to yield the same results for a triangular element that a UMAT routine that computes internal forces when the solution depends on system states and two external fields, including temperature.
8.117.14The Abaqus UMAT interface shall be able to yield the same results for a triangular element that a UEL routine that computes internal forces, calling a UMAT routine, when the solution depends on system states and two external fields, including temperature.
8.117.15The Abaqus UMAT interface shall be able to yield the same results for a triangular element that a UEL routine that computes internal forces, calling a UMAT routine, when the solution depends on system states and two external fields, including temperature, which vary significantly with space.
8.117.16The Abaqus UEL interface shall be able to yield the same results for a triangular element that a UMAT routine that computes internal forces when the solution depends on system states and two external fields, including temperature, which vary significantly with space. Therefore, the system shall properly handle state and external fields whose values depend on element location.

solid_mechanics: Uexternaldb
8.118.1
The Abaqus UEXTERNALDB interface shall execute the UEXTERNALDB routine at the
1. beginning of the simulation.
2. beginning of each MOOSE time step.
3. end of each MOOSE time step.
4. end of the simulation.
8.118.2The Abaqus UEXTERNALDB interface shall interface with a general user object to obtain the proper loading step and pass it to the user routine.
8.118.3The Abaqus UEXTERNALDB interface shall provide the Abaqus API to obtain the simulation output directory.
8.118.4The Abaqus UEXTERNALDB interface shall provide the Abaqus API to obtain the simulation job name.
8.118.5
The Abaqus UEXTERNALDB interface shall provide Abaqus compatible message routine supporting the
1. info level.
2. warning level.
3. non-fatal error level.
4. fatal error level.
8.118.6
The Abaqus UEXTERNALDB interface shall provide Abaqus compatible mutex management
1. supporting single threaded execution.
2. supporting at least 2 threads.
3. supporting at least 4 threads.
8.118.7
The Abaqus UEXTERNALDB interface shall provide the Abaqus API for accessing shared memory objects
1. supporting single threaded execution.
2. supporting two threaded execution.
3. supporting four threaded execution.

solid_mechanics: Umat
8.119.1The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with support for stateful properties
8.119.2The system shall raise an error if a UMAT that requires incremental quantities is being used with a total strain formulation
8.119.3The system shall provide an interface to use Abaqus UMAT materials written in Fortran77 as constitutive models, with support for finite strain elastic material models
8.119.4The system shall provide an interface to use Abaqus UMAT materials written in Fortran77 as constitutive models, with support for small strain elastic material models
8.119.5The system shall provide an interface to use Abaqus UMAT materials written in C/C++ as constitutive models, with support for finite strain elastic material models
8.119.6The Abaqus UMAT interface shall produce the same results as the built-in MOOSE material models for finite strain elasticity
8.119.7The Abaqus UMAT interface shall produce the same results as the built-in MOOSE material models for small strain elasticity
8.119.8The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with support for incremental strain elastic material models
8.119.9The Abaqus UMAT interface shall produce the correct stresses for the hyperelasticNeo-Hookean model under shear loading. The results shall not depend on the time step size.Reference computation with fine time stepping.
8.119.10The Abaqus UMAT interface shall produce the correct stresses for the hyperelasticNeo-Hookean model under shear loading. The results shall not depend on the time step size.Comparison with coarse time stepping.
8.119.11The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with the ability to provide different parameters for different blocks using the same interface plug-in.
8.119.12The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with the ability to provide different parameters for different blocks using distinct interface plug-ins.
8.119.13The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with the ability to provide different parameters for different blocks using distinct interface plug-ins in multi processor runs.
8.119.14The system shall be able to simulate a multiblock solid mechanics test in which the blocks are meshed contiguously using the UMAT interface.
8.119.15The system shall be able to simulate a multiblock solid mechanics test in which the blocks are meshed contiguously.
8.119.16
8.119.17
8.119.18The Abaqus UMAT interface shall produce the same results (stress, strain, displacements) as the built-in MOOSE capabilities for a mechanical problem with an external field (not temperature) affecting material behavior
8.119.19The system shall avoid regression in a simple mechanical problem where a strain field modifies the stiffness of the material through CompositeElasticityTensor. This test also serves as a reference for UMAT external field verification
8.119.20The Abaqus UMAT interface shall produce the same results (stress, strain, displacements) as the built-in MOOSE capabilities for a mechanical problem with an external field (not temperature) step increment affecting material behavior
8.119.21The system shall avoid regression in a simple mechanical problem where a real number representing the strain increment modifies the stiffness of the material through CompositeElasticityTensor. This test also serves as a reference for UMAT external field increment verification
8.119.22The Abaqus UMAT interface shall produce the same results (stress, strain, displacements) as the built-in MOOSE capabilities for a mechanical problem with two external fields (not temperature) affecting material behavior
8.119.23The system shall avoid regression in a simple mechanical problem where two strain fields modifies the stiffness of the material through CompositeElasticityTensor. This test also serves as a reference for UMAT external field verification
8.119.24The Abaqus UMAT interface shall produce the same results (stress, strain, displacements) as the built-in MOOSE capabilities for a mechanical problem with two external material properties affecting material behavior
8.119.25The system shall avoid regression in a simple mechanical problem where two strain fields modifies the stiffness of the material through CompositeElasticityTensor. This test also serves as a reference for UMAT external material property verification
8.119.26The Abaqus UMAT interface shall pass correct values for STRAN, DSTRAN, TIME, CMNAME, NDI, NSHR, NTENS, COORDS, DROT, CELENT, DFGRD0, DFGRD1, NOEL, NPT, and KINC, when a single element is subjected to axial loading. Values are checked against verified references
8.119.27The Abaqus UMAT interface shall produce the same results (stress, strain, displacements) as the built-in MOOSE capabilities for a mechanical problem with two external fields (not temperature) step increment affecting material behavior and shear deformation
8.119.28The system shall avoid regression in a simple mechanical problem where two strain increments modify the stiffness of the material through CompositeElasticityTensor. This test also serves as a reference for UMAT behavior in the presence of shear deformation.
8.119.29The Abaqus UMAT interface shall pass the correct values for STRAN, DSTRAN, TIME, CMNAME, NDI, NSHR, NTENS, COORDS, DROT, CELENT, DFGRD0, DFGRD1, NOEL, NPT, and KINC, when shear deformation (causing nontrivial rotation kinematics) is present.
8.119.30The Abaqus UMAT interface shall pass the correct values for the deformation gradient DFGRD1 and the rotation increment DROT, accounting for the fact, that the calling C++ code is row majorand the called Fortran code is column major.
8.119.31The Abaqus UMAT interface shall pass the values for STRAN, DSTRAN, TIME, CMNAME, NDI, NSHR, NTENS, COORDS, DROT, CELENT, DFGRD0, DFGRD1, NOEL, NPT, and KINC, when shear deformation is present and large deformation kinematics is true, which captures variables or properties that may rely on the displaced mesh, such as COORDS and CELENT.
8.119.32The Abaqus UMAT interface shall pass correct values into a C UMAT routine for STRAN, DSTRAN, TIME, CMNAME, NDI, NSHR, NTENS, COORDS, DROT, CELENT, DFGRD0, DFGRD1, NOEL, NPT, and KINC, when a single element is subjected to axial loading. Values are checked against verified references
8.119.33The system shall generate identical numerical results regardless of whether the UMAT interface is written in Fortran or C. Generation of reference results.
8.119.34The system shall generate identical numerical results regardless of whether the UMAT interface is written in Fortran or C. Verification.
8.119.35The system shall allow for the correct use of Eigen matrices in C UMAT functions.
8.119.36UMAT shall generate the correct shear deformation and stress when a single element is sheared via a Dirichlet boundary condition on the XY plane with anisotropic shear stiffness; results must coincide with those of MOOSE.
8.119.37UMAT shall generate the correct shear deformation and stress when a single element is sheared via a Dirichlet boundary condition on the YZ plane with anisotropic shear stiffness; results must coincide with those of MOOSE.
8.119.38The system shall generate reference results of a sheared finite element on the XY plane.
8.119.39The system shall generate reference results of a sheared finite element on the YZ plane.
8.119.40The system shall generate a system Jacobian from the UMAT routines using large strain kinematics that is of the same order as the one generated by MOOSE using Rashid strain incrementation when the shear stiffness is anisotropic.
8.119.41The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with support for finite strain elastic material models and temperature dependence with the ability to switch boundary conditions according to the controls system.
8.119.42The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with support for finite strain elastic material models and temperature dependence with the ability to switch boundary conditions according to the controls system using loading steps from a general user object.
8.119.43The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with support for finite strain elastic material models and temperature dependence with the ability to switch boundary conditions according to the controls system using loading steps from a general user object, to which the user only provides the step durations.
8.119.44The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with support for finite strain elastic material models and temperature dependence with the ability to switch boundary conditions according to the controls system using loading steps from a general user object, to which the user only provides the number of steps and the total time interval.
8.119.45The system shall ensure that solution synchronization times are enforced by default when employing a user object to determine the steps of the simulation.
8.119.46The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with support for finite strain elastic material models and temperature dependence with the ability to switch boundary conditions according to the controls system using loading steps from a general user object.
8.119.47The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with support for finite strain elastic material models and temperature dependence
8.119.48The Abaqus UMAT interface shall produce the same results as the built-in MOOSE material models for finite strain elasticity and temperature dependence
8.119.49The system shall provide an interface to use Abaqus UMAT materials as constitutive models, with support for finite strain elastic material models and temperature dependence, include the temperature step increment
8.119.50The Abaqus UMAT interface shall produce the same results as the built-in MOOSE material models for finite strain elasticity and temperature dependence, including the temperature step increment
8.119.51The system shall provide an interface to use Abaqus UMAT materials that allows the user to control the time step increment within the UMAT routine and combine it with MOOSE native time step controls, including cutback and growth factors, and soft terminations.

solid_mechanics: Visco
8.120.1
8.120.2
8.120.3
8.120.4
8.120.5
8.120.6
8.120.7

solid_mechanics: Volumetric Deform Grad
8.121.1The ComputeDeformGradBasedStress class shall correctly compute the stress based on the lagrangian strain.
8.121.2The ComputeDeformGradBasedStress class shall correctly compute the stress from lagrangian strain when volumetric locking correction is used.
8.121.3The ComputeVolumeDeformGrad and the VolumeDeformGradCorrectedStress classes shall correctly compute the volumetric deformation gradient, total deformation gradient and transform the stress from previous configuration to the current configuration.
8.121.4The ComputeVolumeDeformGrad and the VolumeDeformGradCorrectedStress classes shall correctly compute the volumetric deformation gradient, total deformation gradient and transform the stress from previous configuration to the current configuration when volumetric locking correction is used.

solid_mechanics: Volumetric Eigenstrain
8.122.1The system shall correctly compute an eigenstrain tensor that results in a solution that exactly recovers the specified volumetric expansion, and the reported volumetric strain computed by RankTwoScalarAux shall match the prescribed volumetric strain.
8.122.2The system shall correctly compute an eigenstrain using AD tensor that results in a solution that exactly recovers the specified volumetric expansion, and the reported volumetric strain computed by RankTwoScalarAux shall match the prescribed volumetric strain.
8.122.3The volumetric strain computed using RankTwoScalarAux for a unit cube with imposed displacements shall be identical to that obtained by imposing an eigenstrain that causes the same deformation of that model.

solid_mechanics: Volumetric Locking Verification
8.123.1The mechanics system shall correctly model the deformation of a 2D membrane with nearly incompressible material when volumetric locking correction is set to true.
8.123.2The mechanics system shall correctly model the locking behavior of a 2D membrane with nearly incompressible material when volumetric locking correction is set to false.

solid_mechanics: Weak Plane Shear
8.124.1The system shall not allow a weak plane shear friction angle less than the dilation angle.
8.124.2The system shall not allow a zero length transverse direction vector
8.124.3The 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
8.124.4
For the case of small deformations the system shall
1. return a value of shear of 1 and a normal stress value of 0 when the apex of the shear envelope is (0,2);
2. utlize the small smoother parameter under both pure tension and pure shear to solve;
3. return the stress state to the yield surface under random loading;
4. utilize the cap smoothing functionality under both pure tension and pure shear to solve the simulation.
8.124.5
For large deformation cases the system shall
1. produce a solution without plasticity after the mesh is rotated by 90 degrees;
2. return the stress state to the yield surface under pure tensile loading conditions after the mesh is rotated by 45 degrees;
3. return the stress state to the yield surface under pseudo-random loading directions;
4. return the stress state to the yield surface under pseudo-random loading while using the cap-tip yeild surface scheme.
8.124.6
Under small deformation loading selected to induce cohesion hardening the system shall
1. with repeated stretch loading;
2. for pure tension and then pure shear to map out the yield function;
3. with linearly increasing loading and an exponential hardening model;
4. with a linearly increasing load and a cubic hardening model.
8.124.7The system shall return the stress state to the yield surface under random large deformation loads.

solid_mechanics: Weak Plane Tensile
8.125.1The system shall not allow a negative tensile strength value.
8.125.2
With a single element stretched by 1e-6m, for small deformation, the system shall
1. return the stress state to the yield surface when the applied stress is above the yield surface tensile strength;
2. return the stress state to the specified vector yield surface;
3. not calculate a plastic solution because the applied stress is lower than the tension cutoff for the yield surface.
8.125.3
Under large deformation loading the system shall
1. calculate a solution without plasticity after the mesh is rotated by 90 degrees;
2. return the stress state to the yield surface under pure tensile loading after the mesh is rotated by 45 degrees.
8.125.4
Under small deformation loading selected to induce hardening the system shall
1. the deformation shall produce increments of plastic strain in response to incremental loading steps;
2. evolve the plastic internal state variable for an exponential hardening model;
3. return the stress state to the evolving yield surface when loaded solely in the z-direction for an exponential hardening model;
4. return the stress state to the evolving yield surface when loaded solely in the z-direction for an cubic hardening model;

Usability Requirements

No requirements of this type exist for this application, beyond those of its dependencies.

Performance Requirements

No requirements of this type exist for this application, beyond those of its dependencies.

System Interfaces

solid_mechanics: Neml2
8.84.1The system shall be capable of capturing errors occurring during the evaluation of a NEML2 model.
8.84.2The system shall throw an error upon instantiation of NEML2 interface objects if NEML2 is not available.
8.84.3
The system shall be capable of modularly interfacing with the NEML2 library to perform constitutive update
1. to solve an elasticity problem;
2. to solve a perfect viscoplasticity problem;
3. to solve a viscoplasticity problem with isotropic hardening;
4. to solve a viscoplasticity problem with kinematic hardening;
5. to solve a viscoplasticity problem with both isotropic and kinematic hardening;
6. to solve a viscoplasticity problem with non-associative Chaboche hardening;
7. to solve a viscoplasticity problem using radial return;
8. to solve a rate-independent problem with perfect plasticity;
9. to solve a rate-independent problem with isotropic hardening;
10. to solve a rate-independent problem with kinematic hardening;
11. to solve a rate-independent problem with both isotropic and kinematic hardening.
8.84.8
The system shall be capable of interfacing with the NEML2 library to perform constitutive update
1. to solve an elasticity problem;
2. to solve a perfect viscoplasticity problem;
3. to solve a viscoplasticity problem with isotropic hardening;
4. to solve a viscoplasticity problem with kinematic hardening;
5. to solve a viscoplasticity problem with both isotropic and kinematic hardening;
6. to solve a viscoplasticity problem with non-associative Chaboche hardening;
7. to solve a viscoplasticity problem using radial return;
8. to solve a rate-independent problem with perfect plasticity;
9. to solve a rate-independent problem with isotropic hardening;
10. to solve a rate-independent problem with kinematic hardening;
11. to solve a rate-independent problem with both isotropic and kinematic hardening.
12. in a batched fashion.

System Operations

Human System Integration Requirements

The Solid Mechanics module is command line driven and conforms to all standard terminal behaviors. Specific human system interaction accommodations shall be a function of the end-user's terminal. MOOSE (and therefore the Solid Mechanics module) does support optional coloring within the terminal's ability to display color, which may be disabled.

Maintainability

The latest working version (defined as the version that passes all tests in the current regression test suite) shall be publicly available at all times through the repository host provider.
Flaws identified in the system shall be reported and tracked in a ticket or issue based system. The technical lead will determine the severity and priority of all reported issues and assign resources at their discretion to resolve identified issues.
The software maintainers will entertain all proposed changes to the system in a timely manner (within two business days).
The core software in its entirety will be made available under the terms of a designated software license. These license terms are outlined in the LICENSE file alongside the Solid Mechanics module source code. As a MOOSE physics module, the license for the Solid Mechanics module is identical to that of the framework - that is, the GNU Lesser General Public License (LGPL) version 2.1 license.

Reliability

The regression test suite will cover at least 86% of all lines of code within the Solid Mechanics module at all times. Known regressions will be recorded and tracked (see Maintainability) to an independent and satisfactory resolution.

System Modes and States

MOOSE applications normally run in normal execution mode when an input file is supplied. However, there are a few other modes that can be triggered with various command line flags as indicated here:

Command Line Flag	Description of mode
`-i <input_file>`	Normal execution mode
`--split-mesh <splits>`	Read the mesh block splitting the mesh into two or more pieces for use in a subsequent run
`--use-split`	(implies -i flag) Execute the simulation but use pre-split mesh files instead of the mesh from the input file
`--yaml`	Output all object descriptions and available parameters in YAML format
`--json`	Output all object descriptions and available parameters in JSON format
`--syntax`	Output all registered syntax
`--registry`	Output all known objects and actions
`--registry-hit`	Output all known objects and actions in HIT format
`--mesh-only` (implies -i flag)	Run only the mesh related tasks and output the final mesh that would be used for the simulation
`--start-in-debugger <debugger>`	Start the simulation attached to the supplied debugger

note

The list of system-modes may not be extensive as the system is designed to be extendable to end-user applications. The complete list of command line options for applications can be obtained by running the executable with zero arguments. See the command line usage.

Physical Characteristics

The Solid Mechanics module is software only with no associated physical media. See System Requirements for a description of the minimum required hardware necessary for running the Solid Mechanics module.

Environmental Conditions

Not Applicable

System Security

MOOSE-based applications such as the Solid Mechanics module have no requirements or special needs related to system security. The software is designed to run completely in user-space with no elevated privileges required nor recommended.

Information Management

The core framework and all modules in their entirety will be made publicly available on an appropriate repository hosting site. Day-to-day backups and security services will be provided by the hosting service. More information about MOOSE backups of the public repository on INL-hosted services can be found on the following page: GitHub Backups

Polices and Regulations

MOOSE-based applications must comply with all export control restrictions.

System Life Cycle Sustainment

MOOSE-based development follows various agile methods. The system is continuously built and deployed in a piecemeal fashion since objects within the system are more or less independent. Every new object requires a test, which in turn requires an associated requirement and design description. The Solid Mechanics module development team follows the NQA-1 standards.

Packaging, Handling, Shipping and Transportation

No special requirements are needed for packaging or shipping any media containing MOOSE and Solid Mechanics module source code. However, some MOOSE-based applications that use the Solid Mechanics module may be export-controlled, in which case all export control restrictions must be adhered to when packaging and shipping media.

Verification

The regression test suite will employ several verification tests using comparison against known analytical solutions, the method of manufactured solutions, and convergence rate analysis.

Introduction
System Overview
References
Definitions and Acronyms
System Requirements
Verification