Ibssmfra 14 Asm 351

Nomenclature
Symboles Nomination Unités

𝝈𝒄 le tenseur symétrique des contraintes de Cauchy. N/m2
𝒇 le vecteur des répartitions volumiques de forces.
𝒅𝟐 𝑼 l’accélération du point 𝑴, dont le champ des
𝚼 𝑴 = 𝟐 m/S2
𝒅𝒕 déplacements.
𝑰 Matrice unité.
(𝝀 , 𝝁) Coefficient de Lamé.
𝝂 Coefficient de Poisson.
𝑬 Module de Young. GPa
𝛔 vecteur des contraintes. N/m2
𝐃 matrice d’élasticité.
𝛆 vecteur des déformations N/m2
𝑨 l’aire de la section droite. m2
𝑮 le module d’élasticité transversal.
𝒌𝒚 le coefficient d’aire cisaillée.
𝑻𝒚 l’effort tranchant. N
𝑴 masse de système. Kg
𝝆𝒇 masse volumique du fluide. Kg/m3
𝝆 masse volumique du solide. Kg/m3
𝒏 lanormale extérieure unitaire.
𝚪𝐒𝐅 le domaine interaction fluide-structure. m2
𝒖 la vitesse du fluide. m/s
𝒑 la pression Pa
𝑰 la matrice identité.
𝝈 le tenseur des contraintes visqueuses. N/m2
𝚽 potentiel du vitesse.
𝚪𝐋 surface libre m2
Ω𝑭 domaine du fluide. m2
Ω𝑺 domaine du solide. m2
∫ intégrale.
{ } matrice colonne : vecteur
[ ] matrice (utilisée aussi pour les références).
𝒗 𝒛, 𝒕 la flèche. m
𝑵𝒔 fonction de forme.
𝒂𝒆 déplacement nodaux. m
𝑲𝒆 matrice de rigidité élémentaire. Kg/S2
𝑴𝒆 matrice de masse élémentaire. Kg
𝜞∞ frontière à l’infini.
Page 1
Nomenclature
𝜞𝑹 frontière fini.
𝚽 𝒏 vecteur des potentiels nodaux élémentaires
𝒑 vecteur des pressions latérales des poutres Pa
𝑭𝑺 matrice globale de couplage fluide-structure.
𝑻
𝑭𝑺 matrice globale transposée de couplage fluide-structure.
𝑴𝒂 matrice de masse ajoutée Kg
matrice du système en tenant compte des conditions de
[𝑯𝒓 ] pression nulles sur la surface libre
𝒓𝟎 , 𝒛𝟎 les coordonnés du centre de l’élément quadratique𝑸𝟒 m
Page 2
Liste des tableaux
Tableau VI-1 :Convergence des fréquences (Hz) en fonction du nombre d’éléments……59
Tableau VI-2 :Valeurs des fréquences naturelles (Hz)………………………………...…..60
Tableau VI-3a :Variation des fréquences pour une même valeur d’inertie pour des sections
différentes…………………………………………………………………………………..………..61
Tableau VI-3b :Variation des fréquences pour une même valeur de surface pour des sections
différentes……………………………………………………………………………………62
Tableau VI-4 : Variation des fréquences en fonction du rapport L/r pour une poutre de
section
circulaire……………………………………………………………………………………....64
Tableau VI-5 :Variation des fréquences en fonction du nombre de masse Rm=rhos/rhof….66
Tableau VI-6 : Variation des fréquences pour différents matériaux………………………...67
Tableau VI-7 :Variation des fréquences pour différents matériaux………………………...68

Liste des figures
Figure I-1 : couplage structure-fluide et couplage structure fluide-fluide………………………...7
FigureI-2 :Solide élastique immergé dans un fluide……………………………………………...7
Figure I-3 : Couplage SPH/FEM…………………………………………………………………..8
Figure I-4 :Couplage SPH/FEM/FEM………………………………………………………….…8
Figure I-5 : Couplage FEM/FEM………………………………………………………………….9
Figure II.1 : Elément de poutre dans le plan XY…………………………………………………..12
Figure II.2:Schématisation générale du problème structure……………………………………..15
Figure II.3: Déplacement de la section droite d’abscisse x……………………………...……….16
Figure II.4:Hypothèse de Bernoulli...…………………………………………...…..……………18
FigureII.5:Schématisation générale du problème structure……………………………….............25
Figure II.6: Schématisation générale du problème structure-fluide………………………………25
Figure III.1:Schématisation générale d’un problème de couplage mécanique fluide-structure.....28
Figure III.2:Schématisation générale du problème fluide-structure……...………………………29
Figure III.3:Principe de maillage du problème couplé dans le plan (𝑟, 𝑧)…….………………..36
Figure IV.1:Poutre vertical immergé dans un fluide…………………………………………......42
Figure IV.2:élément poutre à deux nœuds………………………………………………………..44
Figure IV.3: maillage du domaine fluide………………………………………………………….47
Figure V.1:Schéma de maillage et de numérotation de la structure………………………………52
Figure VI-1 : poutreconsole partiellement immergée…………………………………………….58
Figure VI -2 : courbes de convergence…………………………………………………………....59
Figure VI -3a : Variation des fréquences en fonction du rapport des surfaces pour différentes
sections pour la même valeur d’inertie, (a) mode 1, (b) mode 2…………………………………..61
Figure VI -3b : Variation des fréquences en fonction du rapport d’inertie des différentes sections
pour la même valeur de surface, (a) mode 1, (b) mode 2………………………………………….63
Figure VI- 4 :Tracés des graphes des quatre modes pour les rapports L/r= 2, 2.2, 4, 5, 10 (a :
mode1, b : mode2, c : mode3, d : mode4)………………………………………………………..…65
Figure VI-5: variation des fréquences en fonction du rapport de masse volumique de l’acier / masse
volumique du cuivre et laiton………………………………………………………………..……...67
Figure VI-6:Variation des fréquences en fonction du rapport de masse (Rm=rhos/rhof) et du rapport

hauteur d’immersion sur la longueur de la poutre H/L (a) :mode1,(b) :mode……………..69
Sommaire
Introduction générale………………………………………………………………….1
Chapitre I: Revubibliographique….…………………………………………………4
Chapitre II : Formulation mathématique
II.1Introduction……………………………………………………..............................11
II.2Equations de solide………………………………………………………………..11
II-2-1Equations locales…..........................................................................................12
II-2-2 Conditions aux limites……...……………………………………………...…15
II-3 Théorie des poutres……………………………………………………….………16

II- 3-1 Généralités…………………………………………………………………16
II- 3-2 Equations d’équilibre……………………………………………………...17
II.4 Mécanique des fluides …..………………………………………………………..18
II.5 Couplage…………….…………………………………………………………….24
Chapitre III :Modélisatuon par élément finis et couplage
III.1 Introduction…………………………………………………………………….28
III.2 Méthode des éléments finis pour l’analyse vibratoire d’une structure couplée en
présence d’un fluide incompressible ………………………………………………....28
III.3 Les éléments finis en mécanique du solide ……………………………………...30
III.4 Méthode des éléments finis en mécanique des fluides.(Forme faible de la
méthode des résidus pondérés)………………………………………………………...33
III.5 Couplage éléments finis-éléments finis………………………………………….35
Chapitre IV:Cas d'une poutre émergée dans un fluide
IV.1Etude des vibrations d’une poutre immergée dans un fluide…………………….42
IV.1.1Equation du problème couplé ………………………………………………43
IV.1.1.1Equations du comportement de poutre……………………………….....43

IV.1.1.2Equations du fluide……………………………………………………….43
IV.2Modélisation numérique et couplage………………………………………...…44
IV.2.1 Discrétisation de la poutre……………………………………………………44
IV.2.2 Discrétisation du milieu fluide……………………………………………….45
Chapitre V:Organisation de la programmation
V.1 Introduction……………………………………………………………………....51
V-2Description du programme………………………………………………………..51
V-3 Résolution………………………………………………………………………....53
V-4 Organigramme……………………………………………………………….........55
Chapitre VI :Résultat et interprétation
VI.1 Introduction………………………………………………………………………57
VI-2 Description de la poutre………………………………………………………….57
VI-3 Validation du programme……………………………………………...…58
VI-4 Influence des paramètres géométriques et physiques sur les fréquences propres
de vibration de la poutre………………………………………………………………60
Conclusion…………………………………………………………............................70
Références bibliographiques……………………………………………….......... 71
Résumé
L’étude présentée dans ce mémoire concerne l’analyse du comportement dynamique d’une

poutre immergée totalement ou partiellement dans le fluide.
Le mouvement de la poutre dans le fluide conduit a une interaction fluide-structure qui donne
une augmentation de la masse du système dite masse ajoutée.
Les lois de comportement ont été établies par la relation déplacements-contraintes pour la
poutre et les équations de conservation de la masse et de la quantité du mouvement pour le
fluide avec la condition de non glissement au niveau de l’interface.
Le couplage définit par l’interaction fluide-structure nécessite la caractérisation de condition

aux limites qui décrivent le dialogue entre le solide et le fluide au niveau de l’interface. La
modélisation de la poutre et le fluide est faite par la méthode de l’élémentfinis.
Les équations différentielles du mouvement de la structure sans contact avec le fluide ont été
établies. La détermination des équations de couplage nous ont permis de calculer la masse
ajoutée générée par le fluide, qui nous a permisd’aboutir à un système d’équations
différentielles qui caractérise le comportement dynamique des deux milieux.
Les fréquences propres de la poutre immergée correspondant aux quatres premiers modes sont
calculées à l’aide d’un programme élaboré sous Matlab.
Après avoir étudié la convergence et validé le programme avec plusieurs littératures, on a

étudié plusieurs exemples, cette étude nous a permis de déterminer l’influence des paramètres
physiques et géométriques de la poutre immergée.
Parmi ces cas d’étude, nous décrivons la variation de la longueur d’immersion de la poutre
dans le fluide, des rapports de la longueur sur le diamètre de la poutre du rapport de la masse
volumique du solide sur cellede fluide, longueur de la poutre sur la hauteur d’immersion et en
même temps rapport masse volumique du solide sur la masse volumique du fluide et en fin le
rapport d’inertie des différentes sections pour une même surface, puis le rapport de surface de
différentes section pour une même inertie.
Mots clés : éléments finis, couplage, interaction, masse ajoutée, équation différentiel, poutre
immergée.
Abstract
The study presented in this memory relates to the analysis of the dynamic behavior of a beam
immersed completely or partially in the fluid.
The movement of the beam in the fluid leads has an interaction fluid-structure which gives an
increase in the mass of the system known as added mass.
The laws of behavior were established by the relation displacement-constraints for the beam
and the conservation equations of the mass and the quantity of the movement for the fluid
with the condition of nonslip on the level of the interface.
The coupling defines by the interaction fluid-structure requires the characterization of

boundary condition which describes the dialogue between the solid and the fluid on the level
of the interface. The modeling of the beam and the fluid are made by the method the finite
element.
The differential equations of the movement of the structure without contact with the fluid
were established. The determination of the equations of coupling enabled us to calculate the
added mass generated by the fluid, which enabled us to lead to a system of differential
equations which characterizes the dynamic behavior of the two mediums.
The Eigen frequencies of the immersed beam corresponding to the fuor first modes are
calculated using a program worked out under Matlab.
After having studied convergence and having validated the program with several literatures,
one studied several examples, this study enabled us to determine the influence of the physical
and geometrical parameters of the immersed beam.
Among these cases of study, we describe the variation in the length of immersion of the beam
in the fluid, the reports/ratios length on the diameter of the beam of the ratio of the density of
the solid on that of fluid, length of the beam on the height of immersion and at the same time
density report/ratio of the solid on the density of the fluid and in end the ratio of inertia of the
various sections for the same surface, then the report/ratio of surface of different section for
the same inertia.
Key words:finite elements, coupling, interaction, added mass, differential equation,

immersed beam.
‫يهخص‬
‫انذساسخ انز‪ ٙ‬لذيذ ف‪ْ ٙ‬زِ انًزكشح رزؼهك ثزحه‪ٛ‬م انسهٕن انذ‪ُٚ‬بي‪ٛ‬ك‪ ٙ‬نؼزجخ يغًٕسح خزئ‪ٛ‬ب أ كه‪ٛ‬ب ف‪ ٙ‬سبئم‪.‬‬
‫حشكخ انؼزجخ ف‪ ٙ‬انسبئم رؤد٘ انٗ رفبػم ػزجخ– سبئم انز٘ ‪ٚ‬ؼط‪ُٛ‬ب رزا‪ٚ‬ذ ف‪ ٙ‬كزهخ انُظبو ٔ انز‪ ٙ‬رسًٗ انكزهخ انًضبفخ‬
‫لٕاَ‪ ٍٛ‬انسهٕن انحشك‪ ٙ‬يؼطبح ثؼاللخ االخٓبد انًزؼهمخ ثزُمم انؼزجخ ٔ يؼبدنخ حفظ انكزهخ ٔ كً‪ٛ‬خ انسبئم يغ ػذو ٔخٕد‬
‫حشكخ االَزالق ػهٗ يسزٕٖ انٕاخٓخ‪.‬‬
‫االلزشاٌ انًؼشف ثبنزفبػم سبئم‪-‬ػزجخ ‪ٚ‬زطهت رٕص‪ٛ‬ف انششٔط انحذ‪ٚ‬خ انز‪ ٙ‬رصف انحٕاس ث‪ ٍٛ‬انصهت ٔ انسبئم ػهٗ‬
‫يسزٕٖ انٕاخٓخ‪.‬‬
‫‪ٚ‬زى ٔضغ ًَبرج انسبئم ٔ انصهت ثٕاسطخ طش‪ٚ‬مخ انؼُبصش انًزُبْ‪ٛ‬خ‪.‬‬
‫لذ اَشئذ انًؼبدالد انزفبضه‪ٛ‬خ نحشكخ انصهت دٌٔ ارصبنّ يغ انسبئم‪.‬‬
‫االلزشاٌ ‪ٚ‬سًح نُب ثزحذ‪ٚ‬ذ انًؼبدالد نحسبة انكزهخ انًضبفخ انُبردخ ػٍ انسٕائم االيش انز٘ ‪ٚ‬سًح نُب ثزحم‪ٛ‬ك َظبو يٍ‬
‫انًؼبدالد انزفبضه‪ٛ‬خ انز‪ ٙ‬رً‪ٛ‬ز انسهٕن انذ‪ُٚ‬بي‪ٛ‬ك‪ ٙ‬نهٕسط‪.ٍٛ‬‬
‫انزشدداد انطج‪ٛ‬ؼ‪ٛ‬خ نهؼزجخ انًغًٕسح رٕافك االًَبط االسثؼخ االٔنٗ انًحسٕثخ ثًسبػذح ثشَبيح *يبطالة*‬
‫ثؼذ دساسخ انزمبسة ٔ رطج‪ٛ‬ك انجشَبيح دسسُب ايثهخ ػذ‪ٚ‬ذح ٔ ْزا سًح نُب ثبنزؼشف ػهٗ ربث‪ٛ‬ش انًؼبيالد انف‪ٛ‬ز‪ٚ‬بئ‪ٛ‬خ ٔ‬
‫انُٓذس‪ٛ‬خ نؼزجخ يغًٕسح‪.‬‬
‫يٍ ث‪ ٍٛ‬انحبالد انًذسٔسخ َحذد طٕل اغًبس انؼزجخ ف‪ ٙ‬انسبئم ‪َ.‬ست انطٕل ػهٗ لطش انؼزجخ ٔ َست انكزهخ انحدً‪ٛ‬خ‬
‫ػهٗ َظ‪ٛ‬شرٓب نهسبئم‪.‬طٕل انؼزجخ ػهٗ اسرفبع االغًبس ٔ ف‪َ ٙ‬فس انٕلذ َست انكزم انحدً‪ٛ‬خ نهصهت ػهٗ انكزهخ‬
‫انحدً‪ٛ‬خ نهسبئم ٔ ف‪ ٙ‬االخ‪ٛ‬ش َسجخ انؼطبنخ الشكبل يخزهفخ نُفس انًسبحخ ثى َسجخ انًسبحخ نًخزهف االشكبل نُفس‬
‫انؼطبنخ‪.‬‬
‫انًفبر‪ٛ‬ح‬
‫انؼُبصش انًزُبْ‪ٛ‬خ ‪.‬االلزشاٌ‪.‬انزفبػم‪.‬انكزهخ انًضبفخ‪.‬انًؼبدالد انزفبضه‪ٛ‬خ‪.‬انؼزجخ انًغًٕسح‪.‬‬

INTRODUCTION GENERALE
Introduction générale
L’interaction fluide-structure est un domaine vaste il concerne un ensemble important de

systèmes mécaniques ou de phénomènes naturels. Ces phénomènes apparaissent Lorsque on
est en présence de deux milieux différents ; c’est dans le domaine naval ou le milieu fluide
correspond à l’eau, que nous avons une large gamme de cas d’interaction allant du sous marin
aux bateaux. Si le fluide considéré est de l’air, nous sommes dans le domaine de
l’aérodynamique
Ces différents problèmes sont donc complexes à étudier. Certaines simplifications peuvent
cependant être mises en œuvre. Dans notre cas il s’agit d’étudier les caractéristiques
dynamiques d’une poutre partiellement ou totalement immergée soumise à des vibrations. Le
fluide considéré est un fluide parfait incompressible et irrotationnel.
Dés l’aube de l’industrialisation du monde moderne, le problème de l’interaction fluide-

structure est apparut dans de très nombreuses applications telle que l’aviation, les
constructions navales et les off-shores etc. Cette branche hybride entre la mécanique des
fluides et celle des structures permet de quantifier les interactions fluide-structure.
Les phénomènes sont d’autant plus significatifs aujourd’hui qu’il est nécessaire et important
de pouvoir les modéliser correctement dés la phase avant projet.
Ces phénomènes d’interaction entre les fluides et la structure sont nombreux. Ils peuvent être
classés par la nature des échanges entre les deux milieux :
 Couplage thermique : le fluide et le solide échangent de la chaleur.
 Couplage chimique : le solide et le fluide réagissent entre eux.
 Couplage mécanique : le solide et le fluide échangent des énergies mécaniques.
