Code_Aster ®
Version
5.0
Titrate:
TPLP01 Domaine in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.05.001-A Page:
1/8
Organization (S): EDF/AMA, EDF/UTO/LOCATED, Delta CAD
Handbook of Validation
V4.05 booklet: Stationary thermics of the plane structures
V4.05.001 document
TPLP01 - Domaine 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.
Handbook of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPLP01 Domaine in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.05.001-A Page:
2/8
1
Problem of reference
1.1 Geometry
0.8
F
E
0.4
0.2
0.8
C
D
0.2
With
B
0.4
Dimensions in meters
1.2
Properties of material
=
1. W/m.°C
thermal conductivity
1.3
Boundary conditions and loadings
· side [AF] Température imposed Tp = 10°C,
· side [OF] Température imposed Tp = 0°C,
· side [AB], [BC], [CD], [EFF], flow = 0.
1.4 Conditions
initial
Without object.
Handbook of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPLP01 Domaine in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.05.001-A Page:
3/8
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.
Handbook of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPLP01 Domaine in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.05.001-A Page:
4/8
3 Modeling
With
3.1
Characteristics of modeling
PLAN (QUAD8)
y
F
E
N23
N37
N45
Limiting conditions:
N21
N35
N43
- dimensioned AF:
T=10°C
- dimensioned:
T=0°C
- dimensioned AB, BC, CD, EFF: =0
N19
N33
N41
C
D
N17
N31
X
N15
N29
With
B
3.2
Characteristics of the grid
A number of nodes:
53
A number of meshs and types: 12 QUAD8
3.3 Functionalities
tested
Commands
AFFE_MODELE
THERMIQUE
PLAN
TOUT
AFFE_CHAR_THER
TEMP_IMPO
THER_LINEAIRE
EXCIT
CHARGE
RECU_CHAMP
NUME_ORDRE
Handbook of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPLP01 Domaine in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.05.001-A Page:
5/8
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
Obstruction memory:
8 megawords
Time CPU To use: 2.25 seconds
Handbook of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPLP01 Domaine in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.05.001-A Page:
6/8
5 Modeling
B
5.1
Characteristics of modeling
3D (HEXA8)
y
L
N19
N30
N35
K
Limiting conditions:
N29
N36
F
N20
E
- face AGLF:
T=10°C
N18
N27
N34
- face DJKE:
T=0°C
- faces ABHG, BHIC:
=0
N17
N28
N33
- faces CDJI, FEKL:
=0
N16
I
N26
N31
J
N15
N25
N32
N13
C N24
D
Thickness = 0.2 m
N14
N23
G
N11
N21
H
N12
N22
X
With
B
Z
5.2
Characteristics of the grid
A number of nodes:
42
A number of meshs and types: 12 HEXA8
5.3 Functionalities
tested
Commands
AFFE_MODELE
THERMIQUE
3D
TOUT
AFFE_CHAR_THER
TEMP_IMPO
THER_LINEAIRE
EXCIT
CHARGE
RECU_CHAMP
NUME_ORDRE
Handbook of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPLP01 Domaine in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.05.001-A Page:
7/8
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
Obstruction memory:
8 megawords
Time CPU To use: 2.00 seconds
Handbook of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPLP01 Domaine in L with geometrical singularity
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.05.001-A Page:
8/8
7
Summary of the results
Two modelings give results whose certain values exceed the fixed tolerance
initially (1%):
· for modeling A (PLAN 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 grid
(more important refinement in the zone of the geometrical singularity).
The results are regarded as acceptable taking into account modelings carried out.
Handbook of Validation
V4.05 booklet: Stationary thermics of the plane structures
HT-66/02/001/A
Outline document