Température d'équilibre d'une paroi multicouche

Considérons un problème de l'état stable de flux de chaleur dans une plaque d'épaisseur h avec la conductivité thermique k, dont la surface est maintenu à une température t1 et t2  (voyez figure).

La variation de température le long de l'épaisseur de la plaque est définie par la relation:

Donc, le flux de chaleur en tout point est égale à:

Maintenant, supposons que la plaque est un composite, qui est, est constitué de n couches d'épaisseurs h1,h2,...,hn et les coefficients de conductivité thermique k1,k2,...,kn, respectivement. Ensuite, le flux de chaleur pour chaque couche f1,f2,...,fn peut être trouvée à partir de la formule:

fi = - ki (ti+1 - ti) / hi = (ti - ti+1) / Ri ,
Ri = hi / ki ,
ti < ti+1 , i = 1,2,...,n

Laissons les couches ont idéale contact thermique à travers les interfaces, puis le flux de chaleur sera continue lorsqu'on passe d'une couche à l'autre, et pour ce problème particulier, ce sera la même en tout point (c'est-à, f1 = f2 = ... = fn = f). La variation de température entre les surfaces opposées extérieures de l'ensemble de plaque de composite est égale à la somme des variations de température dans chaque couche unique:

(t1 - t2) + (t2 - t3) + ... + (ti - ti+1) + ...+ (tn - tn+1) = t1 - tn+1

Ensuite:

t1 - tn+1 = f1 R1 + f2 R2 + ...+fn Rn = f ( R1 + R2 + ...+ Rn ),

f = ( t1 - tn+1 ) / (R1 + R2 + ...+ Rn) .

Nous laisser utiliser les données suivantes: nombre de couches n = 3, la longueur et la largeur de chaque couche est 0.5 m et 0.3 m , respectivement, les épaisseurs des couches h1,h2,h3 sont égales à 0.007 m, 0.01 m et 0.003 m . Températures appliquées t1 et t4 sont égales à 273.15 K (ou 0 oC) et 373.15 K (ou 100 oC) respectivement.

Les coefficients de conductivité thermique: k1 = 200 W / m.K, k2 = 390 W / m.K, k3 = 43 W / m.K .

Le modèle éléments finis avec des températures appliquées

Donc, f = - 7.6682E+005 W / m2, t2 = -R1 f + t1 = 299.9887 K, t3 = - (R1 + R2) f + t1 = 319.6508 K.

Après avoir effectué le calcul à l'aide de AutoFEM, les résultats suivants sont obtenus:

Tableau 1. Paramètres de maillage éléments finis

Type d'élément fini

Nombre de nœuds

Nombre d'éléments finis

tétraèdre linéaire

792

3063

Tableau 2. Résultat "Temperature"

Surface Sij de séparation des couches i et j

Solution numérique
Temperature T*, K

Solution analytique
Temperature T, К

Erreur δ = 100%* |T* - T| / |T|

S12

2.999887E+02

2.999887E+02

0.01

S23

3.196508E+02

3.196508E+02

0.01

Tableau 3. Résultat "Flux de chaleur"

Solution numérique
Flux de chaleur f *, W/m2

Solution analytique
Flux de chaleur f , W/m2

Erreur δ = 100%* | f* - f | / | f |

7.66821375E+005

7.66821372E+005

3.91E-007

 

 

*Les résultats des tests numériques dépendent du maillage éléments finis et peuvent différer légèrement de celles indiquées dans le tableau.

 

En savoir plus sur AutoFEM Thermal Analysis

Retour au sommaire

autofem.com