Analyse modale d’une structure industrielle avec prise

en compte du couplage fluide/structure

Jean-François Sigrist, Christian Lainé, Bernard Peseux

To cite this version:

Jean-François Sigrist, Christian Lainé, Bernard Peseux. Analyse modale d’une structure industrielle
avec prise en compte du couplage fluide/structure. Mécanique and Industries / Mécanique et Indus-
tries, EDP Sciences, 2005, 6 (5), pp.553 - 563. <10.1051/meca:2005067>. <hal-01377405>

HAL Id: hal-01377405

[Link]
Submitted on 7 Oct 2016

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est

archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents
entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non,
lished or not. The documents may come from émanant des établissements d’enseignement et de
teaching and research institutions in France or recherche français ou étrangers, des laboratoires
abroad, or from public or private research centers. publics ou privés.

Distributed under a Creative Commons ttribution�ttribution 4.0 International License

 1

Analyse modale d’une structure industrielle avec prise en compte

du couplage fluide/structure
Jean-François Sigrist1,a , Christian Lainé1 et Bernard Peseux2
1
DCN Propulsion, Service Scientifique et Technique, 44620 La Montagne, France
2
École Centrale de Nantes, Laboratoire Mécanique Matériaux, 1 rue de La Noë, 44321 Nantes Cedex, France

Résumé – Nous présentons une analyse modale d’une structure industrielle couplée avec un fluide, en utili-
sant les techniques numériques de calculs couplés fluide/structure. Compte tenu de la nature axisymétrique
de la géométrie et de la nature non axisymétrique des équations de couplage, la modélisation du problème
est réalisée au moyen d’éléments finis axisymétriques développés en série de Fourier. Un code de calcul
est implanté dans Matlab pour permettre l’analyse modale de la structure. Différentes formulations du
problème sont comparées ; les résultats de calcul Matlab sont comparés avec les résultats de calculs ob-
tenus avec le code généraliste Ansys. Les développements mis en œuvre pour cet exemple seront à terme
intégrés dans le code Ansys pour l’étude de problèmes couplés en pression/déplacement avec développement
en série de Fourier.

Mots clés : Analyse modale / interactions fluide/structure / éléments finis axisymétriques en série de
Fourier / formulations symétrique et non symétrique

Abstract – Modal analysis of an industrial structure with fluid-structure interaction. This

paper gives an example of numerical methods with finite element coupling for the numerical study of an
industrial fluid/structure coupled problem. The problem is solved by coupling a finite element discretization
of both fluid and structure domain. The coupled system is described in terms of pressure and displacement
for the fluid and structure problem. The unknown degrees of freedom are expanded in a Fourier serie. A
numerical code is developed in Matlab in order to perform the modal analysis of the coupled problem.
Numerical calculations are performed and compared with the Ansys code in order to validate our numerical
developments.

Key words: Modal Analysis / fluid-structure inetraction / harmonic axisymmetric finite elements /
symmetric and non-symmetric formulations

Introduction Le présent article se propose d’illustrer une applica-

tion de ces techniques de calcul à l’étude dynamique d’une
La simulation numérique de problèmes couplés structure représentative de l’ensemble panier et cuve d’un
fluide/structure a fait l’objet de nombreuses études réacteur de propulsion navale, avec prise en compte du
numériques qui ont permis le développement de méthodes couplage fluide/structure.
de calcul, consistant à réaliser un couplage entre une Le modèle simplifié de la structure étudiée est
discrétisation de type éléments finis du domaine struc- représenté par la figure 1. On adopte une modélisation
ture et une discrétisation de type équation intégrale ou axisymétrique de l’ensemble et on distingue quatre prin-
éléments finis du problème fluide [1, 2]. Ces techniques cipaux éléments :
de calcul sont développées dans des laboratoires de re-
cherche, et appliquées par exemple pour réaliser l’étude – le corps du panier (1) : il s’agit de la partie
de structures immergées, avec prise en compte ou non de supérieure de la structure intérieure du réacteur,
phénomènes de surface libre dans le fluide [3–10], mais qui contient les éléments combustibles. La géométrie
leur utilisation courante en bureau d’étude est cependant est essentiellement cylindrique. Les caractéristiques
assez rare dans le cadre de projets industriels [11, 12]. géométriques et physiques sont notées E1 (module
d’Young), ρ1 (masse volumique), ν1 (coefficient de
a
Auteur correspondant : [Link]@[Link] Poisson) et e1 (épaisseur) ;
 2

Fig. 1. Géométrie industrielle étudiée. Représentation axisymétrique de l’ensemble panier et cuve d’un réacteur de propulsion
navale.

– le fond du panier (2) : il s’agit de la partie la structure, les inconnues du problème couplé en série de
inférieure de la structure intérieure du réacteur. La Fourier sont développées selon la relation [13, 14] :
géométrie est plus complexe et fait apparaı̂tre deux
parties sphériques. Les grandeurs caractéristiques sont    0  n 
 ur   ur   ur cos (nθ) 
notées E2 , ρ2 , ν2 et e2 ; 
 
 
 
 
 


     
– la cuve (3) : il s’agit de l’enveloppe externe du  uθ 
 
 u0θ 
 n=+∞


 u n
