diff --git a/ProcessLib/LIE/SmallDeformation/Tests.cmake b/ProcessLib/LIE/SmallDeformation/Tests.cmake
index 1120e08a0664cd97f76c4f0489f60cb5ebdc986a..5ecdf5b5b003f1d6f3e6e224cd355b5b1e3a98f2 100644
--- a/ProcessLib/LIE/SmallDeformation/Tests.cmake
+++ b/ProcessLib/LIE/SmallDeformation/Tests.cmake
@@ -79,3 +79,20 @@ AddTest(
     expected_single_joint_displacement_controlled_pcs_0_ts_10_t_1.000000.vtu single_joint_displacement_controlled_pcs_0_ts_10_t_1.000000.vtu displacement_jump1 displacement_jump1 1e-16 1e-16
     expected_single_joint_displacement_controlled_pcs_0_ts_10_t_1.000000.vtu single_joint_displacement_controlled_pcs_0_ts_10_t_1.000000.vtu aperture aperture 1e-16 1e-16
 )
+
+AddTest(
+    NAME LIE_M_cohesive_zone_modeI
+    PATH LIE/Mechanics
+    EXECUTABLE ogs
+    EXECUTABLE_ARGS cohesive_zone_load_path.prj
+    WRAPPER time
+    TESTER vtkdiff
+    REQUIREMENTS NOT OGS_USE_MPI
+    DIFF_DATA
+    GLOB cohesive_zone_load_path_pcs_0_ts_*.vtu displacement displacement 1e-16 0
+    GLOB cohesive_zone_load_path_pcs_0_ts_*.vtu displacement_jump1 displacement_jump1 1e-16 0
+    GLOB cohesive_zone_load_path_pcs_0_ts_*.vtu sigma_xx sigma_xx 1e-6 0
+    GLOB cohesive_zone_load_path_pcs_0_ts_*.vtu sigma_yy sigma_yy 1e-6 0
+    GLOB cohesive_zone_load_path_pcs_0_ts_*.vtu sigma_zz sigma_zz 1e-6 0
+    GLOB cohesive_zone_load_path_pcs_0_ts_*.vtu sigma_xy sigma_xy 1e-6 0
+)
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.gml b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.gml
new file mode 100644
index 0000000000000000000000000000000000000000..3df1a09712144daec6fdc0745d4460736968d4fd
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:a208a8c803f6b4a0e80074418247a1ff327260261742cd059fc164b082150e52
+size 890
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.prj b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.prj
new file mode 100644
index 0000000000000000000000000000000000000000..e0970482970da5a29a9f2d10916076f008aae635
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.prj
@@ -0,0 +1,219 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <mesh>cohesive_zone_load_path.vtu</mesh>
+    <geometry>cohesive_zone_load_path.gml</geometry>
+    <processes>
+        <process>
+            <name>SD</name>
+            <type>SMALL_DEFORMATION_WITH_LIE</type>
+            <integration_order>2</integration_order>
+            <dimension>2</dimension>
+            <process_variables>
+                <process_variable>displacement</process_variable>
+                <process_variable>displacement_jump1</process_variable>
+            </process_variables>
+            <constitutive_relation>
+                <type>LinearElasticIsotropic</type>
+                <youngs_modulus>E</youngs_modulus>
+                <poissons_ratio>nu</poissons_ratio>
+            </constitutive_relation>
+            <fracture_model>
+                <type>CohesiveZoneModeI</type>
+                <normal_stiffness>Kn</normal_stiffness>
+                <shear_stiffness>Ks</shear_stiffness>
+                <fracture_toughness>Gc</fracture_toughness>
+                <peak_normal_traction>t_np</peak_normal_traction>
+                <penalty_aperture_cutoff>4e-6</penalty_aperture_cutoff>
+                <tension_cutoff>1</tension_cutoff>
+            </fracture_model>
+            <fracture_properties>
+                <material_id>0</material_id>
+                <initial_aperture>aperture0</initial_aperture>
+            </fracture_properties>
+            <secondary_variables>
+                <secondary_variable type="static" internal_name="sigma_xx" output_name="sigma_xx"/>
+                <secondary_variable type="static" internal_name="sigma_yy" output_name="sigma_yy"/>
+                <secondary_variable type="static" internal_name="sigma_zz" output_name="sigma_zz"/>
+                <secondary_variable type="static" internal_name="sigma_xy" output_name="sigma_xy"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="SD">
+                <nonlinear_solver>basic_newton</nonlinear_solver>
+                <convergence_criterion>
+                    <type>DeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <abstol>1e-14</abstol>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <output>
+                    <variables>
+                        <variable>displacement</variable>
+                        <variable>displacement_jump1</variable>
+                        <variable>sigma_xx</variable>
+                        <variable>sigma_yy</variable>
+                        <variable>sigma_zz</variable>
+                        <variable>sigma_xy</variable>
+                    </variables>
+                </output>
+                <time_stepping>
+                    <type>FixedTimeStepping</type>
+                    <t_initial>0</t_initial>
+                    <t_end>100</t_end>
+                    <timesteps>
+                        <pair>
+                            <repeat>100</repeat>
+                            <delta_t>1.0</delta_t>
+                        </pair>
+                    </timesteps>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>cohesive_zone_load_path</prefix>
+            <timesteps>
+                <pair>
+                    <repeat>10</repeat>
+                    <each_steps>10</each_steps>
+                </pair>
+            </timesteps>
+        </output>
+    </time_loop>
+    <parameters>
+        <parameter>
+            <name>E</name>
+            <type>Constant</type>
+            <value>1e14</value>
+        </parameter>
+        <parameter>
+            <name>nu</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>Kn</name>
+            <type>Constant</type>
+            <value>1e12</value>
+        </parameter>
+        <parameter>
+            <name>Ks</name>
+            <type>Constant</type>
+            <value>1e12</value>
+        </parameter>
+        <parameter>
+            <name>Gc</name>
+            <type>Constant</type>
+            <value>200</value>
+        </parameter>
+        <parameter>
+            <name>t_np</name>
+            <type>Constant</type>
+            <value>1e7</value>
+        </parameter>
+        <parameter>
+            <name>displacement0</name>
+            <type>Constant</type>
+            <values>0 0</values>
+        </parameter>
+        <parameter>
+            <name>aperture0</name>
+            <type>Constant</type>
+            <value>1e-5</value>
+        </parameter>
+        <parameter>
+            <name>displacement_unit_top</name>
+            <type>Constant</type>
+            <value>5e-7</value>
+        </parameter>
+        <parameter>
+            <name>displacement_unit_bottom</name>
+            <type>Constant</type>
+            <value>-5e-7</value>
+        </parameter>
+        <parameter>
+            <name>displacement_top</name>
+            <type>CurveScaled</type>
+            <curve>timeRamp_displacement</curve>
+            <parameter>displacement_unit_top</parameter>
+        </parameter>
+        <parameter>
+            <name>displacement_bottom</name>
+            <type>CurveScaled</type>
+            <curve>timeRamp_displacement</curve>
+            <parameter>displacement_unit_bottom</parameter>
+        </parameter>
+        <parameter>
+            <name>zero_u</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+    </parameters>
+    <curves>
+        <curve>
+            <name>timeRamp_displacement</name>
+            <coords>0 20 40 50 60 80 100</coords>
+            <values>0 20  0 -8  0 20  40</values>
+        </curve>
+    </curves>
+    <process_variables>
+        <process_variable>
+            <name>displacement</name>
+            <components>2</components>
+            <order>1</order>
+            <initial_condition>displacement0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>TwoElementsWithInterface</geometrical_set>
+                    <geometry>PLY_TOP</geometry>
+                    <type>Dirichlet</type>
+                    <component>1</component>
+                    <parameter>displacement_top</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <geometrical_set>TwoElementsWithInterface</geometrical_set>
+                    <geometry>PLY_BOTTOM</geometry>
+                    <type>Dirichlet</type>
+                    <component>1</component>
+                    <parameter>displacement_bottom</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <geometrical_set>TwoElementsWithInterface</geometrical_set>
+                    <geometry>PLY_RIGHT</geometry>
+                    <type>Dirichlet</type>
+                    <component>0</component>
+                    <parameter>zero_u</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+        </process_variable>
+        <process_variable>
+            <name>displacement_jump1</name>
+            <components>2</components>
+            <order>1</order>
+            <initial_condition>displacement0</initial_condition>
+            <boundary_conditions>
+            </boundary_conditions>
+        </process_variable>
+    </process_variables>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_newton</name>
+            <type>Newton</type>
+            <max_iter>100</max_iter>
+            <linear_solver>general_linear_solver</linear_solver>
+        </nonlinear_solver>
+    </nonlinear_solvers>
+    <linear_solvers>
+        <linear_solver>
+            <name>general_linear_solver</name>
+            <eigen>
+                <solver_type>SparseLU</solver_type>
+                <scaling>1</scaling>
+            </eigen>
+        </linear_solver>
+    </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..0d37bfd14e7bb42237eca6bdf690ae78c662222d
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:935a92068d70978a79faaffb8c578e817c55f82bad36df3384ad83c849ca913d
+size 1436
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_100_t_100.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_100_t_100.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..43ebc2efb6f167adb50d35e56a742cc91102cc28
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_100_t_100.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:a669d9b6d5702da69c87b70fee0db873d2fb8742c294fa7f1e0fda58d6afd42c
+size 4843
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_10_t_10.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_10_t_10.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..a05efb97f2df09d949803592f9cc2b128043e640
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_10_t_10.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:76bfb0967635791d150215eea1bfbd0f451278e9da77a06ac0b279bc7306f385
+size 4891
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_20_t_20.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_20_t_20.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..69e7c43f774abcf4341b8ca459a3abefcb43241e
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_20_t_20.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:1606c52ce0a2978f1b98e2997041df0206e4dfca6b63e3a46bfb50cb7c4a0db3
+size 4860
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_30_t_30.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_30_t_30.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..b5d883d1dd0aebf0f0fc675646714923c14a62fd
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_30_t_30.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:d4d90f41815a912240d2534a6e35e40d257c2b5f49da60686a1df12409697bda
+size 4887
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_40_t_40.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_40_t_40.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..ab2f2a4fa845fea3714ff58cee9367963f38686a
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_40_t_40.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f0b05467b511d32ea1fc90f09b140125f3844122424c0a5e778075fa2a9d2b06
+size 4947
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_50_t_50.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_50_t_50.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..c40abbcaf92e2be6b93ef44079c6237ba419254c
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_50_t_50.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:535798480751b38ace96c594ebed0ed8d9a1e7841195a8b3282da30203a405d5
+size 4862
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_60_t_60.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_60_t_60.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d6604ca2534c3f16a8d3c6cacce2ca5095992fc9
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_60_t_60.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:26c72f6fd82d4dbe20dece6ec63277f69e2951842364e6df05ac7e1a51d8a0c5
+size 4943
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_70_t_70.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_70_t_70.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..4a6812b23e0a823bc3ec66f5785556a239359339
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_70_t_70.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:5f2052e25a91f55ba3d44d90969b942ebbb6ba7362b70f502e9f53bfc3e408b5
+size 4917
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_80_t_80.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_80_t_80.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..ffb26b9c870704e8ed83b7c90b5021f7a699633e
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_80_t_80.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:0fe6d62008cb5fdd863cd1009bd5347bccd4c85cdca9a74646fd8cef9e28e9e1
+size 4845
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_90_t_90.000000.vtu b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_90_t_90.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..067d2b4cdd394cafb698eb9fb8abf258bc1b8dcd
--- /dev/null
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path_pcs_0_ts_90_t_90.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:1d84297a8335acbf31a7114bbfa586bbe3ce20ce331c725b19b744bbb21beb3b
+size 4921