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 œuvre 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).

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.

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}_{\mathrm{2 }}={10}^{7}{\mathrm{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).

Masse du cylindre 1: \({m}_{1}=2.041{10}^{6}\mathrm{kg}\)

Le cylindre étant rigide, son mouvement peut être modélisé par un système masse-ressort à un degré de liberté:

Achacune de ses extrémités, le cylindre est relié à deux ressorts en parallèle: la raideur équivalente de chacun est \(k’=2{k}_{2}\)

La fréquence propre vautalors :

, soit:

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:

Masse du cylindre 2: \({m}_{2}=3,674{10}^{6}\mathrm{kg}\)

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}{\mathrm{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:

, soit:

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é:

Les fréquences propres de ce système sont données par la formule exacte [bib2]:

,

soit

\({f}_{1}=0.497\mathrm{Hz}\) et \({f}_{2}=5.263\mathrm{Hz}\) .

Les deux modes propres admettent, pour valeur numérique:

.

Calcul de la matrice de masse ajoutée

Potentiels fluides

Reprenant [bib1], on établit que:

La forme de la matrice de masse ajoutée, dans cette configuration, est:

Avec:

Le coefficient de couplage inertiel

est considéré comme négligeable devant les coefficients d’auto-masse ajoutée

et

. Les fréquences propres du système ne dépendent, en première approximation, que de ces deux dernièrs coefficients.

Résultats de référence

Résultat analytique.

Références bibliographique

• 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

• BLEVINS R.D. : Formulas for Natural frequency and mode shape, Ed. Krieger

Modélisation A

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

Valeurs testées

	Fréquences analytique en air ( \(\mathrm{Hz}\) )
Premier mode en air	0.497
Deuxième mode en air	5.263
	Masse ajoutée théorique ( \(\mathrm{kg}\) )
\({M}^{11}\)	1.753 106
\({M}^{22}\)	2.675 106
	Fréquences analytique des modes en eau ( \(\mathrm{Hz}\) )
Premier mode en eau	0.365
Deuxième mode en eau	4.004

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.