Projet de Fin dEtude Juin 2009

Spcialit Gnie Civil

Modlisation de linteraction fluide/structure

[10]

Auteur : SCHOUBRENNER Mathieu

Elve ingnieur de 5me anne
Modlisation de linteraction fluide/structure

2
 Modlisation de linteraction fluide/structure

Remerciements

Je tiens tout dabord remercier M.Marc LEFRANC pour son accueil au

sein de Force Technology, mais aussi pour son enseignement sur lHistoire de la
Norvge.

Je remercie aussi messieurs Bogdan IWANOSKI et Rune GLADS pour

leurs aides prcieuses tout au long de ce projet.

Je tiens galement remercier, M. Pierre REGENASS, pour ses conseils

aviss et Mme Saida MOUHOUBI pour mavoir transmis la proposition de sujet
de PFE.

Je remercie lensemble des employs de Force Technology Norway pour

leur accueil au cours de ses 20 semaines de stage, et pour mavoir notamment
permis de participer lHolmenkollstafetten.

Je tiens finalement remercier YMERI Denion pour son aide de

traducteur au cours de ses 20 semaines.

3
 Modlisation de linteraction fluide/structure

Abstract

This project focuses on the capacity of LS-Dyna to compute fluid-structure

interaction. The main concern in this study is the coupling between the fluid and
the structure. In fact there are many types of coupling with relevant advantages
and disadvantages. The aim of this survey is, the analyze, the modelisation and
simulation of fluid structure interaction with the software LS-Dyna. LS-Dyna is not
developed in this goal.
Structure displacements are not negligible and its real impact on the fluid.
The fluid is considered compressible and Newtonian, represented in an Eulerian
configuration and governed by Navier Stokes equations. The structure is
assumed rigid or elastic, homogenous and isotropic, represented in a Lagrangian
configuration and governed by elasticity equation. An ALE method permits to
solve fluid equation in a mobile field fitting with structure deformation.

To compare and analyze LS-dyna results, we will relate to different

theories and other experiments
The aim of this study is the test of LS-Dyna software. To achieve this
project, different stages have been fulfilled.
-understanding of the problem and of its theory
-discovery of the LS-Dyna software
-working out of modeling principle and creation of different models
-analyze and check of the results

4
 Modlisation de linteraction fluide/structure

Rsum
Dans ce projet de fin dtudes, nous nous sommes intresss la
capacit de calcul du logiciel LS-Dyna dans le domaine dinteraction fluide
structure. La difficult principale dans ce type de problme est le couplage entre
la structure et le fluide. En effet, diffrents types de couplages existent avec leurs
avantages et inconvnients. Cette tude a pour objet d'analyser, de modliser et
de simuler diffrents problmes d'interaction fluide-structure avec le logiciel LS-
Dyna qui na pas t dvelopp a priori pour solutionner ce genre de problmes.
Les dplacements de la structure ne sont pas ngligeables et induisent
une relle dformation du domaine fluide. Le fluide est suppos compressible et
newtonien ; il est reprsent dans une configuration d'Euler par les quations de
Navier-Stokes. La structure est suppose rigide ou lastique, homogne et
isotrope ; elle est dcrite dans une configuration de Lagrange par les quations
de l'lasticit linarises. Une mthode ALE permet de rsoudre les quations
du problme fluide dans un domaine mobile s'adaptant aux dformations de la
structure.

Pour comparer les rsultats obtenus par le logiciel LS-Dyna au cours des
diffrentes simulations nous nous rapporterons diffrentes thories, mais aussi
des donnes exprimentales.

Le sujet de ce projet de fin dtude est de tester le logiciel LS-Dyna, et

pour parvenir cet objectif diffrentes tapes ont donc t suivi.
-comprhension du problme et de la thorie qui sy rapporte
-prise en main du logiciel LS-Dyna
-laboration dun principe de modlisation et ralisation de
diffrents modles.
-Analyse et comparaison des rsultats.

5
 Modlisation de linteraction fluide/structure

Sommaire
Remerciements.....................................................................................................3
Abstract.................................................................................................................4
Rsum ................................................................................................................5
Sommaire .............................................................................................................6
Liste des notations et symboles ............................................................................8

1 Introduction..................................................................................................9

2 Introduction gnrale ................................................................................10

2.1 Prsentation de linteraction fluide/structure.........................................10
2.2 Les enjeux industriels ...........................................................................10
2.2.1 Gnie civil......................................................................................10
2.2.2 Gnie maritime ..............................................................................12

3 Prsentation du logiciel LS-Dyna.............................................................13

3.1 LS-Dyna................................................................................................13
3.2 Organisation des commandes de programmation (keyword) ...............14
3.3 Fonctionnement thorique de LS-DYNA...............................................15

4 Modlisation du systme ..........................................................................17

4.1 Modlisation de linteraction fluide/structure.........................................17
4.2 Lagrangien-Eulerien Arbitraire (ALE)....................................................17
4.2.1 Formulation Eulerienne .................................................................17
4.2.2 Formulation Lagrangienne.............................................................18
4.2.3 Formulation ALE............................................................................18
4.3 Couplage fluide/structure......................................................................22
4.3.1 Gnralits ....................................................................................22
4.3.2 Procdure du couplage .................................................................22
4.3.3 Conditions et algorithme de couplage ...........................................23
4.4 Couplage par pnalit ..........................................................................24
4.5 Description des lments finis utiliss dans le modle fluide structure 25
4.5.1 Description gomtrique du maillage de la structure.....................25
4.5.2 Description de la mthode.............................................................26
4.5.3 Description du maillage du fluide...................................................26
4.6 Hourglass [annexe C] ...........................................................................27
4.6.1 Elment sous intgrs ...................................................................27
4.6.2 Mode Hourglass ............................................................................27

5 Vrification de la modlisation du fluide (hydrostatique) ......................29

5.1 Equation dtat .....................................................................................29
5.2 Densit du fluide...................................................................................30
5.3 Pression hydrostatique .........................................................................32
5.4 Test de ladvection ...............................................................................35

6
 Modlisation de linteraction fluide/structure

5.5 Efforts de leau .....................................................................................37

5.5.1 Rsultat .........................................................................................38
5.5.2 Conclusion.....................................................................................40

6 Vrification du comportement de la structure (cas statique) ................41

6.1 Dfinition de la plaque ..........................................................................41
6.2 Courbe du matriau ..............................................................................42
6.3 Thorie des plaques .............................................................................42
6.3.1 Introduction....................................................................................42
6.3.2 Equation diffrentielle ....................................................................43
6.4 1ere vrification ......................................................................................44
6.4.1 Modle...........................................................................................44
6.4.2 Rsultats .......................................................................................46
6.5 2eme vrification.....................................................................................46
6.5.1 Modle...........................................................................................47
6.5.2 Rsultats thorique........................................................................47
6.5.3 Rsultats et comparaison ..............................................................48
6.5.4 Plaque mince avec grand dplacement.........................................50
6.6 Etat de contrainte de la plaque .............................................................51
6.6.1 Tenseur des contraintes ................................................................51
6.6.2 Contraintes principales ..................................................................52
6.6.3 Pression isostatique ......................................................................56

7 Comportement dune structure sous leffet dynamique du fluide ........58

7.1 Dambreak exprience...........................................................................58
7.1.1 Prsentation de lexprience.........................................................58
7.1.2 Rsultats et comparaison avec lexprience .................................59
7.2 Autres simulations ................................................................................62
7.2.1 Force dimpact ...............................................................................63
7.2.2 Evolution de la force en fonction de la plaque ...............................66
7.3 Observation du couplage par pnalits ................................................69
8 Conclusion .................................................................................................70

Tableau des figures ............................................................................................72

Bibliographie .......................................................................................................74

Annexes .............................................................................................................75
A. Equations de Navier Stokes.........................................................................76
B. Algorithme dadvection de Van Leer............................................................77
C. Viscosit de Hourglass ................................................................................79
D. Ecoulement du fluide ...................................................................................82
E. Programme ..................................................................................................83
F. Contraintes de Von Mises............................................................................87

7
 Modlisation de linteraction fluide/structure

Liste des notations et symboles

-v : vitesse du fluide [m/s].

-p : pression [Pa].
- : oprateur gradient.
-t : variable temps.
- : densit [kg/m3].
- : domaine physique.
- : frontire du domaine physique.
-, , : coordonnes cartsiennes pour le passage llment de rfrence.
-x, y, z : coordonnes cartsiennes.
-g : rapport des chaleurs massiques.
-m : taux relatif de la fraction molaire.
-ref : rfrence.
-w : flche [m].
-D : rigidit flexionnelle de la poutre [N.m]
-F : force [N].
-M : moment [N.m].
-E : module dYoung [Pa]
- : contrainte [Pa]
- u : coefficient de poisson
- h, t : paisseur [m]
-ALE, Arbitraire Lagrange Euler

8
 Modlisation de linteraction fluide/structure

1 Introduction
La dtermination des efforts exercs sur une structure par un fluide est
dpendante de lcoulement de celui-ci et des dformations de la structure. On
parle de systmes coupls et dinteraction fluide/structure lorsque la structure ne
peut pas tre considre comme un lment gomtrique (cas des structures
rigides pouvant tre assimiles un lment gomtrique, o on ne parle plus
de couplage). Le couplage apparat systmatiquement ds quune structure
vibrante/lastique est immerg dans un fluide au repos ou non
Dans ltat actuel des standards de calcul, linteraction entre fluide et
structure nest tenue en compte que dans le cas de structure de grande
dimension (thorie de la diffraction). Pour les structures de petite taille
linteraction au sens de couplage est nglige, les effets du chargement dun
fluide devant tre considrs indpendamment de la taille de la structure.
Cependant la capacit des logiciels actuels est telle quil est possible de
simuler lcoulement du fluide en tenant compte de la prsence de la structure et
des dformations de celle-ci sous laction du fluide.

On se propose ici de tester le logiciel LS-DYNA dans loptique de valider

sa capacit rpondre ce type de problmes. En effet, la socit Force
Technology souhaite savoir si elle peut se servir de ce logiciel lors de phases de
dimensionnements de structures o le problme dinteraction fluide structure
apparait. En sachant que ce logiciel na pas t dvelopp dans cette optique et
que sa fonction principale et la simulation de crash.
On calculera les efforts sur la structure dans lhypothse de structure
lastique, soumis un fluide en mouvement. Aussi une premire tude en phase
statique nous permettra dtre une premire tape dans le test de LS-Dyna.
Ds lors les rsultats obtenus seront confronts avec les rsultats
exprimentaux, ou avec les diffrentes thories appropries. En effet pour les
problmes hydrostatiques nous nous ramnerons diffrentes thories,
notamment la thorie des plaques de Timoshenko. Dans le cas dun fluide en
mouvement nous nous rapporterons notamment une exprience ralise au
sein de luniversit de Groningen aux Pays-Bas.

Lobjet de ce rapport est de prsenter les rsultats des diffrents tests

raliss avec le logiciel LS-Dyna concernant des problmes dinteraction fluide
structure. Il permettra aussi dnoncer les bases de la thorie dinteraction fluide
structure

9
 Modlisation de linteraction fluide/structure

2 Introduction gnrale

2.1 Prsentation de linteraction fluide/structure

La simulation numrique de phnomnes coupls a connu un essor

constant ces dernires annes. Ce dveloppement est d en particulier aux
succs prcdents de la simulation numrique en gnral, mais aussi
l'accroissement permanent des performances des calculateurs.
Parmi ces phnomnes coupls se trouvent les interactions fluide-
structure. Elles mettent en jeu une structure toujours mobile, rigide ou
dformable, et un fluide liquide ou gazeux, en coulement autour ou contre une
partie de la structure. Ces phnomnes sont dits coupls, parce que l'volution
de chacun des deux lments dpend de celle de l'autre. Ainsi par exemple, la
forme de la voile d'un bateau (en rgime permanent) dpend de l'coulement de
l'air autour de celle-ci. Rciproquement, cet coulement dpend de la forme de la
voile.
On pourrait citer un trs grand nombre d'exemples o la gomtrie de la
structure agit sur lcoulement et les efforts du fluide sur la structure. Parmi
ceux-ci, on peut exhiber entre autres les phnomnes hydro-lastiques (fluide en
phase liquide) : coulements autour d'un navire, d'un sous-marin, d'une digue
dans un port ou de piles de pont, coulements liquides l'intrieur de conduites,
mouvements de liquides dans un rservoir, etc... On distingue galement les
phnomnes aro-lastiques o le fluide est en phase gazeuse : coulements
autour des vhicules ariens (avions, missiles, etc...) et terrestres (trains
grande vitesse, automobiles, etc...), influence du vent sur les constructions
souples (ponts suspendus, rfrigrants de centrale nuclaire, etc...).

