Herezh++

Bonjour Frank,
effectivement il y avait un certain nombre de points à regarder.
1) au niveau de la déformation logarithmique que tu utilises (ce n'est peut-être pas une bonne idée, car pour l'instant la matrice tangente n'est pas très précise, mais ce n'est pas une erreur), donc il y avait un pb d'initialisation de dérivée. Dans le cas implicite, ou purement explicite, c'était ok. Mais en relaxation dynamique, en fait on est en explicite, tout en voulant de temps en temps calculer la raideur globale (pour le calcul de masse) ce qui correspond à une partie qui relève habituellement d'un calcul type implicite ... donc ce cas particulier d'une prise en compte d'un ddl externe (ici la température) pour une succession implicite explicite implicite, etc. n'était pas parfaitement pris en compte -> c'est a priori ok maintenant
Je recompile une nouvelle version dans ce sens: 6.805
2) un second pb vient de ton choix d'évolution de la température et de la loi associée. Là ce n'est pas un pb d'Herezh++ mais un pb de compréhension du fct d'Herezh.
En relaxation dynamique, dans le cas que tu as utilisé, c-a-d utilisation de la raideur réelle pour l'estimation de la masse virtuelle, le calcul démarre par la construction de la matrice de raideur (calcul associé au type implicite).
Or dans ton cas, au début du calcul, la température est nulle ! et pour une température nulle et bien compte tenu de la loi thermodépendante que tu utilises, le module d'Young est nulle !! d'où une masse nulle, d'où une divergence !!

Donc ce qu'il faut faire:
On impose une température initiale non nulle: et ensuite pendant le calcul, si on demande un recalcul de la raideur, ce sera la température en cours (= la température initiale + la température que l'on impose en chargement) qui sera utilisée.

Par exemple dans mon cas j'ai modifié ta mise en donnée dans le .info avec une température initiale de 10 °

1. nom du maillage | Ref noeud | Blocages |
   #---------------------------------------------------------------------
   N_to 'TEMP = 10 ' # température initiale

initialisation -----------
#---------------------------------------------------------------------

La gestion des ddl (ou variables) qui ne sont pas des déplacements est différentes des ddl méca classique. Les ddl de position sont par défaut initialisées aux positions données par le maillage, alors que les autres ddl sont initialisées à 0 !!! à moins de spécifier une initialisation particulière.
NB: c'est expliqué dans la doc, mais ce n'est pas toujours évident de s'en souvenir.

bon... à tester avec la nouvelle version.

Merci Gérard pour cette réponse. Je comprends bien l'importance de l'initialisation pour avoir une raideur non nulle au démarrage. Pour la déformation logarithmique c'était juste pour ce cas test pour vérifier par rapport à un calcul analytique.

J'ai ajouté l'initialisation :

initialisation
N_to 'TEMP = 10 ' # température initiale

et j'ai la même erreur mais je ne suis pas sûr d'avoir placé cette instruction au bon endroit dans le jeu de données (voir le fichier joint).
As-tu fais le test chez toi ?

Quand tu écris :"dans le cas que tu as utilisé, c-a-d utilisation de la raideur réelle pour l'estimation de la masse virtuelle" est-ce que cela signifie qu'il y a d'autres façon de faire vauier la raideur en fonction d'un champ thermique ?

Frank