HEC - Cours - MEPA 2014

2014
HYDROSTATIQUE ET
HYDRAULIQUE EN CHARGE
Formation Mastère Eau Potable et Assainissement
José VAZQUEZ et Matthieu DUSFRESNE

Dernière mise à jour : Septembre 2014
José VAZQUEZ et Matthieu DUFRESNE
Avant propos
L’hydraulique est incontournable dans le domaine de l’environnement. En effet, elle a une
place déterminante dans la compréhension, l’analyse et le diagnostic des réseaux
d’adduction d’eau potable, des stations de traitement, des réseaux d’assainissement et des
rivières. De plus, le contrôle de ces systèmes nécessite une instrumentation qui oblige le
concepteur et l’exploitant à une connaissance poussée du fonctionnement hydraulique de
ces ouvrages.
D’un point de vue réglementaire, la directive 2000/60/CE du Parlement européen établit

un cadre pour une politique communautaire dans le domaine de l’eau. Elle incite les Etats
membres (dont évidemment la France) à protéger et restaurer la qualité de leurs ressources
en eau afin de parvenir à un bon état chimique et écologique. L’eau est donc une
préoccupation majeure dans notre civilisation.
L’objectif de cet ouvrage destiné aux techniciens et ingénieurs est de fournir les bases
nécessaires à la compréhension et au calcul des phénomènes présents en hydraulique
appliquée au génie de l’eau et de l’environnement. Chaque notion d’hydraulique est
ponctuée par une série d’exercices permettant d’illustrer les concepts présentés. Les
exemples sont issus d’ouvrages hydrauliques existant en réseau. Les techniques de calcul
qui sont associées à la résolution des équations mises en œuvre sont élaborées dans un
souci d’efficacité.
Cet ouvrage est composé de plusieurs chapitres qui sont décrits dans les paragraphes
suivants.
Avant d’entrer directement dans la description mathématique et la quantification des

phénomènes hydrauliques rencontrés en eau potable, en assainissement et en rivière, le
premier chapitre est consacré à la description physique et phénoménologique des
différents types d’écoulements ainsi que des forces en présence. Elle permet aux non-
initiés de découvrir l’hydraulique de façon qualitative par le biais de l’expérimentation.
Le deuxième chapitre s’intéresse aux fluides au repos. Par exemple, la plupart des
dispositifs de contrôle de débit tels que les vannes autorégulées nécessitent un mécanisme
pouvant être approché par le biais de l’hydrostatique. Cette partie est donc consacrée à
l’action de l’eau sur les parois et sur les corps immergés.
Le chapitre suivant est dédié à l’hydraulique en charge. Ce type d’écoulement est le plus
souvent rencontré dans les réseaux d’adduction d’eau potable et parfois en assainissement.
Après un rappel des équations de Bernoulli, le paragraphe suivant s’intéresse aux pertes de
charge linéaires et singulières. On aborde ensuite le tracé de la ligne piézométrique et de la
ligne de charge qui permettent de caractériser le fonctionnement d’un réseau en terme de
sur ou sous pression. Un dernier sous-chapitre traite du coup de bélier avec les techniques
de calcul associées.
Bibliographie
AGHTM : Les stations de pompage d’eau, Editions Tec et Doc (2000).
BERTRAND-KRAJEWSKI J.L., Mesures en hydrologie urbaine et assainissement, éd.
Tec et doc, ed. 2000.
BONNIN, J., Ecoulement des fluides dans les tuyauteries. Techniques de l'ingénieur, A
738, 1-22, ed. 1983.
CARLIER M. : Hydraulique générale et appliquée, Editions Eyrolles (1972).
COMOLET R., Mécanique expérimentale des fluides, Masson, ed.1982.
IDEL’CIK I.E., Mémento des pertes de charge. Eyrolles (1986).
FRELIN, M., Coups de bélier. Techniques de l'ingénieur, BM 4 176, 1-27, (2002).
GRAF W. H., ALTINAKAR M. S. : Hydrodynamique : Une introduction, Traité de
Génie Civil, Ecole polytechnique fédérale 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 Génie Civil, Ecole polytechnique
fédérale de Lausanne, Presse polytechnique et universitaire romanes (1993).
GRAF W. H., ALTINAKAR M. S. : Hydraulique fluviale : écoulement non permanent
et phénomènes de transport, Tome 2, Traité de Génie Civil, Ecole polytechnique
fédérale de Lausanne, Presse polytechnique et universitaire romanes (1996).
HAGER W. H. : Wastewater hydraulics theory and practice, Springer, ed. 1999.
LENCASTRE A. : Hydraulique générale, Editions Eyrolles (1996).
LESIEUR M. : La turbulence, Presses Universitaires de Grenoble, Ed. 1994.
PERNES, P., Hydraulique unidimensionnelle - Partie 1 - Analyse dimensionelle et
similitudes - Généralités sur les écoulements unidimensionnels - Ecoulements en
charge - Ecoulements à surface libre, Cemagref Editions 2003.
PERNES, P., Hydraulique unidimensionnelle - Partie 2 - Coups de bélier et phénomène
d'oscillation en masse - Pompes centrifuges, Cemagref Editions 2004.
PERNES, P., Mécanique des milieux continus déformables - Application à la mécanique
des liquides parfaits et des liquides newtoniens, Cemagref Editions 2003.
SCHIESTEL R. : Modélisation et simulation des écoulements turbulents, Editions
Hermès (1993).
SINNIGER R.O., HAGER W. H. : Constructions hydrauliques : Ecoulements
stationnaires, Traité de Génie Civil, Ecole polytechnique fédérale de Lausanne,
Presse polytechnique et universitaire romanes (1989).
VIOLET P.L., CHABARD J.P., Mécanique des fluides appliquée, Presses des ponts et
chaussées, ed. 1998.

Sommaire
CHAPITRE I : CARACTERISTIQUES DES

ECOULEMENTS ......................................................................... 6
1. - ECOULEMENTS EN CHARGE ...................................................................................................... 6
1.1. - REGIMES D’ECOULEMENTS....................................................................................................... 6
1.2. - DISTRIBUTION DES VITESSES DANS UN TUYAU RECTILIGNE .................................................... 7
1.3. - ELARGISSEMENT ET RETRECISSEMENT DANS UNE CONDUITE .................................................. 8
1.4. - SORTIE D’UN RESERVOIR .......................................................................................................... 9
1.5. - COUDE ...................................................................................................................................... 9
1.6. - JONCTION ET BIFURCATION .................................................................................................... 10
1.7. - OBSTACLES DANS UN ECOULEMENT EN CHARGE ................................................................... 10
1.8. - LE COUP DE BELIER ................................................................................................................ 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 D’INFLUENCE AVAL .................................................................................................. 14
2.4. - COURBES DE REMOUS ............................................................................................................. 15
2.5. - LE RESSAUT ............................................................................................................................ 15
2.6. - LES SEUILS ET DEVERSOIRS .................................................................................................... 15
3. - PROPRIETES DES LIQUIDES...................................................................................................... 17
3.1. - MASSE VOLUMIQUE ................................................................................................................ 17
3.2. - POIDS SPECIFIQUE ................................................................................................................... 17
3.3. - VISCOSITE............................................................................................................................... 17
3.4. - PRESSION DE VAPEUR SATURANTE ......................................................................................... 19
3.5. - TENSION SUPERFICIELLE ........................................................................................................ 20
4. - CARACTERISATION DES FORCES DANS UN ECOULEMENT ..................................................... 21
4.1. - LES FORCES ............................................................................................................................ 21
4.2. - IMPORTANCE DES DIFFERENTES FORCES ................................................................................ 22
5. - LES EQUATIONS DE BASE ......................................................................................................... 23
6. - OUTILS MATHEMATIQUES ....................................................................................................... 24
CHAPITRE II : HYDROSTATIQUE ........................................ 25

1. - EQUATIONS DE L’HYDROSTATIQUE ........................................................................................ 25
2. - VARIATION DE LA PRESSION DANS UN FLUIDE INCOMPRESSIBLE ......................................... 26
3. - VARIATION DE LA PRESSION DANS UN FLUIDE COMPRESSIBLE ............................................ 28
4. - FORCES HYDROSTATIQUES SUR LES PAROIS.......................................................................... 28
4.1. - PAROI PLANE EN POSITION INCLINEE...................................................................................... 29
4.2. - PAROI A SURFACE GAUCHE..................................................................................................... 31
5. - FORCES HYDROSTATIQUES SUR DES CORPS IMMERGES ........................................................ 32
5.1. - FORCES D’ARCHIMEDE ........................................................................................................... 32
5.2. - EQUILIBRE DES CORPS IMMERGES .......................................................................................... 32
CHAPITRE III : L’HYDRAULIQUE EN CHARGE............... 34
1. - L’EQUATION DE CONTINUITE .................................................................................................. 34
2. - EQUATIONS DE BERNOULLI..................................................................................................... 36
2.1. - CAS DES FLUIDES PARFAITS.................................................................................................... 36
2.2. - ECRITURE DE BERNOULLI POUR LES FLUIDES REELS ............................................................. 49
3. - EVALUATION DES PERTES DE CHARGE ................................................................................... 51
3.1. - LES PERTES DE CHARGE LINEAIRE .......................................................................................... 51
3.2. - LES PERTES DE CHARGE SINGULIERE...................................................................................... 61
3.3. - EXEMPLE DE TRACE DE PERTE DE CHARGE ............................................................................ 67
4. - POSITION DE LA LIGNE PIEZOMETRIQUE ............................................................................... 68
5. - LES CONSTRUCTIONS GRAPHIQUES ........................................................................................ 70
6. - LES COUPS DE BELIER .............................................................................................................. 72
6.1. - INTRODUCTION ....................................................................................................................... 72
6.2. - FORMULATION MATHEMATIQUE ............................................................................................ 72
6.3. - COUP DE BELIER DE MASSE .................................................................................................... 73
6.4. - COUP DE BELIER D’ONDES ...................................................................................................... 75
6.5. - PROTECTIONS ANTI-BELIER .................................................................................................... 86

Chapitre I :
CARACTERISTIQUES DES ECOULEMENTS
L’hydraulique est l’étude des écoulements.

On distingue deux types d’écoulements :
 les écoulements en charge, dans lesquels l’eau remplit complètement la
canalisation, c’est le cas notamment des réseaux d’eau potable,
 les écoulements à surface libre (interface entre l’eau et l’air), c’est le cas des
rivières et des réseaux d’assainissement.
1. - ECOULEMENTS EN CHARGE
1.1. - Régimes d’écoulements
Le régime d’un écoulement se caractérise par la fluctuation temporelle des vitesses
et des pressions au sein de la veine liquide.
Expérience :
Un premier réservoir d’eau de niveau constant est vidangé par un tuyau. Une vanne
placée à l’extrémité du tuyau permet de faire varier le débit Q (m3/s). Un deuxième tuyau
est placé à l’intérieur du réservoir. Il contient un colorant et permet d’obtenir un mince
filet fluide coloré au centre du tuyau.
Mesure de la
Injection d’un colorant vitesse en un point
Sens de Vanne de
l’écoulement régulation du débit
Réservoir
Zone d’observation
Quand la vitesse est très faible (quelques millimètres par seconde) le filet coloré reste bien
défini, rectiligne et parallèle à l’axe du tuyau. Le régime est dit laminaire. L’écoulement
laminaire est rare dans le domaine de l’hydraulique de l’eau potable et de l’assainissement,
toutefois il n’est pas inexistant.
Les figurent suivantes montrent un écoulement laminaire dans une veine liquide et
l’évolution de la vitesse en un point (vitesse ponctuelle) en fonction du temps.
1. Ecoulements en charge 7
Evolution de la vitesse ponctuelle en fonction du

temps en régime laminaire
Sens de l’écoulement
0.0014
0.0012
Vitesse (m/s)
0.001
0.0008
0.0006
Injection Filet de colorant 0.0004
de colorant 0.0002
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Temps (s)
Quand la vitesse est plus élevée, le filet devient ondulé et très instable. Il se mélange
rapidement au fluide ambiant. Des tourbillons de différentes tailles apparaissent. Le
régime est dit turbulent.
Evolution de la vitesse ponctuelle en fonction du
temps en régime turbulent
0.8
0.7
0.6
Vitesse (m/s)
0.5
0.4
0.3
Injection Filet de colorant

0.2
de colorant 0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Temps (s)
La turbulence se caractérise donc par la création de tourbillons. ils mélangent les matières
dissoutes dans l’eau, comme par exemple le chlore dans un réseau d’eau potable ou le rejet
d’une station de traitement des eaux usées dans une rivière. La mise en place d’un
agitateur dans un bassin crée de la turbulence et ainsi il tend à homogénéiser les matières
dissoutes.
1.2. - Distribution des vitesses dans un tuyau rectiligne

En fonction du régime d’écoulement (laminaire ou turbulent) les forces dans
l’écoulement (elles seront détaillées dans un prochain chapitre) sont différentes. Cela a
pour effet une distribution de la vitesse ponctuelle moyennée dans le temps (que l’on
appelle vitesse moyenne temporelle) à l’intérieur du tuyau qui est différente suivant le
régime d’écoulement.
La connaissance du profil de vitesse dans une conduite permet de calculer le débit. La
plupart des capteurs en réseau mesurent la vitesse moyenne temporelle en quelques points
et reconstituent le profil complet de la vitesse pour en déduire le débit.
Le graphique suivant compare la forme du profil de la vitesse en régimes laminaire et
turbulent dans une conduite circulaire. En abscisse est représenté le rapport entre la
distance à la paroi (r) et le rayon de la conduite  R  . L’ordonnée représente le
 r 