2.2 Les enjeux industriels

2.2.1 Gnie civil

Un des premiers grands secteurs o les simulations exprimentales (et

numriques par la suite) d'interactions fluide-structure sont apparues est le Gnie
Civil, notamment pour les structures trs lances (ponts suspendus) et
faisceaux de cbles. Le dveloppement de ce domaine d'application est li
l'effondrement du pont de Tacoma (1940). Ce pont, conu dans les annes 1930
avec un budget rduit aprs la crise de 1929, tait trs lanc et fin pour
l'poque. Sa structure rvolutionnaire trs lgre lui permettait nanmoins de ne
pas tre susceptible de connatre le mme sort que le Tay Bridge: le 28
dcembre 1879, la trave centrale du Tay bridge s'est effondre Dundee, tuant

10
 Modlisation de linteraction fluide/structure

75 personnes, sous l'effet "statique" d'un vent perpendiculaire de force 10-

11(vitesse de vent de 89 117km/h).

Fig. 1: Pont de Tacoma [3]

Cependant, la matrise des effets statiques n'est pas suffisante. Pour le

pont de Tacoma ses concepteurs en ont fait la dure exprience. Sous l'effet d'un
vent constant de 68 km/h, celui-ci s'est effondr aprs avoir connu de
nombreuses oscillations d'amplitudes croissantes (certains haubans ont alors fini
par casser, entrainant la ruine du pont). Il est intressant de noter que ce pont
tait rput pour tre trs mobile, longtemps avant sa destruction. Si cet
effondrement n'tait pas survenu il est vident que la fatigue des lments
structurels du pont serait apparue de manire trs prmature.
A cette poque, les tudes taient uniquement de nature exprimentale.
En soufflerie, on valuait diffrentes caractristiques arolastiques de profils-
types (sections droites, surfaces polies et angles non mousss). Lorsque seule
la simulation exprimentale est utilise, le processus de dfinition de la section
de pont est long et coteux. En effet, pour obtenir des rsultats fiables, il faut
construire des maquettes ralistes, gnralement une chelle de l'ordre du
1/50me, des ponts envisags (sachant que, de toute faon, les expriences en
soufflerie exprimentale ne sont pas parfaites et ne peuvent pas respecter toutes
les similitudes rgissent le phnomne physique. Aussi les paramtres sans
dimension sont propres chaques modlisations car dpendant de la
modlisation du phnomne physique).
Avec l'accroissement de la puissance des ordinateurs, la simulation
numrique est de plus en plus utilise pour l'analyse arolastique de sections
de ponts. La soufflerie numrique permet en tout cas de rduire
considrablement les cots et les dures des tudes prliminaires pour des
projets trs lancs ou trs grands. Ces rsultats numriques ne reprsentent
que des simulations lmentaires autour de profils de ponts. Pour l'instant, ces
simulations ne peuvent absolument pas se substituer aux simulations en
soufflerie exprimentale, mme si celles-ci sont extrmement coteuses. En effet,
des simulations numriques compltes ( plusieurs chelles gomtriques) et
tridimensionnelles autour de structures compltes, restent toujours dlicates.
Cela nanmoins volue rapidement. Il faut aussi dire que, le jour o celles-ci

11
 Modlisation de linteraction fluide/structure

seront informatiquement plus aises, la simulation globale aura progress, car

les simulations en soufflerie exprimentale comportent de nombreux dfauts
(mesures, capteurs, erreurs de similitude, etc.) A l'heure actuelle, les simulations
numriques sont capables d'acclrer le design grossier des sections. Pour le
design fin, il faut toujours au moins une validation exprimentale incluant le
maximum de dtails.

2.2.2 Gnie maritime

Dans le domaine du Gnie Maritime, de nombreuses tudes numriques

concernent diffrents aspects des interactions fluide-structure. La plupart des
tudes font jouer un rle particulier la surface libre de la mer. Une des
sollicitations les plus importantes pour des constructions portuaires ou des
structures pour l'extraction off-shore du ptrole est simplement la houle (ondes
sans dplacement de matire la surface libre).
Pour les installations ptrolires, il existe une trs grande diversit de
structures envisages. Certains projets comportent des plateformes fixes
rigidement ou non au sol sous-marin, des barges flottantes maintenues par des
cbles, des pipe-lines verticaux permettant de ramener le ptrole vers la surface,
ventuellement munis de flotteurs intermdiaires. Il existe aussi des projets de
jacket de trs grande taille
Pour toutes ces configurations, les sollicitations de la houle en surface
sont trs dimensionnantes, mais il existe d'autres aspects fluide-structure ne
pas ngliger, comme les courants sous-marins o l'coulement du fluide ramen
la surface dans les pipe-lines verticaux.

Fig. 2: Varit de plate-formes ptrolires [3]

12
 Modlisation de linteraction fluide/structure

3 Prsentation du logiciel LS-Dyna

3.1 LS-Dyna

LS-DYNA est un logiciel de calcul par lments finis (FEM) qui permet
d'analyser des comportements physiques hautement non-linaires avec des
mthodes numriques dites explicites ou implicites. On trouve typiquement des
applications avec de trs grandes dformations dans un laps de temps trs court
comme par exemple un crash test. LS-DYNA s'est impos au fil du temps
comme un outil de calcul trs puissant, ceci grce ses trs nombreuses
fonctionnalits. Il permet le calcul de problmes 2D ou 3D des plus simples aux
plus complexes. LS-DYNA compte de nombreux utilisateurs dans les domaines
de l'automobile, l'aronautique, le gnie maritime et civil.
Ces fonctionnalits donnent l'utilisateur la possibilit de calculer un large
ventail de cas comme la simulation de crash test, de la scurit des passagers
(industrie automobile, navale, aronautique, ferroviaire), du formage des mtaux
et du calcul de chute (industrie de biens de consommation). LS-DYNA permet
galement le calcul non-linaire de flambage ou de rsistance sismique.
Le calcul coupl fluide-structure permet de simuler l'interaction entre une
structure et un fluide comme l'aquaplaning, la chute d'un rcipient rempli de
fluide ou leffort de la houle sur une structure offshore.
LS-DYNA peut tre utilis sur pratiquement toutes les plates-
formes: UNIX (SGI, HP, Compaq, SUN, CRAY, etc.), Windows PC et LINUX.
LS-DYNA peut exploiter les performances d'un ordinateur multi-processeur
(SMP; Shared Memory) ou d'un cluster (MPP; Massively Parallel), il est
compltement vectoris.

Fig. 3: organigramme LS-Dyna [7]

13
 Modlisation de linteraction fluide/structure

3.2 Organisation des commandes de programmation (keyword)

La figure ci-dessous relate la philosophie gnrale, lorganisation et la

relation entre les diffrentes entits de programmation (keyword). Par exemple
les donnes inclues par le keyword ELEMENT, sont la fois : son identit (EID),
lappartenance de llment une partie (PID) et lidentit des nuds qui le
composent (N1, N2, N3, N4). Les coordonnes cartsiennes de chaques nuds
sont dfinies par le keyword NODE. Aussi les diffrentes parties du systme sont
dfinies par le keyword PART. On attribue ainsi chaques parties son identit
(PID), lidentit de sa section (SID), lidentit du matriau qui la compose (MID)
ou encore lidentit de lquation dtat qui la gouverne (EOSID)
Cest seulement par la suite que lon dfinit rellement les
caractristiques physiques et mathmatiques du systme. Ainsi avec le keyword
SECTION_SCHELL on dfinit les rgles dintgrations, la formulation de
llment, son paisseur... Le keyword EOS permet de dfinir lquation dtat.

Fig. 4 Organisation des commandes de programmation [7]

Identit des keywords

NID = identit du nud, EID = identit de llment

Dfinition des keywords

Par exemple : *NODE
10101, 0, 0, 0
Le nud 10101 se trouve aux coordonnes 0, 0, 0.

Et ainsi de suite on dfinit les diffrents keywords du programme, comme la

section, le matriau, lquation dtat

14
 Modlisation de linteraction fluide/structure

3.3 Fonctionnement thorique de LS-DYNA

A partir des quations de Navier Stokes [annexe A], conservation de la

masse, conservation de quantit de mouvement, conservation de lnergie,
lquation constitutive dun fluide newtonien et les conditions aux limites [forme
forte], on obtient la forme faible qui est intgrable.
Pour mettre jour les variables partir des quations prcdentes, le pas
de temps est divis en deux cycles de calcul: un cycle Lagrangien qui dtermine
les inconnues physiques du problme et un cycle dadvection qui ramne la
solution physique sur une configuration ALE (Arbitraire Lagrange Euler) mobile
(maillage mobile de vitesse arbitraire vale). Le calcul Lagrangien implique que le
maillage se dforme avec le fluide et peut se distordre. La seconde tape permet
dviter cette distorsion en dplaant le maillage de manire prserver la
rgularit de ce dernier. Le mouvement du maillage est indpendant du
mouvement du fluide ce qui entrane des flux de masse, de quantit de
mouvement et dnergie.

Cycle Lagrangien

1) Connaissant le champ de contraintes, de pression, les forces de Hourglass et

la viscosit de choc linstant tn dans chaques lments, les forces aux nuds
sont calcules. On en dduit les acclrations en divisant ces forces par la
masse nodale.
2) Les acclrations sont intgres pour donner les vitesses linstant t(n+1)/2.
3) Les vitesses sont intgres leur tour pour fournir le champ de dplacements
linstant t n+1 mettant ainsi jour la position Lagrangienne du maillage.
4) La viscosit de choc artificielle et la viscosit de Hourglass sont calcules
partir de v(n+1)/2
5) Lnergie interne est mise jour.
6) La pression est dtermine partir de lquation dtat en utilisant la masse
volumique et lnergie interne par unit de volume.
7) Un nouveau pas de temps est calcul

Cycle dadvection (ou Eulerienne)

Remarque: le cycle Eulerien projette, sur le nouveau maillage, la solution

physique dtermine dans le cycle Lagrangien en procdant par les sous-tapes
suivantes:

8) Un nouveau maillage est gnr en relaxant le maillage Lagrangien dans un

calcul ALE.
9) Les variables centres dans llment (contrainte, nergie interne,...) sont
projetes par un algorithme dadvection de type Van Leer.
10) Les vitesses aux nuds sont projetes leur tour sur le nouveau maillage
Un nouveau maillage, dcal dun demi-pas dans les 3 directions, de manire
ce que les centres des mailles soient les nuds du maillage original, est cre.

15
 Modlisation de linteraction fluide/structure

Schma dexplication

n tn+1
t

1,2 Fig. 5 cycle Lagrangien 4,5,6,7

tn+1 tn+1

8,9,10

Fig. 6 cycle Eulerien ou ALE

16
 Modlisation de linteraction fluide/structure

4 Modlisation du systme

4.1 Modlisation de linteraction fluide/structure

Le couplage fluide-structure apparat systmatiquement ds que lon a