θ sin (nθ) 

= +
réacteur. Elle est modélisée en première analyse d’une 
 uz 
 
 u0z 
 
 unz cos (nθ) 

partie cylindrique et d’une partie sphérique. Les gran- 
 
 
 
 n=1 
 


p   
 p0   pn cos (nθ) 
 
deurs caractéristiques sont notées E3 , ρ3 , ν3 et e3 ;
– le fluide primaire (4) : il est contenu entre le panier et  m 
la cuve. On note ρF sa masse volumique et c la célérité  ur sin (mθ) 

 

des ondes sonores dans ce milieu. 
 m 

m=+∞
  −u θ cos (mθ) 
+ (1)
Parmi la grande variété de modélisations géométriques 
 um
z sin (mθ) 

et de représentations des inconnues du problème couplé, m=1 
 


 p sin (mθ)
m 

nous privilégions une approche compatible avec les exi-
gences d’analyse en bureau d’étude : simplicité du modèle
géométrique, interprétation aisée des modes propres. Ce développement fait apparaı̂tre des termes axi-
L’analyse est basée sur une discrétisation de type symétriques (ordre 0), des termes symétriques (ordre n ≥
éléments finis pour la structure, couplée avec une 1) et des termes antisymétriques (ordre m ≥ 1). La sollici-
discrétisation de type éléments finis pour le fluide. La tation dynamique appliquée correspond à un chargement
modélisation adopte une formulation : dans une direction transverse à la structure de révolution,
– en déplacement pour la structure, les inconnues du nous nous s’intéressons donc à la détermination des modes
problème structure sont les déplacements, notés ur , symétriques d’ordre 1 (i.e. en cos θ, n = 1).
uθ et uz en coordonnées cylindriques (approche La discrétisation du problème couplé est réalisée
lagrangienne) ; dans un code développé dans Matlab [15], utilisant un
– en pression pour le fluide, l’inconnue du problème élément coque axisymétrique1 et un élément acoustique
fluide est le champ de pression p (approche 1
eulérienne). Compte tenu des épaisseurs de la structure industrielle, la
modélisation par éléments coque atteint ses limites. L’objet
Compte tenu de l’hypothèse de modélisation axi- du présent travail étant de mettre en évidence le couplage des
symétrique de la géométrie et de la nature non axi- modes de structure par le fluide, l’utilisation d’un modèle de
symétrique de la sollicitation dynamique appliquée sur coque élastique reste cependant justifiée dans ce contexte.
 3

axisymétrique. Les calculs seront comparés avec des

calculs éléments finis réalisés avec le code commercial
Ansys2 [16], utilisant deux types de modèles :
– modèle bidimensionnel axisymétrique mettant en
œuvre des éléments finis axisymétriques développés en
série de Fourier et formulés en déplacement ;
– modèle tridimensionnel mettant en œuvre des
éléments finis volumiques formulés en pression/
déplacement.
Le choix d’une représentation en pression/
déplacement du problème couplé, avec développement
des inconnues en série de Fourier, présente les avantages
suivants et répond aux exigences d’une analyse pour
application en bureau d’étude :
Fig. 2. Élément fini coque.
– représentation simple du problème fluide par la va-
riable pression ; Tableau 1. Premières fréquences propres du panier dans le
– non singularité des matrices de masse et de raideur du vide.
problème fluide pour les harmoniques n ≥ 1 et m ≥ 1 ;
– absence de modes propres non physique pouvant pa- Fréquence (Hz) Calcul Ansys Calcul Matlab ε (%)
f1 44,69 44,70 0,02
rasiter l’analyse modale.
f2 140,67 140,77 0,07
L’analyse modale conduite sur le problème industriel f3 245,71 245,72 0,00
est réalisée en considérant la structure élastique (panier et f4 284,25 284,62 0,13
cuve) seule, le fluide acoustique seul et enfin le problème f5 322,1 322,80 0,21
couplé élasto-acoustique. Dans ce dernier cas, différentes
formulations sont étudiées. Tableau 2. Premières fréquences propres de la cuve dans le
vide.
Fréquence (Hz) Calcul Ansys Calcul Matlab ε (%)
1 Étude de la structure dans le vide f1 78,71 78,708 –0,01
f2 260,65 260,92 0,10
f3 431,58 424,88 –1,55
Nous commençons l’étude par une analyse modale sur
f4 480,22 480,64 0,09
chaque élément de la structure (panier et cuve) dans le
f5 555,87 553,28 0,47
vide. La discrétisation du problème utilise des éléments
coque à deux nœuds et quatre degrés de liberté par nœud
(voir Fig. 2). L’analyse modale de la structure consiste à résoudre
Le champ de déformations est donné par le modèle le problème aux valeurs propres symétrique suivant :
de Sanders-Koiters, qui permet de décrire le comporte-
ment en flexion de la coque élastique [18]. La relation KS − ω 2 MS × U(ω) = 0 (3)
déplacements/déformations est alors la suivante :
avec KS et MS les matrices du système structure. Le
 
 ∂uz  problème aux valeurs propres est symétrique, nous uti-

 − 


 ∂z 
 lisons l’algorithme de Cholesky [19, 20] pour extraire les

 


 1 ∂u u 

