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