une structure dformable dans un fluide en mouvement ou une structure en
mouvement dans un fluide initialement au repos. Il est d lacclration des
particules du fluide (effet dinertie) ou encore aux effets visqueux dpendants de
la vitesse du fluide au voisinage de la structure. La variation de pression qui en
rsulte induit sur la structure une force dinertie qui agit comme une masse
ajoute fixe au systme mcanique et se dplaant solidairement avec lui. En
prsence dun fluide visqueux, un amortissement est induit par la viscosit et
sajoute lamortissement propre de la structure.
Avec le dveloppement des codes de calcul en mcanique des fluides et
laugmentation des ressources informatiques, la simulation numrique devient
une alternative intressante et complmentaire pour ltude des couplages
fluide-structure. En effet on peut envisager de reproduire numriquement les
couplages fluide-structure tudis exprimentalement et analytiquement. On est
alors confront aux problmatiques suivantes : assurer la conservation de
lnergie lors de la transmission des donnes entre les codes et grer linterface
mobile o schangent les donnes et o se cre le couplage. Les couplages
fluide-structure soprent au niveau de linterface fluide-structure. En thorie on
suppose que lon a continuit des contraintes et des vitesses linterface entre le
fluide et la structure. La difficult est de trouver une formulation adapte pour la
rsolution des quations fluide-structure qui tient compte des parois mobiles et
permet de simuler des coulements fortement convectifs. La formulation
Lagrangien-Eulerien Arbitraire (ALE) est une solution qui permet de grer des
maillages mobiles avec des dplacements.

4.2 Lagrangien-Eulerien Arbitraire (ALE)

4.2.1 Formulation Eulerienne

La description Eulerienne tudie le passage dune certaine quantit de

matire travers une rgion donne de lespace. Le systme de rfrence
associ est fixe et est appel Systme de Rfrence Spatial (SRS). Dans cette
description, lobservateur se place en un point fixe et tudie comment varient les
caractristiques de la particule passant par ce point linstant considr. Cette
formulation est gnralement celle qui est adopte en mcanique des fluides.
Dans le cadre de la mthode des lments finis, la formulation Eulerienne
impose un maillage fixe dans lespace. Ce maillage dfinit un volume de contrle
invariable au cours du temps et une particule matrielle donne peut migrer dun
lment fini vers un autre. Certaines particules matrielles peuvent mme quitter

17
 Modlisation de linteraction fluide/structure

le domaine dfini initialement. Dans ce cas toutes les informations ultrieures

relatives ces particules sont perdues.
Lutilisation de formulations Euleriennes fait apparatre des difficults lors
de ltude de problmes frontires variables. Les frontires du corps tudi ne
concident gnralement ni partout ni chaques instants avec le bord dun
lment. Il en rsulte, de grandes difficults pour prendre en compte les
conditions aux limites matrielles, ainsi quune dfinition peu prcise des
frontires matrielles lorsque celles-ci voluent au cours du temps.
En revanche, les formulations Euleriennes autorisent dimportantes
distorsions de la matire. Ceci explique pourquoi elles sont largement utilises
en mcanique des fluides ou dans le cas de la mise en forme des matriaux pour
ltude des processus stationnaires.

4.2.2 Formulation Lagrangienne

La description Lagrangienne tudie le mouvement du corps considr

comme une certaine quantit de matire se dplaant dans lespace partir
dune rfrence connue. Le systme de rfrence, appel Systme de
Rfrence Matriel (SRM), est donc attach la matire. Dans ce cas,
lobservateur peut tre considr comme attach une particule matrielle
donne et il tudie lvolution temporelle des proprits de cette particule au
cours de son mouvement. Cette formulation est classiquement adopte en
mcanique des solides.
Dans le cadre de la mthode des lments finis, la formulation
Lagrangienne impose chaques nuds de rester en concidence avec la mme
particule matrielle tout au long du processus de dformation. Cette proprit est
intressante dun point de vue numrique puisque, lors de lintgration des lois
de comportement, les termes convectifs napparaissent pas dans les quations.
Dautre part, un lment donn contient toujours le mme ensemble de
particules matrielles et est donc entran par la dformation. En particulier, un
nud situ la frontire du corps reste sur cette frontire au cours de la
simulation et, dans ce cas, le suivi dune frontire mobile ne prsente pas de
difficults.
En revanche, en prsence de grandes dformations, le maillage peut tre
fortement distordu et entraner une imprcision des rsultats voire limpossibilit
de conduire la simulation terme. Des oprations de ractualisation de maillage
(remaillage) sont alors indispensables, oprations coteuses en temps (homme
et machine).

4.2.3 Formulation ALE

4.2.3.1 Introduction

La formulation ALE est particulirement bien adapte ltude des

problmes prsentant de grandes dformations et plus gnralement ltude

18
 Modlisation de linteraction fluide/structure

des problmes coupls. Elle permet de combiner les avantages de la

reprsentation Lagrangienne (bonne dfinition des frontires variables, facilit
dimposer des conditions aux limites matrielles) et ceux de la formulation
Eulerienne (possibilit de prise en compte de grandes distorsions). Le principe
de cette formulation mixte est le dveloppement de maillages se dplaant
taux indpendant (bien que non nul) de celui de la matire. En consquence, la
matire se dplace par rapport au maillage, comme dans une formulation
Eulerienne pure, mais la forme extrieure des lments est contrle par les
conditions aux limites du problme, comme dans une formulation Lagrangienne.

Fig. 7: volution du maillage

La figure illustre les difficults associes aux formulations Eulerienne et

Lagrangienne ainsi que lintrt des formulations de type ALE.

4.2.3.2 Mise en uvre des formulations ALE

Ce paragraphe prsente les tapes principales de la mise en uvre des

formulations ALE. Pour plus de dtails, il convient de se reporter, par exemple,
aux travaux de Ponthot et de Souli.

Relation fondamentale

Pour la formulation ALE, trois configurations sont distingues :

- La configuration initiale C0 o chaque particule est associe un nud.
- La configuration Ct. Cest la configuration aprs volution du domaine
(dformations, mouvement) pour laquelle lquilibre est recherch.
- La configuration de rfrence Cref . Cest dans cette configuration que les
nuds sont positionns et associs des particules.

Pour la formulation Lagrangienne, Cref et C0 sont identiques alors que pour

la formulation Eulerienne, ce sont Cref et Ct qui sont identiques. En ce qui
concerne la formulation ALE, Cref est choisie de la faon la plus judicieuse
possible. Dans ce cas, lors du processus dvolution, la particule est transporte
par une application alors que le nud est transport par une application .
Dans la configuration Ct, le nud nest alors plus associ la mme particule
r
que dans C0. Nous avons donc les variables Euleriennes ( x , t ) qui sont

19
 Modlisation de linteraction fluide/structure

r
exprimes, en fonction des coordonnes de la configuration initiale ( X , t ) tel que
r r
x = ( X , t ) et en fonction des coordonnes de la configuration de rfrence
r r r
( , t ) tel que x = ( , t ) . Ainsi donne le mouvement du materiau et
reprsente le dplacement du domaine ALE.

Fig. 8: notation en fonction de la formulation

Lors de la mise en quation, on est amen calculer la drive par

rapport au temps de champs tels que ceux relatifs la vitesse, la masse
volumique. Le choix du rfrentiel conditionne alors la forme des expressions
des drives. Soit f une fonction reprsentant un champ. Lorsque le systme de
rfrence est le systme de rfrence spatial (formulation Eulerienne), la
fonction est note f ( x, t ) . Si le systme de rfrence est matriel (formulation
~
Lagrangienne), la fonction est note f ( X , t ) . Enfin, pour un systme de grille
(formulation ALE), la fonction est note f ( , t )
il vient donc :
~
f ( x, t ) = f ( X , t ) = f ( , t )
Eq: 4-1

En calculant les drives de chaque fonction par rapport au temps, on obtient :

- Dans le cas dune formulation Eulerienne :

Eq: 4-1

- Dans le cas dune formulation Lagrangienne :

Eq: 4-2
- Dans le cas dune formulation ALE :

Eq: 4-3

20
 Modlisation de linteraction fluide/structure

O :
f est le gradient de la fonction
f est le gradient de la fonction f
u : Vitesse de la particule par rapport au systme de rfrence spatial
w : Vitesse du maillage par rapport au systme de rfrence spatial
c : Vitesse de la particule par rapport au systme de rfrence de grille (vitesse
convective)
Lexpression (ALE) permet de retrouver le cas Eulerien () et le cas
Lagrangien (u = w).
En combinant les quations Lagrangiennes, Euleriennes et ALE, on
obtient la relation suivante,

Eq: 4-4

r
x f
r r
x = ( , t)
r r
x = ( X ,t)

r
f
r
X f
r r
= o 1 ( , t )
Fig. 9 fluide dans les trois configurations

r
x s

r r
x = ( X ,t)

r
X s

Fig. 10 structure dans la configuration initiale et courante

21
 Modlisation de linteraction fluide/structure

4.3 Couplage fluide/structure

4.3.1 Gnralits

On distingue gnralement trois types d'approches pour modliser le

couplage :
-l'approche simplifie consiste essayer de dcoupler les quations de
chaque milieu en exprimant un systme en fonction de lautre. Dans certains cas,
il est possible de trouver une expression analytique reprsentant laction dun
systme sur lautre. Ainsi, les forces arodynamiques s'exerant sur un profil
d'aile peuvent sous certaines hypothses simplificatrices tre exprimes en
fonction de la dformation du profil. Cette pratique a t couramment employe,
notamment pour comprendre le flottement des ailes d'avion.
-l'approche dcale ou partitionne permet de rsoudre l'aide de codes
de calculs ddis la physique du fluide d'une part et du solide d'autre part. Pour
chaque milieu, un code de calcul spcifique est utilis et la difficult consiste
faire transiter l'information d'un code l'autre. Par ailleurs, les codes ne
fonctionnent pas simultanment mais successivement, ce qui introduit une erreur
numrique supplmentaire qui peut tre fatidique pour la simulation.
-l'approche monolithique reprsente thoriquement la solution optimale
puisque le fluide et le solide sont traits par un mme code de calcul. Toutefois,
lorsque les gomtries ou la physique du problme traiter deviennent
complexes ce type de mthode n'est plus envisageable puisque chaque milieu
(fluide ou solide) ncessite des procdures de calcul numrique spcifiques.

Chaque mthode prsente des avantages et des inconvnients, mais

dans la pratique, l'approche dcale est la plupart du temps utilise. En effet,
l'approche simplifie est insuffisante dans la majorit des cas puisque la
rtroaction d'un milieu sur l'autre est mal prise en compte. Le dernier type
d'approche est en gnral peu envisageable pour traiter des configurations
complexes, d'autant plus que le formalisme employ pour modliser fluides et
solides diffre gnralement : la formulation Lagrangienne et prfrentiellement
employe pour dcrire les solides, tandis que les fluides sont plutt traits dans
une approche Eulrienne. L'approche dcale constitue donc une bonne
alternative, d'autant plus qu'elle permet de rutiliser les procdures de calcul
spcifiques chaque milieu. La difficult rside dans le traitement de l'interaction
qui ncessite de dvelopper des algorithmes de couplage performants.

4.3.2 Procdure du couplage

Laction de couplage fluide/structure doit ncessairement commencer par

une phase de mise en place dalgorithme de couplage entre trois oprateurs
principaux : un solveur de lcoulement du fluide, un solveur pour la structure et
une procdure de mise jour du maillage fluide.

22
 Modlisation de linteraction fluide/structure

4.3.3 Conditions et algorithme de couplage

Le couplage entre le fluide et la structure intervient au niveau de linterface

entre les deux domaines. Les conditions du couplage sont donc imposes sur la
frontire c. Soit c la frontire dinterface entre le fluide et la structure tel
que c = s f . Les conditions de couplage consistent en une opration de
transfert de champs entre le fluide et la structure, et concerne le transfert des
contraintes de pression au solide, qui en se dformant transfert les champs de
dplacements au fluide.

Ces conditions doivent vrifier.

-la continuit des vitesses
-la continuit des efforts

Les tapes de transfert des donnes se droulent de la manire suivante :

1. le code fluide transmet les pressions, nodales de linterface, au code
structure.
2. le code structure calcule la solution et remet jour les coordonnes des
nuds.
3. la structure transmet les nouvelles positions des nuds de linterface
4. le code fluide met jour le maillage, en tenant compte des nouvelles
positions des nuds, et calcule la solution (pressions)
5. retour ltape 1.

Fig. 11 transfert des donnes entre deux codes

23
 Modlisation de linteraction fluide/structure