pourcentage du rapport entre la vitesse maximale du profil et la vitesse ponctuelle V   .
R
On constate que la vitesse maximale est au centre de la canalisation. Par contre, le profil

8 I- Caractéristiques des écoulements
de vitesse turbulent varie beaucoup plus qu’en laminaire au voisinage de la paroi. Cette
zone de fort gradient de vitesse est appelée couche limite.
Evolution de la vitesse dans une conduite circulaire
Centre de la conduite
100
Profil de vitesse
90
Turbulent
Vitesse adimensionnelle V/Vmax en %
80
70
60
50
Profil de vitesse
Laminaire
40
30
20
Zone de fort
10 gradient de vitesse
0
0 0,2 0,4 0,6 0,8 1 0,8 0,6 0,4 0,2 0
Distance à la paroi r/R
1.3. - Elargissement et rétrécissement dans une conduite

La transition entre deux conduites de diamètre différent pour un écoulement en
charge provoque une répartition transversale des vitesses longitudinales totalement
différentes de celles vues précédemment.
On constate une zone ou les veines liquides proches de la paroi se décollent sur une
longueur L. Dans cette zone, on observe des recirculations à l’origine de perturbations
importantes dans l’écoulement.
Zone de décollement
de veine liquide
L
Rétrécissement
Elargissement
1.4. - Sortie d’un réservoir

La sortie d’un réservoir provoque également une modification de la répartition de
la vitesse et une zone de décollement 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 présence de cette zone a pour effet de diminuer légèrement la section de
passage de l’écoulement.
Zone de
recirculation
Sens de
Resserrement des
lignes de courants l’écoulement
1.5. - Coude
En raison de la courbure de la conduite, un mouvement hélicoïdal des lignes de
courant peut s’établir. La présence de ce mouvement de rotation persiste sur une longueur
en aval du coude (à peu près 50 fois le diamètre de la conduite). Lorsque la courbure est
importante, des zones de recirculation peuvent apparaître.
Le mécanisme interne de l’écoulement au travers d’un coude n’est pas encore bien connu,
malgré le grand nombre de chercheurs qui ont analysé cet élément important.
A-A
Zone de A
recirculation
Mouvement
hélicoïdal
Evolution des lignes de courant dans un coude circulaire

1.6. - Jonction et bifurcation

L’écoulement dans une jonction provoque une zone de séparation et une zone de
mélange. L’écoulement de la branche dans laquelle la vitesse est la plus faible est entraîné
par l’eau provenant de la branche ayant une vitesse plus élevée. Ce phénomène est à
l’origine de décélérations et d’accélérations de l’écoulement dans les deux branches.
Zone de mélange
Zone de
séparation
Jonction
L’écoulement dans une bifurcation se comporte de manière légèrement différente du fait

de la dérivation latérale. Une zone de séparation apparaît également dans la branche
latérale.
1.7. - Obstacles dans un écoulement en charge

Les figures suivantes montrent l’évolution des tourbillons à l’aval d’un obstacle
dans un écoulement en fonction de l’intensité de la vitesse.
Sens de
l’écoulement
Obstacle dans un écoulement. L’écoulement est à faible L’augmentation légère de la

vitesse. On observe l’apparition vitesse entraîne un allongement
de deux tourbillons symétriques des deux tourbillons.
attachés à l’obstacle. C’est le
début de la turbulence.
A vitesse plus importante les L’augmentation de la vitesse a

tourbillon commencent à pour effet d’accroître la
osciller, puis se détachent. fréquence des oscillations.
Les vannes à opercule ou les vannes papillons créent un obstacle dans les conduites. Ce
sont des organes de contrôle du débit.
Vanne à opercule Vanne papillon
1.8. - Le coup de bélier

Le coup de bélier est un ensemble de phénomènes hydrauliques complexes
provenant de la modification brutale du champ de vitesse et de pression dans une conduite.
L’importance technique des coups de bélier est extrêmement grande. En effet, l’amplitude
et la soudaineté des surpressions et dépressions qui sont mises en jeu peuvent détériorer la
conduite ou les organes qui y sont branchées. Il est nécessaire par conséquent d’empêcher
la production de telles perturbations ou du moins de les atténuer.
Le coup de bélier peut se produire, par exemple, dans le cas de la fermeture brutale d’une
vanne.
Vanne Considérons, dans un premier temps, un
Sens de ouverte
Réservoir l’écoulement réservoir se vidangeant par une vanne en
position ouverte. La vitesse dans la
Vitesse Uo canalisation est Uo.
Entrée Sortie
Vanne A un instant donné, la fermeture brutale de

fermée la vanne stoppe instantanément le fluide au
Réservoir Vitesse Uo
voisinage de celle-ci. Compte tenu de son
inertie, l’eau pousse fortement sur la vanne.
U=0 Une surpression importante se crée dans la
Zone de zone où la vitesse s’annule.
surpression
Onde de Gonflement Vanne A l’amont de la vanne, l’eau se comprime

choc de la conduite fermée
Réservoir Uo
faiblement et la canalisation gonfle. Tranche
U=0
par tranche les particules liquides
Déplacement transforment la vitesse en surpression. On
de l’onde observe ainsi la formation et la progression
d’une onde de choc. Celle-ci remonte de la
Zone de
surpression vanne où elle a pris naissance vers le
réservoir.
La vitesse de l’onde peut avoisiner les 1000
m/s en fonction du type de matériau de la
conduite.

Gonflement Vanne L’onde de choc atteint le réservoir. La

de la conduite fermée
Réservoir conduite entière est en surpression et la
U=0
vitesse de l’eau est nulle. Le réservoir est à
niveau d’eau constant. Sa pression ne
pouvant quasiment pas varier, le réservoir
Pression imposée Zone de
par le réservoir surpression
impose donc la pression à l’entrée de la
conduite.
Onde de Vanne La pression dans le réservoir étant beaucoup

choc fermée
Réservoir Uo
plus faible que la surpression dans la
U=0
conduite, la canalisation, au voisinage du
réservoir, va retrouver sa forme initiale. La
Déplacement
de l’onde diminution du diamètre de la conduite va
chasser l’eau vers le réservoir à la vitesse
Uo. De proche en proche se phénomène se
produit et crée ainsi une onde de choc se
déplaçant cette fois-ci vers la vanne.
U=0 Dès que l’onde atteint la vanne, la conduite

Réservoir Vitesse Uo entière retrouve quasiment sa pression
d’origine et sa vitesse initiale inversée. La
vanne étant toujours fermée, la vitesse de
Zone de l’eau à son voisinage est nulle.
dépression
Onde de
choc
Vanne Compte tenu de son inertie, l’eau au
fermée voisinage de la vanne entre en dépression et
Réservoir Uo
U=0 la conduite diminue son diamètre. Une onde
de choc en dépression se crée et se déplace
Déplacement
de l’onde Zone de vers le réservoir.
dépression
Gonflement Vanne Comme précédemment, le réservoir impose

de la conduite fermée
Réservoir la pression au niveau de l’entrée de la
U=0
conduite.
Pression imposée
par le réservoir
Onde de Vanne La conduite retrouve alors sa forme initiale.

choc fermée
Réservoir Uo
Le volume d’eau manquant est pris dans le
U=0
réservoir et la conduite retrouve sa vitesse
initiale Uo.
Déplacement
de l’onde
U=0 L’onde arrive au niveau de la vanne et on

Réservoir Vitesse Uo retrouve la configuration initiale.
2. Ecoulements à surface libre 13
2. - ECOULEMENTS A SURFACE LIBRE

2.1. - La surface libre
La surface libre est l’interface entre l’air et l’eau. La pression y est égale le plus
souvent à la pression atmosphérique.
Les écoulements dans les canaux naturels (rivière) 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 répartition des iso-vitesses dans le cas d’un régime
turbulent. Dans certains cas, la vitesse peut-être maximale un peu au-dessous de la surface
libre.
V
Vmax
1,0
0,75 1,0
0,5
0,25 Vitesse 0,95
maximale 0,9
0,85
Comme pour les écoulements en charge, nous avons également un régime d’écoulement
laminaire et turbulent. (L’écoulement laminaire est rare en hydraulique).

2.3. - Notion d’influence aval

Expérience : l’onde de gravité
Considérons un canal à pente nulle, de section et de hauteur d’eau constante. Le
fluide est au repos (vitesse nulle U=0). A un instant t, on perturbe la surface libre du canal.
Il se développe deux ondes se déplaçant à la même célérité c. Ces ondes sont appelées
ondes de gravité.
c c
U=0 U0 U=0 U0 U=0
Expérience : notion d’influence aval

On utilise cette fois-ci un canal à section transversale, pente, hauteur (h) et débit
(Q) constants. On crée une perturbation grâce à une vanne que l’on ferme et que l’on
ouvre très rapidement.
Au niveau de la surface libre, il se crée

Aval
deux ondes :
Amont  une onde se propageant vers l’aval,
c’ < 0  une onde se propageant vers
c’’ > 0
l’amont quand la vitesse (U) dans le
Q
canal est inférieure à la vitesse de
U<c l’onde de gravité (c).
Quand la vitesse (U) dans le canal est

Aval
supérieure à la vitesse de l’onde de gravité
Amont (c), alors deux ondes se propagent vers
c’ > 0 l’aval avec des vitesses différentes.
c’’ > 0
Q
U : vitesse de l’écoulement
U>c c : célérité des ondes
c’ : vitesse de l’onde amont
c’’ : vitesse de l’onde aval
De cette expérience, on en déduit que :

 dans le cas où la vitesse du fluide est supérieure à la vitesse de l’onde, l’amont
n’est pas influencé par les conditions hydrauliques à l’aval (régime torrentiel);
 dans le cas contraire, on a une remontée de l’onde qui va perturber l’amont (régime
fluvial), ce phénomène est appelé influence aval.
2. Ecoulements à surface libre 15
2.4. - Courbes de remous

Par rapport à l’écoulement en charge, un écoulement à surface libre a une difficulté
supplémentaire qui est la détermination de la position de la surface libre par rapport au
fond du canal (tirant d’eau). Celle-ci est variable en fonction des caractéristiques du fluide
et de l’écoulement. La figure suivante montre un exemple de forme de surface libre.
Sens de
l’écoulement
Seuil
2.5. - Le ressaut
Le ressaut hydraulique se caractérise par une variation importante et croissante de
la hauteur d’eau de l’amont vers l’aval du phénomène sur une courte distance. Dans la
plupart des cas, une agitation importante de la surface libre permet rapidement de localiser
le phénomène, comme par exemple dans le cas d’une ressaut fort.
Le ressaut hydraulique est l’un des phénomènes les plus complexes en hydraulique à
surface libre. Les connaissances actuelles sur le ressaut hydraulique ne sont pas encore
suffisamment étendues pour que l’écoulement interne soit parfaitement compris.
Ressaut fort Ressaut ondulé
2.6. - Les seuils et déversoirs

Le seuil crée un obstacle dans un canal, qui oblige le tirant d’eau à augmenter et
donc l’eau à passer par dessus.
Dans le cas d’un seuil dénoyé, l’eau chute à l’aval du seuil.
Dans le cas ou le tirant d’eau à l’aval de l’ouvrage est important, la chute d’eau ne peut
plus avoir lieu. Dans ces circonstances, le seuil est dit noyé.
Seuil dénoyé Seuil très faiblement noyé

Seuil faiblement noyé Seuil complètement noyé
En fonction de la forme du seuil et de la vitesse de l’écoulement, il peut apparaître une

zone de dépression à l’aval du seuil.
Zone de dépression
Le déversoir est un ouvrage de bifurcation qui permet un partage des débits dans
deux canaux. Par rapport à une simple bifurcation, où les débits sont partagés quelle que
soit la hauteur d’eau, dans un déversoir, le déversement n’a lieu que si la hauteur du fluide
atteint la hauteur de la crête déversante.
Déversoir latéral
Conduite Conduite
Amont Aval
Vue de dessus
Conduite
Déversée
AM ONT AVAL
M ilieu n a tu re l
R essa u t
h y d r a u liq u e
3. Propriétés des liquides 17
3. - PROPRIETES DES LIQUIDES
Dans l’établissement des principes de l’hydraulique, certaines propriétés des fluides jouent
un rôle important, d’autres seulement un rôle mineur ou aucun rôle du tout. En
hydrostatique (fluide au repos) c’est le poids spécifique qui est la propriété la plus
importante, tandis qu’en hydrodynamique (fluide en mouvement), la densité et la viscosité
sont des propriétés dominantes. La pression de vapeur prend de l’importance quand
interviennent des basses pressions, le liquide en question contient des bulles de vapeur,
c’est le phénomène de cavitation. La tension de surface influe sur les conditions statiques
et dynamiques dans les conduits très étroits, c’est le phénomène de capillarité.
3.1. - Masse volumique

Masse M
La masse volumique (ρ) est le rapport :    [Kg/m3]
Volume Vol
Pour les liquides, le volume est pratiquement insensible aux variations de pression et, dans
la majorité des cas, il augmente faiblement quand la température augmente, l’eau faisant
exception à cette règle en dessous de 4°C.
eau = 1000 kg/m3
mercure = 13546 kg/m3
air sec = 1,205 kg/m3
T (C°) 4 6 10 20 40 60 80 100
eau (kg/m3) 1000,03 1000,06 1000,30 1001,80 1007,85 1017,08 1029,02 1043,16
Attention : Contrairement aux liquides, les gaz sont fortement compressibles. La variation
de masse volumique dépend de la température et de la pression :   f  p, T  .
3.2. - Poids spécifique

