diff --git a/examples/Validation/beam_extension_compression/beam.xml b/examples/Validation/beam_extension_compression/beam.xml
new file mode 100644
index 00000000000..76d18779ce6
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/beam.xml
@@ -0,0 +1,41 @@
+<Node name="Group">
+
+    <Node name="plugins">
+        <RequiredPlugin pluginName="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedProjectiveConstraint] -->
+        <RequiredPlugin pluginName="Sofa.Component.Engine.Select"/> <!-- Needed to use components [BoxROI] -->
+        <RequiredPlugin pluginName="Sofa.Component.LinearSolver.Direct"/> <!-- Needed to use components [SparseLDLSolver] -->
+        <RequiredPlugin pluginName="Sofa.Component.LinearSolver.Iterative"/> <!-- Needed to use components [CGLinearSolver] -->
+        <RequiredPlugin pluginName="Sofa.Component.LinearSolver.Ordering"/> <!-- Needed to use components [NaturalOrderingMethod] -->
+        <RequiredPlugin pluginName="Sofa.Component.LinearSystem"/> <!-- Needed to use components [ConstantSparsityPatternSystem] -->
+        <RequiredPlugin pluginName="Sofa.Component.Mass"/> <!-- Needed to use components [MeshMatrixMass] -->
+        <RequiredPlugin pluginName="Sofa.Component.ODESolver.Backward"/> <!-- Needed to use components [EulerImplicitSolver] -->
+        <RequiredPlugin pluginName="Sofa.Component.SolidMechanics.FEM.Elastic"/>
+        <RequiredPlugin pluginName="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
+        <RequiredPlugin pluginName="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [TetrahedronSetGeometryAlgorithms,TetrahedronSetTopologyContainer,TetrahedronSetTopologyModifier] -->
+        <RequiredPlugin pluginName="Sofa.Component.Topology.Container.Grid"/> <!-- Needed to use components [RegularGridTopology] -->
+        <RequiredPlugin pluginName="Sofa.Component.Topology.Mapping"/> <!-- Needed to use components [Hexa2TetraTopologicalMapping] -->
+        <RequiredPlugin pluginName="Sofa.Component.Visual"/> <!-- Needed to use components [LineAxis,VisualGrid,VisualStyle] -->
+        <RequiredPlugin pluginName="Sofa.GL.Component.Rendering3D"/> <!-- Needed to use components [OglSceneFrame] -->
+    </Node>
+
+    <DefaultAnimationLoop/>
+    <VisualStyle displayFlags="hideBehaviorModels showForceFields" />
+
+    <VisualGrid size="0.1"/>
+    <LineAxis size="0.1"/>
+
+    <EulerImplicitSolver name="backward_Euler" rayleighStiffness="0.1" rayleighMass="0.1" />
+
+    <RegularGridTopology name="grid" min="-0.05 -0.01 -0.01" max="0.05 0.01 0.01" n="20 5 5"/>
+    <MechanicalObject template="Vec3" name="state" showObject="true"/>
+
+    <NodalMassDensity property="1100"/>
+    <FEMMass template="Vec3,Hexahedron"/>
+
+    <BoxROI template="Vec3" name="box_roi_left" box="-0.051 -0.011 -0.011   -0.049 0.011 0.011" drawBoxes="1" />
+    <LinearMovementProjectiveConstraint template="Vec3" keyTimes="0 1 3 4 1e12" movements="0 0 0  -0.02 0 0  0.02 0 0  0 0 0  0 0 0" indices="@box_roi_left.indices" />
+
+    <BoxROI template="Vec3" name="box_roi_right" box="0.049 -0.011 -0.011   0.051 0.011 0.011" drawBoxes="1" />
+    <LinearMovementProjectiveConstraint template="Vec3" keyTimes="0 1 3 4 1e12" movements="0 0 0  0.02 0 0  -0.02 0 0  0 0 0  0 0 0" indices="@box_roi_right.indices" />
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/hexahedron/corotational/parallel/assembled/CorotationalFEMForceField.scn b/examples/Validation/beam_extension_compression/hexahedron/corotational/parallel/assembled/CorotationalFEMForceField.scn
new file mode 100644
index 00000000000..9b530ebb9b2
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/hexahedron/corotational/parallel/assembled/CorotationalFEMForceField.scn
@@ -0,0 +1,15 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
+    <NaturalOrderingMethod/>
+    <SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>
+
+    <CorotationalFEMForceField name="FEM" template="Vec3,Hexahedron"
+                               youngModulus="2e6" poissonRatio="0.45" topology="@grid"
+                               computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
+    <VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
+                    colorMap="green yellow orange red purple #221C35"/>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/hexahedron/corotational/parallel/matrixfree/CorotationalFEMForceField.scn b/examples/Validation/beam_extension_compression/hexahedron/corotational/parallel/matrixfree/CorotationalFEMForceField.scn
new file mode 100644
index 00000000000..373a72cb45c
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/hexahedron/corotational/parallel/matrixfree/CorotationalFEMForceField.scn
@@ -0,0 +1,13 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />
+
+    <CorotationalFEMForceField name="FEM" template="Vec3,Hexahedron"
+                               youngModulus="2e6" poissonRatio="0.45" topology="@grid"
+                               computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
+    <VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
+                    colorMap="green yellow orange red purple #221C35"/>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/hexahedron/corotational/sequential/assembled/CorotationalFEMForceField.scn b/examples/Validation/beam_extension_compression/hexahedron/corotational/sequential/assembled/CorotationalFEMForceField.scn
new file mode 100644
index 00000000000..90d6513a214
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/hexahedron/corotational/sequential/assembled/CorotationalFEMForceField.scn
@@ -0,0 +1,15 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
+    <NaturalOrderingMethod/>
+    <SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>
+
+    <CorotationalFEMForceField name="FEM" template="Vec3,Hexahedron"
+                               youngModulus="2e6" poissonRatio="0.45" topology="@grid"
+                               computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
+    <VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
+                    colorMap="green yellow orange red purple #221C35"/>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/hexahedron/corotational/sequential/matrixfree/CorotationalFEMForceField.scn b/examples/Validation/beam_extension_compression/hexahedron/corotational/sequential/matrixfree/CorotationalFEMForceField.scn
new file mode 100644
index 00000000000..5dc05bb8c77
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/hexahedron/corotational/sequential/matrixfree/CorotationalFEMForceField.scn
@@ -0,0 +1,13 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />
+
+    <CorotationalFEMForceField name="FEM" template="Vec3,Hexahedron"
+                               youngModulus="2e6" poissonRatio="0.45" topology="@grid"
+                               computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
+    <VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
+                    colorMap="green yellow orange red purple #221C35"/>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/hexahedron/linear/parallel/assembled/LinearSmallStrainFEMForceField.scn b/examples/Validation/beam_extension_compression/hexahedron/linear/parallel/assembled/LinearSmallStrainFEMForceField.scn
new file mode 100644
index 00000000000..3f7dbb1b531
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/hexahedron/linear/parallel/assembled/LinearSmallStrainFEMForceField.scn
@@ -0,0 +1,15 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
+    <NaturalOrderingMethod/>
+    <SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>
+
+    <LinearSmallStrainFEMForceField name="FEM" template="Vec3,Hexahedron"
+                                    youngModulus="2e6" poissonRatio="0.45" topology="@grid"
+                                    computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
+    <VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
+                    colorMap="green yellow orange red purple #221C35"/>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/hexahedron/linear/parallel/matrixfree/LinearSmallStrainFEMForceField.scn b/examples/Validation/beam_extension_compression/hexahedron/linear/parallel/matrixfree/LinearSmallStrainFEMForceField.scn
new file mode 100644
index 00000000000..af5d15c64ee
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/hexahedron/linear/parallel/matrixfree/LinearSmallStrainFEMForceField.scn
@@ -0,0 +1,13 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />
+
+    <LinearSmallStrainFEMForceField name="FEM" template="Vec3,Hexahedron"
+                                           youngModulus="2e6" poissonRatio="0.45" topology="@grid"
+                                           computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
+    <VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
+                    colorMap="green yellow orange red purple #221C35"/>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/hexahedron/linear/sequential/assembled/LinearSmallStrainFEMForceField.scn b/examples/Validation/beam_extension_compression/hexahedron/linear/sequential/assembled/LinearSmallStrainFEMForceField.scn
new file mode 100644
index 00000000000..ed37723b776
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/hexahedron/linear/sequential/assembled/LinearSmallStrainFEMForceField.scn
@@ -0,0 +1,15 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
+    <NaturalOrderingMethod/>
+    <SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>
+
+    <LinearSmallStrainFEMForceField name="FEM" template="Vec3,Hexahedron"
+                                           youngModulus="2e6" poissonRatio="0.45" topology="@grid"
+                                           computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
+    <VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
+                    colorMap="green yellow orange red purple #221C35"/>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/hexahedron/linear/sequential/matrixfree/LinearSmallStrainFEMForceField.scn b/examples/Validation/beam_extension_compression/hexahedron/linear/sequential/matrixfree/LinearSmallStrainFEMForceField.scn
new file mode 100644
index 00000000000..dd108fd4f4c
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/hexahedron/linear/sequential/matrixfree/LinearSmallStrainFEMForceField.scn
@@ -0,0 +1,13 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />
+
+    <LinearSmallStrainFEMForceField name="FEM" template="Vec3,Hexahedron"
+                                           youngModulus="2e6" poissonRatio="0.45" topology="@grid"
+                                           computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
+    <VonMisesStress template="Vec3,Hexahedron" name="stress" topology="@grid" stressEvaluator="@FEM"
+                    colorMap="green yellow orange red purple #221C35"/>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/tetrahedron/corotational/parallel/assembled/CorotationalFEMForceField.scn b/examples/Validation/beam_extension_compression/tetrahedron/corotational/parallel/assembled/CorotationalFEMForceField.scn
new file mode 100644
index 00000000000..ac3ee5c0d14
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/tetrahedron/corotational/parallel/assembled/CorotationalFEMForceField.scn
@@ -0,0 +1,22 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
+    <NaturalOrderingMethod/>
+    <SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>
+
+    <Node name="fem">
+        <TetrahedronSetTopologyContainer name="Tetra_topo"/>
+        <TetrahedronSetTopologyModifier name="Modifier" />
+        <TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
+        <Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>
+
+        <CorotationalFEMForceField name="FEM" template="Vec3,Tetrahedron"
+                                   youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
+                                   computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
+        <VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
+                        colorMap="green yellow orange red purple #221C35"/>
+    </Node>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/tetrahedron/corotational/parallel/matrixfree/CorotationalFEMForceField.scn b/examples/Validation/beam_extension_compression/tetrahedron/corotational/parallel/matrixfree/CorotationalFEMForceField.scn
new file mode 100644
index 00000000000..8fac16de66f
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/tetrahedron/corotational/parallel/matrixfree/CorotationalFEMForceField.scn
@@ -0,0 +1,20 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />
+
+    <Node name="fem">
+        <TetrahedronSetTopologyContainer name="Tetra_topo"/>
+        <TetrahedronSetTopologyModifier name="Modifier" />
+        <TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
+        <Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>
+
+        <CorotationalFEMForceField name="FEM" template="Vec3,Tetrahedron"
+                                   youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
+                                   computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
+        <VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
+                        colorMap="green yellow orange red purple #221C35"/>
+    </Node>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/tetrahedron/corotational/sequential/assembled/CorotationalFEMForceField.scn b/examples/Validation/beam_extension_compression/tetrahedron/corotational/sequential/assembled/CorotationalFEMForceField.scn
new file mode 100644
index 00000000000..673671ecffc
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/tetrahedron/corotational/sequential/assembled/CorotationalFEMForceField.scn
@@ -0,0 +1,22 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
+    <NaturalOrderingMethod/>
+    <SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>
+
+    <Node name="fem">
+        <TetrahedronSetTopologyContainer name="Tetra_topo"/>
+        <TetrahedronSetTopologyModifier name="Modifier" />
+        <TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
+        <Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>
+
+        <CorotationalFEMForceField name="FEM" template="Vec3,Tetrahedron"
+                                   youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
+                                   computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
+        <VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
+                        colorMap="green yellow orange red purple #221C35"/>
+    </Node>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/tetrahedron/corotational/sequential/matrixfree/CorotationalFEMForceField.scn b/examples/Validation/beam_extension_compression/tetrahedron/corotational/sequential/matrixfree/CorotationalFEMForceField.scn
new file mode 100644
index 00000000000..ef839934ced
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/tetrahedron/corotational/sequential/matrixfree/CorotationalFEMForceField.scn
@@ -0,0 +1,20 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />
+
+    <Node name="fem">
+        <TetrahedronSetTopologyContainer name="Tetra_topo"/>
+        <TetrahedronSetTopologyModifier name="Modifier" />
+        <TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
+        <Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>
+
+        <CorotationalFEMForceField name="FEM" template="Vec3,Tetrahedron"
+                                   youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
+                                   computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
+        <VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
+                        colorMap="green yellow orange red purple #221C35"/>
+    </Node>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/tetrahedron/linear/parallel/assembled/LinearSmallStrainFEMForceField.scn b/examples/Validation/beam_extension_compression/tetrahedron/linear/parallel/assembled/LinearSmallStrainFEMForceField.scn
new file mode 100644
index 00000000000..c67ccfb20e4
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/tetrahedron/linear/parallel/assembled/LinearSmallStrainFEMForceField.scn
@@ -0,0 +1,22 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
+    <NaturalOrderingMethod/>
+    <SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>
+
+    <Node name="fem">
+        <TetrahedronSetTopologyContainer name="Tetra_topo"/>
+        <TetrahedronSetTopologyModifier name="Modifier" />
+        <TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
+        <Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>
+
+        <LinearSmallStrainFEMForceField name="FEM" template="Vec3,Tetrahedron"
+                                               youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
+                                               computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
+        <VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
+                        colorMap="green yellow orange red purple #221C35"/>
+    </Node>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/tetrahedron/linear/parallel/matrixfree/LinearSmallStrainFEMForceField.scn b/examples/Validation/beam_extension_compression/tetrahedron/linear/parallel/matrixfree/LinearSmallStrainFEMForceField.scn
new file mode 100644
index 00000000000..d376254908b
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/tetrahedron/linear/parallel/matrixfree/LinearSmallStrainFEMForceField.scn
@@ -0,0 +1,20 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />
+
+    <Node name="fem">
+        <TetrahedronSetTopologyContainer name="Tetra_topo"/>
+        <TetrahedronSetTopologyModifier name="Modifier" />
+        <TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
+        <Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>
+
+        <LinearSmallStrainFEMForceField name="FEM" template="Vec3,Tetrahedron"
+                                               youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
+                                               computeForceStrategy="parallel" computeForceDerivStrategy="parallel"/>
+        <VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
+                        colorMap="green yellow orange red purple #221C35"/>
+    </Node>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/tetrahedron/linear/sequential/assembled/LinearSmallStrainFEMForceField.scn b/examples/Validation/beam_extension_compression/tetrahedron/linear/sequential/assembled/LinearSmallStrainFEMForceField.scn
new file mode 100644
index 00000000000..45ea1d2071f
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/tetrahedron/linear/sequential/assembled/LinearSmallStrainFEMForceField.scn
@@ -0,0 +1,22 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <ConstantSparsityPatternSystem template="CompressedRowSparseMatrix" name="A" checkIndices="false"/>
+    <NaturalOrderingMethod/>
+    <SparseLDLSolver name="linear_solver" template="CompressedRowSparseMatrix"/>
+
+    <Node name="fem">
+        <TetrahedronSetTopologyContainer name="Tetra_topo"/>
+        <TetrahedronSetTopologyModifier name="Modifier" />
+        <TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
+        <Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>
+
+        <LinearSmallStrainFEMForceField name="FEM" template="Vec3,Tetrahedron"
+                                               youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
+                                               computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
+        <VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
+                        colorMap="green yellow orange red purple #221C35"/>
+    </Node>
+
+</Node>
diff --git a/examples/Validation/beam_extension_compression/tetrahedron/linear/sequential/matrixfree/LinearSmallStrainFEMForceField.scn b/examples/Validation/beam_extension_compression/tetrahedron/linear/sequential/matrixfree/LinearSmallStrainFEMForceField.scn
new file mode 100644
index 00000000000..5549770ea73
--- /dev/null
+++ b/examples/Validation/beam_extension_compression/tetrahedron/linear/sequential/matrixfree/LinearSmallStrainFEMForceField.scn
@@ -0,0 +1,20 @@
+<?xml version="1.0"?>
+<Node name="root" dt="0.01" gravity="0 0 0">
+    <include href="../../../../beam.xml"/>
+
+    <CGLinearSolver iterations="250" name="linear_solver" tolerance="1.0e-12" threshold="1.0e-12" />
+
+    <Node name="fem">
+        <TetrahedronSetTopologyContainer name="Tetra_topo"/>
+        <TetrahedronSetTopologyModifier name="Modifier" />
+        <TetrahedronSetGeometryAlgorithms template="Vec3" name="GeomAlgo" drawTetrahedra="false"/>
+        <Hexa2TetraTopologicalMapping input="@grid" output="@Tetra_topo" swapping="true"/>
+
+        <LinearSmallStrainFEMForceField name="FEM" template="Vec3,Tetrahedron"
+                                               youngModulus="2e6" poissonRatio="0.45" topology="@Tetra_topo"
+                                               computeForceStrategy="sequenced" computeForceDerivStrategy="sequenced"/>
+        <VonMisesStress template="Vec3,Tetrahedron" name="stress" topology="@Tetra_topo" stressEvaluator="@FEM"
+                        colorMap="green yellow orange red purple #221C35"/>
+    </Node>
+
+</Node>