Fig. 12 algorithme du coupleur

La procdure de couplage est base sur un algorithme explicite. A chaque

itration tn, aprs convergence du code fluide, on obtient la pression
correspondante Pn, qui sera transmise au code structure. Pn sera considre
comme une sollicitation extrieure aux frontires de la structure.
Aprs convergence du code structure, on extrait les nouvelles positions xn
et yn, puis on actualise la gomtrie pour un calcul du code fluide.

4.4 Couplage par pnalit

Le contact entre le fluide et la structure impose en premier lieu des

conditions de vitesses selon la normale la surface de contact qui traduisent la
non-pntration du fluide dans la structure. Dautre part comme les frottements
du fluide sur la structure sont trs faibles pour un fluide rel, les contraintes
tangentielles dans le fluide et sur la structure au niveau de la zone de contact
sont ngliges.
La condition de non pntration est traduite approximativement en
introduisant entre le nud et llment en contact un ressort caractris par sa
raideur, afin de repousser les nuds qui viendraient en contact, grce la
cration defforts fictifs de pnalits. Cette dernire peut tre constante ou
dpendre de la pntration estime du nud (quivalente un amortissement
visqueux), selon les algorithmes de calcul utilises. La traverse des lments
par les nuds peut tre anticipe par limposition dun intervalle (dit gap) de
dtection de contact, ce qui vite des variations trop brutales des forces
appliques aux nuds

24
 Modlisation de linteraction fluide/structure

Fig. 13 schma du contact par pnalisation

Dans tous les cas, la gestion du contact faite par le logiciel LS-Dyna se
compose dune phase de recherche des pntrations ventuelles puis dune
phase de rduction des pntrations par la mise en place de ressorts fictifs entre
le nud pntrant et la facette traverse. En outre, afin de respecter le principe
mcanique de laction et de la raction qui sy opposent. Celles-ci sont fonction,
entre autres, de la position du point de contact par rapport aux quatre nuds de
llment impact.

La principale difficult dans le cas dimpact de leau repose sur le fait que
les raideurs des lments du fluide et de la structure sont trs diffrentes, celle
de linterface tant gnralement une moyenne des raideurs des lments en
contact.

4.5 Description des lments finis utiliss dans le modle fluide

structure

4.5.1 Description gomtrique du maillage de la structure

La structure est discrtise par des lments 2D, de types quadrilatres

(quadrangles). La forme rgulire des diffrentes structures qui seront tudies
nimplique pas que nous devons utiliser des lments de type triangle plus
rompu aux contraintes dune structure irrgulire

Fig. 14 lment de discrtisation de la structure

25
 Modlisation de linteraction fluide/structure

Fig. 15 maillage de lune des structures

Le maillage est raffin au niveau de leau sur la structure pour une

meilleure prcision.

4.5.2 Description de la mthode

-Maillage
-Construction de la matrice de raideur de chaque lment [Ke]
-Assemblage de la matrice globale [K]
-Construction du vecteur chargement {F}
-limination de certains degrs de libert (si besoin)
-Dtermination des dplacements: {U}= [K]-1{F}
-Calcul des quantits drives du dplacement

4.5.3 Description du maillage du fluide

Le fluide est discrtis par des lments 3D, de types tridres. Cet
lment sadapte particulirement bien aux distorsions du fluide et son

26
 Modlisation de linteraction fluide/structure

mouvement. Cependant ces lments ne sont pas reconnus pour tre les
meilleurs dun point de vue de la prcision des rsultats.

Fig. 16 discrtisation du fluide

4.6 Hourglass [annexe C]

4.6.1 Elment sous intgrs

Pour parvenir des performances maximales sur les temps de calculs, les
lments disponibles dans les codes de calcul comme LS-Dyna sont dits sous
intgrs. Cest dire quils possdent bien moins de points de Gauss que de
nuds

Intgration complte: 4 points

dans le plan milieu

Sous-intgration: 1 point
dans le plan milieu
Fig. 17 lment sous intgr

4.6.2 Mode Hourglass

Pour les lments briques comme pour les lments quadrangles, le

dsavantage dune intgration avec un seul point dintgration est la ncessit

27
 Modlisation de linteraction fluide/structure

de contrler les modes nergie nulle appels aussi modes de Hourglass. Ces
modes impliquent des dformations nulles aux points de quadratures, et lorsque
les vitesses des nuds opposs diagonalement sont gales, alors le taux de
dformations sannule aux point de quadrature.
Ce mode entrane alors un travail nul do le nom de mode nergie nulle.
Linterprtation gomtrique de ce phnomne est que le plan parallle au
gradient de vitesses passant par le point de quadrature ne se dforme pas dans
un mode de Hourglass. La raideur de ce mode est, en gnral, plus petite que la
plus basse des raideurs physiques et ainsi la frquence associe est aussi plus
petite que la plus petite frquence du problme physique. On observe alors des
oscillations perturbant les grandeurs physiques. De plus, le Hourglass est une
instabilit spatiale visible dans la dformation des mailles lors dun calcul
Lagrangien.
Une technique pour contrer ces modes indsirables est lutilisation de
viscosit artificielle ou encore lutilisation dun crouissage lastique capable de
stopper la formation de ces modes de dformation mais ayant un effet
ngligeable sur les modes globaux de dformations stables.

28
 Modlisation de linteraction fluide/structure

5 Vrification de la modlisation du fluide

(hydrostatique)

5.1 Equation dtat

Les quations dtats des diffrents fluides sont des polynmes de

lnergie interne (E) par unit de volume et du taux relatif de la fraction
ref
volumique = o ref est la masse volumique de rfrence. La
ref
pression dpend de lnergie par la relation dtat suivante.
E
P = A1 ( ) + A2 ( )
V
Eq. 5-1
o A1 ( ) = C 0 + C1 + C 2 2
Eq. 5-2
et A2 ( ) = C 3 + C 4 + C 5 + C 6
2 3

Eq. 5-3

o Ci sont des constantes.

Remarque : est frquemment utilis comme une mesure de la compression

(>0), ou de lexpansion (<0).

Dans un fluide, le coefficient de compressibilit K et la densit

permettent de dterminer la vitesse du son c dans ce fluide selon la formule
K
c=

Les quations dtats pour leau et pour lair sont respectivement dfinis
par :

- P = C1 avec Ci=0 pour i=0,26 et C1=ref.c2 =2.2GPa (C1=K). Ce

coefficient augmente avec la pression)

Remarque : Lorsque lon considre un fluide incompressible cela signifie que le

coefficient de compressibilit est trs grand, il est donc trs difficile dapporter
suffisamment dnergie pour comprimer le fluide et avoir une variation de
volume notoire.

-P= ( 1) E avec Ci=0 pour i=0,,3,6 et C4= C5= 1 o est le
ref
rapport des chaleurs spcifiques.

29
 Modlisation de linteraction fluide/structure

Il suffit donc de dterminer E et V linstant initial (t=0) pour pouvoir

dfinir compltement la pression au cours du temps
-E0= ref.Cv.T lnergie interne initial par volume de rfrence
-V0=1 volume relatif initial

Dfinition de leau

ref
ref.c2

Fig. 18: programmation de lquation dtat de leau

E0 V0

5.2 Densit du fluide

La simulation consiste vrifier ici les diffrentes caractristiques

physiques de leau. Ainsi dans un premier temps nous vrifierons que la densit
de leau est bien celle souhait, savoir 998,21kg/m3 20C.
Soit de leau sur 55cm lche dans une cuve remplie de vide initialement.
Daprs les figures ci-dessous, on remarque que lon a bien une densit
de fluide gale la valeur rfrence prcdemment

30
 Modlisation de linteraction fluide/structure

vide

eau

Fig. 19 maillage du systme

Fig. 20 Contour de la densit du fluide

Fig. 21: iso surfaces de la densit du fluide

31
 Modlisation de linteraction fluide/structure

5.3 Pression hydrostatique

Nous vrifions ici que la pression hydrostatique obtenue sur le fond de la

cuve deau est bien identique la thorie, savoir :
Pfond = * g * h = 998,21 * 9,81 * 0,55 == 5,38kPa
Eq. 5-4

Plaque rigide (lment gomtrique)

Eau: h=0,55m

Fig. 22 modle pression hydrostatique

Nous testons diffrents types de chargement de la gravit sur leau. Dans

un premier test nous lappliquerons de faon instantane. Par la suite nous
appliquerons un chargement de la gravit pousant les contours dune sinusode
plus ou moins directe (voir figure ci-dessous).
diffrents types de chargements

12

10 chargement 1s
chargement 2s
chargement 0,5s
chargement instantan
8
gravit [m.s ]
-2

0
0 0,5 1 1,5 2 2,5 3 3,5 4
temps [s]

Fig. 23 diffrents types de chargements

32
 Modlisation de linteraction fluide/structure

Remarque : Pour des questions de fiabilits et de prcision, nous avons somm

les pressions hydrostatiques sur 15 lments du fond de la cuve comme nous le
montrent les figures ci-dessous.

Fig. 24 points de mesure de la pression hydrostatique

Les rsultats obtenus devront donc tre compars 15x5,38kPa soit

80,76kPa.

Test 1 : gravit instantane

Rsultat obtenu incohrent, une

application instantane de la gravit
nest pas envisageable. En effet
nous nobtenons pas les 80,76kPa
thoriques escompts. Leffet
dynamique dun tel chargement est
trop important.

Fig. 25 test gravit instantane

33
 Modlisation de linteraction fluide/structure

Test 2 : gravit applique suivant une courbe sinusodale, mais avec un temps
de chargement variant de 0,5s 2s

Avec un temps de chargement de

0,5s on remarque que ce dernier
nest pas totalement complt. Nous
obtenons une moyenne de pression
quivalente 75kPa environ. De
plus nous avons un effet dynamique
important car la pression nest pas
stabilise. La valeur de la pression
hydrostatique oscille autour dune
valeur proche de 6,25% de la valeur
thorique.

Fig. 26: test temps de chargement 0,5s

Avec un temps de chargement de 1s on

obtient une pression hydrostatique
proche de la thorie avec des effets
dynamiques limits. En effet on obtient
une valeur moyenne de 78kPa, soit un
rsultat proche 2,5% de la thorie

Fig. 27: test temps de chargement 1s

34
 Modlisation de linteraction fluide/structure

Rsultat similaire au prcdent.

Fig. 28: test temps de chargement 2s

Conclusion : Une application instantane de la gravit nest pas envisageable

dans notre cas dtude. En effet les rsultats obtenus divergent de la valeur
souhaite. Aussi un temps de chargement trop court engendrera des effets
dynamiques indsirables, et ne permettra pas au chargement de parfaitement se
drouler. De mme un temps de chargement trop long nest pas souhaitable car
il augmente le temps de calcul. En effet, notre souhait sera de librer leau une
fois que son chargement est complt, donc il faut trouver le meilleur compromis
possible entre stabilit et rapidit.

Cest dans cet tat desprit que nous allons tester le cycle dadvection.

5.4 Test de ladvection

Le cycle dadvection permet de projeter les variables centres dans un

lment sur le nouveau maillage ALE. Ce cycle ce produit aprs le cycle
Lagrangien dcrit prcdemment, et en fonction de la tendance du maillage
Lagrangien ce distordre on le remet plus ou moins souvent jour (maillage
ALE).
Les rsultats obtenus dans le paragraphe prcdent lon t pour une
thorie de Van Leer [Annexe C] du second ordre remis jour aprs chaque cycle
Lagrangien. On va ainsi tester sparment linfluence de la remise jour (tous
les 3 cycles) et de la mthode employe (mthode donor cells, littralement
cellule donneuse, qui est de type du premier ordre).

35
 Modlisation de linteraction fluide/structure

Test 1 : remise jour du maillage tous les 3 cycles (temps de chargement 2s)

Avec un nombre de 3 cycles entre

chaque advection on obtient un
rsultat semblable, proche du
prcdent mais avec un gain de
temps de 46 %. Cependant ce
rsultat est relativiser tant
donn que nous sommes dans une
phase statique. Cela nous permet
nanmoins de nous sensibiliser
la notion de temps dans les
simulations numriques.