Il représente la force de gravité agissant sur la masse par unité de volume :
 = g. [N/m3]
4 3
eau = 10 N/m
3.3. - Viscosité
La viscosité d’un fluide en mouvement est la propriété qui exprime sa résistance à une
force tangentielle.
Expérience :
y
L
Répartition de la vitesse
entre deux plaques en
régime laminaire

Considérons un fluide placé entre deux plaques planes, parallèles, distantes de L et

horizontales. L’une est fixe et l’autre est en mouvement uniforme de vitesse Uo. Pour
générer une vitesse de la plaque supérieure (surface A), il faut exercer une force F.
Cette force est la résultante des forces de frottements visqueux.

L’expérimentation permet de déduire 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 :
F U U
 0 
A L y
2
 [N.s/m ] est appelé viscosité dynamique ou absolue.
Le viscosimètre :
On considère deux cylindres coaxiaux séparés par un intervalle e dont l’espace entre eux
est rempli par un liquide. On fait tourner le cylindre extérieur à vitesse constante () et on
maintient fixe le cylindre intérieur.
Vitesse de
Cylindre fixe
rotation : 
Variation de la
vitesse de du
liquide suivant
l’épaisseur e
Niveau r
de l’eau
h
e
Cylindre en
rotation
Le fluide en contact avec le cylindre extérieur va y adhérer et par conséquent va être

animé de la vitesse V du cylindre extérieur. Le fluide en contact avec le cylindre fixe aura
une vitesse nulle. La viscosité fait naître une force de frottement que l’on mesure par le
couple M.
Les expériences ont montré que :
 si e est faible par rapport au rayon intérieur r, la courbe représentative 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 définit un deuxième coefficient de viscosité, le coefficient de viscosité cinématique :

 [m2/s]

eau à 20°C = 10-3 N.s/m2 mercure = 1,554.10-3 N.s/m2 air = 18,5.10-6 N.s/m2
eau à 20°C = 10-6 m2/s mercure = 0,1147.10-6 m2/s air = 15,6.10-6 m2/s
3. Propriétés des liquides 19
Température (°C) 0 5 10 20 30 40 50 60 70 80 90 100
Viscosité dynamique
(N.s/m2) 10‐3 1,787 1,519 1,307 1,002 0,798 0,653 0,547 0,467 0,404 0,355 0,315 0,282
Viscosité cinématique
(m2/s) 10‐6 1,787 1,519 1,307 1,004 0,801 0,658 0,553 0,475 0,413 0,365 0,326 0,294
On appelle fluide parfait un fluide dont la viscosité serait nulle (fluide inexistant dans la
nature). La viscosité existe dès qu’il y a mouvement relatif entre particules, que ce soit en
régime laminaire ou turbulent.
3.4. - Pression de vapeur saturante

L’ébullition est un phénomène de changement d’état, dans lequel le liquide passe à
l’état de vapeur. Tous les liquides ont tendance à s’évaporer ; la phase liquide se
transforme en phase gazeuse. Au cours de cette transformation, les molécules de vapeur
exercent une pression appelée pression de vapeur saturante. Dans le cas de l’eau, la
pression de vapeur (ps) croît avec une augmentation de la température (T).
La pression de vapeur saturante pour l’eau est donnée par la relation empirique suivante :
2795
log10  ps   22.435   3.868log10  T  273.15  avec ps en Pa et T en °Celsius.
T  273.15
T (C°) Ps (Pa) Ps(mCE)
1 645.2 0.07
5 855.1 0.09
10 1200.9 0.12
15 1664.9 0.17
20 2280.0 0.23
25 3086.1 0.31
30 4131.3 0.42
35 5472.8 0.56
40 7177.9 0.73
45 9325.1 0.95
50 12005.8 1.22
55 15324.5 1.56
60 19400.6 1.98
65 24369.3 2.48
70 30382.4 3.10
75 37609.6 3.83
80 46239.2 4.71
85 56478.6 5.76
90 68555.6 6.99
95 82718.4 8.43
100 99236.0 10.12
Si, à température constante, on abaisse la pression à la surface d’un liquide, ce dernier se

met à bouillir lorsqu’on atteint la pression de vapeur saturante correspondant à cette

température. Dans l’écoulement des liquides, il peut arriver que la pression en certains
points devienne inférieure à la pression de vapeur saturante. Le liquide entre alors
localement en ébullition et des bulles de vapeur apparaissent au sein même de
l’écoulement. Ce phénomène, appelé cavitation, est le plus souvent nuisible pour les
installations où il se produit (canalisation, pompes, turbine…). Les variations de volume
lors du changement d’état sont telles qu’il se produit au sein du fluide de véritables
explosions de bulles au moment de la vaporisation et de violentes implosions, lors de la
condensation.
3.5. - Tension superficielle
Une molécule liquide au repos est soumise aux forces d’attractions que les molécules
voisines exercent sur elle. Une molécule à la surface libre d’un liquide ou à la surface de
séparation de deux liquides non miscibles n’est plus soumise à l’action de forces
symétriques, puisqu’elle n’est plus entourée symétriquement par d’autres molécules de
même nature. Ainsi la résultante des forces moléculaires n’est plus nulle. La surface de
séparation se comporte comme une membrane tendue.
La force d’attraction tangentielle à la surface nécessaire pour arracher des particules
agissant le long d’un segment de longueur unitaire est appelée tension superficielle.
Surface libre
Molécules
Forces entre molécules

au niveau de la surface libre
Forces entre molécules

au sein du fluide
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 à 20°C : 0,0724 N/m.
4. caractérisation des forces dans un écoulement 21
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 d’accélération pure :
Elles proviennent de la variation de la vitesse (V) de la masse d’une fluide
(M) dans le temps.
V
Faccélération pure  M
t
Prenons par exemple deux réservoirs à la même hauteur, dont l’un est vide et
l’autre plein, reliés par une conduite de diamètre constant, horizontal et muni
d’une vanne. A l’ouverture de la vanne, il se produit un écoulement. La
variation de la vitesse dans le temps  V t  crée au sein de l’écoulement
une force d’accélération pure.
- Les forces d’accélération convective :
Elles proviennent de la variation de la vitesse (Vx, Vy, Vz) dans l’espace
(repère [x, y, z]).
 V V V 
Faccélération convective  M  .Vx  .Vy  .Vz 
 x y z 
Prenons une conduite dont l’écoulement ne varie pas dans le temps (
V / t  0 ). L’écoulement étant permanent, le débit est identique en tout
point de la canalisation. Or, si la surface A est supérieure à la surface B alors
la vitesse en B (Vb) est supérieur à la vitesse en A (Va).
Cette variation de vitesse va engendrer une accélération qui va générer une
force d’accélération convective.
A B C D
Régime Régime Régime
accéléré uniforme accéléré
Va Vb Vc Vd
Va  Vb Vc  Vd
0 0
Lab Lcd
Vb  Vc
0
L bc
Lab Lbc Lcd
Convergent Divergent

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
F
à la surface (A) d’un fluide : p 
A
- Les forces de frottement de viscosité :
Nous avons vu précédemment qu’un fluide, dont les particules sont en
mouvement relatif, génèrent des forces de frottement dues à la viscosité. Le
U
force de frottement s’écrit : F  A
y
- Les forces générées par la turbulence :
La turbulence décrite au premier chapitre joue un rôle majeur dans
l’écoulement des fluides. La turbulence à tendance à « freiner »
l’écoulement. Une façon de les représenter mathématiquement consiste à les
assimiler à des forces de frottement, ce qui est faux compte tenu de la nature
même de la turbulence.
4.2. - Importance des différentes forces

En général, l’hydraulicien doit s’occuper de l’effet de la force dominante. Dans la
plupart des problèmes d’écoulement des fluides, la pesanteur, la viscosité et l’élasticité
sont prépondérantes, mais pas toujours simultanément.
L’importance relative des différentes forces agissant sur un liquide est calculée par des
nombres adimensionnels représentant les rapports entre ces forces. L’analyse
dimensionnelle permet de simplifier ces rapports. Les différentes dimensions utilisées
sont :
L : longueur,
T : temps,
ρ : masse volumique.
4.2.1. - Le Nombre d’Euler

Le nombre de Euler est le rapport entre les forces d’inertie (Ma) et les forces de pression
(pA).
Ma L3 (L / T 2 ) V2 M : la masse,
 2
  a : l’accélération du fluide,
pA pL p
p : la pression,
A : la surface d’application de la pression,
V : vitesse caractéristique de l’écoulement.
Ce rapport est utilisé lorsqu’on calcul les forces (statiques et dynamiques) de l’eau sur une
paroi ou un ouvrage.
4.2.2. - Le Nombre de Reynolds

Le Nombre de Reynolds (Re) est le rapport entre les forces d’inertie (Ma) et les forces de
viscosité (τA).
Re= 
 
Ma L3 (L / T 2 )

VL τ : la force de frottement par unité de surface,
A  A : la surface d’application de la force de frottement,
 L/T 2
 L L : longueur caractéristique de l’écoulement.
 L 
4. caractérisation des forces dans un écoulement 23
Ce rapport permet de distinguer le régime laminaire et turbulent.

Si Re > valeur dépendant de l’ouvrage étudié alors turbulent
Si Re < valeur dépendant de l’ouvrage étudié alors laminaire
Exemples :
Ecoulement dans une Ecoulement entre deux Ecoulement autour d’une
conduite circulaire : plaques parallèles : sphère :
Si Re > 2500 alors turbulent Si Re < 500 alors laminaire Si Re < 1 alors laminaire
Si Re < 2000 alors laminaire Sinon turbulent Sinon turbulent
L : diamètre de la conduite. L : distance entre les deux L : diamètre de la sphère.
plaques.
4.2.3. - Le Nombre de Froude

Le nombre de Froude est le rapport entre les forces d’inertie (Ma) et les forces de
pesanteur (Mg).

 
Ma L3 (L / T 2 ) L2 V 2 V 2 g : l’accélération de la pesanteur,
 
Mg  
L3 g  
L3 g Lg L : longueur caractéristique de l’écoulement,
V : vitesse caractéristique de l’écoulement.
Ce rapport permet de distinguer le régime 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 d’inertie (Ma) et les forces d’élasticité
(EA).
Ma

 
L3 ( L / T 2 ) V 2

E : module d’élasticité.
EA EL2 E
4.2.5. - Le Nombre de Weber

Le nombre de Weber est le rapport entre les forces d’inertie (Ma) et les forces de tension
superficielle (σL).
Ma

 
L3 ( L / T 2 ) LV 2

σ : tension de surface.
L L 
5. - LES EQUATIONS DE BASE

La détermination des caractéristiques d’un écoulement consiste à rechercher la
pression et la vitesse en tous points. Pour cela, il faut écrire des équations d’équilibre (ou
de conservation) entre les différentes forces agissant sur le fluide.
Ces équations d’équilibre sont :
- l’équation de continuité ou équation de conservation de la masse (m) du fluide.
dm
0
dt
- l’équation de la quantité de mouvement qui traduit l’équilibre entre la somme
des forces extérieures (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)
 F
dt

6. - OUTILS MATHEMATIQUES
Elément de volume : dv = dx.dy.dz

dz
dy
dx

Dérivée partielle :
x
P P P P
Dérivée totale : dP  dt  dx  dy  dz
t x y z
dP P P x P y P z
   
dt t x t y t z t
 f 
 x 
  f 
Gradient d’un scalaire : grad ( f )  grad ( f )   
 y 
 f 
 
 z 
 Vx Vx Vx 
 z 
Vx   x y
Vy Vy Vy 
Gradient d’un vecteur : grad(V)  grad  Vy   
 x y z 
 Vz   Vz Vz Vz 
 
 x y z 
Vx Vy Vz
Divergence d’un vecteur : div ( V )   
x y z
  Vz Vy 
 x   y  z 
   Vx   Vx Vz 
 
Rotationnel : RotV       Vy    
 y   z x 
    Vz    Vy  Vx 
   x  y 
 z   
Chapitre II :
HYDROSTATIQUE
L’hydrostatique étudie les conditions d’équilibre des liquides au repos. Ce chapitre

aborde l’étude de la répartition de la pression, notamment en fonction de la distance
verticale, ainsi que les forces qui en résultent.
1. - EQUATIONS DE L’HYDROSTATIQUE
Nous avons vu au chapitre précédent que l’ensemble des forces agissant sur un
fluide sont de deux natures : les forces de volume et les forces de surface. En reprenant les
résultats acquis dans ce chapitre, nous allons détaillé l’ensemble de ses forces.
Considérons dans un réservoir un fluide au repos, dont on extrait un petit parallélépipède
d’eau d’axe vertical z. Soit p la pression en son centre. Il est soumis aux forces verticales
suivantes :
 p dz 
p  dx.dy
 z 2 
.g.dx.dy.dz
p(x,y,z)
dz
z
y  p dz 
p  dx.dy
Réservoir  z 2 
x
Concernant les forces de volume, il n’en existe qu’une seule la force de pesanteur. Elle
s’écrit de la façon suivante :
Fpesanteur  g.  dx.dy.dz 
Les forces d’inertie n’existe pas puisque le fluide est au repos (vitesse nulle).
Concernant les forces de surface, la pression agit sur la face supérieure et inférieure de
l’élément. Ces forces s’écrivent de la façon suivante :
 p dz 
Force de pression sur le surface inférieure : Fpression inf.   p   dx.dy
 z 2 
 p dz 
