Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

COURS Hydraulique Generale MEPA 2010

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 85

ECOLE NATIONALE DU GENIE DE LEAU ET

DE LENVIRONNEMENT DE STRASBOURG


HYDRAULIQUE
GENERALE




Rservoir deau potable de Strasbourg

FORMATION CES/MASTERE
Eau potable et assainissement

Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains)
Avant propos

Lhydraulique est incontournable dans le domaine de lenvironnement. En effet, elle a une
place dterminante dans la comprhension, lanalyse et le diagnostic des rseaux
dadduction deau potable, des stations de traitement, des rseaux dassainissement et des
rivires. De plus, le contrle de ces systmes ncessite une instrumentation qui oblige le
concepteur et lexploitant une connaissance pousse du fonctionnement hydraulique de
ces ouvrages.
Dun point de vue rglementaire, la directive 2000/60/CE du Parlement europen tablit
un cadre pour une politique communautaire dans le domaine de leau. Elle incite les Etats
membres (dont videmment la France) protger et restaurer la qualit de leurs ressources
en eau afin de parvenir un bon tat chimique et cologique. Leau est donc une
proccupation majeure dans notre civilisation.

Lobjectif de cet ouvrage destin aux techniciens et ingnieurs est de fournir les bases
ncessaires la comprhension et au calcul des phnomnes prsents en hydraulique
applique au gnie de leau et de lenvironnement. Chaque notion dhydraulique est
ponctue par une srie dexercices permettant dillustrer les concepts prsents. Les
exemples sont issus douvrages hydrauliques existant en rseau. Les techniques de calcul
qui sont associes la rsolution des quations mises en uvre sont labores dans un
souci defficacit.

Cet ouvrage est compos de plusieurs chapitres qui sont dcrits dans les paragraphes
suivants.

Avant dentrer directement dans la description mathmatique et la quantification des
phnomnes hydrauliques rencontrs en eau potable, en assainissement et en rivire, le
premier chapitre est consacr la description physique et phnomnologique des
diffrents types dcoulements ainsi que des forces en prsence. Elle permet aux non-
initis de dcouvrir lhydraulique de faon qualitative par le biais de lexprimentation.

Le deuxime chapitre sintresse aux fluides au repos. Par exemple, la plupart des
dispositifs de contrle de dbit tels que les vannes autorgules ncessitent un mcanisme
pouvant tre approch par le biais de lhydrostatique. Cette partie est donc consacre
laction de leau sur les parois et sur les corps immergs.

Le chapitre suivant est ddi lhydraulique en charge. Ce type dcoulement est le plus
souvent rencontr dans les rseaux dadduction deau potable et parfois en assainissement.
Aprs un rappel des quations de Bernoulli, le paragraphe suivant sintresse aux pertes de
charge linaires et singulires. On aborde ensuite le trac de la ligne pizomtrique et de la
ligne de charge qui permettent de caractriser le fonctionnement dun rseau en terme de
sur ou sous pression. Un dernier sous-chapitre traite du coup de blier avec les techniques
de calcul associes.

Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Bibliographie



AGHTM : Les stations de pompage deau, Editions Tec et Doc (2000).

BERTRAND-KRAJEWSKI J.L., Mesures en hydrologie urbaine et assainissement, d.
Tec et doc, ed. 2000.

CARLIER M. : Hydraulique gnrale et applique, Editions Eyrolles (1972).

COMOLET R., Mcanique exprimentale des fluides, Masson, ed.1982.

GRAF W. H., ALTINAKAR M. S. : Hydrodynamique : Une introduction, Trait de
Gnie Civil, Ecole polytechnique fdrale de Lausanne, Presse
polytechnique et universitaire romanes (1995).

GRAF W. H., ALTINAKAR M. S. : Hydraulique fluviale : coulement permanent
uniforme et non uniforme, Tome 1, Trait de Gnie Civil, Ecole
polytechnique fdrale de Lausanne, Presse polytechnique et
universitaire romanes (1993).

GRAF W. H., ALTINAKAR M. S. : Hydraulique fluviale : coulement non permanent
et phnomnes de transport, Tome 2, Trait de Gnie Civil, Ecole
polytechnique fdrale de Lausanne, Presse polytechnique et
universitaire romanes (1996).

HAGER W. H. : Wastewater hydraulics theory and practice, Springer, ed. 1999.

LENCASTRE A. : Hydraulique gnrale, Editions Eyrolles (1996).

LESIEUR M. : La turbulence, Presses Universitaires de Grenoble, Ed. 1994.

SCHIESTEL R. : Modlisation et simulation des coulements turbulents, Editions
Herms (1993).

SINNIGER R.O., HAGER W. H. : Constructions hydrauliques : Ecoulements
stationnaires, Trait de Gnie Civil, Ecole polytechnique fdrale de
Lausanne, Presse polytechnique et universitaire romanes (1989).

VIOLET P.L., CHABARD J.P., Mcanique des fluides applique, Presses des ponts et
chausses, ed. 1998.




Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Sommaire

CHAPITRE I : CARACTERISTIQUES DES ECOULEMENTS ..................................................... 6
1. - ECOULEMENTS EN CHARGE ..................................................................................................................... 6
1.1. - Rgimes dcoulements .................................................................................................................. 6
1.2. - Distribution des vitesses dans un tuyau rectiligne ......................................................................... 7
1.3. - Elargissement et rtrcissement dans une conduite ....................................................................... 8
1.4. - Sortie dun rservoir ...................................................................................................................... 9
1.5. - Coude ............................................................................................................................................. 9
1.6. - Jonction et bifurcation ................................................................................................................. 10
1.7. - Obstacles dans un coulement en charge .................................................................................... 10
1.8. - Le coup de blier .......................................................................................................................... 11
2. - ECOULEMENTS A SURFACE LIBRE .......................................................................................................... 13
2.1. - La surface libre ............................................................................................................................ 13
2.2. - Distribution des vitesses dans un canal ....................................................................................... 13
2.3. - Notion dinfluence aval ................................................................................................................ 14
2.4. - Courbes de remous ...................................................................................................................... 15
2.5. - Le ressaut ..................................................................................................................................... 15
2.6. - Les seuils et dversoirs ................................................................................................................ 15
3. - PROPRIETES DES LIQUIDES .................................................................................................................... 17
3.1. - Masse volumique .......................................................................................................................... 17
3.2. - Poids spcifique ........................................................................................................................... 17
3.3. - Viscosit ....................................................................................................................................... 17
3.4. - Pression de vapeur saturante ....................................................................................................... 19
3.5. - Tension superficielle .................................................................................................................... 19
4. - CARACTERISATION DES FORCES DANS UN ECOULEMENT ...................................................................... 20
4.1. - Les forces ..................................................................................................................................... 20
4.1.1. - Les forces de volumes ........................................................................................................................... 20
4.1.2. - Les forces de surfaces ............................................................................................................................ 21
4.2. - Importance des diffrentes forces ................................................................................................ 21
4.2.1. - Le Nombre dEuler ................................................................................................................................ 21
4.2.2. - Le Nombre de Reynolds ........................................................................................................................ 21
4.2.3. - Le Nombre de Froude ............................................................................................................................ 22
4.2.4. - Le Nombre de Cauchy ........................................................................................................................... 22
4.2.5. - Le Nombre de Weber ............................................................................................................................ 22
5. - LES EQUATIONS DE BASE ...................................................................................................................... 22
6. - OUTILS MATHEMATIQUES ..................................................................................................................... 23

CHAPITRE II : HYDROSTATIQUE ....................................................................................... 24
1. - EQUATIONS DE LHYDROSTATIQUE ....................................................................................................... 24
2. - VARIATION DE LA PRESSION DANS UN FLUIDE INCOMPRESSIBLE ........................................................... 25
3. - VARIATION DE LA PRESSION DANS UN FLUIDE COMPRESSIBLE ............................................................. 27
4. - FORCES HYDROSTATIQUES SUR LES PAROIS .......................................................................................... 27
4.1. - Paroi plane en position incline .................................................................................................. 28
4.2. - Paroi surface gauche ................................................................................................................ 30
5. - FORCES HYDROSTATIQUES SUR DES CORPS IMMERGES ......................................................................... 31
5.1. - Forces dArchimde ..................................................................................................................... 31
5.2. - Equilibre des corps immergs ...................................................................................................... 31

Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
CHAPITRE III : LHYDRAULIQUE EN CHARGE ............................................................. 33
1. - LEQUATION DE CONTINUITE ................................................................................................................ 33
2. - EQUATIONS DE BERNOULLI .................................................................................................................. 35
2.1. - Cas des fluides parfaits ................................................................................................................ 35
2.1.1. - Equations dEuler .................................................................................................................................. 35
2.1.2. - Equations de Bernoulli .......................................................................................................................... 38
2.1.3. - Reprsentation graphique ...................................................................................................................... 39
2.1.4. - Interprtation nergtique ...................................................................................................................... 39
2.2. - Ecriture de Bernoulli pour les fluides rels ................................................................................. 47
3. - EVALUATION DES PERTES DE CHARGE .................................................................................................. 48
3.1. - Les pertes de charge linaire ....................................................................................................... 48
3.1.1. - Profil des vitesses dans une section circulaire en charge ....................................................................... 48
3.1.2. - Expression de la perte de charge linaire ............................................................................................... 51
3.2. - Les pertes de charge singulire ................................................................................................... 57
3.3. - Exemple de trac de perte de charge ........................................................................................... 63
4. - POSITION DE LA LIGNE PIEZOMETRIQUE ................................................................................................ 64
5. - LES CONSTRUCTIONS GRAPHIQUES ....................................................................................................... 66
6. - LE COUP DE BELIER ............................................................................................................................... 68
6.1. - Introduction ................................................................................................................................. 68
6.2. - Formulation mathmatique .......................................................................................................... 71
6.3. - Coup de blier de masse .............................................................................................................. 71
6.4. - Coup de blier dondes ................................................................................................................ 73
6.4.1. - Conservation de la masse pour un coulement dans lequel se propage une onde .................................. 73
6.4.2. - Clrit des ondes .................................................................................................................................. 74
6.4.3. - Equations aux caractristiques ............................................................................................................... 75
6.4.4. - Application la fermeture instantane dune vanne (sans pertes de charge dans la conduite) .............. 77
6.4.5. - Application la fermeture progressive dune vanne (sans pertes de charge dans la conduite) .............. 81
6.4.6. - Application la fermeture instantane dune vanne (prise en compte simplifie des pertes de charge) 82
6.5. - Protections anti-blier ................................................................................................................. 83
6.6. - Rfrences bibliographiques sur les coups de blier ................................................................... 85


Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Chapitre I :
CARACTERISTIQUES DES ECOULEMENTS

Lhydraulique est ltude des coulements.
On distingue deux types dcoulements :
les coulements en charge, dans lesquels leau remplit compltement la
canalisation, cest le cas notamment des rseaux deau potable,
les coulements surface libre (interface entre leau et lair), cest le cas des
rivires et des rseaux dassainissement.

1. - ECOULEMENTS EN CHARGE
1.1. - Rgimes dcoulements
Le rgime dun coulement se caractrise par la fluctuation temporelle des vitesses
et des pressions au sein de la veine liquide.

Exprience :
Un premier rservoir deau de niveau constant est vidang par un tuyau. Une vanne
place lextrmit du tuyau permet de faire varier le dbit Q (m
3
/s). Un deuxime tuyau
est plac lintrieur du rservoir. Il contient un colorant et permet dobtenir un mince
filet fluide color au centre du tuyau.

Injection dun colorant
Vanne de
rgulation du dbit
Sens de
lcoulement
Zone dobservation
Rservoir
Mesure de la
vitesse en un point

Quand la vitesse est trs faible (quelques millimtres par seconde) le filet color reste bien
dfini, rectiligne et parallle laxe du tuyau. Le rgime est dit laminaire. Lcoulement
laminaire est rare dans le domaine de lhydraulique de leau potable et de lassainissement,
toutefois il nest pas inexistant.
1. Ecoulements en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
7
Les figurent suivantes montrent un coulement laminaire dans une veine liquide et
lvolution de la vitesse en un point (vitesse ponctuelle) en fonction du temps.


Sens de lcoulement
Filet de colorant
Injection
de colorant

Evolution de la vitesse ponctuelle en fonction du
temps en rgime laminaire
0
0.0002
0.0004
0.0006
0.0008
0.001
0.0012
0.0014
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Temps (s)
V
i
t
e
s
s
e

(
m
/
s
)


Quand la vitesse est plus leve, le filet devient ondul et trs instable. Il se mlange
rapidement au fluide ambiant. Des tourbillons de diffrentes tailles apparaissent. Le
rgime est dit turbulent.



Sens de lcoulement
Filet de colorant
Injection
de colorant

Evolution de la vitesse ponctuelle en fonction du
temps en rgime turbulent
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Temps (s)
V
i
t
e
s
s
e

(
m
/
s
)


La turbulence se caractrise donc par la cration de tourbillons. ils mlangent les matires
dissoutes dans leau, comme par exemple le chlore dans un rseau deau potable ou le rejet
dune station de traitement des eaux uses dans une rivire. La mise en place dun
agitateur dans un bassin cre de la turbulence et ainsi il tend homogniser les matires
dissoutes.
1.2. - Distribution des vitesses dans un tuyau rectiligne
En fonction du rgime dcoulement (laminaire ou turbulent) les forces dans
lcoulement (elles seront dtailles dans un prochain chapitre) sont diffrentes. Cela a
pour effet une distribution de la vitesse ponctuelle moyenne dans le temps (que lon
appelle vitesse moyenne temporelle) lintrieur du tuyau qui est diffrente suivant le
rgime dcoulement.
La connaissance du profil de vitesse dans une conduite permet de calculer le dbit. La
plupart des capteurs en rseau mesurent la vitesse moyenne temporelle en quelques points
et reconstituent le profil complet de la vitesse pour en dduire le dbit.
Le graphique suivant compare la forme du profil de la vitesse en rgimes laminaire et
turbulent dans une conduite circulaire. En abscisse est reprsent le rapport entre la
distance la paroi (r) et le rayon de la conduite( ) R . Lordonne reprsente le
I- Caractristiques des coulements
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
8
pourcentage du rapport entre la vitesse maximale du profil et la vitesse ponctuelle
r
V
R
| |
|
\
.
On constate que la vitesse maximale est au centre de la canalisation. Par contre, le profil
de vitesse turbulent varie beaucoup plus quen laminaire au voisinage de la paroi. Cette
zone de fort gradient de vitesse est appele couche limite.