modes propres du système.

 −
θ
+
r 
 Les premières fréquences propres du panier et de la

 r ∂θ r 


 
 cuve, calculées avec le code développé dans Matlab, sont

 ∂u ∂u 


 −
1 z r 
 données dans les tableaux 1 et 2 et comparées avec

 + 

r ∂θ ∂z les fréquences propres calculées avec Ansys, en utilisant
{ε} = 2
∂ ur (2) un modèle de coque élastique axisymétrique en série de

 − 


 
 Fourier (éléments shell61 [16]). Les déformées modales

 ∂z 2 


 2 
 correspondantes dans le plan (r, z) sont données par la fi-

 1 ∂ ur ∂uθ 


 + 
 gure 3 ; on considère les modes symétriques d’harmonique

 r 2 ∂θ 2 ∂θ 


 
 n = 1, la sollicitation dynamique sur la structure étant en
2
 ∂ 2
u 3 ∂u 1 ∂u 


 −
r
−
θ z 
 cos (θ). Les résultats donnés par le code développé dans
r ∂z∂θ 4 ∂z
+
4r ∂θ 
2 Matlab sont identiques aux résultats de calcul Ansys, à
la fois en termes de valeurs des fréquences propres que de
2 déformées modales.
Un précédent travail a permis de valider, par comparaison
calcul/essai, l’utilisation du code Ansys pour l’analyse modale La figure 4 donne les courbes de convergence des
de structures élastiques couplées avec un fluide acoustique [17]. premières fréquences propres du panier et de la cuve pour
 4

Fig. 3. Premiers modes propres d’ordre n = 1 du panier (a) et de la cuve (b) dans le vide obtenus avec le code de calcul
développé dans Matlab.

Fig. 4. Courbes de convergence des premières fréquences

propres du panier et de la cuve dans le vide.

un nombre N d’éléments, soit un nombre total de degrés

de liberté de 4N (en ordonnée, on porte ε l’écart relatif Fig. 5. Élément fini fluide.
entre les résultats du calcul Matlab et le résultat du cal-
cul de référence avec Ansys avec les éléments shell61).
La convergence est obtenue assez rapidement pour les
modèles étudiés, ce qui permettra de travailler avec des
modèles de taille raisonnable dans le cas de l’analyse mo-
dale du système avec couplage fluide/structure.

2 Étude de la cavité acoustique seule

Nous réalisons ensuite une analyse modale de la ca-
vité fluide délimitée par le panier et la cuve, considérés
comme rigides. La discrétisation du domaine fluide utilise
des éléments finis linéaires à quatre nœuds et un degré de
liberté par nœud (voir Fig. 5).
Le maillage de l’ensemble est visualisé par la figure 6.
L’analyse modale de la cavité acoustique consiste
à résoudre le problème aux valeurs propres symétrique
suivant :
KF − ω 2 MF × P(ω) = 0 (4)
avec KF et MF les matrices du système fluide. Nous utili- Fig. 6. Maillage de la cavité acoustique délimitée par le panier
sons les mêmes méthodes d’extraction de valeurs propres et la cuve dans Matlab.
que pour le problème structure.
 5