Force de pression sur le surface supérieure : Fpression sup.    p   dx.dy
 z 2 

26 II- Hydrostatique
Les forces de viscosité et de turbulence n’existent pas puisqu’il n’y a pas de vitesse
relative entre les particules de fluide.
L’équation de l’hydrostatique est déterminé en écrivant l’équilibre de l’ensemble des

forces :  F  Forces d' inertie . 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 
p
soit : -.g  0
z
On peut écrire de façon analogue les équations d’équilibre dans les autres directions :
p p
0; 0
x y
Ces trois équations montrent que la pression est indépendante de x et de y, c’est-à-dire que
la pression ne varie pas dans les directions x et y ou encore qu’elle est constante dans un
plan horizontale. Cela est vérifié tant que l’on reste dans un même fluide (ρ constante). La
pression ne dépend que de z, ce qu’on écrit :
p
 -.g ou dp=-.g.dz
z
2. - VARIATION DE LA PRESSION DANS UN FLUIDE

INCOMPRESSIBLE
Pour un fluide incompressible ( = Cte), l’intégration de l’équation ci-dessus
s’écrit : p + .g.z = Cte.
On écrit fréquemment : p* = p + .g.z = Cte.

On nomme p* l’énergie potentielle par unité de volume.
p* p
Une autre écriture consiste à diviser l’équation par .g :   z  Cste ,
g g
p*
On remarque que est homogène à z c’est-à-dire à une longueur ce qui offre un moyen
g
p*
pratique pour la représenter graphiquement. est appelée hauteur piézométrique.
g
Dans ce cas, la pression p est mesurée en hauteur de colonne du liquide  .
Unités de mesure :
L’unité légale est le pascal : 1 Pa = 1N/m2
Un multiple du pascal est le bar : 1 bar = 105 Pa
Il existe d’autres unités plus pratiques :
- Le mètre de colonne d’eau : 1,0mCE  1,0mCE.eau .g  9,81.103 Pa  0,098bar
105
1,0bar  105 Pa   10,19mCE
eau .g
2. Variation de la pression dans un fluide incompressible 27
- Le millimètre de mercure : 1,0mmHg  103  mercure .g  1,33.102 Pa

105
1.0bar  105 Pa  .1000  752,5mmHg
mercure .g
En hydraulique des réseaux urbains, la pression est le plus souvent exprimée en mètres de
colonne d’eau et en bars.
Pression absolue et pression relative :

La pression absolue est définie par rapport à la pression dans le vide qui correspond à la
pression nulle. On en déduit donc que la pression minimale possible est zéro.
La pression relative se définit par rapport à une référence que l’on choisi le plus souvent
égale à la pression atmosphérique. Cela consiste finalement à faire une translation du
repère des pressions. La pression nulle est donc équivalente à la pression atmosphérique
(pa). La pression minimale correspond donc à : -pa (pression atmosphérique négative).
Prenons par exemple un réservoir où la surface libre est à la pression atmosphérique (pa).
En écrivant l’équation de l’hydrostatique par rapport à un plan de référence, on a :
p + .g.z = Cte
Entre le point 1 et 2, on a :
p1 + .g.z1 = p2 + .g.z2 = Cte,
soit : p1 = p2 + .g.(z2 -z1) = pa + .g.h, avec pa=105 Pa.
La pression en 1 est mesurée en pression absolue.
La figure suivante montre la répartition des pressions suivant la hauteur.
Pression
atmosphérique : Pa
2
h + =
1
p1=gh pa pa+p1
z1 z2
Plan de référence
Dans la pratique, on préfère souvent mesurer par rapport à la pression atmosphérique (pa).
Dans ce cas, la pression au point 1 devient : p1’ = p1 - pa = .g.(z2 -z1) = .g.h. La pression
en 1 est mesurée en pression relative par rapport à la pression atmosphérique.
La plupart des instruments de mesure fournit une pression relative appelée également
pression manométrique.

3. - VARIATION DE LA PRESSION DANS UN FLUIDE

COMPRESSIBLE
Contrairement aux liquides, les gaz sont fortement compressibles. La variation de la
masse volumique ρ dépend de pression p et de la température T :   f  p, T  . S’il s’agit
p
d’un gaz parfait, nous avons :  rT , r étant une constante. On peut ainsi définir à

p0
l’altitude z0 :  rT0 , ce qui nous permet de calculer la constante r.
0
dp dp p T0
En reprenant l’équation de l’hydrostatique : -.g  on a :  0 g .
dz dz p0 T
Si on se place à température constante T=T0, l’équation peut ainsi s’écrire :
dp g  g 
  0 dz => p=p0 exp   0 z 
p p0  p0 
Dans le cas de l’air, p0=1,013 105 Pa et 0=1.205 kg/m3, l’équation précédente permet de
caractériser l’évolution de la pression dans l’air en fonction de l’altitude :
 g   z 
p=p0 exp   0 z   1, 013.105 exp   
 p0   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 caractériser la pression du fluide sur une
surface. Pour cela, on a besoin de :
- l’intensité : la pression dépend de la profondeur d’eau h. Elle est calculée par la
relation : p = .g.h,
- la zone d’application : la pression s’applique sur une surface (ds),
- la direction : la pression est toujours perpendiculaire à la surface d’application.
p
h
ds surface gauche
surface d’application
de la pression
Le calcul des forces hydrostatiques sur une surface quelconque plongée dans l’eau,
consiste à déterminer les trois caractéristiques suivantes :
4. Forces hydrostatiques sur les parois 29
- l’intensité de la force qui s’applique sur la surface ds : dF = p.ds = .g.h.ds,

- le point d’application de la force,
- la direction.
Dans le but de fournir des résultats facilement applicables, on distingue les surfaces planes
et les surfaces gauches.
4.1. - Paroi plane en position inclinée

On s’intéresse aux surfaces planes de forme quelconque entièrement immergée dans l’eau.
La figure suivante représente à gauche la surface immergée et à droite une vue A-A de
cette surface. On définit un repère (x,y) dont l’axe (x) est sur la surface libre et (y) dirigé
vers le bas et passant par la surface plane. Le point G(xG,yG) est le centre de gravité de la
section. On définit le repère (ξ, η) comme étant une translation du repère (x,y) centré en G.
L’intensité de la force résultante agissant sur la surface S est définie par :
F   dF  .g  hds
S S
L’intégration de cette équation s’écrit :

F = .g.hG.S,
hG : hauteur d’eau du centre de gravité de la paroi immergée,
S : surface de la paroi immergée.
x
h=y.sin(α)
hP hG
Vue A-A x
y
yG Surface S
ds yP
Ixx
F
ds
ξ
G
P
y
η Iξξ
Vue A-A y
Le point d’application de la force résultante des pressions P(xp, yp) est appelé : centre de
pression ou de poussée.
La position de ce point est définie par la position du barycentre des surfaces élémentaires
(ds) pondérées par la pression sur chaque surface, ce qui revient à calculer le moment
équivalent des forces de pression, c’est-à-dire :

 x.dF  x .F
S
p
 y.dF  y .F
S
p
Dans le grande majorité des cas les surfaces sont symétriques par rapport à l’axe η, ce qui
revient à dire que : xp = xG.
La deuxième intégrale s’écrit :
S y ds
2
I
yp   yG 
 yds
S
y G .S
Iξξ représente l’inertie de la section suivant les axes ξξ.

Le tableau suivant fournit le centre de gravité, la surface et l’inertie pour quelques formes
de surface plane.
2h h
v ; v 
v 3 3
h bh bh 3
ξ G ξ S ; I 
v’ 2 36
b
h h
v ; v 
2 2
v bh 3
h ξ G ξ S  bh ; I  
12
v’
b
b h  2B  b  h  B  2b 
v   ; v   
3 B b  3 B b 
h 3  B2  4Bb  b 2 
v
h
h
ξ G ξ S   B  b  ; I  
v’ 2 36  B  b 
d a 3 d a 3
v  ; v  
2 2 2 2
v 2 2
G ξ d 3 3a 3 5d 4 5a 4 3
d ξ S  ; I   
2 2 48 3 16
v’
a
4. Forces hydrostatiques sur les parois 31
 4sin 3    
v v  R 1  
ξ G ξ  3  2  sin  2   
v’  
φ
R v  R 1  cos      v
R2
S
2
 2  sin  2 
R 4 1  cos  2  
3
R4
I   4  sin  4   9 2  sin 2
16   
v  R ; v  R
v
R 4
G ξ S  R 2 ; I  
ξ 4
v’
4.2. - Paroi à surface gauche

La première méthode permettant le calcul de la force résultante F sur une paroi gauche
consiste à décomposer la force élémentaire dF suivant les axes x et y :
dF  dFx .x  dFy .y
L’évaluation des projections de dF suivant x et y, peut se faire de la manière suivante :
dFx  cos().p.ds  cos()..g.h.ds
dFy  sin().p.ds  sin()..g.h.ds
h1
1 3
h2
dFx
2 dFy
dF α
y
x
L’intégration dFx et dFy sur toute la surface de l’élément courbe permet d’évaluer le force
résultante F. Compte tenu de la surface courbe, l’angle α est variable, ce qui complique le
calcul de l’intégrale.

Une deuxième méthode consiste à isoler un volume de fluide et à faire l’équilibre des
forces extérieures agissant sur ce volume.
Dans l’exemple suivant, le volume de fluide isolé est composé d’un ensemble de surfaces
planes horizontales et verticales et de la surface gauche (23). Le choix des surfaces planes
se justifie par l’utilisation des relations précédentes.
h1 Fv
p1 p2
1 3
h2 W
Fx
2
p3
y FR
x
p3
En faisant l’équilibre des forces suivant l’horizontale, on en déduit que la composante
horizontale de la force hydrostatique (FR) est donnée par Fx. La composante verticale est la
somme de la force Fv et du poids de l’eau W. En faisant la composition vectorielle des
forces Fx et Fv, on en déduit FR. Il suffit d’écrire ensuite le moment des forces par rapport
un point quelconque pour localiser la position de FR.
5. - FORCES HYDROSTATIQUES SUR DES CORPS IMMERGES

