TP Couplage Thermo-Mécanique-Électrique

Chauffage d’un flan rectangulaire par effet Joule

Mécatronique 4ème année − S8
Hervé LAURENT

L’objectif de ce TP est de simuler la chauffe par effet Joule d’un flan rectangulaire sur Abaqus
[1] (éventuellement également sur Comsol MultiPhysics® [2]). Le moyen de chauffage est
présenté sur la figure 1.

Il permet de chauffer un flan rectangulaire d’un acier trempant développé par ArcelorMittal
(appelé Usibor® 1500P) en quelques dizaines de secondes. Les résultats des deux codes seront
comparés avec les températures expérimentales du fichier en annexe.

Figure 1 – Schéma du prototype de chauffage d’un flan rectangulaire et photo d’un essai de
chauffe.

1
1 Couplage thermo-mécanique-électrique
Le couplage étudié sera ici un couplage fort entre la partie électrique-thermique et mécanique,
c’est-à-dire que le courant électrique va chauffer la tôle par effet Joule qui va ensuite modifier le
champ de contrainte dans la tôle. Il sera donc possible de simuler une "sollicitation" mécanique
de la tôle après cette chauffe comme par exemple, faire un essai de traction ou d’emboutissage.

L’équation de la thermique s’écrira dans ce cas :

∂T
ρCp = div [k · grad(T )] + QJoule (1)
∂t
où :
- T (⃗x, t) est la température fonction de la position ⃗x et du temps t,
- ρ la masse volumique,
- Cp la chaleur spécifique,
- k la conductivité thermique (supposé isotrope).
QJoule est la source de chaleur induite par l’effet Joule donnée par la conservation de la charge
électrique :  
QJoule = div J⃗ = σe . grad ⃗ (V ) (2)
∼

où :
- σe est le tenseur de conductivité électrique (supposé isotrope),
∼

- J⃗ le vecteur de densité électrique,

- V le potentiel électrique.

2
2 Modélisation éléments finis
Seul un quart de la géométrie sera représenté en utilisant les conditions de symétrie du problème
(voir figure 2).

Figure 2 – Conditions aux limites sur le 1/4 d’un flan rectangulaire chauffé par le passage
d’un courant électrique (longueur, largeur et épaisseur du flan : L=252,5 ; b=40 et h=1,5 mm).

Des éléments hexaèdres de type Q3D8 pour « Thermal Electrical Structural » seront à utiliser.
Ils prennent donc en compte les degrés de liberté thermique (ddl 11), mécanique (ddl de 1 à 6)
et électrique (ddl 9).
Les calculs se feront en unité SI : m, N, Pa. Il est aussi préférable d’utiliser les K. Il faudra
attribuer les constantes physiques suivantes au modèle d’Abaqus en cliquant sur Model-1 et
en utilisant « Model-Edit Attributes » et en renseignant la température absolue à 0K et la
constante de Stephan-Boltzman à : σ=5,67 × 10−8 W/m2 /K4

3 Matériaux
3.1 Usibor® 1500P
L’ensemble des propriétés de cet acier subit une forte variation aux alentours de 750 °C, ceci est
dû à un changement d’état métallurgique (austénitisation). Ces données étant confidentielles,
on supposera des propriétés constantes qui sont données dans le tableau 1.

ρ [kg m−3 ] k [W m−1 K−1 ] Cp [J kg−1 K−1 ] σe [Ω m] E [Pa] ν

7940. 40. 450. 4,2 × 106 2,11 × 1011 0.3

Table 1 – Propriétés de l’Usibor à température ambiante.

3
3.2 Cuivre
Les propriétés des électrodes en cuivre sont celles issues de la base de données de COMSOL
Multiphysics® , elles sont présentées dans le tableau 2.

ρ [kg m−3 ] k [W m−1 K−1 ] Cp [J kg−1 K−1 ] σe [Ω m] E [Pa] ν

8700. 400. 385. 6 × 107 1,24 × 1011 0.3

Table 2 – Propriétés du cuivre à température ambiante.

4 Conditions aux limites pour la thermique

Les conditions aux limites thermiques et électriques utilisées pour effectuer cette simulation
sont présentées sur la figure 2.

4.1 Convection
Les pertes convectives seront décrites par la loi de Newton :

−⃗n . φ
⃗ conv = h (Tair − T ) (3)

où ⃗n est le vecteur normal à la surface et φ

⃗ conv est le vecteur flux de chaleur par convection.

Les coefficients d’échange convectif appliqués dans le modèle hinf , hsup , hcote et helec sont sup-
posés constants (paramètres à ajuster) :
- hinf =5,3 W m−2 K−1 pour la face inférieure du flan
- hsup =10,6 W m−2 K−1 pour la face supérieure du flan
- hcote =20,8 W m−2 K−1 pour les faces latérales du flan
- helec =7 W m−2 K−1 pour les faces de l’électrode

4.2 Rayonnement
Les échanges radiatifs avec l’environnement sont pris en compte par la loi de Stefan-Boltzman :
4
− T4 (4)


−⃗n . φ
⃗ rad = ε . σb Tenv

avec une température d’environnement Tenv =20 °C et où φ

⃗ rad est le vecteur flux de chaleur
par rayonnement.

