Code_Aster ®
Version
5.0
Titrate:
Thermal TTLP300 Transfert in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.300-A Page:
1/6
Organization (S): EDF/AMA, EDF/UTO/LOCATED, Delta CAD
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
V4.23.300 document
TTLP300 - Thermal Transfert in a bar
metal orthotropic
Summary:
This test results from the validation independent of version 3 in linear transitory thermics.
It is about a problem plane 2D represented by only one modeling (plane).
The functionalities tested are as follows:
· plane thermal element,
· orthotropic material,
· transitory algorithm of thermics,
· limiting condition: convection.
The interest of the test lies in the taking into account of an orthotropic material.
The results are compared with a solution based on a graphic estimate.
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
Thermal TTLP300 Transfert in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.300-A Page:
2/6
1
Problem of reference
1.1 Geometry
y
has
has
= 50.8 X 10-3m have
B = 25.4 X 10-3m
D
C
B
X
With
B
(0) =
B
X, y, C
H Tf
1.2
Properties of material
X = 34.614 W/m °C conductivity thermal along axis X
y = 6.237 W/m °C
thermal conductivity along the axis y
CP = 37.719 J/kg. °C
specific heat
= 6407.38
kg/m3 masses
voluminal
1.3
Boundary conditions and loadings
Convection:
· H = 1362.71 W/m2 °C,
· Tf = 37.78 °C.
1.4 Conditions
initial
T (X, y, T = 0) = 260° C
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
Thermal TTLP300 Transfert in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.300-A Page:
3/6
2
Reference solution
2.1
Method of calculation used for the reference solution
The original reference solution given in the book [bib1] is based on a graphic estimate.
This reference is quoted in the handbook of checking of ANSYS [bib2]
2.2
Results of reference
Temperature at points A B C D at the moment T = 3s
2.3
Uncertainty on the solution
Unknown factor, it was not possible to get the original reference (delivers old, more published).
2.4 References
bibliographical
[1]
Schneider, P.J., “Conduction Heat Transfer”, Addison-Wesley Publishing Co., Inc. Reading,
Farmhouse., 2nd Printing, 1957.
[2]
ANSYS: “Checking manual”, 1st edition, June 1, 1976
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
Thermal TTLP300 Transfert in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.300-A Page:
4/6
3 Modeling
With
3.1
Characteristics of modeling
PLAN (QUAD8)
Y
N82
C
Limiting conditions:
- dimensioned AB and DA: = 0
- dimensioned BC and CD: Text = 37.78 °C
H = 1362.71 W/m ² °C
N113
D
B
y
N1
X
= 30°
X
AN26
3.2
Characteristics of the grid
A number of nodes:
113
A number of meshs and types: 30 QUAD8
3.3 Functionalities
tested
Commands
AFFE_MODELE
THERMIQUE
PLAN
TOUT
DEFI_MATERIAU
THER_ORTH
AFFE_CARA_ELEM
MASSIF
AFFE_CHAR_THER_F
ECHANGE
THER_LINEAIRE
TEMP_INIT
VALE
EXCIT
CARA_ELEM
RECU_CHAMP
INST
3.4 Remarks
The discretization in step of time is as follows:
10 steps for [0. , 1.D-4] is T = 1.D-5
9 steps for [1.D-4, 1.D-3] are T = 1.D-4
9 steps for [1.D-3, 1.D-2] are T = 1.D-3
9 steps for [1.D-2, 1.D-1] are T = 1.D-2
9 steps for [1.D-1, 1.D0] are T = 1.D-1
20 steps for [1.D+0, 3.D0] are T = 1.D-1
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
Thermal TTLP300 Transfert in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.300-A Page:
5/6
4
Results of modeling A
4.1 Values
tested
Identification Reference
Aster
% difference
Tolerance
moment T = 3 S
Points T (°C)
With (N26)
237.50
238.95
0.611
5%
B (N113)
137.22
140.71
2.541
5%
C (N82)
65.98
66.19
0.318
5%
D (N1)
94.44
93.30
- 1.206
5%
4.2 Parameters
of execution
Version: 5.03
Machine: SGI - ORIGIN 2000 - R12000
Obstruction memory:
8 megawords
Time CPU To use: 4.70 seconds
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
Thermal TTLP300 Transfert in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.300-A Page:
6/6
5
Summary of the results
The results obtained are satisfactory, the maximum change is of 2.5% lower than the fixed tolerance
initially (5%) (the reference solution is obtained graphically).
This test made it possible to test in linear transient, modeling PLAN, the commands:
· DEFI_MATERIAU associated with key word THER_ORTH, allowing to define the characteristics
of an orthotropic material,
· AFFE_CARA_ELEM associated with key word MASSIF, allowing to define the axes of orthotropism.
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
Outline document