Fig. 29: test advection mise jour tous les 3 cycles de calcul

Test 2 : mthode dadvection du premier ordre

Avec une mthode dadvection du

premier ordre on retrouve les mmes
rsultats quobtenus prcdemment.
Une mthode du premier ordre peut
donc convenir mais mme remarque
que prcdemment, nous sommes
dans un cas particulier o leau est
statique.

Fig. 30: test advection du premier ordre

36
 Modlisation de linteraction fluide/structure

Conclusion : Nous pouvons choisir de rduire le temps de calcul en augmentant

le nombre de cycles entre chaque advection ou encore en changeant le type
dadvection, mais chaque cas dtude dtient sa propre vrit.

5.5 Efforts de leau

Nous vrifions ici que leffort de leau sur la plaque engendr par la
gravit et la densit de leau est gale la valeur thorique, savoir :

b h
F = * g * zdydz
0 0
Eq. 5-5
b h
F = dy * g * zdz
0 0
Eq. 5-6
h
1
F = [ y ]0 * * gz 2
b

2 0
Eq. 5-7
1
F = * gbh 2
2
Eq. 5-8

Avec h=0,55m et b=0,15m on obtient :

1
F= * 998,21 * 9,81 * 0,15 * 0,55 2 = 222,17 N
2
Eq. 5-9

Plaque: matriau lasto-plastique,

encastre sur ses bords

Cuve deau de 0,55m*0,15m

Fig. 31 modle effort hydrostatique

37
 Modlisation de linteraction fluide/structure

Pour vrifier que leffort de leau est reproduit sur la plaque nous allons
sommer lensemble des efforts selon laxe x sur ses limites extrieures.

Points de mesures : les nuds

des bords encastrs

Fig. 32 pression hydrostatique

5.5.1 Rsultat

Fig. 33 Somme des forces sur le contour de la plaque

Nous obtenons bien une valeur proche des 222,17kN thoriques,

cependant nous allons maintenant nous intresser la rpartition des
contraintes sur la plaque. En effet la somme des forces ne nous donne quune
valeur sur lintensit mais ne nous garantie pas que le modle de fluide
dvelopp agit comme leau.

38
 Modlisation de linteraction fluide/structure

Plaque dpaisseur 1mm

et de largeur 0,35m

0,55m

0,175m

Fig. 34 contrainte du leau suivant x

Nous remarquons que la zone la plus sollicite se trouve au 2/3 de la

surface libre de leau. Nous avons une contrainte de 14,75kPa ce niveau.
Cette zone correspond au centre de gravit du triangle de la pression
hydrostatique. Visuellement il apparaitrait donc que la rpartition des contraintes
soit cohrente.

Fig. 35 Somme des contraintes de la plaque suivant x

Thoriquement la somme des contraintes suivant x est gale :

xx = F / S xx
xx = 222,17 /(0,001 * 0,35) = 0,63MPa
Eq. 5-10

39
 Modlisation de linteraction fluide/structure

Daprs les courbes prcdentes nous avons une somme des contraintes
de 0,65MPa, soit un cart relatif de 3,18% avec la valeur thorique. On peut
donc en conclure que la rpartition et lintensit de leffort de leau sur la plaque
soient cohrentes avec la thorie.

5.5.2 Conclusion

Nous pouvons conclure que le modle de fluide ici dvelopp et test

regroupe les caractristiques physiques de leau. A savoir sa densit, pression
hydrostatique, rpartition des efforts
Aussi le couplage entre le fluide (leau) et la structure (plaque, lasto-
plastique) semble fonctionner. Pour tester plus en profondeur si le couplage
fluide/structure est cohrent nous allons vrifier le comportement de la plaque
grce la thorie des plaques de Timoshenko.

40
 Modlisation de linteraction fluide/structure

6 Vrification du comportement de la structure (cas

statique)
Nous conservons le mme modle que celui de la figure : Fig. 31 modle
effort .

6.1 Dfinition de la plaque

Fig. 36 programmation du matriau de la structure

Dfinition de lpaisseur (t) de la plaque

Limite de plasticit
[Pa]

Dfinition de la courbe du Coefficient de Poisson

matriau
Module dYoung [Pa]

Masse volumique [kg/m3]

Nous avons donc dfini une plaque en acier S235, de masse volumique
7850kg/m3, de module dYoung (E).210000MPa, de coefficient de poisson 0,30
et donc de limite de plasticit 355MPa.

41
 Modlisation de linteraction fluide/structure

6.2 Courbe du matriau

proprit du matriau

400
350
contrainte [Mpa]

300
250
200
150
100
50
0
0 0,02 0,04 0,06 0,08 0,1 0,12 0,14 0,16
dplacement relatif
Fig. 37 courbe du matriau

6.3 Thorie des plaques

6.3.1 Introduction

La capacit de flexion dune plaque dpend davantage de son paisseur

que de ses autres dimensions. Aussi nous pouvons distinguer trois types de
plaques : (1) les plaques minces avec de petits dplacements, (2) des plaques
minces avec de grands dplacements, (3) les plaques paisses.
Lorsque que lon parle de petit ou de grand dplacement, on sous entend
des dplacements petit ou grand en comparaison de lpaisseur.

6.3.1.1 Dmarche

Fig. 38 dmarche de la thorie des plaques

42
 Modlisation de linteraction fluide/structure

6.3.1.2 Hypothses

La thorie des plaques minces suppose que :

-le plan moyen (quivalent de la courbe moyenne des poutres) est
initialement plan ;
-le feuillet moyen (quivalent de la fibre neutre des poutres) ne subit pas
de dformation dans son plan ; on ne considre que le dplacement transversal
w des points du feuillet moyen ;
-les sections normales au feuillet moyen restent normales lors de la
dformation ; en consquence, on peut ngliger le cisaillement ;
-l'paisseur est faible ; en consquence, les contraintes dans le sens de
l'paisseur sont supposes nulles ;
-on reste en petites dformations.

6.3.2 Equation diffrentielle

-Daprs la thorie des plaques de Timoshenko nous avons :

4w 4w 4w q
+ 2 + =
x 4 x 2 y 2 y 4 D
Eq. 6-1

Avec w la flche
q le chargement
D (voir ci-dessous)

-Introduction du coefficient D

Fig. 39 flexion de plaque

43
 Modlisation de linteraction fluide/structure

Daprs la loi de Hooke, on peut crire les dplacements relatifs x et y en

fonction des contraintes x et y :

x y
x =
E E
Eq. 6-2
y x
y = =0
E E
Eq. 6-3
Do
(1 2 ) x
x =
E
Eq. 6-4
Et
E x Ez d 2 w
x = =
1 2 1 2 dx 2
Eq. 6-5
Nous obtenons par intgration le moment M :
h h
2 2
Ez d 2 w Eh 3 d 2w
M = x zdz = dz =
h1 12(1 2 ) dx 2
2
h dx 2

2 2
Eq. 6-6
On introduit la notation
Eh 3
=D
12(1 2 )
Eq. 6-7
La quantit D[N.m], prend la place de la quantit EI dans le cas de la
thorie des poutres, et est appele rigidit flexionnelle de la plaque.

6.4 1ere vrification

Remarque : Daprs Eq. 6-1 nous remarquons que dans des cas similaires de
plaques (mmes conditions aux limites) et de chargement, le dplacement sera
proportionnel D, donc lpaisseur au cube.

6.4.1 Modle

Soit le modle de la figure : Fig. 32 pression hydrostatique. Etant donn quil

nexiste aucune bibliographie donnant la flche thorique pour ce type de plaque
et quil serait fastidieux de tenter de la calculer, nous allons seulement nous
rapporter la remarque prcdente. En effet nous allons dans le cas de

44
 Modlisation de linteraction fluide/structure

diffrentes paisseurs, vrifier que la flche de la plaque varie

proportionnellement avec t3 (paisseur de la plaque au cube)

Fig. 40 points de mesure 1ere vrification

Pour obtenir la flche maximale de la plaque, nous faisons la moyenne

des quatre rsultats obtenus aux nuds de llment plac au centre de la zone
soumise la pression hydrostatique.

Nous avons choisi quatre plaques encastres sur leurs contours

extrieurs conformment au modle rappel ci-dessus. Ces quatre plaques ont
des paisseurs diffrentes variant de 2mm 10mm. Voir figure ci-dessous

t=2mm t=3,5mm t=5mm t=10mm

Fig. 41 dplacements des diffrentes plaques dans le 1er cas

45
 Modlisation de linteraction fluide/structure

6.4.2 Rsultats

Fig. 42 graphique proportionnalit flche-paisseur

Les rsultats sont prsents sous la forme dun graphe deux ordonnes.
A droite nous avons la valeur de D en [N.m] et gauche la valeur de la flche en
[m]. Ce type de graphe permet de mettre en vidence la correspondance entre la
valeur de D et la flche, ce qui rejoint la thorie des plaques de Timoshenko.
Aussi cela semble prouver que le couplage entre le fluide et la structure
fonctionne, et valide en partie la modlisation de la structure :
Cependant lorsque lon augmente lpaisseur de la plaque, les effets
dynamiques deviennent plus prpondrants (oscillations de la plaque 10mm). En
effet, tant donn que pour cette plaque la flche est nettement plus faible, les
effets dynamiques du calcul auront plus dimportance.

6.5 2eme vrification

Remarque : Nous allons nous ramener une des tudes menes par
Timoshenko, savoir : une plaque subissant la pression hydrostatique mais
encastre sur le contour de la zone o cette pression agit.. Nous pourrons ainsi
comparer les flches obtenues aux flches thoriques.

46
 Modlisation de linteraction fluide/structure

6.5.1 Modle

Fig. 43 modlisation plaque encastre sur le contour de la pression

vide

b=0,275m

Plaque encastre

a=0,55m

Repre plac au centre de la plaque

q0=5,38kPa
6.5.2 Rsultats thorique

Remarque : Les rsultats thoriques sont exprims dans le cas de la thorie des
plaques minces avec de petits dplacements. On verra plus-tard que celle-ci a
son importance.

Soit le chargement :
q0 q0 z
q0 = +
2 2 a
Eq. 6-8

Nous pouvons dcomposer le chargement en deux parties. Une partie

chargement uniforme et une partie chargement triangulaire. Cela permet
dutiliser les rsultats de Timoshenko dans le cas dune plaque encastre sur
ses quatre cots soumise un chargement surfacique constant, dont le rsultat
est le suivant :
qa 4
La flche au centre est gale w = 0.0014 donc dans notre cas
D
q0 a 4
w1 = 0.0014
2D
Eq. 6-9

47
 Modlisation de linteraction fluide/structure

On sait aussi daprs la thorie des plaques, quune plaque encastre sur
les cots (a) et simplement appuye sur les cots (b), soumise un chargement
q z
triangulaire 0 (=combinaison des rsultats pour ce type de plaque, pour un
2 a
q z
chargement uniforme gale 0 + pour un chargement triangulaire q 0 ), sa
2 a
flche qui est gale

m
+1
2q 0 a 4
(1) 2 2 + m tanh m my 1 my mz
w2 = 5
D m = 2, 4, 6... m 5
(2
cosh m
cosh
a
+
cosh m
sinh
a
) sin
a
Eq. 6-10

Consquemment pour liminer la courbure le long des cots simplement

appuys, on leurs applique un moment (CL= conditions aux limites) tel que :

m 1

my
(M z ) z = a = (1) 2
E m cos
2 1, 3, 5.. b
m 1

mz
(M ) y y= b = (1) 2
Fm cos
2 1, 3, 5.. b
Eq. 6-11

On obtient aprs calcul, dans le cas de calcul o b/a=0.5, une flche au

centre de la plaque de :

w = w1 + w2 (+CL)
q0 a 4
w = 0,000080
D
Eq. 6-12

6.5.3 Rsultats et comparaison

Nous choisissons de travailler sur trois paisseurs de plaques diffrentes,
t=1mm, t=2mm et t=5mm.

Do les rsultats thoriques suivant :

Epaisseur t[mm] D[N.m] w[mm]
1 19,23 2,05
2 153,85 0,256
5 2403,85 0,0164

