diff --git a/MaterialLib/FractureModels/CohesiveZoneModeI.h b/MaterialLib/FractureModels/CohesiveZoneModeI.h
index 8bbfc83acae6ece64948883fc090f0919a2faa61..23931d7e16f37da6498a0c113825744adcd626c1 100644
--- a/MaterialLib/FractureModels/CohesiveZoneModeI.h
+++ b/MaterialLib/FractureModels/CohesiveZoneModeI.h
@@ -49,6 +49,11 @@ struct MaterialPropertiesParameters
     double fracture_opening_at_residual_traction(double const t,
                                                  X const& x) const
     {
+        if (peak_normal_traction(t, x)[0] == 0.0)
+        {
+            return 0.0;
+        }
+
         return 2 * fracture_toughness(t, x)[0] / peak_normal_traction(t, x)[0];
     }