5.1. - Forces d’Archimède
Supposons qu’une surface fermée formant un corps solide de poids volumique s, de
volume total V et de volume immergé Vimmergé, se trouve immergée entièrement ou
partiellement (Vimmergé  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 : Vimmergé. 
(Appelées forces d’Archimède ou forces de portance)
La force d’Archimède est appliquée au centre de gravité du liquide déplacé (centre de
poussée). Pour un corps plein de poids volumique homogène et entièrement immergé, le
centre de poussée est confondu avec le centre de gravité du corps. Il n’en est pas de même
pour les corps flottants.
5.2. - Equilibre des corps immergés

Un corps est en équilibre si le poids W et la force d’Archimède sont égaux, opposés et
situés sur la même ligne verticale. Dans le cas contraire, il en résulte un mouvement.
La stabilité peut se définir de la façon suivante : si on incline un corps d’un angle par
rapport à la verticale, le corps est soumis à un couple de redressements qui le fait tourner
jusqu’à ce qu’il revienne à sa position initiale.
4. Forces hydrostatiques sur des corps immergés 33
L’instabilité est donc définie par un couple qui tend à augmenter l’inclinaison.
G
P
G
P
Déséquilibre
volontaire
G
Déséquilibre P Mouvement
volontaire Mouvement naturel de la
naturel de la quille
quille
G P
Equilibre stable Equilibre instable
Déséquilibre Mouvement Déséquilibre Mouvement

volontaire naturel volontaire naturel
Dans cet exemple, on constate que la position d’équilibre stable est vérifiée pour un angle
Θ’, et que la position instable correspond à un angle Θ’’ qui a la particularité : Θ’< Θ’’.
On en déduit qu’il existe donc un angle limite Θlimite de basculement entre l’état stable et
instable.

Chapitre III :
L’HYDRAULIQUE EN CHARGE
La dynamique des fluides consiste à étudier le mouvement des particules fluides soumises
à un système de forces. Bien souvent, on commence par l’étude les fluides fictifs dit
« fluides parfaits ». Ils ont la particularité de ne pas avoir de viscosité et de ne pas
développer de la turbulence. Ils permettent d’établir l’équation de Bernoulli facilement.
Les fluides réels engendrent des forces de frottement dues à la turbulence et à la viscosité.
La présence de ces forces induit une perte de charge (énergie) qui est une transformation
irréversible de l’énergie mécanique en énergie thermique.
Ce chapitre aborde, dans un premier temps, l’équation de continuité et surtout
l’établissement de l’équation de Bernoulli. Dans un deuxième temps, l’évaluation des
pertes de charges ainsi que les méthodes de calcul des réseaux hydrauliques en charge sont
exposées.
1. - L’EQUATION DE CONTINUITE
Cette équation exprime le principe de conservation de la masse : la variation de
masse de fluide d’un élément de volume dv pendant un temps dt est égale à la masse de
fluide entrante dans ce volume déduite de la masse de fluide sortante.
z
(x,y,z) dz
y
dy
x dx
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
1-Equation de Continuité 35
  
Après un temps dt dans ce même volume, la masse est égale à :    dt dx.dy.dz
 t 

On constate donc une variation de masse de : dt.dx.dy.dz
t
Variation de masse de fluide pendant la durée dt entre ce qui entre et ce qui sort :
On défini le débit massique par : q massique suivant x  .Q volumique suivant x  .u.dy.dz
   q massique suivant x  dx     q massique suivant x  dx 

 q massique suivant x  .  q massique suivant x  q massique suivant x  . 
 x 2   x 2 
   
1 2
La masse de fluide entrante par la face 1 (suivant x) pendant le temps dt est :

   q massique suivant x  dx     u.dydz  dx 
 q massique suivant x  .  .dt   u.dydz  .  .dt
 x 2   x 2 
 
La masse de fluide sortante par la face 2 (suivant x) pendant le temps dt est :
   q massique suivant x  dx     u.dydz  dx 
 q massique suivant x  .  .dt   u.dydz  .  .dt
 x 2   x 2 
 
La différence de masse entre ce qui entre et ce qui sort dans le volume dv est donc suivant
u
x:  dx.dy.dz.dt
x
u
v w
On a de même suivant y et z avec V   v  :  dx.dy.dz.dt ;  dx.dy.dz.dt
y z
 w 
En écrivant que la variation de masse de fluide d’un élément 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
C’est l’équation de continuité d’un fluide conservatif.
Cas des fluides incompressibles :


Si le fluide est en mouvement permanent (  0 ) la masse volume  est indépendante du
t
temps et si le fluide est incompressible,  est indépendant de x, y et z, l’équation de
continuité se réduit à :
div(V) = 0

36 III- Hydraulique en charge
L’interprétation physique de cette équation est la suivante : les débits Q entrant et sortant à
travers un volume quelconque et rempli du fluide doivent être égaux.
Exemple de conduite à section variable :

(Intégration de div(V) sur le volume et passage à une intégrale de surface par Ostogradski)
U1.S1 = U2.S2
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
développer de la turbulence. Ils permettent d’établir l’équation de Bernoulli facilement.
2.1.1. - Equations d’Euler
 p dz 
p   dx.dy
 z 2 
.g.dx.dy.dz
p(x,y,z)
 p dz 
p   dx.dy
 z 2 
En raisonnant, dans un premier temps, suivant la verticale (z), les forces qui agissent sur
cet élément de volume dv = dS.dz, sont :
 Les forces de volumes :
- les forces de pesanteur provenant de la gravité : - .g.dv
2-Equations de Bernoulli 37
dw  w w x w y w z 
- Les forces d’inertie :  dv      dv
dt  t x t y t z t 
w
- les forces d’accélération pure : dv
t
- les forces d’accélération convective : .grad(w).V.dv
 Les forces de surfaces :
 p dz 
- les forces de pression : Surface inférieure :  p   dx.dy
 z 2 
 p dz 
: Surface supérieure :   p   dx.dy
 z 2 
- les forces de frottement de viscosité :0
- les forces de frottement de turbulence :0
En écrivant l’équilibre de l’ensemble des forces :  F  Forces d' inertie

 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 l’équation 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
V 1
soit :  grad (V ).V  grad (p  .g.z)  0
t 
Cette équation est appelée l’équation d’Euler.
Recherche d’une solution

Les caractéristiques du mouvement d’un fluide incompressible dépendent de quatre
inconnues : Vx, Vy, Vz et de p.
 L’équation d’Euler fournit trois équations.
 L’équation de conservation de la masse fournit une équation : div(V) = 0
 Pour obtenir une solution, il faut également définir un nombre suffisant de conditions
aux limites :
- à la surface libre : p(x, y, z, t) = Cste = pa.
- à la paroi : V doit être parallèle à la paroi.

Utilisation des équations d’Euler pour le calcul des forces hydrauliques sur une
surface
Le théorème d’Euler permet de calculer la réaction de l’eau sur un élément. En raisonnant,

suivant un tube de courant en régime permanent, les forces qui agissent sur cet élément
sont :
- Les forces de volumes : - les forces de pesanteur provenant de la gravité : .Vvol .g

- les forces d’inertie :
- les forces d’accélération pure :0
- les forces d’accélération convective :
V
V . s V.dv  S .V.(V.n ext ).ds  j .V j.(V j.n ext j ).S j
- Les forces de surfaces : - les forces de pression sur les surfaces Si :  p .S
i
i i
- les forces de frottement de viscosité :0

- les forces de frottement de turbulence :0
- l’action du fluide sur l’élément :R
En écrivant l’équilibre de l’ensemble des forces :  F  Forces d' inertie

.Vvol.g   pi .Si  R   .V j.(V j.n ext j ).S j
i j
.Vvol .g  p1 .S1 .n ext 1  p 2 .S 2 .n ext 2  R  .V1 .S1 .V 1  .V2 .S 2 .V 2

2.1.2. - Equations de Bernoulli

Le système d’équations intrinsèques consiste à écrire les équations d’Euler en régime
permanent (  t  0 ) dans un repère particulier. Ce repère 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 :
 
  dV dV  ds
V  Vs et  sV
dt dt dt
  
d s d s ds n
avec :  .  V
dt ds dt R 

R : rayon de courbure et n le vecteur perpendiculaire à s .
V 1  
V  (.g.h  p ) suivant s
s  s
V 1  
V.   (.g.h  p) suivant n
R  n
Avec, V : l’intensité de la vitesse d’une particule fluide.
 En restant sur la ligne de courant on peut intégrer la première équation et l’on

obtient l’équation de Bernoulli :
V2
  .g.h  p  Cste
2
L’équation de Bernoulli est valable en tout point du fluide incompressible en
mouvement permanent.
 En se plaçant sur une ligne normale à la ligne de courant et un rayon de courbure

infini (ligne droite), on obtient :
.g.h  p  Cste

2.1.3. - Représentation graphique

Les différents termes de l’équation de Bernoulli peuvent être écrits en quantités
homogènes à des hauteurs de liquide de poids volumique  :
V2 p pt W E
 h    Cste = 
2g   mg g
hauteur due cote du hauteur due charge totale
à la vitesse point à la pression
Hauteur piézométrique
En terme de pression :
V2 W
  .g.h  p  pt  Cste =
2 Vvolume
Pression pression statique pression totale
dynamique
2.1.4. - Interprétation énergétique
L’équation de Bernoulli peut donner lieu à une interprétation énergétique :

V2
  p*  pt  Cste
2
Energie Energie Energie
cinétique potentielle mécanique totale
L’équation mécanique totale contenue dans un volume unitaire est donc le travail
mécanique total que la particule est susceptible d’accomplir. Ainsi, on peut dire que
l’équation de Bernoulli traduit la conservation de l’énergie mécanique totale par unité de
volume au cours du mouvement permanent.
Dans un écoulement, l’énergie mécanique totale par unité de volume de fluide, peut être
modifiée d’une section à l’autre en introduisant par exemple dans le circuit une machine
hydraulique.
Ainsi l’expression de l’équation de Bernoulli s’écrit :
V12 p1 V22 p E
 h1 + =  h2 + 2 
2g  2g  g
Si l’échange d’énergie se fait des parois de la machine vers le fluide nous avons affaire à
une pompe, si au contraire, l’échange d’énergie se fait du fluide au parois de la machine,
nous avons affaire à une turbine.
L’équation d’énergie est modifiée par le terme E, qui représente l’augmentation par une
pompe ou la diminution par une turbine de l’énergie mécanique totale par unité de masse
de liquide en mouvement.
2.1.4.1. - Turbine
Exemple d’installation :
Travail
La puissance est définie par : Puissance 
Temps

Pour une turbine la puissance hydraulique est dissipée par :

  V 2 
 Volume.   p *
Variation d'énergie E  2 
Pt    
Unité de temps t t

Volume.  V 2


 p *   g.Q.H
 
t  2 
E 
Pt  .Q. t   .Q.H 
g
Types de turbines :
Turbine Pelton Turbine Francis

Turbine Kaplan
Choix d’une turbine :

2.1.4.2. - Pompes
Exemple d’installation :
Principe de montage :
Types de pompes :
Choix d’une pompe :
Travail
La puissance est définie par : Puissance 
Temps
 Ep 
Pour une pompe la puissance hydraulique est fournie par : Pp  .Q.   .Q.H 
 g 

Exemple de courbe caractéristique pour une pompe multicellulaire à étage (de 1 à 12

étages) :
H= Charge totale en sortie – Charge totale en entrée
Caractéristiques par étage :
Pression
absolue
Définition du NPSH « Net Positive Suction Head » (charge nette à l’aspiration)
La pompe transmet une vitesse au fluide par l’intermédiaire de la roue ou de l’hélice. La

vitesse du fluide qui entre dans la roue augmente et par conséquent la pression diminue à
l’entrée engendrant une aspiration et le maintien de l’amorçage. La rotation à l’intérieur de
la pompe augmente la vitesse du fluide tandis que la force centrifuge qui le comprime sur
la périphérie augmente sa pression. Lorsqu’à la suite de survitesses locales des zones de
sous-pression apparaissent sur les bords d’attaque des aubes des roues de pompes
centrifuges, la pression descend en dessous de la pression de vapeur saturante du liquide,
ce phénomène entraîne la formation de bulles de vapeur. L’eau utilisée dans les
installations industrielles contient de l’air dissous dans une proportion atteignant 25mg par
litre. Lorsque l’eau est soumise à une pression décroissante, on observe le dégagement de
l’air dissous quand la pression atteint la pression de saturation et une vaporisation quand la
pression atteint la pression de vapeur. Un phénomène de dégazage puis de vaporisation
partielle du liquide vont apparaître à l’intérieur de la pompe. Les micro-bulles formées
sont transportées par le fluide dans les zones de plus haute pression où elles implosent au
contact des pièces mobiles. Un tel phénomène s’accompagne de variations énormes de
pression (plusieurs milliers de bars) localisées sur de très petites surfaces. Il en résulte une
érosion intense des surfaces métalliques et une chute de pression dans toute l’installation.
Le critère qui sert à éviter la cavitation dans une pompe est le NPSH.
NPSH= Charge totale à l’entrée de la pompe – pression de vapeur saturante

V 2 entrée pompe Pentrée pompe

NPSH    h vapeur saturante
2g g
(Le plan de référence est pris au niveau de l’entrée de la pompe)
Cette caractéristique est fournie par l’installateur dans le cas ou l’on parle de NPSH
disponible. Cette expression est indépendante de la pompe utilisée et ne dépend que de
l’installation (Longueur et diamètre d’aspiration, pertes de charge le long de cette
canalisation, hauteur géométrique d’aspiration).
On parle de NPSH requis pour la charge nette à l’aspiration fournie par le constructeur de
la pompe. Cette relation est indépendante de l’installation de la pompe et ne dépend que
de ce qui se passe dans la pompe.
Pour que l’installation puisse correctement fonctionner il faut : NPSHdisponible > NPSHrequis
Caractéristiques d'une pompe à différentes vitesses de rotation
Les pompes centrifuges vérifient les lois de similitude qui à partir d'une courbe
caractéristique établie pour une vitesse de rotation N de la roue de la pompe permettent
d'obtenir la caractéristique pour une vitesse de rotation N' quelconque.
Si on connaît pour une vitesse N1, le débit Q1, la hauteur manométrique totale Ht1 et la
puissance absorbée P1, on sait qu'il existe deux courbes caractéristiques (Ht en fonction de
Q et P en fonction de Q) pour la vitesse N2 tels que :
;
On peut ainsi reconstruire point par point les caractéristiques pour la vitesse de rotation N2
en prenant des points différents des caractéristiques établies pour la vitesse N1.
2.2. - Ecriture de Bernoulli pour les fluides réels
Nous avons vu que pour le cas d’un fluide réel et en régime permanent, d’autres forces
interviennent, notamment les forces dues au frottement, qui font apparaître une dissipation
de l’énergie mécanique en énergie thermique.
Le long d’une ligne de courant, l’énergie mécanique diminue du fait de l’augmentation de

l’énergie thermique ; on appelle ce phénomène la perte de charge due aux frottements dans
un liquide.
Après intégration entre deux points sur la même ligne de courant l’équation de Bernoulli
s’écrit :
2 2
V1 p V2 p
 1  h1   2  h 2  j1  2
2g  2g 
j1 2 : est la perte de charge entre 1 et 2 en (m).

Toute la difficulté consiste à déterminer cette perte de charge.

Une représentation graphique peut être visualisée par la figure suivante :
3- Evaluation des pertes de charge 51
3. - EVALUATION DES PERTES DE CHARGE

L’écoulement d’un fluide réel dans une conduite représente une des applications
classiques de l’hydrodynamique théorique et expérimentale. La maîtrise de ces
écoulements pour le calcul du fonctionnement des réseaux d’eau potable est indispensable.
On distingue deux types de perte de charge :
- la perte de charge linéaire représentant l’énergie perdue entre les deux points,
- la perte de charge singulière qui intervient lorsque l’écoulement uniforme est
localement perturbé.
3.1. - Les pertes de charge linéaire
3.1.1. - Profil des vitesses dans une section circulaire en charge

Dans le cas d’un écoulement permanent   0 , le profil des vitesses dans une section est
t
dépendant de l’intensité de la vitesse et de la géométrie de l’ouvrage. En effet, nous avons
vu que la turbulence apparaît quand la vitesse dépasse un seuil (expérience de Reynolds).
Ceci peut être représenté par la figure suivante :
Dans la zone laminaire seules les forces de frottement interviennent. Dans la zone
turbulente les forces de turbulence deviennent prépondérantes et une sous-couche
visqueuse très mince apparaît.
Dans le cas des conduites en charge où l’écoulement 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ù l’écoulement 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 très faible
et le profil est aplati en son centre. La zone de la couche limite est de très faible épaisseur
et proche de la paroi. Le gradient de vitesse est très important et les forces de viscosité
sont donc très importantes dans cette couche. Une approche mathématique exacte n’est
pas possible dans ce type d’écoulement complexe. Les relations qui seront construite sont
toutes d’origines empiriques.
Deux types d’écoulements turbulents dans les conduites réelles (présentant des aspérités)
existent :
- conduite hydraulique lisse : l’écoulement est séparé de la paroi par un film
laminaire,
- conduite hydraulique rugueuse : la sous-couche laminaire disparaît et la
turbulence arrive jusqu’à la paroi.
Décrivons le profil de l’écoulement le long d’un diamètre de conduite hydraulique lisse.