Evolution de la vitesse dans une conduite circulaire
0
10
20
30
40
50
60
70
80
90
100
0 0,2 0,4 0,6 0,8 1 0,8 0,6 0,4 0,2 0
Distance la paroi r/R
V
i
t
e
s
s
e

a
d
i
m
e
n
s
i
o
n
n
e
l
l
e


V
/
V
m
a
x

e
n

%
Centre de la conduite
Profil de vitesse
Laminaire
Profil de vitesse
Turbulent
Zone de fort
gradient de vitesse


1.3. - Elargissement et rtrcissement dans une conduite
La transition entre deux conduites de diamtre diffrent pour un coulement en
charge provoque une rpartition transversale des vitesses longitudinales totalement
diffrentes de celles vues prcdemment.
On constate une zone ou les veines liquides proches de la paroi se dcollent sur une
longueur L. Dans cette zone, on observe des recirculations lorigine de perturbations
importantes dans lcoulement.



L

Elargissement

Zone de dcollement
de veine liquide
L

Rtrcissement

1. Ecoulements en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
9
1.4. - Sortie dun rservoir
La sortie dun rservoir provoque galement une modification de la rpartition de
la vitesse et une zone de dcollement de la veine liquide. On remarque que les lignes de
courant (ligne enveloppe du champ de vitesse) se resserrent au passage de la zone de
recirculation. La prsence de cette zone a pour effet de diminuer lgrement la section de
passage de lcoulement.


Sens de
lcoulement
Zone de
recirculation
Resserrement des
lignes de courants

1.5. - Coude
En raison de la courbure de la conduite, un mouvement hlicodal des lignes de
courant peut stablir. La prsence de ce mouvement de rotation persiste sur une longueur
en aval du coude ( peu prs 50 fois le diamtre de la conduite). Lorsque la courbure est
importante, des zones de recirculation peuvent apparatre.
Le mcanisme interne de lcoulement au travers dun coude nest pas encore bien connu,
malgr le grand nombre de chercheurs qui ont analys cet lment important.


Zone de
recirculation
A
A
A-A
Mouvement
hlicodal

Evolution des lignes de courant dans un coude circulaire

I- Caractristiques des coulements
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
10
1.6. - Jonction et bifurcation
Lcoulement dans une jonction provoque une zone de sparation et une zone de
mlange. Lcoulement de la branche dans laquelle la vitesse est la plus faible est entran
par leau provenant de la branche ayant une vitesse plus leve. Ce phnomne est
lorigine de dclrations et dacclrations de lcoulement dans les deux branches.

Jonction
Zone de
sparation
Zone de mlange


Lcoulement dans une bifurcation se comporte de manire lgrement diffrente du fait
de la drivation latrale. Une zone de sparation apparat galement dans la branche
latrale.
1.7. - Obstacles dans un coulement en charge
Les figures suivantes montrent lvolution des tourbillons laval dun obstacle
dans un coulement en fonction de lintensit de la vitesse.

Obstacle dans un coulement.

Sens de
lcoulement

Lcoulement est faible
vitesse. On observe lapparition
de deux tourbillons symtriques
attachs lobstacle. Cest le
dbut de la turbulence.

Laugmentation lgre de la
vitesse entrane un allongement
des deux tourbillons.

A vitesse plus importante les
tourbillon commencent
osciller, puis se dtachent.

Laugmentation de la vitesse a
pour effet daccrotre la
frquence des oscillations.


Les vannes opercule ou les vannes papillons crent un obstacle dans les conduites. Ce
sont des organes de contrle du dbit.
1. Ecoulements en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
11

Vanne opercule

Vanne papillon

1.8. - Le coup de blier
Le coup de blier est un ensemble de phnomnes hydrauliques complexes
provenant de la modification brutale du champ de vitesse et de pression dans une conduite.
Limportance technique des coups de blier est extrmement grande. En effet, lamplitude
et la soudainet des surpressions et dpressions qui sont mises en jeu peuvent dtriorer la
conduite ou les organes qui y sont branches. Il est ncessaire par consquent dempcher
la production de telles perturbations ou du moins de les attnuer.
Le coup de blier peut se produire, par exemple, dans le cas de la fermeture brutale dune
vanne.

Rservoir
Vanne
ouverte
Sens de
lcoulement
Entre Sortie
Vitesse Uo

Considrons, dans un premier temps, un
rservoir se vidangeant par une vanne en
position ouverte. La vitesse dans la
canalisation est Uo.


Vanne
ferme
Rservoir
U=0
Vitesse Uo
Zone de
surpression

A un instant donn, la fermeture brutale de
la vanne stoppe instantanment le fluide au
voisinage de celle-ci. Compte tenu de son
inertie, leau pousse fortement sur la vanne.
Une surpression importante se cre dans la
zone o la vitesse sannule.


Rservoir
Vanne
ferme
Uo
Gonflement
de la conduite
U=0
Dplacement
de londe
Zone de
surpression
Onde de
choc

A lamont de la vanne, leau se comprime
faiblement et la canalisation gonfle. Tranche
par tranche les particules liquides
transforment la vitesse en surpression. On
observe ainsi la formation et la progression
dune onde de choc. Celle-ci remonte de la
vanne o elle a pris naissance vers le
rservoir.
La vitesse de londe peut avoisiner les 1000
m/s en fonction du type de matriau de la
conduite.


I- Caractristiques des coulements
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
12

Rservoir
Vanne
ferme
Gonflement
de la conduite
U=0
Zone de
surpression
Pression impose
par le rservoir

Londe de choc atteint le rservoir. La
conduite entire est en surpression et la
vitesse de leau est nulle. Le rservoir est
niveau deau constant. Sa pression ne
pouvant quasiment pas varier, le rservoir
impose donc la pression lentre de la
conduite.


Rservoir
Vanne
ferme
Uo
U=0
Dplacement
de londe
Onde de
choc

La pression dans le rservoir tant beaucoup
plus faible que la surpression dans la
conduite, la canalisation, au voisinage du
rservoir, va retrouver sa forme initiale. La
diminution du diamtre de la conduite va
chasser leau vers le rservoir la vitesse
Uo. De proche en proche se phnomne se
produit et cre ainsi une onde de choc se
dplaant cette fois-ci vers la vanne.


Rservoir
U=0
Vitesse Uo
Zone de
dpression

Ds que londe atteint la vanne, la conduite
entire retrouve quasiment sa pression
dorigine et sa vitesse initiale inverse. La
vanne tant toujours ferme, la vitesse de
leau son voisinage est nulle.


Rservoir
Vanne
ferme
Uo
U=0
Dplacement
de londe Zone de
dpression
Onde de
choc

Compte tenu de son inertie, leau au
voisinage de la vanne entre en dpression et
la conduite diminue son diamtre. Une onde
de choc en dpression se cre et se dplace
vers le rservoir.


Rservoir
Vanne
ferme
Gonflement
de la conduite
U=0
Pression impose
par le rservoir

Comme prcdemment, le rservoir impose
la pression au niveau de lentre de la
conduite.



Rservoir
Vanne
ferme
Uo
U=0
Dplacement
de londe
Onde de
choc

La conduite retrouve alors sa forme initiale.
Le volume deau manquant est pris dans le
rservoir et la conduite retrouve sa vitesse
initiale Uo.



Rservoir
U=0
Vitesse Uo

Londe arrive au niveau de la vanne et on
retrouve la configuration initiale.

2. Ecoulements surface libre
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
13
2. - ECOULEMENTS A SURFACE LIBRE
2.1. - La surface libre
La surface libre est linterface entre lair et leau. La pression y est gale le plus
souvent la pression atmosphrique.
Les coulements dans les canaux naturels (rivire) et artificiels (irrigation, assainissement)
sont, dans la plupart des cas, des coulements surface libre.


Surface
libre
Seuil


2.2. - Distribution des vitesses dans un canal
La distribution des vitesses dans une section transversale varie en fonction de la
forme de la section et de la rugosit des parois.
Les figures suivantes montrent la rpartition des iso-vitesses dans le cas dun rgime
turbulent. Dans certains cas, la vitesse peut-tre maximale un peu au-dessous de la surface
libre.

0,25

0,5

0,75

1,0

0,85

0,95

1,0

0,9

Vitesse
maximale
max
V
V

Comme pour les coulements en charge, nous avons galement un rgime dcoulement
laminaire et turbulent. (Lcoulement laminaire est rare en hydraulique).


I- Caractristiques des coulements
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
14
2.3. - Notion dinfluence aval
Exprience : londe de gravit
Considrons un canal pente nulle, de section et de hauteur deau constante. Le
fluide est au repos (vitesse nulle U=0). A un instant t, on perturbe la surface libre du canal.
Il se dveloppe deux ondes se dplaant la mme clrit c. Ces ondes sont appeles
ondes de gravit.

c c
U=0 U0
U=0
U0
U=0



Exprience : notion dinfluence aval
On utilise cette fois-ci un canal section transversale, pente, hauteur (h) et dbit
(Q) constants. On cre une perturbation grce une vanne que lon ferme et que lon
ouvre trs rapidement.


Q
U < c
c < 0
c > 0
Aval
Amont

Au niveau de la surface libre, il se cre
deux ondes :
une onde se propageant vers laval,
une onde se propageant vers lamont
quand la vitesse (U) dans le canal
est infrieure la vitesse de londe
de gravit (c).


Q
U > c
c > 0
c > 0
Aval
Amont

Quand la vitesse (U) dans le canal est
suprieure la vitesse de londe de gravit
(c), alors deux ondes se propagent vers
laval avec des vitesses diffrentes.

U : vitesse de lcoulement
c : clrit des ondes
c : vitesse de londe amont
c : vitesse de londe aval

De cette exprience, on en dduit que :
dans le cas o la vitesse du fluide est suprieure la vitesse de londe, lamont
nest pas influenc par les conditions hydrauliques laval (rgime torrentiel);
dans le cas contraire, on a une remonte de londe qui va perturber lamont (rgime
fluvial), ce phnomne est appel influence aval.
2. Ecoulements surface libre
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
15
2.4. - Courbes de remous
Par rapport lcoulement en charge, un coulement surface libre a une difficult
supplmentaire qui est la dtermination de la position de la surface libre par rapport au
fond du canal (tirant deau). Celle-ci est variable en fonction des caractristiques du fluide
et de lcoulement. La figure suivante montre un exemple de forme de surface libre.

Sens de
lcoulement
Seuil

2.5. - Le ressaut
Le ressaut hydraulique se caractrise par une variation importante et croissante de
la hauteur deau de lamont vers laval du phnomne sur une courte distance. Dans la
plupart des cas, une agitation importante de la surface libre permet rapidement de localiser
le phnomne, comme par exemple dans le cas dune ressaut fort.
Le ressaut hydraulique est lun des phnomnes les plus complexes en hydraulique
surface libre. Les connaissances actuelles sur le ressaut hydraulique ne sont pas encore
suffisamment tendues pour que lcoulement interne soit parfaitement compris.

Sens de lcoulement

Ressaut fort

Sens de lcoulement

Ressaut ondul
2.6. - Les seuils et dversoirs
Le seuil cre un obstacle dans un canal, qui oblige le tirant deau augmenter et
donc leau passer par dessus.
Dans le cas dun seuil dnoy, leau chute laval du seuil.
Dans le cas ou le tirant deau laval de louvrage est important, la chute deau ne peut
plus avoir lieu. Dans ces circonstances, le seuil est dit noy.


Seuil dnoy

Seuil trs faiblement noy
I- Caractristiques des coulements
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
16

Seuil faiblement noy

Seuil compltement noy

En fonction de la forme du seuil et de la vitesse de lcoulement, il peut apparatre une
zone de dpression laval du seuil.
Zone de dpression



Le dversoir est un ouvrage de bifurcation qui permet un partage des dbits dans
deux canaux. Par rapport une simple bifurcation, o les dbits sont partags quelle que
soit la hauteur deau, dans un dversoir, le dversement na lieu que si la hauteur du fluide
atteint la hauteur de la crte dversante.

Dversoir latral
Conduite
Amont
Conduite
Aval
Conduite
Dverse
Vue de dessus


AMONT
Milieu naturel
AVAL
Ressaut
hydraulique

3. Proprits des liquides
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
17
3. - PROPRIETES DES LIQUIDES

Dans ltablissement des principes de lhydraulique, certaines proprits des fluides jouent
un rle important, dautres seulement un rle mineur ou aucun rle du tout. En
hydrostatique (fluide au repos) cest le poids spcifique qui est la proprit la plus
importante, tandis quen hydrodynamique (fluide en mouvement), la densit et la viscosit
sont des proprits dominantes. La pression de vapeur prend de limportance quand
interviennent des basses pressions, le liquide en question contient des bulles de vapeur,
cest le phnomne de cavitation. La tension de surface influe sur les conditions statiques
et dynamiques dans les conduits trs troits, cest le phnomne de capillarit.
3.1. - Masse volumique
La masse volumique () est le rapport :
ol
Masse M
Volume V
= = [Kg/m
3
]
Pour les liquides, le volume est pratiquement insensible aux variations de pression et, dans
la majorit des cas, il augmente faiblement quand la temprature augmente, leau faisant
exception cette rgle en dessous de 4C.

eau
= 1000 kg/m
3

mercure
= 13546 kg/m
3

air sec
= 1,205 kg/m
3

Attention : Contrairement aux liquides, les gaz sont fortement compressibles. La variation
de masse volumique dpend de la temprature et de la pression : ( ) f p, T = .
3.2. - Poids spcifique
Il reprsente la force de gravit agissant sur la masse par unit de volume :
= g. [N/m
3
]

eau
= 10
4
N/m
3

3.3. - Viscosit
La viscosit dun fluide en mouvement est la proprit qui exprime sa rsistance une
force tangentielle.
Exprience :
Considrons un fluide plac entre deux plaques planes, parallles, distantes de L et
horizontales. Lune est fixe et lautre est en mouvement uniforme de vitesse Uo. Pour
gnrer une vitesse de la plaque suprieure (surface A), il faut exercer une force F.

Cette force est la rsultante des forces de frottements visqueux.

0
L
y

Rpartition de la vitesse
entre deux plaques en
rgime laminaire

Uo Uo

I- Caractristiques des coulements
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
18
Lexprimentation permet de dduire une proportionnalit entre le rapport de la force F et
la surface A avec le rapport entre la vitesse Uo et la longueur L telle que :
0
U F U
A L y

= =


