Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
1/10

Organisme(s) : EDF/AMA














Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
Document : V4.01.001





TPLA01 - Cylindre creux infini en équilibre
thermique





Résumé :

Thermique stationnaire linéaire.

Modèle axisymétrique ; 3 modélisations.

Solution analytique.

Intérêt du test :

· tous les éléments axisymétriques : triangles et quadrangles, degrés 1 et 2,
· conditions aux limites variées : échange, température imposée, flux imposé,
· validation partielle de la matrice de "masse" thermique car on fait un "faux" transitoire.

Les résultats ne sont pas affectés par la distorsion des mailles h/l = 40.

Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
2/10


1
Problème de référence

1.1 Géométrie

z
Ri
Re
Rayon intérieur
Ri = 0.300 m
Rayon extérieur
Re = 0.35 m
Point F

r = 0.32 m
r
z
D
C
E
F
A
B
r


1.2
Propriétés de matériaux

= 1 W/m °C

Cp = 2. J/m3 °C (chaleur volumique)

1.3
Conditions aux limites et chargements

· [DC] [AB] :
= 0 W/m2
· [EA] :
T = Ti = 100°C
· [ED] :
= i = 1729.9091 W/m2 (flux rentrant)
· [CB] : échange
h = he = 500 W/m2 °C

T = Te = 17.03444 °C

1.4 Conditions
initiales

Pour faire ce calcul stationnaire, on fait un calcul transitoire pour lequel les conditions aux limites sont
constantes dans le temps. Ceci permet de tester les calculs élémentaires de masse intervenant dans
le premier membre ainsi que le deuxième membre.
Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
3/10


2
Solution de référence

2.1
Méthode de calcul utilisée pour la solution de référence

r
T(r) = Ti + Log
Ri

T - T
T
e
i
i : température en peau " interne"
=
;
R
T


e
e : température en peau " externe"
Log


Ri


T

avec : les flux "radiaux"
sur les parois du cylindre sont :
r


i
= + ,
R

i


e
= +
.
R

e



2.2
Résultats de référence

Températures et flux aux points A, B, D, F.



2.3
Incertitude sur la solution

Solution analytique.

Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
4/10


3 Modélisation
A

3.1
Caractéristiques de la modélisation

axis (TRIA3, QUAD4)


y
l = 2.5 10-3
D
C
x y

h = 0.1
A 0.3 0.0 N1
B 0.35
0.0
N41
E
G
D
0.30 0.10 N43
E
0.30 0.05 N2
F 0.32
0.0
N17
j
0
x
i
A
F
B



3.2
Caractéristiques du maillage

Nombre de noeuds : 63.

Nombre de mailles et types : 40 TRIA3, 20 QUAD4


3.3 Fonctionnalités
testées

Commandes



AFFE_MODELE
THERMIQUE
AXIS

THER_LINEARE
TEMP_INIT
STATIONNAIRE

EXCIT
CHARGE

AFFE_CHAR_THER
TEMP_IMPO
FLUX_REP
ECHANGE


Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
5/10


4
Résultats de la modélisation A

4.1 Valeurs
testées


Identification Référence
Aster %
différence
T(A) 100 100.000 0.00
T(B) 20 20.002 0.01
T(F) 66.506
66.507 0.00
T(D) 100 99.933 ­0.07
(A) 1729.91 1722.62 ­0.42
(B) 1482.78 1487.98
0.35
(D) 1729.91 1723.26 ­0.38
(F) 1621.79 1628.17
0.39
(F) 1621.79 1615.50 ­0.39

Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
6/10


5 Modélisation
B

5.1
Caractéristiques de la modélisation

axis (TRIA6, QUAD8)


y
l = 2.5 10-3
D
C
x y

h = 0.1
A 0.3 0.0 N180
B 0.35
0.0
N10
E
D
0.30 0.10 N178
E
0.30 0.05 N183
F 0.32
0.0
N112
0
x
A
F
B



5.2
Caractéristiques du maillage

Nombre de noeuds : 185.

Nombre de mailles et types : 40 TRIA6, 20 QUAD8


5.3 Fonctionnalités
testées

Commandes



AFFE_MODELE
THERMIQUE
AXIS

THER_LINEARE
TEMP_INIT
STATIONNAIRE

EXCIT
CHARGE

AFFE_CHAR_THER
TEMP_IMPO
FLUX_REP
ECHANGE


Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
7/10


6
Résultats de la modélisation B

6.1 Valeurs
testées


Identification Référence
Aster %
différence
T(A) 100 100.00
0.00
T(B) 20 20.00 0.00
T(F) 66.506
66.506 0.00
T(D) 100 100.00
0.00
(A) 1729.91 1729.89
0.00
(B) 1482.78 1482.77
0.00
(D) 1729.91 1729.88
0.00
(F) 1621.79 1621.77
0.00
(F) 1621.79 1621.78
0.00


Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
8/10


7 Modélisation
C

7.1
Caractéristiques de la modélisation

axis (TRIA6, QUAD9)


y
l = 2.5 10-3
D
C
x y

h = 0.1
A 0.3 0.0 N199
B 0.35
0.0
N10
E
D
0.30 0.10 N197
E
0.30 0.05 N203
F 0.32
0.0
N124
0
x
A
F
B



7.2
Caractéristiques du maillage

Nombre de noeuds : 205.

Nombre de mailles et types : 40 TRIA6, 20 QUAD9


7.3 Fonctionnalités
testées

Commandes



AFFE_MODELE
THERMIQUE
AXIS

THER_LINEARE
TEMP_INIT
STATIONNAIRE

EXCIT
CHARGE

AFFE_CHAR_THER
TEMP_IMPO
FLUX_REP
ECHANGE


Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
9/10


8
Résultats de la modélisation C

8.1 Valeurs
testées


Identification Référence
Aster %
différence
T(A) 100 100.00
0.00
T(B) 20 20.00 0.00
T(F) 66.506
66.506 0.00
T(D) 100 100.00
0.00
(A) 1729.91 1729.89
0.00
(B) 1482.78 1482.77
0.00
(D) 1729.91 1729.88
0.00
(F) 1621.79 1621.77
0.00
(F) 1621.79 1621.77
0.00


Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Code_Aster ®
Version
5.0

Titre :

TPLA01 - Cylindre creux infini en équilibre thermique


Date :
30/08/02
Auteur(s) :
J. PELLET Clé
:
V4.01.001-C Page :
10/10


9
Synthèse des résultats

Ce problème est correctement résolu :

· avec les différents types d'éléments quel que soit le degré d'interpolation,
· n'est pas affecté par la forme des éléments h/l = 40.
Manuel de Validation
Fascicule V4.01 : Thermique stationnaire des structures axisymétriques
HT-66/02/001/A

Document Outline