Temperatura di una parete multistrato

Consideriamo un problema di stato stazionario come il flusso di calore in una piastra di spessore h con conducibilità termica k, le cui superfici hanno differenti temperature t1 e t2 (si veda figura).

La variazione di temperatura lungo lo spessore della piastra h è definita dalla relazione:

Quindi, il flusso di calore in ogni punto è uguale a:

Ora, supponiamo che la piastra è un composito, che è, costituito da n strati con spessore h1,h2,...,hn e coefficienti di conducibilità termica k1,k2,...,kn, rispettivamente. Quindi, il flusso di calore per ogni strato f1,f2,...,fn si ricava dalla seguente formula:

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

Se gli strati hanno un contatto termico ideale tra le interfacce, allora il flusso di calore sarà continuo quando si passa da uno strato all'altro, e per questo problema particolare, sarà la stessa in ogni punto (cioè, f1 = f2 = ... = fn = f). La variazione di temperatura tra le superfici esterne opposte l'intera piastra composita è uguale alla somma delle variazioni di temperatura in ogni singolo strato:

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

Poi:

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

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

Usiamo i seguenti dati: numero di strati n = 3, lunghezza 0.5 m e larghezza 0.3 m, spessori degli strati h1,h2,h3 sono pari a 0.007 m, 0.01 m e 0.003 m. Le temperature applicate t1 e t4 sono uguali a 273.15 K (o 0 oC) e 373.15 K (o 100 oC) rispettivamente.

Coefficienti di conducibilità termica: k1 = 200 W/m.K, k2 = 390 W/m.K, k3 = 43 W/m.K .

Il modello ad elementi finiti con temperature applicate

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

Dopo aver effettuato il calcolo con l'aiuto di AutoFEM, i questi risultati sono ottenuti:

Tabella 1. I parametri della maglia ad elementi finiti

Tipo di elemento finito

Numero di nodi

Numero di elementi finiti

tetraedro lineare

792

3063

Tabella 2. Risultato "Temperatura"

Superficie Sij di separazione di strati i e j

Soluzione numerica
Temperatura T*, K

Soluzione analitica
Temperatura T, К

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

S12

2.999887E+02

2.999887E+02

0.01

S23

3.196508E+02

3.196508E+02

0.01

Tabella 3. Risultato "Flusso di calore"

Soluzione numerica
Flusso di calore f *, W/m2

Soluzione analitica
Flusso di calore f , W/m2

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

7.66821375E+005

7.66821372E+005

3.91E-007

 

 

* I risultati dei test numerici dipendono dalla maglia ad elementi finiti e possono differire leggermente da quelli indicati nella tabella.

 

 

 

Torna ai contenuti

autofem.com