diff --git a/Parametres/EnteteParaGlob.h b/Parametres/EnteteParaGlob.h
index 7192f68..1a07356 100644
--- a/Parametres/EnteteParaGlob.h
+++ b/Parametres/EnteteParaGlob.h
@@ -41,7 +41,7 @@
 	   EnumLangue ParaGlob::langueHZ = FRANCAIS; // langue utilisée pour les entrées sorties
     int ParaGlob::nbComposantesTenseur = 1; // nombre de composantes par defaut a 1
     int ParaGlob::nivImpression = 2; // niveau d'impression
-    string ParaGlob::nbVersion = "7.030" ;   // numéro de version du logiciel
+    string ParaGlob::nbVersion = "7.031" ;   // numéro de version du logiciel
     string ParaGlob::NbVersionsurfichier = ""; // numéro de version lue en entrée fichier
     int ParaGlob::nb_diggit_double_calcul= 17; // nombre de chiffre significatifs utilisé pour 
                                         // l'affichage des double précision pour l'archivage
diff --git a/comportement/Hypo_elastique/Hypo_hooke2D_C.cc b/comportement/Hypo_elastique/Hypo_hooke2D_C.cc
index 14b0f0c..f2e727d 100644
--- a/comportement/Hypo_elastique/Hypo_hooke2D_C.cc
+++ b/comportement/Hypo_elastique/Hypo_hooke2D_C.cc
@@ -1546,8 +1546,8 @@ void Hypo_hooke2D_C::Calcul_SigmaHH (TenseurHH& sigHH_t,TenseurBB& DepsBB_,DdlEl
 //       save_resul.eps33 = Isigma * untier / Kc_use - log_var_surf;
        // modif 2 mai 2024: dans le cas d'une loi hypo élastique on utilise la forme incrémentale
        // delta eps_3^3 = -nu/(1.-nu) (delta eps_1^1+delta eps_2^2)
-         double DepsBH_33 =  -nu/(1.-nu) * DepsBH.Trace();//trace 2D
-       save_resul.epsBH33 += DepsBH_33;
+       double DepsBH_33 =  -nu/(1.-nu) * DepsBH.Trace();//trace 2D
+       save_resul.epsBH33 = DepsBH_33 * deltat + save_resul.epsBH33_t;
       };
 
     // traitement des énergies
@@ -1940,8 +1940,8 @@ void Hypo_hooke2D_C::Calcul_DsigmaHH_tdt (TenseurHH& sigHH_t,TenseurBB& DepsBB_,
        //       save_resul.eps33 = Isigma * untier / Kc_use - log_var_surf;
        // modif 2 mai 2024: dans le cas d'une loi hypo élastique on utilise la forme incrémentale
        // delta eps_3^3 = -nu/(1.-nu) (delta eps_1^1+delta eps_2^2)
-         double DepsBH_33 =  -nu/(1.-nu) * DepsBH.Trace();//trace 2D
-       save_resul.epsBH33 += DepsBH_33;
+       double DepsBH_33 =  -nu/(1.-nu) * DepsBH.Trace();//trace 2D
+       save_resul.epsBH33 = DepsBH_33 * deltat + save_resul.epsBH33_t;
       };
 
 
diff --git a/comportement/lois_combinees/LoiContraintesPlanes.cc b/comportement/lois_combinees/LoiContraintesPlanes.cc
index 734f702..b8feb60 100644
--- a/comportement/lois_combinees/LoiContraintesPlanes.cc
+++ b/comportement/lois_combinees/LoiContraintesPlanes.cc
@@ -1534,7 +1534,7 @@ void LoiContraintesPlanes::Grandeur_particuliere
                       break;
                  case DEFORMATION_LOGARITHMIQUE : case DEF_CUMUL_CORROTATIONNEL : case DEFORMATION_CUMU_LOGARITHMIQUE :
                      // cas d'une def logarithmique ou une approximation
-                      {  tyTQ(1+(*idecal)) = exp(hsurh0);(*idecal)++;
+                      {  tyTQ(1+(*idecal)) = log(hsurh0);(*idecal)++;
                       };
                       break;
                  default :