Code_Aster
®
Version
5.0
Titrate:
TPLP01 Field in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER
Key
:
V4.05.001-A
Page:
1/8
Manual of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
Organization (S):
EDF/AMA, EDF/UTO/LOCATED, Delta CAD
Manual of Validation
V4.05 booklet: Stationary thermics of the plane structures
V4.05.001 document
TPLP01 - Field in L with singularity
geometrical
Summary:
This test results from the validation independent of version 3 in linear stationary thermics.
It is about a problem plane 2D presenting two modelings, one planes, the second voluminal one.
The objective is to validate, in the presence of a geometrical singularity, the plane thermal elements and 3D with
for boundary condition an imposed temperature.
The results are compared with those provided by VPCS.
Code_Aster
®
Version
5.0
Titrate:
TPLP01 Field in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER
Key
:
V4.05.001-A
Page:
2/8
Manual of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
1
Problem of reference
1.1 Geometry
0.8
0.8
0.4
0.4
With
B
C
D
E
F
0.2
0.2
Dimensions in meters
1.2
Properties of material
=
1. W/m.°C
thermal conductivity
1.3
Boundary conditions and loadings
·
side [AF] imposed Temperature T
p
= 10°C,
·
side [OF] imposed Temperature T
p
= 0°C,
·
side [AB], [BC], [CD], [EFF], flow = 0.
1.4 Conditions
initial
Without object.
Code_Aster
®
Version
5.0
Titrate:
TPLP01 Field in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER
Key
:
V4.05.001-A
Page:
3/8
Manual of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
2
Reference solution
2.1
Method of calculation used for the reference solution
The reference solution is that given in card TPLP01/89 of guide VPCS.
2.2
Results of reference
Temperature at the points of a squaring of with dimensions 0.2m X 0.2m.
2.3
Uncertainty on the solution
Analytical solution.
2.4 References
bibliographical
[1]
Guide validation of the software packages of structural analysis. French company of
Mechanics, AFNOR 1990 ISBN 2-12-486611-7
[2]
G.T. Symm-, National Physical Laboratory Division off Numerical Analysis and Computing,
treatment off singularities in solution off Laplace' S equation by integral year equation method,
NPL Carryforward NAC 31, Jan 1973.
Code_Aster
®
Version
5.0
Titrate:
TPLP01 Field in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER
Key
:
V4.05.001-A
Page:
4/8
Manual of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
3 Modeling
With
3.1
Characteristics of modeling
PLAN (QUAD8)
With
B
C
D
F
E
Limiting conditions:
- dimensioned AF:
T=10°C
- dimensioned:
T=0°C
- dimensioned AB, BC, CD, EFF:
=0
X
y
N19
N21
N23
N17
N15
N33
N35
N37
N31
N29
N41
N43
N45
3.2
Characteristics of the mesh
A number of nodes:
53
A number of meshs and types: 12 QUAD8
3.3 Functionalities
tested
Controls
AFFE_MODELE
THERMICS
PLAN
ALL
AFFE_CHAR_THER
TEMP_IMPO
THER_LINEAIRE
EXCIT
CHARGE
RECU_CHAMP
NUME_ORDRE
Code_Aster
®
Version
5.0
Titrate:
TPLP01 Field in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER
Key
:
V4.05.001-A
Page:
5/8
Manual of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
4
Results of modeling A
4.1 Values
tested
Identification Reference
Aster difference
tolerance
Temperature (°C)
X = 0.2 y = 0.0 (N15)
9.316
9.283
-
0.357
1%
y = 0.2 (N17)
9.001
9.108
1.186
1%
y = 0.4 (N19)
8.514
8.519
0.054
1%
y = 0.6 (N21)
8.018
8.015
-
0.037
1%
y = 0.8 (N23)
7.869
7.883
0.177
1%
X = 0.4
y = 0.0 (N29)
9.009
8.961
-
0.529
1%
y = 0.2 (N31)
8.640
8.669
0.334
1%
y = 0.4 (N33)
6.667
6.667
-
0.005
1%
y = 0.6 (N35)
5.680
5.666
-
0.254
1%
y = 0.8 (N37)
5.495
5.519
0.443
1%
X = 0.6
y = 0.4 (N41)
2.972
2.963
-
0.310
1%
y = 0.6 (N43)
2.881
2.877
-
0.132
1%
y = 0.8 (N45)
2.816
2.834
0.650
1%
4.2 Parameters
of execution
Version: 5.03
Machine: SGI - ORIGIN 2000 - R12000
Overall dimension memory:
8 megawords
Time CPU To use: 2.25 seconds
Code_Aster
®
Version
5.0
Titrate:
TPLP01 Field in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER
Key
:
V4.05.001-A
Page:
6/8
Manual of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
5 Modeling
B
5.1
Characteristics of modeling
3D (HEXA8)
With
B
C
D
F
E
Limiting conditions:
- face AGLF:
T=10°C
- face DJKE:
T=0°C
- faces ABHG, BHIC:
=0
- faces CDJI, FEKL:
=0
X
y
Thickness = 0.2 m
G
H
I
J
L
K
Z
N12
N14
N15
N17
N20
N11
N13
N16
N18
N19
N22
N23
N21
N24
N25
N28
N29
N26
N27
N30
N32
N33
N36
N31
N34
N35
5.2
Characteristics of the mesh
A number of nodes:
42
A number of meshs and types: 12 HEXA8
5.3 Functionalities
tested
Controls
AFFE_MODELE
THERMICS
3D
ALL
AFFE_CHAR_THER
TEMP_IMPO
THER_LINEAIRE
EXCIT
CHARGE
RECU_CHAMP
NUME_ORDRE
Code_Aster
®
Version
5.0
Titrate:
TPLP01 Field in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER
Key
:
V4.05.001-A
Page:
7/8
Manual of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
6
Results of modeling B
6.1 Values
tested
Identification Reference
Aster difference
tolerance
Temperature (°C)
X = 0.2 y = 0.0 (N11)
9.316
9.294
-
0.239
1%
(N12)
9.316
9.294
-
0.239
1%
y = 0.2 (N13)
9.001
9.015
0.160
1%
(N14)
9.001
9.015
0.160
1%
y = 0.4 (N15)
8.514
8.505
-
0.108
1%
(N16)
8.514
8.505
-
0.108
1%
y = 0.6 (N17)
8.018
8.025
0.093
1%
(N18)
8.018
8.025
0.093
1%
y = 0.8 (N19)
7.869
7.861
-
0.096
1%
(N20)
7.869
7.861
-
0.096
1%
X = 0.4
y = 0.0 (N21)
9.009
8.996
-
0.139
1%
(N22)
9.009
8.996
-
0.139
1%
y = 0.2 (N23)
8.640
8.661
0.247
1%
(N24)
8.640
8.661
0.247
1%
y = 0.4 (N25)
6.667
6.667
-
0.005
1%
(N26)
6.667
6.667
-
0.005
1%
y = 0.6 (N27)
5.680
5.669
-
0.188
1%
(N28)
5.680
5.669
-
0.188
1%
y = 0.8 (N29)
5.495
5.502
0.123
1%
(N30)
5.495
5.502
0.123
1%
X = 0.6
y = 0.4 (N31)
2.972
2.990
0.621
1%
(N32)
2.972
2.990
0.621
1%
y = 0.6 (N33)
2.881
2.959
2.712
1%
(N34)
2.881
2.959
2.712
1%
y = 0.8 (N35)
2.816
2.845
1.024
1%
(N36)
2.816
2.845
1.024
1%
6.2 Parameters
of execution
Version: 5.03
Machine: SGI - ORIGIN 2000 - R12000
Overall dimension memory:
8 megawords
Time CPU To use: 2.00 seconds
Code_Aster
®
Version
5.0
Titrate:
TPLP01 Field in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER
Key
:
V4.05.001-A
Page:
8/8
Manual of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
7
Summary of the results
Two modelings give results whose certain values exceed the fixed tolerance
initially (1%):
·
for modeling A (PLANE with meshs QUAD8), the maximum change is 1.19%
(going beyond for only one value out of the 13 tested),
·
for modeling B (3D with meshs HEXA8), the maximum change is 2.7%
(going beyond for two values out of 26 tested).
The modeling of the geometrical singularity (presence of an important heat gradient close to
singularity) is represented better with quadratic elements (modeling A).
For two modelings, the precision should be improved by using a finer mesh
(more important refinement in the area of the geometrical singularity).
The results are regarded as acceptable taking into account modelings carried out.