On distingue habituellement deux régions : une région centrale dite externe qui correspond
à plus de 90% de l’écoulement et une région voisine de la paroi dite interne.
Dans la région centrale, très aplatie, la vitesse est indépendante des aspérités de la paroi et
dépend peu de la viscosité. L’écoulement est formé de tourbillons importants, allongés,
contenant des tourbillons plus petits. L’intensité de la turbulence est à peu très constante
dans ce domaine.
Dans la zone interne, voisine de la paroi, la vitesse passe d’une valeur élevé à zéro sur une
distance relativement faible. Le fort gradient de vitesse qui en résulte impose au fluide des
contraintes de cisaillement élevées et donc d’intenses frottements visqueux. On comprend

donc que la dissipation d’énergie qui entraîne une partie de la perte de charge ait lieu
principalement dans cette zone. En revanche, l’intensité élevée 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 résulte
des propriétés utilisables pour les transferts thermiques.
La turbulence est mesurée par le nombre de Reynolds. Dans le cas des conduites
circulaires le domaine laminaire et turbulent est différencié par une valeur de Reynolds de
:
La variation du débit ou de la vitesse dans un tuyau de diamètre constant provoque une

perte de charge hr. Cette dernière varie en fonction de la vitesse à peu près linéairement
quand l’écoulement reste laminaire et quadratiquement quand l’écoulement est turbulent.
Entre les deux types d’écoulement il y a une zone de transition, c’est-à-dire une zone
d’incertitude où il faut faire preuve de beaucoup de réserve.
3.1.2. - Expression de la perte de charge linéaire
Compte tenu des difficultés pour résoudre l’équation de Navier-Stokes, la perte de charge
sera traduite par une équation empirique du type :
L U2 H U2
H   ; J 
D 2g L 2gD
 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 diamètre
U la vitesse moyenne
J la perte de charge part unité de longueur

3.1.2.1. - En régime laminaire Re < 2000

En régime laminaire, seules les forces de viscosité interviennent. La vitesse est très faible
et l’état de surface de la paroi n’intervient pas dans le calcul du coefficient.
64
On détermine ainsi :  
Re
32 U
Pour une conduite circulaire, on a : J 
g D2
3.1.2.2. - En régime turbulent Re > 4000

Une surface rugueuse implique un état de surface telles que ses irrégularités ont une action
directe sur les forces de frottement, c’est-à-dire qu’au delà d’un certain degré de rugosité
la configuration géométrique de la surface exerce une influence bien déterminée sur
l’écoulement. Dans le cas des écoulements en charge appliqués au domaine de l’ingénieur,
la rugosité va jouer un rôle très important.
Une surface rugueuse peut être considérée comme étant constituée par une série de
protubérances élémentaires. Elles sont caractérisées par une certaine hauteur k (désignée
parfois par )
Cette rugosité peut être mesurée par rapport au diamètre de la conduite. On parle alors de
k 
rugosité relative : 
D D
Si la rugosité est ondulée le film laminaire pourra, dans une certaine mesure, se modeler
sur les ondulations. Si au contraire elle présente des arêtes vives, le film sera aisément
écorché puis déchiré par les aspérités, l’influence de la viscosité dans la couche limite sera
diminuée et la turbulence fortement augmentée. Ceci se traduit par une vitesse moyenne
réduite.
Les surfaces auxquelles on a affaire dans la pratique (béton, acier, fonte) ne présentent pas
des caractéristiques de rugosité uniforme. En effet, les protubérances de la surface sont
inégales et irrégulièrement distribuées. Elles proviennent du rivetage, des recouvrements,
des joints, des défauts sur la paroi, de corrosion, d’incrustations...
On définit ainsi deux types de rugosité : lisse et rugueuse. La différence entre les deux
dépend du nombre de Reynolds.
Les tableaux suivants donnent quelques indications sur la valeur de .


Le calcul de la perte de charge linéaire va donc prendre en compte le nombre de Reynolds

et la rugosité de la conduite. Le choix se fait en fonction de la turbulence c’est-à-dire de la
valeur du nombre de Reynolds.
Régime turbulent en tuyau lisse :

.Re. 
4000 < Re et  12 : Dans cette zone, la turbulence est encore modérée.
D
L’épaisseur de la sous-couche limite est suffisante pour englober toutes les aspérités de la
conduite qui se comporte dès lors comme un tuyau lisse.
Deux expressions empiriques sont souvent utilisées:

Formule de Blasius :   (100 Re)0.25
1  Re  
Formule implicite de Von Karman :  2. log10  

  2 . 51 
Régime turbulent en tuyau rugueux :

.Re.  k 
 200 : La turbulence devient très importante et  ne dépend plus que de 
D D D
1   
Formule de Nikuradse :  2. log10  
  3,71D 
Régime turbulent en tuyau lisse, rugueux et en zone de transition :

Tous les résultats de ces travaux sont résumés dans la formule de Colebrook qui s’est
fortement inspiré des résultats de Von Karman et de Nikuradse.
1  2.51  
Formule de Colebrook :  2. log10    pour 108 > Re > 4000
  Re  3.7 D 
Cette relation implicite est difficile à exploiter analytiquement et est le plus souvent
représentée sur un graphique.
f
ks/D

On pourra utiliser les tables de Lechapt et Calmon en 1965 en remplaçant la formule de

Colebrook par une formule approchée de la forme : J = L.QM.D-N où L, M et N sont des
constantes pour une rugosité donnée. Pour des valeurs de vitesse comprises en 0.4m/s et
2m/s cette formule conduit à des écarts relatifs inférieurs à 3% par rapport à la formule de
Colebrook.
J mm / m  L
Qm / s
3 M
DmN
k [mm] 2 1 0,5 0,25 0,1 0,05 0,025
L 1.863 1.601 1.4 1.16 1.1 1.049 1.01
M 2 1.975 1.96 1.93 1.89 1.86 1.84
N 5.33 5.25 5.19 5.11 5.01 4.93 4.88
La formule d’Hazen et Williams est également souvent utilisée :

Q1.85
J  10.69 1.85
; Coefficient de Hazen et Williams fonction de k.
CHW .D4.87
3.2. - Les pertes de charge singulière
La perte de charge singulière, localisée dans une section de la conduite, est provoquée par
un changement de direction et d’intensité de la vitesse (voir premier chapitre).
L’écoulement uniforme est perturbé et devient localement un écoulement non uniforme.
La turbulence joue un rôle considérable, alors que les forces de viscosité sont
négligeables. La perte de charge n’a donc lieu qu’en régime turbulent.
Une telle non-uniformité de la vitesse peut être provoquée par :
- un branchement de section de la conduite,
- un changement de direction (coude),
- un branchement ou raccordement,
- un dispositif de mesure et contrôle de débit...
Comme pour les pertes de charge linéaire, les pertes de charges singulières se traduisent
par la relation :
V2
H  K
2g
K est fonction des caractéristiques géométriques et du nombre de Reynolds.
La valeur de K est donnée pour les différents cas les plus classiques dans les tableaux
suivants :

Raccordement d’une conduite avec un grand réservoir

Départ
Sans saillie à l’intérieur du réservoir, avec raccordement à angles vifs
K  0 .5
Sans saillie à l’intérieur du réservoir, avec raccordement à angles vifs, ajutage débitant à
gueule bée
K 1
Avec saillie à l’intérieur du réservoir
K 1
Pour une saillie dont la longueur est comprise

entre 1 et 2 fois le diamètre
Sans saillie à l’intérieur du réservoir, avec raccordement de profil arrondi
K  0.05
Cette valeur est une moyenne, elle dépend du

profil de l’arrondi.
Arrivée
S2
K 1
S1 2 2
 S  1 S 
K  1  1    1 
 S 2  9  S2 
S2>>S1 => K=1
Coudes
Arrondi
 7

  d 2  
K   0.131  1.847  
  2.r   90
 
 en degrés
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 indépendant du diamètre
 22.5 30 45 60 90
K 0.07 0.11 0.24 0.47 1.13

Tés
Branchement de prise à 90° de même diamètre et à angles vifs
2
V
H r  K r t
2g
2
Vt
H b  K b
2g
Qb/Qt Qb/Qt 0 0.2 0.4 0.6 0.8 Qb/Qt 1

Kr 0.40 0.26 0.15 0.06 0.02 0.00
Kb 1.00 1.01 1.05 1.15 1.32 1.45
Branchement d’amenée à 90° de même diamètre et à angles vifs

2
V
H r  K r t
2g
2
V
H b  K b t
2g
Qb/Qt Qb/Qt 0 0.1 0.2 0.4 0.8 Qb/Qt 1

Kr 0 0.16 0.27 0.46 0.60 0.55
Kb -0.60 -0.37 -0.18 0.26 0.94 1.20
Cônes
Convergent
La perte de charge est négligeable.
Divergent
2

   
1.25
  D 2 
K  3,2. tg   .1   1  
  2    D2  
 
Si il y a décollement et le comportement
est identique à celui d’un élargissement brusque.

2
 D1  3 6 8 10 12 14 16 20 24 30 40
 
 D2 
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 diamètre

Rétrécissement
  D 2 
K  0,5.1   2  
  D1  
 
D2 0.01 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
D1
K 0.500 0.495 0.480 0.455 0.420 0.375 0.320 0.255 0.180 0.095
Elargissement
2 2
 S  1 S 
K  1  1    1 
 S 2  9  S2 
D1 0.01 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
D2
K 1.000 0.980 0.922 0.829 0.708 0.569 0.424 0.287 0.175 0.109

Appareils de robinetterie
Vanne opercule
p 1 2 3 4 5 6 7
D 8 8 8 8 8 8 8
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.3. - Exemple de tracé de perte de charge
Considérons une conduite de longueur L12. On dispose en 3 d’un rétrécissement brusque.
Supposons deux conduites horizontales avec une liaison du type élargissement ou

rétrécissement brusque.

4. - POSITION DE LA LIGNE PIEZOMETRIQUE