Les valeurs d’émissivité seront dans un premier temps constantes, puis l’influence des émissivités
suivantes seront analysées :
- émissivité de ε = 0, 35 pour les faces supérieures et inférieures du flan.
- émissivité de ε = 0, 8 pour les faces latérales du flan.
- émissivité de ε = 0, 65 pour les électrodes en cuivre.
- Pas d’échanges radiatifs pour les faces en vis-à-vis et/ou espaces confinés

4
4.3 Température imposée
Les bases des électrodes étant loin de la source de chaleur, une température Telec =20 °C
(température de l’air ambiant) est imposée.

4.4 Isolation thermique

Des flux de chaleur nuls sont imposés au niveau des symétries.

5 Conditions aux limites pour la partie électrique

5.1 Isolation électrique
Des conditions d’isolation électrique (flux nuls) sont imposées sur toutes les faces extérieures,
car en contact avec de l’air.

5.2 Potentiel électrique/densité de courant

Un potentiel nul est imposé au centre du flan (potentiel de référence). Comme pour l’essai
expérimental, le chauffage du flan est contrôlé par l’intensité imposée. En conséquence, une
densité de courant sera imposée numériquement sur les cotés des électrodes sous la forme d’un
fichier (*AMPLITUDE) donné en annexe.
**
** LOADS
**
*INCLUDE,INPUT=Intensite_d.inp
** Name: DensiteCourant Type: Surface current
*Dsecurrent, amplitude=INTENSITE
** Calculer la surface equivalente de l’electrode
Electrode-SUP-1.Surf-Cote-Electrode,CS,250.
Le coefficient multiplicateur de « 250 » correspond ici à l’inverse de la valeur de l’aire de la
surface où est appliquée le courant.

6 Résistance de Contact Thermique/Electrique

Aux interfaces électrode/flan, les résistances de contact électrique (RCE) et des résistances
de contact thermiques (RCT) ne peuvent être négligées. Elles dépendent de la pression, de la
température et des matériaux en présence.

6.1 Résistance de Contact Thermique

Le transfert de chaleur entre 2 surfaces en contact sera simulé par :
−⃗n . φ
⃗ cont = RCT · (T1 − T2 ) (5)
où φ
⃗ cont est le vecteur flux de chaleur, RCT le coefficient de Résistance de Contact Thermique
et T1 et T2 les températures entre les surfaces en contact.

5
6.2 Résistance de Contact Électrique
La densité électrique de chaleur entre 2 surfaces en contact sera simulée par :

−⃗n . I⃗cont = RCE · (V1 − V2 ) (6)

où I⃗cont est le vecteur densité électrique, RCE le coefficient de Résistance de Contact Electrique
et V1 et V2 le potentiel électrique entre les surfaces en contact.

Utiliser, dans un premier temps, les valeurs suivantes pour ces résistances :
RCT =1/15 × 10−4 =666 W m−2 K−1 et RCE=1/100 × 10−8 =1 × 106 Ω−1 m−2 :

**
** INTERACTION PROPERTIES
**
*Surface Interaction, name=Int-Contact-Meca-Electric
1.,
*Friction
0.,
*Gap Conductance
666., 0.
666., 1.e+06
*Gap Electrical Conductance
1.e+06, 0.
1.e+06, 1.e+06

Ensuite, il s’agira de les faire évoluer pour tendre vers la solution expérimentale.

6
7 Travail à réaliser
La figure 3 donne la position des thermocouples (TC1 à TC6) pendant les essais.

Figure 3 – Positions des thermocouples sur un flan rectangulaire.

1. Comparer les températures expérimentales et numériques sur les thermocouples (TC1,

TC2 et TC3)
2. Regarder l’influence des RCE/RCT sur les évolutions des températures.
3. Regarder l’influence des coefficients thermiques/matériaux sur ces évolutions.
4. Pour les plus courageux, comparer les résultats entre Abaqus et Comsol sur cet exemple.

7
A Input file template
Un exemple de mise en données de calcul thermo-mécanique-électrique Abaqus est proposée
ci-dessous :

*HEADING
...
** Specify the coupled thermal-electrical-structural element type
*ELEMENT, TYPE=Q3D8
...
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0., STEFAN BOLTZMANN=5.67E-08
*STEP
*COUPLED TEMPERATURE-DISPLACEMENT, ELECTRICAL
Data line to define incrementation
*BOUNDARY
Data lines to define nonzero boundary conditions on displacement,
temperature or electrical potential degrees of freedom
*CFLUX and/or *CFILM and/or
*CRADIATE and/or *DFLUX and/or
*DSFLUX and/or *FILM and/or
*SFILM and/or *RADIATE and/or
*SRADIATE
Data lines to define thermal loads
*CLOAD and/or *DLOAD and/or *DSLOAD
Data lines to define mechanical loads
*CECURRENT
Data lines to define concentrated currents
*DECURRENT and/or *DSECURRENT
Data lines to define distributed current densities
*FIELD
Data lines to define field variable values
*END STEP

Références
[1] Abaqus. ABAQUS Documentation Version 2020, 2020. http ://abaqus-doc.irdl.fr/.

[2] Comsol MultiPhysics. Comsol Multyphysics® Documentation, 2020.