Code_Aster
®
Version
5.0
Titrate:
Thermal TTLP300 Transfer in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND,
E. SCREWS, F. LEBOUVIER
Key
:
V4.23.300-A
Page:
1/6
Manual of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
Organization (S):
EDF/AMA, EDF/UTO/LOCATED, Delta CAD
Manual of Validation
V4.23 booklet: Transitory thermics of the plane systems
V4.23.300 document
TTLP300 - Thermal transfer 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.
Code_Aster
®
Version
5.0
Titrate:
Thermal TTLP300 Transfer in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND,
E. SCREWS, F. LEBOUVIER
Key
:
V4.23.300-A
Page:
2/6
Manual of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
1
Problem of reference
1.1 Geometry
has
has
B
B
X
,
y
,
, C
(0) =
y
X
H T
F
= 50.8 X 10 have
- 3
m
B = 25.4 X 10
- 3
m
C
B
With
D
1.2
Properties of material
X
= 34.614 W/m °C
thermal conductivity along axis X
y
= 6.237 W/m °C
thermal conductivity along the axis y
C
p
= 37.719
J/kg. °C
specific heat
= 6407.38
kg/m
3
mass
voluminal
1.3
Boundary conditions and loadings
Convection:
·
H = 1362.71 W/m
2
°C,
·
T
F
= 37.78 °C.
1.4 Conditions
initial
T (X, y, T = 0) = 260° C
Code_Aster
®
Version
5.0
Titrate:
Thermal TTLP300 Transfer in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND,
E. SCREWS, F. LEBOUVIER
Key
:
V4.23.300-A
Page:
3/6
Manual of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
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 manual 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”, 1
St
edition, June 1, 1976
Code_Aster
®
Version
5.0
Titrate:
Thermal TTLP300 Transfer in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND,
E. SCREWS, F. LEBOUVIER
Key
:
V4.23.300-A
Page:
4/6
Manual of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
3 Modeling
With
3.1
Characteristics of modeling
PLAN (QUAD8)
Y
X
C
B
With
D
Limiting conditions:
- dimensioned AB and DA:
= 0
- dimensioned BC and CD: T
ext.
= 37.78 °C
H = 1362.71 W/m ² °C
= 30°
X
y
N1
N26
N113
N82
3.2
Characteristics of the mesh
A number of nodes:
113
A number of meshs and types: 30 QUAD8
3.3 Functionalities
tested
Controls
AFFE_MODELE
THERMICS
PLAN
ALL
DEFI_MATERIAU
THER_ORTH
AFFE_CARA_ELEM
SOLID MASS
AFFE_CHAR_THER_F
EXCHANGE
THER_LINEAIRE
TEMP_INIT
VALE
EXCIT
CARA_ELEM
RECU_CHAMP
INST
3.4 Remarks
The discretization in pitch of time is as follows:
10 pitches for [0. , 1.D-4] is
T = 1.D-5
9 pitches for [1.D-4, 1.D-3] are
T = 1.D-4
9 pitches for [1.D-3, 1.D-2] are
T = 1.D-3
9 pitches for [1.D-2, 1.D-1] are
T = 1.D-2
9 pitches for [1.D-1, 1.D0] are
T = 1.D-1
20 pitches for [1.D+0, 3.D0] are
T = 1.D-1
Code_Aster
®
Version
5.0
Titrate:
Thermal TTLP300 Transfer in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND,
E. SCREWS, F. LEBOUVIER
Key
:
V4.23.300-A
Page:
5/6
Manual of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
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
Overall dimension memory:
8 megawords
Time CPU To use: 4.70 seconds
Code_Aster
®
Version
5.0
Titrate:
Thermal TTLP300 Transfer in an orthotropic metal bar
Date:
20/09/02
Author (S):
C. DURAND,
E. SCREWS, F. LEBOUVIER
Key
:
V4.23.300-A
Page:
6/6
Manual of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A
5
Summary of the results
The results obtained are satisfactory, the maximum change is 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 controls:
·
DEFI_MATERIAU
associated the key word
THER_ORTH
, allowing to define the characteristics
of an orthotropic material,
·
AFFE_CARA_ELEM
associated the key word
SOLID MASS
, allowing to define the axes of orthotropism.