Herezh++

Gérard,

Hugo a relevé une valeur étrange de Cos3hpi_sig dans un calcul de traction sur élément membrane en contrainte plane. Dans le cas contrainte plane, normalement, Cos3phi_sig devrait être correct puisque SIG33=0.

Pourtant, il trouve cos(3phi) = 1.53
C'est ce que je trouve également dans l'exemple pièce jointe sur un élément quadrangle.

J'ajouterais que le calcul des mises (deformation et contrainte) sont bons au sens 2D (donc en considérant composante 33 = 0).
Mais que les calculs de Lode ne sont pas bons (toujours au sens 2D), ni sur SIG, ni sur EPS.

Pour exemple dans le détail, je trouve au dernier incrément :
SIG11 = 2.603228571430e+01
SIG22 = -2.573789813879e-03
SIG12 = 0.000000000000e+00
contrainte_mises = 2.603357270462e+01
Cos3phi_sig = 1.533389226440e+00

Pour SIG33 = 0, on doit obtenir :
contrainte_mises = 26.033572 => OK
Cos3phi_sig = 0.9999999670 => NOT OK

Et pour les déformations :
EPS11 = 8.677685950413e-02
EPS22 = -2.604086501368e-02
EPS12 = 0.000000000000e+00
def_duale_mises = 6.820913748080e-02
Cos3phi_eps = 8.987084506902e-01

Pour EPS33 = 0, on doit obtenir :
def_duale_mises = 0.068209 => OK
Cos3phi_eps = 0.7858437489 => NOT OK

sachant que pour SIG et EPS, je calcule cos(3phi) de la même manière :
cos(3phi) = 27/2 * J3 / sqrt(-3 * J2)**3

avec :
d = déviateur
J2 = -1/2 trace( d . d )
J3 = det( d )