Notre étude fait partie de ce dernier type d’interaction, nous étudions le comportement
élastique d’une structure « poutre » immergée. Il faut noter que la résolution de ce genre de
problème de couplage dépend du milieu que l’on veut privilégier (le fluide ou bien la
structure), le transfert d’informations se fera du milieu secondaire vers le milieu primaire en
d’autres termes des échanges d’énergie mécanique se produisent dans les deux sens :
Page 1
- Le fluide exerce des efforts sur la structure en mouvement, ce qui modifie la

dynamique de la structure.
- La structure impose des déplacements sur le fluide et modifie ainsi les caractéristiques
de l’écoulement du fluide.
Ce couplage s’opère au niveau de l’interface, il induit sur la structure une certaine force
d’inertie et en utilisant le principe de la continuité des contraintes et des vitesses à
l’interface fluide-structure, on montrera que l’effet du fluide sur le solide génère une masse
supplémentaire que l’on appelle « masse ajoutée » qui est liée au système mécanique et se
déplace solidairement avec lui. Autrement elle sera additionnée à la masse du solide. Il est
nécessaire d’établir une formulation adaptée pour la résolution des équations différentielles
régissant le système structure-fluide, en utilisant des méthodes numériques. Notre choix fut
porté sur le couplage « élément finis- éléments finis» ,et finalement on est ramené a résoudre
le système des équations différentielles aux valeurs propres.
Dans le présent travail la structure sera considérée comme milieu primaire, ainsi les variables
relatives au fluide seront éliminées au profit des variables de la structure vibrante.
Notre travail se divise en quatre parties principales:
Dans le premier chapitre, seront exposées les équations physiques générales régissant la
dynamique d’un corps solide, suivi d’un rappel sur la théorie des poutres, ainsi que les lois de
la mécanique des fluides. Par la suite nous poserons les hypothèses de petits déplacements et
petites déformations.
Dans le second chapitre, nous procédons à La formulation du problème vibratoire d’une

poutre en établissant le système des équations différentielles. Les matrices de masse et de
rigidité sont calculées par intégration exacte en se basant sur la théorie classique et en utilisant
la méthode des éléments finis.
La méthode de couplage établie nous a permit de calculer la masse ajoutée due au fluide. En
additionnant cette grandeur à la masse du solide dans le système d’équations différentielles
relatives à la poutre non immergée, nous obtenons un nouveau système d’équations
différentielles permettant la caractérisation dynamique de la poutre immergée dans le fluide,
aboutissant finalement à un problème vibratoire aux valeurs propres.
Page 2
Dans le chapitre quatre, nous avons présenté une description du programme élaboré sous
MATLAB.
Et enfin dans le dernier chapitre nous avons validé notre programme en montrant que
l’élément poutre utilisé converge. Puis nous avons étudié l’influence des paramètres
géométriques et de masses à savoir le rapport rayon sur longueur, hauteur d’immersion sur la
longueur, et rapport de la masse volumique du solide sur la masse volumique du fluide, sur le
comportement dynamique de la poutre.
Nous terminons par une conclusion générale et des perspectives futures à ce travail.
Page 3
REVUES BIBLIOGRAPHIQUES
CHAPITRE I
Couplage fluide-structure, domaine de la mécanique aux multiples finalités, regroupe

deux grandes types de problèmes : ceux ou le fluide est en écoulement par rapport à la
structure et ceux ou fluide et structure évoluentautour d’une configuration statique
initiale.
L’interaction structure-fluide apparait systématiquement dès qu’une structure se trouve en

contact avec un fluide.
Le mouvement des particules du fluide, au voisinage de la structure génère une variation

de pression.
L’étude de couplage a enrichi considérablement la littérature scientifique et technique

depuis plus d’un demi-siècle et ce parallèlement au développement et à l’élaboration des
grandes défis technologique qui ont marqué cette période (conquêtede l’espace,
développementdes moyens de transport aériens…).
La paternité de la théorie des poutres est attribuée à Galilée, mais des études récentes
indiquent que Léonard de Vinci l'aurait précédé. Léonard avait supposé que la
déformation variait de manière linéaire en s'éloignant de la surface neutre, le coefficient
de proportionnalité étant la courbure, mais il ne put finaliser ses calculs car il n'avait pas
imaginé la loi de Hooke. De son côté, Galilée était parti sur une hypothèse incorrecte (il
supposait que la contrainte était répartie uniformément en flexion), et c'est Antoine Parent
qui obtint la distribution correcte.
Ce sont Leonhard EULER et Jacques Bernoulli qui émirent la première théorie utile vers
1750, tandis que Daniel Bernoulli, le neveu du précédent, écrivit l'équation différentielle
pour l'analyse vibratoire. À cette époque, le génie mécanique n'était pas considéré comme
une science, et l'on ne considérait pas que les travaux d'une académie des mathématiques
puissent avoir des applications pratiques, et l'on continua à bâtir les ponts et les bâtiments
de manière empirique.
Ce n'est qu'au XIXe siècle, avec la Tour Eiffel et les grandes roues, que l'on démontra la
validité de la théorie à grande échelle.
Une formulation exacte du problème de poutre était étudié en premier lieu en terme
d’équations générales d’élasticité par Pochhammar(1876)et Chree(1889)[1].
Page 4
CHAPITRE I
La théorie de poutre de Rayleigh (1877) [ 2 ] fournit une amélioration marginale sur la

théorie d'Euler Bernoulli en incluant les effets de la rotation de la section transversale. En
conséquence, elle corrige partiellement la surestimation des fréquences normales dans le
modèle d'Euler Bernoulli. Cependant, les fréquences normales sont encore surestimées.
Les premiers investigateurs incluent Davies (1937) [3], qui a étudié les effets de l'inertie
rotatoire sur une poutre encastré -libre.
Timoshenko (1921, 1922) [4, 5] a proposé une théorie de poutre qui ajoute l'effet du
cisaillement comme l'effet de la rotation a la poutre d'Euler Bernoulli. Le modèle de
Timoshenko est une amélioration importante pour les poutres non-minces et pour des
réponses à haute fréquence où le cisaillement ou les effets rotatoires ne sont pas
négligeables. Après Timoshenko, plusieurs auteurs ont obtenu les équations de fréquence.
Pour le problème interaction fluide-structure il y’avait deux approches théoriques la
théorie de Wagner qui dérive de celle de VonKarmanet la théorie non linéaire des
similitudes Dobrovol’skaya,Pierson,Cumberbatch,.
Bien qu’il soit difficile d’étendre ces théories à des cas très complexes, elles donnent des
résultatsde références sur lesquels les résultats des méthodes numériques peuvent
compter pour leurs validations.
Pour les méthodes numériques, les difficultés résident dans le suivi de la surface libre et
dans la gestion de l’interaction entre la formulation fluide et la formulation structure.
Pour cette dernière, une formulation Lagrangienne en élément finis est suffisante pour
bien la modéliser [6]. Depuis la fin des années 50, le phénomène de ballotement dans les
réservoirs desengins spatiaux et aériens préoccupe le National AdvisoryComittee for
Aeronautics(NACA) américain que dés 1958 effectue des compagnes d’essais pour
évaluer l’importance de ce phénomène sur la stabilité en vol des avions [7].
En1962,Abranson qui fait figure de poinnierdans l’étude du ballottement, parallèlement
les Russes : Rapoport,Moiseyev et Rumyantrev effectuent des travaux traduits et publiés
en 1968, concernant la stabilité des structures contenant un fluide [8].
Parmi les précurseurs de l’analyse des phénomènes hydro-élastique, citons par exemple
Zienkieviez et al en 1968,pour ses études concernant les barrages avec l’une des
premières introductions au concept de matière de masse ajoutée[9].Belytshko pour ses
travaux relatifs au domaine nucléaire dés 1977.Citons également les méthodes
Page 5
CHAPITRE I
numériques développées par Coppolino en 1976 dans le logiciel MSC/NASTRAN. Pour

le calcul hydro-élastique par éléments finis des lancers axisymétriques.
En 1966, Tong propose une première formulation variationelle symétrique adaptée aux
méthodes éléments finis, comportait cependant certaines approximations.
Kruszewski a obtenu les trois modes antisymétriques d’une poutre à porte-faux et trois
modes antisymétriques et symétriques d'une poutre libre.Traill-Nash et collier ont donné
un traitement théorique assez complet aussi bien que des résultats expérimentaux pour le
cas d'une poutre uniforme.
Nous comptons utiliser la méthode des éléments finis(MEF) pour l’analyse de la structure
et de déterminer la masse ajoutée de fluide pour nous permettre de réaliser le couplage
.Cette méthode a été utilisée et développé par Peseux [10].
Jean-François Sigrist, Christian Lainé et Bernard Peseux [11],[12] utilise une analyse
modale d’une structure industrielle couplée avec un fluide, avec les techniques
numériques de calculs couples fluide/structure. La nature de la structure et axisymétrique
en géométrie et, la modélisation du problème est réalisée par les éléments finis
axisymétriques développes en série de Fourier par la discrétisation en éléments finis.
La modélisation est effectuée par un code de calcul implanté dans Matlab pour permettre
l’analyse modale de la structure, la comparaison se fait par ANSYS à des résultats de
calcul Matlab avec les résultats de calculs obtenus avec le code généraliste Ansys, pour
l’étude de problèmes couples en pression/déplacement avec développement en série de
Fourier, a partir de ca trouver les valeurs propre par Matlab et le code de calcule Ansys.
Myung Jo Jhung,WalTaekimet Yong Ho Ryu [13] ont déterminé la fiabilité de
cylindres ou de réservoirs avec anneaux rempli de liquide, sur les caractéristiques
modales en considérant l’effet d'interaction fluide-structure. Ils utilisent la méthode des
éléments finis pour la structure et la méthode de masse ajoutée à la masse de la structure
pour simplifier la résolution du problème et obtenir sous forme de donnée le déplacement
et la fréquence naturelle de remplissage du cylindre.
Le couplage entre un fluide et une structure est la base de beaucoup de phénomènes.
Parmi c'est phénomènes, nous pouvons citer la déformation de la structure, qui peut
Page 6
CHAPITRE I
trouver dans beaucoup de domaines, même dans le corps humain qui est une structure
élastique complexe. On trouve plusieurs définitions sur le couplage fluide-structure, soit
mathématique ou physique, c'est-à-dire la formulation algébrique ou physique sur les
bases des données du problème à résoudre.
Le couplage est l’étude de l’interaction crée au niveau de l’interface entre deux propriétés
physiques différentes ou de même propriété physique, comme solide-solide, fluide-
structure, fluide-fluide.
Les figures suivantes montrent les différents cas de couplage fluide structure qu'on peut
trouver.
Figure I-1: couplage structure-fluide et couplage structure- fluide – fluide.
Figure I-2: Solide élastique immergé dans un fluide.
Page 7
CHAPITRE I
Parmi les méthodes utilisées pour la prise en compte du couplage, nous pouvons citer la
méthode par pénalisation.
La méthode SPH (SmoothedParticleHydrodynamics), qui est une méthode par
pénalisation se caractérise par la masse des particules et par deux grandeurs de
discrétisation, la distance entre les particules, diamètre de la particule.
Le couplage SPH/FEM repose sur une méthode de contact par pénalisation. Les
particules fluides sont les "nœuds" esclaves. La particule pénètre dans un élément
structure d’une distance d (voir Figure 3), la force appliquée est proportionnelle à cette
distance: 𝑓 = 𝑘 × 𝑑dans laquelle k est la raideur par pénalisation.
Figure I-3: Couplage SPH/FEM.
Ce couplage est utilisé pour gérer l’interaction fluide/structure entre un maillage

Lagrangien, pour la structure, et un maillage Eulérien pour le fluide. Le maillage
Lagrangien étant immergé dans la grille Eulérienne.
Figure I-4: Couplage SPH/FEM/FEM.
La distance relative entre les positions Lagrangiennes des nœuds mais la distance relative
entre une particule fluide et un point interne de l’élément de la structure. Les quantités
Page 8
CHAPITRE I
aux points internes telles que les forces d’interaction, les masses ou les vitesses sont
reliées aux quantités nodales par les fonctions de forme NI.
Parmi les travaux nous avons aussi la méthode de FEM/FEM par pénalisation qui se
caractérise par une discrétisation par éléments finis de la structure et du fluide et
l’interface qui est caractérisée par une raideur entre deux nœud sous forme de contact
sans frottement ni glissement (voir Figure .5).
FigureI- 5: Couplage FEM/FEM.
La méthode des multiplicateurs de Lagrange est plus naturelle que la méthode des
pénalités, mais nécessite un temps de calcul plus élevé.
L’idée générale est de trouver, en respectant la contrainte spatiale issue des positions des
nœuds esclaves et maîtres, la force de contact vérifiant au mieux l‘équation de
conservation de quantité de mouvement.
Page 9
Page 10
FORMULATION MATHEMATIQUE
CHAPITRE II
II-1. Introduction
L’étude du mouvement d’un corps fait introduire divers paramètres qui le décrivent,
d’où la nécessité d’examiner la conservation de certaines quantités pendant ce
mouvement, comme la masse, la quantité de mouvement et le moment cinétique. Dans
ce contexte, et face à un problème typique de la mécanique des solides, on s’intéresse
au calcul des déplacements, des déformations et des contraintes. La contrainte est le
paramètre qui caractérise l’interaction mécanique d’un corps avec son environnement.
Cette contrainte est la cause du mouvement, d’où la nécessité d’une formulation
mathématique. Plusieurs théories ont été développées et sont actuellement utilisées en
ingénierie. Chacune de ces approches vise à modéliser certains aspects spécifiques du
comportement des matériaux.
La première partie de ce chapitre est consacrée à la description mathématique des

équations générales qui traduisent le comportement du solide, représenté par les deux
modèles d’états de contraintes planes ou de déformations planes, ces équations sont
appliquées dans les cas des poutres. La deuxième part ie est consacrée aux équations
régissant le milieu fluide, suivie d’une troisième consacrée au couplage des deux
milieux.
II-2. Equations de solide
II.2.1 Equations locales :

On prend l’hypothèse des petites perturbations en élasticité linéaire [14], on se place
dans un repère galiléen ou absolu. A chaque instant t, pour toute partie 𝜴𝒔 d’un
système matériel 𝑺, le torseur dynamique de 𝜴𝒔 est égal au torseur des forces
extérieures s’exerçant sur 𝜴𝒔 .
 les équations d’équilibre s’écrivent :
𝒅𝒊𝒗𝝈𝒄 + 𝒇 = 𝝆𝒔 𝜸 𝒅𝒂𝒏𝒔 𝜴𝒔 (2-1)
Où :
o 𝝈𝒄 est le tenseur symétrique des contraintes de Cauchy,
Page 11
CHAPITRE II
o 𝒇 est le vecteur des répartitions volumiques de forces,

𝒅𝟐 𝑼
o 𝚼 𝑴 = est l’accélération du point 𝑴.
𝒅𝒕𝟐
dont le champ des déplacements est:
𝑼 𝑴 = 𝒖 𝑴 𝒆𝟏 + 𝒗 𝑴 𝒆𝟐 + 𝒘 𝑴 𝒆𝟑 (2-2)
z,w x,u
y,v u1 u2
v1 v2
θ1 θ2
e2 o x,u
e3 e1
θ
z,w
Figure II-1 :Elément de poutre dans le plan XY
 Lois de comportement
Page 12
CHAPITRE II
Le tenseur des déformations associé est donné par [14]
𝟏 𝒕
𝜺 = 𝜺𝒊𝒋 (𝑼) = 𝒈𝒓𝒂𝒅𝑼 𝑴 + 𝒈𝒓𝒂𝒅 𝑼 𝑴 (2-3)
𝟐
Où le gradient du déplacement est un tenseur symétrique.
Les déformations de la structure induisent un champ de contraintes dans la structure,

représentées par le tenseur des contraintes 𝝈𝒊𝒋 (𝑼), l’hypothèse du comportement
élastique de la structure permet d’écrire la relation entre les déformations et les
contraintes sous la forme suivante (Loi de Hooke) :
(2-4)
𝝈𝒊𝒋 𝑼 = 𝝀𝜺𝒌𝒌 𝑼 𝜹𝒊𝒋 + 𝟐𝝁𝜺𝒊𝒋 (𝑼)
D’où les équations définissant la loi de comportement
−𝝀 𝝈 𝝂 𝟏+𝝂 (2-5a)
𝜺= 𝒕𝒓𝝈𝑰 + = − 𝒕𝒓𝝈𝑰 + 𝝈
𝟐𝝁 𝟐𝝀 + 𝟐𝝁 𝟐𝝁 𝑬 𝑬
Ou son inverse
𝝂𝑬 𝑬 (2-5b)
𝝈 = 𝝀𝒕𝒓𝜺𝑰 + 𝟐𝝁𝜺 = 𝒕𝒓𝜺𝑰 + 𝜺
𝟏 + 𝝂 𝟏 + 𝟐𝝂 𝟏+𝝂
𝑰 : Matrice unité
(𝝀 , 𝝁): Coefficient de Lamé,
𝝂 Coefficient de Poisson
𝑬 : Coefficient de Young
Le problème plan à deux dimensions est déduit à partir du problème à trois dimensions
donné par [15].
Toutes les variables sont indépendantes du troisième axe et ne sont fonction que des
deux autres coordonnées 𝒙, 𝒚; on parle alors d’état de contraintes planes ou de
déformations planes.
Page 13
CHAPITRE II
Dans les deux cas de figure, états de contraintes planes [15] et de déformations planes,
nous pouvons tirer profit des situations les caractérisant pour rendre les calculs plus
faciles.
Etat de déformations planes:
Le modèle s’applique à un corps de section constante quelconque, infiniment long,

