Herezh++

oui, effectivement cela ne fonctionne pas avec une loi de Mooney Rivlin et toutes les lois qui utilisent les mêmes invariants.
Pour l'instant les seules lois hyper-élastiques qui peuvent fonctionner avec les contraintes planes sont les lois qui utilisent les invariants: V = volume finale / volume initiale, Qeps = intensité du déviateur de déformation, cos(3*phi), phi étant l'angle de phase (ou de Lode) dans le plan déviatoire. Donc les lois :
1) ISOHYPER3DFAVIER3
2) ISOHYPER3DORGEAS1
3) ISOHYPER3DORGEAS2
4) ISOHYPERBULK3
5) ISOHYPERBULK_GENE

Toutes ces lois utilisent la déformation d'Almansi et la convexité des potentiels n'est réellement assurée que jusqu'à 10-15 % environ, au dela on peut avoir plusieurs solutions et donc des divergences.

Les autres lois vont également être disponibles dès que j'aurai fait les modifs nécessaires pour le calcul de l'opérateur tangent d_sig / d_eps en curviligne, sachant que cela existe déjà en coordonnées absolues.

Donc pour des premiers tests le plus simple est d'utiliser les lois hyper citées plus haut. Elles permettent de rendre compte de comportement déjà très non linéaire.

Affaire à suivre,
@bientôt