Tableau 3. Premières fréquences propres de la cavité acous- Tableau 4. Fréquences propres du système panier et cuve.
tique délimitée par le panier et la cuve. Comparaison entre un calcul sans eau, avec masse d’eau re-
portée et avec masse d’eau ajoutée (calcul couplé avec le code
Fréquence (Hz) Calcul Ansys Calcul Matlab ε (%)
Ansys par couplage éléments finis/éléments finis en formula-
f1 224,58 224,57 0,00
tion déplacement/déplacement avec représentation des incon-
f2 307,96 307,04 –0,30
nues en série de Fourier).
f3 419,06 406,12 –3,09
f4 479,80 482,52 0,57 Fréquence (Hz) Sans eau Masse d’eau Masse d’eau
f5 669,69 669,87 0,93 reportée ajoutée
f1 44,70 43,38 26,03
f2 78,71 75,55 65,52
Les premières fréquences propres de la cavité (pour f3 140,77 136,56 102,94
l’harmonique n = 1), calculées avec le code développé f4 245,71 238,27 162,30
dans Matlab, sont données dans le tableau 3 et com- f5 260,65 250,19 180,32
parées avec les fréquences propres calculées avec Ansys,
utilisant les éléments fluides acoustiques formulés en
Ceci met en évidence l’intérêt de mettre en œuvre un
déplacement3 et développés en série de Fourier (éléments
calcul couplé fluide/structure pour l’analyse modale de la
fluid81 du code [16]). Les formes modales correspon-
structure industrielle.
dantes (représentées en iso-valeurs de pression) dans le
Les résultats des calculs couplés sont obte-
plan (r, z) sont données par la figure 7. Les résultats
nus avec le code Ansys, dans une formulation
donnés par le code développé dans Matlab sont identiques
déplacement/déplacement. Cette formulation possède
aux résultats de calcul Ansys, à la fois en termes de va-
l’inconvénient de générer des modes non physiques,
leurs des fréquences propres que de déformées modales.
qui peuvent perturber l’analyse modale : ces modes
Le problème fluide comporte environ N = 350 degrés de
correspondent aux mouvements de cisaillement du fluide
liberté.
sans changement de volume, la structure restant fixe. On
lui préfère, pour une application en bureau d’étude, une
formulation en pression/déplacement (non symétrique)
3 Étude de la structure couplée avec le fluide ou en pression/potentiel des déplacements/déplacement
(symétrique), qui ne présente pas ces inconvénients [21].
3.1 Approche simplifiée Cependant, le code Ansys ne permet pas de traiter, pour
ce type de formulation, des problèmes axisymétriques
avec un développement des inconnues en série de
L’approche simplifiée du problème consiste à utili- Fourier [16]. Nous développons donc un code dans le
ser un modèle de masse d’eau reportée pour prendre en progiciel Matlab permettant de réaliser l’analyse modale
compte la présence du fluide : la masse d’eau présente de la structure étudiée, en utilisant les éléments finis
entre le panier et la cuve est reportée pour moitié entre présentés aux sections précédentes.
les deux structures, ce qui permet de respecter le bilan
masse total de l’ensemble structure et fluide. Cet artifice
de calcul ne permet pas de bien décrire les effets du fluide 3.2 Approche (( couplage fort )), formulation
sur la structure : l’utilisation du modèle de masse d’eau non symétrique
reportée conduit à sous-estimer très largement dans ce cas
l’influence du fluide sur la structure, comme le montrent La technique de couplage éléments finis/éléments finis
les résultats du tableau 4. Compte tenu du confinement, consiste à discrétiser les termes d’échanges entre le fluide
la masse d’eau ajoutée (qui traduit l’effet physique du et la structure ; on calcule ainsi une matrice d’interac-
fluide sur la structure) est nettement plus importante que tion fluide/structure en couplant les degrés de liberté en
la masse d’eau reportée : la baisse de fréquence propre est déplacement de la structure et le degré de liberté pression
moins importante avec le modèle de masse d’eau reportée du fluide (Fig. 8).
qu’avec le modèle de masse d’eau ajoutée [22]. D’autre La matrice de couplage est calculée par assemblage des
part, les effets de couplage des modes de cuve et de panier matrices élémentaires de couplage ri , pour un couplage
par le fluide ne sont pas pris en compte dans un modèle défini par la figure 8. Chaque matrice élémentaire ri est
de masse d’eau reportée. calculée à partir des fonctions de forme des deux éléments,
selon la relation :
3
Dans Ansys, les éléments fluides à formulation en zi+1
déplacement sont obtenus par modification des éléments so-
lides. La différence entre les éléments fluides et solides est prin- ri = {N F (z)} ni  [NS (z)] dz (5)
cipalement dans la loi de comportement : la raideur en cisaille- zi
ment est proche de zéro et les termes associés aux contraintes
sont ramenés à la valeur du module du fluide. On obtient ainsi avec {N F } le vecteur des fonctions de forme de l’élément
un élément qui ne peut reprendre le cisaillement, mais dont les fluide, [NS ] la matrice des fonctions de forme de l’élément
composantes de contrainte normale sont proches de la pression structure et ni  le vecteur des composantes de la normale
instantanée [16, 21]. locale sortante du milieu fluide.
 6

Fig. 7. Premiers modes propres de la cavité acoustique délimitée par le panier et la cuve.

Tableau 5. Validation de l’implantation de l’algorithme non

symétrique de Lanczos dans le code de calcul développé
dans Matlab. Comparaison sur l’exemple théorique proposé
dans [24].
Valeurs propres théoriques Valeurs propres calculées
0,10000 + 0,00000 i 0,1000 + 0,0000 i
4,00000 ± 1,00000 i 4,0000 ± 1,0000 i
2,40000 ± 7,32000 i 2,39974 ± 7,31897 i
7,00000 + 0,00000 i 6,99999 + 0,00000 i
7,00000 + 0,00000 i 7,00088 + 0,00000 i
8,00000 ± 127,00000 i 7,99997 ± 127,00001 i
11,00000 + 0,00000 i 11,00005 + 0,00000 i
130,000000 + 0,00000 i 129,99998 + 0,00000 i
Fig. 8. Couplage éléments finis/éléments finis. −2000,00000 + 0,00000 i −1999,99995 + 0,00000 i

Le problème aux valeurs propres à résoudre dans le cas L’algorithme de Lanczos non symétrique est implanté
du couplage de la géométrie industrielle prend la forme dans le code de calcul développé dans Matlab, et validé
générale suivante : par comparaison avec le cas théorique proposé dans [24]
et rappelé en annexe ; les résultats des deux calculs sont
 
 
KC 0 R T  UC (ω)  donnés dans le tableau 5.
C  
    On conduit une analyse modale pour un problème cou-
 0 KP R T
P  × UP (ω) =
    plant le fluide avec l’ensemble panier et cuve.