On remarque ds lors que la plaque dpaisseur 1mm, aurait selon la

thorie des plaques minces avec petit dplacement, une flche en son centre
deux fois suprieure son paisseur. Nous ferons donc une tude particulire
de ce cas par la suite.

48
 Modlisation de linteraction fluide/structure

t=2mm t=5mm

Fig. 44 dplacements des plaques dans le 2eme cas

Fig. 45 comparaison thorie/ls-dyna pour t= 2mm

Fig. 46 comparaison thorie/ls-dyna pour t=5mm

Nous retrouvons des rsultats trs proches de la thorie. Cela confirme la

modlisation de la structure, ainsi que le couplage fluide structure utilis par le

49
 Modlisation de linteraction fluide/structure

logiciel LS-Dyna. La prochaine investigation concernera le cas dun fluide en

mouvement impactant une structure lastique.
Mais avant nous avons vu que la plaque de 1 mm ne semble plus
rpondre aux exigences dune plaque mince avec de petits dplacements.

6.5.4 Plaque mince avec grand dplacement

Timoshenko a tabli un lien entre la thorie des plaques minces avec

petits dplacements et celle avec grands dplacements comme nous le montre
le graphique ci-dessous.

Fig. 47 Correspondance entre les deux thories

En effet les courbes linaires reprsentent la thorie des petits

dplacements, tandis que les courbes incurves reprsentent celle des grands
dplacements. On a vu que si lon considrait la premire thorie pour la plaque
t=1mm on obtenait w/h=2. Or si lon utilise la seconde thorie on obtient pour le
mme cas de figure w/h=1,16

1,16

Fig. 48 rapport w/h

50
 Modlisation de linteraction fluide/structure

Rsultats et comparaison

Epaisseur t[mm] w(petit dep.)[mm] w(grand dep.)[mm]

1 2,05 1,16

Fig. 49 comparaison thorie/ls-dyna pour t=1mm

Nous pouvons donc constater que les rsultats obtenus grce au logiciel
LS-Dyna sont semblables ceux obtenus par la thorie des plaques minces
avec petit ou grand dplacement. Cela confirme nos premires impressions de la
capacit du logiciel rpondre a un problme de couplage fluide/structure dans
le cas particulier de pression hydrostatique.

6.6 Etat de contrainte de la plaque

Remarque : Nous utilisons la modlisation dfinie Fig. 43 modlisation plaque encastre

sur le contour de la pression.

6.6.1 Tenseur des contraintes

Le tenseur des contraintes est une reprsentation utilise en mcanique

des milieux continus pour reprsenter l'tat de contrainte, c'est--dire les forces
surfaciques (parfois appeles efforts) mises en jeu entre les portions dformes
du milieu. Le terme a t introduit par Cauchy vers 1822.

51
 Modlisation de linteraction fluide/structure

Comme les forces surfaciques sont dfinies pour chaque surface coupant
le milieu, le tenseur est dfini en chaque point du solide, ou localement. L'tat de
contrainte du solide est donc reprsent par un champ de tenseur. On parle
aussi de ce fait de champ de contrainte.
Dans le cadre de l'lasticit linaire, le champ de contrainte est reli au
champ de dformation par la loi de Hooke gnralise, c'est--dire que l'on peut
crire l'quation tensorielle (et non algbrique) = E.

XX XY XZ

T = YX YY YZ
ZY ZZ
ZX

Lorsque le solide est l'quilibre (les contraintes et dformations restent

constantes dans le temps, il n'y a pas d'onde) et si l'on peut ngliger les forces
volumiques (en particulier le poids), alors le tenseur des contraintes est
symtrique, c'est--dire que : ij = ji

Fig. 50 points de mesure des contraintes

6.6.2 Contraintes principales

On peut trouver une base de l'espace tel que le tenseur des contraintes
soit une matrice diagonale :

1 0 0

T = 0 2 0
0 0 3

52
 Modlisation de linteraction fluide/structure

Cependant nous sommes dans un cas particulier ou lon peut considrer

que nous sommes en contraintes planes, dans ce cas on peut crire le tenseur

0
T = 1
0 2

Daprs la reprsentation du cercle de mohr on retrouve rapidement les

formules nous donnant les diffrentes valeurs des contraintes principales.

(ZZ , YZ)

3 2 1

(YY , YZ)

Fig. 51 cercle de mohr

A partir des contraintes donnes par LS-Dyna nous allons comparer les
rsultats donns par les formules thoriques pour les contraintes principales et
les valeurs de celles-ci donnes par LS-Dyna.

1ere contrainte principale

-Formule 1 : hypothses contraintes planes

+ ZZ 1
1 = YY + ( ZZ YY ) 2 + 4 YZ
2 2
Eq. 6-13
-Formule 2 :
XX + ZZ 1
1 = + ( ZZ XX ) 2 + 4 XZ
2 2
Eq. 6-14

53
 Modlisation de linteraction fluide/structure

-Comparaison

Fig. 52 : 1ere contrainte principale

2eme contrainte principale

-Formule 1 : hypothses contraintes planes

+ ZZ 1
2 = YY + ( YY ZZ ) 2 + 4 YZ
2 2
Eq. 6-15
-Formule 2 :
YY + XX 1
2 = + ( YY XX ) 2 + 4 YX
2 2
Eq. 6-16
-Comparaison

Fig. 53 : 2eme contrainte principale

54
 Modlisation de linteraction fluide/structure

3eme contrainte principale

-Formule 2 :
XX + YY 1
3 = + ( XX YY ) 2 + 4 XY
2 2
Eq. 6-17
-Formule 2 :
XX + ZZ 1
3 = + ( XX ZZ ) 2 + 4 ZX
2 2
Eq. 6-18
-Comparaison

Fig. 54 :3eme contrainte principale

Conclusion

Nous pouvons remarquer des diffrences notables entre les valeurs que
lon est sens obtenir selon la thorie, et les valeurs obtenues par LS-Dyna. Tout
en sachant que les valeurs de bases sont les mmes, savoir les contraintes
donnes par LS-Dyna.
Cependant les valeurs du logiciel se rapprochent plus des valeurs
donnes pour des contraintes planes. Cela semble indiquer que les contraintes
suivant lpaisseur de la plaque sont peu fiables (voir Fig. 54 :3eme contrainte
principale). Mais elles ne sont pas responsables des diffrences car ngligeables
par rapport aux autres contraintes. Et cest contradictoire avec la vrification des
contraintes effectues prcdemment.
La diffrence entre la thorie et LS-Dyna est aussi visible pour les
contraintes de Von-Mises, Tresca (voir graphe ANNEXE F)
Il est cependant difficile de conclure quoi que se soit, cependant si les
contraintes releves sont bonnes cela signifierait un disfonctionnement de LS-
Dyna dans le calcul des contraintes principales et autres.

55
 Modlisation de linteraction fluide/structure

6.6.3 Pression isostatique

La pression isostatique p est dfinie comme le tiers de la trace de la

matrice, c'est--dire comme la moyenne des termes diagonaux
Sa valeur ne dpend pas du repre x y z utilis, c'est une gnralisation
de la notion de pression hydrostatique dans les liquides.

-Formule 1 (avec les contraintes donnes par LS-Dyna)

XX + YY + ZZ
p=
3
Eq. 6-19
-Formule 2 (avec les contraintes principales calcules de faon thorique partir
des contraintes donnes par Ls-Dyna)

1 + 2 + 3
p=
3
Eq. 6-20
-Formule 3 (avec les contraintes principales calcules par Ls-Dyna)

1 + 2 + 3
p=
3
Eq. 6-21
La comparaison sur la pression isostatique nous permet aussi de mettre
en vidence la trace du tenseur de contrainte qui doit rester invariante. Ainsi cela
nous permettra de vrifier si les calculs effectus prcdemment par le logiciel
LS-Dyna font en sorte de conserver la trace du tenseur de contrainte.

Comparaison

56
 Modlisation de linteraction fluide/structure

Conclusion

Les trois formules nous donnent des valeurs identiques ce qui prouvent
que le calcul des contraintes principales effectu par LS-Dyna fait en sorte de
conserver la trace du tenseur de contrainte.
En revanche la pression releve par LS-Dyna est trs loign des
rsultats que lon t en droit dattendre. Encore une fois certaines donnes
fournies par LS-Dyna semblent suspectes.

57
 Modlisation de linteraction fluide/structure

7 Comportement dune structure sous leffet

dynamique du fluide

7.1 Dambreak exprience

7.1.1 Prsentation de lexprience

Lexprience Dambreak peut tre perue comme un modle simple

dcoulement deau sur le pont dun navire. Les dimensions du modle
consistent en un rservoir de 3,22m de longueur par 1m de large et par 1m de
haut. La partie droite du rservoir est dans un premier temps ferme par une
porte guillotine. Derrire cette porte, 0,55m de hauteur deau est prte
scouler ds que la porte sera leve. A lintrieur du rservoir est place une
structure reprsentant le modle lchelle dun container sur le pont dun
bateau.
Pendant lexprience, des mesures sur les forces, la pression et la
hauteur de leau sont effectues aux diffrents points suivant ce modle.

Points de mesure du
niveau deau

Points de mesure de la pression

numrots de 1 7.
Fig. 55 Points de mesures

Lexprience est ralise avec et sans le container. On pourra dans un

premier temps modliser un rservoir rigide puis lastique sans la prsence dun
container, et ainsi vrifier la hauteur deau au point 1, ou encore la pression sur
le mur du fond du rservoir

58
 Modlisation de linteraction fluide/structure

Fig. 56 modlisation de lexprience

Fig. 57 arrt sur image exprience et conflow [9]

7.1.2 Rsultats et comparaison avec lexprience

7.1.2.1 Ecoulement

Pour valider la modlisation du fluide lcoulement est primordiale. En

effet lcoulement relie les paramtres physiques comme la vitesse, la
viscosit
Comme nous pouvons le voir sur la courbe ci-dessous la hauteur deau
au centre du rservoir (1) lors de la simulation se rapproche des niveaux deau
obtenus lors des essais en laboratoire. On peut donc en conclure que
lcoulement de leau lors de la simulation se rapproche suffisamment dun
coulement deau rel dans une cuve pour tre considr comme raliste. Au-
del de la hauteur de leau, cest bien la vitesse du dplacement de leau qui est
vrifie [hauteur en fonction du temps].

59
 Modlisation de linteraction fluide/structure

niveau deau [1]

0,60
experience 1
experience 2
ls-dyna
0,50

0,40
niveau [m]

0,30

0,20

0,10

0,00
0,00 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00 4,50
temps [s]

Fig. 58 volution de la hauteur [1] deau

7.1.2.2 Pression

pression[1] sur mur du fond

14,00
ls-dyna rigide
experience 1
12,00 experience 2
ls-dyna lastique

10,00

8,00
pression [kPa]

6,00

4,00

2,00

0,00
0,00 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00 4,50

-2,00
temps [s]

Fig. 59 volution de la pression 3cm de hauteur sur le mur du fond de la cuve

60
 Modlisation de linteraction fluide/structure

Pression [3] sur mur du fond

5,00
experience 1
4,50 experience 2
ls-dyna lastique
4,00 ls-dyna rigide

3,50

3,00
pression [kPa]

2,50

2,00

1,50

1,00

0,50

0,00
0,00 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00 4,50
-0,50
temps [s]

Fig. 60 volution de la pression 18cm de hauteur sur le mur du fond de la cuve

Pression [5] sur mur du fond

4,50
ls-dyna rigide
experience 1
4,00
experience 2
ls-dyna lastique
3,50

3,00

2,50
pression [kPa]

2,00

1,50

1,00

0,50

0,00
0,00 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00 4,50

-0,50
temps [s]

Fig. 61 volution de la pression 33cm de hauteur sur le mur du fond de la cuve

Nous pouvons observer sur les diffrentes courbes ci-dessus, lvolution

de la pression en diffrents points du mur du fond de la cuve, dans le cas de
lexprience, de structure rigide et lastique.

61
 Modlisation de linteraction fluide/structure

