Code_Aster
®
Version
7.2
Titre :
FDLV110 - Calcul de masse ajoutée sur modes obtenus
Date :
01/03/04
Auteur(s) :
N. GREFFET, G. ROUSSEAU
Clé
:
V8.01.110-A
Page :
1/8
Manuel de Validation
Fascicule V8.01 : Fluide
HT-66/04/005/A
Organisme(s) :
EDF-R&D/AMA, EDF-DPN/UTO
Manuel de Validation
Fascicule V8.01 : Fluide
Document : V8.01.110
FDLV110 - Calcul de masse ajoutée sur modes
obtenus par sous-structuration
Résumé :
Ce test du domaine de l'analyse modale et de l'interaction fluide-structure met en oeuvre le calcul de masse
ajoutée sur une structure constituée de trois cylindres concentriques séparés par deux anneaux de fluide (eau)
dont on suppose le comportement régi par la théorie potentielle (fluide parfait, incompressible au repos). Le
modèle est tridimensionnel pour l'eau. La structure est représentée par des éléments de type coque mince
dans la modélisation A (la structure est rigide dans la solution de référence). Celle-ci est caractérisée par deux
modes propres évalués par sous-structuration dynamique, avec interface de type CRAIG-BAMPTON.
L'intérêt du test réside dans l'utilisation de la fonctionnalité `
NOEUD_DOUBLE
` de l'opérateur `
CALC_MATR_AJOU
`.
Cette fonctionnalité permet de calculer les effets de masse ajoutée d'une structure représentée par un maillage
surfacique (sans épaisseur) qui se trouve baignée dans un fluide. Les fluides choisis dans ce cas-test sont de
densités différentes de part et d'autre du cylindre intermédiaire (eau à températures différentes).
Code_Aster
®
Version
7.2
Titre :
FDLV110 - Calcul de masse ajoutée sur modes obtenus
Date :
01/03/04
Auteur(s) :
N. GREFFET, G. ROUSSEAU
Clé
:
V8.01.110-A
Page :
2/8
Manuel de Validation
Fascicule V8.01 : Fluide
HT-66/04/005/A
1
Problème de référence
1.1 Géométrie
z
y
x
fluide
L
surface fluide
cylindre milieu
ressorts
cylindre externe
R
f
R
i
R
e
k
1
k
2
L=50m ; R
f
=1m ; R
i
=
5
3
m ; R
e
=3m ; k
1
=10
9
N.m
1
; k
2
=0.5 10
7
N.m
1
;
fluide
=1000 kg.m
3
;
s
=7800 kg.m
3
; épaisseur de la coque : 50 cm.
1.2
Propriétés des matériaux
Fluide : masse volumique
1
= 1000 kg.m
3
;
2
= 750 kg.m
3
.
Structure :
s
= 7800 kg/m
3
;
E
= 2.1 10
11
Pa ;
= 0.3 (acier).
1.3
Conditions aux limites et chargements
Le cylindre extérieur est d'une part relié à un bâti fixe via les quatre ressorts de raideur unitaire k
1
,
relié d'autre part au cylindre milieu par quatre ressorts de raideur unitaire k
2
. Les deux structures sont
rigides dans cette solution de référence.
Code_Aster
®
Version
7.2
Titre :
FDLV110 - Calcul de masse ajoutée sur modes obtenus
Date :
01/03/04
Auteur(s) :
N. GREFFET, G. ROUSSEAU
Clé
:
V8.01.110-A
Page :
3/8
Manuel de Validation
Fascicule V8.01 : Fluide
HT-66/04/005/A
2
Solution de référence
On calcule les modes propres du système après avoir vérifié ceux de chaque sous-structure. On
évalue ensuite la masse ajoutée sur les modes en air.
2.1
Décomposition en sous-structures
Première sous-structure : cylindre intermédiaire
La première sous-structure est constituée du cylindre intermédiaire et de quatre ressorts de raideur
k
2
= 10
7
N.m
1
. Ces ressorts sont encastrés à l'interface avec le cylindre externe qui constitue la
deuxième sous-structure (interface de type CRAIG-BAMPTON).
R
i
k
2
Masse du cylindre 1 : m
1
= 2.041 10
6
kg
Le cylindre étant rigide, son mouvement peut être modélisé par un système masse-ressort à un degré
de liberté :
m
1
Achacune de ses extrémités, le cylindre est relié à deux ressorts en parallèle : la raideur équivalente
de chacun est k
'
= 2k
2
La fréquence propre vaut alors :
1
2
1
4
2
1
2
2
1
m
k
m
k
f
=
=
, soit :
=
f
Hz
0 705
.
Deuxième sous-structure : cylindre extérieur
La deuxième sous-structure est le cylindre externe relié d'une part à l'interface par les mêmes
ressorts, d'autre part à un bâti fixe :
R
e
k
1
k
Masse du cylindre 2 : m
2
=3,674 10
6
kg
Code_Aster
®
Version
7.2
Titre :
FDLV110 - Calcul de masse ajoutée sur modes obtenus
Date :
01/03/04
Auteur(s) :
N. GREFFET, G. ROUSSEAU
Clé
:
V8.01.110-A
Page :
4/8
Manuel de Validation
Fascicule V8.01 : Fluide
HT-66/04/005/A
La raideur équivalente d'une attache de ce cylindre par le système de ressorts en série k
1
et k
2
valant
9,9 10
6
N.m
1
(quatre attaches du même type relient en parallèle le cylindre à un encastrement), la
fréquence propre est donnée par :
(
)
1
1
2
1
2
4
2
1
m
k
k
k
k
f
+
=
, soit :
Hz
.
f
522
0
=
N.B. : le troisième cylindre (cylindre intérieur) n'a pas été modélisé dans notre cas car il s'agint d'un
cylindre fixe. Il constitue donc une paroi fixe du domaine fluide.
Modes en air de la structure complète (cylindre intermédiaire et cylindre externe)
C'est un système à deux degrés de liberté :
m
2
m
1
k
3
=4k
1
k
4
=4k
2
Les fréquences propres de ce système sont données par la formule exacte [bib2] :
f
k
m
k
m
k
m
k
m
k
m
k
m
k k
m m
i
=
+
+
±
+
+
-
1
2
4
32
3
2
4
2
4
1
3
2
4
2
4
1
2
3 4
1
2
,
soit
Hz
et
Hz
263
5
497
0
2
1
.
f
.
f
=
=
.
Les deux modes propres admettent, pour valeur numérique :
X
X
X
X
1
2 1
3
1
2 2
3
1
5 10
9 10
1
=
= -
-
-
.
.
et
.
2.2
Calcul de la matrice de masse ajoutée
Potentiels fluides
Reprenant
[bib1], on établit que :
( )
( )
( )
(
)
-
+
-
+
+
-
+
×
-
=
-
+
×
-
+
+
-
+
=
-
-
2
2
2
2
2
2
2
2
2
2
2
2
2
2
3
2
1
2
2
2
2
2
3
2
2
2
2
2
2
2
2
2
1
1
;
.
9
.
5
;
i
e
e
i
e
i
e
i
e
f
i
f
i
i
i
e
e
i
e
i
e
i
e
f
i
f
i
i
R
R
R
R
R
R
R
R
R
R
R
R
R
R
e
R
R
R
R
R
e
R
R
R
R
R
R
R
R
R
et
Code_Aster
®
Version
7.2
Titre :
FDLV110 - Calcul de masse ajoutée sur modes obtenus
Date :
01/03/04
Auteur(s) :
N. GREFFET, G. ROUSSEAU
Clé
:
V8.01.110-A
Page :
5/8
Manuel de Validation
Fascicule V8.01 : Fluide
HT-66/04/005/A
La forme de la matrice de masse ajoutée, dans cette configuration, est :
=
22
21
12
11
a
a
a
a
a
M
M
M
M
M
Avec :
.
,
,
kg
kg
kg
15318
)
10
.
9
(
)
10
.
5
(
10
.
676
,
2
10
.
753
,
1
2
2
2
2
2
3
2
2
2
2
2
2
2
2
2
3
12
6
2
2
2
2
2
22
6
2
2
2
2
2
2
2
2
2
11
-
=
-
+
×
-
-
+
+
-
+
×
=
=
-
+
=
=
-
+
+
-
+
=
-
-
i
e
e
i
e
i
e
e
i
f
i
f
i
i
a
i
e
e
i
e
a
i
e
e
i
f
i
f
i
i
a
R
R
R
R
R
R
R
R
R
R
R
R
R
R
L
M
R
R
R
R
R
L
M
R
R
R
R
R
R
R
R
R
L
M
Le coefficient de couplage inertiel
12
a
M
est considéré comme négligeable devant les coefficients
d'auto-masse ajoutée
11
a
M
et
22
a
M
. Les fréquences propres du système ne dépendent, en première
approximation, que de ces deux dernièrs coefficients.
2.3
Résultats de référence
Résultat analytique.
2.4 Références
bibliographique
[1]
ROUSSEAU G., LUU H.T. : Masse, amortissement et raideur ajoutés pour une structure
vibrante placée dans un écoulement potentiel - Bibliographie et implantation dans le
Code_Aster - HP-61/95/064
[2]
BLEVINS R.D. : Formulas for Natural frequency and mode shape, Ed. Krieger
Code_Aster
®
Version
7.2
Titre :
FDLV110 - Calcul de masse ajoutée sur modes obtenus
Date :
01/03/04
Auteur(s) :
N. GREFFET, G. ROUSSEAU
Clé
:
V8.01.110-A
Page :
6/8
Manuel de Validation
Fascicule V8.01 : Fluide
HT-66/04/005/A
3 Modélisation
A
3.1
Caractéristiques de la modélisation
Pour le système 3D sur lequel on calcule les coefficients ajoutés :
Cylindre :
2400 mailles QUAD4
éléments de coques MEDKQU4
12 mailles SEG2
éléments ressorts MECA_DIS_T_L
Fluide :
3600 mailles QUAD4
éléments thermiques THER_FACE4
sur les surfaces cylindriques
7200 mailles HEXA8
éléments thermiques THER_HEXA8
dans le volume annulaire fluide
3.2 Fonctionnalités
testées
Commandes
CALC_MATR_AJOU OPTION
`NOEUD_DOUBLE'
`MASS_AJOU'
Code_Aster
®
Version
7.2
Titre :
FDLV110 - Calcul de masse ajoutée sur modes obtenus
Date :
01/03/04
Auteur(s) :
N. GREFFET, G. ROUSSEAU
Clé
:
V8.01.110-A
Page :
7/8
Manuel de Validation
Fascicule V8.01 : Fluide
HT-66/04/005/A
4
Résultats de la modélisation A
4.1 Valeurs
testées
Fréquences analytique
en air (Hz)
Fréquences calculées
en air par Aster (Hz)
Ecart (%)
Premier mode en air
0.497 0.496
0.2
Deuxième mode en air
5.263 5.147
0.2
Masse ajoutée théorique
(kg)
Masse ajoutée calculée
(kg)
Ecart (%)
M
11
1.753 10
6
1.741
10
6
0.6
M
22
2.675 10
6
2.567
10
6
4.0
Fréquences analytique
des modes en eau (Hz)
Fréquences en eau
calculées par Aster (Hz)
Ecart (%)
Premier mode en eau
0.365
0.364
0.2
Deuxième mode en eau
4.004 4.061
1.0
Code_Aster
®
Version
7.2
Titre :
FDLV110 - Calcul de masse ajoutée sur modes obtenus
Date :
01/03/04
Auteur(s) :
N. GREFFET, G. ROUSSEAU
Clé
:
V8.01.110-A
Page :
8/8
Manuel de Validation
Fascicule V8.01 : Fluide
HT-66/04/005/A
5
Synthèse des résultats
Le calcul de masse ajoutée sur des modes estimés par sous-structuration est satisfaisant. Ceci a
permis de valider l'option `
NOEUD_DOUBLE
' de la commande `
CALC_MATR_AJOU
'. L'écart observé sur
le deuxième coefficient de masse ajoutée s'explique par la discrétisation du deuxième cylindre. Le
nombre d'éléments est un peu insuffisant pour calculer de manière exacte l'intégrale du champ de
pression sur la structure.