0 0 KF 
 P(ω) 
 Le tableau 6 compare les valeurs des premières
    valeurs propres du système complet (panier et cuve)
MC 0 0  UC (ω)  couplé avec le fluide dans le cas d’un calcul Matlab
 
    et d’un calcul Ansys. Le calcul de référence Ansys est
2
ω  0 M P 0  × UP (ω) (6)
   obtenu par couplage éléments finis/éléments finis en
−ρF RC −ρF RP MF  P(ω) 
  formulation déplacement/déplacement et représentation
des inconnues en série de Fourier, le calcul Matlab uti-
avec MC , KC les matrices masse et raideur de la cuve, lise le couplage éléments finis/éléments finis en formu-
MP , KP les matrices masse et raideur du panier, MF , lation pression/déplacement et représentation des incon-
KF les matrices du fluide et RC , RP les matrices de cou- nues en série de Fourier, selon les principes exposés
plage de la cuve et du panier avec le fluide. Les matrices précédemment.
masse et raideur du problème couplé sont creuses comme La figure 10 donne une représentation des déformées
le montre la figure 9 (la représentation proposée met en modales des premiers modes d’ordre n = 1 du système
évidence les éléments non nuls des matrices de masse et couplé panier, cuve et fluide primaire.
de raideur du système couplé). Les résultats de calcul obtenus avec les deux méthodes
Le système obtenu après couplage est non symétrique : sont proches, les écarts étant toutefois plus importants
nous utilisons alors une adaptation de l’algorithme de que pour le modèle de panier ou de cuve seuls (écarts de
Lanczos [23, 24] pour traiter le problème de l’extraction l’ordre de 4 % en moyenne). L’écart devient assez signi-
des éléments propres du problème non symétrique (7). ficatif pour les fréquences au-delà du cinquième mode ;
 7

Fig. 9. Structure des matrices masse et raideur pour le problème couplé panier, cuve et fluide en formulation (u, p).

Fig. 10. Premiers modes propres de l’ensemble panier et cuve couplé avec le fluide.

Tableau 6. Première fréquences propres du système couplé. en fonction de la taille N du problème matriciel (ε est
Comparaison des résultats de calcul Matlab et Anys. l’écart relatif entre les résultats du calcul Matlab et les
Fréquence (Hz) Calcul Ansys Calcul Matlab ε (%) résultats du calcul de référence avec Ansys en formula-
f1 26,66 26,03 –2,36 tion déplacement/déplacement) et la figure 12 compare
f2 64,02 65,52 2,50 les temps de calcul pour l’extraction des valeurs propres
f3 108,33 102,94 –4,98 du système couplé en utilisant un PC standard et une sta-
f4 158,24 162,30 2,57 tion Unix avec les mêmes versions de Matlab. L’utilisation
f5 201,23 180,32 10,39 de Matlab dans l’environnement Unix est nécessaire pour
pouvoir traiter un problème de grande dimension (à partir
de N ≈ 500).
nous atteignons les limites de l’outil Matlab pour le trai-
tement de ce type de problèmes : la taille des problèmes
devient trop importante pour pouvoir obtenir une bonne 3.3 Approche (( couplage fort )), formulation
précision de calcul avec les fonctions de base de Matlab4 . symétrique
Néanmoins, les résultats donnés par le code Matlab res-
tent très satisfaisants pour un calcul des premiers modes. La formulation du problème couplé utilisant la des-
La figure 11 donne les courbes de convergence des cription (u, p) (déplacement/pression) conduit à un
trois premières fréquences propres du système couplé problème aux valeurs propres non symétriques :
       
4 KS R T U(ω) MS 0 U(ω)
Une solution consisterait à développer certaines fonctions × = ω2 ×
du programme sous forme de sous routines C++ ou Fortran, 0 KF P(ω) −ρF RT MF P(ω)
compilées et intégrées au code Matlab [15]. (7)
 8

(interface avec la structure), les équations du problème

fluide sont alors :
∂2ϕ p
ρF + 2 = 0 dans ΩF (9a)
∂xi ∂xj c
p ω2
= ϕ dans ΩF (9b)
ρ F c2 c2
∂ϕ
· nj = 0 sur ∂ΩFn (9c)
∂xj
p = 0 sur ∂ΩFo (9d)

Fig. 11. Courbes de convergence des premières fréquences

Les équations de couplage entre les deux problèmes
propres du panier et de la cuve couplés avec le fluide. expriment la continuité de la composante normale du ten-
seur des contraintes et la continuité de la composante nor-
male du déplacement :

σij (u) · nj = ϕρF ω 2 ni sur ΓFS (10a)

∂ϕ
· nj = uj · nj sur ΓSF (10b)
∂xj

La formulation variationnelle du problème spectral

ainsi posé est alors la suivante : trouver u(ω), p(ω),
ϕ(ω) dans VS × VF × HF et ω réel tels que pour tout δu,
δp, δϕ dans VS × VF × HF , on ait [27–29] :

kS (u(ω), δu) + kF (p(ω), δp) =

