Herez++_test

suite à notre conversation d'hier, qq remarques sur les algos de Newton utilisés en CP et CP2, en relation avec le pilotage du critère plis.
- Dans une première étape, systématiquement tout l'incrément de charge est appliqué.
- Puis s'il y a divergence, on sub-divise l'incrément.
Après réflexion, j'ai l'impression que c'est correct de garder cette solution comme solution par défaut. Car le pas de chargement est imposé par l'algo global et si tout ce passe bien, il n'y a pas de raison de subdiviser dès le début.
Ceci étant, je vais quand même introduire la possibilité de considérer un pas maxi différent de 100% de l'incrément par défaut. ça peut-être intéressant par exemple si on a une divergence et que l'on fait un restart en voulant examiner spécifiquement le comportement local sans changer l'incrément global...
à suivre ...

au niveau du calcul de la variation de volume je vais intégrer explicitement la prise en compte d'une déformation logarithmique cumulée qui sera explicitement attribuée aux lois hypo élastique ce qui permettra d'être affranchi du passage par le test relatif aux def Almansi
Cela a des répercussions à différents endroits donc il faut que je prenne des précautions mais j'ai l'impression que c'est la bonne solution et cela ira dans le sens de ce que je prévois de faire en cours d'année relativement à la prise en compte des mesures de type log.
bon... c'est démarré ... à suivre

ok. Et quels sont les changements de paramètres que tu as fait ?

pour aller plus loin, j'ai mis en pièce jointe un exemple de loi orthotrope E=200, NU=0.8 dans le plan de la membrane et NU=0.19 dans l'épaisseur. Dans le .info, il y a 3 lois :
- MAT_sans_pli : loi orthotrope hypoélastique sans critère pli => pas de souci avec cette loi
- MAT_avec_pli : idem avec pli => nombreux messages d'erreur de dépassement 0.5
(j'ai mis le .log de ce calcul dans l'archive pièce jointe)
- MAT_avec_pli_almansi : même paramètres et avec pli mais Almansi => tourne correctement avec un petit passage compliqué avec message erreur.

Graphiquement, ça donne ça (toujours dans un contexte où le critère pli doit donner la loi de Hooke sigma = 200 * espilon)

- Loi hypoélastique MAT_avec_pli :
ça part bien au début, puis le même genre de non linéarité qu'en isotrope, le tout calé à peu près sur la loi de Hooke au sens log

- Loi almansi MAT_avec_pli_almansi :
c'était pour comparer/vérifier. Comme en isotrope, le comportement est très calé sur la loi de Hooke au sens Almansi, mais cette fois, un petit passage compliqué que l'on voit sur la courbe (décroche un peu puis se recale sur le résultat théorique)