Herezh++

J'essaye d'utiliser la loi HH dans un calcul BSO.

Après le déploiement, lorsque l'on applique les CLL avec la loi, HH Herezh se met en erreur "loi de comportement" dès la 1ere itération.
Pour rendre progressive l'application des conditions de symétrie cyclique je passe par une phase de contact. Le déroulement est :

inc1 : déploiement, loi 1 x ISOELAS2D_C + 0 x LC2+pli
inc2 : même chargement et changement loi -> 0.001 x ISOELAS2D_C + 0.99 x LC2+pli
inc3 : activation contact
...

Cela fonction bien pour les lois IE et HE. Avec ces 2 loi l'incrément inc2 demande 2 itérations seulement. Avec une loi HH (Julien fin 2019) il faut 500 itérations. C'est ok.

Lorsque l'on applique le contact avec la loi HH (inc3) le comportement de Herezh est étrange. Les forces de contact sont bien présentes (F_N max) mais les noeuds ne bougent pas ! L'énergie cinétique reste quasi nulle et max_deltaX= 0.

Le paramètre PENALISATION_PENETRATION permet de moduler les forces normales mais rien ne bouge. Que ce passe t-il ???

Ci-joint un jeu de données "minimal". Les incréments I et 2 durent moins de 2 min.

Merci de jeter un oeil sur ce comportement étrange...

Frnk