Herezh++

Bonjour Gérard,

ci-joint un calcul avec un pb de sortie Gmsh. Le calcul converge jusqu'à la fin. Le dernier incrément convergé est le n° 73. Mais on dirait qu'il y a un bug dans la dernière sortie Gmsh. Les 2 derniers $Nodedata des fichiers .pos ont tous les deux un numéro d'incrément égal à 72.

Pas grave en soi, mais ça met en échec certains post-traitements hors Herezh (par exemple : impossible d'utiliser hz_visuGmsh.pl).

version Herezh : 7.056 sur mac

Le contexte est un calcul axisymétrique de compression d'un joint torique dans un logement avec prise en compte du contact. La compression est exercée par un quadrangle.
Le bug est assez difficile à reproduire. Je n'ai pas essayé de faire un exemple un peu plus léger mais le calcul n'est pas trop long du fait de lois HYPOELAS assez simple et éléments triangle linéaire (environ 2min sur ma machine).

(ps : le résultat est assez joli, avec utilisation de TYPE_PENALISATION_PENETRATION 12, mais ce n'est pas vraiment le sujet du ticket).

Je n'ai pas vraiment de piste pour cerner le problème :
- il y a la présence du critère d'arrêt : min(Res,Res/MaX(Reaction_et_PVExterne))_verification_ddouble
Peut-être que le modifieur "_verification_ddouble" a un rôle dans le bug.

- La présence de contact est peut-être en cause aussi.