ω 2 × [mS (u(ω), δu) − mF (ϕ(ω), δϕ) + r̃(u(ω), δϕ)

Fig. 12. Temps de calcul CPU pour l’extraction des valeurs
propres du système couplé avec deux plateformes de calcul. + r̃∗ (ϕ(ω), δu) + m̃F (p(ω), δϕ) + m̃∗F (ϕ(ω), δp)] (11)

Les espaces fonctionnels sont  respectivement

 [30] :
VS = u ∈ H 1 (ΩS )3 , u|∂ΩSo = 0 , VF = ϕ ∈ H 1 (ΩF )
 
L’utilisation d’une formulation du problème couplé et HF = p ∈ L2 (ΩF ), p|∂ΩFo = 0 ; les formes bi-
à l’aide d’une description (u, p, ϕ) (déplacement, pres- linéaires définies
 sur ces espaces sont données  par
sion, potentiel des déplacements [13]) permet d’obtenir kS (u, δu) = σij (u) · ·εij (δu)dΩ, mS (u, δu) = ρS u ·
un problème symétrique5. Les équations continues du ΩS
 pδp  ΩS
pδp
problème couplé sont alors les suivantes. ΩS désignant δudΩ, mF (p, δp) = c2 dΩ kF (p, δp) = ρF c2 dΩ,
le domaine structure, de frontières disjointes ∂ΩSn (ef-  pδϕ ΩF
 ΩF
fort imposé), ∂ΩSo (déplacement imposé) et ΓFS (inter- m̃F (p, δϕ) = c2 dΩ, r̃(u, δϕ) = ρF u · nδϕdΓ et
ΩF Γ
face avec le fluide), les équations du problème structure  pδϕ ∗
sont alors : m̃F (p, δϕ) = c2 dΩ ; m̃F et r̃∗ sont les opérateurs
ΩF
transposés de m̃F et r̃.
∂σij (u) La discrétisation de la formulation variationnelle (11)
ρS ω 2 u i + = 0 dans ΩS (8a) conduit au problème aux valeurs propres suivant :
∂xij
σij (u) · nSj = 0 sur ∂ΩSn (8b)  
 
KS 0 0  U(ω) 
 
ui = 0 sur ∂ΩSo (8c)    
 0 K F 0  × P(ω) =
   
0 0 0  Φ(ω) 
 
ΩF désignant le domaine fluide, de frontières ∂ΩFn
(déplacement imposé), ∂ΩFo (pression imposée) et ΓSF    
MS 0 R̃  U(ω) 
   

  P(ω)
ω2  0 0 M̃F × (12)
   
5
R̃T M̃T   Φ(ω) 

F −M F
D’autres formulations symétriques sont possibles, voir par
exemple [25, 26].
 9

Fig. 13. Structure des matrices masse et raideur pour le problème couplé panier, cuve et fluide en formulation (u, p, ϕ).

Tableau 7. Premières fréquences propres du problème couplé,

modes d’ordre n = 1. Comparaison des formulations (u, p)
et (u, p, ϕ).
Fréquence (Hz) Formulation Formulation ε (%)
(u, p) (u, p, ϕ)
f1 26,03 25,99 0,15
f2 65,62 64,72 1,37
f3 102,94 103,98 –1,01
f4 162,30 162,17 0,08
f5 180,32 183,70 –1,87

Après élimination de la variable Φ, le problème aux

valeurs propres (12) est écrit sous la forme symétrique Fig. 14. Premières fréquence propre du système couplé en
suivante : formulation (u, p, ϕ) et conditionnement en norme −1 de la
    matrice KF en fonction du nombre de degrés de liberté N du
KS 0 U(ω) problème.
  × =
0 K F P(ω)
    La formulation (u, p, ϕ) fait intervenir la ma-
−1 −1
MS + R̃M F R̃T R̃M F M̃T
F  U(ω)  trice K−1F : le calcul des valeurs propres du système couplé
ω2  −1 −1
× est donc influencé par la précision de l’inversion matri-
M̃F M F R̃T M̃F M F M̃T  P(ω) 
F cielle. Dans Matlab, nous utilisons la fonction inv, qui
(13) réalise une inversion matricielle à partir de l’algorithme de
Cholesky ou l’algorithme de Gauss selon la nature définie
En utilisant les mêmes fonctions de forme pour p et ϕ, positive de la matrice à inverser [20]. La figure 14 donne
en supposant les caractéristiques fluide c et ρF constantes, l’évolution de la première fréquence propre du système
nous avons : couplé en fonction de la taille du problème et l’évolution
MF du conditionnement de la matrice KF pour le problème
K F = R̃ = ρF R M F = ρF KF M F = MF
ρF de l’inversion matricielle et met en évidence la corrélation
(14) entre les deux évolutions.
Le problème (13) est symétrique mais la matrice de Si le fluide est considéré incompressible, c → +∞
masse n’est plus creuse, comme le montre la figure 13 et mF (., .) ≡ 0. Le problème couplé se réduit du point
(représentation des termes non nuls des matrices). de vue structure au problème symétrique suivant, faisant
Le tableau 7 compare les résultats de calculs obtenus intervenir l’opérateur de masse ajoutée [13] :
avec le code Matlab en utilisant les formulations (u, p) et
(u, p, ϕ). Les écarts entre les deux types de calcul sont de ω 2 (MS + MH ) × U(ω) = KS × U(ω) (15)
l’ordre de 1 %. Les modes propres obtenus à partir de la
avec :
formulation (u, p, ϕ) possèdent les propriétés d’orthogo-
MH = ρF RT K−1
F R (16)
nalité vis-à-vis des matrices masse et raideur du système
couplé, ce qui permet l’application de méthodes modales L’opérateur MH étant symétrique, le problème aux
pour les calculs dynamiques [31]. valeurs propres (15) l’est également.
 10