Considérons une conduite reliant deux réservoirs. La ligne piézométrique correspondant
aux pression relatives est représentée approximativement par la droite AA’ (On a négligé
les pertes de charges dues à l’entrée et à la sortie des réservoirs. La ligne piézométrique
BB’ correspond aux pressions absolues (Pa/ = 10.33m).
Si la conduite toute entière est située au dessous de AA’, la pression dépasse la pression
atmosphérique. Cette hypothèse correspond à une situation normale. Il faut prévoir des
ventouses aux points les plus élevés pour la sortie de l’air accumulé et des décharges de
fond N et N’’ pour la vidange et le nettoyage.
Si la conduite passe au-dessus de la ligne piézométrique AA’, la partie du tronçon au-

dessus de AA’ est en dépression. En général, on doit éviter les zones en dépression ; en
effet, une fente éventuelle facilitera l’entrée de corps étrangers et peut contaminer l’eau.
La pose d’une ventouse normale est contre indiquée, étant donné qu’elle permettrait
4- Position de la ligne piézométrique 69
l’entrée de l’air et provoquerait en conséquence, une réduction de débit. Dans la majorité

des cas, on évite ce type de tracé. Si toutefois le tracé
Si la conduite s’élève au-dessus de la ligne horizontale qui passe par A, il n’y aura
écoulement que si toute la conduite a été remplie d’eau au préalable (effets de
siphonnage).
Si la forme de la conduite est toujours située au-dessous de l’horizontale passant par A,

mais dépassant BB’, l’écoulement se fait sans siphonnage. Cependant, la pression est nulle
pour la partie située au-dessus de BB’ : il y a cavitation. L’écoulement se fait avec une
section partiellement pleine et le débit diminue. Les dégagements d’air et de vapeur d’eau
dans les régions de forte dépression rendent l’écoulement irrégulier.

Si la conduite a un tronçon au-dessus de l’horizontale passant par A, il faudrait amorcer le

siphon, l’écoulement se ferait en section partiellement pleine et de façon très irrégulière.
Si la conduite dépasse la cote B, il est impossible d’amorcer l’écoulement.
5. - LES CONSTRUCTIONS GRAPHIQUES
Ces constructions permettent de déterminer le point de fonctionnement d’un réseau (débit,

pression) en utilisant les courbes de fonctionnement et les courbes caractéristiques des
pompes.
Prenons l’exemple suivant :
5- Les constructions graphiques 71
Déterminons dans un premier temps la caractéristique du réseau en fonction du débit :

CR(Q).
CR(Q)  H Sortie pompe  H Entrée pompe
Cette caractéristique doit être déterminée en passant uniqiement par le réseau.
Déterminons, dans un premier temps, la charge totale (en mCE) au point B, HB, par
application du théorème de Bernouilli.
H1  H B  PdC1B
2 2
p V p V
h1  1  1  h B  B  B  j1B
g 2g g 2g
H B  h1  j1 B
Déterminons, dans un deuxième temps, la charge totale (en mCE) au point C, HC, par
application du théorème de Bernouilli.
H C  H 2  PdCC 2
2 2
p V p V
h C  C  C  h 2  2  2  jC 2
g 2g g 2g
H C  h 2  jC  2
On peut ensuite représenter la caractéristique du réseau :
CR(Q)  H C  H B   h 2  h1     PdCaspiration  PdC refoulement  .
Il suffit de tracer la caractéristique de la pompe HMT en fonction du débit (donnée fournie
par le constructeur) : HMT(Q).
Point de fonctionnement d'une instalation
80
70
Point de
fonctionnement de
60 l'installation
50
Courbe
caractéristique de la
H (mCE)
pompe : HMT(Q)
40
30
20
Courbe
10 Hauteur caractéristique du
géométrique de réseau : CR(Q)
l'installation
0
0 10 20 30 40 50 60 70 80
Débit (l/s)
Nous pouvons ainsi déterminer le point de fonctionnement de la pompe et de l’installation.

6. - LES COUPS DE BELIER

6.1. - Introduction
Toute modification du régime permanent de l’écoulement dans un réseau de conduites
(démarrage ou arrêt d’une pompe, fermeture d’une vanne, etc.) est susceptible de
provoquer des phénomènes se caractérisant par des variations brutales de la vitesse et de la
pression ; on parle de « coups de bélier ». Les conséquences néfastes des coups de bélier
sont la rupture de canalisations ou encore la détérioration d’appareils traversés par
l’écoulement. L’ingénieur doit donc être capable de prévoir ces phénomènes hautement
instationnaires et de déterminer leur amplitude. Si besoin, un dispositif de protection doit
être installé.
On distingue le coup de bélier de masse (parfois appelé « oscillation en masse ») et le
coup de bélier d’ondes (appelé simplement « coup de bélier » dans le langage courant).
Dans le premier cas, l’écoulement répond à la modification du régime d’écoulement
uniquement par des transferts de masse ; la masse volumique du fluide est inchangée et la
canalisation ne se déforme pas. Dans le second cas, il y répond par des transferts de masse,
une modification de la masse volumique du fluide et éventuellement par une déformation
de la canalisation.
L’objectif de cette partie est de faire comprendre ces phénomènes complexes au moyen
d’une description mathématique simplifiée.
6.2. - Formulation mathématique

Les deux équations de base de l’hydraulique, à savoir la loi de conservation de la masse
(équation 1) et la loi de la quantité de mouvement (équation 2) vont permettre de
caractériser le coup de bélier.
soit : 0 Première équation de Saint Venant (1)

6- Le coup de bélier 73
avec : l'équation devient :
.
∗
soit : .
t : le temps,
x : l’abscisse curviligne suivant la canalisation (pas nécessairement horizontale),
ρ : la masse volumique de l’eau,
S : la section de l’écoulement (section intérieure de la conduite),
U : la vitesse permanente en moyenne uniformisée sur la section,
P : la pression relative par rapport à la pression atmosphérique,
g : l’accélération gravitationnelle,
j : le coefficient adimensionnel de perte de charge,
i : la pente de la canalisation.
6.3. - Coup de bélier de masse

Dans le cas du coup de bélier de masse, les modifications du régime permanent de
l’écoulement sont suffisamment faibles pour que la masse volumique du fluide ne soit pas
modifiée et que la canalisation ne soit pas déformée. Dans ce cas, l’équation de
conservation de la masse s’écrit comme suit :
0 0
(4)
La précédente é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)
L’équation de la quantité de mouvement s’écrit alors :

∗
1
.
(6)
Appliquons cette dernière équation au calcul des oscillations de niveau d’eau dans une
cheminée d’équilibre (protection anti-bélier), comme illustré sur la Figure 1. Cet exemple
a été proposé par Frelin (Frelin 2002).
Figure 1. Installation avec cheminée d'équilibre (Frelin 2002).
Ici, une conduite circulaire forcée entre une retenue d’eau de section importante et une
installation de turbinage est équipée d’une cheminée d’équilibre de section circulaire.
Cette conduite est le siège d’un écoulement permanent jusqu’à l’instant t = 0. A cet
instant, la vanne située en amont de la turbine est fermée instantanément.
On considère les données suivantes: L = 4000 m, d = 1.5 m, D = 3 m, z0 – z2 = 20 m,
Q0 = 1 m3/s (débit volumique).
On suppose que la masse volumique de l’eau reste constante et que la conduite est
complètement indéformable. Déterminons, sous l’hypothèse d’absence de pertes de charge
(J = 0), la période des oscillations du niveau d’eau dans la cheminée ainsi que le niveau
maximal.
Pour cela, intégrons l’équation (6) entre les abscisses curvilignes des points 1 et 2. Nous
obtenons l’expression suivante :
∆ ∗
0 (7)
Par ailleurs, compte tenu des sections importantes du réservoir et de la cheminée par
rapport à la section de la canalisation, l’application de l’équation de Bernoulli sans pertes
de charge conduit aux équations suivantes :
2 2
p U1 p* U
z0  z1  1   1  (8)
 g 2g  g 2g
2 2
p2 U 2 p2* U
z0  s  z 2     (9)
 g 2g  g 2g
La vitesse est en effet identique au point 1 et au point 2 dans la mesure où nous avons fait
les hypothèses d’un fluide incompressible et d’une conduite indéformable (équation 5).
L’utilisation de ces équations aboutit à l’expression suivante pour la loi de la quantité de
mouvement :
dU
L  gs  0 (10)
dt
Le débit dans la conduite étant identique à celui dans la cheminée, la loi de conservation
de la masse (équation 4) peut s’exprimer, quelque soit le temps t, comme suit :
d2  D 2 ds
U (11)
4 4 dt
La dérivée de la vitesse par rapport au temps s’écrit alors :
2
dU  D  d 2 s
  (12)
dt  d  dt 2
En utilisant cette expression dans l’équation (10), nous aboutissons à l’équation
différentielle suivante, où la seule variable est le temps t :
2 2
D d s
  L 2  gs  0 (13)
 d  dt
Elle admet comme solution la forme générique suivante, où C1 et C2 sont des constantes :
d g  d g 
s  C1 cos  t   C2 sin  t  (14)
D L  D L 
ds 4Q0
Sachant que s et valent respectivement 0 et à l’instant t = 0 (conditions
dt  D2
initiales), nous obtenons le niveau d’eau dans la cheminée en fonction du temps.
4Q0 L d g 
s sin  t  (15)
 Dd g D L 
D L
Ainsi, la période des oscillations vaut 2 , soit 4 minutes et 14 secondes. La
d g
4Q0 L
hauteur maximale dans la cheminée vaut quant à elle z0  z2  , soit 25.7 m.
 Dd g
6.4. - Coup de bélier d’ondes

Les équations du coup de bélier ont été définies précédemment :
Conservation de la masse : 0
∗
Conservation de la quantité de mouvement : .
6.4.1. - Mécanisme de propagation d'une onde

Si dans un milieu continu, au repos, on provoque le déplacement rapide d’un élément on
constate que les éléments voisins se déplacent et agissent à leur tour sur les particules en
contact ; ainsi, le déplacement se propage, de proche en proche, avec une vitesse finie dans
tout le milieu.
Si ce déplacement est de faible amplitude (ce qui sera le cas dans l’étude des coups de
bélier), la vitesse de propagation est appelée célérité du son et on la désigne
habituellement par a.

Dans un fluide, cette propagation s’effectue dans une direction confondue avec celle du
déplacement. Elle s’accompagne, à cette même célérité, d’une variation de pression, de
masse volumique et de température. Ce domaine mobile, par rapport au fluide constitue
une onde.
Pour schématiser plus simplement ce phénomène de propagation d’ondes, considérons une
longue tuyauterie rectiligne, absolument indéformable, de section constante, et contenant
un liquide.
Supposons que l’une de ses extrémités comporte un piston mobile et que l’autre soit
fermée. Dans l’éventualité hypothétique d’un liquide incompressible, le déplacement du
piston se transmettrait instantanément à la totalité du fluide contenu dans la tuyauterie. Il
ne pourrait pas se déplacer sans que la pression qui règne dans le liquide devienne
théoriquement égale à l’infini.
Mais en réalité, comme nous l’avons précédemment évoqué, les liquides ne sont pas
rigoureusement incompressibles. Si l’on déplace brusquement le piston d’une petite
quantité, il apparaîtra immédiatement une augmentation de pression dans la couche de
liquide en contact avec lui. Cette couche, n’étant plus en équilibre avec les tranches plus
éloignées, va se détendre en comprimant à son tour les couches voisines et ainsi de suite
avec une vitesse de propagation égale à a, puisque la tuyauterie est supposée
indéformable.
6.4.2. - Célérité des ondes : cas des canalisations indéformables

Pour déterminer la célérité a considérons, tout d’abord, l’écoulement permanent d’un
liquide parfait contenu à l’intérieur d’une canalisation rectiligne, rigoureusement
indéformable et de section constante S. Si nous provoquons alors une perturbation, il y
aura au passage de l’onde une discontinuité de la pression, de la masse volumique, de la
température et également de la vitesse.
Pendant le temps dt, la masse de fluide concernée par le parcours de cette onde est :
. . .
et elle subit une accélération . Pour cette masse de fluide, l’équation de la dynamique
s’écrit :
. . . . . . soit : . .
À la date t1, cette onde élémentaire se trouve à une abscisse x1 dans la canalisation et, à t2,
elle se trouvera dans une position x2.
t1
t2
Au temps t1, toutes les particules de fluide comprises entre les sections A1 et A2 ont
conservé leur masse volumique ρ et leur vitesse V. Ce qui n’est pas le cas, pour cette
même masse de fluide contenue à la date t2 entre les sections A3 et A2, où ces grandeurs
sont respectivement devenues .
Avant le passage de l’onde, la masse de fluide comprise entre A1 et A2 était :
. . .
Après le passage de l’onde, cette même masse de fluide est comprise entre les sections
A3 et A2 et elle devient :
En égalisant ces expressions et en négligeant les termes du second ordre, on a :

dρ ρd
En utilisant l'équation précédente, on a :
Définition du coefficient de compressibilité :

La compressibilité d'un corps représente la variation de volume du corps en
réponse à une variation de pression :
1
En faisant l'expérience à masse

constante :
.
0 . .
Soit : => => =>

ε = 2050MPa pour l'eau ce qui donne a=1432m/s

Dans ce cas la conservation de la masse peut s'écrire :

Conservation de la masse : 0
avec 0 => canalisation indéformable et , on a : 0 soit :
. 0
6.4.3. - Célérité des ondes : cas des canalisations déformables
P.2.R.Long=σθ.2.e.Long
Dans ce cas la conservation de la masse peut s'écrire :

0 soit 0
1
0
1
Par analogie avec l'équation précédente, la célérité des ondes s'écrit :

1
.
1
Pour les conduites façonnées en matériau peu élastique, la célérité des ondes est en général
comprise entre 700 et 1300 m/s ; elle peut descendre jusqu’à quelques dizaines de m/s
dans le cas de conduites élastiques, façonnées par exemple en caoutchouc ou en plastique
(Bonnin 1983). En considérant un ordre de grandeur de 1000 m/s pour la célérité du coup
de bélier d’ondes, on comprend la rapidité du phénomène dans des conduites peu
élastiques par rapport au coup de bélier de masse.
6.4.4. - Equations des caractéristiques
∗
Le système de Saint Venant peut s'écrire en utilisant : car z est indépendant du
temps :
∗
Conservation de la masse : . 0 (1)
∗
Conservation de la quantité de mouvement : . (2)
∗ ∗
en faisant : 1 2 : .
∗ ∗
en faisant : 2 1 : .
en se déplaçant suivant les trajectoires : et , le système

précédent devient :
∗ ∗
1
. . .
∗ ∗
1
. . .
∗
1
.
∗
1
.
La célérité des ondes étant généralement très grande (de l’ordre de 1,000 m/s) devant la
vitesse moyenne uniformisée (de l’ordre de 1 m/s), nous pouvons approximer U  a par a
et U  a par – a. Considérons de plus qu’il n’y a pas de pertes de charge (j = 0). Enfin,
utilisons comme variables le débit Q transitant dans la canalisation et la hauteur
p*
piézométrique y définie comme la quantité . Dans ces conditions, nous pouvons
g
écrire :

D  aQ  dx
  y   0 avec a (28)
Dt  gS  dt
D  aQ  dx
  y   0 avec  a (29)
Dt  gS  dt
La première équation signifie qu’un observateur se déplaçant à la vitesse a verra la
aQ
quantité  y constante tout au long de son déplacement ; la seconde, qu’un
gS
aQ
observateur se déplaçant à la vitesse – a verra la quantité   y constante tout au long
gS
de son déplacement.
La méthode de l’épure Schnyder – Bergeron consiste à résoudre graphiquement le
phénomène de coup de bélier en tirant profit du résultat précédent. Plus précisément, cette
méthode permet de déterminer la pression et le débit en tout point de la conduite et à tout
instant sans déterminer la nature exacte de l’onde mais simplement en considérant des
« observateurs » se déplaçant aux vitesses a et – a. Nous allons l’utiliser pour quelques
exemples simples afin de mettre en évidence le phénomène de coup de bélier d’ondes.
6.4.5. - Application à la fermeture instantanée d’une vanne (sans pertes de charge dans la
conduite)
Considérons une conduite reliant un réservoir à charge constante et une vanne, comme
illustré sur la Figure 2. Alors que l’écoulement y était permanent (débit Q0) jusqu’à
l’instant t = 0, la vanne est instantanément fermée à l’instant t = 0.
Figure 2. Fermeture d’une vanne (Frelin 2002).
Pour des raisons de facilité de compréhension, on raisonnera à partir du temps

