Herezh++

Gérard,

voici un ticket qui peut attendre. Et je pense que ça peut être en lien avec ce ticket (erreur de sortie maple avec loi HYPER_EXTERNE_W sur la partie sphérique du potentiel) :
https://herezh.irdl.fr/issues/368

Si on applique LOI_CONTRAINTES_PLANES sur une loi 3D hyperélastique dont la partie sphérique est définie via HYPER_EXTERNE_W, la loi de comportement renvoie un message d'erreur du style :

*** pb dans le calcul de la nouvelle epaisseur:  coef de compressibilite moyen = 0.00302637 est inferieur a delta la trace de sigma = 0.00907927 ce qui va donner une epaisseur finale negative !!  on continue sans changer l'epaisseur  ... 
 QuadraMemb::CalEpaisseurMoyenne_et_vol_pti(...) 

(ps : je viens de constater que d'après les valeurs indiquées, on a un coef 3 entre 0.00302637 et 0.00907927. Il y a peut-être quelque chose à creuser de ce côté. Confusion entre K et 3K pour les lois HYPER_EXTERNE_W ? De mon côté, je suis sûr d'avoir un coef de compressibilité tangent Kt égal à environ 2000, soit K3/3 avec K3=6000 comme indiqué ci-dessous. Avec Kt défini au sens thermodynamique :
Kt = V d(-p)/dV = V d^2 W / d V^2 (W : potentiel config initiale)
)

Dans l'exemple joint, j'ai construit 3 lois dont la partie sphérique est de type ISOHYPERBULK3, K3=6000 (je ne détaille pas la partie déviatorique, ça n'a pas d'importance).
Voici les lois :
- MAT_hyperbulk : partie sphérique via loi ISOHYPERBULK3
- MAT_hyper_W : partie sphérique identique à ISOHYPERBULK3 en passant par une loi HYPER_EXTERNE_W
- MAT_hyper_W_GENE : partie sphérique identique à ISOHYPERBULK3 en passant par une loi ISOHYPERBULK_GENE

Le cas test est un quadrangle dont j'impose tous les déplacements sur 3 de ses noeuds et je laisse un noeud libre selon X et Y.
Chaque loi a son propre fichier .info. J'ai mis le .log du calcul et le .maple obtenu.

Le fichier gnu trace les résultats suivants au premier point d'intégration :

1) EPS11 - contrainte hydrostatique => on voit bien que ce sont les mêmes lois :

2) temps - V=v/v0 => ça donne l'impression d'une erreur de sortie maple (colonne 12 = VOLUME_PTI) :

Mais au final, erreur de sortie maple ou pas, il y a le message d'erreur sur l'épaisseur, qui n'apparait pas avec la loi ISOHYPERBULK3 et ISOHYPERBULK_GENE.
Et Frank a rencontré des problèmes de convergence avec des lois HYPER_EXTERNE_W sur la partie hydro (test d'une version ISOHYPERBULK_GENE en cours, à voir si c'est mieux).

Pour l'instant, je laisse de côté HYPER_EXTERNE_W pour la partie sphérique et vais me focaliser sur des lois ISOHYPERBULK_GENE, en espérant que ça passe mieux dans les calculs.
Donc pas d'urgence...

merci
Julien