La structure relle tait ralise en plexiglas dpaisseur 2cm, pouvant

donc tre considre comme rigide.
Les rsultats donns par LS-Dyna ont t filtr pour obtenir les courbes
dallures de la pression sur le mur du fond. Cela permet une meilleure lisibilit au
dtriment dune perte d informations sur linstantanit des pressions et sur
linstant.

Les rsultats obtenus par LS-Dyna ne collent pas parfaitement

lexprience mais ne sont pas non plus aux antipodes. Dans le cas dune
structure rigide (structure qui doit coller le mieux lexprience) les courbes de
pressions obtenues ont tendance reflter une moyenne en lissant les extrema
des pressions obtenues par lexprience.
Au contraire, dans le cas dune structure lastique les courbes de
pressions obtenues ont tendances englober les maxima, et malgr le
filtre/lissage elles refltent mieux linstantanit de leffort.
Dans le cas de simulation avec structure souple, les rsultats obtenus
semblent donc plus sensibles aux diffrents paramtres rentrant en jeu.
Cette mme tude sera poursuivie avec un maillage de la structure plus
affin. (non prsente dans ce rapport)

7.2 Autres simulations

Soit une plaque de 3mm en acier S235 (Fig. 37 courbe du matriau) place
un mtre dune cuve deau de hauteur deau 55cm. On libre leau aprs 3s le
temps pour celle-ci datteindre ses proprits physiques.

Fig. 62 modle effet dynamique de leau sur la structure

62
 Modlisation de linteraction fluide/structure

7.2.1 Force dimpact

7.2.1.1 Thorie

Lors de limpact de leau sur la plaque on peut considrer de faon plus ou

moins intuitive que lon est dans le cas dun jet deau venant vitesse constante
et section constante impacter la plaque. Cette considration nous permettra
ainsi destimer la valeur de leffort que nous devons obtenir lors de cet impact.
Cependant la difficult de cibler correctement le moment/temps de limpact ainsi
que la section deau impactant la plaque peut rendre ce calcul difficile et erron.
Daprs la thorie, la force horizontale si lon nglige les frottements est
donne par la relation suivante :

F = vQ
Eq. 7-1
Avec Q = vS
F = v S
Eq. 7-2

7.2.1.2 Vitesse et surface dimpact

Fig. 63 Vitesse dimpact de leau sur la structure

La vitesse du fluide est de 4m/s au moment de limpact sur la plaque.

63
 Modlisation de linteraction fluide/structure

Fig. 64 impact de leau sur la structure

Fig. 65 section deau impactant la structure

7.2.1.3 Force thorique

On estime quune hauteur deau de 5cm vient impacter la plaque, ce qui

en raison des dimensions de la plaque nous donne une surface de 137,5cm. On
a vu prcdemment que la vitesse dimpact est estime 4m/s. On obtient donc
une force dimpact de :

F = 998,21 * 4 * 137,5.10 4
Eq. 7-3

F = 219,6kN
Eq. 7-4

64
 Modlisation de linteraction fluide/structure

7.2.1.4 Force Ls dyna

Fig. 66 volution de leffort de leau sur la plaque

Leffort engendr par le premier impact de leau sur la plaque est

reprsent par la premire droite de la courbe dans lintervalle de temps [3,3s ;
3,55s]. En effet le lissage de la courbe fait que lon nobtient pas une lecture de
force instantane mme si leffort est lui impulsionnel.

Fig. 67 Force dimpact (LS-Dyna)

Nous obtenons une force dimpact selon LS-Dyna quivalente 220kN.

65
 Modlisation de linteraction fluide/structure

7.2.1.5 Conclusion

Nous avons obtenu un rsultat thorique et informatique (LS-Dyna)

rsolument proche lun de lautre, savoir 219,6kN pour une valeur avoisinant
les 225kN. Le transfert des efforts du fluide sur la plaque semble donc cohrent
et tmoigne du bon fonctionnement de lalgorithme de couplage.

7.2.2 Evolution de la force en fonction de la plaque

Nous allons dans cette partie tudier la force dimpact de leau en fonction
de diffrents critres (largeur de la plaque, ouverture dans la plaque).

7.2.2.1 Diffrents cas de plaques

-1ere plaque : voir tude ci-dessus

-2eme plaque : on prend la mme modlisation que ci-dessus, cependant on

ralise une fente de 0,10m*1,27m

Fig. 68 Plaque muni dune ouverture

66
 Modlisation de linteraction fluide/structure

-3eme plaque : on rduit la largeur de la plaque prcdente de moiti

Fig. 69 Plaque de largeur rduite

7.2.2.2 Rsultats

Nous considrerons que la vitesse dimpact et que la hauteur de la

section deau impactant restent inchanges.

-1ere plaque :

Nous utiliserons ce rsultat

comme rfrence de comparaison

Fig. 70 force dimpact sur la 1re structure

67
 Modlisation de linteraction fluide/structure

-2eme plaque :

Nous pouvons faire ici deux

remarques:
-La force dimpact reste inchange
ceux qui rvle un dfaut dans notre
modle (fente trop troite par rapport
aux lments de discrtisation du
fluide par exemple dou dcoulement
travers la plaque dans un premier
temps). Etant donn la taille de la fente
leffort devrait tre rduit au minimum
de 36%
-La force maximale est quant elle
rduite de plus de 50%, ce qui est le
rsultat de la diminution de la section
(36%) et de la forme de la courbe de
pression (voir Fig. 72 courbe de
pression sur une plaque), une fois Fig. 71 force dimpact sur la 2me structure
que lcoulement travers la plaque
est tabli.

Fig. 72 courbe de pression sur une plaque

-3eme plaque :

Nous retrouvons un rsultat

relativement prvisible. En effet la
section de la plaque tant rduite de
moiti cela implique que la force en est
rduite dautant. Si lon regarde les
points caractristiques, force dimpact
de 115kN ou la force maximale 250kN,
on retrouve bien des forces deux fois
moindre que dans notre premier cas
dtude.

Fig. 73 force dimpact sur la 3me structure

68
 Modlisation de linteraction fluide/structure

7.3 Observation du couplage par pnalits

Le champ des forces dinteraction est valu et appliqu en respectant le

principe de laction et de la raction. La force dinteraction est estime par une
mthode de pnalisation consistant permettre une lgre pntration entre les
deux maillages fluide (ALE) et structure (Lagrange). Une force de rappel
proportionnelle cette distance est alors value par lquation suivante et elle
est applique, par symtrie, sur le fluide et la structure.

F = -k * d
Eq. 7-5

o F, d et k sont respectivement la force, la pntration et une raideur. Cette

dernire est lobjet du problme.

Fig. 74 fuite a linterface fluide structure

Si la raideur est trop faible, les interpntrations deviendront

inacceptables et le fluide traversera la structure. En thorie, la raideur devrait
tre trs grande afin de limiter les pntrations. Cependant, si la raideur est trop
importante, la force de rappel trop grande devant les forces physiques
dstabilisera le calcul.

69
 Modlisation de linteraction fluide/structure

8 Conclusion

Ce projet de fin dtudes donne un premier aperu de la capacit du

logiciel LS-Dyna de rpondre de problme dinteraction fluide structure (ndlr : le
logiciel LS-Dyna est avant tout un logiciel de simulation de crash et donc pas
forcment en clin de solutionner ce type de problme). En effet je nai pas la
prtention de pouvoir affirmer que ce logiciel peut solutionner tous les types de
problmes faisant appel linteraction fluide structure.

Pour ce faire nous avons mis au point une stratgie qui sest droul en
deux parties. Aprs avoir dvelopp un modle de fluide convenable nous avons
tudi son impact dans le cas de cas simple de pression hydrostatique et dans
celui plus complexe dun fluide en mouvement.
Ltude dans le cas dune pression hydrostatique nous a rvl que les
rsultats obtenus par la simulation avec LS-Dyna collaient parfaitement aux
thories (thorie des plaques de Timoshenko, contraintes induites par la
pression hydrostatiques). Cependant nous avons pu relever que certaines
valeurs donnes par le logiciel, dans ce cas dtude, paraissaient errones. En
effet les valeurs des contraintes principales, de la pression isostatique ou encore
du critre de Von Mises divergeaient notoirement des valeurs que lon t en
droit dattendre.
Ltude dans le cas dun fluide en mouvement, en se rapportant une
exprience ralis prcdemment luniversit de Groningen, nous a permis de
comparer les rsultats obtenus dans le cas de structure souple ou rigide. Ainsi
nous avons pu constater que dans le cas de structure lastique les contraintes
sont amplifies, et cela pour un mme coulement, en comparaison avec une
structure rigide. Cela prouve bien que nous ne pouvons pas ngliger les
interactions fluide/structure lors des phases de dimensionnement.
Nous avons aussi vrifi les efforts obtenus lors de limpact de leau sur
diffrentes plaques souples. Ceux-ci nous ont conforts dans le comportement
du logiciel LS-Dyna car les rsultats obtenus semblaient cohrents.

Nous pouvons donc dire que ce projet met en vidence le bon

fonctionnement du couplage dans le cas simple de structures soumises une
pression hydrostatique. Dans le cas de fluide en mouvement les rsultats sont
plus difficiles analyser. On a pu cependant mettre en vidence les diffrences
entre structures rigides ou souples soumises au mme coulement. Aussi nous
avons pu montrer la cohrence des forces gnres par limpact dun fluide sur
une structure souple.

70
 Modlisation de linteraction fluide/structure

Plus personnellement ce projet ma permis dacqurir des comptences

nouvelles, que ce soit dans le domaine des lments finis, de la mcanique des
fluides mais surtout dans le domaine de linteraction fluide/structure.
En effet ce projet ma permis de dcouvrir les bases de la thorie de
linteraction fluide/structure, un domaine en pleine expansion preuve en est le
nombre de thses effectues ces dernires annes sur le sujet. Aussi il serait
illusoire de croire quen 20 semaines jai pu assimiler toute cette thorie.
Ce projet a galement t pour moi loccasion de mettre en uvre de
nombreuses connaissances acquises tout au long de ma scolarit lINSA
Strasbourg. En effet jai t amen utiliser la mthode des lments finis pour
les phases de comprhension du fonctionnement du logiciel, lanalyse des
structures et la mcanique des fluides lorsquil a fallu analyser les diffrents
rsultats.

71
 Modlisation de linteraction fluide/structure

Tableau des figures

Fig. 1: Pont de Tacoma [3] ............................................................................................11
Fig. 2: Varit de plate-formes ptrolires [3] ................................................................12
Fig. 3: organigramme LS-Dyna [7].................................................................................13
Fig. 4 Organisation des commandes de programmation [7]...........................................14
Fig. 5 cycle Lagrangien .................................................................................................16
Fig. 6 cycle Eulerien ou ALE .........................................................................................16
Fig. 7: volution du maillage..........................................................................................19
Fig. 8: notation en fonction de la formulation .................................................................20
Fig. 9 fluide dans les trois configurations.......................................................................21
Fig. 10 structure dans la configuration initiale et courante .............................................21
Fig. 11 transfert des donnes entre deux codes............................................................23
Fig. 12 algorithme du coupleur ......................................................................................24
Fig. 13 schma du contact par pnalisation ..................................................................25
Fig. 14 lment de discrtisation de la structure............................................................25
Fig. 15 maillage de lune des structures ........................................................................26
Fig. 16 discrtisation du fluide .......................................................................................27
Fig. 17 lment sous intgr .........................................................................................27
Fig. 18: programmation de lquation dtat de leau ....................................................30
Fig. 19 maillage du systme..........................................................................................31
Fig. 20 Contour de la densit du fluide ..........................................................................31
Fig. 21: iso surfaces de la densit du fluide...................................................................31
Fig. 22 modle pression hydrostatique..........................................................................32
Fig. 23 diffrents types de chargements........................................................................32
Fig. 24 points de mesure de la pression hydrostatique..................................................33
Fig. 25 test gravit instantane .....................................................................................33
Fig. 26: test temps de chargement 0,5s.........................................................................34
Fig. 27: test temps de chargement 1s ...........................................................................34
Fig. 28: test temps de chargement 2s ...........................................................................35
Fig. 29: test advection mise jour tous les 3 cycles de calcul .......................................36
Fig. 30: test advection du premier ordre ........................................................................36
Fig. 31 modle effort hydrostatique ...............................................................................37
Fig. 32 pression hydrostatique ......................................................................................38
Fig. 33 Somme des forces sur le contour de la plaque ..................................................38
Fig. 34 contrainte du leau suivant x...........................................................................39
Fig. 35 Somme des contraintes de la plaque suivant x..................................................39
Fig. 36 programmation du matriau de la structure .......................................................41
Fig. 37 courbe du matriau............................................................................................42
Fig. 38 dmarche de la thorie des plaques..................................................................42
Fig. 39 flexion de plaque ...............................................................................................43
Fig. 40 points de mesure 1ere vrification .....................................................................45
Fig. 41 dplacements des diffrentes plaques dans le 1er cas.......................................45
Fig. 42 graphique proportionnalit flche-paisseur ......................................................46
Fig. 43 modlisation plaque encastre sur le contour de la pression.............................47
Fig. 44 dplacements des plaques dans le 2eme cas ...................................................49
Fig. 45 comparaison thorie/ls-dyna pour t= 2mm.........................................................49
Fig. 46 comparaison thorie/ls-dyna pour t=5mm..........................................................49