conventionnellement d’axe noté 𝒛, encastré à ses deux extrémités. Le corps est soumis
à un chargement orthogonal à 𝒛, uniforme sur toute sa longueur. Dans ces conditions,
aucune dilatation suivant l’axe 𝒛 n’est possible : chaque tranche ne se déforme que
dans son plan. Dans ce cas, le déplacement dans la direction 𝒛 peut être considéré
comme nul. Les dérivées des déplacements dans cette même direction sont nulles :
𝜺𝒛𝒛 = 𝜺𝒙𝒛 = 𝜺𝒚𝒛 = 𝜺𝒛𝒙 = 𝜺𝒛𝒚 = 𝟎
Par contre il existe des contraintes 𝝈𝒛𝒛 qui contrebalancent l’effet de Poisson c'est-à-
dire la tendance à la dilatation suivant l’axe 𝒛. la contrainte qui retient chaque tranche
dans le plan vaut:
𝝈𝒛𝒛 = 𝝂(𝝈𝒙𝒙 + 𝝈𝒚𝒚 )
Etat de contraintes planes:
Une structure plane mince est en état de contraintes planes si les charges sont
appliquées dans son plan. L’épaisseur de la structure est toujours très petite par rapport
aux deux autres dimensions et est symétrique de part et d’autre du plan (𝒙, 𝒚), appelée
surface neutre.
Les composantes du vecteur contraintes 𝝈 associé au tenseur des contraintes de

Cauchy 𝝈𝒄 sont :
𝑻 (2-6)
𝝈 = 𝝈𝒙𝒙 , 𝝈𝒚𝒚 , 𝝈𝒛𝒛 , 𝝈𝒙𝒚, 𝝈𝒚𝒛, 𝝈𝒛𝒙,
Page 14
CHAPITRE II
Ainsi les contraintes 𝝈𝒛𝒛 = 𝝈𝒙𝒛 = 𝝈𝒚𝒛 = 𝝈𝒛𝒙 = 𝝈𝒛𝒚 deviennent nulles. La déformation
𝜺𝒛𝒛 n’est pas nulle par effet de poisson. Elle peut être calculée en utilisant l’hypothèse
𝝈𝒛𝒛 = 𝟎
𝑬
𝜺𝒛𝒛 = 𝟏−𝝂𝟐 𝜺𝒙𝒙 − 𝜺𝒚𝒚 (2-7)
La loi de Hooke reliant le vecteur de contraintes au vecteur de déformations est:
𝛔 = 𝐃 𝛆 (2-8)
Ces deux approches sont basées sur des hypothèses fondamentalement différentes et
ne s’appliquent donc pas aux mêmes structures. Elles sont déterminées à partir des
équations (2-5a) et (2-5b).
II.2.2. Conditions aux limites:
La résolution des équations régissant le comportement d’une structure nécessite

l’introduction de conditions aux limites de la structure [16].
Pour un solide occupant un domaine Ω𝑺 dans ℝ𝟑 de frontière𝝏Ω𝑺, on peut diviser la

surface 𝝏Ω 𝑺 en deux parties, tel que 𝝏Ω𝑺 = 𝝏Ω𝑺𝟎 𝑼 𝝏Ω𝑺𝝈 (Figure (I-1)).
Figure (II-2) : Schématisation

générale du problème structure
Le mouvement de la structure est décrit ainsi par un champ de déplacement 𝑼 = (𝑼𝒊 )

défini en tout point 𝒙 du domaine ΩS. Les déformations sont représentées par le
tenseur des déformations noté 𝜺 𝑼 = 𝜺𝒊𝒋 (𝑼) et les contraintes par le tenseur 𝝈𝒊𝒋 :
𝝈 𝑼 = 𝝈𝒊𝒋 (𝑼) .
Page 15
CHAPITRE II
Avec les hypothèses précédentes, la relation fondamentale de la dynamique pour le

milieu continu s’écrit :
𝝏𝝈𝒊,𝒋 (𝑼)
−𝝎𝟐 𝝆𝑺 𝑼𝒊 − =𝟎 𝑑𝑎𝑛𝑠 Ω𝑺 (2-9)
𝝏𝒙𝒊
Avec la relation liant déplacement/déformation et la loi de comportement

déformation/contraintes données respectivement par (2-3) et (2-4), l’équation locale
précédente est assortie des conditions aux limites suivantes :
𝒖𝒊 = 𝟎 𝒔𝒖𝒓 𝝏𝜴𝒔𝟎 (2-10)

𝝈𝒊𝒋 𝑼 𝒏𝑺𝒋 = 𝟎 𝒔𝒖𝒓 𝝏𝜴𝒔𝝈
(2-11)
Qui traduisent le fait que la structure est bloquée sur 𝝏Ω𝑺𝟎 et libre d’effort sur 𝝏Ω𝑺𝝈 .
II-3. Théorie des poutres

II.3.1. Généralités :
Par définition une poutre est un solide engendré par une surface plane 𝑨 qui peut être
constante ou non la ligne moyenne est appelée axe neutre de la poutre [17].
Figure (II-3) : Déplacement

de la section droite d’abscisse
x
Compte tenu des hypothèses précédentes, le champ de déplacements s’écrit
𝒖 𝒙, 𝒚 = −𝒚𝜽𝒛(𝒙)
𝒗 𝒙, 𝒚 = 𝒗𝑮 𝒙 = 𝒗(𝒙) (2-12)
Page 16
CHAPITRE II
Les déformations :
𝝏𝒖 𝝏𝜽𝒛
𝜺𝒙𝒙 = = −𝒚
𝝏𝒙 𝝏𝒙 (2-13)
𝝏𝒖 𝝏𝒗 𝝏𝒗
𝜸𝒙𝒚 = + = −𝜽𝒛 +
𝝏𝒚 𝝏𝒙 𝝏𝒙
Les contraintes
𝝏𝜽𝒛
𝝈𝒙𝒙 = 𝑬𝜺𝒙𝒙 = −𝑬𝒚 (2-14)
𝝏𝒙
𝝏𝒗
𝝈𝒙𝒚 = 𝑮𝒌𝒚 𝜸𝒙𝒚 = 𝑮𝒌𝒚 − 𝜽𝒛 (2-15)
𝝏𝒙
 𝐴 est l’aire de la section droite,
 𝐺 = 𝐸 2(1 + 𝑣) est le module d’élasticité transversal,
 𝒌𝒚 est le coefficient d’aire cisaillée. Ce coefficient traduit le fait que le
cisaillement n’est pas uniforme dans la section, Il est défini par
𝑻𝟐
𝒚
= ∫𝑨 𝝈𝟐𝒙𝒚 𝒅𝑨 (2-16)
𝑨𝒌𝒚
Avec 𝑻𝒚 représente l’effort tranchant déduit de la théorie d’élasticité.
La contrainte normale dans la section droite est égale à
𝒚
𝝈𝒙𝒙 = − 𝑴𝒇𝒛 (2-17)
𝒍𝒛
II.3.2. Equations d’équilibre
La poutre porte une force et un couple répartis d’intensité linéique respectivement 𝑷𝒚
et 𝒎𝒄 .
L’équilibre d’un morceau de poutre compris entre les sections droites d’abscisses 𝒙 et
𝒙 + 𝒅𝒙 est régi par les équations d’équilibre données par les expressions (2-18) et
(2-19)
𝝏𝑻𝒚
+ 𝒑𝒚 = 𝝆𝑨𝒗 (2-18)
𝝏𝒙
𝝏𝑴𝑭𝒚 (2-19)
+ 𝑻𝒚 +𝒎𝒛 = 𝝆𝑰𝒛 𝜽𝒛
𝝏𝒙
Page 17
CHAPITRE II
Ces équations peuvent s’écrire en fonction des déplacements comme suit
𝝏 𝝏𝒗
𝑮𝑨𝒌𝒚 − 𝜽𝒛 + 𝒑𝒚 = 𝝆𝑨𝒗 (2-20)
𝝏𝒙 𝝏𝒙
𝝏 𝝏𝜽𝒛 𝝏𝒗 (2-21)
𝑬𝑰𝒛 + 𝑮𝑨𝒌𝒚 − 𝜽𝒛 +𝒎𝒛 = 𝝆𝑰𝒛 𝜽𝒛
𝝏𝒙 𝝏𝒙 𝝏𝒙
Ce modèle est dit modèle de Timoshenko. Si la poutre est longue, on admet

l’hypothèse de Bernoulli: la section droite reste normale à la fibre moyenne après
déformation (Figure (II-3)), d’où la relation cinématique
𝝏𝒗 (2-22)
𝜽𝒛 =
𝝏𝒙
Figure (II-4) : Hypothèse de

Bernoulli
II.4. Mécaniques des fluides
Les équations générales de la dynamique des fluides sont présentées dans ce chapitre.
Ces équations résultent des lois de conservation de la physique qui sont : conservation
de la masse, conservation de la quantité de mouvement (seconde loi de Newton) et
conservation de l’énergie (premier principe de la thermodynamique).
Conservation de la masse
Page 18
CHAPITRE II
Le principe de la conservation de la masse fait partie du fondement de la mécanique

classique. Il exprime la constance de la masse d’un système matériel, suivi dans son
mouvement.
La masse d’un système est donnée par

(2-23)
𝑴= 𝝆𝒇 𝒅𝜴
𝜴
Selon le principe énoncé, on a
𝑫𝑴 (2-24)
=𝟎
𝑫𝒕
En appliquant le théorème de transport, il résulte que
𝑫𝑴 𝝏
= 𝝏𝒕 ∫𝛀 𝝆𝒇 𝒅𝛀 + ∫𝚪 𝝆𝒇 𝒖 𝒏 𝒅𝚪 (2-25)
𝑫𝒕
Où 𝒏 est la normale extérieure unitaire à la frontière du domaine Γ et 𝑢 est la vitesse

du fluide.
En appliquant le théorème de la divergence, on écrit
𝝏𝝆𝒇 (2-26)
+ 𝛁. 𝝆𝒇 𝑼 𝒅𝛀 = 𝟎
𝛀 𝝏𝒕
L’intégrant doit être nul puisque le volume 𝛀 est arbitraire. On à donc
𝝏𝝆𝒇
+ 𝛁. 𝝆𝒇 𝑼 = 𝟎 (2-27)
𝝏𝒕
Conservation de la quantité de mouvement

Le principe de la quantité de mouvement traduit la seconde loi de Newton, appliqué à
un système matériel : pour tout système matériel, le taux de variation par rapport au
temps du tenseur des quantités de mouvement est égal au tenseur de forces extérieures
appliquées sur ce système.
D’après la mécanique classique, le tenseur des quantités de mouvement relatif à un tel

système, est exprimé par
(2-28)
𝑰= 𝝆𝒇 𝑼𝒅𝛀
𝛀
Page 19
CHAPITRE II
Les forces extérieures agissant sur le volume Ω sont les forces de volumes et les forces
de surface
(2-29)
𝑭= 𝝆𝒇 𝒇𝒅𝛀 + −𝒑 𝑰 + 𝝈 . 𝒏𝒅𝚪
𝛀 𝚪
Le principe énoncé exprime que

𝑫𝑰
= 𝑭 (2-30)
𝑫𝒕
En appliquant le théorème de transport, il résulte que

𝝏
𝝏𝒕 𝛀
∫ 𝝆𝒇 𝑼𝒅𝛀 + ∫𝚪 𝝆𝒇 𝒖(𝑼 𝒏)𝒅𝚪 = ∫𝛀 𝝆𝒇 𝒇𝒅𝛀 + ∫𝚪 −𝒑 𝑰. 𝒏𝒅𝚪 +
∫𝚪 𝝈. 𝒏𝒅𝚪 (2-31)
𝒑 est la pression, 𝑰 est la matrice identité, 𝝈 est le tenseur des contraintes visqueuses.
En appliquant le théorème de la divergence, on obtient
𝝏𝝆𝒇 𝑼 (2-32)
+ 𝛁. 𝝆𝒇 𝑼 × 𝑼 + 𝛁𝒑 + 𝛁. 𝝈 − 𝝆𝒇 𝒇 𝒅𝛀 = 𝟎
𝛀 𝝏𝒕
L’intégrant doit être nul puisque le volume Ω est arbitraire. On a donc
𝝏𝝆𝒇 𝑼
+ 𝛁. 𝝆𝒇 𝑼 × 𝑼 + 𝛁𝒑 + 𝛁. 𝝈 − 𝝆𝒇 𝒇 (2.33)
𝝏𝒕
Le tenseur des contraintes visqueuses est donné par

𝒕 (2-34)
𝝈 = 𝝁 𝛁𝑼 + 𝛁𝑼
En introduisant l’opérateur rotationnel l’équation de conservation de la quantité de

mouvement s’écrit
𝝏𝑼 𝟏
𝝆𝒇 + 𝟐 𝒈𝒓𝒂𝒅𝑼𝟐 + 𝒓𝒐𝒕𝑼 ∧ 𝑼 = 𝒇 + 𝒅𝒊𝒗𝝈𝑪 (2-35)
𝝏𝒕
Le tenseur des contraintes s’exprime en fonction de la pression p, de la viscosité 𝝁 et

du tenseur gradient des taux de déformations.
𝝈𝑪 = −𝒑𝑰 + 𝟐𝝁𝑻𝑫 (2-36)
Page 20
CHAPITRE II
Dans le cadre de notre étude nous prenons les hypothèses simplificatrices suivantes :
L’écoulement considéré est un fluide parfait, i rotationnel et incompressible
 Fluide parfait
Dans un fluide parfait la viscosité est nulle, donc l’équation (2-36) devient (2-37)
𝝈𝑪 = −𝒑𝑰 (2-37)
 Mouvement irrotationnel
Le champ de vitesse de l’écoulement est irrotationnel quand
(2-38)
𝒓𝒐𝒕𝑼 = 𝟎
Cette équation donne la condition nécessaire et suffisante pour que le potentiel de
vitesse 𝚽 puisse exister : 𝑼 = 𝒈𝒓𝒂𝒅𝚽 , et par conséquent un écoulement

irrotationnel est appelé écoulement potentiel.
Et dans ces conditions l’équation d’Euler devient
𝝏𝑼 𝟏
𝝆𝒇 + 𝟐 𝒈𝒓𝒂𝒅𝑼𝟐 = 𝒇 + 𝒅𝒊𝒗𝝈𝑪 = 𝒇 + 𝒈𝒓𝒂𝒅𝒑 (𝟐-39)
𝝏𝒕
 Les forces volumiques dérivent d’un potentiel
L’action du champ de pesanteur conduit à :
𝒇 = −𝝆𝒇 𝒈𝒌 = −𝒈𝒓𝒂𝒅(𝝆𝒇 𝒈𝒛) (2-40)

Et après intégration, l’équation d’Euler dans le domaine fluide donne l’équation de
Bernoulli
𝝏𝚽 𝟏 𝒑
+ 𝟐 𝑼𝟐 + 𝒈𝒛 + 𝝆 = 𝑪(𝒕) (2-41)
𝝏𝒕 𝒇
Les hypothèses de petits mouvements du fluide
Donc en négligeant 𝒖𝟐 devant les autres termes, l’équation de Bernoulli conduit à
L’équation de Lagrange linéaire qui en choisissant 𝑪(𝒕) = 𝟎 est
Page 21
CHAPITRE II
𝝏𝚽
𝒑 = 𝝆𝒇 (2-41)
𝝏𝒕
 Hypothèse de compressibilité
 Si le fluide est incompressible, implique une vitesse de propagation des
perturbations de pression faibles dans l’écoulement (ondes acoustiques) d’où les
équations de l’écoulement potentiel décrit par les variables 𝑼, 𝝓 𝑒𝑡 𝒑 et l’équation de

continuité se transforme en équation de Laplace
(2-42)
𝚫𝚽 = 𝟎 𝒅𝒂𝒏𝒔 𝛀𝒇
L’écoulement étudié dans ce travail est supposé évoluer de manière iso-volumique.
Cela signifie que l’on peut négliger les effets de dilatation thermique dus aux
variations de température, ainsi que les variations de volume liées à des compressions
détentes associées à des changements de pression.
Pour simplifier, les champs de températures seront supposés uniformes et constants.
L’échauffement résultant des frottements étant négligeable, cela impose à la
température des frontières d’être uniforme et indépendante du temps.
Par suite, puisque ces écoulement ne sont le siège d’aucun frottement ni échange de
chaleur, ils sont isentropiques.
Les effets de variations des pressions sur la masse volumique peuvent alors etre
négligés.
L’incompressibilité impose 𝒅𝒊𝒗𝑼 = 𝟎 (2-43)
Avec la conservation de masse cela conduit à
(2-44)
𝝏𝝆𝒇
=𝟎
𝝏𝒕
Qui signifie que la masse volumique d’une particule fluide suivie pendant une période
de temps dt reste constante.
Bien que cette contrainte n’interdise pas à 𝝆 de varier d’une particule à une autre, on
supposera que l’écoulement provient d’une région de l’espace constituant soit un
réservoir, soit un domaine infini à l’amont, où la masse volumique est uniforme.
Ainsi 𝝆 aura la même valeur pour toutes les particules fluides, et donc sera donnée,
uniforme et indépendante du temps.
Page 22
CHAPITRE II
Le principe fondamental de la dynamique est alors traduit par l’équation (2-38) dans
laquelle on néglige les termes de frottement, et qui simplifie en l’équation d’EULER :
𝒅𝑼 (2-45)
𝝆𝒇 = 𝝆𝒇 𝒈 − 𝒈𝒓𝒂𝒅𝒑
𝒅𝒕
Elle correspond à une loi de comportement réduite à
𝝈𝒊𝒋 𝒙, 𝒕 = −𝒑(𝒙, 𝒕)𝜹𝒊𝒋 (2-46)
Qui exprime la contrainte
𝑻𝒊 𝒙, 𝒕, 𝒏 = −𝒑 𝒙, 𝒕 𝜹𝒊𝒋𝒏𝒋 = −𝒑 𝒙, 𝒕 𝒏𝒊 (2-47)
Qui est normale à cet élément, par conséquent les efforts tangentiels de contact
(frottements) sont exclus
La masse volumique étant connue on dispose alors des équations suivantes
𝒅𝑼 (2-48)
𝝆𝒇 = 𝝆𝒇 𝒈 − 𝒈𝒓𝒂𝒅𝒑
𝒅𝒕
C'est-à-dire d’un système de quatre équations indépendantes, pour les quatre

