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