[N.s/m
2
] est appel viscosit dynamique ou absolue.

Le viscosimtre :
On considre deux cylindres coaxiaux spars par un intervalle e dont lespace entre eux
est rempli par un liquide. On fait tourner le cylindre extrieur vitesse constante () et on
maintient fixe le cylindre intrieur.

Niveau
de leau
Cylindre en
rotation
Cylindre fixe
Vitesse de
rotation :
e
r
h
Variation de la
vitesse de du
liquide suivant
lpaisseur e

Le fluide en contact avec le cylindre extrieur va y adhrer et par consquent va tre
anim de la vitesse V du cylindre extrieur. Le fluide en contact avec le cylindre fixe aura
une vitesse nulle. La viscosit fait natre une force de frottement que lon mesure par le
couple M.
Les expriences ont montr que :
si e est faible par rapport au rayon intrieur r, la courbe reprsentative de la
variation de la vitesse entre r et r+e est une droite,
le couple (M) varie proportionnellement la vitesse et on a :
( ) 2. .r.h .r
M V
e

=
On dfinit un deuxime coefficient de viscosit, le coefficient de viscosit cinmatique :

= [m
2
/s]

eau 20C
= 10
-3
N.s/m
2

eau 20C
= 10
-6
m
2
/s

mercure
= 1,554.10
-3
N.s/m
2

mercure
= 0,1147.10
-6
m
2
/s

air
= 18,5.10
-6
N.s/m
2

air
= 15,6.10
-6
m
2
/s

On appelle fluide parfait un fluide dont la viscosit serait nulle (fluide inexistant dans la
nature). La viscosit existe ds quil y a mouvement relatif entre particules, que ce soit en
rgime laminaire ou turbulent.
3. Proprits des liquides
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
19
3.4. - Pression de vapeur saturante
Lbullition est un phnomne de changement dtat, dans lequel le liquide passe
ltat de vapeur. Tous les liquides ont tendance svaporer ; la phase liquide se
transforme en phase gazeuse. Au cours de cette transformation, les molcules de vapeur
exercent une pression appele pression de vapeur saturante. Dans le cas de leau, la
pression de vapeur (p
s
) crot avec une augmentation de la temprature (T).
La pression de vapeur saturante pour leau est donne par la relation empirique suivante :
( ) ( )
10 s 10
2795
log p 22.435 3.868log T 273.15
T 273.15
= +
+
avec p
s
en Pa et T en Celsius.

Si, temprature constante, on abaisse la pression la surface dun liquide, ce dernier se
met bouillir lorsquon atteint la pression de vapeur saturante correspondant cette
temprature. Dans lcoulement des liquides, il peut arriver que la pression en certains
points devienne infrieure la pression de vapeur saturante. Le liquide entre alors
localement en bullition et des bulles de vapeur apparaissent au sein mme de
lcoulement. Ce phnomne, appel cavitation, est le plus souvent nuisible pour les
installations o il se produit (canalisation, pompes, turbine). Les variations de volume
lors du changement dtat sont telles quil se produit au sein du fluide de vritables
explosions de bulles au moment de la vaporisation et de violentes implosions, lors de la
condensation.

3.5. - Tension superficielle

Une molcule liquide au repos est soumise aux forces dattractions que les molcules
voisines exercent sur elle. Une molcule la surface libre dun liquide ou la surface de
sparation de deux liquides non miscibles nest plus soumise laction de forces
symtriques, puisquelle nest plus entoure symtriquement par dautres molcules de
mme nature. Ainsi la rsultante des forces molculaires nest plus nulle. La surface de
sparation se comporte comme une membrane tendue.
La force dattraction tangentielle la surface ncessaire pour arracher des particules
agissant le long dun segment de longueur unitaire est appele tension superficielle.

Forces entre molcules
au sein du fluide
Forces entre molcules
au niveau de la surface libre
Surface libre
Molcules

Les effets de tension superficielle ne sont pas importants dans les coulements en eau
potable ou en assainissement et ne sont donc pas pris en compte.
Tension de surface : Air eau 20C : 0,0724 N/m.
I- Caractristiques des coulements
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
20
4. - CARACTERISATION DES FORCES DANS UN ECOULEMENT
4.1. - Les forces
Les forces qui agissent sur un volume fini de fluide sont de deux types :
- Les forces de volumes,
- Les forces de surfaces.
4.1.1. - Les forces de volumes
Elles se composent des forces suivantes :
- Les forces de pesanteur provenant de la gravit.
- Les forces dacclration pure :
Elles proviennent de la variation de la vitesse (V) de la masse dune fluide
(M) dans le temps.
acclration pure
V
F M
t


Prenons par exemple deux rservoirs la mme hauteur, dont lun est vide et
lautre plein, relis par une conduite de diamtre constant, horizontal et muni
dune vanne. A louverture de la vanne, il se produit un coulement. La
variation de la vitesse dans le temps ( ) V t cre au sein de lcoulement
une force dacclration pure.
- Les forces dacclration convective :
Elles proviennent de la variation de la vitesse (V
x
, V
y
, V
z
) dans lespace
(repre [x, y, z]).
acclration convective x y z
V V V
F M .V .V .V
x y z
| |
= + +
|

\

Prenons une conduite dont lcoulement ne varie pas dans le temps
( 0 t / V = ). Lcoulement tant permanent, le dbit est identique en tout
point de la canalisation. Or, si la surface A est suprieure la surface B alors
la vitesse en B (V
b
) est suprieur la vitesse en A (V
a
).
Cette variation de vitesse va engendrer une acclration qui va gnrer une
force dacclration convective.

V
a
V
b

V
d

V
c

Rgime
acclr
Rgime
acclr
Rgime
uniforme
L
ab

Convergent
L
cd

Divergent
L
bc

A B C D
a b
ab
V V
0
L

>
b c
bc
V V
0
L

=
c d
cd
V V
0
L

<


4. caractrisation des forces dans un coulement
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
21
4.1.2. - Les forces de surfaces
Elles se composent des forces suivantes :
- Les forces de pression :
La pression (p) est le rapport entre une force F agissant perpendiculairement
la surface (A) dun fluide :
F
p
A
=
- Les forces de frottement de viscosit :
Nous avons vu prcdemment quun fluide, dont les particules sont en
mouvement relatif, gnrent des forces de frottement dues la viscosit. Le
force de frottement scrit :
U
F A
y


