4Méthodologie pour la simulation. Code_Aster Simulation numérique du soudage
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
12/24
Clé :
U2.03.05
Révision :
4036
4 Méthodologie pour la simulation
Les retours d’expérience cités dans cette partie sont essentiellement issus des trois études suivantes :
•
Une tuyauterie en 316L chanfreinée soudée par le procédé TIG [2], avec du métal d’apport en treize passes (fiche CEA/EDF/FRAMATOME 3488) ; les 2 premières passes de cet exemple constituent un cas-test Code_Aster [V7.42.100]) ;
•
Une plaque chanfreinée en 316L soudée par le procédé TIG [4], avec du métal d’apport en deux passes de soudage (fiche CEA/EDF 2425) ; les différents fichiers de commande rattachés à cette étude sont stockés dans la base d’études de Code_Aster (SERVICE/BDD
Etudes) ;
•
Une plaque en 316L non chanfreinée, avec création d’un ligne de fusion par le procédé TIG
[3], mais sans métal d’apport (thèse INSA de Lyon de L. Depradeux soutenue en 2004).
Par souci de simplicité, nous les nommerons par la suite fiche 3488, fiche 2425 et thèse L.D 2004.
4.1
Type de modélisation : 2D ou 3D ?
En toute rigueur, le procédé de soudage est strictement tridimensionnel, l’apport de chaleur et
éventuellement de matière étant mobile et de vitesse souvent constante. Les champs de températures et mécaniques engendrés sont donc tridimensionnels et transitoires.
Il est donc préférable de réaliser une analyse transitoire tridimensionnelle.
Cependant, ce choix n’étant pas toujours possible du fait de la complexité, du temps calcul et de la capacité mémoire requis par des calculs 3D, on a souvent recourt à une modélisation 2D en simulation numérique du soudage (notamment pour du soudage multi-passe avec un grand nombre de passes).
Cas d’une modélisation 2D
•
Dans la majorité des cas, les simulations 2D considèrent une section transverse, perpendiculaire à l’avancée de la torche de soudage. Ces simulations sont réalisées avec l’hypothèse des déformations planes ;
•
Des calculs 2D axisymétriques sont également adoptés pour le cas du soudage de conduites cylindriques.
Le point faible des deux modélisations précédentes, c’est qu’on néglige l’effet de vitesse du soudage et qu’on suppose que le cordon (ou la chaleur) est déposé(e) simultanément sur toute la longueur de la pièce à souder.
•
D’autres choix, moins classiques, sont également envisageables : on peut choisir par exemple de mailler une section longitudinale (dans le sens du soudage) ou bien, si on considère que les températures sont constantes dans l’épaisseur, on peut choisir de mailler le plan moyen de la plaque. Ces simulations sont alors réalisées avec l’hypothèse des contraintes planes.
Retour d’expérience
•
Fiche 3488 : Le calcul du tube chanfreiné est réalisé en axisymétrique. Des comparaisons entre cette approche 2D axisymétrique et une approche 3D ont été réalisées dans [1] et montrent la bonne représentativité de l’approche 2D, même si celle-ci ne tient pas compte des effets de vitesse.
Thèse L.D 2004 : Sur la plaque non chanfreinée, plusieurs modélisations sont testées :
Calcul 3D thermo-mécanique complet ;
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
13/24
Clé :
U2.03.05
Révision :
4036
Calcul 2D en déformation plane où on modélise une section perpendiculaire à l’avancée de la torche ;
•
Calcul 2D en contrainte plane où on modélise le plan moyen de la plaque.
On constate, dans ce cas de figure, que le calcul 2D en déformation plane ne reproduit pas bien les résultats obtenus sur le calcul 3D, en terme de déplacements et de contraintes résiduelles, excepté pour ces dernières dans la zone centrale de la plaque. Quant au calcul 2D en contrainte plane, il donne des résultats très proches en contraintes résiduelles du calcul 3D.
4.2
Le maillage
Conseils généraux
•
Comme toute étude numérique, il est toujours difficile de donner précisément des valeurs sur la densité du maillage ; cela dépend beaucoup du problème physique considéré. Dans le cas du soudage, le maillage doit être suffisamment fin autour de la source de chaleur (cordons de soudure et zone fondue) pour appréhender correctement les forts gradients thermiques et mécaniques dans cette zone.
•
Dans le cas d’une modélisation thermique quasi-stationnaire en repère mobile (cf. § 3.3.2), la rapidité des calculs autorise une densité de maillage plus importante, notamment autours de la source.
•
Dans le cas de soudage sur aciers avec transformations de phases métallurgiques, il est nécessaire de mailler de façon suffisamment fine la ZAT. Ceci est d’autant plus important que le maxima de contraintes est généralement atteint en périphérie de ZAT, en zone partiellement austénitisée. La densité de maillage dans cette zone est donc primordiale. Les simulations de soudage sur un acier qui présente des transformations de phases nécessitent donc à priori des maillages plus volumineux en nombres d’éléments par rapport à des simulations de soudage sur aciers sans transformations (à étude identique).
•
Pour des calculs mécaniques 3D, il est très important de déraffiner le plus possible le maillage mécanique lorsqu’on s’éloigne de la zone soudée. En effet, la forte non-linéarité du problème
(et les temps de calculs qui vont avec) et l’aspect transitoire (nécessitant de nombreux pas de temps) limitent en pratique le nombre de nœuds du modèle.
•
Comme pour n’importe quel problème, il est nécessaire de prévoir dès l’étape de maillage les groupes destinés à recevoir les conditions aux limites thermiques et mécaniques (trajet de la source, surface d’échange…), ainsi que les zones destinées au post-traitement. Ceci est d’autant plus vrai pour le maillage thermique où le chargement (source chaleur) et les conditions aux limites (convection et rayonnement) se déplacent au fur et à mesure de la dépose des cordons.
•
Pour la zone fondue, il est souvent difficile de la représenter car, en l’absence de macrographie, sa forme est inconnue.
Retour d’expérience
•
Fiche 3488 : Les cordons de soudure peuvent être représentés géométriquement et maillés de façon plus ou moins complexes. Il y a 3 choix possibles, en allant du plus compliqué vers le plus simple :
On peut choisir de respecter à la fois le volume et la forme de la passe. La forme des cordons étant courbe, on devra mailler des surfaces à bords courbes, donc utiliser des
éléments finis au moins de degré 2 ;
On respecte seulement le volume de la passe ;
On respecte au mieux le volume de la passe.
Des comparaisons ont été faites dans [2] avec ces trois possibilités. Il s’avère que les résultats des calculs mécaniques diffèrent très peu d’un maillage à l’autre (à thermique équivalente). En conclusion, un maillage courbe n’apporte rien de significatif au niveau des résultats par rapport à un maillage polygonal. Il importe en revanche de respecter approximativement le volume des cordons déposés.
Manuel d'utilisation Fascicule u2.03 : Thermo-mécanique
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
14/24
Clé :
U2.03.05
Révision :
4036
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
15/24
Clé :
U2.03.05
Révision :
4036
4.3
Calcul thermique
Le point critique pour la partie thermique concerne la modélisation de l’apport de chaleur.
4.3.1 Introduction
Dans Code_Aster , le calcul de l’évolution thermique s’effectue par la résolution non linéaire de l’équation de la chaleur (opérateur THER_NON_LINE [R5.02.02]) dans le volume de la pièce, étant donné une condition initiale et des conditions aux limites thermiques sur les frontières. La résolution est transitoire, la source de chaleur se déplace sur le maillage.
La diffusion de la chaleur est traitée par une formulation enthalpique, qui est l’intégrale de la chaleur
T
spécifique sur la température
T =
∫
0
C p
u du
. On peut donc fournir, soit la conductivité et la chaleur spécifique
C p
en fonction de la température, soit, et c’est préférable (cf. remarque cidessous), la conductivité et l’enthalpie en fonction de la température.
Remarque sur le choix d’une formulation enthalpique
•
A la température de fusion, la capacité thermique (qui est l’énergie à fournir pour élever la température du corps) subit une discontinuité qui se traduit par la chaleur latente de fusion, qui représente l’énergie à fournir pour franchir la température de changement de phase. Dans le cas d’un mélange de constituant (c’est le cas des alliages), la fusion s’étale entre les températures de liquidus et de solidus
Tl
et
Ts
. Une formulation « enthalpique » de l’équation de la chaleur est utile pour la prise en compte des phénomènes de chaleur latente de changement de phases. En effet cette formulation évite de « rater » la transformation, ce qui risque d’arriver pour la formulation en capacité si les pas de temps ne sont pas suffisamment petits.
4.3.2 Gestion de l’ajout de matière (les cordons de soudure)
Description des différentes méthodes
Il existe plusieurs méthodes pour prendre en compte, dans le modèle thermique transitoire, l’ajout successif de matière lors des différentes passes de soudage :
•
La première méthode consiste à ne considérer qu’un seul modèle (au sens de Code_Aster défini par la commande AFFE_MODELE) contenant toutes les passes et à « désactiver » artificiellement les cordons non encore déposés en leur imposant une conductivité thermique nulle (
10 −
5
W / m° C
en pratique).
•
La seconde méthode consiste à emboîter les modèles thermiques (au sens de Code_Aster) : un maillage global de tous les cordons est utilisé, mais le modèle thermique n’est affecté que sur la partie du maillage correspondant au métal déjà déposé (lors de l’étape AFFE_MODELE). Les cordons non encore déposés ne sont pas inclus dans le modèle. Dans ce cas, il est nécessaire d’avoir autant de modèles thermiques qu’il y a de cordons déposés. Les différents modèles sont emboîtés les uns dans les autres, c’est-à-dire que le modèle thermique I contient les passes de
à
1
I
et le modèle thermique
I 1
contient les passes
1
à
I
et la passe
I 1
. Cette méthode est plus propre que la précédente dans la mesure où les éléments non déposés ne perturbent pas le calcul de la passe en cours, puisqu’ils ne sont pas inclus dans le modèle. En revanche, ceci pose un problème lors de l’enchaînement des calculs thermiques, les champs de température du modèle
I
n’étant pas définis en tous les nœuds du modèle correspondant à la passe suivante
I 1
. Il faut donc effectuer un prolongement des champs calculés d’un modèle à l’autre.
Pour ce faire, on procède en 4 étapes :
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
16/24
Clé :
U2.03.05
Révision :
4036
On commence par créer un champ de température à l’ambiante ( par la commande CREA_CHAMP (opération ‘AFFE’),
T20
) sur tout le maillage
On étend le premier champ de température calculé à la passe I (numéro d’ordre 1) en complétant par
T20
sur les nouvelles mailles (CREA_CHAMP opérations ‘EXTR’ puis
‘ASSE’),
On stocke ce champ dans une nouvelle structure de données de type evol_ther par la commande CREA_RESU,
⇒
on fait une boucle sur les numéros d’ordre restants et on répète les opérations 2 et 3 pour chaque numéro d’ordre en enrichissant la structure de données créée en 3 (mot-clé reuse de
CREA_RESU).
•
Enfin, la troisième méthode consiste, sur un seul modèle thermique (au sens de Code_Aster), à dédoubler les nœuds à l’interface de chaque cordon, ce qu’il faut prévoir dans l’étape de construction du maillage. Lors de l’affectation du chargement thermique (AFFE_CHAR_THER), on n’affecte que la partie du modèle qui correspond au métal déjà présent. Ainsi, les cordons non déposés ne voient pas le chargement, et le dédoublement des nœuds empêche le transfert de chaleur à l’interface des cordons non déposé (une température égale à l’ambiante est appliquée dans les cordons non encore déposé). Lors de l’ajout d’un nouveau cordon, on recolle les nœuds correspondant en appliquant l’égalité des températures aux nœuds dédoublés par la commande
LIAISON_GROUP. Cette méthode est certainement la plus fiable mais la construction du maillage devient fastidieuse si le nombre de passes est grand.
Retour d’expérience
•
L’avantage de la première méthode (un seul modèle et conductivité presque nulle pour le cordon non encore déposé) est sa simplicité de mise en œuvre. En revanche, elle peut conduire à des oscillations numériques de la température dues à la discontinuité de la conductivité aux interfaces entre les cordons, comme constaté dans la fiche 3488 [2].
Néanmoins, cette méthode est souvent choisie du fait de sa simplicité.
•
Les trois méthodes ont été testées sur la plaque 2425 [4] et comparées à des mesures expérimentales de température. Les trois calculs conduisent à des résultats sensiblement identiques aux points de mesures.
4.3.3 Modélisation de la source de chaleur
La modélisation véritablement complète du problème thermique nécessiterait la prise en compte des transferts de chaleur électro-thermo-fluides dans l’arc électrique, avec prise en compte des phénomènes électromagnétiques, et la modélisation des mouvements convectifs dans le bain fondu, en interaction avec le plasma de couverture, et les transferts de chaleur dans la partie solide. Bien qu’une littérature abondante existe en ce qui concerne la prise en compte de la modélisation de l’arc et du bain de fusion dans la simulation, nous nous placerons par la suite dans une optique de modélisation avec un code de calcul industriel, qui ne permet pas de modéliser les phénomènes présents dans l’arc et le bain fondu. Ceux-ci sont remplacés par la définition d’une source de chaleur de forme adéquate. Seule la conduction de la chaleur est donc modélisée, c’est pourquoi il ne sera pas fait allusion à la modélisation de l’arc et du bain par la suite.
Description des deux méthodes possibles
•
La première méthode consiste à imposer des cycles de température à la matière que l’on dépose. Ces températures imposées peuvent être appliquées soit sur le seul cordon déposé, soit sur l’ensemble cordon déposé plus zone fondue.
Par raison de simplicité, la fonction température imposée
T imp
est choisie, le plus souvent, constante en espace mais variable dans le temps, soit
T imp
=
T imp
x , y , z ;t=T
imp
t
;
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
17/24
Clé :
U2.03.05
Révision :
4036
De plus, on impose souvent cette température uniquement sur le cordon déposé : en effet, les formes et tailles de zones fondues sont rarement connues, et même si elles le sont, cela rend fastidieux la construction du maillage, surtout si le nombre de cordons est élevé ;
La méthode en température est conseillée dans le cas de modélisations 2D ou bien lorsque aucune donnée n’est disponible sur les paramètres de soudage. La fonction
T imp
t
a souvent la forme suivante :
T imp max
tT
{
t=
T
T t max
1
éventuellement,
ini
−
T t max
−
T min
t−t
3
−
t
2
2
T max
si t≤t
1 si t
1
≤
t≤t
2 si t
1
≤
t≤t
3
On prescrit une montée de la température initiale température éventuellement supérieure
T max
d’une durée
t
1
(la montée en température prescrite peut être linéaire ou non linéaire, exponentielle par exemple), puis un maintien à cette température pendant une durée (
t
2
−
t
1
T ini
à la température de fusion ou une
), et enfin éventuellement, on peut prescrire le début du refroidissement jusqu’à la température
T min
.
Les paramètres
t
1,
t
2,
t
3,
T min
peuvent être identifiés si des mesures de température existent ;
La suite du refroidissement s’effectue avec des échanges par convection et rayonnement.
•
La seconde méthode consiste à imposer un flux de chaleur surfacique (
volumique (
J / sm
3
) au cordon de soudure modélisé.
J / sm
2
) ou
Elle est préférentiellement utilisée lorsque l’énergie de soudage délivrée est connue : c’est le cas pour le soudage à l’arc, pour lequel on peut estimer la puissance délivrée UI (avec
U
la tension et
I
l’intensité du courant). De cette puissance délivrée, seul une fraction
UI
pénètre effectivement dans la pièce et participe à l’échauffement et la fusion des matériaux. Pour le procédé TIG, par exemple, le paramètre
est de l’ordre de 0,6 à 0,9 ;
Il est alors possible de répartir cette puissance, en surface ou en volume (ou les deux), sur le maillage. La répartition de ce flux de chaleur peut-être extrêmement simple (constante en espace dans le métal déposé), ou plus élaborée (modèle double-ellipsoïde à répartition gaussienne élaboré par J. Goldak). Le calage des paramètres de la source de chaleur retenue nécessite des mesures de température.
Par la suite, nous détaillons la démarche à suivre dans le cas d’une modélisation par flux de chaleur, en distinguant les simulations 3D et 2D. On suppose qu’il existe des données disponibles en température.
Cas d’une modélisation 3D avec flux de chaleur
•
Nous donnons ci-dessous quelques formes possibles pour la représentation spatiale des flux volumiques ou surfaciques :
Exemple d’un flux volumique, constant suivant soudage, ici
x
.
Q
V
=
Q
V
x , y , z ;t=Q
V
x ;t =
{
Q max
−
Q
L/ 2
max
L/2
x y
et
z
et variable dans la direction de si 0≤x≤
x−L si
L
2
L
2
≤
x ≤L
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
18/24
Clé :
U2.03.05
Révision :
4036
Exemple d’un flux surfacique et d’intensité maximale
Q max
.
Q
S
=
Q
S
x , y ; t
cylindrique à base circulaire de rayon
R
•
Lorsque des mesures de températures sont disponibles, les paramètres de cet apport de chaleur (tailles, valeurs de puissance maximale injectée) peuvent être calés de manière à reproduire au mieux, par calcul, les cycles thermiques mesurés.
•
Cette étape de calage peut être effectuée facilement et avec de faibles temps de calculs par l’intermédiaire d’un calcul quasi-stationnaire en repère mobile (THER_NON_LINE_MO). Cela suppose que le trajet de la source est rectiligne uniforme, à vitesse constante. Dans ce cas de figure, on suppose qu’un régime permanent s’établit, et l’équation de la chaleur écrite dans ce repère devient indépendante du temps. On obtient donc le champ de température dans ce repère (fonction de l’espace uniquement). On peut éventuellement récupérer la dimension temporelle par un changement de variable (
X =x – V.t
). On peut donc rapidement (un seul pas de temps) identifier les paramètres de la source, qui sont ensuite réutilisés dans le cas du calcul transitoire.
•
Une fois l’identification réalisée et si la puissance du procédé de soudage est connue, on pourra en déduire la valeur du paramètre
dans la puissance nette (afin de vérifier la cohérence de la valeur obtenue par rapport au procédé). En effet, l’intégrale de la densité volumique (respectivement surfacique) de flux de chaleur sur le volume choisi
(respectivement la surface choisie) doit être égale à la puissance nette, soit
∫
V
Q
V
x , y , z ;t dv= UI et
S
∫
Q
S
x , y , z ;t ds= UI
Dans les deux exemples présentés ci-dessus, cela donne :
Exemple du flux volumique :
S
∫
L
0
Q
V
x ;t dx=S∗
L
2
∗
Q max
=
UI
où
S
est la surface du cordon déposée perpendiculaire à la direction du soudage.
Exemple du flux surfacique :
∫
S
Q
S
x , y , z ;t ds=∗R
2
∗
Q max
=
UI
.
•
Le flux de chaleur étant identifié spatialement, il faut maintenant mettre en œuvre le déplacement de cette source sur le maillage pour le calcul thermique transitoire. On définit ainsi une fonction dans la direction de soudage, ici
x
, et du temps par l’opérateur
DEFI_NAPPE : un flux fonction du temps est appliqué successivement sur les mailles situées sous le trajet de la torche de soudage, en considérant une montée à valeur maximale, suivit d’une descente à zéro. La descente du flux sur une maille correspond à la montée du flux sur la maille adjacente suivante. Ainsi les sources sont « allumées » successivement le long de l’axe de soudage lors du passage de la torche, ce qui correspond à une source de chaleur d’intensité constante qui voyage dans l’espace le long de l’axe de soudage.
Cas d’une modélisation 2D avec flux de chaleur
•
Dans le cas d’une modélisation en 2D, l’approche est différente car il faut choisir un flux
(volumique ou surfacique) qui soit variable dans le temps et ceci afin de prendre en compte, dans le calcul 2D, l’approche et l’éloignement de la torche. Pour la répartition spatiale de ce flux, on choisit souvent un flux constant en espace. Un exemple de flux surfacique est donné ci-dessous :
Q r
=
Q r
x , y , z ;t =Q
r
x ;t=
Q
{
Q t
−
Q max max
1
max t
t−t
2
t
3
−
t
2
si t≤t
1 si t
1
≤
t≤t
2 si t
2
≤
t≤t
3
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
19/24
Clé :
U2.03.05
Révision :
4036
•
Là encore, il faudra identifier, à partir de mesures de température, les paramètres du flux choisi (dans notre exemple, les instants
t
1
,
t
2
et
t
3
ainsi que la valeur
Q max
).
•
Si l’on connaît la puissance délivrée par le procédé, on peut également en déduire le
:
t
∫
S
{
∫
Q r
t ds
}
=
S
∫
t
Q r
t dt=
UI
V
⇔
t
∫
Q r
t dt=
UI
SV
ce qui fait :
t
∫
Q r
t dt=Q
max
∗
1
2
t
3
t
2
−
t
1
=
UI
SV
où
V
est la vitesse de soudage connue et
S
la surface du cordon déposé connue. Pour affiner l’identification des paramètres de la source, on peut utiliser le paramètre un calcul 3D quasi-stationnaire en repère mobile.
Retour d’expérience
trouvé par
•
Thèse INSA L.D 2004 :
Pour la modélisation 3D transitoire de la plaque, plusieurs formes de flux surfaciques ont
été testées : triangulaire dans le sens du soudage, cylindrique à base circulaire et gaussienne. Pour ces trois flux, la répartition en température s’est avérée à peu près identique.
Pour la modélisation 2D (déformation plane), trois représentations ont été testées : température imposée, flux surfacique et flux volumique avec identification des paramètres sur des mesures. Seule la modélisation avec température imposée a conduit la fusion du métal déposé. Avec un flux volumique, la zone fondue n’a pas été reproduite.
•
Fiche 3488 : L’approche en température imposée a été utilisée sur le calcul axisymétrique du tube. Les résultats de cette simulation montrent que les calculs surestiment la zone fondue et les maxima des cycles thermiques. L’approche en température imposée est donc trop
énergétique et c’est l’approche en flux qui est conseillée pour cette étude.
4.3.4 Conditions aux limites thermiques avec l’environnement
Les conditions aux limites thermiques sont modélisées par des échanges convectifs et radiatifs de la pièce soudée avec l’environnement, qui s’écrivent :
−
∂
T
∂
n
=
hT −T
air
T
4
−
T
4
air
où
n
est la normale extérieure,
h
le coefficient d’échange convectif,
la constante de Stefan-
Boltzman et
l’émissivité du matériau.
Remarques
•
Dans le cas du soudage de pièces relativement épaisses, les valeurs prises pour le coefficient d’échange et l’émissivité n’ont souvent qu’une influence secondaire sur les calculs thermiques : en effet, le flux de chaleur échangé avec l’air ambiant pendant le chauffage et au refroidissement est très faible, comparativement au flux de chaleur par conduction dans la pièce provenant de l’apport de chaleur.
•
A proximité immédiate de la source, le flux émis par rayonnement est non négligeable dans la réalité, mais les inconnues correspondantes peuvent être intégrées dans la définition de la taille et la répartition de la source de chaleur.
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
20/24
Clé :
U2.03.05
Révision :
4036
•
Le problème thermique se ramène donc le plus souvent à un problème purement conductif et la non prise en compte des échanges radiatifs et convectifs avec l’environnement n’est pas essentielle.
4.3.5 Caractéristiques thermophysiques
Les propriétés thermophysiques intervenant dans l’équation de la chaleur sont fonction de la température. Il s’agit de la conductivité thermique
T
et de la chaleur spécifique
C p
ou de l’enthalpie.
Remarques
•
En ce qui concerne la conductivité thermique, il est parfois d’usage de l’augmenter artificiellement pour les températures supérieures à la température de fusion, de façon à
« prendre en compte » les phénomènes convectifs à l’intérieur du bain, et à homogénéiser la température de celui-ci. Dans la littérature, la conductivité est ainsi augmentée d’un facteur 2,
10 ou 100 selon les auteurs (c’est ce qui est fait dans le tube 3488 avec un facteur 100). En pratique, il est souvent raisonnable de prendre une conductivité constante à partir de la température de fusion, et on estime qu’il n’est pas nécessaire de l’augmenter artificiellement,
étant donnée que les inconnues concernant le bain de fusion (phénomènes convectifs, etc…) sont englobées dans la définition d’une source de chaleur ad hoc.
•
Lors du soudage, la chaleur latente de fusion (cf. remarque du paragraphe 3.3.1) est faible par rapport à l’apport de chaleur provenant du procédé, et il n’est pas forcément utile de la prendre en compte. En outre, cette chaleur latente de fusion peut également être incluse implicitement dans la définition de l’apport de chaleur.
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
4.4
Calcul métallurgique
Version default
Date :
10/08/2010
Page :
21/24
Clé :
U2.03.05
Révision :
4036
Le point critique pour la partie métallurgique concerne, essentiellement, la connaissance de données expérimentales afin d’identifier correctement les modèles de comportement métallurgiques.
4.4.1 Introduction
Dans cette partie, on ne s’intéresse qu’au matériau type acier, en particulier le 16MND5, qui peut comporter cinq phases métallurgiques différentes : la ferrite, la perlite, la bainite, la martensite, dites phase froide ou
, et l’austénite, dite phase chaude ou
.
Lorsqu’on chauffe un matériau, les phases se transforment en phase . Quand on refroidit le matériau, l’austénite se transforme, selon la vitesse de refroidissement, en ferrite et/ou perlite et/ou bainite et/ou martensite. Il faut donc définir pour le chauffage la cinétique de transformation
et pour le refroidissement la nature et la cinétique des transformations possibles .
Les phases métallurgiques en présence dépendent donc de la température et de la vitesse de température. Le calcul de ces phases est réalisé en post-traitement d’un calcul mécanique, par la commande CALC_META. Il faut préciser dans cette commande l’état métallurgique initial.
4.4.2 Modèles de comportement au chauffage et au refroidissement
Actuellement dans Code_Aster, pour des matériaux type acier, il n’existe qu’un seul modèle de comportement permettant de calculer, à chaque pas de temps, les proportions des différentes phases.
Néanmoins, ce modèle est différent au chauffage et au refroidissement. On trouvera dans
[R4.04.01] l’expression détaillée des cinétiques au chauffage et au refroidissement.
4.4.3 Interaction thermique => Métallurgique
Comme nous l’avons déjà précisé auparavant, il n’y a pas de couplage entre les calculs thermique et métallurgique. Or, les propriétés thermo-physiques (conductivité, chaleur spécifique ou enthalpie) du point matériel dépendent des proportions des différentes phases en présence : si ceci n’est en général pas trop préjudiciable vis à vis des prédictions mécaniques, cela peut l’être pour une prévision fine de la métallurgie finale (mais ce n’est pas, en général, l’objectif d’une simulation numérique du soudage).
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
22/24
Clé :
U2.03.05
Révision :
4036
4.5
Calcul mécanique
Les points critiques pour la partie mécanique concernent le choix de la loi de comportement et l’identification des paramètres de cette loi.
4.5.1 Principe du calcul mécanique
Le calcul mécanique s’effectue par résolution des équations d’équilibre (opérateur STAT_NON_LINE), en tenant compte de la dépendance des propriétés mécaniques vis à vis de la température et,
éventuellement des proportions des phases métallurgiques.
4.5.2 Gestion de l’ajout de matière (les cordons de soudure)
En mécanique, il existe deux méthodes pour gérer l’ajout de matière, équivalentes dans le principe à celles de la thermique.
•
On peut considérer un seul modèle mécanique comportant tous les cordons, où les nœuds à l’interface entre cordons sont dédoublés. Lors de l’ajout d’un nouveau cordon, on affecte des conditions à l’interface portant sur l’incrément des déplacements (LIAISON_DDL). Tout comme en thermique, cette méthode est certainement la plus fiable mais la construction du maillage devient fastidieuse si le nombre de passes est élevé. On n’a pas de retour d’expérience récent sur cette méthode.
•
La méthode la plus couramment utilisée en mécanique est celle où on a un seul modèle mécanique pour toutes les passes, les éléments non encore présents étant désactivés artificiellement par affectation d’un module d’Young « quasi-nul » (
E=10
−
11
. E 20 ° C
en pratique)). La valeur de module d’Young à affecter dans le métal non déposé résulte d’un compromis : si la valeur est trop élevée, des contraintes seront générées dans le métal non encore déposé, mais si elle est trop faible, cela peut entraîner des problèmes de convergence.
4.5.3 Loi de comportement
Dans Code_Aster, il existe différentes lois de comportement pour décrire le comportement non linéaire d’un matériau (écrouissages isotrope, cinématique, mixte, effet visqueux,….). L’utilisateur trouvera dans le document [U4.51.11] l’inventaire de ces modèles. Étant donnés que ces modèles ne sont pas spécifiques aux seules simulations du soudage, nous ne les décrirons pas.
Remarque
Néanmoins, nous pouvons citer les modèles mécaniques avec effet des transformations structurales, qui ont été développés pour les aciers, tel le 16MND5 et spécialement dans le cadre des activités soudage.
Ces modèles, qui sont décrits en détail dans les documents [R4.04.02] et [R4.04.03], permettent de modéliser les phénomènes suivants : comportement plastique ou comportement visqueux,
écrouissage isotrope linéaire ou non linéaire ou écrouissage cinématique linéaire, plasticité de transformation, restauration d’écrouissage d’origine métallurgique, restauration d’écrouissage d’origine visqueuse. On peut réaliser un calcul en petites déformations mais également en grandes déformations, soit avec l’option PETIT_REAC, soit avec l’option SIMO_MIEHE (les modèles avec écrouissage cinématique n’existent pas avec SIMO_MIEHE).
Les deux questions importantes au niveau du choix de la loi de comportement portent sur :
•
Doit-on prendre en compte les effets visqueux du fait des hautes températures engendrées lors du soudage ?
•
Quel type d’écrouissage doit-on considérer (isotrope, cinématique ou mixte) ? Cette question est importante en soudage car ce procédé implique des cycles de traction et compression ;
Manuel d'utilisation Fascicule u2.03 : Thermo-mécanique
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
23/24
Clé :
U2.03.05
Révision :
4036
Les effets visqueux sont souvent négligés en raison essentiellement du manque de données expérimentales. Un modèle plastique peut alors être suffisant si son identification est réalisée à partir d’essais où la vitesse de chargement est proche de celle rencontrée dans la structure étudiée [6].
Sinon, il est préférable de prendre un modèle visqueux, surtout si l’on souhaite simuler un chargement de type détensionnement.
Pour le type d’écrouissage, on a constaté dans [3], [4] et [6] qu’un modèle à écrouissage isotrope conduit en général à un niveau final de contrainte très élevé (comparé à des mesures). Le niveau de contraintes maximal final sera d’autant plus élevé que les cycles de plastification en tractioncompression seront nombreux. C’est pourquoi plus le nombre de passes modélisé est élevé, plus les contraintes résiduelles prévues par une modélisation avec écrouissage isotrope seront élevées. On peut atténuer ce résultat en utilisant un écrouissage isotrope avec de la restauration d’origine visqueuse, qui va compenser l’effet de l’écrouissage isotrope. Au contraire, un modèle à écrouissage cinématique a tendance à sous-estimer les contraintes résiduelles. Le comportement réel du matériau est souvent une combinaison des deux écrouissages.
4.5.4 Remarque sur les grandes déformations
La prise en compte des grandes déformations n’est en général pas indispensable pour une simulation du soudage sur structures épaisses. En revanche, ii les structures sont minces (fortes distorsions), il est nécessaire de les prendre en compte.
Dans Code_Aster, il existe deux possibilités pour tenir compte des grandes déformations dans l’opérateur STAT_NON_LINE :
•
L’option PETIT_REAC adapté lorsque les rotations sont petites ;
•
L’option SIMO_MIEHE pour une formulation exacte des grandes déformations, mais restreinte aux modèles à écrouissage isotrope.
4.5.5 Conditions aux limites
Retour d’expérience
•
Fiche 3488 : La modélisation axisymétrique du soudage sur le tube suppose implicitement à tort que le soudage a lieu simultanément sur toute la circonférence du tube, donc que la température s’élève partout dans le chanfrein. Dans la réalité, la source de chaleur progresse vers une partie de structure restée froide, qui bride obligatoirement la zone soudée. La pièce, au niveau de la source de chaleur, ne peut donc se dilater librement. Cet effet d’autobridage doit s’estomper quand la chaleur diffuse et disparaître au cours de la phase de refroidissement. Pour remédier à ce problème, on peut imposer un bridage axial au tube, uniquement dans la phase de chauffage. On empêche ainsi le tube de se dilater librement au chauffage, par contre il est libre de se déformer au refroidissement.
4.6
Temps de calcul
On constate (en version STA9), aux travers de différentes études de soudage menées ces dernières années, que l’essentiel du temps calcul est consommé dans le calcul mécanique. Pour donner un ordre de grandeur, nous donnons deux exemples ci-dessous :
•
Cas test HTNA100A qui considère un tube chanfreiné soudé en deux passes avec procédé TIG
(fiche dite 3488) :
Nombre de nœuds = 3632
Modélisation : axisymétrique
Temps CPU du calcul thermique = 307s
Temps CPU du calcul métallurgique = 25s
Temps CPU du calcul mécanique = 5169s
Manuel d'utilisation Fascicule u2.03 : Thermo-mécanique
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Notice d’utilisation pour la simulation numérique [...]
Responsable :
Jean ANGLES
Version default
Date :
10/08/2010
Page :
24/24
Clé :
U2.03.05
Révision :
4036
•
Etude sur une plaque chanfreinée soudée en deux passes avec procédé TIG (fiche dite 2425)
Nombre de nœuds = 11600
Modélisation : 3D
Temps CPU du calcul thermique = 951s
Temps CPU du calcul mécanique = 64601s
Remarques
•
Le calcul métallurgique éventuel ne coûte rien en temps calcul car il s’agit uniquement d’intégrer les modèles de comportement ; il n’y a pas d’équations d’équilibre à résoudre.
•
Le temps mécanique pourrait être certainement amélioré si on considérait (lors de la résolution des équations d’équilibre mécanique par l’opérateur STAT_NON_LINE), dans la phase de prédiction, la dépendance du module d’YOUNG et peut-être du coefficient de dilation thermique, vis à vis de la température. Actuellement, seule est prise en compte, dans le second membre, la déformation thermique (en fait la dérivée). Il existe une fiche d’évolution dans Code_Aster qui trace ce problème (pas inhérent au soudage).
•
La partie mécanique étant un problème fortement non linéaire, il est souvent nécessaire d’utiliser la recherche linéaire dans le module STAT_NON_LINE afin de faciliter la convergence. De plus, un découpage adapté du pas de temps en fonction de l’évolution des températures permet de faciliter la convergence.
Manuel d'utilisation
Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.03 : Thermo-mécanique
