Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
1/14
Organisme(s) : EDF/MTI/MMN, DeltaCAD
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
Document : V3.01.101
SSLL101 - Tuyauterie : Problème de HOVGAARD
Résumé :
Il s'agit d'un test élastique linéaire, en statique, d'une tuyauterie tridimensionnelle non coplanaire comportant
des coudes. Il existe un test en dynamique de la même structure (SDLX02) [V2.05.002].
On teste les éléments POU_D_T, POU_C_T, TUYAU (SEG3 et SEG4) et TUYAU_6M (SEG3) par l'intermédiaire
de cinq modélisations, chaque coude est représenté par :
·
20 éléments POU_D_T dans la modélisation A,
·
2 éléments POU_C_T dans la modélisation B,
·
28 éléments TUYAU (SEG3) dans la modélisation C,
·
28 éléments TUYAU_6M (SEG3) dans la modélisation D,
·
28 éléments TUYAU (SEG4) dans la modélisation E.
Les chargements sont :
·
pesanteur,
·
dilatation thermique,
·
forces nodales.
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
2/14
1
Problème de référence
1.1 Géométrie
La géométrie, ainsi que les points de modélisation sont représentés sur la figure suivante :
0.
A
1.828
0.922
2.75
z
0.922
B
=
0.922
=
0.
4 5
=
3
6
=
7
=
A
8
9
10
B
3.69
11
2
12
y
13
=
14
=
=
15
=
1.96
=
x
1
longueurs données en mètres
·
diamètre extérieur du tuyau :
0.185 m
·
épaisseur du tuyau :
6.12 mm
·
rayon de courbure des coudes :
0.922 m
·
tuyauterie pleine d'eau
1.2
Propriétés de matériaux
E = 1.658 E+11 Pa
= 0.3
= 13404.10 kg/m3
= 0.1288 E4 /C°
1.3
Conditions aux limites et chargements
·
Points 1 et 15 encastrés,
·
Chargement :
1) Pesanteur selon z,
2) Elévation uniforme de température de 472.22 C°,
3) Forces
nodales.
Noeuds
2
3
4 - 10
5 - 9
6 - 7 - 8
11
12 - 13
14
Fz (N)
327.654
102.5145 222.687
176.580
624.897 788.724
214.839
117.720
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
3/14
2
Solution de référence
2.1
Méthode de calcul utilisée pour la solution de référence
Les solutions de référence retenues pour vérifier les modélisations Code_Aster sont les suivantes :
· pour les modélisations POUTRE : comparaison avec les codes : POUX, ADL et TITUS-T [bib1],
utilisant une modélisation de type poutre,
· pour les modélisations TUYAUX : comparaison avec le code ABAQUS, utilisant une
modélisation de type tuyau. Le nombre de mode de Fourier (M) utilisé lors du calcul de la
référence est identique à celui utilisé lors des calculs avec le Code_Aster.
2.2
Résultats de référence
Cas de
Déplacement
Modélisation Poutre
Modélisation
Modélisation
Chargement
au point 3
(POUX,ADL,TITUS)
Tuyau: M=3
Tuyau: M=6
(ABAQUS)
(ABAQUS)
Poids propre
DX 0.1658E3
0.16517E3
0.16512E3
DY 0.2040E4 0.13870E4
0.13946E4
DZ 0.8010E5 0.80376E5
0.80369E5
DX 0.1651E3
0.16445E3
0.16441E3
Force nodale
DY 0.2080E4 0.14245E4
0.14320E4
DZ 0.9516E5 0.10047E4
0.10047E4
DX 6.1418E3
6.3277E3
6.3236E3
Dilatation
DY 13.090E3
13.092E3 13.093E3
DZ 16.799E3
16.798E3 16.798E3
2.3
Incertitude sur la solution
2%
2.4 Références
bibliographiques
[1]
M.W. KELLOG Co. Design of Piping Systems. New York, 1956 - Problème n°5.9
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
4/14
3 Modélisation
A
3.1
Caractéristiques de la modélisation
Les éléments courbes sont modélisés par des éléments droits.
Un demi élément courbe est modélisé par 20 éléments droits.
2.75
=
=
N400
=
N500 N600
=
N300
N700
=
A
N800
N900
N1000
B
3.69
N1100
N200
N1200
y
N1300
=
N1400
=
=
N1500
=
1.96
=
x
N100
Longueurs données en mètres
3.2
Caractéristiques du maillage
Nombre de noeuds : 93
Nombre de mailles et type : 92 POU_D_T
3.3 Fonctionnalités
testées
Commandes
AFFE_CHAR_MECA PESANTEUR
TEMP_CALCULEE
AFFE_MODELE 'MECANIQUE' 'POU_D_T'
TOUT
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
5/14
4
Résultats de la modélisation A
4.1 Valeurs
testées
Identification Déplacement
Référence
poutre
Aster %
Poids propre
N300 DX
0.1658E3
0.1658E3
0.0
DY 0.2040E4
0.2039E4
0.02
DZ 0.8010E5
0.8010E5
0.0
N300 DX
0.1651E3
0.1651E3
0.04
Force nodale
DY 0.2080E4
0.2080E4
0.01
DZ 0.9516E5
0.9516E5
0.004
N300 DX
6.1418E3
6.1413E3
0.007
Dilatation
DY 13.090E3
13.091E3 0.012
DZ 16.799E3
16.799E3 0.003
4.2 Remarque
Les différences entre les résultats Aster et la solution de référence poutre sont toutes inférieures à
0.04%
4.3 Paramètres
d'exécution
Version : 5.6
Machine : Origin 2000
Encombrement mémoire :
16 Mo
Temps CPU User :
3 secondes
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
6/14
5 Modélisation
B
5.1
Caractéristiques de la modélisation
Les éléments courbes sont modélisés par des éléments "POU_C_T".
2.75
=
=
N400
=
N500 N600
=
N300
N700
=
A
N800
N900
N1000
B
3.69
N1100
N200
N1200
y
N1300
=
N1400
=
=
N1500
=
1.96
=
x
N100
Longueurs données en mètres
5.2
Caractéristiques du maillage
Nombre de noeuds : 15
Nombre de mailles et type :
10 POU_D_T
4
POU_C_T
5.3 Fonctionnalités
testées
Commandes
AFFE_CHAR_MECA PESANTEUR
TEMP_CALCULEE
AFFE_MODELE 'MECANIQUE'
'POU_D_T'
TOUT
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
7/14
6
Résultats de la modélisation B
6.1 Valeurs
testées
Identification Déplacement
Référence
poutre Aster %
Poids propre
N300 DX
0.1658E3
0.1658E4
0.017
DY 0.2040E4
0.2053E5
0.65
DZ 0.8010E5
0.8010E6
0.006
N300 DX
0.1651E3
0.1652E4
0.04
Force nodale
DY 0.2080E4
0.2080E5
0.02
DZ 0.9516E5
0.9516E6
0.002
N300 DX
6.1418E3
6.1404E3
0.02
Dilatation
DY 13.090E3
13.090E2 0.005
DZ 16.799E3
16.799E2
0.003
6.2 Remarque
Les différences entre les résultats Aster et les résultats de référence poutre sont toutes inférieures à
0.02% sauf pour DY en poids propre où l'écart est de 0.65%
6.3 Paramètres
d'exécution
Version : 5.6
Machine : Origin 2000
Encombrement mémoire :
16 Mo
Temps CPU User :
3 secondes
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
8/14
7 Modélisation
C
7.1
Caractéristiques de la modélisation
5
E
Modélisation TUYAU (SEG3)
D
4
z
A
F
B
y
5
10
G
4
Conditions aux limites : Points C et H
H
- DDL de Poutre : DX = DY = DZ = DRX = DRY = DRZ = 0
- DDL de Coque : UIm = VIm = WIm = 0 (m=2,3)
C
UOm = VOm = WOm = 0 (m=2,3)
x
WI1 = WO1 = WO = 0
7.2
Caractéristiques du maillage
Nombre de noeuds :
57
Nombre de mailles et type :
28 SEG3
7.3 Fonctionnalités
testées
Commandes
AFFE_MODELE AFFE
MODELISATION='TUYAU'
AFFE_CARA_ELEM
POUTRE:(SECTION:'CERCLE')
ORIENTATION:(CARA:'GENE_TUYAU'
VALE : (X Y Z)
AFFE_CHAR_MECA PESANTEUR
TEMP_CALCULEE
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
9/14
8
Résultats de la modélisation C
8.1 Valeurs
testées
Identification Déplacement
Référence
tuyau Aster %
(M = 3)
Poids propre
Point D DX
0.16517E3
0.1636E3
0.93
DY
0.13870E4
0.1251E4
9.80
DZ
0.80376E5
0.8018E5
0.24
Point D DX
0.16445E3
0.1629E3
0.94
Force nodale
DY
0.14245E4
0.1288E4
9.61
DZ
0.10047E4
0.1003E4
0.20
Point D DX
6.3277E3
6.4534E3
1.99
Dilatation
DY
13.092E3
13.103E3
0.08
DZ
16.798E3
16.880 E3
0.49
8.2 Remarque
Les résultats obtenus avec le Code_Aster sont similaires à ceux d'ABAQUS par des éléments tuyaux
sauf pour le déplacement DY (poids propre et force nodale) où l'écart est de l'ordre de 10%.
8.3 Paramètres
d'exécution
Version : 5.6
Machine : SGI-Origin2000 R12000
Encombrement mémoire : 16 Mo
Temps CPU User : 4,7 secondes
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
10/14
9 Modélisation
D
9.1
Caractéristiques de la modélisation
5
E
Modélisation TUYAU_6M (SEG3)
D
4
z
A
F
B
y
5
10
G
4
Conditions aux limites : Points C et H
H
- DDL de Poutre : DX = DY = DZ = DRX = DRY = DRZ = 0
- DDL de Coque : UIm = VIm = WIm = 0 (m=2,6)
C
UOm = VOm = WOm = 0 (m=2,6)
x
WI1 = WO1 = WO = 0
Longueurs données en mètres
9.2
Caractéristiques du maillage
Nombre de noeuds :
57
Nombre de mailles et type :
28 SEG3
9.3 Fonctionnalités
testées
Commandes
AFFE_MODELE AFFE
MODELISATION='TUYAU_6M'
AFFE_CARA_ELEM
POUTRE:(SECTION:'CERCLE')
ORIENTATION:(CARA:'GENE_TUYAU'
VALE : (X Y Z)
AFFE_CHAR_MECA PESANTEUR
TEMP_CALCULEE
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
11/14
10 Résultats de la modélisation D
10.1 Valeurs
testées
Identification Déplacement
Référence
tuyau Aster %
(M = 6)
Poids propre
Point D DX
0.16512E3
0.1636E3
0.93
DY
0.13946E4
0.1258E4
9.78
DZ
0.80369E5
0.8018E5
0.24
Point D DX
0.16441E3
0.1629E3
0.94
Force nodale
DY
0.14320E4
0.1295E4
9.58
DZ
0.10047E4
0.1003E4
0.21
Point D DX
6.3236E3
6.4495E3
1.99
Dilatation
DY
13.093E3
13.104E3
0.08
DZ
16.798E3
16.880E3
0.49
10.2 Remarque
Les résultats obtenus avec le Code_Aster sont similaires à ceux d'ABAQUS pour des éléments tuyaux
sauf pour le déplacement DY (poids propre et force nodale) où l'écart est de l'ordre de 10%.
10.3 Paramètres
d'exécution
Version : 5.6
Machine : SGI-Origin2000 R12000
Encombrement mémoire : 16 Mo
Temps CPU User : 10,56 secondes
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
12/14
11 Modélisation
E
11.1 Caractéristiques de la modélisation
5
E
Modélisation TUYAU (SEG4)
D
4
z
A
F
B
y
5
10
G
4
Conditions aux limites : Points C et H
H
- DDL de Poutre : DX = DY = DZ = DRX = DRY = DRZ = 0
- DDL de Coque : UIm = VIm = WIm = 0 (m=2,3)
C
UOm = VOm = WOm = 0 (m=2,3)
x
WI1 = WO1 = WO = 0
11.2 Caractéristiques du maillage
Nombre de noeuds :
85
Nombre de mailles et type :
28 SEG4
11.3 Fonctionnalités
testées
Commandes
CREA_MAILLAGE MODI_MAILLE
OPTION:'SEG3_4'
AFFE_MODELE AFFE
MODELISATION='TUYAU'
AFFE_CARA_ELEM
POUTRE:(SECTION:'CERCLE')
ORIENTATION:(CARA:'GENE_TUYAU'
VALE : (X Y Z)
AFFE_CHAR_MECA PESANTEUR
TEMP_CALCULEE
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
13/14
12 Résultats de la modélisation E
12.1 Valeurs
testées
Identification Déplacement
Référence
tuyau Aster %
(M = 3)
Poids propre
Point D DX
0.16517E3
0.1648E3
0.22
DY
0.13870E4
0.1321E4
4.73
DZ
0.80376E5
0.8024E5
0.18
Point D DX
0.16445E3
0.1638E3
0.37
Force nodale
DY
0.14245E4
0.1400E4
1.74
DZ
0.10047E4
0.9997E5
0.50
Point D DX
6.3277E3
6.329E3
0.02
Dilatation
DY
13.092E3
13.105E3
0.10
DZ
16.798E3
16.843E3
0.27
12.2 Remarques
Le maillage en SEG4 est obtenu à partir d'un maillage SEG3 avec la commande CREA_MAILLAGE,
MODI_MAILLE avec l'option 'SEG3_4'. Il est important que le noeud milieu du SEG3 soit bien au
milieu, le Code_Aster vérifie cette condition avec une tolérance.
Les résultats obtenus avec le Code_Aster sont similaires à ceux d'ABAQUS avec des éléments tuyaux
sauf pour le déplacement DY (poids propre et force nodale) ou l'écart est de l'ordre de 5% et 2%.
12.3 Paramètres
d'exécution
Version : 5.6
Machine : SGI-Origin2000 R12000
Encombrement mémoire : 16 Mo
Temps CPU User : 7,62 secondes
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Code_Aster ®
Version
5.0
Titre :
SSLL101 - Tuyauterie : problème de HOVGAARD
Date :
29/10/01
Auteur(s) :
J.M. PROIX, F. LEBOUVIER Clé
:
V3.01.101-C Page :
14/14
13 Synthèse des résultats
Modélisation poutre :
Les résultats sont similaires à la solution de référence (modélisation poutre : moyenne de résultats de
3 codes) aussi bien pour la modélisation A, où chaque coude est discrétisé par 20 éléments droits,
POU_D_T, que pour la modélisation B, où l'on utilise les éléments courbe POU_C_T. On note
simplement, dans ce cas, un écart un peu plus important en une valeur de déplacement (0.65%).
Modélisation tuyau :
Les résultats Code_Aster sont similaires à ceux d'ABAQUS (pour des éléments tuyaux), sauf pour le
déplacement DY et pour les chargements poids propre et forces nodales où l'écart avec la solution de
référence est plus important avec les mailles SEG3 (10%) qu'avec les mailles SEG4 (5%).
Le chargement de dilatation donne des résultats similaires.
Ce cas-test permet de tester une tuyauterie non coplanaire.
Manuel de Validation
Fascicule V3.01 : Statique linéaire des structures linéiques
HI-75/01/010/A
Document Outline