inconnues que sont les trois composantes 𝑼 𝒙, 𝒕 du champ de vitesse et le champ de
pression𝒑 𝒙, 𝒕 .
Pour exploiter ce système, il est a priori nécessaire de lui associer les conditions aux
limites suivantes
Condition de surface libre sur 𝚪𝐋
Compte tenu des hypothèses retenues pour l’étude du fluide donnera
𝝏𝟐 𝚽 𝝏𝚽
 + 𝒈 𝝏𝒛 = 𝟎 sur le plan 𝒛 = 𝟎 (surface libre au repos), lorsque les
𝝏𝒕𝟐
mouvements de surface libre sont de faibles amplitudes, et que le fluide
Page 23
CHAPITRE II
incompressible est étudié dans la gamme des basses fréquences des

ondes de gravité telles que les fréquences de houle)
 𝚽 = 𝟎 sur le plan 𝒛 = 𝟎 dans le cas d’un fluide incompressible en
hautes fréquences ou bien dans le cas d’un fluide compressible.
𝝏𝚽
 𝚽= = 𝟎 sur le plan 𝒛 = 𝟎 dans le cas d’un fluide incompressible
𝝏𝒏
en très basses fréquence. Dans ce cas, la surface libre est équivalente à
une paroi rigide.
Conditions à l’infini :
Dans le cas d’une étude en milieu non borné, la condition à l’infini affichée est la
condition de Sommerfeld (non réflexion des ondes à l’infini)
𝝏𝚽 𝟏 𝝏𝚽
𝐥𝐢𝐦𝒓→∞ 𝒓𝜶 ( + ) = 𝟎 𝜶 = 𝟎. 𝟓 𝒆𝒏 𝟐𝑫 𝒆𝒕 𝜶 = 𝟏 𝒆𝒏 𝟑𝑫 (2-49)
𝝏𝒓 𝒄 𝝏𝒕
Et :
𝐥𝐢𝐦𝒓→∞ 𝚽 = 𝚽𝟏 𝚽𝟏 potentiel incident, dû à la houle par exemple ou nul

dans le cas où il n’y a pas d’écoulement.
Condition de glissement sur les parois du solide
La condition de Dirichlet imposée, implique que les particules fluide, ne pénètrent pas
les parois de solides, c'est-à-dire qu’il y’a continuité des vitesses normales des
particules fluides et solides
𝝏𝚽
𝒗𝒔 ∙ 𝒏 = = 𝒈𝒓𝒂𝒅𝚽 ∙ 𝒏 𝒔𝒖𝒓 𝚪𝐒𝐅 (2-50)
𝝏𝒏
Et 𝒗𝒔 ∙ 𝒏 = 𝑼 ∙ 𝒏
Où 𝒏 est la normale extérieure au domaine structure, donc intérieure au domaine

fluide, cette condition peut également être écrite en termes de pression
𝝏𝒑 (2-51)
−𝝆𝒇 𝒚𝒔 ∙ 𝒏 = 𝒔𝒖𝒓 𝚪𝐒𝐅
𝝏𝒏
II.5. Couplage
L’interaction structure-fluide exprime le problème mécanique du contact d’un corps

solide et d’un fluide. Cette étape essentielle de notre problème nécessite la
Page 24
CHAPITRE II
caractérisation de conditions aux limites qui décrivent le dialogue entre le solide et le

fluide au niveau de l’interface.
Figure (II-5) :
Schématisation générale du
problème structure-fluide
A l’interface, figure (2-5) le lien entre vecteurs normaux est défini de telle manière
que 𝒏 = 𝒏𝒔 = −𝒏𝒇 (2-52)
On rappelle que le problème étudié est celui d’une structure élastique couplée avec un
fluide incompressible, parfait et i rotationnel, la situation étant décrite par la figure
(2-16).
La structure occupe le domaine Ω𝑺 , de frontière 𝝏Ω𝑺 ; le fluide occupe le
domaineΩ𝑭 , de frontière 𝝏Ω𝑭 . les deux domaines sont en contact sur 𝜞𝑺𝑭 =
𝝏Ω𝑺 ∩ 𝝏Ω𝑭. Les éléments de frontière 𝝏Ω𝑺𝟎 , 𝝏Ω𝑺𝝈 , 𝝏Ω𝑭𝟎 , et 𝝏Ω𝑭 désignent
respectivement, les portions de frontière du domaine de la structure à déplacement
imposé, à effort imposé, et les portions de frontière du domaine fluide à pression et
gradient de pression imposé.
Figure (II-6) :
Schématisation générale du
problème structure-fluide
Page 25
CHAPITRE II
Le couplage des deux milieux est modélisé par les équations qui expriment la
continuité des composantes normales du tenseur des contraintes et du champ
d’accélération à l’interface structure-fluide.
La condition équation (2-47), peut également être écrite en termes de pression
𝝈𝒊𝒋 𝑼 𝒏𝒔𝒋 = 𝒑𝒏𝒊 sur 𝜞𝑺𝑭 (2-53)
𝝏𝒑
𝒏𝒋 = 𝝎𝟐 𝝆𝑭 𝑼𝒋 𝒏𝒋 sur 𝜞𝑺𝑭 (2-54)
𝝏𝒙𝒊
La relation (2-53) traduit le fait que le fluide exerce des efforts de pression sur la
structure au milieu de l’interface structure-fluide, dans la direction normale sur
l’interface.
La relation (2-54) traduit le fait que la structure impose une accélération au fluide au
niveau de l’interface structure-fluide dans la direction normale à l’interface.
Le problème consiste à la recherche de fonctions scalaires satisfaisant l’équation de

Laplace pour les hypothèses citées ci-dessus lorsque le fluide incompressible est limité
par une surface libre déformable (onde de gravité).
Page 26
Page27
MODELISATION PAR ELEMENTS FINIS ET COUPLAGE
CHAPITRE III
III.1. Introduction
Les phénomènes liés au couplage fluide-structure apparaissent à des degrés dive

d’importante, pour toute structure en contact avec un fluide, on entend par problème de
couplage fluide-structure mécanique tout problème dans lequel ;
 Les milieux solides et fluides échangent des énergies d’origine mécanique.

 Les phases solides et fluides sont distinctes.
 Les milieux solides et fluides n’échangent pas de masse.
Le solide et le fluide occupant des domaines de l’espace note

respectivement𝛀𝐒 et𝛀𝐅 sont couplés dans le sens suivant au niveau de la zone de contact
entre le fluide et la structure (ou interface fluide-structure𝚪), des échanges d’énergie

mécanique se produisent dans les deux sens.[23 ; 43].
 Le fluide exerce des efforts sur la structure en mouvementce qui modifie la

dynamique de la structure.
 La structure impose des déplacements sur le fluide et modifie ainsi les
caractéristiques de l’écoulement fluide.
Forces
Structure
Fluide Déplacement, déformation,
Vitesse, pression et contrainte
Déplacements
Interface fluide-structure
Figure: (III-1) Schématisation générale d’un problème de couplage mécanique
fluide-structure
III.2. Méthode des éléments finis pour l’analyse vibratoire d’une

structure couplée en présence d’un fluide incompressible
On s’intéresse à la description des vibrations linéaires d’une structure en présence d’un

fluide incompressible.
Page 28
CHAPITRE III
On commence par rappeler les équations de l’élasticité linéaire décrivant les petits
mouvements d’une structure autour d’une position de référence.
les notions de base dynamique des structuressont proposés dans les ouvrages de
références [24,25, 26].
On présente ensuite la méthode des éléments finis utilisée pour résoudre le problème
structure dans les ouvrages des références [30,28, 29].
On s’intéresse à la description des vibrations linéaires d’une structure couplée avec un
fluide incompressible.
On commence par rappeler les équations décrivant le problème couplé, les termes de
couplage traduisant l’action du fluide sur la structure (effort de pression imposé au
niveau de l’interface fluide-structure) et l’action de la structure sur le fluide
(accélération normale imposée au niveau de l’interface fluide-structure).
On expose ensuite la méthode de résolution éléments finis du problème, utilisant un

couplage entre la discrétisation éléments finis de la structure.
La discrétisation de termes de couplage conduit à la définition d’un opérateur de

couplage fluide-structure, pris d’un opérateur de masse ajoutée.
 Contexte d’étude et hypothèses générales
Le problème étudié est celui d’une structure élastique linéaire couplée avec un fluide
incompressible, la situation étant décrite par la figure (III-2)
𝛛𝛀𝐅𝛑
𝛛𝛀𝐒𝟎
𝛀𝐒
𝚪 (𝑬, 𝑽, 𝝆𝑺 ) 𝐧𝐒
𝛀𝐅 𝛈 𝐮
(𝑪, 𝝆𝑭 ) 𝛛𝛀𝐒𝛔
𝐏
Figure: (III-2) Schématisation générale du problème fluide-
structure.
𝛛𝛀𝐅𝟎 𝐮
𝐧𝐅
𝚪
Page 29
𝛈
CHAPITRE III
La structure occupe le domaine𝛀𝐒 , de frontière𝝏𝛀𝐒
Le fluide occupe le domaine𝛀𝐅 ,de frontière𝝏𝛀𝐅.
Les deux domaines sont en contact sur𝚪 = 𝝏𝛀𝐒 ⋂𝝏𝛀𝐅 .
Les éléments de frontière𝛛𝛀𝐒𝟎 , 𝛛𝛀𝐒𝛔 , 𝛛𝛀𝐅𝟎 et 𝛛𝛀𝐅𝛑 désignent respectivement
comme dans les sections précédentes.
III-3 Les éléments Finis en mécanique du solide
La relation numérique des problèmes posés par la formulation utilisé de façon classique
Formulation général
Nous rappelons seulement dans ce paragraphe les étapes de la discrétisation par élément
finis de la forme variationnelle « principe des travaux virtuels » en élasto-
dynamique[27, 29, 32].
Le domaineΩ𝒔 est divisé en sous-domaines élémentairesΩ𝒆 de formes simples sur
lesquels il est aisé de définir, à partir de l’approximation polynomiale généralisée, une
approximation nodale (interpolation) de la forme :
{𝒖} = [𝑵]{𝒂}𝒆 (3-1)
Ou
 {𝒖} : Vecteur dont les composantes sont celles qui caractérisent le
vecteurdes déplacementsu ⃗en𝑴(𝒙)deΩ𝒆
 {𝒂}𝒆 : Vecteur des déplacements nodaux du domaine élémentaireΩ𝒆
 [𝑵] : matrice des fonctions d’interpolation
Le calcul du tenseur des petites déformationsE permet d’exprimer les composantes du
vecteur des déformations{𝝃}en fonction des déplacements nodaux sous la forme :
{𝝃} = [𝑩] {𝒂}𝒆 (3-2)
La loi de Hooke généralisée, en remplaçant {𝝃}par sa valeur s’écrit
{𝝈} = [𝑫] [𝑩] {𝒂}𝒆 (3-3)
Page 30
CHAPITRE III
La forme quadratique associée à l’énergie de déformation sur l’élément fini est
(3-4)
𝟐𝑬𝒅𝒆 = ∫Ω 𝒆 {𝝃}’ {𝝈}𝒅𝑽
Soit donc en remplaçant{𝝃}et{𝝈}par leurs valeurs précédentes :
𝟐𝑬𝒅𝒆 = < 𝒂 > ∫ Ω𝒆

𝑩 ’ 𝑫 𝑩 𝒅𝑽 𝒂 𝒆 (3-5)
On en déduit la forme générale de la matrice raideur d’un élément finiΩ𝒆 dans la base
des déplacements nodaux {𝒂}𝒆
[𝑲]𝒆 = ∫Ω 𝒆 [𝑩]𝒕 [𝑫][𝑩]𝒅𝑽 (3-5)
De même la forme quadratique associé a l’énergie cinétique de l’élément permet de

définir la matrice masse de l’élément finiΩ𝒆 dans la basse des déplacements
nodaux{𝒂}𝒆 . telle que
𝟐𝑬𝒄𝒆 = ∫ 𝒖𝟐 𝒑𝒅𝑽 = < 𝒂 >𝒆 ∫Ω 𝒆 [𝑵]𝒕 𝝆[𝑵]𝒅𝑽 {𝒂}𝒆 (3-6)

Ω𝒆
(3-7)
𝒕
[𝒎]𝒆 = ∫Ω 𝒆 [𝑵] 𝝆[𝑵]𝒅𝑽
Et enfin le calcul du travail virtuel des forces extérieures conduit à la définition de

forces généralisées qui peuvent être de différentes natures
 Forces généralisées associées aux répartitions volumiques de forces {𝒒𝒗 }
𝜹𝑻𝒗𝒆 = ∫Ω 𝒆 < 𝛿𝒂 >𝒆 {𝒒𝒗 }𝒅𝑽 = < 𝜹𝒂 >𝒆 {𝝋𝒗 }𝒆 (3-7)
Soit
{𝝋𝒗 }𝒆 = ∫Ω𝒆 [𝑵]𝒕 {𝒒𝒗 }𝒅𝑽 (3-8)
 Forces généralisées associées aux répartitions surfaciques de forces{𝒒𝒔}
∫𝝏Ω 𝒆 < 𝛿𝒂 > 𝒆 {𝒒𝒔 }𝒅𝑽 = < 𝜹𝒂 > 𝒆 {𝝋𝒔 }𝒆 (3-9)
Page 31
CHAPITRE III
Soit
(3-10)
{𝝋𝒔}𝒆 = ∫Ω 𝒆 [𝑵]𝒕 {𝒒𝒔 }𝒅𝑽
 Forces généralisées associées aux forces appliquées aux nœuds de l’élément

(efforts de liaisons)
𝜹𝑻𝒔𝒆 =< 𝜹𝒂 >𝒆 {𝝋𝒍 }𝒆 (3-11)
 Forces généralisées associées à une dissipation d’énergie. Soit si {𝑭𝒅 }est la «

force élémentaire dissipative », telle que
{𝑭𝒅 } = −𝒃{𝒖𝒑 } (3-12)
Le travail virtuel de cette force s’écrit
𝜹𝑻𝒅𝒆 = −∫Ω 𝒆 < 𝛿𝒂 > {𝑭𝒅 }𝒅𝑽 = < 𝜹𝒂 >𝒆 {𝝋𝒅 }𝒆 (3-13)
{𝝋𝒅 }𝒆 = − ∫Ω 𝒆 [𝑵]𝒕 𝒃[𝑵]𝒅𝑽 {𝒂}𝒆 (3-14)
Ce qui conduit à la définition d’une matrice amortissement[𝒃]𝒆 telle que
[𝒃]𝒆 = ∫Ω 𝒆 [𝑵]𝒕𝒃[𝑵]𝒅𝑽 (3-15)
En prenant les variations par rapport à chaque variable généralisée (écriture des
équations de Lagrange) on obtient l’équation matricielle caractérisant le mouvement de
l’élément fini
[𝒎]𝒆 {𝒂}𝒆 + [𝒃]𝒆 {𝒂}𝒆 + [𝑲𝒔 ]{𝒂}𝒆 = {𝝋𝟏 }𝒆 + {𝝋𝒗 }𝒆 + {𝝋𝒔 }𝒆 (3-16)
L’équation matricielle de la structure assemblée est obtenue en additionnant les formes

intégrales élémentaires et en prenant les variations des quantités additionnées est
[𝑴𝒔 ]{𝒂} + [𝑩𝒔 ]{𝒂} + [𝑲𝒔 ]{𝒂} = {𝑭𝒆𝒙𝒕 } + {𝝋𝒗 }𝒂𝒔𝒔 + {𝝋𝒔 }𝒂𝒔𝒔 (3-17)
Avec{𝑭𝒆𝒙𝒕 }vecteur des forces extérieures directement appliquées aux nœuds de la

structure.
Page 32
CHAPITRE III
III.4. Méthode des éléments finis en mécanique des

fluides.(Forme faible de la méthode des résidus pondérés)
Les équations des problèmes sont obtenues par linéarisation des équations générales
d’écoulement du fluide [23,34, 43].
Les équationslinéaires des fluide peuvent être décrites en utilisant une formulation en
déplacement 𝝁𝑭 [36,37]Le potentiel des vitesses[38]en potentiel des déplacements
𝝋[39,40]ou en pression 𝒑[41,42].
Nous considérons des fluides avec les hypothèses précédentes.
Ils sont caractérisés par les équations caractéristiques [43]
𝟏 𝝏𝟐 𝝋
∆𝝋 − = 𝟎 𝒅𝒂𝒏𝒔 𝛀𝑭
𝑪𝟐 𝝏𝒕𝟐
𝝋 = 𝟎, 𝒐𝒖
𝟐
𝝏𝝋 𝟏 𝝏 𝝋 𝒆𝒏 𝒊𝒏𝒄𝒐𝒎𝒑𝒓𝒆𝒔𝒔𝒊𝒃𝒍𝒆
+ 𝟐
=𝟎
𝝏𝒛 𝒈 𝝏𝒕
𝝏𝝋
= 𝑽⃗𝑺 ∙ 𝒏⃗ 𝒔𝒖𝒓 𝒍𝒆𝒔 𝒊𝒏𝒕𝒆𝒓𝒇𝒂𝒄𝒆𝒔 𝒇𝒍𝒖𝒊𝒅𝒆 − 𝒔𝒐𝒍𝒊𝒅𝒆 𝚪𝐒𝐅
𝝏𝒛
Sur la surface libre𝚪𝐋 (plan𝒛 = 𝟎)
Dans un premier temps, on considère que le domaine𝛀𝑭 est borné.
Il est divisé en sous-domaines élémentaires𝛀𝒆 sur lesquels le potentiel des vitesses𝝋est
interpolé en fonction des potentiels nodaux 𝚽 𝒆
𝚽= 𝑵 𝚽 𝒆
Pour simplifier la mise en œuvre de la formulation éléments finis, cette relation est
écrite en considérant que l’approximation nodale est réalisée sur le domaine𝛀𝑭 complet
sous la forme
 𝚽 vecteur des potentiels aux nœuds du domaine fluide complet𝛀𝑭
 𝑵⃗𝒊 qui s’identifie à la fonction d’interpolation𝑵𝒊 calculée au nœud𝒊d’un
