Herezh++

Gérard,

j'obtiens des valeurs erronées dans les torseurs de réaction et la grandeur EPAISSEUR_FINALE avec une loi HYPER_EXTERNE_W.
Plus précisément, l'erreur apparait quand on utilise HYPER_EXTERNE_W pour définir la partie sphérique du potentiel.

Je suis dans un contexte de LOI_CONTRAINTES_PLANES, je ne sais pas si ça a une importance.

Dans l'exemple joint (traction sur un quadrangle), j'ai défini une loi HYPER_EXTERNE_W dont la partie déviatorique est exactement égale à une loi mooney-rivlin si C10=C01.
Et sur la partie sphérique, en se plaçant en quasi-incompressible, on a un comportement du style :
pression hydro = K * ln(V).
Les constantes, les fonctions nD et la loi de comportement sont dans le répertoire MAT_HYPER/ (constantes.hz_courbe, fonction_nD.hz_courbe et loi.hz_loi)

Dans le fichier MAT_HYPER/loi.hz_loi, on peut choisir :
1) partie sphérique : soit ISOHYPERBULK, soit un HYPER_EXTERNE_W très similaire (en quasi-incompressible, on a dans les 2 cas une loi telle que : pression hydro = 2000 * ln(V)
2) partie déviatorique : soit MOONEY_RIVLIN_3D, soit son exacte copie HYPER_EXTERNE_W si C10=C01=5, pour un module de cisaillement de l'ordre de 2G = C10+C01 = 10

En testant un peu toutes les combinaisons sphérique/déviatorique, je m'aperçois d'un bug si la partie sphérique est définie avec HYPER_EXTERNE_W (aucun pb avec la partie déviatorique).
Herezh écrit alors dans le .maple des valeurs erronées pour le torseur selon la direction X de traction (composante Rx colonne 2) et pour l'épaisseur finale (EPAISSEUR_FINALE colonne 20, j'ai testé aussi EPAISSEUR_MOY_FINALE).

Le fichier gnuplot gnu trace 3 résultats :
1) courbe EPS11 - SIG11 => quelque soit la combinaison de choix sphérique/déviatorique, on a la courbe suivante (exemple ici dans le cas avec ISOHYPERBULK/3K=6000) :

2) concernant la force de réaction, je trace 3 résultats différents : soit directement la valeur Rx du .maple, soit la contrainte multipliée par la section initiale (1 x 0.038), soit la contrainte multipliée par la section actuelle (évaluée en fonction de la déformation 22 + isotropie)
Par ce test, j'en déduis que Herezh calcule très bien la contrainte, ainsi que toutes les grandeurs (force réac, épaisseur, etc...). Ce n'est qu'un pb de sortie maple.

Par exemple, dans le cas ISOHYPERBULK + déviatorique HYPER_EXTERNE_W, on a :

mais si on remplace ISOHYPERBULK par le sphérique HYPER_EXTERNE_W, on a :

La sortie Rx (courbe rouge) est largement au-dessus du résultat attendu. Et pourtant, en zoomant avec le fichier gnu, on verrait que les 2 autres résultats sont très similaires au cas ISOHYPERBULK => ce n'est qu'un problème de sortie

3) épaisseur :
tout à fait logique dans le cas ISOHYPERBULK (diminution d'épaisseur depuis la valeur ini de 0.038) :

erronée avec sphérique HYPER_EXTERNE_W :

Il reste à comprendre pourquoi :
- est-ce lié au fait que la partie sphérique de HYPER_EXTERNE_W dépend de III_B au lieu de V ?
- ou simple bug qui traine (un décalage de colonne, etc... ?)

merci,
Julien