Herezh++

oui, cela parait un peu étrange, je vais essayer de regarder cela rapidement ...
En tout cas, je pense que tes tests sont pertinents.

L'analyse est bonne, mais en fait dans Herezh, pour les potentiels hyperélastiques, je tiens bien compte de la variation de volume due à la dilatation.
Il faut que l'info relatifs à la variation de volume soit disponible dans la mise en données.
Donc, a priori on devrait avoir quelque chose de correct ...
Bon ... je vais regarder ton fichier

dans le .info d'UMAT, j'ai essayé d'ajouter un mot-clé du style :
""""
dilatation_thermique
E_tout 1
""""
mais ça ne change rien

Je ne comprends pas pourquoi il faut indiquer une loi thermo dans le .info d'UMAT.

Le principe est que c'est Abaqus qui gère la dilatation thermique. ça pose un souci si il faut indiquer à Herezh une loi thermo dans le .info d'UMAT :

D'une part, cette loi thermo serait légèrement en désaccord avec Abaqus. Abaqus gère la dilatation en envoyant uniquement la partie méca de la déformation dans l'UMAT. Il retranche donc une partie thermique calculée en déformation logarithmique. Tandis que Herezh calcule une dilatation en déformation d'Almansi. En petite déf, on pourrait croire que c'est pas une grande différence. Mais dans le cas d'un gros module de compressibilité, ça va se voir.

D'autre part, il y a pas mal d'options dans Abaqus pour la dilatation thermique. ça va être compliqué et peu robuste de retranscrire ça en Herezh.

De plus, ça me mettrait en difficulté pour les travaux avec Benoit. Je leur ai créé un script qui prépare des fichiers d'UMAT automatiquement. Si en plus il faut gérer cette dilatation, bonjour l'usine.

Ma réponse concernait le fonctionnement actuel d'Herezh pour lequel la variation de volume n'est pas calculée via le tenseur de déformation, mais directement via les coordonnées.
L'utilisation directe du fonctionnement actuel d'Herezh nécessite bien d'indiquer une loi thermo-physique qui contient un paramètre de dilatation.
Je vais donc corriger le bug pour que cela fonctionne en l'état.
Ensuite si on veut un processus particulier pour la connexion avec Abaqus pour la loi hyperbulk, il faudra que je réfléchisse a implanter une méthodologie particulière qui n'est pas immédiate.
Dans un premier temps, pour suppléer à la loi hyperbulk, dont l'objectif est de retraduire le comportement volumique, l'utilisation de la partie sphérique d'une loi hypoélastique thermo-dépendante est peut-être la bonne solution. J'ai l'impression que le résultat serait équivalent ?

ok, je me suis enflammé.

la méthode la plus simple serait de calculer la variation de volume en faisant l'exponentielle de la trace du tenseur de déformation méca (STRAN + DSTRAN). A faire uniquement dans le cas Umat que Herezh repèrerait grâce au fait d'utiliser un élément de type POINT CONSTANT...ou alors créer un nouveau mot-clé du style :

materiau_Umat
E_tout 1

oui, on pourrait effectivement calculer la variation de volume de cette manière. Mais la difficulté n'est pas là en fait. Dans le cas des potentiels qui utilisent directement V, il me faut également la variation de V par rapport aux invariants de la déformation. Cette variation est différente suivant que l'on utilise la mesure d'Almansi ou la mesure logarithmique (cumulée ou non). J'ai fait des développements théoriques à ce niveau mais aucune implémentation pour l'instant. Il faut donc que je revois les implications d'un nouveau choix de calcul de V, en particulier pour le calcul de l'opérateur tangent, et ensuite que j'implémente cela d'où ma suggestion.