élément.
Page 33
CHAPITRE III
En remplaçant 𝚽par son approximation dans l’équation locale précédente, on obtient un
résidu que l’on pondère par des fonctions de pondération𝜳. Considérons la forme
intégrale𝑾 𝚽 = 𝑾𝟏 𝚽 + 𝑾𝟐 𝚽
Avec
𝟏 𝝏𝟐 𝝋
𝑾𝟏 𝚽 = ∫𝜴𝑭 𝜳∆𝚽𝒅𝑽 𝒆𝒕 𝑾𝟐 𝚽 = − ∫𝜴𝑭 𝜳 𝟐 𝟐 𝒅𝑽
𝑪 𝝏𝒕
En appliquant la formule deGreen(intégration par parties généralisée) à la forme
intégrale𝑾𝟏 𝝋 onobtient
∫𝜴𝑭 𝜳∆𝚽𝒅𝑽 + ∫𝜴𝑭 𝒈𝒓𝒂𝒅⃗𝜳 ∙ 𝒈𝒓𝒂𝒅⃗𝚽𝒅𝑽 = ∫𝝏𝜴𝑭 𝜳𝒈𝒓𝒂𝒅⃗𝚽 ∙ 𝒏⃗ 𝒅𝑰

Donc la forme intégrale s’écrit
𝝏𝟐 𝝋
𝑾 𝚽 = −∫𝜴𝑭 𝜳 𝟐 𝒅𝑽 − +∫𝜴𝑭 𝒈𝒓𝒂𝒅⃗𝜳 ∙ 𝒈𝒓𝒂𝒅⃗𝚽𝒅𝑽
𝝏𝒕
+ ∫𝝏𝜴𝑭 𝜳𝒈𝒓𝒂𝒅⃗𝚽𝝋 𝒏⃗𝒅𝑰
La formulation faible du problème fluide se traduit par
La recherche du potentiel, donc le vecteur des potentiels nodaux 𝚽 , qui annulera
forme intégrale𝑾 𝚽 , quelque soit les fonctions de pondération 𝜳cinématiquement
admissible (nulles sur la partie de frontière 𝝏𝜴𝑭 ou le potentiel 𝚽est connu) et tel que
𝝏𝝋
sur la partie de frontière complémentaires𝝏𝜴𝟐 , soit donné.
𝝏𝒏
Soit compte tenu de l’orientation de la normale choisie pour écrire la condition de
glissement, il vient
𝝏𝚽
𝒈𝒓𝒂𝒅⃗𝚽 ∙ 𝒏⃗ = = −𝑽⃗𝑺 𝒏⃗ 𝒔𝒖𝒓 𝚪𝑺𝑭 (𝚪𝑺𝑭 ⊂ 𝝏𝜴𝟐 )
𝝏𝒏
𝝏𝚽 𝟏 𝝏𝟐 𝚽
=− (𝒐𝒏𝒅𝒆 𝒅𝒆 𝒈𝒓𝒂𝒗𝒊𝒕é 𝒆𝒏 𝒊𝒏𝒄𝒐𝒎𝒑𝒓𝒆𝒔𝒔𝒊𝒃𝒍𝒆)
𝝏𝒏 𝒈 𝝏𝒕𝟐
𝝏𝚽
= 𝟎 (𝒕𝒓è𝒔 𝒃𝒂𝒔𝒔𝒆𝒔 𝒇𝒓é𝒒𝒖𝒆𝒏𝒄𝒆𝒔 𝒆𝒏 𝒊𝒏𝒄𝒐𝒎𝒑𝒓𝒆𝒔𝒔𝒊𝒃𝒍𝒆)
𝝏𝒏
Donc la condition de surface libre sera prise en compte par des conditions sur 𝝏𝜴𝟏 ou
bien par des conditions sur 𝝏𝜴𝟐 .
Page 34
CHAPITRE III
La principale difficulté dans l’analyse du comportement d’un fluide par élément fini en
domaine non borné est la prise en compte correcte des conditions à l’infini.
Il existe plusieurs manières de détourner cette difficulté
 Ramener la frontière à distance fini et afficher sur cette frontière des

conditions similaires le comportement du fluide à l’infini.
 Développer des éléments finis particuliers : les éléments infinis.
 Développer une formulation mixte élément fini-équations intégrales et
comme pour la condition de surface libre, suivant le type d’approximation,
les conditions à l’infini, seront prises en compte comme des conditions
sur𝝏𝜴𝟏 ou sur𝝏𝜴𝟐 .
Une méthode courante guidant le choix des fonctions de pondération est la méthode de
Galerkin(qui s’identifie pour les problèmes de mécanique du solide au principe des
travaux virtuels) ; elle consiste à développer les fonctions𝜳sur les mêmes fonctions de
base que celle utilisées pour l’interpolation de 𝚽
Dans notre cas on n’a pas tenu compte de la compressibilité du fluide
Dans ce cas la matrice d’amortissement est nulle.
III.5. Couplage éléments finis-éléments finis
La résolution du problème couplé peut être conduite avec une méthode éléments finis,
on utilise un couplage entre
 Des éléments finis poutre à deux nœuds, à deux degrés de liberté par nœud
 Des éléments finis fluide axisymétrique en série de Fourier à quatre
nœuds, à un degré de liberté par nœud permettant par maillage du pavé
𝑹. 𝑹′ × 𝟎. 𝑯 la résolution du problème fluide :
Le principe du maillage du problème dans le plan(𝒓, 𝒛)peut être
schématisé par la figure suivante
On utilise un maillage à𝑵éléments finis de la poutre et un maillage à 𝑴 × 𝑵𝟎 éléments
finis du fluide avec𝑵𝟎 𝑵 = 𝑯 𝑳 (3-17)
Page 35
CHAPITRE III
𝒛
Elément fini Elément fini
Structure Fluide
𝑯
𝑳
𝑵 + 𝟏 Nœuds 𝑵𝟎 + 𝟏
𝟎 𝒓
Nœuds
𝑹
𝑹′
𝑴 + 𝟏 Nœuds
Figure (III-3)
Principe de maillage du problème couplé dans le plan (𝒓, 𝒛)
Cette analyse du comportement de la structure et du fluide par la méthode des éléments

fini est d’une mise en œuvre classique et aisé pour traiter les problèmes intérieurs et elle
a fait l’objet de nombreux travaux et communications.
Nous prenons ici la formulation générale conduisant à la construction des deux

systèmes d’équations linéaires couplés, obtenus par discrétisation en éléments finis des
deux milieux. Afin de ne pas alourdir inutilement les écritures matricielles, nous
supposerons qu’à l’interface fluide-structure, les discrétisations en éléments finis des
deux milieux sont identiques. Nous ne ferons donc pas de différences entre la frontière
mouillée d’un élément fini structure(𝚪𝐒𝐅 )𝒆 etd’un élément fini fluide en contact avec la
structure. Cette frontière élémentaire commune sera notée simplementҐ𝒔𝒇𝒆 . Le cas ou
les maillages structure et fluide à l’interface sont différents ne présente pas de difficultés
théoriques particulières, seule la mise en œuvre pratique de la méthode est un peu plus
délicate.
Page 36
CHAPITRE III
L’approximation nodale par éléments finis sur chaque élément solide s’écrit
{𝒖} = [𝑵𝒔 ]{𝒂}𝒆 (3-18)
Elle conduit donc conformément au calcul classique des quantités élémentaires énergies
cinétique, de déformation et travaux virtuels. Pour le calcul des forces généralisées
associées à la répartition surfacique, il intervient deux types de répartitions :
 des répartitions connues d’origine mécanique

 des répartitions inconnues due au chargement du fluide : les pressions
dynamiques sur Ґ𝒔𝑭 la partie « mouillée » de𝝏Ω𝑺𝝈 .
Donc pour les éléments finis ayant une frontièreҐ𝒆 appartenant àҐ𝒔𝑭 le travail virtuel des
forces dues au fluide (forces de pression) s’écrit :
𝜹𝑻𝒑𝒉 = ∫𝚪𝐬𝐟𝐞 𝜹𝒖⃗(−𝒑)𝒏⃗𝒅𝑺 (3-19)
Et en tenant compte de l’interpolation nodale, la force généralisée associée a alors pour

expression
(3-20)
𝛗𝐡 𝐞 = ∫𝚪𝐅 𝐍𝐬 𝐧 −𝐩 𝐝𝐒
𝐞
De même l’interpolation du potentiel des vitesses en fonction des variables nodales

élémentaires s’écrit
𝚽 = 𝑵𝑭 𝚽 𝒆 (3-21)
Et l’équation de Lagrange permet d’exprimer la pression dynamique en tout point d’un

élément fini fluide
𝑷 = −𝝆𝒇 𝑵𝑭 𝚽 𝒆 (3-22)
Et en particulier àl’interface fluide –structure, la force généralisée équivalente à la

pression dynamique à pour expression
Page 37
CHAPITRE III
𝝋𝒉 𝒆 = 𝝆𝒇 ∫𝚪𝐬𝐟𝐞 𝑵𝑻𝒔 𝒏𝑵𝒇 𝒅𝑺 𝚽 𝒆

(3-23)
Que l’on peut écrire
(3-24)
𝝋𝒉 𝒆 = 𝝆𝒇 𝒔𝒇 𝒆 𝚽 = − 𝒔𝒇 𝒆 𝒑 𝒆
(3-25)
Avec[𝒔𝒇]𝒆 = ∫Ґ 𝑵𝑻𝒔 𝒏𝑵𝒇 𝒅𝑺
𝒔𝒇𝒆
Matrice de couplage structure-fluide élémentaire, et{𝒑}𝒆 vecteur des pressions nodales

élémentaires.
En procédant d’une manière analogue pour l’analyse du fluide, on peut exprimer la

vitesse normale et un point quelconque de l’interface appartenant à un élément fini par
𝑽𝒏 = 𝑽⃗𝒔 . 𝒏⃗ = {𝒏}𝒕 𝑵𝒔 Ґ𝒔𝒇𝒆 {𝒂}𝒆 (3-26)
Donc les termes(𝒈𝒊)𝒆 calculés sur𝒔𝒇𝒆 € Ґ𝑺𝑭 s’écrivent
{𝒈}𝒆 = −∫𝚪𝐬𝐟𝐞 𝑵𝑻𝒇 {𝒏}𝒕 𝑵𝒔 𝒅𝑺 {𝒂}𝒆 (3-27)
Soit
(3-28)
{𝒈}𝒆 = − [𝑭𝒔]𝒆 {𝒂}𝒆 = −[𝑺𝑭]𝒕𝒆 {𝒂}𝒆
Donc après assemblage des formes intégrales élémentaire discrétisées, pour le fluide et
pour la structure, nous obtenons les deux systèmes d’équations linéaires couplés
suivant
𝑴𝒔 𝒂 + 𝑩𝒔 𝒂 + 𝑲𝒔 𝒂 = 𝑭𝒆𝒙𝒕 + 𝝆𝒇 [𝑺𝑭] 𝝋
(3-29)
𝑴𝒇 𝝋 + 𝑩𝒇 𝝋 + 𝑯𝒇 𝝋 = 𝑮𝒆𝒙𝒕 − [𝑭𝑺] 𝒂
Page 38
CHAPITRE III
Avec 𝑭𝒆𝒙𝒕 vecteur des forces généralisées associées à des forces d’excitation de la
structure, volumique, surfacique ou ponctuelles connues et 𝑮𝒆𝒙𝒕 un second membre du

système matriciel caractérisant le comportement du fluide associé à des conditions de
glissement sur des parois ou la vitesse normale est connue.
La matrice de couplage fluide-structure assemblée[𝑭𝑺]ou[𝑺𝑭] = [𝑭𝑺]𝒕 ne fait

correspondre que les inconnues nodales du fluide et de la structure attachées à
l’interfaceҐ𝑺𝑭 .
Dans le cas général la résolution de ce système se fera en groupant les inconnues

nodales 𝒂 et 𝝋 dans un même vecteur des inconnues nodales 𝒚 tel que
[𝑴] 𝒚 + [𝑩] 𝒚 + [𝑯] 𝒚 = 𝑸 (3-30)
Avec
[𝑴𝒔 ] (3-31)
[𝑴] =
[𝑴𝒇 ]
[𝑩𝒔 ]−𝝆𝒇 [𝑺𝑭] (3-32)

[𝑩] =
[𝑭𝑺] [𝑩𝒇 ]
𝑲𝒔 (3-33)
𝑯 =
𝑯𝒇
{𝑭𝒆𝒙𝒕 } (3-34)
𝑸 =
{𝑮𝒆𝒙𝒕 }
La résolution de ce système est simplifiée si on considère que le fluide est

incompressible et en se plaçant dans une gamme de fréquence telle que l’on puisse
considérer que la condition de surface libre= 𝟎. Dans ces conditions l’équation
matricielle traduisant le comportement fluide s’écrit
𝑯𝒇 𝝋 = 𝑮𝒆𝒙𝒕 − 𝑭𝒔 𝒂 (3-35)
Et il est alors possible d’éliminer les variables nodales fluides du système précédent et
dans ces conditions le système devient
Page 39
CHAPITRE III
([𝑴𝒔 ] + [𝑴𝒂 ]) {𝒂} + [𝑩𝒔 ]{𝒂} + [𝑲𝒔 ]{𝒂} = {𝑭𝒈𝒕 } (3-36)
Avec
[𝑴𝒂 ] = 𝝆𝒇 [𝑭𝒔 ]𝒕 [𝑯𝒇 ]−𝟏 [𝑭𝒔 ]Matrice « masse ajoutée » (3-37)
Et
{𝑭𝒈𝒕 } = {𝑭𝒆𝒙𝒕 } − 𝝆𝒇 [𝑯𝒇 ]−𝟏 {𝑮𝒆𝒙𝒕 } (3-38)
Comme nous l’avons déjà signalé, si le milieu fluide n’est pas borné, il est nécessaire de
mettre en œuvre des procédures spéciales permettant de tenir compte de la frontière à
l’infini. Ces procédures étant plus ou moins bien adaptées aux types de problèmes à
résoudre. Par exemple pour un calcul de fréquence et modes propres d’une structure
immergée en milieu infini, on peut tout en gardant une bonne précision étudier la
structure en milieu borné.
Page 40
Page 41
CAS D’UNE POUTRE IMMERGEE DANS UN FLUIDE
CHAPITRE IV
IV.1. Etude des vibrations d’une poutre immergée dans un fluide

Pour traiter ce problème, nous mettons en œuvre les différentes méthodes de résolution qui
ont été décrites précédemment.
On considère une poutre verticale d’axe𝒁, de section circulaire de rayon 𝒓 = 𝒂,encastrée à

une extrémité, libre à l’autre, immergée dans les autres directions.
étudie les vibrations de flexion de poutre dans le plan(𝒓, 𝒛), on suppose que𝒉 ≤ 𝒍 .
𝒁 𝒁
𝒗 𝒛, 𝒕 Vue de face
𝒉 𝒍
Vue de dessus
Figure IV.1 poutre vertical immergée dans un fluide
Page 42
CHAPITRE IV
IV.1.1. Equation du problème couplé

IV.1.1.1. Equations du comportement de poutre
On se place dans le cadre classique de l’hypothèse de Bernoulli, ce qui conduit au système

d’équations suivantes :
𝑬𝑰𝒗,𝒛𝒛𝒛𝒛+ 𝝆𝒔 𝑺𝒗 = 𝒇𝒉 𝒛, 𝒕 𝒑𝒐𝒖𝒓 𝒛 ∈ 𝟎, 𝒍
𝑴𝒇 = 𝑬𝒍𝒗,𝒛𝒛 𝒆𝒕 𝑻 = −𝑬𝒍𝒗,𝒛𝒛𝒛
𝒂𝒗𝒆𝒄 𝒍𝒆𝒔 𝒄𝒐𝒏𝒅𝒊𝒕𝒊𝒐𝒏𝒔 𝒍𝒊𝒎𝒊𝒕𝒆𝒔 𝒔𝒖𝒊𝒗𝒂𝒏𝒕𝒆𝒔 𝒒𝒖𝒊 𝒄𝒐𝒓𝒓𝒆𝒔𝒑𝒐𝒏𝒅𝒆𝒏𝒕 à 𝒖𝒏𝒆 𝒑𝒐𝒖𝒕𝒓𝒆 𝒆𝒏𝒄𝒂𝒔𝒕𝒓é𝒆 − 𝒍𝒊𝒃𝒓𝒆
𝒗 𝟎, 𝒕 = 𝟎
𝒗,𝒛 𝟎, 𝒕 = 𝟎
𝒗,𝒛𝒛 𝒍, 𝒕 = 𝟎 (4-1)
𝒗,𝒛𝒛𝒛 𝒍. 𝒕 = 𝟎
𝒗(𝒛, 𝒕)est la flèche suivant l’axe 𝒙, 𝑬le module d’Young du matériau, 𝝆𝒔 sa masse
volumique, 𝑰le moment quadratique de la section droite et𝑺l’aire de la section.
𝒇𝒉 (𝒛, 𝒕)est la résultante linéique suivant𝒙des pressions :
𝟐𝝅
𝒇𝒉 𝒛, 𝒕 = − 𝟎
𝒑 𝒂, 𝜽, 𝒛, 𝒕 𝒂𝒄𝒐𝒔𝜽𝒅𝜽 𝒑𝒐𝒖𝒓 𝒛 ∈ 𝟎, 𝒉 (4-2)
𝒇𝒉 𝒛, 𝒕 = 𝟎 𝒑𝒐𝒖𝒓 ∈ 𝟎, 𝒍
IV.1.1.2. Equations du fluide
En supposant le fluide parfait et incompressible, et l’écoulementirrotationnel, on aboutit au