Tableau 8. Premières fréquences propres du problème couplé. La discrétisation du problème couplé a mis en œuvre
Comparaison entre un calcul couplé compressible et incom- des éléments finis axisymétriques développés en série de
pressible. Fourier :
Fréquence (Hz) Fluide Fluide – élément fini structure : coque élastique formulé en
compressible incompressible déplacement ;
f1 26,03 25,99 – élément fini fluide : fluide acoustique formulé en
f2 65,62 64,82
pression.
f3 102,94 104,23
f4 162,30 179,25 Un code de calcul dédié a été programmé sous Matlab
et réalise le maillage du problème, l’assemblage des ma-
Tableau 9. Calculs de la première fréquence propre pour trices élémentaires, la détermination des éléments propres
différentes situations de calcul (panier seul, cuve seule, fluide du système et la visualisation des modes propres.
seul, système couplé). Formulation fluide en pression, formu- Les comparaisons entre les résultats de cal-
lation structure en déplacement. cul Matlab et Ansys (utilisant un couplage
déplacement/déplacement) ont permis de valider les
Modèle Calcul Calcul ε (%) développements réalisés dans Matlab et la capacité du
Ansys 3D Matlab 2D
code élémentaire à traiter des géométries simples mais
Panier 41,25 44,70 8,36
représentatives de cas industriels.
Cuve 78,64 78,70 0,08
Fluide 224,68 224,68 –0,05 À terme, l’élément fluide ainsi développé et son
Panier, cuve et fluide 25,94 26,03 0,35 couplage avec des éléments structure sera implanté
dans le code Ansys pour permettre une analyse de
géométries complexes, en utilisant la formulation pres-
sion/déplacement dans le cas de géométries de révolution
Le tableau 8 donne les fréquences propres du système avec termes de chargement non axisymétriques.
couplé calculé dans le cas d’un fluide incompressible et
d’un fluide compressible : les fréquences obtenues et les
déformées modales du système structure sont identiques
dans les deux cas, ce qui montre que les effets de com- Annexe : Cas élémentaire de validation
pressibilité n’interviennent pas sur les premiers modes du du programme de calcul des modes propres
système couplé. Les fréquences propres du système couplé d’un problème généralisé non symétrique
avec fluide incompressible sont suffisamment éloignées des
premières fréquences propres de la cavité acoustique, ce
qui permet de négliger les effets acoustiques pour les pre- Le programme d’extraction des valeurs propres pour
miers modes [31]. un problème généralisé non symétrique utilisé dans le
code développé dans Matlab est une adaptation de l’al-
Nous comparons enfin les résultats de calcul obte- gorithme de Lanczos [23]. Le principe de l’algorithme
nus avec un modèle tridimensionnel Ansys utilisant les est détaillé dans [24] ; cet algorithme est implanté dans
éléments structure shell63 (éléments coques linéaires bi- le code de calcul Matlab pour traiter le cas industriel
dimensionnels [16]) et les éléments fluide fluid30 (éléments présenté, la validation de la fonction Matlab ainsi ob-
fluides acoustiques linéaires tridimensionnels [21]) et un tenue est faite à partir du problème présenté dans [24],
modèle bidimensionnel axisymétrique utilisant le code correspondant aux matrices 12 × 12 suivantes :
de calcul Matlab. Le tableau 9 donne une comparaison
des valeurs du premier mode calculé pour le panier, la
cuve, le fluide seuls ou couplés. Les résultats de calcul
sont équivalents dans les deux approches, mais les temps
de calcul sont très importants dans le cas d’un modèle
tridimensionnel, particulièrement pour le système couplé
compte tenu de la non symétrie du problème aux valeurs
propres.

Conclusion

La présente étude a porté sur l’analyse modale d’une

structure industrielle, en prenant en compte le couplage
fluide/structure. Une simplification de la géométrie du
problème réel nous a permis de nous placer dans le cas
d’un problème à géométrie de révolution, présentant un
chargement non axisymétrique.
 11

Références [14] O.C. Zienkiewicz, R.L. Taylor, The Finite Element

