From 25e6f768bad351517a57191557d72c61516c5e25 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?G=C3=A9rard=20Rio?= Date: Wed, 15 May 2024 17:47:27 +0200 Subject: [PATCH] V 7.031 : corr boulette sur Hypo_hooke2D_C (cal epsBH33) --- Parametres/EnteteParaGlob.h | 2 +- comportement/Hypo_elastique/Hypo_hooke2D_C.cc | 8 ++++---- comportement/lois_combinees/LoiContraintesPlanes.cc | 2 +- 3 files changed, 6 insertions(+), 6 deletions(-) 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 :