diff --git a/MaterialLib/SolidModels/MFront/MohrCoulombAbboSloan.mfront b/MaterialLib/SolidModels/MFront/MohrCoulombAbboSloan.mfront
index f9f7bbf7e87e19859975e04a88541899e8efc1b0..97d5cfa1b82b55d45d85c454ec6544a557477f94 100644
--- a/MaterialLib/SolidModels/MFront/MohrCoulombAbboSloan.mfront
+++ b/MaterialLib/SolidModels/MFront/MohrCoulombAbboSloan.mfront
@@ -82,6 +82,11 @@ a.setEntryName("TensionCutOffParameter");
     cos_3_lodeT = cos(3. * lodeT);
     sin_3_lodeT = sin(3. * lodeT);
     tan_3_lodeT = tan(3. * lodeT);
+
+    // Compute initial elastic strain
+    const auto S = invert(D);
+    eel = S*sig;
+
     // elastic prediction
     const auto sig_el = computeElasticPrediction();
     const auto s_el = deviator(sig_el);