Method, Basic Formulation and Linear Problems, 4th edi-
[1] J. Makerle, Finite and Boundary Element Techniques in tion, Mac Graw Hill, 1989
Acoustics, a Bibliography (1990-1992), Finite Element [15] M. Mokhtari, A. Mesbah, Apprendre et maı̂triser Matlab,
Analysis Design 15 (1994) 263–272 Springer Verlag, 1997
[2] J. Makerle, Fluid-Structure Interaction Problems, Finite [16] P. Khonke, ANSYS Theory Reference, Swanson Analysis
Element Approach and Boundary Elements Approaches, System, 1986
A Bibliography, Finite Elements in Analysis and Design [17] J.F. Sigrist, C. Lainé, B. Peseux, Calcul couplé
31 (1999) 231–240 fluide/structure avec Ansys : étude numérique et
[3] G. Delhommeau, Contribution à l’étude théorique et à expérimentale d’une plaque élastique en contact avec un
la résolution numérique de problème de résistance de fluide lourd compressible, J. Physique France 12 (2002)
vagues, Thèse de Doctorat, École nationale supérieure de Pr11-137–144
mécanique de Nantes, 1978 [18] D.G. Ashwell, R.H. Gallagher, Finite Elements for Thin
[4] G. Delhommeau, Les problèmes de diffraction-radiation Shells and Curved Members, Wiley, 1976
et de résistance de vague : étude théorique et résolution [19] P.G. Ciarlet, Introduction à l’analyse numérique matri-
numérique par la méthode des singularités, Thèse de doc- cielle et à l’optimisation, Masson, 1988
torat d’État, École nationale supérieure de mécanique de [20] G.H. Golub, C.F. Van Loan, Matrix Computation, Johns
Nantes, 1987 Hopkins University Press, 1984
[5] H. Oudin, Étude du comportement hydroélastique [21] D.B. Woyak, Acoustic and Fluid Structure Interaction,
de structures marines par une formulation mixte Swanson Analysis System, 1995
équations intégrale-éléments finis, Thèse de doctorat, [22] J.F. Sigrist, C. Lainé, D. Lemoine, B. Peseux, Choice
École nationale supérieure de mécanique de Nantes, and limits of a linear fluid model for the numerical
1986 study in dynamic fluid structure interaction problem,
[6] B. Peseux, Contribution à l’étude de structures partielle- Pressure Vessel and Piping, Cleveland, 20–24 July 2003,
ment ou totalement immergées en matériau homogène ou 454, pp. 87–93
en composite, Thèse de doctorat d’État, École nationale [23] F. Chatelin, Valeurs propres de matrices, Masson, 1988
supérieure de mécaniques de Nantes, 1989 [24] C. Rajakumar, C.G. Rogers, The Lanczos Algorithm
[7] B. Peseux, G. Delhommeau, Fluid-Structure Interaction: Applied to Unsymmetric Generalized Eigenvalue
Influence of Compressibility on the Dynamic Response Problem, Int. J. Numerical Methods in Engineering 32
of Immersed Structures, Aéronautique Conférence (1991) 1009–1026
DGLR/AIAA, Strasbourg, 1993 [25] G.C. Everstine, Letter to the Editor: A Symmetric
[8] B. Peseux, F. Bérot, Vibro-Acoustic Behavior of Potential Formulation for Fluid-Structure Interaction, J.
Submersed Cylindrical Shells: Analytical Formu- Sound and Vibration 79 (1981) 157–190
lation and Numerical Model, J. Fluids and Structures 12 [26] G. Sandberg, P.A. Goransson, A Symmetric Finite
(1998) 959–1003 Element Formulation for Acoustic Fluid-Structure
[9] B. Peseux, A. Piard, Dynamic Response of Submersed Interaction Analysis, J. Sound and Vibration 123 (1988)
Structures in Bonded Domain with Gravity Waves, 507–515
Pressure Vessel and Piping, San Diego, 26–30 July 1998, [27] J. Boujot, Interaction fluide structure en régime transi-
336, pp. 187–196 toire, La Recherche Aérospatiale 3 (1984) 203–209
[10] J.P. Quevat, Couplage de la méthode des équations du second ordre correspondant à d i v e r s m o dèles de
intégrales et de la méthode des éléments finis pour la couplage fluide structure, Compte Rendu à l’Académie des
résolution d’un problème d’hydroélasticité, Thèse de doc- Sciences, I, 301 (1985) 391–394
torat, École nationale supérieure de mécanique de Nantes, [28] J. Boujot, Etude d’une famille de problèmes d’évolution
1979 [29] J. Boujot, Mathematical Formulation of Fluid Structure
[11] J. Makerle, Finite Element in the Analysis of Interaction Problems, Mathematical Modelling and
Pressure Vessel and Piping, a Bibliography (1976- Numerical Analysis 21 (1987) 239–260
1996), Int. J. Pressure Vessel and Piping 69 (1996) [30] P.A. Raviart, J.M. Thomas, Introduction à l’analyse
279–339 numérique des équations aux dérivées partielles, Masson,
[12] J. Makerle, Earthquake Analysis of Structures: FEM and 1983
BEM Approaches, a Bibliography (1995-1998), Finite [31] R.J. Gibert, Vibration des structures, Interaction avec
Element Analysis Design 32 (1999) 113–124 les fluides, Sources d’excitation aléatoires, Collection de
[13] H.J.P. Morand, R. Ohayon, Fluid Structure Interaction, la Direction des Études et Recherches d’Électricité de
Wiley & Sons, 1995 France, 69 (1986)