adimensionnalisé par la durée nécessaire à l’onde pour parcourir la longueur de la
conduite :
t
 (30)
La
Par exemple, l’onde générée en V (vanne) à μ = 0 arrivera en R (réservoir) à μ = 1.
L’écoulement permanent consiste en un débit Q0 et une charge z0 – z2 dans toute la
2
U
conduite (pas de pertes de charge). En supposant le terme négligeable dans
2g
l’expression de la charge, nous pouvons affirmer que la hauteur piézométrique y0 en

régime permanent vaut z0 – z2 dans toute la conduite (pas de pertes de charge).
Une fois la vanne fermée, les conditions aux limites sont un débit nul à la vanne et une
hauteur piézométrique égale à y0 au niveau de l’entrée du réservoir.
Figure 3. Epure pour une fermeture instantanée de la vanne (Frelin 2002).
Plaçons-nous en μ = – 1+ au niveau de la vanne (c’est-à-dire une unité de temps

adimensionnel avant sa fermeture). Le signe « + » signifie qu’on se situe à un instant
infiniment petit après le temps μ = – 1 ; cette notation permet d’éviter toute ambigüité due
dx
à la nature discontinue du coup de bélier d’ondes. En se déplaçant selon l’équation a
dt
nous arriverons à la vanne en μ = 0+, c’est-à-dire juste après la fermeture de la vanne. En
μ = – 1+, les conditions du régime permanent règnent : le débit et la hauteur piézométrique
valent respectivement Q0 et y0 au niveau du réservoir. En μ = 0+ à la vanne, la quantité
aQ
 y étant conservée, nous pouvons écrire la relation suivante :
gS
 aQ   aQ  aQ0
  y    y   y0 (31)
 gS  V 0  gS  R 1 gS
Le débit étant nul à la vanne, nous pouvons déterminer la hauteur piézométrique à la
vanne juste après sa fermeture :
aQ0
yV 0   y0 (32)
gS
Cette expression montre que la fermeture instantanée de la vanne est responsable d’une
aQ0
augmentation brutale de la hauteur piézométrique de . Le point V0+ est le point haut
gS
de l’épure (Figure 3). La droite reliant les points R–1 et V0+ a pour équation :
+
aQ aQ0
y  y0 (33)
gS gS
dx
Déplaçons-nous à présent depuis la vanne selon l’équation  a . Arrivant au réservoir
dt
aQ
en μ = 1+, nous voyons conservée la quantité   y.
gS

 aQ   aQ  aQ0
  y    y   y0 (34)
 gS  R1  gS  V 0 gS
La hauteur piézométrique à l’entrée du réservoir valant y0, le débit à l’entrée du réservoir
vaut :
gS  aQ0 
QR1     y0  y0   Q0 (35)
a  gS 
Ainsi, du fait des fortes pressions régnant dans la canalisation, le débit s’inverse : un
écoulement prend place de la conduite vers le réservoir. La droite reliant les points V0+ et
R1+ a pour équation :
aQ aQ
  y  y0  0 (36)
gS gS
dx
Gagnons à présent la vanne en nous déplaçant selon l’équation  a . La quantité
dt
aQ
 y étant conservée, nous pouvons déterminer la pression à la vanne en μ = 2+ :
gS
aQ
yV 2  y0  0 (37)
gS
Ceci signifie que la vanne subit une dépression. La droite reliant les points R1+ et V2+ a
pour équation :
aQ aQ
 y  y0  0 (38)
gS gS
dx
En nous déplaçant vers le réservoir selon l’équation  a , nous pouvons exprimer le
dt
aQ
débit au niveau du réservoir en μ = 3+ en conservant la quantité  y :
gS
QR 3  Q0 (39)
La droite reliant les points V2+ et R3+ a pour équation :
aQ aQ
  y  y0  0 (40)
gS gS
En poursuivant les allers-retours entre le réservoir et la vanne et en procédant de la même
façon à partir de différents points à différents instants, nous pouvons déterminer le débit et
la hauteur piézométrique en tout point de la conduite et à tout instant. Le phénomène
périodique ainsi mis en évidence est illustré sur la Figure 4. La Figure 5 présente quant à
elle la variation de pression au niveau de la vanne ainsi que la variation de débit au niveau
de l’entrée du réservoir. Nous constatons sur cette dernière figure que le coup de bélier
d’ondes a ici une période égale à 4μ, c’est-à-dire 16 s en considérant les données du calcul
effectué plus haut pour le coup de bélier de masse dans une conduite équipée d’une
cheminée d’équilibre et en supposant une célérité de 1,000 m/s. Nous pouvons ainsi
constater la fréquence beaucoup plus importante du coup de bélier d’ondes par rapport au
coup de bélier de masse.
Figure 4. Décomposition du coup de bélier d’ondes due à la fermeture d’une vanne

(Frelin 2002).
Figure 5. Variations de la pression à la vanne et du débit à l’entrée du réservoir (Frelin

2002).
Des figures similaires à la Figure 5 pourraient être tracées en chaque point de la conduite.
La résolution de l’épure Schnyder – Bergeron nous a permis de mettre en évidence les
aQ
pressions maximale et minimale subie par la conduite, respectivement y0  0 et
gS
aQ0
y0  .
gS

Il peut arriver que la pression absolue minimale atteigne la pression de vapeur saturante de
l’eau (0.239 mCE à 20°C). Dans ce cas, le liquide se vaporise ; une poche de cavitation se
forme. La conduite doit être suffisamment résistante pour ne pas imploser sous l’action de
la pression atmosphérique sur sa surface extérieure. La démarche décrite ci-dessus, qui
suppose un fluide sous phase liquide et qui n’est donc plus valable, doit être adaptée
(Pernès 2004).
6.4.6. - Application à la fermeture progressive d’une vanne (sans pertes de charge dans la
conduite)
Considérons à présent une fermeture lente de la vanne, toujours sans considérer de pertes
de charge. Utilisons pour cela l’équation 41 qui constitue un exemple possible de loi de
fermeture de vanne.
2
 Q 
 vanne   0   (41)
1 t T 
Ici, ξvanne est la perte de charge singulière due à la vanne (en mCE) ; ξ0, la perte de charge
singulière de la vanne lorsqu’elle est complètement ouverte ; T, le temps total de
fermeture. L’épure illustrée sur la Figure 6 se construit de la même façon que pour
l’exemple précédent, à ceci près qu’il faut considérer au moyen de l’équation 41 un débit à
la vanne non nul pour tous les instants compris entre le début de la fermeture et la
fermeture complète.
Nous pouvons constater sur cette figure qu’une fermeture lente atténue l’amplitude du
phénomène, aussi bien pour les surpressions que pour les dépressions.
Figure 6. Epure pour une fermeture progressive de la vanne (Frelin 2002).
6.4.7. - Application à la fermeture instantanée d’une vanne (prise en compte simplifiée des
pertes de charge)
Considérons à nouveau la fermeture instantanée de la vanne mais cette fois-ci en
considérant d’une façon simplifiée les pertes de charge dans la conduite. Supposons ainsi
que la perte de charge dans la conduite soit localisée à l’entrée du réservoir (perte de
charge singulière entre les points R et W infiniment proches). Pour construire l’épure,
nous considérerons des allers-retours non pas entre V et R mais entre V et W. La hauteur
piézométrique en W correspond à la hauteur piézométrique en R à laquelle on retranche la
perte de charge correspondant au débit. Dans le cas des débits négatifs, remarquons que la
hauteur piézométrique est plus importante en W qu’en R. L’épure ainsi tracée, illustrée sur
la Figure 7, converge vers son point central, c’est-à-dire un débit nul et la hauteur
piézométrique régnant dans le réservoir : les pertes de charge ayant atténué le coup de
bélier, il n’y a plus d’écoulement.

Figure 7. Epure pour une fermeture instantanée de la vanne avec une prise en compte
simplifiée des pertes de charge (Frelin 2002).
6.5. - Protections anti-bélier

On distingue les anti-béliers spécifiques des anti-béliers non spécifiques (Bonnin 1983).
Parmi les premiers, les vannes à fermeture lente permettent d’empêcher les coups de bélier
dus à une fermeture brutale de vanne ; ces dispositifs sont en revanche inefficaces contre
les arrêts brutaux de pompes, toujours susceptibles de se produire (panne du moteur,
coupure d’alimentation, manœuvre involontaire, etc.).
Toujours parmi les anti-béliers spécifiques, les volants d’inertie sont des dispositifs
capables d’emmagasiner de l’énergie pendant le démarrage de la machine et de la fournir à
la pompe dès la coupure de son alimentation : l’arbre de la pompe continue ainsi à tourner
pendant quelques secondes à quelques dizaines de secondes.
Les anti-béliers non spécifiques ne cherchent pas à empêcher les coups de bélier mais à
atténuer leurs effets. Dans la mesure où ils ne protègent pas entre le lieu de génération du
coup de bélier et leur localisation, ils doivent être installés à des endroits bien spécifiques.
Une soupape est ainsi un dispositif qui s’ouvre à la pression atmosphérique lorsque la
pression dans la conduite devient supérieure à une valeur seuil. Ce faisant, la soupape
laisse échapper un certain débit. Utile contre les surpressions, ce dispositif est en revanche
complètement inefficace contre les dépressions.
Pour atténuer aussi bien les surpressions que les dépressions, une méthode consiste à
tempérer les variations de pression dans la conduite par des variations de masse dans un
réservoir associé. C’est le principe de fonctionnement des cheminées d’équilibres
(réservoir à l’air libre), illustrées sur la Figure 8, et des ballons anti-bélier (réservoir sous
pression), illustrés quant à eux sur la Figure 9. En fonctionnement permanent, le réservoir
est partiellement rempli. Lorsque la conduite est le siège d’une surpression, le réservoir
accepte de l’eau de la conduite de façon à tempérer la surpression. Au contraire, lorsque la
conduite est le siège d’une dépression, le réservoir apporte de l’eau à la conduite de façon
à tempérer la dépression.
Figure 8. Cheminée d’équilibre (Bonnin 1983).
Figure 9. Ballon anti-bélier (Bonnin 1983).

HEC - Cours - MEPA 2014

Transféré par

Droits d'auteur :

Formats disponibles

HEC - Cours - MEPA 2014

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

HEC - Cours - MEPA 2014

Transféré par

Droits d'auteur :

Formats disponibles

2014

Formation Mastère Eau Potable et Assainissement

José VAZQUEZ et Matthieu DUSFRESNE

D’un point de vue réglementaire, la directive 2000/60/CE du Parlement européen établit

Avant d’entrer directement dans la description mathématique et la quantification des

José VAZQUEZ et Matthieu DUFRESNE

CHAPITRE I : CARACTERISTIQUES DES

CHAPITRE II : HYDROSTATIQUE ........................................ 25

José VAZQUEZ et Matthieu DUFRESNE

CARACTERISTIQUES DES ECOULEMENTS

L’hydraulique est l’étude des écoulements.

Evolution de la vitesse ponctuelle en fonction du

Injection Filet de colorant

1.2. - Distribution des vitesses dans un tuyau rectiligne

José VAZQUEZ et Matthieu DUFRESNE

Evolution de la vitesse dans une conduite circulaire

1.3. - Elargissement et rétrécissement dans une conduite

1.4. - Sortie d’un réservoir

Evolution des lignes de courant dans un coude circulaire

José VAZQUEZ et Matthieu DUFRESNE

1.6. - Jonction et bifurcation

L’écoulement dans une bifurcation se comporte de manière légèrement différente du fait

1.7. - Obstacles dans un écoulement en charge

Obstacle dans un écoulement. L’écoulement est à faible L’augmentation légère de la

A vitesse plus importante les L’augmentation de la vitesse a

Vanne à opercule Vanne papillon

1.8. - Le coup de bélier

Vanne A un instant donné, la fermeture brutale de

Onde de Gonflement Vanne A l’amont de la vanne, l’eau se comprime

José VAZQUEZ et Matthieu DUFRESNE

Gonflement Vanne L’onde de choc atteint le réservoir. La

Onde de Vanne La pression dans le réservoir étant beaucoup

U=0 Dès que l’onde atteint la vanne, la conduite

Gonflement Vanne Comme précédemment, le réservoir impose

Onde de Vanne La conduite retrouve alors sa forme initiale.

U=0 L’onde arrive au niveau de la vanne et on

2. - ECOULEMENTS A SURFACE LIBRE

2.2. - Distribution des vitesses dans un canal

José VAZQUEZ et Matthieu DUFRESNE

2.3. - Notion d’influence aval

U=0 U0 U=0 U0 U=0

Expérience : notion d’influence aval

Au niveau de la surface libre, il se crée

Quand la vitesse (U) dans le canal est

De cette expérience, on en déduit que :

2.4. - Courbes de remous

Ressaut fort Ressaut ondulé

2.6. - Les seuils et déversoirs

Seuil dénoyé Seuil très faiblement noyé

José VAZQUEZ et Matthieu DUFRESNE

Seuil faiblement noyé Seuil complètement noyé

En fonction de la forme du seuil et de la vitesse de l’écoulement, il peut apparaître une

3. - PROPRIETES DES LIQUIDES

3.1. - Masse volumique

3.2. - Poids spécifique

José VAZQUEZ et Matthieu DUFRESNE

Considérons un fluide placé entre deux plaques planes, parallèles, distantes de L et

Cette force est la résultante des forces de frottements visqueux.

Le fluide en contact avec le cylindre extérieur va y adhérer et par conséquent va être