diff --git a/MaterialLib/FractureModels/MohrCoulomb.cpp b/MaterialLib/FractureModels/MohrCoulomb.cpp
index 1aad3c972347ab9e1649cbace8cdcb0814d172b2..1341cef13ab3dadde21a8887b4aa402f70008994 100644
--- a/MaterialLib/FractureModels/MohrCoulomb.cpp
+++ b/MaterialLib/FractureModels/MohrCoulomb.cpp
@@ -99,7 +99,7 @@ void MohrCoulomb<DisplacementDim>::computeConstitutiveRelation(
         sigma.noalias() = Ke * (w - w_prev);
 
         sigma.coeffRef(index_ns) *=
-            logPenalty(aperture0, aperture, _penalty_aperture_cutoff);
+            logPenaltyDerivative(aperture0, aperture, _penalty_aperture_cutoff);
         sigma.noalias() += sigma_prev;
     }
 
@@ -179,8 +179,8 @@ void MohrCoulomb<DisplacementDim>::computeConstitutiveRelation(
 
             sigma.noalias() = Ke * (w - w_prev - state.w_p + state.w_p_prev);
 
-            sigma.coeffRef(index_ns) *=
-                logPenalty(aperture0, aperture, _penalty_aperture_cutoff);
+            sigma.coeffRef(index_ns) *= logPenaltyDerivative(
+                aperture0, aperture, _penalty_aperture_cutoff);
             sigma.noalias() += sigma_prev;
         };