Herezh++

Ajout fichier matériaux

Frank,
j'ai essayé d'analyser le problème et j'ai quelques conclusions.
1) il y a bien eu une modification dans la loi critère qui peut expliquer que les résultats soient différents avec la version actuelle. Cependant, il s'agit bien d'une amélioration et non d'un bug (a priori). Du coup, les résultats divergents que l'on observe ici montrent sans doute que la méthodologie globale n'est pas assez robuste.
2) En examinant les messages d'erreur et en debug, j'ai pu constater qu'il y avait à partir d'un certain moment une énorme variation d'épaisseur. Je pense que c'est du à l'algorithme de Newton qui applique la contrainte de sig33 = 0. Ce n'est même pas au niveau des plis. Une explication est le fait que la cinématique appliquée au travers de l'algo de relaxation dynamique, introduit de fortes oscillations locales de déformations qui mettent certaine fois en défaut l'algo de Newton.
Je pense (mais à vérifier) que l'amortissement visqueux serait moins impactant à ce niveau.
Ceci étant, une solution pour l'amortissement cinétique, que l'on a déjà évoquée est d'avoir une loi de base sans algo de newton interne, qui soit utilisée jusqu'à une certaine stabilité de la structure et seulement ensuite l'introduction de la loi finale.

Une solution lourde et pas pratique, serait de stabiliser la géométrie avec une première loi. Sauvegarder la forme. Repartir de la forme avec la loi définitive.
-> à mon avis ce n'est pas la bonne voie.

Une autre solution serait d'utiliser une loi des mélanges, pondérée par la précision actuelle sur le résidu, qui permettrait d'entrer deux lois, la première simple, la seconde de plis. Au début du calcul, c'est la première loi qui est utilisée et lorsque la précision est suffisante, on bascule sur la seconde loi.
Un premier intérêt est qu'avec la loi des mélanges on peut spécifier que tant que la pondération sur la seconde loi = 0, la seconde loi n'est pas activée. Du coup le calcul n'utilise en fait qu'une seule loi ce qui permet de garantir une rapidité.

Un second intérêt est que cette technique peut se dérouler dans un même incrément de manière automatique. De plus, ceci peut-être vrai sur plusieurs incréments...

bon... mais actuellement la pondération disponible sur la loi des mélanges, n'intègre pas les variables globales ... donc il faut que j'introduise cela dans le code.

Je pense que c'est la bonne solution. J'espère pouvoir mettre cela en place rapidement (mais en dehors des vacances !).

à suivre ...

J'ai une 2e situation de non régression entre ces 2 version de Herezh qui me parait plus suspecte. Le cas joint (que j'espère complet !) comporte 2 incréments de calcul, il s'agit du problème d'équilibre du BSO. Il converge rapidement (< 14min) avec la version 7.771, c'est donc déjà un petit succès !
Avec la version > 7.771 le premier incrément de charge donne des résultats identiques (exactement même fichier log). Par contre il diverge à la première itération de l'incrément 2.
Est-ce lié au calcul de la masse en début d'incrément ?
Merci

Frank

Gérard,

J'ai un cas plus simple (3 fichiers et quelques dizaines de secondes de calcul) qui met en évidence ce problème de divergence qui apparait à partir de la version 7.772. La divergence intervient systématiquement au 2e incrément même si il ne se passe rien !

Ce problème assez handicapant car je ne peux pas bénéficier des toutes les évolutions que tu a apportées à Herezh depuis et celles à venir, et je ne vois pas comment le contourner...

Merci

Frank,
une initialisation mal placée à un changement d'incrément !!
Sur ton dernier test c'est ok. Je mets la nouvelle version 6.782, dans la partie test, car je n'ai pas eu le temps de passer tous les tests de validation.
J'espère que c'est ok,
@+