- Les forces gnres par la turbulence :
La turbulence dcrite au premier chapitre joue un rle majeur dans
lcoulement des fluides. La turbulence tendance freiner
lcoulement. Une faon de les reprsenter mathmatiquement consiste les
assimiler des forces de frottement, ce qui est faux compte tenu de la nature
mme de la turbulence.
4.2. - Importance des diffrentes forces
En gnral, lhydraulicien doit soccuper de leffet de la force dominante. Dans la
plupart des problmes dcoulement des fluides, la pesanteur, la viscosit et llasticit
sont prpondrantes, mais pas toujours simultanment.
Limportance relative des diffrentes forces agissant sur un liquide est calcule par des
nombres adimensionnels reprsentant les rapports entre ces forces. Lanalyse
dimensionnelle permet de simplifier ces rapports. Les diffrentes dimensions utilises
sont :
L : longueur,
T : temps,
: masse volumique.
4.2.1. - Le Nombre dEuler
Le nombre de Euler est le rapport entre les forces dinertie (Ma) et les forces de pression
(pA).
( )
p
V
pL
) T / L ( L
pA
Ma
2
2
2 3
=

=
M : la masse,
a : lacclration du fluide,
p : la pression,
A : la surface dapplication de la pression,
V : vitesse caractristique de lcoulement.
Ce rapport est utilis lorsquon calcul les forces (statiques et dynamiques) de leau sur une
paroi ou un ouvrage.
4.2.2. - Le Nombre de Reynolds
Le Nombre de Reynolds (Re) est le rapport entre les forces dinertie (Ma) et les forces de
viscosit (A).
Re=
( )

=
|

\
|

VL
L
L
T / L
) T / L ( L
A
Ma
2
2 3

: la force de frottement par unit de surface,
A : la surface dapplication de la force de frottement,
L : longueur caractristique de lcoulement.

I- Caractristiques des coulements
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
22
Ce rapport permet de distinguer le rgime laminaire et turbulent.
Si Re > valeur dpendant de louvrage tudi alors turbulent
Si Re < valeur dpendant de louvrage tudi alors laminaire

Exemples :
Ecoulement dans une
conduite circulaire :
Si Re > 2500 alors turbulent
Si Re < 2000 alors laminaire
L : diamtre de la conduite.
Ecoulement entre deux
plaques parallles :
Si Re < 500 alors laminaire
Sinon turbulent
L : distance entre les deux
plaques.
Ecoulement autour dune
sphre :
Si Re < 1 alors laminaire
Sinon turbulent
L : diamtre de la sphre.
4.2.3. - Le Nombre de Froude
Le nombre de Froude est le rapport entre les forces dinertie (Ma) et les forces de
pesanteur (Mg).
( )
( ) ( ) Lg
V
g L
V L
g L
) T / L ( L
Mg
Ma
2
3
2 2
3
2 3
=

=
g : lacclration de la pesanteur,
L : longueur caractristique de lcoulement,
V : vitesse caractristique de lcoulement.
Ce rapport permet de distinguer le rgime fluvial ou torrentiel que nous avons vu au 2.3.
4.2.4. - Le Nombre de Cauchy
Le nombre de Cauchy est le rapport entre les forces dinertie (Ma) et les forces dlasticit
(EA).
( )
E
V
EL
) T / L ( L
EA
Ma
2
2
2 3

=
E : module dlasticit.

4.2.5. - Le Nombre de Weber
Le nombre de Weber est le rapport entre les forces dinertie (Ma) et les forces de tension
superficielle (L).
( )

2 2 3
LV
L
) T / L ( L
L
Ma

: tension de surface.

5. - LES EQUATIONS DE BASE
La dtermination des caractristiques dun coulement consiste rechercher la
pression et la vitesse en tous points. Pour cela, il faut crire des quations dquilibre (ou
de conservation) entre les diffrentes forces agissant sur le fluide.
Ces quations dquilibre sont :
- lquation de continuit ou quation de conservation de la masse (m) du fluide.
dm
0
dt
=

- lquation de la quantit de mouvement qui traduit lquilibre entre la somme
des forces extrieures (F) qui exercent une influence sur le fluide et le taux de
variation de la quantit de mouvement du fluide de masse m.
d m v
dt
F
( . )
=


Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
6. - OUTILS MATHEMATIQUES

Elment de volume : dv = dx.dy.dz


dx
dy
dz

Drive partielle :
x


Drive totale : dz
z
P
dy
y
P
dx
x
P
dt
t
P
dP

=
dP P P x P y P z
dt t x t y t z t

= + + +


Gradient dun scalaire :
(
(
(
(
(
(

= =

z
f
y
f
x
f
) f ( grad ) f ( grad
Gradient dun vecteur :
(
(
(
(
(
(
(

=
(
(
(

=
z
Vz
y
Vz
x
Vz
z
Vy
y
Vy
x
Vy
z
Vx
y
Vx
x
Vx
Vz
Vy
Vx
grad ) V ( grad
Divergence dun vecteur :
z
Vz
y
Vy
x
Vx
) V ( div

=
Rotationnel : ( )
(
(
(
(
(
(

=
(
(
(

(
(
(
(
(
(

=
y
Vx
x
Vy
x
Vz
z
Vx
z
Vy
y
Vz
Vz
Vy
Vx
z
y
x
V Rot

II- Hydrostatique
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
24
Chapitre II :
HYDROSTATIQUE

Lhydrostatique tudie les conditions dquilibre des liquides au repos. Ce chapitre
aborde ltude de la rpartition de la pression, notamment en fonction de la distance
verticale, ainsi que les forces qui en rsultent.
1. - EQUATIONS DE LHYDROSTATIQUE
Nous avons vu au chapitre prcdent que lensemble des forces agissant sur un
fluide sont de deux natures : les forces de volume et les forces de surface. En reprenant les
rsultats acquis dans ce chapitre, nous allons dtaill lensemble de ses forces.
Considrons dans un rservoir un fluide au repos, dont on extrait un petit paralllpipde
deau daxe vertical z. Soit p la pression en son centre. Il est soumis aux forces verticales
suivantes :

p(x,y,z)
x
y
z
dz
Rservoir
p dz
p dx.dy
z 2
| |

\

p dz
p dx.dy
z 2
| |
+
|

\

.g.dx.dy.dz

Concernant les forces de volume, il nen existe quune seule la force de pesanteur. Elle
scrit de la faon suivante :
( )
pesanteur
F g. dx.dy.dz =
Les forces dinertie nexiste pas puisque le fluide est au repos (vitesse nulle).

Concernant les forces de surface, la pression agit sur la face suprieure et infrieure de
llment. Ces forces scrivent de la faon suivante :
Force de pression sur le surface infrieure :
pression inf.
p dz
F p dx.dy
z 2

| |
=
|

\

Force de pression sur le surface suprieure :
pression sup.
p dz
F p dx.dy
z 2

| |
= +
|

\


Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Les forces de viscosit et de turbulence nexistent pas puisquil ny a pas de vitesse
relative entre les particules de fluide.

Lquation de lhydrostatique est dtermin en crivant lquilibre de lensemble des
forces : inertie d' Forces F =

. En projetant cette quation suivant la verticale, on a :



p dz p dz
.g.(dx.dy.dz) p dx.dy p dx.dy 0
z 2 z 2

| | | |
+ + =
| |

\ \

soit :
p
- .g 0
z


On peut crire de faon analogue les quations dquilibre dans les autres directions :
p p
0 ; 0
x y

= =


Ces trois quations montrent que la pression est indpendante de x et de y, cest--dire que
la pression ne varie pas dans les directions x et y ou encore quelle est constante dans un
plan horizontale. Cela est vrifi tant que lon reste dans un mme fluide ( constante). La
pression ne dpend que de z, ce quon crit :
p
- .g ou dp=- .g.dz
z




2. - VARIATION DE LA PRESSION DANS UN FLUIDE
INCOMPRESSIBLE
Pour un fluide incompressible ( = Cte), lintgration de lquation ci-dessus
scrit : p + .g.z = Cte.

On crit frquemment : p* = p + .g.z = Cte.
On nomme p* lnergie potentielle par unit de volume.
Une autre criture consiste diviser lquation par .g : Cste z
g
p
g
* p
= +

,
On remarque que
g
* p

est homogne z cest--dire une longueur ce qui offre un moyen


pratique pour la reprsenter graphiquement.
g
* p

est appele hauteur pizomtrique.


Dans ce cas, la pression p est mesure en hauteur de colonne du liquide .

Units de mesure :
Lunit lgale est le pascal : 1 Pa = 1N/m
2
Un multiple du pascal est le bar : 1 bar = 10
5
Pa
Il existe dautres units plus pratiques :
- Le mtre de colonne deau :
3
eau
1, 0mCE 1, 0mCE. .g 9,81.10 Pa 0,098bar = =
5
5
eau
10
1, 0bar 10 Pa 10,19mCE
.g
= =


II- Hydrostatique
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
26
- Le millimtre de mercure :
3 2
mercure
1, 0mmHg 10 .g 1, 33.10 Pa

=
5
5
mercure
10
1.0bar 10 Pa .1000 752,5mmHg
.g
= =


En hydraulique des rseaux urbains, la pression est le plus souvent exprime en mtres de
colonne deau et en bars.

Pression absolue et pression relative :
La pression absolue est dfinie par rapport la pression dans le vide qui correspond la
pression nulle. On en dduit donc que la pression minimale possible est zro.
La pression relative se dfinit par rapport une rfrence que lon choisi le plus souvent
gale la pression atmosphrique. Cela consiste finalement faire une translation du
repre des pressions. La pression nulle est donc quivalente la pression atmosphrique
(p
a
). La pression minimale correspond donc : -pa (pression atmosphrique ngative).

Prenons par exemple un rservoir o la surface libre est la pression atmosphrique (p
a
).
En crivant lquation de lhydrostatique par rapport un plan de rfrence, on a :
p + .g.z = Cte
Entre le point 1 et 2, on a :
p
1
+ .g.z
1
= p
2
+ .g.z
2
= Cte,
soit : p
1
= p
2
+ .g.(z
2
-z
1
) = p
a
+ .g.h, avec p
a
=10
5
Pa.
La pression en 1 est mesure en pression absolue.
La figure suivante montre la rpartition des pressions suivant la hauteur.

Pression
atmosphrique : Pa
1
2
h
+ =
p
1
=gh p
a
p
a
+p
1
Plan de rfrence
z
1
z
2


Dans la pratique, on prfre souvent mesurer par rapport la pression atmosphrique (p
a
).
Dans ce cas, la pression au point 1 devient : p
1
= p
1
- p
a
= .g.(z
2
-z
1
) = .g.h. La pression
en 1 est mesure en pression relative par rapport la pression atmosphrique.
La plupart des instruments de mesure fournit une pression relative appele galement
pression manomtrique.
4. Forces hydrostatiques sur les parois
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
27

3. - VARIATION DE LA PRESSION DANS UN FLUIDE
COMPRESSIBLE
Contrairement aux liquides, les gaz sont fortement compressibles. La variation de la
masse volumique dpend de pression p et de la temprature T : ( ) f p, T = . Sil sagit
dun gaz parfait, nous avons :
p
rT =

, r tant une constante. On peut ainsi dfinir


laltitude z
0
:
0
0
0
p
rT =

, ce qui nous permet de calculer la constante r.


En reprenant lquation de lhydrostatique :
dp
- .g
dz
= on a :
0
0
0
T dp p
g
dz p T
= .
Si on se place temprature constante T=T
0
, lquation peut ainsi scrire :
0 0
0
0 0
g g dp
dz => p=p exp z
p p p
| |
=
|
\


Dans le cas de lair, p
0
=1,013 10
5
Pa et
0
=1.205 kg/m
3
, lquation prcdente permet de
caractriser lvolution de la pression dans lair en fonction de laltitude :
5 0
0
0
g z
p=p exp z 1, 013.10 exp
p 8570
| | | |
=
| |
\
\


4. - FORCES HYDROSTATIQUES SUR LES PAROIS

Les forces hydrostatiques sur une surface proviennent des forces de pressions du
fluide agissant sur cette surface.
Il convient, donc dans un premier temps, de caractriser la pression du fluide sur une
surface. Pour cela, on a besoin de :
- lintensit : la pression dpend de la profondeur deau h. Elle est calcule par la
relation : p = .g.h,
- la zone dapplication : la pression sapplique sur une surface (ds),
- la direction : la pression est toujours perpendiculaire la surface dapplication.


p
h
ds
surface gauche
surface dapplication
de la pression

Le calcul des forces hydrostatiques sur une surface quelconque plonge dans leau,
consiste dterminer les trois caractristiques suivantes :
II- Hydrostatique
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
28
- lintensit de la force qui sapplique sur la surface ds : dF = p.ds = .g.h.ds,
- le point dapplication de la force,
- la direction.
Dans le but de fournir des rsultats facilement applicables, on distingue les surfaces planes
et les surfaces gauches.
4.1. - Paroi plane en position incline
On sintresse aux surfaces planes de forme quelconque entirement immerge dans leau.
La figure suivante reprsente gauche la surface immerge et droite une vue A-A de
cette surface. On dfinit un repre (x,y) dont laxe (x) est sur la surface libre et (y) dirig
vers le bas et passant par la surface plane. Le point G(x
G
,y
G
) est le centre de gravit de la
section. On dfinit le repre (, ) comme tant une translation du repre (x,y) centr en G.
Lintensit de la force rsultante agissant sur la surface S est dfinie par :

= =
S S
hds g . dF F
Lintgration de cette quation scrit :
F = .g.h
G
.S,
h
G
: hauteur deau du centre de gravit de la paroi immerge,
S : surface de la paroi immerge.

G
P
y
x
y
P
y
G
y


h
P
h
G
h=y.sin()

y
x
ds
F

Vue A-A

I
xx
I

Vue A-A

Surface S
ds

Le point dapplication de la force rsultante des pressions P(x
p
, y
p
) est appel : centre de
pression ou de pousse.
La position de ce point est dfinie par la position du barycentre des surfaces lmentaires
(ds) pondres par la pression sur chaque surface, ce qui revient calculer le moment
quivalent des forces de pression, cest--dire :
4. Forces hydrostatiques sur les parois
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
29
p
S
p
S
x.dF x .F
y.dF y .F
=
=


Dans le grande majorit des cas les surfaces sont symtriques par rapport laxe , ce qui
revient dire que : x
p
= x
G
.
La deuxime intgrale scrit :
2
S
p G
G
S
y ds
I
y y
y .S
yds

= = +

reprsente linertie de la section suivant les axes .


Le tableau suivant fournit le centre de gravit, la surface et linertie pour quelques formes
de surface plane.

h
v
v
b
G


3
2h h
v ; v
3 3
bh bh
S ; I
2 36

= =
= =



h
v
v
b
G


3
h h
v ; v
2 2
bh
S bh ; I
12

= =
= =


h
v
v
B
G

b

( )
( )
( )
3 2 2
h 2B b h B 2b
v ; v
3 B b 3 B b
h B 4Bb b
h
S B b ; I
2 36 B b

+ +
| | | |
= =
| |
+ +
\ \
+ +
= + =
+



d
v
v
a
G


2 2 4 4
d a 3 d a 3
v ; v
2 2 2 2
d 3 3a 3 5d 5a 3
S ; I
2 2 16 48 3

= = = =
= = = =

II- Hydrostatique
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
30


v
v
R
G


( )
( ) ( )
( ) ( )
( ) ( )
( ) ( )
( ) ( )
( ) ( )
3
2
3
4 4
4sin
v R 1
3 2 sin 2
v R 1 cos v
R
S 2 sin 2
2
1 cos 2
R R
I 4 sin 4
16 9 2 sin 2

| |

= |
|

\
=
=

=




v
v
G


4
2
v R ; v R
R
S R ; I
4

= =

= =


4.2. - Paroi surface gauche
La premire mthode permettant le calcul de la force rsultante F sur une paroi gauche
consiste dcomposer la force lmentaire dF suivant les axes x et y :
x y
dF dF .x dF .y = +
Lvaluation des projections de dF suivant x et y, peut se faire de la manire suivante :
x
y
dF cos( ).p.ds cos( ). .g.h.ds
dF sin( ).p.ds sin( ). .g.h.ds
= =
= =


1
2
3
h
1
h
2
dF

x

y

dF
y
dF
x



Lintgration
x
dF et
y
dF sur toute la surface de llment courbe permet dvaluer le force
rsultante F. Compte tenu de la surface courbe, langle est variable, ce qui complique le
calcul de lintgrale.
4. Forces hydrostatiques sur les parois
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
31

Une deuxime mthode consiste isoler un volume de fluide et faire lquilibre des
forces extrieures agissant sur ce volume.
Dans lexemple suivant, le volume de fluide isol est compos dun ensemble de surfaces
planes horizontales et verticales et de la surface gauche (23). Le choix des surfaces planes
se justifie par lutilisation des relations prcdentes.

1
2
3
h
1
h
2
p
1
p
2

p
3

p
3

W
F
x
F
v
F
R
x

y


En faisant lquilibre des forces suivant lhorizontale, on en dduit que la composante
horizontale de la force hydrostatique (F
R
) est donne par F
x
. La composante verticale est la
somme de la force F
v
et du poids de leau W. En faisant la composition vectorielle des
forces F
x
et F
v
, on en dduit F
R
. Il suffit dcrire ensuite le moment des forces par rapport
un point quelconque pour localiser la position de F
R
.

5. - FORCES HYDROSTATIQUES SUR DES CORPS IMMERGES
5.1. - Forces dArchimde
Supposons quune surface ferme formant un corps solide de poids volumique
s
, de
volume total V et de volume immerg V
immerg
, se trouve immerge entirement ou
partiellement (V
immerg
V) dans un liquide au repos de poids volumique . Les forces
verticales qui agissent sur le corps sont :
- les forces de pesanteur : V.
s

- les forces de pression du liquide : V
immerg
.
(Appeles forces dArchimde ou forces de portance)
La force dArchimde est applique au centre de gravit du liquide dplac (centre de
pousse). Pour un corps plein de poids volumique homogne et entirement immerg, le
centre de pousse est confondu avec le centre de gravit du corps. Il nen est pas de mme
pour les corps flottants.

5.2. - Equilibre des corps immergs
Un corps est en quilibre si le poids W et la force dArchimde sont gaux, opposs et
situs sur la mme ligne verticale. Dans le cas contraire, il en rsulte un mouvement.
La stabilit peut se dfinir de la faon suivante : si on incline un corps dun angle par
rapport la verticale, le corps est soumis un couple de redressements qui le fait tourner
jusqu ce quil revienne sa position initiale.
II- Hydrostatique
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
32
Linstabilit est donc dfinie par un couple qui tend augmenter linclinaison.


G
P
G
P
G
P Dsquilibre
volontaire
Mouvement
naturel de la
quille
G
P
Mouvement
naturel de la
quille
Dsquilibre
volontaire
Equilibre stable
Equilibre instable



Dsquilibre
volontaire
Mouvement
naturel
Dsquilibre
volontaire
Mouvement
naturel

Dans cet exemple, on constate que la position dquilibre stable est vrifie pour un angle
, et que la position instable correspond un angle qui a la particularit : < .
On en dduit quil existe donc un angle limite
limite
de basculement entre ltat stable et
instable.



Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Chapitre III :
LHYDRAULIQUE EN CHARGE



La dynamique des fluides consiste tudier le mouvement des particules fluides soumises
un systme de forces. Bien souvent, on commence par ltude les fluides fictifs dit
fluides parfaits . Ils ont la particularit de ne pas avoir de viscosit et de ne pas
dvelopper de la turbulence. Ils permettent dtablir lquation de Bernoulli facilement.
Les fluides rels engendrent des forces de frottement dues la turbulence et la viscosit.
La prsence de ces forces induit une perte de charge (nergie) qui est une transformation
irrversible de lnergie mcanique en nergie thermique.
Ce chapitre aborde, dans un premier temps, lquation de continuit et surtout
ltablissement de lquation de Bernoulli. Dans un deuxime temps, lvaluation des
pertes de charges ainsi que les mthodes de calcul des rseaux hydrauliques en charge sont
exposes.

1. - LEQUATION DE CONTINUITE
Cette quation exprime le principe de conservation de la masse : la variation de
masse de fluide dun lment de volume dv pendant un temps dt est gale la masse de
fluide entrante dans ce volume dduite de la masse de fluide sortante.


dx
dz
dy
x
y
z
(x,y,z)
Variation de masse entre t et t+dt
La masse de fluide contenue dans le volume dv=dx.dy.dz est gale au temps t :
.dx.dy.dz
34 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Aprs un temps dt dans ce mme volume, la masse est gale : dz . dy . dx dt
t
|

\
|

+
On constate donc une variation de masse de : dz . dy . dx . dt
t



Variation de masse de fluide pendant la dure dt entre ce qui entre et ce qui sort :
On dfini le dbit massique par :
massique suivant x volumique suivant x
q .Q .u.dy.dz = =

massique suivant x
q
( )
massique suivant x
massique suivant x
q
dx
q .
x 2
| |

|
|
\

1 2
( )
massique suivant x
massique suivant x
q
dx
q .
x 2
| |

| +
|
\

La masse de fluide entrante par la face 1 (suivant x) pendant le temps dt est :
( ) ( ) massique suivant x
massique suivant x
q
u.dydz
dx dx
q . .dt u.dydz . .dt
x 2 x 2
| |

| |
| =
|
|

\
\

La masse de fluide sortante par la face 2 (suivant x) pendant le temps dt est :
( ) ( ) massique suivant x
massique suivant x
q
u.dydz
dx dx
q . .dt u.dydz . .dt
x 2 x 2
| |

| |
| + = +
|
|

\
\

La diffrence de masse entre ce qui entre et ce qui sort dans le volume dv est donc suivant
x : dt . dz . dy . dx
x
u


On a de mme suivant y et z
u
avec V v
w
(
(
=
(
(

: dt . dz . dy . dx
y
v

; dt . dz . dy . dx
z
w


En crivant que la variation de masse de fluide dun lment de volume dv pendant un
temps dt est gale la masse de fluide entrante dans ce volume moins la masse de fluide
sortante, on a :
( u) ( v) ( w)

t x y z
div( V) 0
t

=

+ =


Cest lquation de continuit dun fluide conservatif.


Cas des fluides incompressibles :
Si le fluide est en mouvement permanent ( 0
t
=

) la masse volume est indpendante du


temps et si le fluide est incompressible, est indpendant de x, y et z, lquation de
continuit se rduit :
div(V) = 0

1-Equation de Continuit
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
35
Linterprtation physique de cette quation est la suivante : les dbits Q entrant et sortant
travers un volume quelconque et rempli du fluide doivent tre gaux.

Exemple de conduite section variable :
(Intgration de div(V) sur le volume et passage une intgrale de surface par Ostogradski)


U
1
.S
1
= U
2
.S
2



2. - EQUATIONS DE BERNOULLI
2.1. - Cas des fluides parfaits
Les fluides parfaits ont donc la particularit de ne pas avoir de viscosit et de ne pas
dvelopper de la turbulence. Ils permettent dtablir lquation de Bernoulli facilement.
2.1.1. - Equations dEuler

p(x,y,z)
p dz
p dx.dy
z 2

| |

\

p dz
p dx.dy
z 2

| |
+
|

\

.g.dx.dy.dz


En raisonnant, dans un premier temps, suivant la verticale (z), les forces qui agissent sur
cet lment de volume dv = dS.dz, sont :
Les forces de volumes :
- les forces de pesanteur provenant de la gravit : - .g.dv
36 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
- Les forces dinertie : dv
t
z
z
w
t
y
y
w
t
x
x
w
t
w
dv
dt
dw
|
|

\
|

=
- les forces dacclration pure :
w
dv
t


- les forces dacclration convective : .grad(w).V.dv
Les forces de surfaces :
- les forces de pression : Surface infrieure :
p dz
p dx.dy
z 2
| |

\

: Surface suprieure :
p dz
p dx.dy
z 2
| |
+
|

\

- les forces de frottement de viscosit : 0
- les forces de frottement de turbulence : 0

En crivant lquilibre de lensemble des forces : inertie d' Forces F =


p dz p dz w
.g.(dx.dy.dz) p dx.dy p dx.dy dS.dz .grad(w).V.dS.dz
z 2 z 2 t
p w
.g .grad(w).V
z t
| | | |
+ + = +
| |

\ \

= +



En crivant lquation suivant les trois directions :
1 p u u u u
u v w
x t x y z
1 p v v v v
u v w
y t x y z
1 p w w w w
g u v w
z t x y z

= + + +


= + + +


= + + +



soit : 0 ) z . g . p ( grad
1
V ). V ( grad
t
V
= +

+ +


Cette quation est appele lquation dEuler.


Recherche dune solution
Les caractristiques du mouvement dun fluide incompressible dpendent de quatre
inconnues : V
x
, V
y
, V
z
et de p.
Lquation dEuler fournit trois quations.
Lquation de conservation de la masse fournit une quation : div(V) = 0
Pour obtenir une solution, il faut galement dfinir un nombre suffisant de conditions
aux limites :
- la surface libre : p(x, y, z, t) = Cste = p
a
.
- la paroi : V doit tre parallle la paroi.
2-Equations de Bernoulli
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
37
Utilisation des quations dEuler pour le calcul des forces hydrauliques sur une
surface

Le thorme dEuler permet de calculer la raction de leau sur un lment. En raisonnant,
suivant un tube de courant en rgime permanent, les forces qui agissent sur cet lment
sont :

- Les forces de volumes : - les forces de pesanteur provenant de la gravit : g . V .
vol

- les forces dinertie :
- les forces dacclration pure : 0
- les forces dacclration convective :


= =

j
j j ext j j
V S
ext
S ). n . V .( V . ds ). n . V .( V . dv . V
s
V
.
- Les forces de surfaces : - les forces de pression sur les surfaces S
i
:

i
i
i
S . p
- les forces de frottement de viscosit : 0
- les forces de frottement de turbulence : 0
- laction du fluide sur llment : R

En crivant lquilibre de lensemble des forces : inertie d' Forces F =




= +
j
j j ext j j
i
i
i
vol
S ). n . V .( V . R S . p g . V .



2 2 2 1 1 1 2 ext 2 2 1 ext 1 1 vol
V . S . V . V . S . V . R n . S . p n . S . p g . V . + =

38 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
2.1.2. - Equations de Bernoulli
Le systme dquations intrinsques consiste crire les quations dEuler en rgime
permanent ( 0 t = ) dans un repre particulier. Ce repre est constitu par les lignes de
courant pour le vecteur t et par le vecteur n tel que v n.

en appelant s

le vecteur unitaire de la tangente la trajectoire, on a :


s V V

= et
dt
s d
V s
dt
dV
dt
V d

+ =
avec : V
R
n
dt
ds
.
ds
s d
dt
s d

= =
R : rayon de courbure et n

le vecteur perpendiculaire s

.

n suivant ) p h . g . (
n
1
R
V
. V
s suivant ) p h . g . (
s
1
s
V
V

=
+



Avec, V : lintensit de la vitesse dune particule fluide.

En restant sur la ligne de courant on peut intgrer la premire quation et lon
obtient lquation de Bernoulli :
Cste p h . g .
2
V
2
= + +
Lquation de Bernoulli est valable en tout point du fluide incompressible en
mouvement permanent.

En se plaant sur une ligne normale la ligne de courant et un rayon de courbure
infini (ligne droite), on obtient :
Cste p h . g . = +
2-Equations de Bernoulli
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
39
2.1.3. - Reprsentation graphique
Les diffrents termes de lquation de Bernoulli peuvent tre crits en quantits
homognes des hauteurs de liquide de poids volumique :
g
E
mg
W
= Cste
p

p
h
g 2
V
t
2
= =

+ +

hauteur due cote du hauteur due charge totale
la vitesse point la pression

Hauteur pizomtrique




En terme de pression :
volume
t
2
V
W
= Cste p p .g.h
2
V
= = + +
Pression pression statique pression totale
dynamique
2.1.4. - Interprtation nergtique

Lquation de Bernoulli peut donner lieu une interprtation nergtique :
Cste p * p
2
V
t
2
= = +

Energie Energie Energie
cintique potentielle mcanique totale


40 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Lquation mcanique totale contenue dans un volume unitaire est donc le travail
mcanique total que la particule est susceptible daccomplir. Ainsi, on peut dire que
lquation de Bernoulli traduit la conservation de lnergie mcanique totale par unit de
volume au cours du mouvement permanent.
Dans un coulement, lnergie mcanique totale par unit de volume de fluide, peut tre
modifie dune section lautre en introduisant par exemple dans le circuit une machine
hydraulique.
Ainsi lexpression de lquation de Bernoulli scrit :
g
E p
+ h
2g
V
=
p
+ h
2g
V
2
2
2
2 1
1
2
1

+
Si lchange dnergie se fait des parois de la machine vers le fluide nous avons affaire
une pompe, si au contraire, lchange dnergie se fait du fluide au parois de la machine,
nous avons affaire une turbine.
Lquation dnergie est modifie par le terme E, qui reprsente laugmentation par une
pompe ou la diminution par une turbine de lnergie mcanique totale par unit de masse
de liquide en mouvement.
2.1.4.1. - Turbine

Exemple dinstallation :


La puissance est dfinie par :
Temps
Travail
Puissance =
2-Equations de Bernoulli
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
41

Pour une turbine la puissance hydraulique est dissipe par :
[ ]
2
2
. *
2
Variation d'nergie
Unit de temps
.
* . .
2
t
V
Volume p
E
P
t t
Volume
V
p g Q H
t

( | |
+
( |

\
= = =

| |
= + =
|

\

( ) H . Q .
g
E
. Q . P
t
t
=
|
|

\
|
=


Types de turbines :

Turbine Pelton

Turbine Francis

42 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Turbine Kaplan


Choix dune turbine :


2-Equations de Bernoulli
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
43
2.1.4.2. - Pompes
Exemple dinstallation :


Principe de montage :



Types de pompes :

44 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Choix dune pompe :


La puissance est dfinie par :
Temps
Travail
Puissance =


Pour une pompe la puissance hydraulique est fournie par : ( ) H . Q .
g
E
. Q . P
p
p
=
|
|

\
|
=

2-Equations de Bernoulli
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
45
Exemple de courbe caractristique pour une pompe multicellulaire tage (de 1 12
tages) :
H= Charge totale en sortie Charge totale en entre

Caractristiques par tage :

46 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Dfinition du NPSH Net Positive Suction Head (charge nette laspiration)

La pompe transmet une vitesse au fluide par lintermdiaire de la roue ou de lhlice. La
vitesse du fluide qui entre dans la roue augmente et par consquent la pression diminue
lentre engendrant une aspiration et le maintien de lamorage. La rotation lintrieur de
la pompe augmente la vitesse du fluide tandis que la force centrifuge qui le comprime sur
la priphrie augmente sa pression. Lorsqu la suite de survitesses locales des zones de
sous-pression apparaissent sur les bords dattaque des aubes des roues de pompes
centrifuges, la pression descend en dessous de la pression de vapeur saturante du liquide,
ce phnomne entrane la formation de bulles de vapeur. Leau utilise dans les
installations industrielles contient de lair dissous dans une proportion atteignant 25mg par
litre. Lorsque leau est soumise une pression dcroissante, on observe le dgagement de
lair dissous quand la pression atteint la pression de saturation et une vaporisation quand la
pression atteint la pression de vapeur. Un phnomne de dgazage puis de vaporisation
partielle du liquide vont apparatre lintrieur de la pompe. Les micro-bulles formes
sont transportes par le fluide dans les zones de plus haute pression o elles implosent au
contact des pices mobiles. Un tel phnomne saccompagne de variations normes de
pression (plusieurs milliers de bars) localises sur de trs petites surfaces. Il en rsulte une
rosion intense des surfaces mtalliques et une chute de pression dans toute linstallation.

Le critre qui sert viter la cavitation dans une pompe est le NPSH.
NPSH= Charge totale lentre de la pompe pression de vapeur saturante
2
entre pompe
entre pompe
vapeur saturante
P
V
NPSH h
2g g
= +


(Le plan de rfrence est pris au niveau de lentre de la pompe)
Cette caractristique est fournie par linstallateur dans le cas ou lon parle de NPSH
disponible. Cette expression est indpendante de la pompe utilise et ne dpend que de
linstallation (Longueur et diamtre daspiration, pertes de charge le long de cette
canalisation, hauteur gomtrique daspiration).
On parle de NPSH requis pour la charge nette laspiration fournie par le constructeur de
la pompe. Cette relation est indpendante de linstallation de la pompe et ne dpend que
de ce qui se passe dans la pompe.
Pour que linstallation puisse correctement fonctionner il faut : NPSH
disponible
> NPSH
requis
2-Equations de Bernoulli
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
47
2.2. - Ecriture de Bernoulli pour les fluides rels

Nous avons vu que pour le cas dun fluide rel et en rgime permanent, dautres forces
interviennent, notamment les forces dues au frottement, qui font apparatre une dissipation
de lnergie mcanique en nergie thermique.

Le long dune ligne de courant, lnergie mcanique diminue du fait de laugmentation de
lnergie thermique ; on appelle ce phnomne la perte de charge due aux frottements dans
un liquide.
Aprs intgration entre deux points sur la mme ligne de courant lquation de Bernoulli
scrit :
2 1 2
2
2
2
1
1
2
1
j h
p
g 2
V
h
p
g 2
V

+ +

+ = +

+

2 1
j

: est la perte de charge entre 1 et 2 en (m).


Toute la difficult consiste dterminer cette perte de charge.
48 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Une reprsentation graphique peut tre visualise par la figure suivante :




3. - EVALUATION DES PERTES DE CHARGE
Lcoulement dun fluide rel dans une conduite reprsente une des applications
classiques de lhydrodynamique thorique et exprimentale. La matrise de ces
coulements pour le calcul du fonctionnement des rseaux deau potable est indispensable.
On distingue deux types de perte de charge :
- la perte de charge linaire reprsentant lnergie perdue entre les deux points,
- la perte de charge singulire qui intervient lorsque lcoulement uniforme est
localement perturb.
3.1. - Les pertes de charge linaire
3.1.1. - Profil des vitesses dans une section circulaire en charge
Dans le cas dun coulement permanent 0
t
=

, le profil des vitesses dans une section


est dpendant de lintensit de la vitesse et de la gomtrie de louvrage. En effet, nous
avons vu que la turbulence apparat quand la vitesse dpasse un seuil (exprience de
Reynolds). Ceci peut tre reprsent par la figure suivante :
3- Evaluation des pertes de charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
49

Dans la zone laminaire seules les forces de frottement interviennent. Dans la zone
turbulente les forces de turbulence deviennent prpondrantes et une sous-couche
visqueuse trs mince apparat.

Dans le cas des conduites en charge o lcoulement est laminaire, seules les forces de
viscosit interviennent ; le profil des vitesses est parabolique et est donn par la figure
suivante :



Dans le cas o lcoulement est turbulent, le profil des vitesses tend tre uniformis et
est du type logarithmique.



On constate une zone centrale de pleine turbulence o le gradient de vitesse est trs faible
et le profil est aplati en son centre. La zone de la couche limite est de trs faible paisseur
et proche de la paroi. Le gradient de vitesse est trs important et les forces de viscosit
sont donc trs importantes dans cette couche. Une approche mathmatique exacte nest
50 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
pas possible dans ce type dcoulement complexe. Les relations qui seront construite sont
toutes dorigines empiriques.

Deux types dcoulements turbulents dans les conduites relles (prsentant des asprits)
existent :
- conduite hydraulique lisse : lcoulement est spar de la paroi par un film
laminaire,
- conduite hydraulique rugueuse : la sous-couche laminaire disparat et la
turbulence arrive jusqu la paroi.

Dcrivons le profil de lcoulement le long dun diamtre de conduite hydraulique lisse.
On distingue habituellement deux rgions : une rgion centrale dite externe qui correspond
plus de 90% de lcoulement et une rgion voisine de la paroi dite interne.

Dans la rgion centrale, trs aplatie, la vitesse est indpendante des asprits de la paroi et
dpend peu de la viscosit. Lcoulement est form de tourbillons importants, allongs,
contenant des tourbillons plus petits. Lintensit de la turbulence est peu trs constante
dans ce domaine.
Dans la zone interne, voisine de la paroi, la vitesse passe dune valeur lev zro sur une
distance relativement faible. Le fort gradient de vitesse qui en rsulte impose au fluide des
contraintes de cisaillement leves et donc dintenses frottements visqueux. On comprend
donc que la dissipation dnergie qui entrane une partie de la perte de charge ait lieu
principalement dans cette zone. En revanche, lintensit leve de la turbulence se traduit
par un accroissement important des transferts de quantit de mouvement entre particules
de fluide, dans toutes les directions. Ces particules tant porteuses de chaleur, il en rsulte
des proprits utilisables pour les transferts thermiques.

La turbulence est mesure par le nombre de Reynolds. Dans le cas des conduites
circulaires le domaine laminaire et turbulent est diffrenci par une valeur de Reynolds de
:
3- Evaluation des pertes de charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
51

La variation du dbit ou de la vitesse dans un tuyau de diamtre constant provoque une
perte de charge h
r
. Cette dernire varie en fonction de la vitesse peu prs linairement
quand lcoulement reste laminaire et quadratiquement quand lcoulement est turbulent.
Entre les deux types dcoulement il y a une zone de transition, cest--dire une zone
dincertitude o il faut faire preuve de beaucoup de rserve.
3.1.2. - Expression de la perte de charge linaire

Compte tenu des difficults pour rsoudre lquation de Navier-Stokes, la perte de charge
sera traduite par une quation empirique du type :
g 2
U
D
L
H
2
= ;
gD 2
U
L
H
J
2
=

=
est un coefficient de perte de charge. Il est sans dimension et est fonction du nombre de
Reynolds et de la rugosit de la paroi. (Parfois appel f)
L est la longueur de la conduite
D le diamtre
U la vitesse moyenne
J la perte de charge part unit de longueur
3.1.2.1. - En rgime laminaire Re < 2000
En rgime laminaire, seules les forces de viscosit interviennent. La vitesse est trs faible
et ltat de surface de la paroi nintervient pas dans le calcul du coefficient.
On dtermine ainsi :
Re
64
=
Pour une conduite circulaire, on a :
2
D
U
g
32
J

=
52 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
3.1.2.2. - En rgime turbulent Re > 4000
Une surface rugueuse implique un tat de surface telles que ses irrgularits ont une action
directe sur les forces de frottement, cest--dire quau del dun certain degr de rugosit
la configuration gomtrique de la surface exerce une influence bien dtermine sur
lcoulement. Dans le cas des coulements en charge appliqus au domaine de lingnieur,
la rugosit va jouer un rle trs important.
Une surface rugueuse peut tre considre comme tant constitue par une srie de
protubrances lmentaires. Elles sont caractrises par une certaine hauteur k (dsigne
parfois par )

Cette rugosit peut tre mesure par rapport au diamtre de la conduite. On parle alors de
rugosit relative :
D D
k
=
Si la rugosit est ondule le film laminaire pourra, dans une certaine mesure, se modeler
sur les ondulations. Si au contraire elle prsente des artes vives, le film sera aisment
corch puis dchir par les asprits, linfluence de la viscosit dans la couche limite sera
diminue et la turbulence fortement augmente. Ceci se traduit par une vitesse moyenne
rduite.

Les surfaces auxquelles on a affaire dans la pratique (bton, acier, fonte) ne prsentent pas
des caractristiques de rugosit uniforme. En effet, les protubrances de la surface sont
ingales et irrgulirement distribues. Elles proviennent du rivetage, des recouvrements,
des joints, des dfauts sur la paroi, de corrosion, dincrustations...
On dfinit ainsi deux types de rugosit : lisse et rugueuse. La diffrence entre les deux
dpend du nombre de Reynolds.

3- Evaluation des pertes de charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
53
Les tableaux suivants donnent quelques indications sur la valeur de .

54 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES


3- Evaluation des pertes de charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
55
Le calcul de la perte de charge linaire va donc prendre en compte le nombre de Reynolds
et la rugosit de la conduite. Le choix se fait en fonction de la turbulence cest--dire de la
valeur du nombre de Reynolds.


Rgime turbulent en tuyau lisse :
4000 < Re et
. Re.
12
D

: Dans cette zone, la turbulence est encore modre.
Lpaisseur de la sous-couche limite est suffisante pour englober toutes les asprits de la
conduite qui se comporte ds lors comme un tuyau lisse.

Deux expressions empiriques sont souvent utilises:
Formule de Blasius :
25 . 0
Re) 100 (

=
Formule implicite de Von Karman :
|
|

\
|

=
51 . 2
Re
log . 2
1
10



Rgime turbulent en tuyau rugueux :
200
D
. Re .


: La turbulence devient trs importante et ne dpend plus que de
D D
k
=
Formule de Nikuradse :
|
|

\
|
=
D 71 , 3
log . 2
1
10




Rgime turbulent en tuyau lisse, rugueux et en zone de transition :
Tous les rsultats de ces travaux sont rsums dans la formule de Colebrook qui sest
fortement inspir des rsultats de Von Karman et de Nikuradse.
Formule de Colebrook : |

\
|
+

=
D 7 . 3 Re
51 . 2
log . 2
1
10
pour 10
8
> Re > 4000
Cette relation implicite est difficile exploiter analytiquement et est le plus souvent
reprsente sur un graphique.
56 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES

ks/D
f

3- Evaluation des pertes de charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
57
On pourra utiliser les tables de Lechapt et Calmon en 1965 en remplaant la formule de
Colebrook par une formule approche de la forme : J = L.Q
M
.D
-N
o L, M et N sont des
constantes pour une rugosit donne. Pour des valeurs de vitesse comprises en 0.4m/s et
2m/s cette formule conduit des carts relatifs infrieurs 3% par rapport la formule de
Colebrook.

La formule dHazen et Williams est galement souvent utilise :
87 . 4 85 . 1
HW
85 . 1
D . C
Q
69 . 10 J = ; Coefficient de Hazen et Williams fonction de k.


3.2. - Les pertes de charge singulire

La perte de charge singulire, localise dans une section de la conduite, est provoque par
un changement de direction et dintensit de la vitesse (voir premier chapitre).
Lcoulement uniforme est perturb et devient localement un coulement non uniforme.
La turbulence joue un rle considrable, alors que les forces de viscosit sont
ngligeables. La perte de charge na donc lieu quen rgime turbulent.
Une telle non-uniformit de la vitesse peut tre provoque par :
- un branchement de section de la conduite,
- un changement de direction (coude),
- un branchement ou raccordement,
- un dispositif de mesure et contrle de dbit...

Comme pour les pertes de charge linaire, les pertes de charges singulires se traduisent
par la relation :
g 2
V
K H
2
=
K est fonction des caractristiques gomtriques et du nombre de Reynolds.


La valeur de K est donne pour les diffrents cas les plus classiques dans les tableaux
suivants :
58 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Raccordement dune conduite avec un grand rservoir
Dpart
Sans saillie lintrieur du rservoir, avec raccordement angles vifs


5 . 0 K =
Sans saillie lintrieur du rservoir, avec raccordement angles vifs, ajutage dbitant
gueule be


1 K =
Avec saillie lintrieur du rservoir


1 K =

Pour une saillie dont la longueur est comprise
entre 1 et 2 fois le diamtre
Sans saillie lintrieur du rservoir, avec raccordement de profil arrondi


05 . 0 K =

Cette valeur est une moyenne, elle dpend du
profil de larrondi.

3- Evaluation des pertes de charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
59

Arrive
S
1
S
2


1 K =

2
2
1
2
2
1
S
S
9
1
S
S
1 K
|
|

\
|
+
|
|

\
|
=
S
2
>>S
1
=> K=1

Coudes
Arrondi

|
|
|

\
|
|

\
|
+ =
90 r . 2
d
847 . 1 131 . 0 K
2
7

en degrs

r/d => 1 1.5 2 2.5 3
()
11.25 0.037 0.021 0.018 0.017 0.017
22.5 0.074 0.043 0.036 0.034 0.034
30 0.098 0.057 0.048 0.046 0.045
45 0.147 0.085 0.073 0.069 0.067
60 0.196 0.114 0.097 0.092 0.090
90 0.294 0.170 0.145 0.138 0.134
180 0.589 0.341 0.291 0.275 0.269

Brusque


K est indpendant du diamtre

() 22.5 30 45 60 90
K 0.07 0.11 0.24 0.47 1.13

60 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Ts
Branchement de prise 90 de mme diamtre et angles vifs


g 2
V
K H
2
t
r r
=
g 2
V
K H
2
t
b b
=

Q
b
/Q
t Q
b
/Q
t
0 0.2 0.4 0.6 0.8 Q
b
/Q
t
1
K
r
0.40 0.26 0.15 0.06 0.02 0.00
K
b
1.00 1.01 1.05 1.15 1.32 1.45

Branchement damene 90 de mme diamtre et angles vifs

g 2
V
K H
2
t
r r
=
g 2
V
K H
2
t
b b
=

Q
b
/Q
t Q
b
/Q
t
0 0.1 0.2 0.4 0.8 Q
b
/Q
t
1
K
r
0 0.16 0.27 0.46 0.60 0.55
K
b
-0.60 -0.37 -0.18 0.26 0.94 1.20


Cnes
Convergent


La perte de charge est ngligeable.

3- Evaluation des pertes de charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
61
Divergent


2
2
2
1
25 . 1
D
D
1 .
2
tg . 2 , 3 K
|
|

\
|
|
|

\
|

|
|

\
|
|

\
|
=
Si > 20 il y a dcollement et le comportement
est identique celui dun largissement brusque.

2
2
1
D
D
|
|

\
|

3 6 8 10 12 14 16 20 24 30 40
0 0.03 0.08 0.12 0.15 0.19 0.23 0.28 0.37 0.46 0.62 0.90
0.05 0.03 0.07 0.10 0.14 0.17 0.21 0.25 0.33 0.42 0.56 0.82
0.1 0.03 0.06 0.09 0.12 0.16 0.19 0.22 0.30 0.37 0.50 0.73
0.2 0.02 0.05 0.07 0.10 0.12 0.15 0.18 0.23 0.30 0.39 0.58
0.3 0.02 0.04 0.06 0.07 0.09 0.11 0.13 0.18 0.23 0.30 0.44
0.4 0.01 0.03 0.04 0.05 0.07 0.08 0.10 0.13 0.17 0.22 0.33
0.5 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.09 0.12 0.15 0.23
0.6 0.01 0.01 0.02 0.02 0.03 0.04 0.04 0.06 0.07 0.10 0.14


Changement brusque de diamtre
Rtrcissement

|
|

\
|
|
|

\
|
=
2
1
2
D
D
1 . 5 , 0 K

1
2
D
D

0.01 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
K 0.500 0.495 0.480 0.455 0.420 0.375 0.320 0.255 0.180 0.095

Elargissement

2
2
1
2
2
1
S
S
9
1
S
S
1 K
|
|

\
|
+
|
|

\
|
=


2
1
D
D

0.01 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
K 1.000 0.980 0.922 0.829 0.708 0.569 0.424 0.287 0.175 0.109

62 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Appareils de robinetterie
Vanne opercule



D
p

8
1

8
2

8
3

8
4

8
5

8
6

8
7

K 0.07 0.26 0.81 2.06 5.52 17 98







Vanne papillon


5 10 15 20 30 40 45 50 60 70
K 0.24 0.52 0.90 1.5 3.9 11 19 33 118 750


Robinets Boisseau


5 10 15 20 30 40 45 50 55 60
K 0.05 0.29 0.75 1.6 5.5 17 31 53 110 206


Clapet battant


20 30 40 45 50 55 60 65 70 75
K 1.7 3.2 6.6 9.5 14 20 30 42 62 90

3- Evaluation des pertes de charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
63

3.3. - Exemple de trac de perte de charge

Considrons une conduite de longueur L
12
. On dispose en 3 dun rtrcissement brusque.



Supposons deux conduites horizontales avec une liaison du type largissement ou
rtrcissement brusque.

64 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
4. - POSITION DE LA LIGNE PIEZOMETRIQUE
Considrons une conduite reliant deux rservoirs. La ligne pizomtrique correspondant
aux pression relatives est reprsente approximativement par la droite AA (On a nglig
les pertes de charges dues lentre et la sortie des rservoirs. La ligne pizomtrique
BB correspond aux pressions absolues (Pa/ = 10.33m).

Si la conduite toute entire est situe au dessous de AA, la pression dpasse la pression
atmosphrique. Cette hypothse correspond une situation normale. Il faut prvoir des
ventouses aux points les plus levs pour la sortie de lair accumul et des dcharges de
fond N et N pour la vidange et le nettoyage.


Si la conduite passe au-dessus de la ligne pizomtrique AA, la partie du tronon au-
dessus de AA est en dpression. En gnral, on doit viter les zones en dpression ; en
effet, une fente ventuelle facilitera lentre de corps trangers et peut contaminer leau.
La pose dune ventouse normale est contre indique, tant donn quelle permettrait
4- Position de la ligne pizomtrique
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
65
lentre de lair et provoquerait en consquence, une rduction de dbit. Dans la majorit
des cas, on vite ce type de trac. Si toutefois le trac


Si la conduite slve au-dessus de la ligne horizontale qui passe par A, il ny aura
coulement que si toute la conduite a t remplie deau au pralable (effets de
siphonnage).


Si la forme de la conduite est toujours situe au-dessous de lhorizontale passant par A,
mais dpassant BB, lcoulement se fait sans siphonnage. Cependant, la pression est nulle
pour la partie situe au-dessus de BB : il y a cavitation. Lcoulement se fait avec une
section partiellement pleine et le dbit diminue. Les dgagements dair et de vapeur deau
dans les rgions de forte dpression rendent lcoulement irrgulier.

66 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES


Si la conduite a un tronon au-dessus de lhorizontale passant par A, il faudrait amorcer le
siphon, lcoulement se ferait en section partiellement pleine et de faon trs irrgulire.

Si la conduite dpasse la cote B, il est impossible damorcer lcoulement.


5. - LES CONSTRUCTIONS GRAPHIQUES

Ces constructions permettent de dterminer le point de fonctionnement dun rseau (dbit,
pression) en utilisant les courbes de fonctionnement et les courbes caractristiques des
pompes.
Prenons lexemple suivant :

5- Les constructions graphiques
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
67
Dterminons dans un premier temps la caractristique du rseau en fonction du dbit :
CR(Q).
Sortie pompe Entre pompe
CR(Q) H H =
Cette caractristique doit tre dtermine en passant uniqiement par le rseau.
Dterminons, dans un premier temps, la charge totale (en mCE) au point B, H
B
, par
application du thorme de Bernouilli.


B 1 1 B
H h j

=


Dterminons, dans un deuxime temps, la charge totale (en mCE) au point C, H
C
, par
application du thorme de Bernouilli.
C 2 C 2
2 2
C C 2 2
C 2 C 2
H H PdC
p V p V
h h j
g 2g g 2g

= +
+ + = + + +


C 2 C 2
H h j

= +

On peut ensuite reprsenter la caractristique du rseau :
( ) ( )
C B 2 1 aspiration refoulement
CR(Q) H H h h PdC PdC = = + +

.
Il suffit de tracer la caractristique de la pompe H
MT
en fonction du dbit (donne fournie
par le constructeur) : H
MT
(Q).
Point de fonctionnement d'une instalation
0
10
20
30
40
50
60
70
80
0 10 20 30 40 50 60 70 80
Dbit (l/s)
H

(
m
C
E
)
Courbe
caractristique de la
pompe : H
MT
(Q)
Courbe
caractristique du
rseau : CR(Q)
Point de
fonctionnement de
l'installation
Hauteur
gomtrique de
l'installation

Nous pouvons ainsi dterminer le point de fonctionnement de la pompe et de linstallation.
1 B 1 B
2 2
1 1 B B
1 B 1 B
H H PdC
p V p V
h h j
g 2g g 2g

= +
+ + = + + +

68 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
6. - LE COUP DE BELIER
6.1. - Introduction
Limportance technique des coups de blier est extrmement grande. Lamplitude et la
soudainet des surpressions et dpressions qui sont mises en jeu peuvent dtriorer la
conduite ou les organes qui y sont branche. Il est ncessaire par consquent dempcher
la production de telles perturbations ou du moins de les attnuer.

En gnral, on peut considrer les fluides comme incompressibles. Toutefois les milieux
relles prsentent une certaines compressibilits dont il faut tenir compte sil sy produit
de rapide variation de pression (mme trs petites) : ce sont les phnomnes de
compressibilit qui permettent dexpliquer la propagation des ondes lastiques dans les
milieux matriels.

Quand un milieu matriel est affect localement par des variations de pression dorigine
quelconque, la perturbation qui en rsulte se propage dans toutes les directions. Cette
vitesse de propagation n'est pas infinie, mais si les variations de pression sont lentes on
peut admettre quelles se distribuent immdiatement dans tout le milieu (fonctionnement
dun rseau dadduction deau sous pression avec un dmarrage lent dune pompe). Si au
contraire les variations de pression sont rapides, la rpartition de la pression dans le milieu
ne peut pas tre considr comme instantane (dmarrage instantan dune pompe).

Les ondes lastiques se propagent sans modification dans un milieu isotrope infini, mais
sont susceptibles de rflexion et de rfraction quand elles rencontrent une surface sparant
deux milieux diffrents. Deux cas sont envisager :
Quand une onde plane se propageant dans un fluide, rencontre normalement une
surface rigide, elle se rflchit sans changer de signe : ainsi une onde de compression
se rflchit en une onde de compression.
Quand une onde plane se propageant dans un fluide, rencontre normalement une
surface ou la pression reste constante (surface libre dun liquide), il y a rflexion avec
changement de signe : ainsi une onde de compression donne naissance une onde de
dpression.

Dans un fluide parfait les seules ondes lastiques susceptibles de se propager sont des
ondes longitudinales, cest--dire des ondes provoquant des dplacements de particules
normalement aux surfaces dondes.
Vitesse des ondes dans lair : 331m/s.
Vitesse des ondes dans leau : 1 410m/s.

Ainsi, les ondes lastiques longitudinales peuvent se propager dans le fluide contenu dans
une conduite. En faisant abstraction des perturbations provoques par le frottement sur les
parois, on peut admettre que pour une conduite cylindrique, la propagation seffectue par
des ondes planes perpendiculaires laxe de la conduite. Mais la clrit des ondes est
diffrente de la valeur donne prcdemment. En effet, sous laction des variations de
pression provoques par une onde, les parois de la conduite de dforment. Ces
dformations sont petites en gnral, mais il faut en tenir compte dans le calcul de la
vitesse des ondes. Si elle est de lordre de 1000 m/s dans les conduites dacier elle peut
6- Le coup de blier
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
69
descendre quelques centaines de m/s pour les conduites en matire plastique et 15 m/s
pour des conduites en caoutchouc.
Ds que la vanne est ferme (t 0) la ligne pizomtrique remonte jusqu PP une cte
suprieure celle de la ligne de charge. Cette augmentation correspond lnergie
potentielle emmagasine par la compressibilit du fluide. Cette surpression remonte la
conduite et arrive au rservoir en laisant derrire son passage la conduite en surpression.


70 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Arrivant en A, au rservoir, londe subit une rflexion avec changement de signe : elle se
transforme en une onde de dpression. Cette dpression sajoute la surpression laisse
par londe montante prcdente. Puis londe de dpression arrive en B, rencontre la vanne
ferme et y subit une rflexion mais ici sans changement de signe : elle se transforme en
une onde ascendante de dpression.


De tout cela il rsulte quil est ncessaire de calculer les conduites pour quelles rsistent
ces phnomnes : surpressions et dpressions. En particulier, elles devront rsister
lcrasement d la pression atmosphrique dans le cas o les dpressions seraient
suffisantes pour crer la cavitation.
6- Le coup de blier
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
71

On distingue le coup de blier de masse (parfois appel oscillation en masse ) et le
coup de blier dondes (appel simplement coup de blier dans le langage courant).
Dans le premier cas, lcoulement rpond la modification du rgime dcoulement
uniquement par des transferts de masse ; la masse volumique du fluide est inchange et la
canalisation ne se dforme pas. Dans le second cas, il y rpond par des transferts de masse,
une modification de la masse volumique du fluide et ventuellement par une dformation
de la canalisation.
Lobjectif de cette partie est de faire comprendre ces phnomnes complexes au moyen
dune description mathmatique simplifie.
6.2. - Formulation mathmatique
Les deux quations de base de lhydraulique, savoir la loi de conservation de la masse
(quation 1) et la loi de la quantit de mouvement (quation 2), sont crites ci-dessous
pour un coulement unidimensionnel et permanent en moyenne dans le cas dun fluide
compressible et dune canalisation de section dformable (Perns 2003).
( )
( )
0
SU
S
t x

+ =

(1)
*
1 U U p
U gJ
t x x

+ =

(2)
Ici,
*
p p gz = + (3)
Dans ces quations, t est le temps ; x, labscisse curviligne suivant la canalisation (pas
ncessairement horizontale) ; , la masse volumique de leau ; S, la section de
lcoulement (section intrieure de la conduite) ; U , la vitesse permanente en moyenne
uniformise sur la section ; p, la pression relative par rapport la pression atmosphrique ;
g, lacclration gravitationnelle ; z, la cote altimtrique (axe vertical) ; J, le coefficient
adimensionnel de perte de charge.
Le terme de gauche de lquation 2 correspond la drive particulaire,
DU
Dt
, en suivant
un lment de fluide se dplaant la vitesse U .
Les variables x et t de chacune des grandeurs de ces quations ne sont pas indiques pour
des raisons de clart des quations. En toute rigueur, nous devrions crire ( ) ( ) , S x t
plutt que simplement S ,
( )
( ) , SU x t plutt que SU , etc.
6.3. - Coup de blier de masse
Dans le cas du coup de blier de masse, les modifications du rgime permanent de
lcoulement sont suffisamment faibles pour que la masse volumique du fluide ne soit pas
modifie et que la canalisation ne soit pas dforme. Dans ce cas, lquation de
conservation de la masse scrit comme suit :
( )
0
U
S
x

(4)
La prcdente quation signifie que la vitesse dans la canalisation, a priori fonction de x et
de t, est uniquement une fonction du temps t :
( ) ( ) , U x t U t = (5)
72 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Lquation de la quantit de mouvement scrit alors :
*
1
0
dU p
gJ
dt x

+ + =

(6)
Appliquons cette dernire quation au calcul des oscillations de niveau deau dans une
chemine dquilibre (protection anti-blier), comme illustr sur la Figure. . Cet exemple a
t propos par Frelin (Frelin 2002).

Figure. Installation avec chemine d'quilibre (Frelin 2002).
Ici, une conduite circulaire force entre une retenue deau de section importante et une
installation de turbinage est quipe dune chemine dquilibre de section circulaire.
Cette conduite est le sige dun coulement permanent jusqu linstant t = 0. A cet
instant, la vanne situe en amont de la turbine est ferme instantanment.
On considre les donnes suivantes: L = 4000 m, d = 1.5 m, D = 3 m, z
0
z
2
= 20 m,
Q
0
= 1 m
3
/s (dbit volumique).
On suppose que la masse volumique de leau reste constante et que la conduite est
compltement indformable. Dterminons, sous lhypothse dabsence de pertes de charge
(J = 0), la priode des oscillations du niveau deau dans la chemine ainsi que le niveau
maximal.
Pour cela, intgrons lquation (6) entre les abscisses curvilignes des points 1 et 2. Nous
obtenons lexpression suivante :
* *
2 1
1
0
dU
L p p
dt
( + =

(7)
Par ailleurs, compte tenu des sections importantes du rservoir et de la chemine par
rapport la section de la canalisation, lapplication de lquation de Bernoulli sans pertes
de charge conduit aux quations suivantes :
2 2
*
1
1 1
0 1
2 2
p p U U
z z
g g g g
= + + = + (8)
2 2
*
2
2 2
0 2
2 2
p p U U
z s z
g g g g
+ = + + = + (9)
La vitesse est en effet identique au point 1 et au point 2 dans la mesure o nous avons fait
les hypothses dun fluide incompressible et dune conduite indformable (quation 5).
Lutilisation de ces quations aboutit lexpression suivante pour la loi de la quantit de
mouvement :
6- Le coup de blier
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
73
0
dU
L gs
dt
+ = (10)
Le dbit dans la conduite tant identique celui dans la chemine, la loi de conservation
de la masse (quation 4) peut sexprimer, quelque soit le temps t, comme suit :
2 2
4 4
d D ds
U
dt

= (11)
La drive de la vitesse par rapport au temps scrit alors :
2
2
2
dU D d s
dt d dt
| |
=
|
\
(12)
En utilisant cette expression dans lquation (10), nous aboutissons lquation
diffrentielle suivante, o la seule variable est le temps t :
2
2
2
0
D d s
L gs
d dt
| |
+ =
|
\
(13)
Elle admet comme solution la forme gnrique suivante, o C
1
et C
2
sont des constantes :
1 2
cos sin
d g d g
s C t C t
D L D L
| | | |
= +
| |
| |
\ \
(14)
Sachant que s et
ds
dt
valent respectivement 0 et
0
2
4Q
D
linstant t = 0 (conditions
initiales), nous obtenons le niveau deau dans la chemine en fonction du temps.
0
4
sin
Q L d g
s t
Dd g D L
| |
=
|
|
\
(15)
Ainsi, la priode des oscillations vaut 2
D L
d g
, soit 4 minutes et 14 secondes. La
hauteur maximale dans la chemine vaut quant elle
0
0 2
4Q L
z z
Dd g
+ , soit 25.7 m.
6.4. - Coup de blier dondes
Considrons prsent la possible compressibilit du fluide et la possible dformation de la
conduite (uniquement de la section).
6.4.1. - Conservation de la masse pour un coulement dans lequel se propage une onde
Soit un coulement permanent dans une conduite. Suite la fermeture brutale dune
vanne, un coup de blier dondes se produit : une onde de clrit a prend naissance au
niveau de la vanne et se propage vers lamont. En assimilant londe une surface de
discontinuit, la conduite se divise en deux domaines distincts, comme illustr sur la
Figure. Le domaine 1 se situe entre lextrmit aval de la conduite (la vanne ferme situe
sur la gauche de la Figure) et le front de londe. Dans ce domaine, la vitesse de
lcoulement est nulle ; la masse volumique du fluide est + d ; la section de la conduite
est S + dS. Le domaine 2 est constitu de la partie complmentaire de la conduite. Dans la
mesure o le domaine 2 na pas encore t atteint par londe, les conditions suivantes y
rgnent : la vitesse de lcoulement vaut U ; la masse volumique de leau, ; la section de
la conduite, S.
74 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES

Figure. Schma de la propagation dune onde.
En se plaant dans un repre dobservation se dplaant la vitesse, constante, de londe
(Perns 2003), lquation de la conservation de la masse scrit :
( )( )
( )
d S dS a S U a + + = + (16)
En divisant par S et en ngligeant les infiniment petits dordre 2 par rapport aux
infiniment petits dordre 1, il vient lquation 17, qui sera utilise par la suite :
d dS U
S a

+ = (17)
6.4.2. - Clrit des ondes
En milieu infini (non confin), la clrit des ondes a pour expression :
a

= (18)
Dans cette relation, est le coefficient de compressibilit du fluide (environ 210
9
Pa pour
leau dans des conditions usuelles de temprature et de pression pour leau). Retenons
pour la clrit des ondes une valeur voisine de 1,400 m/s. La clrit crot avec la
temprature.
En conduite, la clrit est toujours infrieure cette valeur maximale. Dans le cas dune
canalisation de faible paisseur (e/D << 1), la clrit des ondes a pour expression la
relation suivante (Perns 2004) :
1
a
CD
eE

=
+
(19)
Dans cette quation, le terme rducteur dpend la fois des proprits du fluide et de celle
de la conduite. C est le coefficient dancrage ; D, le diamtre intrieur de la conduite ; e,
lpaisseur de la conduite ; E, le module dYoung du matriau de la conduite. Dautres
expressions applicables aux conduites paisses, aux galeries creuses dans une roche
brute, etc. sont explicites par Perns (Perns 2004).
6- Le coup de blier
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
75
Pour les conduites faonnes en matriau peu lastique, la clrit des ondes est en gnral
comprise entre 700 et 1,300 m/s ; elle peut descendre jusqu quelques dizaines de m/s
dans le cas de conduites lastiques, faonnes par exemple en caoutchouc ou en plastique
(Bonnin 1983). En considrant un ordre de grandeur de 1,000 m/s pour la clrit du coup
de blier dondes, on comprend la rapidit du phnomne dans des conduites peu
lastiques par rapport au coup de blier de masse.
6.4.3. - Equations aux caractristiques
Nous allons chercher exprimer partir de lquation 1 (conservation de la masse) une
expression mettant en jeu la pression toile et la vitesse. Utilisons pour cela les relations
suivantes, o
D
Dt
est la drive particulaire en suivant une particule se dplaant la
vitesse U .
( )
( )
( )
SU U
S
U S
x x x

= +

(20)
( ) ( ) ( ) D S S S
U
Dt t x

= +

(21)
En combinant ces expressions avec la loi de conservation de la masse, nous aboutissons
la relation suivante :
( )
0
D S
U
S
Dt x


+ =

(22)
En divisant cette expression par S, en utilisant lquation 17, puis en multipliant par a, il
vient une expression liant la pression toile et la vitesse :
*
1
0
Dp U
a
a Dt x

+ =

(23)
Ajouter et retrancher membre membre les quations 2 (loi de la quantit de mouvement)
et 23 (loi de conservation de la masse) aboutit aux expressions suivantes :
( ) ( )
* *
1 U U p p
U a U a gJ
t x a t x
( (
+ + + + + =
( (


(24)
( ) ( )
* *
1 U U p p
U a U a gJ
t x a t x
( (
+ + =
( (


(25)
Le premier crochet de lquation 24 peut sinterprter comme
DU
Dt
, drive particulaire de
la vitesse U en suivant une particule se dplaant la vitesse U a + ; le second crochet,
comme
*
Dp
Dt
, drive particulaire de la pression toile p
*
en suivant la mme particule.
De la mme faon, le premier crochet de lquation 25 peut sinterprter comme
DU
Dt
,
drive particulaire de la vitesse U en suivant une particule se dplaant la vitesse
U a ; le second crochet, comme
*
Dp
Dt
, drive particulaire de la pression toile p
*
en
suivant la mme particule. Les quations scrivent alors :
76 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
*
1
avec
DU Dp dx
gJ U a
Dt a Dt dt
+ = = + (26)
*
1
avec
DU Dp dx
gJ U a
Dt a Dt dt
= = (27)
La clrit des ondes tant gnralement trs grande (de lordre de 1,000 m/s) devant la
vitesse moyenne uniformise (de lordre de 1 m/s), nous pouvons approximer U a + par a
et U a par a. Considrons de plus quil ny a pas de pertes de charge (J = 0). Enfin,
utilisons comme variables le dbit Q transitant dans la canalisation et la hauteur
pizomtrique y dfinie comme la quantit
*
p
g
. Dans ces conditions, nous pouvons
crire :
0 avec
D aQ dx
y a
Dt gS dt
| |
+ = =
|
\
(28)
0 avec
D aQ dx
y a
Dt gS dt
| |
+ = =
|
\
(29)
La premire quation signifie quun observateur se dplaant la vitesse a verra la
quantit
aQ
y
gS
+ constante tout au long de son dplacement ; la seconde, quun
observateur se dplaant la vitesse a verra la quantit
aQ
y
gS
+ constante tout au long
de son dplacement.
La mthode de lpure Schnyder Bergeron consiste rsoudre graphiquement le
phnomne de coup de blier en tirant profit du rsultat prcdent. Plus prcisment, cette
mthode permet de dterminer la pression et le dbit en tout point de la conduite et tout
instant sans dterminer la nature exacte de londe mais simplement en considrant des
observateurs se dplaant aux vitesses a et a. Nous allons lutiliser pour quelques
exemples simples afin de mettre en vidence le phnomne de coup de blier dondes.
6- Le coup de blier
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
77
6.4.4. - Application la fermeture instantane dune vanne (sans pertes de charge dans la
conduite)
Considrons une conduite reliant un rservoir charge constante et une vanne, comme
illustr sur la Figure. Alors que lcoulement y tait permanent (dbit Q
0
) jusqu linstant
t = 0, la vanne est instantanment ferme linstant t = 0.

Figure. Fermeture dune vanne (Frelin 2002).
Pour des raisons de facilit de comprhension, on raisonnera partir du temps
adimensionnalis par la dure ncessaire londe pour parcourir la longueur de la
conduite :
t
L a
= (30)
Par exemple, londe gnre en V (vanne) = 0 arrivera en R (rservoir) = 1.
Lcoulement permanent consiste en un dbit Q
0
et une charge z
0
z
2
dans toute la
conduite (pas de pertes de charge). En supposant le terme
2
2
U
g
ngligeable dans
lexpression de la charge, nous pouvons affirmer que la hauteur pizomtrique y
0
en
rgime permanent vaut z
0
z
2
dans toute la conduite (pas de pertes de charge).
Une fois la vanne ferme, les conditions aux limites sont un dbit nul la vanne et une
hauteur pizomtrique gale y
0
au niveau de lentre du rservoir.

Figure. Epure pour une fermeture instantane de la vanne (Frelin 2002).
78 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
Plaons-nous en = 1
+
au niveau de la vanne (cest--dire une unit de temps
adimensionnel avant sa fermeture). Le signe
+
signifie quon se situe un instant
infiniment petit aprs le temps = 1 ; cette notation permet dviter toute ambigit due
la nature discontinue du coup de blier dondes. En se dplaant selon lquation
dx
a
dt
= , nous arriverons la vanne en = 0
+
, cest--dire juste aprs la fermeture de la
vanne. En = 1
+
, les conditions du rgime permanent rgnent : le dbit et la hauteur
pizomtrique valent respectivement Q
0
et y
0
au niveau du rservoir. En = 0
+
la vanne,
la quantit
aQ
y
gS
+ tant conserve, nous pouvons crire la relation suivante :
0
0
V0 R 1
aQ aQ aQ
y y y
gS gS gS
+ +

| | | |
+ = + = +
| |
\ \
(31)
Le dbit tant nul la vanne, nous pouvons dterminer la hauteur pizomtrique la
vanne juste aprs sa fermeture :
0
0
0 V
aQ
y y
gS
+
= + (32)
Cette expression montre que la fermeture instantane de la vanne est responsable dune
augmentation brutale de la hauteur pizomtrique de
0
aQ
gS
. Le point V0
+
est le point haut
de lpure (Figure). La droite reliant les points R1
+
et V0
+
a pour quation :
0
0
aQ aQ
y y
gS gS
+ = + (33)
Dplaons-nous prsent depuis la vanne selon lquation
dx
a
dt
= . Arrivant au rservoir
en = 1
+
, nous voyons conserve la quantit
aQ
y
gS
+ .
0
0
R1 V0
aQ aQ aQ
y y y
gS gS gS
+ +
| | | |
+ = + = +
| |
\ \
(34)
La hauteur pizomtrique lentre du rservoir valant y
0
, le dbit lentre du rservoir
vaut :
0
0 0 0
1 R
aQ gS
Q y y Q
a gS
+
(
= + =
(

(35)
Ainsi, du fait des fortes pressions rgnant dans la canalisation, le dbit sinverse : un
coulement prend place de la conduite vers le rservoir. La droite reliant les points V0
+
et
R1
+
a pour quation :
0
0
aQ aQ
y y
gS gS
+ = + (36)
Gagnons prsent la vanne en nous dplaant selon lquation
dx
a
dt
= . La quantit
aQ
y
gS
+ tant conserve, nous pouvons dterminer la pression la vanne en = 2
+
:
0
0
2 V
aQ
y y
gS
+
= (37)
6- Le coup de blier
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
79
Ceci signifie que la vanne subit une dpression. La droite reliant les points R1
+
et V2
+
a
pour quation :
0
0
aQ aQ
y y
gS gS
+ = (38)
En nous dplaant vers le rservoir selon lquation
dx
a
dt
= , nous pouvons exprimer le
dbit au niveau du rservoir en = 3
+
en conservant la quantit
aQ
y
gS
+ :
0
3 R
Q Q
+
= (39)
La droite reliant les points V2
+
et R3
+
a pour quation :
0
0
aQ aQ
y y
gS gS
+ = (40)
En poursuivant les allers-retours entre le rservoir et la vanne et en procdant de la mme
faon partir de diffrents points diffrents instants, nous pouvons dterminer le dbit et
la hauteur pizomtrique en tout point de la conduite et tout instant. Le phnomne
priodique ainsi mis en vidence est illustr sur la Figure. La Figure prsente quant elle
la variation de pression au niveau de la vanne ainsi que la variation de dbit au niveau de
lentre du rservoir. Nous constatons sur cette dernire figure que le coup de blier
dondes a ici une priode gale 4, cest--dire 16 s en considrant les donnes du calcul
effectu plus haut pour le coup de blier de masse dans une conduite quipe dune
chemine dquilibre et en supposant une clrit de 1,000 m/s. Nous pouvons ainsi
constater la frquence beaucoup plus importante du coup de blier dondes par rapport au
coup de blier de masse.
80 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES

Figure. Dcomposition du coup de blier dondes due la fermeture dune vanne
(Frelin 2002).

Figure. Variations de la pression la vanne et du dbit lentre du rservoir (Frelin
2002).
Des figures similaires la Figure pourraient tre traces en chaque point de la conduite.
La rsolution de lpure Schnyder Bergeron nous a permis de mettre en vidence les
pressions maximale et minimale subie par la conduite, respectivement
0
0
aQ
y
gS
+ et
0
0
aQ
y
gS
.
6- Le coup de blier
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
81
Il peut arriver que la pression absolue minimale atteigne la pression de vapeur saturante de
leau (0.239 mCE 20C). Dans ce cas, le liquide se vaporise ; une poche de cavitation se
forme. La conduite doit tre suffisamment rsistante pour ne pas imploser sous laction de
la pression atmosphrique sur sa surface extrieure. La dmarche dcrite ci-dessus, qui
suppose un fluide sous phase liquide et qui nest donc plus valable, doit tre adapte
(Perns 2004).
6.4.5. - Application la fermeture progressive dune vanne (sans pertes de charge dans la
conduite)
Considrons prsent une fermeture lente de la vanne, toujours sans considrer de pertes
de charge. Utilisons pour cela lquation 41 qui constitue un exemple possible de loi de
fermeture de vanne.
2
vanne 0
1
Q
t T

| |
=
|

\
(41)
Ici,
vanne
est la perte de charge singulire due la vanne (en mCE) ;
0
, la perte de charge
singulire de la vanne lorsquelle est compltement ouverte ; T, le temps total de
fermeture. Lpure illustre sur la Figure se construit de la mme faon que pour
lexemple prcdent, ceci prs quil faut considrer au moyen de lquation 41 un dbit
la vanne non nul pour tous les instants compris entre le dbut de la fermeture et la
fermeture complte.
Nous pouvons constater sur cette figure quune fermeture lente attnue lamplitude du
phnomne, aussi bien pour les surpressions que pour les dpressions.
82 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES

Figure. Epure pour une fermeture progressive de la vanne (Frelin 2002).
6.4.6. - Application la fermeture instantane dune vanne (prise en compte simplifie des
pertes de charge)
Considrons nouveau la fermeture instantane de la vanne mais cette fois-ci en
considrant dune faon simplifie les pertes de charge dans la conduite. Supposons ainsi
que la perte de charge dans la conduite soit localise lentre du rservoir (perte de
charge singulire entre les points R et W infiniment proches). Pour construire lpure,
nous considrerons des allers-retours non pas entre V et R mais entre V et W. La hauteur
pizomtrique en W correspond la hauteur pizomtrique en R laquelle on retranche la
perte de charge correspondant au dbit. Dans le cas des dbits ngatifs, remarquons que la
hauteur pizomtrique est plus importante en W quen R. Lpure ainsi trace, illustre sur
la Figure, converge vers son point central, cest--dire un dbit nul et la hauteur
pizomtrique rgnant dans le rservoir : les pertes de charge ayant attnu le coup de
blier, il ny a plus dcoulement.

6- Le coup de blier
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
83

Figure. Epure pour une fermeture instantane de la vanne avec une prise en compte
simplifie des pertes de charge (Frelin 2002).
6.5. - Protections anti-blier
On distingue les anti-bliers spcifiques des anti-bliers non spcifiques (Bonnin 1983).
Parmi les premiers, les vannes fermeture lente permettent dempcher les coups de blier
dus une fermeture brutale de vanne ; ces dispositifs sont en revanche inefficaces contre
les arrts brutaux de pompes, toujours susceptibles de se produire (panne du moteur,
coupure dalimentation, manuvre involontaire, etc.).
Toujours parmi les anti-bliers spcifiques, les volants dinertie sont des dispositifs
capables demmagasiner de lnergie pendant le dmarrage de la machine et de la fournir
la pompe ds la coupure de son alimentation : larbre de la pompe continue ainsi tourner
pendant quelques secondes quelques dizaines de secondes.
Les anti-bliers non spcifiques ne cherchent pas empcher les coups de blier mais
attnuer leurs effets. Dans la mesure o ils ne protgent pas entre le lieu de gnration du
coup de blier et leur localisation, ils doivent tre installs des endroits bien spcifiques.
Une soupape est ainsi un dispositif qui souvre la pression atmosphrique lorsque la
pression dans la conduite devient suprieure une valeur seuil. Ce faisant, la soupape
laisse chapper un certain dbit. Utile contre les surpressions, ce dispositif est en revanche
compltement inefficace contre les dpressions.
Pour attnuer aussi bien les surpressions que les dpressions, une mthode consiste
temprer les variations de pression dans la conduite par des variations de masse dans un
rservoir associ. Cest le principe de fonctionnement des chemines dquilibres
(rservoir lair libre), illustres sur la Figure, et des ballons anti-blier (rservoir sous
pression), illustrs quant eux sur la Figure. En fonctionnement permanent, le rservoir
est partiellement rempli. Lorsque la conduite est le sige dune surpression, le rservoir
accepte de leau de la conduite de faon temprer la surpression. Au contraire, lorsque la
conduite est le sige dune dpression, le rservoir apporte de leau la conduite de faon
temprer la dpression.
84 III- Hydraulique en charge
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES

Figure. Chemine dquilibre (Bonnin 1983).


Figure. Ballon anti-blier (Bonnin 1983).
6- Le coup de blier
Jos VAZQUEZ (Laboratoire Systmes Hydrauliques Urbains) ENGEES
85
6.6. - Rfrences bibliographiques sur les coups de blier
Bonnin, J., 1983. Ecoulement des fluides dans les tuyauteries. Techniques de l'ingnieur,
A 738, 1-22.

Frelin, M., 2002. Coups de blier. Techniques de l'ingnieur, BM 4 176, 1-27.

Perns, P., 2003. Hydraulique unidimensionnelle - Partie 1 - Analyse dimensionelle et
similitudes - Gnralits sur les coulements unidimensionnels - Ecoulements en
charge - Ecoulements surface libre, Cemagref Editions.

Perns, P., 2004. Hydraulique unidimensionnelle - Partie 2 - Coups de blier et
phnomne d'oscillation en masse - Pompes centrifuges, Cemagref Editions.

Perns, P., 2003. Mcanique des milieux continus dformables - Application la
mcanique des liquides parfaits et des liquides newtoniens, Cemagref Editions.

Vous aimerez peut-être aussi