diff --git a/MaterialLib/SolidModels/MFront/Lubby2.mfront b/MaterialLib/SolidModels/MFront/Lubby2.mfront
index e3f7a5d1e74d1fcb22787aeb0467f09c912ec190..a3d0048369d1d8d265e3086d1a41ae4ef0f86058 100644
--- a/MaterialLib/SolidModels/MFront/Lubby2.mfront
+++ b/MaterialLib/SolidModels/MFront/Lubby2.mfront
@@ -45,6 +45,8 @@ epsM.setEntryName("MaxwellStrain");
 @InitLocalVariables
 {
     mu = computeMu(young, nu);
+    // Compute initial elastic strain
+    eel = 1. / (2. * mu) * sig - nu / young * trace(sig) * Stensor::Id();
 }
 
 @Integrator
diff --git a/MaterialLib/SolidModels/MFront/Lubby2mod.mfront b/MaterialLib/SolidModels/MFront/Lubby2mod.mfront
index 6a89e44f30e5947e9989f5c0b7d4d6dbbac239ce..42ee281d0e32029f5c6c57ec36219fed97ffba44 100644
--- a/MaterialLib/SolidModels/MFront/Lubby2mod.mfront
+++ b/MaterialLib/SolidModels/MFront/Lubby2mod.mfront
@@ -46,6 +46,8 @@ epsK.setEntryName("KelvinStrain");
 @InitLocalVariables
 {
     mu = computeMu(young, nu);
+    // Compute initial elastic strain
+    eel = 1. / (2. * mu) * sig - nu / young * trace(sig) * Stensor::Id();
 }
 
 @Integrator