72
 Modlisation de linteraction fluide/structure

Fig. 47 Correspondance entre les deux thories ...........................................................50

Fig. 48 rapport w/h ........................................................................................................50
Fig. 49 comparaison thorie/ls-dyna pour t=1mm..........................................................51
Fig. 50 points de mesure des contraintes ......................................................................52
Fig. 51 cercle de mohr...................................................................................................53
Fig. 52 : 1ere contrainte principale ................................................................................54
Fig. 53 : 2eme contrainte principale...............................................................................54
Fig. 54 :3eme contrainte principale................................................................................55
Fig. 55 Points de mesures.............................................................................................58
Fig. 56 modlisation de lexprience .............................................................................59
Fig. 57 arrt sur image exprience et conflow [9] ..........................................................59
Fig. 58 volution de la hauteur [1] deau .......................................................................60
Fig. 59 volution de la pression 3cm de hauteur sur le mur du fond de la cuve ..........60
Fig. 60 volution de la pression 18cm de hauteur sur le mur du fond de la cuve ........61
Fig. 61 volution de la pression 33cm de hauteur sur le mur du fond de la cuve ........61
Fig. 62 modle effet dynamique de leau sur la structure ..............................................62
Fig. 63 Vitesse dimpact de leau sur la structure..........................................................63
Fig. 64 impact de leau sur la structure..........................................................................64
Fig. 65 section deau impactant la structure ..................................................................64
Fig. 66 volution de leffort de leau sur la plaque .........................................................65
Fig. 67 Force dimpact (LS-Dyna) .................................................................................65
Fig. 68 Plaque muni dune ouverture ............................................................................66
Fig. 69 Plaque de largeur rduite ..................................................................................67
Fig. 70 force dimpact sur la 1re structure ...................................................................67
Fig. 71 force dimpact sur la 2me structure .................................................................68
Fig. 72 courbe de pression sur une plaque....................................................................68
Fig. 73 force dimpact sur la 3me structure .................................................................68
Fig. 74 fuite a linterface fluide structure........................................................................69

73
 Modlisation de linteraction fluide/structure

Bibliographie

[1] AIT MOUDID.L, Couplage Fluide-Structure pour la simulation numrique des

coulements fluides dans une conduite `a parois rigides ou lastiques, en
prsence dobstacles ou non

[2] AQUELLET.N, Modlisation de limpact hydrodynamique par un couplage

fluide-structure.

[3] PIPERNO.S, Interaction Fluide structure

[4] RIBET.H, Modlisation de limpact sur leau dune structure souple par une
mthode lments finis explicite.

[5] TIMOSHENKO et WOINOWSKY-KRIEGER, Theory of plates and shells.

[6] WEMMENHOVE.R, Numerical simulation of two phase flow in offshore

environments.

[7] Numerical simulation of two phase flow in offshore environments, LS-Dyna

Keyword Users Manual

[8] Encyclopdie en ligne Wikipida.org, juin 2009 : http://fr.wikipedia.org/ et

http://en.wikipedia.org/

[9] University of groningen : http://www.math.rug.nl/~veldman/comflo/comflo.html

[10] Site internet de Force Technology, juin 2009 : http://www.force.no/

74
 Modlisation de linteraction fluide/structure

Annexes

75
 Modlisation de linteraction fluide/structure

A. Equations de Navier Stokes

En mcanique des fluides, les quations de Navier-Stokes sont des

quations aux drives partielles non-linaires qui dcrivent le mouvement des
fluides dans l'approximation des milieux continus. Elles gouvernent par exemple
les mouvements de l'air de l'atmosphre, les courants ocaniques, l'coulement
de l'eau dans un tuyau, et de nombreux autres phnomnes d'coulement de
fluides. Elles sont nommes d'aprs deux physiciens du XIXe sicle, Claude
Navier et George Stokes. Il est possible de dmontrer ces quations partir de
l'quation de Boltzmann.

Formule gnrale

La formulation diffrentielle de ces quations est la suivante:

quation de continuit (ou quation de bilan de la masse)

Eq: A-1

quation de bilan de la quantit de mouvement

Eq: A-2

quation de bilan de l'nergie

Eq: A-3

Avec:
t: le temps (unit SI : s) ;
: la masse volumique du fluide (unit SI : kg.m 3) ;
: la vitesse Eulerienne d'une particule fluide (unit SI : m.s 1) ;
p : la pression (unit SI : Pa) ;

: le tenseur des contraintes visqueuses (unit SI : Pa) ;

: la rsultante des forces massiques s'exerant dans le fluide (unit SI : N.kg
1) ;
e : l'nergie totale par unit de masse (unit SI : J.kg 1) ;
: le flux de chaleur perdu par conduction thermique (unit SI : J.m 2.s 1) ;
r : la perte de chaleur volumique due au rayonnement (unit SI : J.m 3.s 1).

76
 Modlisation de linteraction fluide/structure

B. Algorithme dadvection de Van Leer

Lalgorithme dadvection prsent ici est connu aussi sous le nom
dalgorithme MUSCL (Monotone Upwind Schemes for Conservations Laws). La
mthode de la cellule donneuse suppose que g (reprsentant les grandeurs
centrs dans la maille) est constante sur la maille. Van Leer remplace cette
distribution constante par morceaux par une distribution dun ordre plus leve.
Une distribution linaire par morceaux peut tre utilise:

Eq: B-1
Le transport de g est purement gomtrique et la grandeur t nest pas le
temps physique. Ainsi, quand g passe de ltat n n+1, cette grandeur nest pas
avance dans le temps: elle est projete du maillage ltat n sur le nouveau
maillage ltat n+1. Pour assurer la stabilit des schmas dadvection, le
dplacement du maillage dans un pas de temps Dt doit tre infrieur une
dimension caractristique du maillage L.
x est la coordonne dans le volume de llment, xnj+1/2 la coordonne du
centre de la maille et Snj+1/2 la pente de la droite. Pour assurer la monotonie,
cette pente est limite par:

Eq: B-2
avec :

Eq: B-3
Eq: B-4

o snj+1/2 est une approximation de la pente du second ordre obtenue en

interpolant par une parabole les valeurs de g aux centres des 3 mailles
adjacentes:

Eq: B-5

77
 Modlisation de linteraction fluide/structure

Avec

Eq: B-6
Un second minimum peut tre envisag en remplaant Eq.B-3 et Eq.B-4 par ce
qui suit:

Eq: B-7
Eq: B-8
avec

Eq: B-9
o

Eq: B-20
Eq: B-31
et

Eq: B-42
On obtient la nouvelle valeur de g :

Eq: B-53
Cette mthode dadvection est adopte dans les applications numriques.
Les autres variables values au centre de la maille lors du cycle Lagrangien
sont aussi projetes sur le nouveau maillage par lalgorithme dadvection
prsent. Les vitesses ntant pas centres dans la maille mais dtermines au
nuds, il faut modifier cet algorithme pour ladapter ladvection de la quantit
de mouvement.

78
 Modlisation de linteraction fluide/structure

C. Viscosit de Hourglass
Pour les lments brique comme pour les lments quadrangles, le
dsavantage dune intgration avec un seul point dintgration est la ncessit
de contrler les modes nergie nulle appels aussi modes de Hourglass. Ces
modes impliquent des dformations nulles aux points de quadrature:

Eq: C-1

Fig.C1: element brique

or

Eq: C-2

79
 Modlisation de linteraction fluide/structure

et lorsque les vitesses des nuds opposs diagonalement sont gales

pour i= 1, 2,3

Eq: C-3
Alors le taux de dformation dij (=0, =0, =0) sannule au point de
quadrature. Ce mode entrane alors un travail nul do le nom de mode nergie
nulle. Linterprtation gomtrique de ce phnomne est que le plan parallle au
gradient de vitesses passant par le point de quadrature ne se dforme pas dans
un mode de Hourglass. La raideur de ce mode est, en gnral, plus petite que la
plus basse des raideurs physiques et ainsi la frquence associe est aussi plus
petite que la plus petite frquence du problme physique. On observe alors des
oscillations perturbant les grandeurs physiques. De plus, le Hourglass est une
instabilit spatiale visible dans la dformation des mailles lors dun calcul
Lagrangien. Cette instabilit peut se communiquer et se dvelopper sur
lensemble du maillage comme une instabilit dadvection-diffusion Pour pallier
cette instabilit, 12 forces damortissement fi pour =1,4) et i=1, 2, 3 sont
introduites. Pour cela, on dfinit des vecteurs de forme du Hourglass i comme
tant covecteurs des vecteurs de base du Hourglass hi. Les vecteurs fi sont
orthogonaux un mouvement rigide et un champ de vitesses linaire en xiI
pour que les forces damortissement ne perturbent pas la physique de
lcoulement :

Eq: C-4

Eq: C-5
Les taux de dformations du Hourglass gnralises sont donnes par :

Eq: C-6

80
 Modlisation de linteraction fluide/structure

Les 12 forces nodales rsistant aux 12 modes de Hourglass sont alors crites
sous la forme suivante

Eq: C-7

avec ah, une constante dfinie par:

Eq: C-8
Ch est une constante dfinie par lutilisateur entre 0,05 et 0,15.

Fig.C2: mode de Hourglass

81
 Modlisation de linteraction fluide/structure

D. Ecoulement du fluide

t=0s
t=0,5s

t=1,3s t=2,2s

t=3,6s t=5s

82
 Modlisation de linteraction fluide/structure

E. Programme
bdam-master.k

Paramtres Hourglass

Paramtres de contrle pour les

calculs ALE :( eulrien, nombre
de cycle entre deux advections,
mthode d advection)

83
Modlisation de linteraction fluide/structure

Dfinition du chargement de la
gravit sur le fluide

84
Modlisation de linteraction fluide/structure

Inclusion de programmes
secondaires

85
 Modlisation de linteraction fluide/structure

bdam-parts.k

Dfinition du vide

Dfinition de leau

Dfinition de la
plaque retenant
leau

plaque

eau
vide

86
 F. Contraintes de Von Mises

Von Mises Contraintes

4,50E+08

4,00E+08

3,50E+08

3,00E+08 V-M formule1

V-M formule2
Contraintes [Pa]

V-M formule 3
2,50E+08 Ls-Dyna

2,00E+08

1,50E+08

1,00E+08

5,00E+07

0,00E+00
0 1 2 3 4 5
temps [s]

Formule 1 : calcul avec les contraintes principales thoriques dans le cas de

contraintes planes.
1
( 1 2 )
2

Formule 2 : calcul avec les contraintes principales thoriques

1
( 1 2 ) 2 + ( 1 3 ) 2 + ( 2 3 ) 2
2

Formule 3 : calcul avec les contraintes principales LS-Dyna

1
( 1 2 ) 2 + ( 1 3 ) 2 + ( 2 3 ) 2
2