système d’équations suivant :
𝟏 𝝏𝟐 𝜱
∆𝜱 − = 𝟎𝒅𝒂𝒏𝒔𝜴𝒇
𝑪𝟐 𝝏𝒕𝟐
𝝏𝜱
= 𝒗𝒄𝒐𝒔𝜽 𝒑𝒐𝒖𝒓 𝒓 = 𝒂, ∀𝜽 𝒆𝒕 𝒛 ∈ 𝟎, 𝒉 (𝒄𝒐𝒏𝒅𝒊𝒕𝒊𝒐𝒏 𝒔𝒖𝒓 𝒍𝒂 𝒔𝒖𝒓𝒇𝒂𝒄𝒆
𝝏𝒏
𝝏𝜱
=𝟎 𝒑𝒐𝒖𝒓 𝒛 = 𝟎 𝒆𝒕 ∀𝒓 𝒆𝒕 𝜽 (𝒄𝒐𝒏𝒅𝒊𝒕𝒊𝒐𝒏 𝒔𝒖𝒓 𝒍𝒆 𝒇𝒐𝒏𝒅)
𝝏𝒏
𝚽=𝟎 𝒑𝒐𝒖𝒓 𝒛 = 𝒉 𝒆𝒕 ∀𝒓 𝒆𝒕 𝜽 (𝒄𝒐𝒏𝒅𝒊𝒕𝒊𝒐𝒏 𝒅𝒆 𝒔𝒖𝒓𝒇𝒂𝒄𝒆 𝒍𝒊𝒃𝒓𝒆)
𝝏𝜱 𝟏 𝝏𝜱
𝐥𝐢𝐦 + =𝟎 (𝒄𝒐𝒏𝒅𝒊𝒕𝒊𝒐𝒏 𝒅𝒆 𝑺𝒐𝒎𝒎𝒆𝒓𝒇𝒊𝒆𝒍𝒅)
𝒓→∞ 𝝏𝒏 𝑪 𝝏𝒕
𝒗𝒄𝒐𝒔 𝜽est la vitesse normale du point𝑴appartenant à la surface latérale de la poutre.
Page 43
CHAPITRE IV
La pression en un point quelconque du fluide et en particulier sur la surface latérale mouillée

de la poutre s’obtient par la relation de Lagrange linéarisée. Si on connait le potentiel𝚽, on
peut calculer la résultante des efforts hydrodynamique𝒇𝒉 𝒛, 𝒕
IV.2. Modélisation numérique et couplage
Dans cette partie de notre travail nous donnons une méthode numérique afin de résoudre ce
genre de problème couplé en utilisant une méthode donnée par la réf[42]qui nous permet de
déterminer la masse ajoutée.
IV.2.1. Discrétisation de la poutre
La poutre est discrétisée en éléments poutre, c’est un élément à deux nœuds de numéros
courants𝒋et𝒋 + 𝟏de coordonnées𝒛𝒋 , 𝒛𝒋+𝟏 et de longueur𝑳𝒋 . Sur cet élément la
flèche𝒗(𝒛, 𝒕)est interpolée de façon classique en choisissant une approximation cubique

(polynôme de l’Hermite cubique) soit
𝒗 𝒛, 𝒕 = 𝑵𝒔 𝒂 𝒆 (4-3)
Avec :
𝟏
𝒂 𝒆 = 𝒗𝒋 , 𝜽𝒋 , 𝒗𝒋+𝟏 , 𝜽𝒋+𝟏 (4-4)
𝑵𝒔 Sontles fonctions de forme, données par :
𝑵𝒔 = 𝟏 − 𝟑𝒔𝟐 + 𝟐𝒔𝟑 , 𝒍𝒋 𝒔 − 𝟐𝒔𝟐 + 𝒔𝟑 , 𝟑𝒔𝟐 − 𝟐𝒔𝟑 , 𝒍𝒋 (−𝒔𝟐 + 𝒔𝟑 ) (4-5)
𝒛−𝒛𝒋
Où𝒔 = ∈ 𝟎, 𝟏 (4-6)
𝒍𝒋
Figure (IV-2) : Élément poutre à 2nœuds
Page 44
CHAPITRE IV
A partir des équations (2-5) et (2-7), on calcul la matrice élémentaire de rigidité 𝑲 𝒆 et la

matrice masse 𝑴 𝒆 .
𝟏𝟐 − 𝟔𝑳 − 𝟏𝟐 − 𝟔𝑳
𝝆𝒔 × 𝑨 −𝟔𝑳𝟒𝑳² 𝟔𝑳𝟐𝑳²
𝑲 𝒆 =
𝑳𝟑 −𝟏𝟐𝟔𝑳𝟏𝟐𝟔𝑳
−𝟔𝑳𝟐𝑳² 𝟔𝑳𝟒𝑳²
𝟏𝟓𝟔 − 𝟐𝟐𝑳 − 𝟓𝟒 − 𝟏𝟑𝑳

𝑳 −𝟐𝟐𝑳𝟒𝑳² − 𝟏𝟑𝑳 − 𝟑𝑳²
𝑴 𝒆 = 𝝆𝒔 ×
𝟒𝟐𝟎 𝟓𝟒 − 𝟏𝟑𝑳𝟏𝟓𝟔𝟐𝟐𝑳
𝟏𝟑𝑳 − 𝟑𝑳² 𝟐𝟐𝑳𝟒𝑳²
IV.2.2. Discrétisation du milieu fluide
Compte tenu de la géométrie de révolution d’axe𝒛, le potentiel des vitesses est déterminé en
utilisant un développement en série de Fourier donné par les références [42]
∞ ∞
𝚽 𝒓, 𝜽, 𝒛, 𝒕 = 𝚽𝒏 𝒓, 𝒛, 𝒕 + 𝒏=𝟏 𝚽𝒏 𝒓, 𝒛, 𝒕 𝒄𝒐𝒔𝒏𝜽 + 𝒎=𝟏 𝚽𝒏 𝒓, 𝒛, 𝒕 𝒔𝒊𝒏𝒎𝜽 (4-7)
L’harmonique𝚽correspond à un problème axisymétrique
Le domaine fluide est discrétisé par des éléments Q4 (élément quadrilatère à 4 nœuds) avec
un degré de liberté par nœud. La frontièreà l’infinie𝜞∞ est remplacée par une frontière
finie𝜞𝑹 , l’approximation nodale s’écrit :
𝚽 𝒓, 𝜽, 𝒛, 𝒕 =
∞
𝑵𝒇 𝒓, 𝒛 𝚽 𝟎 + 𝒏=𝟏 𝑵𝒇 𝒓, 𝒛 𝒄𝒐𝒔 𝒏𝜽 𝚽 𝒏 + (4-8)
∞
𝒎=𝟏 𝑵𝒇 𝒓, 𝒛 𝒔𝒊𝒏 𝒏𝜽 𝚽 𝒎 𝚽 𝒏 Vecteur des potentiels nodaux élémentaires pour
l’harmonique 𝒌du développement en série de Fourier (harmonique en cosinus ou sinus).
La formulation faible de type Galerkinconduit à une double infinité plus un système

matriciel ; chaque système correspond à une harmonique𝒏 :
𝑯𝒏𝒇 𝚽𝐧 = 𝐆𝐧
Page 45
CHAPITRE IV
Avec𝚽𝐧 vecteurdes potentiels nodaux du domaine fluide complet,𝑯𝒏 𝐧

𝒇 et 𝐆 matrices
obtenues par assemblage des matrices𝒉𝒏 𝐧

𝒇 et 𝐠 élémentaires telles que :
1
𝐡𝐢 𝐧 =𝛑 Ω𝑒 𝑁𝑓 𝑖,𝑟 𝑁𝑓 𝑗 ,𝑟 + 𝑁𝑓 𝑖,𝑧 𝑁𝑓 𝑗 ,𝑧 𝑟𝑑𝑟𝑑𝑧 + 𝜋𝑟 2 Ω𝑒 𝑁𝑓 𝑖 𝑁𝑓 𝑗 𝑟 𝑑𝑟𝑑𝑧𝑠𝑖𝑛 >0
𝐡𝐢 𝐧 = 𝟐𝛑 Ω𝑒 𝑁𝑓 𝑖,𝑟 𝑁𝑓 𝑗 ,𝑟 + 𝑁𝑓 𝑖,𝑧 𝑁𝑓 𝑗 ,𝑧 𝑟𝑑𝑟𝑑𝑧 𝐬𝐢𝐧 = 𝟎 (4-10)
𝐠𝐢 𝐧 = 𝚪𝐞 𝛜𝚪𝐬𝐟 𝑐𝑜𝑠𝑛𝜃 −𝑣 ∙ 𝑛 𝑑Γ𝑜𝑢 𝑔𝑖 𝑛 = 𝚪𝐞 𝛜𝚪𝐬𝐟 𝑠𝑖𝑛𝑛𝜃 −𝑣 ∙ 𝑛 𝑑Γ
L’orthogonalité des fonctions trigonométriques conduit au découplage des différents

systèmes. En effet
𝟎 𝒔𝒊𝒏 ≠ 𝒑
𝟐𝝅
𝟎
𝒄𝒐𝒔 𝒏𝜽𝒄𝒐𝒔𝒑𝜽𝒅𝜽 = 𝟐𝝅𝒔𝒊𝒏 = 𝒑 = 𝟎
𝝅𝒔𝒊𝒏 = 𝒑 ≠ 𝟎
(4-11)
𝟐𝝅
𝟎 𝒔𝒊𝒎 ≠ 𝒒
𝒔𝒊𝒏𝒎𝜽𝒄𝒐𝒔𝒒𝜽𝒅𝜽 =
𝝅𝒔𝒊𝒎 = 𝒒
𝟎
𝟐𝝅 𝟐𝝅 𝟐𝝅
𝒄𝒐𝒔𝒏𝜽𝒔𝒊𝒏𝒒𝜽𝒅𝜽 = 𝒄𝒐𝒔𝒏𝜽𝒅𝜽 = 𝒔𝒊𝒏𝒑𝜽𝒅𝜽 = 𝟎

𝟎 𝟎 𝟎
Pour ce problème particulier de vibrations de poutre, comme nous l’avons déjà indiqué
précédemment, nous avons :
𝒗𝒏 = 𝒗 𝒄𝒐𝒔𝜽 (4-12)
Donc de tous les seconds membres 𝒈𝒊 𝒏 , seul celui correspondant à l’harmonique 𝒏 = 𝟏du
développement en cosinus est non nul et le problème fluide se ramène à la résolution d’un
seul système linéaire
𝑯𝟏𝒇 𝚽𝟏 = 𝐆𝟏 (4-13)
Avec le second membre qui s’obtient par assemblage de termes
𝒈𝒊 𝟏 = 𝜞𝒆𝝐𝜞𝒔𝒇 𝑵𝒇𝒊 𝒓=𝒂 𝒄𝒐𝒔𝜽 −𝒗𝒄𝒐𝒔𝜽 𝒅𝜞 (4-14)
Page 46
CHAPITRE IV
𝒈𝒊 𝟏 = 𝝅𝒂 𝑵𝒇 −𝒗 𝒅𝒛 (4-15)
𝒓=𝒂
𝑲+𝟐
𝑲+𝟏 𝜞𝑹
𝑴
𝒌
𝑳
𝒂 LJ 𝑹
R
R
Figure (IV -3) maillage du
domaine fluide
Les conditions affichées sur𝜞𝑹 sont ou bien les conditions𝝋 = 𝟎ou bien les conditions
𝝏𝚽
=𝟎 (4-16)
𝝏𝒏
Compte tenu de la discrétisation en éléments finis poutre et en écrivant le systèmeen fonction

des pressions nodales, on obtient
𝑯𝑭 𝒏 𝒑 = 𝝅𝝆𝑭 𝒂 𝑭𝑺 𝒂 (4-17)
Avec[𝑭𝑺]assemblage de matrices [𝒇𝒔]𝒆telles que si les maillages des discrétisations poutre

et fluide coïncident, elles s’expriment par
𝑭𝒔𝒊𝒋 = 𝑵𝒇𝒊 𝑵𝒔 𝒋 𝒅𝒛 (4-18)

𝒆 𝒓=𝒂
Le calcul de la pression sur la surface latérale de la poutre et le calcul de la force généralisée

associée conduit au système matriciel qui traduit le comportement de la poutre immergée
𝑻
𝑴𝑺 𝒂 + 𝑲𝑺 𝒂 = −𝝅𝒂 𝑭𝑺 𝒑 (4-19)
Page 47
CHAPITRE IV
En éliminant les variables de pression (résolution du systèmeavec une matrice du

système[𝑯𝒓 ]tenant compte des conditions de pression nulle sur la surface libre), on en déduit
la matrice « masse ajoutée »
𝑴𝒂 = 𝝅𝟐 𝝆𝒇 𝒂𝟐 𝑭𝑺 𝑯𝒓 −𝟏
𝑭𝑺 (4-20)
En choisissant une discrétisation fluide en éléments quadrilatérales à 4 nœuds (modèle

linéaire), à cotés parallèles aux axes 𝒓𝒆𝒕𝒛,avec (𝒓𝟎 , 𝒛𝟎 ) les coordonnées du centre de
l’élément et𝟒𝒂, 𝟐𝒃les longueurs des cotés, l’interpolation du potentiel des vitesses est
donnée par :
𝚽 =< 𝑵𝒇 > 𝒄𝒐𝒔𝜽 𝚽 𝒆 (4-21)
Où les fonctions Nf sont les fonctions linéaires.
Dans ces conditions, la matrice élémentaire 𝒉 = 𝝅( 𝒃 + 𝒄 )s’exprime par
𝟏 𝟏
𝟏 𝟏 (4-22)
𝒃𝒊𝒋 = 𝑵 𝑵 + 𝑵 𝑵 𝒓𝟎 + 𝒂𝝃 𝒂𝒃𝒅𝝃𝒅𝜼
−𝟏 −𝟏 𝒂𝟐 𝒇𝒓 𝒇𝒓 𝒃𝟐 𝒇𝒛 𝒇𝒛
Soit :
(4-23)
Page 48
CHAPITRE IV
Et
𝒂𝒃
𝒄𝒊𝒋 = ∬ 𝑵𝒊 𝑵𝒋 𝒅𝝃𝒅𝜼
𝒓𝟎 + 𝒂𝝃
(4-24)
Soit ;
Et dans ces conditions la matrice élémentaire 𝑭𝒔 𝒆 couplage entre un élément poutre (𝒆) est
un élément fluide rectangulaire tel que𝟐𝒃 = 𝑰𝒆 s’écrit
𝟕𝒃 𝒃𝟐 𝟑𝒃 𝟐𝒃𝟐
−
𝟏𝟎 𝟓 𝟏𝟎 𝟓
𝟎 𝟎 𝟎 𝟎
𝒇𝒔 𝒆 = = 𝒔𝒇 𝒆 (4-25)
𝟎 𝟎 𝟎 𝟎
𝟑𝒃 𝟐𝒃𝟐 𝟕𝒃 𝒃𝟐
−
𝟏𝟎 𝟓 𝟏𝟎 𝟓
Page 49
Page 50
ORGANISATION DE LA PROGRAMMATION
CHAPITRE V
V-1. Introduction :
Il existe une diversité de méthodes numériques, depuis ces dernières décennies, les
sciences des matériaux ont connu une renaissance grâce aux développements des outils
de simulation.
La simulation numérique est devenue une troisième voie de recherche entre la théorie et
l’expérience. A l’heure actuelle, les simulations numériques permettent de traiter la
majeure partie des problèmes physiques ouverts et les modèles utilisés ne cessent
d’évoluer. La possibilité de simuler des phénomènes toujours plus complexes entraîne
toutefois une augmentation du nombre de grandeurs physiques nécessaires à la
caractérisation des matériaux en utilisant des modèles numériques. Les résultats
numériques peuvent ainsi contribuer à l’acceptation ou bien au rejet des théories
analytiques et peuvent également indiquer les directions dans lesquelles de nouvelles
approches doivent être développées.
La résolution du problème de vibration d’un système couplé solide en interaction avec un

fluide, est ramené à résoudre un système d’équations différentielles qui nécessite
l’utilisation de ces méthodes numériques d’où la nécessité d’élaborer un programme de
calcul.
La machine utilisée pour la réalisation du programme de calcul et l’exécution des

exemples choisis est un PC HP Pentium dual- core, 2 Giga de RAM sous Windows
Vista. Le temps d’exécution est très rapide (quelques secondes)
Le langage technique de programmation utilisé dans ce travail est Matlab version

7.0.6.324 (R2008a).
V-2. Description du programme
Le programme élaboré suivant l’organigramme ci-dessous, permet de calculer des

matrices élémentaires de raideur [𝑲𝒔 ], de masse [𝑴𝒔 ] du solide pour chaque élément, la
matrice masse ajoutée [𝑴𝒂 ] du fluide, et les matrices d’assemblage, suivant le maillage
figure (𝑽 − 𝟏)application des conditions aux limites, ainsi que la résolution du système
des équations différentielles du problème couplé aux valeurs propres avec la commande
Matlab« eig » c'est-à-dire eigenvalues et eigenvectors.
Page 51
CHAPITRE V
Figure (V-1) : Schéma de maillage et de numérotation de la structure
Nous avons écris le programme comme un seul bloc (faisant appel à des fonctions)
suivant la structure suivante
 Introduction des données paramètres physiques et géométriques,

 Lecture des coordonnées des nœuds de l’élément fini,
 Connectivité des éléments,
 Initialisation des matrices,
 Calcul des matrices élémentaires,
 Extraction des coordonnées,
 Assemblage et obtention du système global,
 introduction des conditions aux limites,
 résolution du système d’équations différentielles par la commande« eig »qui
donne les valeurs et vecteurs propres et par la suite les fréquences propres du
système (ou modes).
Data input et symboles utilisés :
 Numérotation de tous les nœuds et tous les éléments de maillage
nnel : nombre de nœuds par élément.
ndln : nombre de degrés de liberté par nœud.
ndle : nombre de degrés de libertés par élément.
nnt: nombre de nœuds total.
nelt : nombre d’éléments total.
ndlt : nombre de degrés de libertés total.
 Coordonnées et connectivités :
Coor (i) : les coordonnées de chaque nœud.
connec (i,j) : connectivité entre les éléments.
ie : numéros de l’élément.
Page 52
CHAPITRE V
V-3. Résolution
En Suivant l’ordre de calcul des différents éléments qui composent l’expression des
équations différentielles (2-36) données au chapitre IV:
𝑴𝒔 + 𝑴𝒂 𝒅 + 𝑲𝒔 𝒅 = 𝟎 (5-1)
On est ramené à résoudre le système d’équations différentielles (2-37) données au
chapitre IV :
𝑲 − 𝝎𝟐 𝑴 𝒅 = 𝟎 (5-2)
Avec 𝑴 = 𝑴𝒔 + 𝑴𝒂 et 𝝎𝟐 fréquences du système étudié
La résolution de l’eq. (4-1) se réalise sousMatlab,qui fait appel à des fonctions sous
programmées.
Les matrices 𝑲 et 𝑴 sont deux matrices symétriques de taille 𝒏𝒙𝒏(n= nombre de degrés
de liberté total non nul), on procède à un changement de variable 𝛚𝟐 = λ,nous obtenons
alors un système à valeurs et vecteurs propres sous lafrome :
[K-λM]d=0 (5-3)
On écrit cette équation sous la forme
𝑲𝒅 = 𝝀𝑴𝒅 (5-4)
La méthode utilisée est « la méthode des itérations inverses ».
Le problème est transformé en forme standard avant d’être résolu par la méthode
d’orthogonalisation de Jacobi.
𝑴est défini positive et symétrique une décomposition de Choleskypeut être appliquée :
𝑴 = 𝑳𝑳𝑻
Où 𝑳est une matrice triangulaire inférieure et𝑳𝑻 sa matrice transposée. (5-5)
On introduit la transformation suivante :
𝒅 = 𝑳−𝟏 𝑻 𝒛 (5-6)
En substituant dans l’équation (4-2) on a :
𝑲 𝑳−𝟏 𝑻 𝒛 = 𝝀𝑳𝑳𝑻 𝑳−𝟏 𝑻 𝒛 (5-7)
Page 53
CHAPITRE V
En multipliant cette équation par la matrice inverse 𝐿−1 nous obtenons
𝑳−𝟏 𝑲 𝑳−𝟏 𝑻 𝒛 = 𝝀𝑳−𝟏 𝑳𝑳𝑻 𝑳−𝟏 𝑻 𝒛 (5-8)

On note que 𝑳−𝟏 𝑳 = 𝑳𝑻 𝑳−𝟏 𝑻
=𝑰
L’équation (4-8) est exprimée sous la forme standard :
𝑯𝒛 = 𝝀𝒅 (5-9)
Où :
𝑯 = 𝑳−𝟏 𝑲 𝑳−𝟏 𝑻 (5-10)
L’algorithme de résolution est comme suit :
1. Décomposition de cholesky𝑴 = 𝑳𝑳𝑻

2. Calcul de 𝑳−𝟏 (matrice triangulaire inférieure inverse)
3. Calcul de 𝑯à partir de l’équation (5-10)
4. Résolution du problème de type valeurs-vecteurs propres 𝑯𝒛 = 𝝀𝒛 (en utilisant la
méthode du Jacobi)
5. Retrouver les vecteurs propres du problème original par l’équation (5-6)
𝑿 = (𝑳−𝟏 )𝑻 𝒁 (5-11)
Page 54
CHAPITRE V
V.4. Organigramme
Début
Data input
Entrées des données
Initialisation des matrices
Introduction des conditions

auxlimites
Extraction des degrés de liberté

pour chaque nœud
Formation des matrices raideurs

élémentaires de chaque élément
Formation des matrices masses

élémentaires de chaque élément
Formation des matrices masses

ajoutées élémentaires de chaque
élément
Calcul de masse total
Assemblage des matrices
Application des conditions aux

limites
Résolution
Data output
Résultats (fréquences propres)
Page 55
Page 56
RESULTATS ET INTERPRETATION
CHAPITRE VI
VI-1 Introduction
Plusieurs paramètres peuvent influencer la vibration d’une poutre. Nous allons axer
notre étude sur la variation des paramètres géométriques longueur sur rayon (L/r), et
hauteur d’immersion sur la longueur (H/L) grandeurs illustrées sur la figure (VI-1), la
variation des moments d’inerties pour quelques poutres de sections différentes ainsi
que la variation du nombre de masse (rapport de la masse volumique sur la masse
volumique du fluide), et observer leurs influences sur les caractéristiques dynamiques
de la poutre, mais il faut en premier temps procéder à la validation du programme
élaboré sous Matlab, par une étude de la convergence et une comparaison avec les
résultats de la littérature.
VI-2 Description de la poutre
Les caractéristiques géométriques et matérielles de la plaque étudiée sont définies

comme suit :
 𝑳 : longueur de la poutre
 𝒓 : rayon de poutre
 𝑯: hauteur immergée
 Coefficient de poisson 𝝂 =0.3
 Module de Young E=206 GPA
 ρ= 7830 Kg/m3 : masse volumique du solide
 ρf = 103 Kg/m3 :masse volumique du fluide
 Rapport de surfaces pour différentes sections (section circulaire/section

circulaire, section triangulaire, rectangulaire)λ : 1, 1.08, et 1.42
 Rapport des moments d’inertie (moment d’inertie de la section circulaire/le

moment d’inertie des différentes sections) α : 0.2, 0.5, et 1
 Rapport de la longueur sur le rayon L/r : 2, 2,5, 4, 5 et 10.
 Rapport de masse volumique Rm= rhos/rhof : 7.57, 7.8, 9.28, 9.39
 Rapport de la longueur sur la hauteur immergée H/L: 1/8, 1/5, 1/4, 1/2, et 1.
Page 57
CHAPITRE VI
Figure VI-1 : poutreconsole partiellement immergée
VI-3 Validation du programme
Afin de s’assurer que le programme est valide, une étude de convergence est établie.
Les résultats que nous avons obtenus par la programmation sont comparés avec les
résultats donnés par les références.
Pour cela nous considérons une poutre console figure (VI-1) : nous calculons les
paramètres de fréquences ainsi que les fréquencesde vibration pour quatre modes.
On donne sur le tableau (VI-1) la variation des valeurs des fréquences des quatre
premiers modes en fonction du nombre d’éléments pour une poutre de longueur L=1,
r=0,1 :
Page 58
CHAPITRE VI
Nombre d’éléments
mode 10 15 20 30 40 50 100
1 186 186 186 186 186 186 186
2 1165.7 1165.7 1165.6 1165.6 1165.6 1165.6 1165.6
3 3265.1 3264 3263.9 3263.8 3263.8 3263.8 3263.8
4 6404.9 6397.4 6396.3 6395.9 6395.9 6395.9 6395.9
Tableau VI-1 :convergence des fréquences (Hz) en fonction du nombre d’éléments
Nous constatons sur le tableau (VI-1), plus le maillage est raffiné (nombre d’élément
plus grand et la taille de l’élément est plus petite) plus les valeurs des fréquences
convergent vers une valeur constante ; ces valeurs que nous pouvons comparer aux
valeurs références données par le tableau (VI-2).
6000
Mode1
5000 Mode2
Mode3
Mode4
Fréquences (Hz)
4000
3000
2000
1000
0
10 20 30 40 50 60 70 80 90 100
Nombre d'éléments
Figure VI-2 : courbes de convergence
Le tracé des courbes de convergence nous montre que les valeurs des fréquences
convergent vers une valeur constante en augmentant le nombre d’éléments. Donc
nous concluons que la convergence est vérifiée.
Nous donnons sur le tableau (VI-2) une comparaison de nos résultats aves les
résultats analytiques et les résultats obtenus par calcul éléments finis donnés par la
littérature ref[43] .
Page 59
CHAPITRE VI
Fréquences naturelles (Hz) pour r=0.1m

ε% ε%
par par
Résultats Résultatsanalytiques rapport rapport
Mode
obtenus Réf [43] Résultats aux aux
MEF résultats résultats
Réf [44] analytiques MEF
1 186 185.91 185.91 0.004 0.004
2 1165.6 1165.1 1165.1 0.004 0.004
3 3263.8 3262.6 3262.5 0.003 0.003
4 6395.9 6394.3 6393.6 0.002 0.002
Tableau VI-2: valeurs des fréquences naturelles (Hz)
Où ε% : est l’écart relative (précision) entre les résultats que nous avons obtenu à
l’aide du programme réalisé et ceux desréférences.
Le tableau (VI-3) représente les fréquences naturelles, nous remarquons que les
valeurs des fréquences se rapprochent des valeurs de référence [44], avec une
précision variant de 0.002% à 0.004%,ce qui montre que ces résultats numériques
sont en bonne concordance avec les valeurs analytiques et les valeurs obtenues par la
MEF [44], nous concluons que notre programme est validé.
VI-4 Influence des paramètres géométriques et physiques sur les fréquences

propres de vibration de la poutre
Etude de cas :
Etudions l’influence des différents paramètres géométriques sur les fréquences des
quatre premiers modes.
Pour cela nous considérons une poutre encastrée-libre de section circulaire,

rectangulaire et triangulaire pour le premier cas, et circulaire pour les cas suivants,
nous faisons varier les paramètres géométriques suivants :
Cas1 : variation des fréquences pour une même surface pour différentes inertie et
pour une même inertie pour différentes valeurs de surfaces.
Le tableau (VI-3a) nous donne les variations des fréquences propres de vibration de la
poutre pour différentes sections (circulaire, rectangulaire, et triangulaire), pour une
même valeur d’inertie. Ces valeurs nous permettent de tracer le graphe représenté
par la figure (VI-3a).
Page 60
CHAPITRE VI
sections Rapports des surfaces λ

rectangulaire triangulaire Circulaire
Modes
1.42 1.08 1
MODE1 323 278.2 275.8
Mode2 2032 1754.3 1735.4
Mode3 5751 4973.4 4896.9
Mode4 11432 9911.4 9697.8
Tableau VI-3a : variation des fréquences pour une même valeur d’inertie pour
des sections différentes
MODE1
330
320
Fréquences (Hz)
310
300
290
280
270
1.0 1.1 1.2 1.3 1.4
Rapport de surface λ
(a)
Page 61
CHAPITRE VI
MODE2
2050
2000
Fréquences (Hz) 1950
1900
1850
1800
1750
1700
1.0 1.1 1.2 1.3 1.4
rapports des surfaces λ
(b)
Figure VI-3a : Variation des fréquences en fonction du rapport des surfaces λ

pour différentes sections pour la même valeur d’inertie, (a) mode 1, (b) mode 2
Le tableau (VI-3b) nous donne les variations des fréquences propres de vibration de
la poutre pour différentes sections (circulaire, rectangulaire, et triangulaire), pour une
même valeur de surface. Ces valeurs nous permettent de tracer le graphe représenté
par la figure (VI-3b).
Sections Rapport des inerties α

circulaire rectangle triangulaire
Modes
1 0.5 0.2
Mode1 275.8 404.12 626
Mode2 1735.4 4811.4 3951
Mode3 4896.9 6818.8 1120.1
Mode4 9697.8 12754.3 2232.2
Tableau VI-3b :variation des fréquences pour une même valeur de surface pour
des sections différentes
Page 62
CHAPITRE VI
Nous avons tracé les graphes pour le premier mode et le deuxième mode.
MODE1
650
600
550
Fréquences (Hz)
500
450
400
350
300
250
0.2 0.4 0.6 0.8 1.0
Rapports d'inertieα
(a)
MODE2
4000
3500
Fréquences (Hz)
3000
2500
2000
1500
0.2 0.4 0.6 0.8 1.0
Rapports d'inertie α
(b)
Figure VI-3b : Variation des fréquences en fonction du rapport d’inertie α des

différentes sections pour la même valeur de surface, (a) mode 1, (b) mode 2
Dans les deux cas étudiés, nous remarquons que les fréquences des modes propres
diminuent lorsque le rapport des inerties augmente, et elles augmentent lorsque le
rapport des surfaces augmente cette augmentation des fréquences est due à des valeurs
de raideur plus importantes.
Page 63
CHAPITRE VI
Pour tous les cas suivent nous prenons un rapport de longueur/ rayon (L/r=5) et
nombre d’éléments égal à onze (11 éléments).
Cas 2 : variation des rapports L/r pour le cas d’une poutre de section circulaire
Le tableau (VI-4) nous donne les valeurs des fréquences pour les quatre premiers
modes de vibration en fonction de la variation des rapports la longueur de la poutre/le
rayon (L/r= 2, 2.5, 4, 5 et 10).
Rapport de L/r
2 2.5 4 5 10
Mode1 726 581 363 275.8 145.2
Mode2 4549 3639 2275 1735.4 909.8
Mode3 12738 10190 6369 4896.9 2547.6
Mode4 24962 19970 12481 9697.8 4992.4
Tableau VI- 4:Variation des fréquences en fonction du rapport L/r pour une poutre
de section circulaire
Mode1
800
700
600
Fréquences (Hz)
500
400
300
200
100
2 4 6 8 10
Rapport L/r
(a)
Page 64
CHAPITRE VI
Mode 2
5000
4500
4000
3500
Fréquences (Hz)
3000
2500
2000
1500
1000
500
2 4 6 8 10
Rapport L/r
(b)
Mode3
14000
12000
10000
Fréquences (Hz)
8000
6000
4000
2000
2 4 6 8 10
Rapport L/r
(c)
Mode4
25000
20000
Fréquences (Hz)
15000
10000
5000
2 4 6 8 10
Rapport L/r
(d)
FigureVI- 4 : Tracés des graphes des quatre modes pour les rapports L/r= 2, 2.2,
4,5, 10(a : mode1, b : mode2, c : mode3, d : mode4).
Page 65
CHAPITRE VI
On constate sur la figure (VI-4) que les fréquences diminuent au fur et à mesure que
le rapport L/raugmente. Ceci s’explique par le fait que la surface exposée au fluide
devient plus grande, d’où un effet d’inertie plus grand qui induit une masse ajoutée
plus grande.
Cas 3 : variation du rapport de masse Rm c.à.d. masse volumique du solide / la masse

volumique du fluide (rhos/rhof).
Le tableau (VI-5) nous donne la variation des fréquences pour quatre rapports de
masse Rm différents :7.57, 7.8, 9.28 et 9.39.
rapport de masse Rm= rhos/ rhof

7.57 7.8 9.28 9.39
Modes
Fréquences (Hz)
mode1 275.4 275.8 278 278.2
mode2 1733.1 1735.4 1748.1 1748.9
mode3 4891.3 4896.9 4927.3 4929.2
mode4 9689.3 9697.8 9743.4 9746.2
Tableau VI-5 : variation des fréquences en fonction du nombre de masse

Rm=rhos/rhof
Sur le tableau (VI-5), Nous remarquons que les fréquences de vibration augmentent
légèrement au fur et à mesure que le rapportRmcroit, cette augmentation des
fréquences est due à la diminution de l’inertie c.à.d. une diminution de la masse du
fluide.
cas 4 :variations des fréquences de vibration en fonction de la masse volumique du

matériau.
Pour ce cas nous utilisons trois matériaux différents : l’acier, cuivre et laiton.
La variation des fréquences est donnée par le tableau (VI-6).
Page 66
CHAPITRE VI
Matériaux
Modes
laiton cuivre acier
8918 8390 7800
Mode1 164.2 184.2 275.8
Mode2 1032.4 1158.4 1735.4
Mode3 2910.5 3267.1 4896.9
Mode4 5757.2 6465.8 9697.8
Tableau VI-6 : Variation des fréquences pour différents matériaux
Ces valeurs nous ont permis de tracer les graphes donnés par la figure (V-5), pour les
quatre modes, les fréquences sont données pour trois rapports masse volumique de
l’acier sur la masse volumique de l’acier, du cuivre et du laiton (1, 0.9, et 0.8).
10000
9000
Mode1
8000
Mode2
7000 Mode3
Fréquences (Hz)
Mode4
6000
5000
4000
3000
2000
1000
0
0.88 0.90 0.92 0.94 0.96 0.98 1.00
Rapport rho1/rhoi
Figure VI-5: variation des fréquences en fonction du rapport de masse

volumique de l’acier / masse volumique du cuivre et laiton
Nous remarquons que la valeur des fréquences augmente quand le rapport des masses
volumiques croit. Cette augmentation est due à une diminution de l’inertie qui
implique une masse moins grande.
Cas 5 :variation des fréquences en fonction de la variation du rapport de mase c.à.d.

rapport de la masse volumique du solide sur la masse volumique du fluide
Page 67
CHAPITRE VI
Rm=rhos/rhof et de la variation du rapport longueur de la poutre sur la hauteur

d’immersion H/L.
Pour trois rapports du nombre de masse (masse volumique du solide sur la masse
volumique du fluide : Rm= rhos/rhof=7.5728, 7.8, 9.2857 nous faisons varier le
rapport longueur d’immersion sur la longueur de la poutre H/L =1, 1 /2, 1/4, 1/5, et
1/8. Nous obtenons les valeurs des fréquences données par les tableaux (5-7).
Rapport rhos/rhof=7.5728
Rapport H/L
Modes 1 1/2 1/4 1/5 1/8
Mode1 282.8 289.1 289.9 290 290.1
Mode2 773.3 1814.5 1817.4 1817.9 1818.6
Mode3 4969.8 5084.2 5090.2 5091.3 5092.8
Mode4 9752.7 9967.5 9977.3 9978.9 9981.2
Rapport rhos/rhof =7.8

Rapport H/L
Modes 1 1/2 1/4 1/5 1/8
Mode1 282.8 289.1 289.9 290 290.2
Mode2 1773.3 1814.6 1817.5 1818 1818.6
Mode3 4969.8 5084.6 5090.4 5091.4 5092.8
Mode4 9752.7 9968 9977.5 9979.1 9981.3
Rapport rhos/rhof =9.28

Rapport H/L
Modes 1 1/2 1/4 1/5 1/8
Mode1 284 289.3 290 290.1 290.2
Mode2 1780.5 1815.4 1817.8 1818.2 1818.8
Mode3 4989.3 5086.2 5091.1 5092 5093.2
Mode4 9788.8 9970.7 9978.7 9980 9981.9
Tableau (VI-7) : variation des fréquences en fonction des rapports

d’immersion pour trois valeurs du nombre de masse Rm= rhos/rhof
Page 68
CHAPITRE VI
Les tableaux (VI-7) nous permettent de tracer les graphes de la variation des
fréquences représentés par la figure (VI-7). Nous avons tracé les graphes pour le
premier mode et le quatrième mode.
(a)
(b)
Figure VI-6: Variation des fréquences en fonction du rapport de masse

(Rm=rhos/rhof) et du rapport hauteur d’immersion sur la longueur de la poutre H/L
(a) :mode1,(b) :mode4
Page 69
CHAPITRE VI
Nous remarquons sur le graphe (VI-8) que les fréquences diminuent lorsque le rapport
de la hauteur d’immersion H/L augmente, ces mêmes fréquences diminuent lorsque le
rapport Rmdiminue. Cette augmentation est due à une augmentation de l’inertie,
générée par une masse ajoutée plus importante.
Pour un fluide léger, c.à.d. un rapport Rmgrandnous avons des fréquences plus
importantes.
L’effet de la masse ajoutée est d’autant plus marqué que la structure est plus
immergée et que le fluide est plus dense.
Au cours de cette investigation, et selon les cas étudiés, nous avons relevé la
caractérisation de l’effet inertiel. L’effet de couplage avec le fluide est caractérisé par
le seul effet de masse ajoutée : la baisse ou l’augmentation des fréquences de
vibration est due à l’effet de la masse ajoutée qui dépend des rapports selon les cas à
savoir le rapport longueur sur rayon de la poutre (L/r: : 2, 2,5, 4, 5 et 10, rapport de
masse volumique du solide/masse volumique du fluide (Rm= rhos/rhof : 7.57, 7.8,
9.28, 9.39) et longueur sur la hauteur d’immergée(H/L: 1/8, 1/5, 1/4, 1/2, et 1).
Page 70
CONCLUSION
Conclusion
Ce travail a permis d’étudier les vibrations libres d’une poutre immergée partiellement et
totalement dans un fluide. Nous déterminons les fréquences naturelles associés aux modes
propres.
Le mouvement des particules du fluide, au voisinage de la structure, donne une variation de

pression, on définit alors le couplage entre la structure et le fluide .Il a été nécessaire d’établir
une formation adaptée pour la résolution des équations structure-fluide, en utilisant des
méthodes numériques .Dans notre cas, on a utilisé le couplage « élément finis- éléments
finis ».La modélisation du solide et le fluide sont réalisés par la méthode des éléments finis.
On sait que les méthodes qui définissent l’interaction des deux milieux prennent en compte
la force d’inertie du fluide agissant sur la structure, dans ce cas nous utilisons l’hypothèse de
la masse ajoutée. Elle est déterminée et conditionnée a celle du solide. Le système d’équations
différentielles définit régit le comportement vibratoire de l’ensemble solide-fluide.Aprés avoir
validé le programme élaboré sur Matlab,nous avons étudié plusieurs cas qui nous ont permis
de déterminer l’influence des paramètres géométriques de la poutre sur les caractéristiques
dynamiques que nous pouvons résumer comme suit :
 L’augmentation du rapport longueur sur le diamètre engendre une diminution de la

fréquence c, à, d : une diminution de la rigidité de la structure.
 L’augmentation du rapport masse volumique du solide sur la masse volumique du
fluide engendre une augmentation de la fréquence c, à, d : une diminution de la masse
ajoutée.
 L’augmentation du rapport longueur de la poutre sur la hauteur d’immersion et la
diminution du rapport masse volumique du solide sur la masse volumique du fluide
engendre une diminution des fréquences c, a, d une augmentation de l’inertie générée
par une masse ajoutée plus importante
 L’augmentation du rapport d’inertie de différentes sections pour une même surface,
donne une diminution des fréquences c'est-à-dire : une diminution d’inertie
Dans un travail futur nous proposons d’étudier un fluide compressible, et donner une
comparaison de l’endommagement de structure sous l’influence de l’interaction solide-
structure pour un fluide compressible et incompressible.
Page 70
LISTE DES REFERENCES
Liste des références
[1]A. E. H. LOVE 1927 A Treatise on the Mathematical Theory of Elasticity. New York:
Dover Publications, Inc.
[2]J. W. STRUTT 1877 Theory of Sound. London: Macmillan Publications Co., Inc.
[3] R. M. DAVIES 1937 Philosophical Magazine, 563. The frequency of transverse vibration
of a loaded Fixed-free bar III. The effect of rotatory inertia of the bar.
[4]S. P. TIMOSHENKO 1921 Philosophical Magazine, 744. On the correction for shear of
the differential equation for transverse vibrations of bars of uniform cross-section.
[5]S. P. TIMOSHENKO 1922 Philosophical Magazine, 125. On the transverse vibrations

ofbars of uniform cross-section.
[6]Aquelet Nicolas Modélisation de l'impact hydrodynamique par un couplage fluide-

structure
[7]JEAN-NOEL GENCE. Professeur des Université, Lyon 1.MECANIQUE DES

FLUIDE.Les écoulements des fluides newtoniens. Principaux mécanisme associés aux
mouvements des fluides.Edition ellipse.
[8]A.A.LAKIS et M.P.PAIDOUSSIS ,1971.Free vibration of cylindrical shells partially filled

with liquid .Journal of Sound and Vibration ,V19,N1,P175-205.
[9]MEYERHOFF, W.K., 1970.Added mass of thin rectangular plates calculated from

potential theory. Journal of Ship research 14,100-111.
[10]B.PESEUX AND V. MARTIN. Numerical procedure for fluid structure interaction:

interface element.JM Crolet and R Ohayon editors
[11] Jean-François SIGRIST et Daniel BROC [Investigation of numerical method for modal
analysis of a tube bundle with fluid-structure interaction] ; Service Technique et Scientifique
DCN Propulsion 44620 la montagne, France
[12] Jean-François Sigrist, Christian Lainé et Bernard Peseux Analyse modale d’une structure
industrielle avec prise en compte du couplage fluide / structure] ; Mécanique & Industries 6,
553–563 (2005) AFM, EDP Sciences 2005 DOI: 10.1051 / meca:2005067.
[13]MyungjoJhung, WalTaekim et Yong Horyu[Dynamic characteristics of cylindrical shells
considering fluid-structure interaction]; Safety Research Division, Korea Institute of Nuclear
Safety 19 Guseong-dong, Yuseong-gu,Daejeon,305-338,Korea*Corresponding author. E-
mail: mjj@kins.re.kr Received August 3, 2009, Accepted for Publication September 16, 2009
Page 71
[14]Bratu C., Vibration des grandes structures marines, Paris, 1978, bulletin atma édition.
[15]Zienkiewicz O. eT Taylor R., The Finite Element Method-Solid and Fluid Mechanics:
Dynamics and Non-Linearity, London, McGraw-Hill, 1991, 4e edition vol. II
[16]Peseux B. et Quevat J.P., Linear fluid-structure interaction in unbounded domains, ASME

Pressure Vessels and Piping Conference, pp. 185-188, Nouvelle Orléans, 1985.
[17]Zienkiewicz O. eT Taylor R., The Finite Element Method-Basic Formulation and

Linear Problems, London,McGraw-Hill, 1989, 4e edition vol. I.
[18] Brid, Stewart et Lightfoot 1960, Cebeci et Brid.
[19]Shaw 1977, Shaprio 1953, Kuetheet Chow 1986.
[20] LAWDAN, Land LIFCHITZ, E(1971)-Mécanique des Fluides, Editions de Moscorv.
[21] Oudin H., Peseux B. et Quevat J.P., Dynamicresponse of coupledproblems, pp. 340-347,
Venise, Italie, 1984.
[22] Peseux B. et Martin V., NumericalProcedure for Fluid-Structure interactioin :

Interface elements. Computationan Methods for fluid-structure interaction, 1994,
Crolet/Ohayon-Longman Higher education edition, 78-103p.
[23]Axisa F.,Modélisation des systèmes mécaniques. Interactions fluide structure, Hermès,

2001.
[24] CIARLET P., Elasticité tridimensionnelle, Masson, 1986.
[25] DUVAUT G., Mécanique des milieux continus, Masson, 1990.
[26] SALEN9ON J., Mécanique des milieux continus, les éditions de l’Ecole polytechnique,
2002.
[27]BATHE K., Finite Element Procedures in Engineering Analysis, Prentice-Hall,

Englewood Cliffs, 1982.
[28]CIARLET P., The Finite Element Method for Elliptic Problems, North Holland, 1978.
Page 72
[29] DHATT G., TOUZOT G., Une présentation de la méthode des éléments finis, Maloine,
1984.
[30] SusbielleG. et Bratu C.H., Vagues et ouvrage pétroliers en mer, Paris, 1981,
technipedition, 78-103p.
[31] DHATT G., TOUZOT G., Une présentation de la méthode des éléments finis, Maloine,
1984.
[32] ZIENKIEWICZ O., TAYLOR R., The Finite Element Method Basic Formulation and
Linear Problems (4ème edition), Mac-Graw Hill, 1989.
[33]Axisa F.,Modélisation des systèmes mécaniques. Interactions fluide structure, Hermès,

2001.
[34] GIBERT R., Vibration des structures. Interaction avec les fluides. Sources d’excitation
aléatoires, Eyrolles, 1986.
[35] MORAND H.J., OHAYON R., Fluid-Structure Interaction, John Wiley & Sons, 1995.
[36] BERMUDEZ A., DURAN R., RODRIGUEZ R., « Finite Element Analysis of
Compressible and Incompressible Fluid-Solid System. », Mathematics of Computation,
Vol.67, n°21, p.111-136, 1998.
[37] HAMDI M., OUSSET Y., VERCHERY G., « A Displacement Method for the Analysis
of Vibrations of Coupled Fluid-Structure Systems», International Journal of Numerical
Methods in Engineering, vol.13, p. 139-150, 1978.
[38] Peseux B. et Delhommeau G., Fluid/structure interaction : Influ
[39] BOUJOUT J., « Mathematical Formulation of Fluid-Structure Interaction Problems »,

Mathematical Modeling and Numerical Analysis, vol.21, p. 239-260, 1987.
[40] SANDBERF G., GORANSSON P., «A Symmetric Finite Element Formulation for
Acoustic Fluid-Structure Interaction Analysis», Journal of Sound and Vibration, vol,
123, n°3, p. 507-515, 1988.
[41] MARCUS M., « A Finite Element Method Appiled to the Vibration of Submersed Plates
», Journal of Ship Research, vol.22, p. 94-99, 1978.
Page 73
[42] SIGRIST J.-F., GARREAU S., « Dynamic Analysis of Fluid-Structure Interaction

Problems with Spectral Method Using Pressure-Based Finite Elements », Finite Element
Analysis in Design, vol. 43, n°4, p. 287-300,2007.
[43] Sahi M., Ouset Y. et G.Verchery,Journal of Sound and Vibration, vol. 54, 1981, pp.
473-500.
[44] DE LANGERE E., Fluides et Solides, les éditions de l’Ecole polytechnique, 2001.
[45] Martin V. etPeseux B., Flowchart and software for numerical coupling of dynamic
problems in elastic media, International Journal for Numerical Methods in
Engineering, vol. 37 n°22, 1994, pp. 4285-4308.
[46] MOKHTARI M., MESBAH A., Apprendre et maitriser Matlab, Springer Verlag, 1997.
Page 74

Ibssmfra 14 Asm 351

Transféré par

Droits d'auteur :

Formats disponibles

Ibssmfra 14 Asm 351

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Ibssmfra 14 Asm 351

Transféré par

Droits d'auteur :

Formats disponibles

Nomenclature

Symboles Nomination Unités

Tableau VI-2 :Valeurs des fréquences naturelles (Hz)………………………………...…..60

Tableau VI-5 :Variation des fréquences en fonction du nombre de masse Rm=rhos/rhof….66

Tableau VI-6 : Variation des fréquences pour différents matériaux………………………...67

Tableau VI-7 :Variation des fréquences pour différents matériaux………………………...68

Figure I-1 : couplage structure-fluide et couplage structure fluide-fluide………………………...7

FigureI-2 :Solide élastique immergé dans un fluide……………………………………………...7

Figure I-3 : Couplage SPH/FEM…………………………………………………………………..8

Figure I-4 :Couplage SPH/FEM/FEM………………………………………………………….…8

Figure I-5 : Couplage FEM/FEM………………………………………………………………….9

Figure II.1 : Elément de poutre dans le plan XY…………………………………………………..12

Figure II.2:Schématisation générale du problème structure……………………………………..15

Figure II.3: Déplacement de la section droite d’abscisse x……………………………...……….16

Figure II.4:Hypothèse de Bernoulli...…………………………………………...…..……………18

FigureII.5:Schématisation générale du problème structure……………………………….............25

Figure II.6: Schématisation générale du problème structure-fluide………………………………25

Figure III.1:Schématisation générale d’un problème de couplage mécanique fluide-structure.....28

Figure III.2:Schématisation générale du problème fluide-structure……...………………………29

Figure III.3:Principe de maillage du problème couplé dans le plan (𝑟, 𝑧)…….………………..36

Figure IV.1:Poutre vertical immergé dans un fluide…………………………………………......42

Figure IV.2:élément poutre à deux nœuds………………………………………………………..44

Figure IV.3: maillage du domaine fluide………………………………………………………….47

Figure V.1:Schéma de maillage et de numérotation de la structure………………………………52

Figure VI-1 : poutreconsole partiellement immergée…………………………………………….58

Figure VI -2 : courbes de convergence…………………………………………………………....59

Figure VI-6:Variation des fréquences en fonction du rapport de masse (Rm=rhos/rhof) et du rapport

Chapitre II : Formulation mathématique

II-3 Théorie des poutres……………………………………………………….………16

Chapitre III :Modélisatuon par élément finis et couplage

III.3 Les éléments finis en mécanique du solide ……………………………………...30

III.4 Méthode des éléments finis en mécanique des fluides.(Forme faible de la

méthode des résidus pondérés)………………………………………………………...33

III.5 Couplage éléments finis-éléments finis………………………………………….35

Chapitre IV:Cas d'une poutre émergée dans un fluide

IV.1Etude des vibrations d’une poutre immergée dans un fluide…………………….42

IV.1.1Equation du problème couplé ………………………………………………43

IV.1.1.1Equations du comportement de poutre……………………………….....43

IV.2Modélisation numérique et couplage………………………………………...…44

IV.2.1 Discrétisation de la poutre……………………………………………………44

IV.2.2 Discrétisation du milieu fluide……………………………………………….45

Chapitre V:Organisation de la programmation

VI-3 Validation du programme……………………………………………...…58

L’étude présentée dans ce mémoire concerne l’analyse du comportement dynamique d’une

Le couplage définit par l’interaction fluide-structure nécessite la caractérisation de condition

Après avoir étudié la convergence et validé le programme avec plusieurs littératures, on a

The coupling defines by the interaction fluid-structure requires the characterization of

Key words:finite elements, coupling, interaction, added mass, differential equation,

‫‪ٚ‬زى ٔضغ ًَبرج انسبئم ٔ انصهت ثٕاسطخ طش‪ٚ‬مخ انؼُبصش انًزُبْ‪ٛ‬خ‪.‬‬

‫لذ اَشئذ انًؼبدالد انزفبضه‪ٛ‬خ نحشكخ انصهت دٌٔ ارصبنّ يغ انسبئم‪.‬‬

‫انؼُبصش انًزُبْ‪ٛ‬خ ‪.‬االلزشاٌ‪.‬انزفبػم‪.‬انكزهخ انًضبفخ‪.‬انًؼبدالد انزفبضه‪ٛ‬خ‪.‬انؼزجخ انًغًٕسح‪.‬‬

L’interaction fluide-structure est un domaine vaste il concerne un ensemble important de

Dés l’aube de l’industrialisation du monde moderne, le problème de l’interaction fluide-

 Couplage thermique : le fluide et le solide échangent de la chaleur.

 Couplage chimique : le solide et le fluide réagissent entre eux.

 Couplage mécanique : le solide et le fluide échangent des énergies mécaniques.

- Le fluide exerce des efforts sur la structure en mouvement, ce qui modifie la

Notre travail se divise en quatre parties principales:

Dans le second chapitre, nous procédons à La formulation du problème vibratoire d’une