COURS DHYDRAULIQUE

T2 : ECOULEMENTS A SURFACE
LIBRE

AL. MAR
Juillet 2004

AVERTISSEMENT

Cette partie du cours dhydraulique est destine aux tudiants de la Formation Initiale
dIngnieurs et des Formations Post-Universitaires (Informatique Applique aux Sciences de
leau, Gnie Sanitaire et Environnement, Eau pour lAgriculture et les Communauts) de
1Ecole Inter-Etats dIngnieurs de 1Equipement Rural. Cest pourquoi nous avons tent dy
dvelopper des aspects de lhydraulique gnrale et de lhydraulique applique lirrigation,
aux ouvrages dart et aux barrages. Les aspects de transports solides et de lhydraulique des
cours deau fond mobile sont omis. Les mthodes de rsolution des formules, dj
dveloppes en hydraulique en charge, ne sont pas reprises.
Le chapitre 1 dfinit lhydraulique surface libre, ses spcificits et les principaux lments
gomtriques et hydrauliques des canaux paramtrs. Les diffrents types dcoulements
quon peut rencontrer et les normes sur laffouillement des canaux non revtus y sont
galement abords.
Le chapitre 2 tudie lcoulement uniforme et donne les principaux rsultats sur les tentatives
de dtermination de la rugosit des canaux qui sont les aspects les plus dificiles de
lhydraulique surface libre. Il constitue la partie la plus commune aux diffrentes filires de
IEIER avec le dimensionnement des canaux.
Le chapitre 3 traite des coulements graduellement varis rencontrs autour des singularits
des canaux. Des mthodes simples dintgration de lquation diffrentielle rgissant ces
coulements y sont dveloppes. Les tudiants intresss aux aspects numriques plus
labors peuvent se rfrer lannexe du chapitre sans jamais oublier que la prcision des
calculs dpend plus dune bonne dtermination des pmnres(rugosit, dbit, etc.) que des
mthodes numriques utilises Des aspects de lcoulement dbit variable et ses
applications sont galement abords dans cette annexe.
Le chapitre 4 sur les coulements brusquement varis donne diffrentes applications du
ressaut hydraulique, des dversoirs et des vannes aux domaines de comptence de 1Ecole.
Le chapitre 5 est une synthse des chapitLes 2, 3 et 4 puisquil fait le point sur les types de
lignes deau quon peut rencontrer autour des singularits dans les cas les plus courants. Des
applications rsultant dtudes exprimentales y sont galement abordes.
Enfin des notions sur les coulements non-permanents sont donnes au chapitre 6. Les
applications de ce chapitre sur des cas simples sont donnes et elles sont moins nombreuses
que dans le cas des coulements en charge.

A. L. MAR

Chapitre 1
Gnralits

CHAPITRE 1

GENERALITES
1 ELEMENTS
GEOMETRIQUES ET HYDRAULIQUES D'UN CANAL
1.1 PROFIL EN TRAVERS
1.2 PROFIL EN LONG
2 CLASSIFICATION DES ECOULEMENTS A SURFACE LIBRE
2.1 ECOULEMENTSCONSERVATIFS

2.2 ECOULEMENTNON CONSERVATIFS

3 VITESSES DANS UNE SECTION DE CANAL

3.1 VITESSE MOYENNE
3.2 REPARTITIONDES VITESSES ET MESURES
3.3 VITESSES LIMITES ET FORCES TRACTRICES
4 REGIMED'ECOULEMENT

4.1 NOMBRE DE REYNOLDS

4.2 NOMBRE DE FROUDE

5 PRESSION ET CHARGE MOYENNES DANS UNE SECTION
5.1 PRESSION
5.2 CHARGE

5.3 CHARGE SPECIFIQUE OU ENERGIE SPEClFlQUE

6 REVETEYENTS DES CANAUX
6.1 OBJECTIF DU REVETEMENT

6.2 STABlLlTE DU REVETEMENT

6.3 TYPES DE REVETEMENT ET MODALITE DE CONSTRUCTION.

fi

CHAPITRE 1
GENERALITES

Lcoulement de leau dans une canalisation peut faire soit en charge, soit surface libre. Ces
deux types sont semblables sur beaucoup daspects (les thormes gnraux sont les mmes)
mais la diffrence rside dans lexistence dune surface libre, cest dire une surface en contact
avec latmosphre.
Les coulements surface libre peuvent se rencontrer dans des canaux artificiels ou dans des
canaux naturels (cours deau, chenaux, etc.).
Les proprits des canaux naturels sont gnralement trs irrguiires. Une tude complte du
comportement des cours deau exige des connaissancess dans dautres disciplines telles que
lhydrologie, la gomorphologie, le transport de sdiment, etc. Elle constitue en elle-mme
une discipline appele hydraulique fluviale qui ne sera pas aborde dans ce cours.
Les canaux artificiels sont consruits ou mnags par lhomme pour les besoins divers
(navigation, hydro-lectricit, irrigation et drainage, assainissement pluvial et gout,
vacuateur de crue de barrage, etc.). Ils peuvent tre revetus ou non et leurs proprits
hydrauliques peuvent tre contrls ou apprhends de faon plus prcise dans leurs
conception. On les dsigne par diffrents termes techniques mais leur principe de
fonctionnement sont les mmes : canal (creus dans le sol, trs long, de faible pente
gnralement), aqueduc (suspendu gnralement pour traverser une dpression), chute et
coursier (gnralement avec une forte pente et trs court), gout (conduite non pleine, pour
lvacuation des eaux uss ou des eaux pluviales), foss (assainissement routier et agricole),
etc.

Aussi, les coulements surface libre prsentent plus de difficults que les coulements en
charge parce que les conditions dcoulement sont plus compliques :
La position de la surface libre peut changer avec le temps et lespace ;
Le dbit, la pente et la surface libre du canal sont interdpendants ;
Les rugosits des surfaces en jeu sont moins standardiss et elles varient avec la
profondeur deau ; do une plus grande incertitude quant aux valeurs adopter pour
les calculs ;
0
Les valeurs des coefficients dtermins exprimentalement dans les formules
universelles de perte de charge (Poiseuille, Prandtl-Von Carman, Blasius, Colebrook)
dpendent de la forme du canal.

1. ELEMENTS
GEOMETRIQUES ET HYDRAULIQUES D'UN CANAL
1.1. Profil en travers
Faisons une coupe perpendiculaire la direction d'un coulement surface libre (figure 1-1).
On dfinit les lments suivants :

Section mouille S : c'est l'aire occupe par l'eau et qui est dlimite par la surface
libre et les parois.
Primtre mouille P : c'est la longueur de la ligne de contact entre l'eau et les

parois .
Rayon hydrauliaue RH : c'est le quotient SAP. Pour une section circulaire de diamtre
D, R = D/4.
Diamtre hydrauliaueDH = RH
Hauteur d'eau ou tirant d'eau ou profondeur d'eau y : c'est la distance verticale entre la
surface libre et le fond du canal (le point le plus bas)
Largeur en gueule ou largeur en miroir ou largeur au plan &eau I : c'est la largeur de la
surface libre dans la section mouille.
profondeur moyenne y m : c'est le rapport Sn
largeur moyenne Im : c'est le rapport S4
Profondeur du centre de gravit JJG: c'est la distance verticale entre la surface libre et
le centre de gravit de la section mouille.

A-'P
.,

F i m e 1-1 :Profil en travers d'un coulement sufface libre et dfinitions

On dfinit galement quelques termes techniques avec le profil en travers d'un canal artificiel
(figure 1.2.).
la revanche r : c'est la distance verticale entre le plan d'eau et les berges. Elle peut
varier de 0,lO m pour les petits canaux 1,5 m pour les grands canaux. On peut
utiliser la formule de LACEY (quation 1-1) pour son dimensionnement :

o Q est le dbit du canal en m3/s et r la revanche en m.

le fruit du talus m : c'est le rapport entre la projection horizontale du talus et sa

projection verticale m = cot g @
Le fruit doit tre choisi en fonction de la nature des berges (pour des raisons de
stabilit de la talus) si le canal n'est pas revtu. Le tableau 1-1 donne la valeur de m
adopter pour les diffrentes natures des berges.
Pour les sections les plus courantes, les relations entre les diffrentes lments gomtriques
sont donnes au tableau 1-2.

Nature des berges

Roche dure, maonnerie ordinaire, bton banch

O 114

Roche fissure, maonnerie en pierres sches

1 12

Argile dure, latrite, terre compacte

314

Alluvions compactes

1/1

Bton non banch

5/4

Gros cailloux, enrochements

312

Dblais en terre ordinaire, sable, gravier

211

Terre remanie, remblais, sables fins

2,511 311

Tableau 1- 1 : Valeurs du finiit des berges m en fonction de la nature des parois

1.2 Profil en long

Pente du canal Z : c'est le rapport de la dnivele sur la distance horizontale ;la dnivele tant
compte positivement vers le bas (figure 1-3)

a) profll en dblais

,/

9-

radier ou plafond

b ) profll en remblais

F i m e 1-2 :profil en travers dun canal trapzique et dfinition des termes techniques

a) pente positive

b) pente negative ou contre pente

Figue 1-3 : Profil en long et dfinition de la pente dun canal

2 CLASSIFICATION DES ECOULEMENTSA SURFACE LIBRE

On peut rsumer la classification des coulements surface libre selon la variation de la

profondeur deau y en fonction de labscisse x le long de la conduite et du temps (cest
exactement les mmes classes et les mmes dfinitions quen coulement en charge).

11

ul

..

A
t

la

I
la

+
h>

tr

=lb

I
1

tr

Section

si k=1/2alors

1
-

Ym

k+ir

-T1
k+l

k+2

avec x=4y/T (1)

k

T=T*(+)

b+2(y-r)+7r r

b+2r

(x -4)r2

(T -4)r2+2(b+2r)y

y 2(b+2r)
b

---(i-mm)
T~ r2
4m m

---(l-rn@)
T r2
4m mT

~ ~ - 1-ma)
4 ~ ~ (
2(2mr-2&+T)
sin@

y 2 r et m = cotg 0
0 en radians

(1) : dans le cas o x l l , T + U donne une bonne approximation de P.

3T

13

2.1 Les coulements permanents

Cest un coulement o le dbit Q et la profondeur deau y ne sont pas fonctions du temps t.
On peut rencontrer dans ce cas les types suivants :
2.1.1.Ecoulements conservatifs
Le dbit nest pas fonction de x (il ny a pas de pertes ni dapports latraux).
2.1.1.1Les coulements uniformes

La profondeur deau y nest pas fonction de x

2.1.1.2 Les coulements varis
La profondeur deau y varie en fonction de x
2.1.1.2. ILes coulements graduellements varis

La variation de y en fonction de x est continue et graduelle. La fonction

y(x) est rgulire.
2.1.1.2.2Les coulements brusquement varis
La variation est brutale sur une courte distance. La fonction y(x) nest pas
rgulire.
2.1.2 Ecoulements non conservatifs.
Le dbit Q(x) varie avec labscisse le long de la conduite (canal avec fuites ou apports
latraux).
2.1.2.1 Ecoulement uniforme
Lcoulement ne peut pas tre uniforme si le dbit varie.
2.1.2.2 Les coulements varis
La profondeur deau y varie en fonction de x
2.7.2.2. 7 Les coulements graduellements varis
La variation de y en fonction de x est continue et graduelle. La fonction
y(x) est rgulire.
2.1.2.2.2 Les coulements brusquement v8rfds

La variation est brutale sur une courte distance. La fonction y(x) nest pas
rgulire.

2.2 Les coulements non permanents

2.2.1. Les coulements non permanents unbnnes
(rares sinon inexistants)

2.2.2. Les coulements non permanent varis

la profondeur d'eau y(t,x) dpend du temps t et de l'espace x
2.2.2.1. Les coulements non permanents graduellement varis
L a fonction y(t,x) est rgulire et les variations sont lentes et progressives.

2.2.2.2. Les coulements non permanents brusquement varis

La variation de y se fait sur une courte distance et un court intervalle de temps.
Des exemples d'coulements sont montrs figure 1-4.

3 VITESSES DANS UNE SECTION DE CANAL

3.1. Vitesse moyenne

On dfinit le dbit Q du canal par le volume d'eau traversant une section droite par
unit de temps. La vitesse moyenne U est le quotient du dbit par la section mouille.

3.2 Repartifion des vitesses et mesure

En fait, la vitesse n'est pas constante dans toute la section. Elle est nulle la paroi et maximale
au 113 environ de la profondeur. Des formules exprimentales donnent l'ordre de grandeur de
cette vitesse maximale dans la section.
Formule de Prony : V m a = U
0,82

15

~-

Ecouiements permanents
..

.-.
...-..--.--.

..
y
\.'.

Ecoulements non-permanents

-.
-.
-.

-.>., --

".. \.
2.

2.

-1 - \

Ecoulement uniforme :
Canal prismatique, loin des extrmits

Changemen
t
de
la
surface libre
le
avec

Ecoulement non permanent uniforme

Rare sinon inexistant

Propagation d'onde
graduellement vari

RV = rapidement vari
GV = graduellement vari
-._
...-.

Canal avec pertes latraies :l'coulement uniforme est

impossible

Figures 1- 4 : Exemples d'coulements surface libre

Raz de mare
rapidement vari

Dans une section, le lieu des points d'gaie vitesse est appel isodrome ou isotache. L'allure
gnrale des isodromes est donne la figure 1-5. pour diffrentes formes de section. La
rugosit et la prsence de coudes interviennent gaiement dans la distribution des vitesses.
Les vitesses se mesurent l'aide d'un moulinet dont l'axe doit tre plac dans la direction de
l'coulement. On obtient une quation de la forme.
V=an+b
o a et b sont des coefficients propres du moulinet et n la vitesse
de rotation obtenue avec un compte tour et un chronomtre s'ils
ne sont pas intrgus.

Avec les mesures de vitesse, on peut calculer le dbit par intgration des vitesses.

Canal

foss peu profond

Conduite circulaire

Canal rectangulaire
troit

Cours d'eau naturel

Figure 1-5 : allures des isodromes en fonction de la forme du canal.

3.3 Vitesses limites et forces tractrices

3.3.1 Vitesse minimale dans un canal
Pour viter les dpts des matriaux en suspension dans un canal, on adopte une vitesse
moyenne suprieure la vitesse minimale donne par la formule de Kennedy (1-2).

17

o y est le tirant deau et e, un coefficient dpendant des matriaux

transportes. Le tableau 1-3 donne les valeurs de e pour diffrents types
de matriaux en suspension.

Nature des matriaux en suspension

Limons et sables trs fins

0,40

Sables fins

0,50

Sables

0,63

Sables graviers

0,90

Tableau 1-3 : Valeurs du coefficient e de la formule de Kennedy

On peut viter les vitesses faibles en adoptant la forme du canal pour les faibles dbits (figure
1 - 6).

Figure 1.6 : Les faibles dbits sont vacus par le lit mineur.

3.3.2. Vitesse maximale dans un canal (Stabilit contre lrosion)

Pour un canal non revtu, il ne faut pas que la vitesse moyenne U dpasse une certaine valeur
sinon il risque dy avoir des affouillements du canal. Cette vitesse maximale est fonction de la
nature du terrain et de la profondeur dcoulement. Les tableaux suivants donnent les valeurs
de la vitesse maximale pour diffrents matriaux avec une profondeur dcoulement de 1 m et
les facteurs de correction pour des profondeurs diffrentes.

18

1) Profondeurs deau Y = 1 m. Canaux rectilignes

a) Matriaux non-cohrents
I

Matriau

Vitesse

Diamtre
mm

Vase.. ......................
Sable fin...................
Sable moyen.. ...........
Sable gros ................
Gravier fin... .............
Gravier moyen ..........
Gravier gros.. ...........

cailloux fuis ..............

0,15

0,005
0,05
0,25
1,O0
2,50
5,OO
10,oo
15,OO

15,o
25,O
40,O
75,O

..........
Gros cailloux. .............

0,20
0,30
0,55
0,65

1,20
1,40
1,80
2,40
2,70
3,50
3,90

100,o

150,O
200,o

0,so
1 ,O0
1,20

Trs peu
Peu
TrS
cornpact
compact
compact
compact
avec un indice avec un indice avec un indice avec un indice
de vide de 2,O de vide de 1,2 de vide de 0,6 de vide de 0,3
1,2
0,6
0,3
02

Matriau cohrent du lit

Argiles sableuses (pourcentage de sable infrieur 50 %).. ........................................

Sols avec beaucoup dargiles. ......................
Argiles. ................................................
Argiles trs fines.................................. r . .

0,45
0,40
0,35
0,32

0,85
0,80

1,30
1,25
120

0,70

1,O5

0 9

1,80
1,70
1,65
1,35

2) Facteur de correction pour des profondeurs deau y ze 1 m

Profondeur moyenne - m
Facteur de correction

0,3
0.8

0,5
0.9

1
1

0,75
0.95

1
1

1,o
1.0

1,5
1.1

2,o
-1 1

2,5
1.2

3,O
-1 3

Sinuosit

Rectiligne

Peu sinueux

Moyennement sinueux

Trs sinueux

Facteur de correction

1,O0

0,95

0,87

0,78

3.3.3 Forces tractrices

Une alternative pour le dimensionnement des canaux contre laffouillement cnsiste en
lapproche par les forces tractrices.

Dans un canal infiniment large, la force tractrice au fond est :

5

19

Dans un canai trapzodal, la distribution de la force tractrice a l'allure indique la figure

1-7a. La force maximum au fond est gale :

La force maximum dans les cts est de :

Sur le fond TM se produit au milieu ;sur les cts T'M est place une distance du fond d =
&*y. Les valeurs de KM ,K'M et & sont donnes dans les tableaux 1-3 et 1-4

Fipure 1-7a : Distribution des forces tractrices dans un canal trapzique

bly
O

1
2
3
4
6
8

2/1
KM
O
0,78
0,89
0,94
0,97
0,98
0,99

K'M
0,65
0,73
0,76
0,76
0,77
0,77
0,77

Kd

KM

033

O
0,78
0,89
0,94
0,97
0,98
0,99

092

092

092

312

K'M
0,565
0,695
0,735
0,743
0,750
0,755
0,760

Kd
0,3

0,2

0,2

0,3

Tableau 1-3 : valeurs de KM,K'Met & dans un canal trapzodai

\

KM
O
0,372
0,686
0,870
0,936

KM

O
0,468
0,686
0,870
0,936

Kd .

..

1,0
1,0
1,O
1,O

m
F M

Kd

211
0,650
093

Tableau 1-4 : valeurs de K'Met

3/2
0,565
093

1/1
0,480
095

213
0,375

097

112
0,325
097

dans un canaI triangulaire

La contrainte ainsi calcule ne doit pas dpasser la force tractrice critique T~ dont les valeurs
sont donnes ci-aprs pour des canaux rectilignes. Pour des canaux avec peu de courbes
(terrain faiblement accident) on doit prendre 0,90 des valeurs indiques ;pour une quantitd
moyenne de courbes (terrain moyennement accident) on prend 0,75 et pour des canaux avec
beaucoup de courbes (terrain trs accident) on, prend 0,60.
) Matriaux non cohrents gros

Au fond on prend comme valeur pour le projet :

o d75 est le diamtre auquel correspond, dans la courbe de

composition gradomtrique, 75 %, en poids, de matriaux de diamtre
infrieur.

Sur les cts on prend :

zf0 = K zo
o K est une fonction de l'angle de repos 8 du matriau et de l'angle des
cts avec l'horizontale $ (figure 1-7b)(quation 1-3).

F i m e 1-7b :Dtermination de la force tractrice critique pour les matriaux non cohrents gros

7.1

2) Matriaux non-cohrents fins :zo en kg/mz

1

Diamtre moyen d50 en mm

0,s

O, 1

0,2
Eau claire ..............................
0,o12 0,013
Eau avec peu de sdiments fins
0,024 0,025
1 Eau avec beaucoup de sdiments fins 10,038 0,038

l,o

2,o

5,0

0,015 0,020 0,029 0,068

0,027 0,029 0,039 0,081
0,041 0,044 0,054 0,090

3) Matriaux cohrents :T~ en kg/mz

Nature du lit

Trspeu
compact
1 avec un
indice de
vide de

1
\

Matriau cohrent du lit\

Argiles sableuses (pourcentage de sable
infrieur 50 % )
Sols avec beaucoup d'argiles
Argiles
Argiles trs fins

Peu
compact
1 avec un
indice de
vide de
2,O 1,2 1,2 0,6

0,020

0,077

0,015
0,012
0,010

0,069
0,061
0,047

Compact
avec un
indice de
vide de
0,6 0,3

compact
avec un
indice de
vide de

0,0160

4. REGIMESD'ECOULEMENT

Le rgime ou le comportement de I'coulement est fonction des :

- effets de la viscosit (frottement interne du liquide)

- effets de la gravit
compars aux effets des forces d'inertie.
4.1. nets de la viscosit :Nombre de Reynokfs'
Si le rapport Forces d'inertie/Forces de viscosit est trs grand, la viscosit n'me pas en jour.
L'coulement est turbulent.

Le nombre de Reynolds R

exprime ce rapprt. Si R est rs grand, ce qui est toujours

le cas en coulement surface libre, l'coulement est turbulent. Le nombre de Reynolds

n'intervient pas dans ce cas. Si R < 2000, l'coulement est laminaire, ce qui est rare dans le cas
des coulements surface libre.

4.2. Effets de la, gravit :Nombre de Froude

Leffet de la gravit est reprsent par le nombre de Froude F, dfini par la relation (1-4).

Le nombre de Froude donne lordre de grandeur des deux rapports :

0
Forces dinertie sur Forces de gravit
ou vitesse moyenne de lcoulement sur clrit des petites ondes.
Si F = 1, lcoulement est critique
Si F > 1, lcoulement est supercritique ou torrentiel.
Si F < 1, lcoulement est subcritique ou fluvial.

Le dnominateur c=,/gy m du nombre de Froude reprsente la clrit des petites ondes

gravitaires en eau peu profonde. Ces ondes peuvent rsulter dun changement momentan de
la profondeur deau (obstacle, perturbation, etc.). Les petites ondes peuvent remonter le
courant deau si lcoulement est fluvial et ne le peuvent pas si lcoulemnet est torrentiel.
Une petite perturbation dans un coulement peut donc tre un critre pour distinguer un
couiemnt fluvial dun coulement torrentiel.

5 PRESSION ET CHARGE MOYENNES DANS UNE SECTlON

5.1 Pression
La distribution des pressions est hydrostatique si lcoulement est quasi-rectiligne et quasiparallle (pas dacclration dans le plan de la section mouille). Cela suppose que le canal
soit prismatique ou assimilable un prisme. La pression en un point courant M (figure 1-8) est
donne par la relation suivante :

O q est 1angIe dfinissant la pente du canal

Au fond du canal, la pression vaudra pgy si lon prend les pressions relatives (pression
atmosphrique = O).
Si lcoulement nest pas quasi-rectiligne, leffet de la courbure est de produire une
acclration centrifuge dans la section mouille. La pression diminue alors dans le sens du
centre de courbure des lignes de courant. La pression au fond dun canal dont le profil en long
prsente une courbure est corrige par rapport la rpartition hydrostatique par les relations
(1-5) suivante.

o rM est le rayon de courbure suppos constant dans la section; VM la

vitesse au point M de la section et Q le dbit. Le rayon de courbure sera
positif si lcoulement est concave et ngatif si lcoulement est
convexe. Cette expression est valable sil ny a pas de dcollement de la
lame deau.

Figure 1-8 : Dfinition de la charge moyenne dans,unesection.

5.2 charge hydraulique
Elle est dfinie comme en hydraulique en charge (voir tome i) et son expression est donne
par la relation (1 -6) avec les notations de lhydraulyque surface libre.

o a est le coefficient dnergie cintique qui tient compte de la

rpartition de la vitesse dans la section (peu diffrent de 1 pour les
coulements turbulents rugueux) ;
a est un coefficient qui tient compte de la courbure du profil en long
(gal A 1 si lcoulement est parallle) ;
q est langle dfinissant la pente du c d 1 = tg q ;
zf est la cote du fond du canal ;
y est la profondeur deau.

La cote pizomtrique est dfhe aussi comme en hydraulique en charge et son expression
devient :

La cote de la sudace libre est donne quant elle par la relation suivante :

Sil ny a pas de courbure du profil en long (a=l) et si la pente du canal est faible
(I=tg(q)SlO%), alors la cote pizomtrique est confondue avec la surface libre. Cest ce cas
que nous tudierons pour la suite.
5.3 charge spcitique ou nergie spcifique
La charge spcifiqw dans la section dun canal est la charge hydraulique mesure par rapport
au fond du canal. Elle est reprsente par lquation (1-7) avec les hypothkses faites au
paragraphe 5.2 (cos q =1 ;azl et a=i)

La charge spcifique est une notion importante dans ltude des coulements surface libre en
particulier dans ltude des singularits.
Elle est bien spcifique au fluide (pression et nergie cintique) puisque mobilisable mme si
la pente du fond du canal est horizontale ou adverse, cest dire sans perte dnergie
potentielle globale (et mme gain avec une pente adverse). En effet on peut avoir
transformation dnergie cintique en nergie potentielle et inversement sans variation de la
charge spcifique.

6. REVETEYENT DES CANAUX

6.1. Objectifs du revetement

Le revtement d'un canal peut tre prvu pour diverses raisons dont les principales sont les
suivantes :

- impermabilisation (rduction des pertes par infiltration)

- rduction du fkit des berges (minimiser les excavations)
- adoption de vitesse leve (diminuer les risques d'rosion.

Si l'eau ne transporte pas de

sables et graviers, la notion de vitesse maximale peut tre oublie dans le design)
- rduction de l'entretien (l'herbe ne pousse pas sur certains revtements)
- augmentation du dbit (en diminuant la rugosit, etc.).

Une tude conomique est souvent ncessaire pour comparer le cot du revtement au gain
qu'il produit.
6.2. Stabilit du revtement

La stabilit du revtement est tudier en fonction de la nature du terrain et du type de

revtement.
Il faut tudier la stabilit des fondations (surtout pour les matriaux compressibles ou solubles
et la stabilit au profit (cas limite avec une vidange rapide). On prvoit des filtres si
ncessaires et les joints doivent tre tanches. Certains canaux d'assainissemnts sont munis de
barbacanes pour dissiper les pressions intersticielles qui ont tendances soulever le
revtement en cas de vidange.
6.3. Type de mvtement et modalit de construction
On rencontre le beton (le plus utilis) mais aussi la maonnerie, des briques, de l'asphate, de
l'argile, des feuilles de plastique etc.
De trs nombreuses techniques sont utilises pour la mise en place.

Pour le bton, on a des paisseurs de 10 20 cm environ. La mise en place du bton peut se

faire sans coffrage si m > 3/2. Autrement, il faut un cofhge. Le revtement des berges peut
tre diffrent de celui du fond du canal.

26

Chapitre 2
Ecoulement uniforme

CHAPITRE 2
ECOULEMENT UNIFORME

1 DEFINITIONET PROPRIETES

1.1 DEFTNITION
1.2 PROPRIETES
2 MISE EN EQUATION
2.1 FORMULE DE CHEZY
2.2

FORMULE DE MANNING-STRICKLER

3 CHOIX DU COEFFICIENT DE STRICKLER

3.1 CHOIX A PRIORI EN FONCTION DE LA NATURE DES PAROIS

3.2 CALCUL DE & A PARTIR DUN

JAUGEAGE

3.3 ESTIMATION DE I(S POUR LES MATERIAUX NON COHERENTS

3.4 ESTIMATION DE Ks POUR LES COURS DEAU NATURELS (FORMULE DE C o W m

ET

TABLE DE HORTON)

3.5

CALCUL DE & POUR LES SECTIONS COMPOSEES

4 CALCUL DE LECOULEMENT UNIFORME

4.1 CALCUL DE LA PROFONDEUR NORMALE OU DUN PARAMETRE DUN CANAL

4.2 DIMENSIONNEMENT DES CANAUX AVEC LA SECTION HYDRAULIQUEMENT FAVORABLE
OU SECTION << ECONOMlQUE bb
4.3 DIMENSIONNEMENT DES CANAUX AVEC UNE VITESSE LIMITE IMPOSEE

4.4 DIMENSIONNEMENT DES CANAUX AVEC UNE PENTE LlMITE

4.5 DIMENSIONS OPTIMALES A DONNER AUX SECTIONS DES CANAUX TRAPEZODAUX

DECOUVERTS

28

CHAPITRE 2
ECOULEMENT UNIFORME

1 D E ~ I T I OET
NPROPRIETES
1.1 Dfinitwn

Thoriquement, un coulement est dit uniforme si les filets liquides sont rectilignes, parallles
entre eux et parallles aux parois. Cela voudrait dire que la distribution des vitesses est la
mme dune section lautre.
Dans la pratique, on dira que lcoulement est uniforme si le dbit Q et la vitesse moyenne U
ne varient pas dune section lautre. Cette approximation est beaucoup plus simple et donne
des solutions satisfaisantes beaucoup de problmes pratiques.
Cependant, cette vitesse moyenne U ne doit pas tre trop leve pour entraner de lair et la
formations dondes dans leau et par consquent une instabilit de lcoulement (phnomne
de leau blanche). CRAYA a donn un critre dinstabilit dfini partir du nombre sans
dimension suivant :

o x est lexposant du rayon hydraulique dans la formule de la vitesse

U (2/3 pour la formule de Manning-Strickler et % pour la formule de
Chzy ci-aprs) ;
F, est le nombre de Froude calcul pour des pentes fortes

Fr=d-f--;
a gY CO=?
et y est un facteur de forme du canal dfini par y =l-RH
dP
ds

(y=l pour

une canal trs large et y0 pour une canal trs troit).

Si le nombre Vcr est suprieur 1, lcoulement sera stable car les ondes vont svanouir
rapidement. Dans le cas contraire, lcoulement devient instable.
Dans la cas de lcoulement uniforme, on sous entend galement quil est permanent car les
coulements non permanents uniformes sont rares sinon inexistants. Les forces de gravit et
les forces de fiottement sur la paroi squilibrent car il ny a pas dacclration du
mouvement.
La profondeur deau y correspondant lcoulement uniforme est appele profondeur
normale et note y,,.

29

1.2 Proprits

la vitesse est constante le long dun filet liquide (paralllisme des lignes de courant et
quation de continuit)
* la section mouille est constante (paralllisme des lignes de courant et des parois)
* la vitesse moyenne ne change pas dune section une autre (quation de continuit et
proprit prcdente)
* la rpartition des pressions est hydrostatique dans une section (paralllisme des lignes de
courant)
* les pentes de la surfaces libre, de la ligne dnergie et du fond du canal sont gales
(paralllisme des lignes de courant, des parois et vitesse constante).
On rencontre lcoulement uniforme ou on le suppose :
i ) Dans les canaux prismatiques, pente, section et rugosit constantes quand on se place

loin des extrmits.

2) Dans les tronons rectilignes des cours deau assimilables des canaux prismatiques.

h!T.ISE EN EQUATION

Lquation du mouvement uniforme a dj t tabli pour les coulements en charge

(paragraphe 5 du chapitre 2 du tome 1) en appliquant le thorme des quantits de mouvement
une tranche dx de lcoulement et a conduit lexpression suivante du frottement la paroi
zo en fonction du rayon hydraulique RH et de la perte de charge unitaire J :

Dans le cas de lcoulement uniforme surface libre, la perte de charge unitaire J est
gale B la pente du canal 1. Cest la seule diffrence avec les coulements en charge car la loi
de frottement est toujours donne dautre part par une relation exprimentale avec la vitesse
moyenne U du type suivant :

2.1 Formule de CHEZY

CHEZY est le premier poser que le frottement la paroi tait proportionnel au carr de la
vitesse moyenne :

z,= pg-U2
2-

C
La combinaison de cette relation avec lquation (2-1) conduit la formule de CHEZY
(quation 2-2).

30

u=C

J Z

o C est le coefficient de CHEZY. Il tient compte de la rugosit

du canai, de sa forme et des conditions d'coulement. C'est un
coefficient qui a la dimension d'une longueur la racine carre
par un temps [L" T'].
Tout ce qui suivra concernera donc le
Systme d'units Internationales (SI)

Il existe de nombreuses formules empiriques ou semi-empiriques qui drivent de la formule

de CHEZY en adoptant des expressions plus ou moins complexes de C. (BAZIN, KUTTER,
MANNING-STRICKLER, etc.).
Par exemple BAZIN a pos la formule (2-3) o y est fonction de la rugosit des parois et cette
formule a t longtemps utilise.

(2-3)

1+Y

RH
Dans la suite du cours, seule la formule de Manning sera utilise. C'est une formule monme
et qui est trs utilise cause de sa simplicit (Manning aux USA, Stricker en Europe).

2.2 Formule de ikt4"INGSTRICKLER

MANNING et STRICKLER ont une formulation simple du Coefficient de CHEZY qui
s'exprime dans la systme SI par les quations (2-4) qui ne diffrent que par la notation.
(2-4a)

c = K,

R ~ (STRICKLER)
%

O n est le coefficient de MANNING et K s = i

(2-4b)
est le

coefficient de STRICKLER.
Ces deux formules sont donc identiques et leurs substitutions dans la formule de CHEZY
donnent la vitesse moyenne U (quation 2-5) et le dbit Q (quation 2-6) en fonction de la
rugosit Ks,du rayon hydraulique RH (section S divise par primtre mouill P) et de la
pente du canal 1.

U'KSRH

b7

(2-5)

31

Le coefficient KS a la dimension L 113 T -1. Les valeurs qui sont donnes dans les tables
correspondent au Systme International (m 1/3/s), On trouve des valeurs de Ks allant de 100
(parois trs lisses) 15 (parois trs irrgulires) et il faut remarquer que plus la paroi est
rugueuse, plus K,est faible.

3 CHOIX DU COEFFICIENT DE STMCKLER

La dtermination de la rugosit est assez dlicate mais on dispose de certaines mthodes

destimation du coefficient de STRICKLER Ks.
3.1 Choix priori enfonction de la nature des parois

Il se fait en fonction de la nature des parois. Le tableau 2-2 donne la valeur adopter pour
diffrents types de matriaux.
3.2 Calcul de KS partir dunjaugeage

Si lon connat le dbit, la pente, la section mouille et le primtre mouille partir dun
jaugeage, on peut estimer & dans un bief. On fera ensuite des interpolations ou des
extrapolations.
3.3 Estimation de &pour les matdriaruc non cohrents

La rugosit est fonction du diamtre moyen des particules et elle est donne par les formules
suivantes dans le systme dunits SI :

-1

n = 0,041 djO6
1

n = 0,038 dw6

Le diamtre des grains d, est obtenu partir de la courbe grandomtrique. La notation

signifie que X % en poids des grains ont un diamtre infrieur dx.

32

Tableau 2-2 :Valeur du Ks de la formule de MANNING-STRICKLER

u = K, R~ 2/3 11/2
(d' aprs LENCASTRE)

KS

n=l/Ks

100 90

0,010 0,0111

85

0,0119

Parois lisses : planches avec des joints mal

soigns ;enduit ordinaire ;grs.

80

0,O 125

Bton lisse ; canaux en bton avec des

joints nombreux.

75

0,0134

Maonnerie ordinaire ; "cement-gun" ;

terre exceptionnellement rgulire.

70

0,O 142

Parois rugueuses ; terre irrgulire ; bton

rugueux ou vieux ;maonnerie vieille ou
mai soigne.

60

0,O 167

Parois trs rugueuses ; terre trs irrgulire

avec des herbes ; rivires rgulires en
lit rocheux.

50

0,0200

Terre en mauvais tat ; rivire en lit de

cailloux

40

0,0250

Terre complment l'abandon ;torrents

transportant de gros blocs.

20 15

0,0500 0,0667

Caractristiques

Parois trs lisses : revtements en mortier de

ciment et sable, trs lisses ; planches
rabotes ; tle mtallique sans soudures
saillantes
Mortier liss

33

3.4 Estimation de Kspour les cours d'eau naturels (Formule de COWAN

et table de

HORTON)
3.4.1 Formule de COWAN (empirique)

COWAN a dvelopp une procdure d'estimation du coefficient de MANNING n partir

d'une tude d'une cinquantaine de petits et moyens chenaux (RH5 5 m). Il a retenu
plusieurs facteurs affectant la rugosit n en proposant la formule (2-8)

n = (q,+ni

+ n2 + n3 + q ) m 5

(2-8)

ou
"0 = valeur de base pour un cour d'eau rectiligne uniforme et lisse et pour le matriau ,

envisag.
ni = traduit l'influence des irrgularits de surface
n2 = traduit les variations de forme et de dimension de la section du lit.
n3 = traduit l'influence des obstructions
r q = traduit celle de la vgtation et des conditions d'coulement

m5 = facteur correctif pour tenir compte des mandres

Les valeurs et critres de dtermination des coefficients ni, "2, "3, rq et n5 sont donns ciaprs et au tableau 2-3 d'aprs le CTGREF.

VALEUR de ni - La surface du canal est considre comme :

lisse

: pour des surfaces comparables aux meilleures qu'on

puisse obtenir avec le matriau envisag.

d'irrgularitmodre

: pour des canaux peu ou pas dragus ; pour des cours

d'eau naturels aux berges modrment rodes ou
affaisses.

d'irrgularit srieuse

: pour des cours d'eau naturels ou des canaux aux berges

trs marcageuses, rodes ou affaisses, pour des canaux

de contours et de surface irrguliers taills dans la pierre.

VALEUR de n2 - Variations de forme et de dimension de la section mouille.

- Variation progressive : pas de variation brusque

34

- Variation alterne occasionnelle :

0

les grandes et les petites sections alternent de temps autre ;

la variation de forme cre ventuellement un dplacement du flot principal d'un bord
l'autre.

- Variation alterne frquente :idem que ci-dessus mais plus fiquement.

VALEUR de n3

- Prsence et caractristique des obstructions : dpts d'ordures, racines

apparentes, souches d'arbres, blocs de pierre, troncs d'arbres, tombs ou
enracins, etc.

11 faut faire attention ne pas faire intervenir les facteurs dj considrs plus
haut. Pour juger de l'effet relatif des obstructions, il faut envisager :

jusqu' quel point les obstructions occupent ou rduisent la surface moyenne de la section;
la forme de ces obstructions (les objets anguleux aux formes pointues entranent une
turbulence plus grande que ceux qui ont des formes arrondies, des surfaces polies) ;
la position et l'espacement des obstructions, transversalement et longitudinalement dans le
bief considr.

VALEUR de n4 Influence de la vgtation

Faible pour des :
pousses denses d'herbes flexibles (gazon) dont la hauteur moyenne est de 1/2 1/3 celle
de l'eau (exemple : herbe bleue ou Bermuda) ;
0

jeunes plans d'arbres souples, comme le saule, le tamaris o la profondeur du flot est trois
quatre fois plus grande que la hauteur de la vgtation.
Modre pour des :

herbes o la profondeur d'eau est une deux fois la hauteur de l'herbe ;

herbes plus rsistantes ou jeunes plans d'arbres feuillage peu pais o la profondeur
d'eau est deux trois fois celle de la vgtation ;
buissons peu fournis, tels que les saules de un deux ans, en hiver le long des berges, sans
vgtation importante dans le lit du cours d'eau dont le rayon hydraulique est suprieur
70 cm.
Importante pour :

des herbes telles que la profondeur du flot soit gale la hauteur de la vgtation ;
pendant la saison d'hiver des arbres tels que saule ou coton de huit dix ans, entremls
de quelques buissons, sans feuillage, pour des rivires de rayon hydraulique > 70 cm ;
en priode vgtative des arbres tels que des saules touffus de un an, entremls d'herbes
pleinement dveloppement sur les berges, sans vgtation notable sur le lit de la rivire,
de rayon hydraulique > 70 cm.

35

Trs importante pour :

des herbes telles que la hauteur d'eau soit infrieure la moiti de la hauteur ;
en priode vgtative saules de un an touffus entremls d'herbes en pleine croissance sur
les berges, ou pousses denses de roseaux au fond du lit du cours d'eau dont le rayon
hydraulique n'est nulle part suprieur 3 ou 5 m ;
en priode vgtative des arbres entremls d'herbes et de buissons, tous avec beaucoup
de feuilles, le rayon hydraulique n'tant suprieur 3 ou 5 mtres.

VALEUR de m5
L'influence des mandres dpend du rapport de la longueur des mandres la
longueur vol d'oiseau du bief considr.
Cette influence est dite modre pour des rapports de l'ordre de 1 1,2,
apprciable pour des rapports de l'ordre de 1,2 1,5, importante pour des rapports
suprieurs 1,5.

Remarques
Cette mthode ne se proccupe pas de l'influence des particules en suspension et des
atterrissements. Les valeurs indiqus dans la table ci-aprs ont t tablies partir d'un
chantillon de 40 50 cours d'eau de faible importance. Nanmoins, cette mthode s'applique
des cours d'eau plus larges dont le rayon hydraulique dpasse, disons, 5 m.
Pour un cours d'eau naturel, K sera gnralement compris entre 20 et 40, les cas extrmes
pouvant aller de 10 pour une rivire trs obstrue (arbres en travers, etc.) 60 70 dans le cas
d'un canal btonn.

3.4.2 Table de HORTON

Dans le tableau 2-4, le CTGREF rsume les travaux de plusieurs auteurs qui ont estim les
valeurs du coefficient de MANNING adopter pour diffrents types de canaux et chenaux
naturels. Pour chaque type de revtement, les valeurs des coefficients de MANNING, de
STRICKLER et de rugosit de COLEBROOK sont donnes. Pour un canal qui sera bien
entretenu, on peut adopter les valeurs normales.
Le tableau 2-5 rsume les valeurs de Ks proposes par FRANCE-GABIONS pour les matelas
RENO et les gabions avec divers modes de mise en place.

3.5 Calcul de Ks pour les sections composes

La section mouille peut tre compose de N sous-sections avec une rugosit propre
chacune d'elles (figure 2-1). C'est le cas par exemple d'un cours d'eau en crue avec un lit
mineur et un champs d'inondation (lit majeur). La rugosit du lit majeure est en gnral plus
leve que celle du lit mineur. L'coulement dans le lit mineur agit sur l'coulement dans la

36

section des champs d'inondation en acclrant le courant. Inversement, l'coulement sur le

champ d'inondation a tendance freiner celui du lit mineur.
a peut aussi le cas d'un canal avec des parois latrales diffrentes de celles du fond.
Si la vitesse moyenne Uj dans les sous-sections sont diffrentes (cours d'eau en crue),
la formule de MANNING-STRICKLER peut tre applique sparment chaque
sous-section pour calculer son dbit. Le dbit total du cours d'eau sera gal la
somme des dbits de chaque sous sections (quation 2-9).

Si les vitesses moyennes Ui dans les sous-sections i sont identiques et gales la

vitesse moyenne dans la section, on peut calculer une rugosit quivalente de la
section et calculer le dbit total avec cette rugosit dans la formule de STRICKLER.
Sur la base de cette hypothse, EINSTEIN a tabli la formule (2-10) en crivant la
,
relation entre les sections ( S = ~ S>.

D'autres hypothses ont t galement faites pour le calcul d'une rugosit

quivalente. C'set ainsi qu'on peut poser que la force de rsistance totale sur le
primtre est gale la somme des forces de rsistance dans chaque sous-section et
que les rayons hydrauliques sont constants (quation 2-1 1).

(2-1 1)

37

F i m e 2-1 Section compose et calcul du dbit

VALEURS

CARACTERISTIQUES DU COURS D'EAU

Matriau

Terre
Roche
Gravier fin
Gravier grossier
Irrgularit du fond et des Lisse
parois
Irrgularit lgre
Irrgularit modre
irrgularit srieuse
Variation de la section
Progressive
Alterne occasionnelle
Alterne frquente
Influence des obstructions
Ngligeable
Faible
Apprciable
importance
Vgtation
Faible
Modre
Eleve
Trs leve
Influence des mandres
Modre
Apprciable
Importance

q,

0,020
0,025
0,024
0,028

0,000
0,005

nl-

10,010
0,020
0,000

n2

0,005
0,010
0,000

n3-

m5

0,010
0,020
0,040
0,005
0,010
0,025
0,050
1,100
1,150
1,300

0,015

- 0,015
- 0,030
- 0,060

- 0,010
- 0,025
- 0,050
- 0,0100

38

Tableau 2-4 : Valeurs des coefficients de rugosit en fonction de la nature de la paroi

(daprs CTGREF)

Ks

0,009 h 0,010

100 110

CANAUX ARTIFICIELS, GALERIES OU CONDUITES A

SURFACE LISSE

1. Surface trs lisses et sans saillies ;trac suffisamment droit

1t4

verre neuf et net

pyroline-cuivre
2. Surfxe lisses et sans saillies ;trac suffisamment droit
verre net et rabot

0,011

0,012

80 90

0,5

- mtal soud non peint

- ciment mortier ou bton bien liss, bien soign et sans dbris
- surfaces trs lisses avec courbures moyennes
3, Surface avec lgres asprits
- acier rivet ou peint
- fer forg ou coul

0,013 0,014

70 80

1,5 2

z 65

35

- bois non rabot

- ciment et mortier - finition la iruelie
- bton coffr avec de lacier ou du bois lisse sans dbris et pas de
courbures

- canaux en bton trs lisse avec joints

- tuyau de b i n a g e ordinaire

- gout vitrifi sans saillie

- brique verniss, grs

- asphalte lisse
- moellons dresss avec joints ciments
- surfaces lisses ou trs lisses avec fortes courbures

4. Surface avec asprits moyennes

- mtal incrust
- mtai rivet avec rivets grossiers
- canaux en mtal avec larges saillies vers lintrieur

- bois trs grossier (madriers)

- bton avec bord lisse et fond rugueux

- petit

canal en bton, assez droit et rgulier dont la surface es

recouverte dun lger dpt
- bois ou bton avec dveloppement dalgues et de mousses
gouts avec regards
- drains enterrs avec joint ouvert
terre particulirement rgulire
- canaux avec plafond en sable fin (surfaces non rides)
- surfaces lisses avec courbes excessives

0,015 0,016

Tableau 2-4 (suite) : Valeurs des coefficients de rugosit en fonction de la nature de la paroi
(d'aprs CTGREF)

5 . Surfaces rugueuses

- mtal trs incrust

- bton coul non liss
- bton coul aux cofhges en bois rugueux

Ks

k(mm)

0,O 17
0,O 18

55 60

10

150

15-20

0,022
0,023

A5

30 60

0,024
0,026

40

60 100

- bton trs rugueux ou vieux

- maonnerie vieille ou mal soigne
- canaux en maonnerie moyenne avec joints nombreux ou
nombreuses courbes

- bois ou bton avec dveloppement dense d'algues ou de

mousse

- canaux en terre &s rgulire, tat neuf, bon alignement

- sable moyen

- pierres dresses, joints ciments

6. Surfaces trs rugueuses

canaux en mtal avec trs fortes saillies vers l'intrieur ou
fortes courbures, ou dveloppement de vgtation importante
ou dbris accumuls
canaux en bton avec maonnerie en trs mauvais tat ou
trs grossire
bton coul non liss
canaux trs larges en gravier fin plus sable ou en terre
rgulire meuble, sans dveloppement de vgtation
- radiers pavs
moellons bruts assembls au ciment

0,O 19
0,021

7. Surface rugosit trs importante

lit en gravier fin

- canaux avec dpts ou vgtation

- canaux en terre moyenne, dimensions modres
- moellons bruts grossirement assembls au ciment
8. Surfaces assez grossires

- aqueducs mtalliques section semi-circulaire en tle plisse

- terre en mauvais tat

- graviers moyens
- canaux en terre, petites dimensions
- canaux en terre, plus larges, avec dveloppement de
vgtation ou gros galets ou pierres disperses

- fosss en bon tat

- canaux en terre sinueux, sans vgtation
- blocage ciment

- bton sur roche rgulirement excave

40

Tableau 2-4 (suite) : Valeurs des coefficients de rugosit en fonction de la nature de la paroi
(d'aprs CTGREF)

9. Surface grossires

- excavation rocheuse trs rgulire

- gros gravier
- pierre sche
- canaux en terre, dragus, sans vgtation ou enherbs

Ks

k (mm)

0,027 0,030

35

100 300

0,031 0,035

30

300 700

0,04

25

1O00

0,05

20

0,05
0,08

20
12

0,05

0,07

20
15

0,lO

10

- chenaux d'vacuation de crue, larges et entretenus

- bton sur roche irrgulirement excave

- canaux et fosss avec pierres rugueuses au fond et vgtation sur les
bords
1O. Surface trs grossires

- excavations rocheuses uniformes

- canaux avec dveloppement considrable de vgtation

- chenaux d'vacuation de crues, larges, mais peu entretenus

- blocage sec

- canaux en terre sinueux avec mauvaises herbes plus ou moins

denses ou plantes aquatiques
- canaux en terre sinueux avec fond en terre et berges en
blocage au fond pierreux et berges recouvertes de
mauvaises herbes
11. Surface excessivement grossires

- excavations rocheuses irrgulires

- canaux en terre en trs mauvais tat, trs sinueux avec
pierres rugueuses et vgtation importante

- lits majeurs d'vacuation de crue dgags, mais entretenus

de faon discontinue
12. Divers

- canaux non entretenus, mauvaises herbes et broussailles

coupes
- canaux en excavation avec broussailles
- canaux avec mauvaise herbes denses aussi hautes que la
hauteur de i'coulement
- fond net, broussailles sur les berges
- mme chose avec niveau d'coulement maximum sans
dbordement
- broussailles trs denses, niveau d'eau lev

41

Tableau 2-4 (suite) : Valeurs des coefficients de rugosit en fonction de la nature de la paroi
(d'aprs CTGREF)
n

Ks

0,025 0,03

30 4 0

0,035
0,040
0,0458 0,050
0,048
0,70
0,100
0,15 0,20

30
25
20

597

0,040
0,050

25
20

- pturages sous broussailles

0,030 8 0,035

30 35

- zones cultives, absence de rcoltes

- zones cultives, rcoltes sur pied
- broussailles disperses et mauvaises herbes ou broussailles

0,030
0,035 0,040
0,050

35
25 35
20

0,06 0,07

15

0,lO

10
25
16
10
8

B. RIVIERES NATURELLES
Pour les cours d'eau section suffisamment constante on se
reportera au tableau A.
1. Petit cours d'eau (largeur maximale infrieur 30 m)

a) cours d'eau de plaine

- net, droit, niveau d'eau lev, peu de variation de la section
mouille
- idem, mais pierres et mauvaises herbes plus nombreuses
- net, sinueux avec seuils et mouilles
- idem, mais avec pierres et mauvaises herbes
- idem, mais niveau bas
- cours paresseux, mauvaises herbes, trous d'eau profonds
- nombreuses mauvaises herbes et nombreux trous d'eau
- pentes et fond irrgulier, nombreuses souches, arbres et
buissons, arbres tombs dans la rivire

b) cours d'eau de montagne

Pas de vgtation dans lit, rives escarpes, arbres et
broussailles pour les niveaux levs
- fond en gravier et cailloux, peu de gros galets
- fond avec gros graviers

20
15
10

2. Plaines d'inondation

et quelques arbres en hiver

quelques arbres et broussailles en t ;broussaille moyenne
ou dense en hiver
- broussaille moyenne ou dense en t
- souches d'arbres sans rejet
- souches d'arbres avec rejets durs
- fort de hautes futaies ;peu de broussailles
- idem, avec niveau d'eau atteignant les branches
- souches denses

0,040

0,060
0,lO
0,12
0,15

3. Grands cours d'eau (largeur maximale suprieure 30 mtres)

La valeur de Ksest supneure celle des petits cours d'eau

d'allure analogue, car les rives offtent moins de rsistance
efficace.
- section rgulire sans broussailles

- section irrgulire et rugueuse

D,025 0,040

25 40

3,040 O, 100

1025

42

Tableau 2-5 :Valeurs de Ks pour les gabions (daprs FRANCE GAF3IONS S.A.)
Revtement

I Ke

ar coule directe
3
4
5

7
8
9

10

C m u x revtus en gabions remplis avec des matriaux bien slectionns

et mise en oeuvre soigne
Canaux revtus en gabions remplis avec des matriaux non slectionns
et mise en oeuvre non soigne
Canaux en terre avec entretien nglig.

38
35
33

Remarque : MR veut dire matelas RENO

4 CALCUL DE LECOULEMENT UNIFORME

4.1 Calcul de la profondeur normale ou dun paramtre dun canal

La formule de MANNING-STRICKLER , associe la dfinition de la vitesse moyenne U,,

sera choisie pour la suite des calculs.

US

Pour lcoulement uniforme, la perte de c h g e unitaire J est gale la pente du canal 1.

Comme J est toujours positif, lcoulement uniforme ne peut avoir lieu que si la pente du
canal 1 est strictement positive. Autrement (canai fond plat ou canal contre pente),,il peut
y avoir un coulement vari mais pas un coulement uniforme.
Les paramtres calculer avec ces quations sont de 4 types :
paramtres gomtriques : rayon hydraulique RH@) et surface mouille S@) qui sont
des fonctions de la profondeur deau y (voir tableau 1-1) ;
paramtres hydrauliques : vitesse moyenne U et dbit Q ;

43

paramtre de rugosit : lcoefficient de STRICKLER Ks ;

paramtre de pose : pente du canal 1.
Le seul paramtre pour lequel lquation de MANNING-STRICKLER nest pas explicite
est la profondeur deau y (ou la largeur au plafond b dans le cas dun canal trapzodal).
Les calculs consisteront donc dterminer un des paramtres connaissant les autres mais
seul le calcul de la profondeur normale yn (ou la largeur au plafond b dans le cas dun
canal trapzodal) peut poser des difficults.
Pour le calcul de la profondeur normale, on peut utiliser une des mthodes ci-aprs :
Une mthode graphique consistant par exemple tracer sur papier millimtr la
-

dbitance du canal SR,,,^ en fonction de y point par point. On porte en ordonne la

Q et on lit la valeur de yn recherche en abscisse (figure 2-2).

valeur connue de -

fi

Des abaques du type des figures 2-3 et 2-4 par exemple pour les canaux rapzodaux
et circulaires. Il existe plusieurs abaques de ce type pour les sections paramtres.
0

Des mthodes numriques de rsolution de fonction implicite en y (mthode de

NEWTON, mthode du point fixe, mthode des intervalles, etc.) ou tout simplement
utiliser le solveur dEXCEL ou de la TI92 (voir paragraphe 6 du chapitre 3 du tome 1
sur les coulements en charge).

Le tableau 2-6 rsume les types de calcul avec la formule de MANNING-STRICKLER et

pour les sections courantes en tenant compte des paramtres gomtriques du canal vus au
tableau 1-2
Dans le cas de la section circulaire, on peut utiliser, dans quelques cas particuliers de taux de
remplissage y/D, les formules approches du tableau 2-7.

Exercice i : On veut couler un dbit Q = 10 m3/s dans un canal trapzodal creus dans la
roche (Ks = 40) de pente, 1= 0,001, dont les berges auront pour fivit m = 114.
1. Quelle sera la profondeur normale si on prend b = 5 m ?
2. Quel serait le dbit si on mettait un revtement en bton (Ks = 70) ?
Rponse
1. y = 1,48m.
2. Q = 17.5 m3/s
Exercice 2 : Une conduite circulaire de diamtre D = 1 m est pose avec une pente 1= 0,l 960.
Le coefficient de Strickler est de 70.
1. Pour un dbit Q = 128 l/s, quelle sera la profondeur dcoulement ?
2. Quel sera le dbit si le taux de remplissage est de 80 % ?
3. Avec ce dbit, quelle pente faudrait-il la poser pour avoir un pourcentage de remplissage
de 75 % ?

44

Rponse
1. yn= 0,55 m
2. Q =0,21OYs
3. 1 =0,11 %O

Figure 2-2 : Mthode graphique de calcul de la profondeur normale.

45

Valeurs de y& ou yiJD

..

!3

7D.
N
E:

l-Section

Calcul de Q
direct

Calcul de 1
direct

Calcul de yS
direct

Calcul de y,,
Solveur ou abaques ou

Calcul de b
Solveur ou abaques
ou itrations

Calcul de D
direct

1 8 en radians

l m

.c--c

bK,l/jyf

Q'

10

b2KiYT

Q
b%FyS

Canai infiniment large

-io,1
Y

Tableau 2-6 : )mulaire du calcul de coulement de l'coulement uniforme pour les sections les plus courantes.

48

Taux de
remplissage
Y
D
1

Angle
mouill

Formule approche du dbit

360"

Qp=O,312 K s & 0:

0'95

308"

Qma = 0,327 Ks & D:

0,83

262'6"

Qp=O,312 Ks& Dt

0'75

240"

0'50

180"

QO,SOP
= 0,151 Ks .Ji 0
9

0,45

168,5"

Qo,45~= O , m Ks JI D!

Q o , ~=~0,284
P
Ks JI Df

Tableau 2-7 : Formule approche d'un canal circulaire pour des cas particuliers de
remplissage.

4.2 Dimensionnement des canaux avec la section hydrauliquementfmrable ou section

u conomique N

4.2.1 Dfinition et proprits de la section hydrauliquement favorable

Pour qu'une section soit conomique il faut que simultanment les grandeurs qui
interviennent dans le cot soient minimales. Ces grandeurs sont la surface S (terrassement)
et le primtre P (revtement).
Il en rsulte que pour un dbit et une pente donne, la section ((conomique N ou
hydrauliquement favorable correspond la vitesse la plus grande.
Donc la section ainsi dfinie est la fois celle qui :
0
Donne le dbit le plus grand pour une surface S et une pente 1 donnes (primtre
mouill minimum) ;
0
Ncessite la section S et le primtre mouill P les plus petits pour un dbit Q et une
pente 1 donns (cot minimum) ;
0
Ncessite la pente 1 la plus faible pour une section S et un dbit Q donns (primtre
minimum).

La section qui sera ainsi calcule est la section hydrauliquement favorable (HF) mais en fait
pas forcment la plus conomique, parce que :
0
dans certains cas le cot n'est pas directement li la section S et au primtre mouill
p;
le raisonnement ne tient pas compte de la revanche.
En outre il existe trs souvent des contraintes de construction qui limitent la profondeur des
canaux et empchent d'adopter la section hydrauliquement favorable (variation des cots

49

dexcavation avec la profondeur due par exemple une htrognit du sous-sol telle que la
prsente de nappe phratique ou de roche dure.
Pour une forme de section donne (rectangulaire, trapzodale de f i t m donn), la section
HF est caractrise par une valeur du rapport entre la profondeur y et la largeur du radier b qui
est constante. Les sections HF dune forme donne sont donc semblables.
La forme circulaire est celle qui assure le plus petit primtre pour une section donne dans
un plan. Il en est de mme de la forme semi circulaire qui a le mme rayon hydraulique dans
le demi-plan. Les sections circulaires ou semi circulaires ne sont cependant ralisables quen
matriaux rsistants (galerie en rocher, conduite mtallique, canalisations et canaux suspendus
en bton).

Dans les autres cas, on adopte plutt des sections trapzodales isocles pour faciliter la
construction et lentretien des canaux.
4.2.2 Calcul de la section trapzodale HF

Minimiser la section S et le primtre mouill P, qui sont tous les deux fonctions de la
profondeur deau y et de la largeur du radier b, revient annuler les diffrentielles totales
exactes (y et b variant) et trouver la solution en (dy,db) diffrente de la solution triviale (0,O).
Lannulation du dterminant du systme linaire en dy et db ainsi form est une condition
ncessaire et suffisante qui conduit la relation (2-12) suivante entre la profondeur y et la
largeur au radier b pour une section trapzodale HF.
(2-12)
Pour la suite des calculs, on peut poser n = & b q - m .

La substitution de la relation (2-12) dans les expressions de la section, du primtre mouill et

du rayon hydraulique en fonction de y et b, conduit :

S=Ay2

On constatera que le rayon hydraulique RH est indpendant du fruit de la berge m pour une
section trapzodale HF.
Pour une section trapzodale isocle la formule de Manning-Strickler devient alors explicite
pour toutes les variables (quations 2-13)

50

(2-13)

Plus gnralement, si l'on connat les fruits de la berge m et de la diagonale n d'un trapze
isocle quelconque, on a les relations suivantes :
b=(n-m)y
1=(n+m)y

P= (n+h)y
Pour n et m fixs, les trapzes isocles correspondants s'obtiennent donc par homothtie de
centre B en faisant varier y (figure 2-5a). On dmontre que n est gal ih 3, pour une section
trapzodale isocle HF ( minimum de P et S pour n et y variables).
La section trapzodale HF est circonscrite un demi cercle de diamtre horizontal
concidant avec la surface libre et dont le rayon est la profondeur d'eau y (figure 2-5b).
En effet on peut tablir , partir de OB=AB/2=1/2=dgy=BC , que le triangle OBC est
isocle et que sa hauteur OF=y=OG.

Bercice 3 : Soit couler un dbit Q = 10 m3/s dans un canal de section trapzodale, creus
dans la roche, de pente l%, dont les berges auront pour fruit m = 1/4. On prendra Ks = 40.
1O ) Calculer y et b pour un volume de dblais et un primtre mouill minimaux.
2") Quelle sera alors la vitesse moyenne dans le canal ?

Rponse
1") y=2,10 m et b= 3,20 m
2") U=1,30 m/s
4.2.3 Comparaison de sections trapzodales isocles HF

1 O) Pour un dbit une rugosit Kset une pente I donns, on peut comparer les paramtres
y, S, P et U des sections trapzodales isocles hydrauliquement favorable en fonction de m.
On prendra comme rfrence la section rectangulaire o m
indices avec O).

O (les variables sont alors

51

a) section trapzodale Isocle quelconque

I

b) section trapzodale isocle HF

F i m e 2-5 : Section trapzodale isocle HF

On peut alors tracer les courbes

montres la figure 2-6,

01,

p, y et p dfinies ci-aprs en fonction de m et qui sont

52

On peut remarquer sur la figure 2.6 que :

0

les variations de a,p, y et p sont faibles pour des valeurs de m infrieures A 413 ;
a et p passent par un maximum ;y et p passent par un minimum pour la valeur de m
gale J3 correspondant 8 = 60" ou an demi-hexagone.
3

pour 8 = 36'56' ou m = 4/3, la section trapzodale HF est exactement quivalente au

re-gle
I3.F YO, so, Po
A toute fin pratique, on peut donc utiliser l'approximation de la section rectangulaire HF,
pour m compris entre O et 4/3 (et mme 5/3) o les variations sur les paramtres a,p, y et p
sont infrieures 5 %, afin de dterminer la pointure d'un canal pour les niits des berges
admissibles.

27 Pour y, K, et Z donns,

les valeurs de Q, P, S et U d'une section trapzodale HF de fruit

des berges m, rapportes respectivement celles d'une section rectangulaire HF (m=O)
donnent les paramtres ci-aprs.

p ' =S--a
so

La fonction h=2a' est trace sur la mme figure 2-6 et dans les mmes limites (m compris
entre O et 5/3), les variations des paramtres sont infrieures 13 %. Elle passe par un

J? (demi-hexagone) et m=4/3, elle donne la mme valeurs que m=O.

minimum m=3
4.2.4 Considdrationspratiques

Le trapze de section mouille minimum ne constitue pas toujours la solution la plus

conomique et la meilleure.
D'abord, elle ne donne, coup sr, un volume de dblais minimum que si la sudace du terrain
o il est creus est horizontale ou sub-horizontale (annexe du chapitre 2).

54

Mme si le terrain est horizontal, on doit tenir compte des lments suivants :

I y la profondeur doit tre limite

le prix au mtre cube de dblais, crot gnralement avec la profondeur. C'est le
cas quand les couches de terrain deviennent plus dures en profondeur.
- les infiltrations dans un canal non revtu croissent avec la pression de l'eau sur le
fond, donc avec la profondeur.
- un canal moins profond est plus facile entretenir.
-

D'o il pourrait tre prfrable d'avoir un canai un peu moins profond, quitte avoir un
volume total de dblais plus grand.

20)la vitesse doit tre limite

La vitesse d'coulement doit tre infrieure une limite Umaxpour viter l'rosion. C'est ce
que nous verrons au paragraphe suivant
4.3 Dimenswnnement des canaux avec une vitesse limite impose

Si l'on dsire respecter une vitesse limite, on fait le calcul suivant pour cette vitesse qui
introduit une relation entre la surface S et le rayon hydraulique RH compltant la formule de
Manning-StricMer pour la dtermination de la profondeur d'eau y et la largeur au radier b.
Etant donne Q, 1 et m on veut dterminer b et y en ne dpassant pas la vitesse d'rosion
Umm.

Dans un premier essai, on prendra U = Umm.

La dfinition du dbit permet de tirer S ;

L'quation de Manning-Strickler permet de tirer RH :

D'o l'on obtiendra le systme de deux quations suivarus en y e b :

y(b+my)=S

55

Par la mthode de substitution, ce systme est quivalent :

b-S-2yG

(2-14a)

A y 2-y+s=o
s

(2-14b)

RH

RH

L'qu tion du second degr (2-14b) admettra des ra ines selon 1- valeur du discriminant A :

ou

A : 1--

avec Um dsignant la vitesse de la section

hydrauliquement favorable et le signe I la proportionnalit.

Si A est ngatif, c'est dire si U est suprieure UW,alors il n'y a pas de solution.
Si A est gale zro ou U=UHF, alors il y'a une solution qui est la section
hydrauliquement favorable :

1 s
Y"==%

Si A est positif ou U infrieure Um, on obtient 2 solutions positives en y qui ne sont

pas toujours des solutions physiques du problme :

-s+&
RH

Yi= 22

On doit en effet vrifier que bi et

positifs.

(correspondant respectivement y1 et y2) restent

56

Or b est positif si la condition ci-dessous est vrifie:

y2 vrifiera toujours cette condition mais il n'en est pas de mme pour yl.

Si y1 vrifie cette condition, on aura deux solutions physiques. On devrait retenir la

solution y2, b2 qui prsente la plus faible profondeur. C'est pourquoi que la mthode
s'applique plus souvent l'tude des coulements larges(endiguement de rivire par
exemple)
Dans tous les cas, la solution retenue devra tre pratiquement ralisable
(dimensions pas ((microscopiques)) ni gigantesques ).
Si on ne trouve pas de solutions (A<O) ,on doit diminuer U et recommencer l'essai avec cette
nouvelle valeur.
De toutes les faons, les variations de A sont lentes de part et d'autre de la valeur zro. Si l'on
prend U = 0,9UHF par exemple, la profondeur y1 est gale 0,54y0 et y2 gale 2,06 yo.
On ne doit donc pas s'attendre une rduction substantielle de la vitesse U en jouant sur les
dimensions y et b quand Ks, m, d, 1 et Q sont fixs C'est plutt la pente qu'il faut modifier
pour diminuer la vitesse de faon substantielle.
4.4 Dimensionnement des canaux trapzodaux avec une pente limite

Lorsque l'on dsire respecter une vitesse limite U, on peut tre conduit rechercher une pente
limite soit pour viter les dpts soit pour viter les rosions. Les pentes ainsi calcules
peuvent conduire soit prvoir des chutes soit prvoir des exhaures des exhaures.
4.4.1 cas d'une section hydrauliquement Favorable

Si par exemple on cherche la pente maximale n'entrainant pas l'rosion (U=U,,) et si l'on
choisit une section HF, on introduit 2 relations entre les variables qui compltent la formule
de Manning-Strickler ; ce qui permet de calculer simultanment la profondeur d'eau y, la
largeur au radier b et la pente 1 connaissant le dbit Q, la vitesse U, la rugosit Ks et le fruit
des berges m.

D'aprs la dfinition du dbit et la section HF :

Q=US=A y* U

La formule de Manning-Strickler s'crit alors :

57

La premire relation permet de tirer y :

En fin la section hydrauliquement favorable permet de calculer b :

b = y (A - m)

En substituant, la valeur de y dans la formule de Manning-Strickler, on peut tirer la pente 1 :

8

Mais on avait dj annonc que la section hydrauliquement favorable pourrait conduire des
profondeurs grandes. On peut donc imposer une largeur b plus grand que la valeur donne par
la calcul qui vient dtre effectu et chercher la pente limite.

Exercice 4 : Calculer la pente 1, la largeur au radier b et la profondeur deau y dun canal

trapzodal de section hydrauliquement favorable avec un fhit des berges m=l pour un dbit
Q=1,3 m3/s, une vitesse maximale de U=0.65 m/s, et une rugosit Ks = 100.
Rponse
y = 1,05 m
b = 0,87 m
I = 10-4
4.4.2 cas dune valeur de b ou y impose

La dfinition du dbit et de la section donne :

S =

9
U

y(b+my)

On aura donc une quation du 2 degr en y qui donne une seule solution positive ou une
quation linaire en b :

V=
J

2m

si m ; t O et y = Q sinon
bU

58

La relation de Manning permet de tirer 1 :

D'o 1en fonction de b ou y:

1013
et I=X[~+%]
K i Q413

413
sinon

Exercice 5 : Calculer la pente 1 et la profondeur d'eau y d'un canal trapzodal de largeur au

radier b=l m avec un h i t des berges m=l pour un dbit Q=1,3 m3/s, une vitesse maximale
de W0.65 d s , et une rugosit K, = 100.
Rponse
y = 1,00 m
I = 10-4
4.5 Dimensions optimales donner aulc sections des canaux trapzodaux dcouverts

Le canal ayant porter un certain dbit Q, il s'agit de raliser une section rpondant aux
mieux aux conditions suivantes, parfois contradictoires
1. largeur des emprises rduite le plus possible (conomie sur les achats de terrain) ;

2. Profondeur la moins grande possible (entretien plus facile, infiltrations moins

importantes, conomie sur le prix de revient unitaire des dblais et le cot total des
travaux d'excavation) ;
3. Section la plus rduite possible (pour diminuer le volume des dblais et le cot total
des travaux) ;

59

5. Section teIle que la vitesse U =

ne soit ni trop faible (pour viter les dpts) ni trop

forte (pour dviter l'rosion).

Dans un projet, on cherche gnexaiement si satisfaire tout d'abord les conditions 3 et 4 ;puis
on modifie les dimensions obtenues, s'il y a lieu, pour donner suite aux deux premires
conditions (largeur et surtout profondeur). La cinquime condition peut recevoir satisfaction
en modifiant la pente 1plutt que les dimensions transversales du canal.

Chapitre 3
Ecoulements graduellement varis

CHAPITRE 3
ECOULEMENTS GRADUELLEMENT VARIES

1 DEFINITION D'UN ECOULEMENT GRADUELLEMENT VARIE

1.1 ECOULEMENTS
VARiES
1.2 ECOULEMENTS
GRADUELLEMENT VARIES
2 VARIATIONS DE LA CHARGE SPECIFIQUE Hs
2.1 DEFINITION
2.2 VARIATION DE Hs EN FONCTION DE LA PROFONDEUR Y POUR UN DEBlT Q DONNE
2.3 VARIATION DE Q EN FONCTION DE Y POUR UNE CHARGE SPECIFIQUE Hs DONNEE.

3 REGIMECRITIQUE
3.1 DEFINITION
3.2 CALCUL DE LA PROFONDEUR CRITIQUE Yc
3.3 APPLICATIONS AUX SECTIONS COURANTES
4 PENTES ET REGIME
4.1 PENTE CRITIQUE
4.2 EXPRESSIONDE LA PENTE CRITIQUE
4.3 TYPES DE CANAUX (CLASSIFICATION SELON LA PENTE).

5 COURBES DE REMOUS
5.1 ETUDEQUALITATIVE DES COURBES DE REMOUS
5.2 METHODE
DE CALCUL DES COURBES DE REMOUS

62

CHAPITRE 3
ECOULEMENTS GRADUELLEMENT VARIES

! DEFINITION DUN ECOULEMENT GRADUELLEMENT VARIE

1.1 Ecoulements varis

Les coulements varis se rencontrent dans les rivires irrgulires, dans les singularits des
canaux et dans les zones de transition entre deux coulements uniformes.
Ces coulements peuvent tre permanents (en un point ou section donn, rien ne change dans le
temps). Mais dune section lautre, la profondeur, les pentes de la surface libre et de la ligne
dnergie, la section mouille varient.
Selon lallure du changement on distinguera :
- les coulements graduellement varis
- les coulements brusquement varis.

1.2 Ecoulements graduellement varis

On a une variation lente et continue de la section mouille, des pentes de la ligne deau et de la
ligne dnergie, ainsi que de la profondeur deau dune section lautre.
On admet les hypothses suivantes dans ltude de ces coulements :
0

La courbure des lignes de courant est ngligeable

Ces lignes sont parallles et normales la section ;do la rpartition des pressions est
hydrostatique dans la section

Le coefficient dnergie cintique a est constante

La pente de la ligne dnergie J est donne par la formule de Manning :

avec x labscisse le long du canal,

mais J diffrent de la pente 1du canai.

63

2 VARIATIONS DE LA CHARGE SPEClFlQUE Hs

2.1 Dfinitions
La charge moyenne H dans la section a dj t dfinie au chapitre 1 et elle est donne par la
Relation suivante :

= z

+y +

U2
a2Y

La charge spcifique est la charge moyenne ramene la base de la section et elle peut sexprimer
par lune des expressions suivantes o lon a pris a=l :

La charge H est toujours dcroissante dans le sens de lcoulement ( quand x crot). Il nen est pas
de mme pour la charge spcifique Hs car le signe de dHs/dx dpend de la diffrence de la pente
du canal I et de la pente de la ligne dnergie J. En effet on peut calculer successivement :

dH- -

dx

dHs

dx

+ -dz

dx

dx
dH,=I-J
dx

(3-2)

On remarquera quen rgime uniforme, la pente de la ligne dnergie J qui caractrise la perte
dnergie par frottement est gale la pente du fond du canal 1 et la pente de la ligne deau par
rapport au fond dy/dx est nulle. Do la charge spcifique Hs est constante le long du canal
comme les autres variables hydrauliques et gomtriques.
Avant datteindre le mouvement uniforme, la variation dnergie potentielle du fluide nest pas en
gnralgaie lnergie dissipe par frottement :

Si 1est suprieur J, lnergie cintique crot donc la vitesse V, do la profondeur deau y

diminue.
Si 1 est infrieur J, lnergie cintique diminue donc la vitesse V, do la profondeur
deau y augmente.

64

2.2 Variation de Hs en fonction de la pmfondeur y pour un dbit Q donn

Lquation (3-1) dfinit, pour une section dtermine, une relation entre la charge spcifique Hs,
la profondeur deau y et le dbit Q. En fixant Q=constante, on obtient une courbe Hs(y) qui a
lallure indique dans les figures 3-1 pour y variant entre O et les grandes profondeurs et pour une
section vase vers le haut. Cette courbe est trs commode pour tudier les problmes de
transition entre 2 biefs de canaux prismatique (voir chapitre 5).
Lune des branches est asymptotique laxe des ordonnes et lautre la premire bissectrice
(droite 45). En effet Hs(y) tend vers linfini comme @/2gS2 et y ;quand y tend respectivement
vers O et linfini.
La courbe admet un minimum pour une valeur de la profondeur appele profondeur critique yc
et qui est la solution de lquation (3-3) qui exprime que le nombre de FROUDE est gal 91.

En effet, le minimum est obtenu avec lannulation de drive de Hs par rapport y (quation 3-4)
obtenue en faisant remarquer que la drive de S par rapport y est gaie la largeur en gueule 1.

On peut galement constater sur la courbe de la figure 3-1b :

Avec une charge spcifique donne Hs suprieure la charge spcifique critique Hsc dune
section dtermine, le mme dbit peut scouler sous 2 profondeurs diffrentes :

- lune y correspondant au rgime torrentiel ou supercritique et qui est infrieure

Yc ;

lautre y correspondant au rgime fluvial ou subcritique et qui est suprieure yc.

Ces 2 profondeurs sont dites conjugues avec la charge spcifique.

La charge spcifique critique Hsc est la charge minimale permettant lcoulement du dbit
Q dans la section considre. Lcoulement du dbit Q a lieu en ce moment avec la
profondeur critique. Ce minimum Hsc est dautant plus grand que le dbit est lui-mme
plus grand
Au voisinage du rgime critique Hsc, une lgre variation de Hs provoque des variations
importantes des profondeurs deau y ; cest pourquoi le rgime critique est instable
(ondulation apprciable de la surface libre due de petites irrgularits).

65

--f

torrentiel

_I_____

Y"

-F-e 3-1 : Variation de la charge spcifique Hs en fonction de la profondeur d'eau y dans une
section dtermine.

66

En rgime fluvial, la charge spcifique HS et la profondeur deau y varient dans le mme

sens (simultanment croissantes ou dcroissantes) ; lnergie cintique Q2/2gS2est plus
faible que lnergie potentielle y.
0

En rgime torrentiel, cest le contraire qui se produit : Hs et y varient dans le sens contraire
et lnergie cintique est plus grande que lnergie potentiel y.

2.3 Variation de Q en fonction de y pour une charge spcitque Us donne.

Ltude de la variation du dbit Q en fonction de la profondeur deau y pour une charge
spcifique donne a galement un intrt pratique dans les transitions lac-canal ( voir chapitre 5).
Si lon fixe Hs dans la formule (3-1), on obtient la relation (3-5) entre le dbit Q et la profondeur
deau y dont une illustration est donne la figure 3-2.
(3-5)
Le dbit est nul pour y=O et y=Hs et on obtient sa valeur maximale (fonction continue variant de O
O pour y variant de O Hs) en annulant la drive dQ/dy donne par lquation suivante :

Le maximum de dbit appel dbit critique Qc respecte donc la mme condition (quations 3-2)
dfinissant le rgime critique vu au paragraphe prcdent :

On peut galement constater sur la courbe de la figure 3-2 :

0

Pour une nergie spcifique donne, le mme dbit, sil est infrieur au dbit critique, peut
scouler deux profondeurs diffrentes :

y correspondant au rgime torrentiel ;

y correspondant au rgime fluvial.

Sil est gal au dbit critique, on aura une seule profondeur qui est la profondeur critique.
Ce dbit est le dbit maximum qui peut transiter dans la section considre avec la charge
spcifique fixe.
Si la charge spcifique Hs reste constante, une augmentation du dbit se traduira par une
lvation de la hauteur deau y en rgime torrentiel, et par un abaissement de y en rgime
fluvial.

67

CY
I

fluvial

1-

-- ---

Figure 3-2 : Variation du dbit Q en fonction de la profondeur deau y pour une section
dtermine

3 REGIMECRITIQUE
3.1 Dfinitions

A partir du paragraphe prcdent, on tirer que le rgime critique dune section dtermine dun
canal faible pente, est caractris par plusieurs conditions importantes qui sont :
Lnergie spcifique est minimale pour un dbit donn ;
Le dbit est maximal pour une charge spcifique donne ;
Limpulsion totale ou force spcifique (voir chapitre 4) est minimale pour un dbit donn ;
Lnergie cintique (U2/2g)est gale la moiti de la profondeur moyenne ym;
Le nombre de FROUDE est gale lunit ;
0

La vitesse moyenne U de lcoulement est gale la clrit des petites ondes gravitaires
en eau peu profonde causes par une perturbation locale.

On appelle lments critiques les lments gomtriques de la section (S, P, RH,1, ...) calculs
pour une profondeur deau gale la profondeur critique et qui seront indics par la lettre c.

68

3.2 Calcul de la profondeur critique Yc

la profondeur critique est donne par la rsolution de lune des quations (3-3) o la surface S, la
profondeur moyenne ym, la largeur en gueule 1 et le nombre de FROUDE sont exprims en
fonction de la profondeur deau y :

Plusieurs mthodes de rsolution peuvent tre utilises pour calculer la profondeur critique dune
section dtermine : mthodes graphiques, abaques, mthodes analytiques pour les sections
paramtres.
3.2.1 Mthodes graphiques

On peut par exemple tracer points par points la courbe

en fonction de la profondeur

deau y et chercher labscisse pour la quelle cette courbe est gale

Q/&

(figure 3-3)

Y
Yc

F i m e 3-3 : dtermination de la profondeur critique par la mthode graphique.

3.2.2 Utilisations des abaques

Plusieurs types dabaques de calcul de la profondeur critique sont donns dans la littrature
hydraulique. La figure 3-4 est un exemple dabaque pour le calcul de la profondeur critique dune
section circulaire ou trapzodale :

69

Pour une section circulaire, on porte sur la graduation logarithmique du bas la valeur de

et on lit l'ordonne

Jg D?

de la courbe dnomme <(circular )) (D est le diamtre de

la section).
Pour une section trapzodale, on porte sur la graduation logarithmique du haut la valeur
de la courbe paramtre selon le b i t des berges m (b
de -et on lit l'ordonne
b
by
est la largeur au radier). On peut constater que ces courbes convergent pour les faibles

&

valeurs de A vers la section infiniment large o le fruit des berges n'intervient plus.

3.2.3 mthodes analytiques

Pour les sections courantes paramtres, la substitution dans l'quation (3-3) des valeurs des
lments gomtriques S et 1 en fonction de la profondeur d'eau y (Tableau 1-2) et la rsolution
algbrique conduisent aux formules du tableau 3-1 ci-aprs. Les fonctions implicites de la
profondeur critique yc peuvent tre rsolues avec les solveurs des calculatrices ou Z'EXCEL ou
bien par des mthodes numriques (voir chapitre 3 du tome 1)
3.3 Applications aux sections courantes
Exemple : Calculer la profondeur critique d'un canal rectangulaire ayant un dbit Q=1,3 m3/s ,
une rugosit &=70, une pente I=l% et une largeur au plafond b=l m.
Rponse : yc = 0,56 m
iExemple 2 : Calculer la profondeur critique d'un canal trapzodal ayant un dbit Q=lO m3/s ;
une largeur au plafond b=20 m, une rugosit Ks=70, une pente I=l%oet un fniit des berges m=0,5.
Rponse : yc = 0'29 m
Exemple 3 : Calculer la profondeur critique d'un canal circulaire ayant un dbit Q 4 , 2 m3/s, une
rugosit Ks=70, une pente I=l% et un diamtre D=l m
Rponse : 8, = 119.59' et yc = 0.25 m

70

M
+
l

QTfl

ml N

L$
II

4
II

ru:

e
O
*O
O

.-.(

.-.1

-e(

Valeurs de yJD ou y&

4 PENTES ET REGIME
4.1 Pente critique

Les discussions ci-dessus sur le rgime critique de lcoulement se rapportent principalement

une section particulire le long du canal appele section critique. La pente ne figure pas dans

la relation dterminant la profondeur critique. Donc si une profondeur est critique pour une
dbit donn et une forme de section dtermine, cette profondeur sera critique queue que soit
la pente.
Si le rgime critique existe tout le long du canal ou sur un bief, lcoulement est alors critique
sur ce bief. Puisque la profondeur critique dpend uniquement des lments gomtriques du
canal (surface S, profondeur moyenne ym et largeur en gueule 1) pour un dbit donn, la
profondeur critique doit tre la mme sur le bief et lcoulement est alors uniforme avec cette
profondeur.

La pente pour laquelle un dbit donn scoule une profondeur normale critique est
appele pente critique. Mais un tel coulement est instable comme nous lavons dit au
paragraphe 2.2.
4.2 Expression de la pente critique
Le calcul de la pente critique exige la rsolution du systme des 2 quations constitu de
lquation du rgime critique et de la formule du type Manning-Strickler :
(3-6a)

(3-6b)
On obtient ainsi deux relations entre le dbit Q, la profondeur critique yc ( travers Sc, 1, et
RH^) et la pente critique 1,.
0

Si le dbit est donn, la premire quation permet de dterminer la profondeur

critique yc et par consquent Sc=S(y,) et RH,=RH(Y,)et la deuxime quation permet
de calculer la pente critique (quation 3-7).

Ic=

Q2

KSS,~R~C~

(3-7)

Compte tenu de la relation entre le dbit et la section critique ( quation 3-6a) et de la

dfinition de la profondeur moyenne y,,, la formule (3-7) est quivalente lune des

73

relations (3-8) qui permet de calculer la pente critique directement avec la seule
profondeur critique est connue.
Si la pente critique est donne, le calcul de la profondeur critique se fait par la
rsolution de l'une des quations (3-8) par un solveur ou par une mthode numrique
ou graphique ;ensuite le dbit est dtermin par l'une des relations (3-6).

rc=
Ic=

(3-8a)

Y K i RH C ~
nc

(3-8b)

On peut tracer la courbe Qc=f(Ic) qui est une fonction implicite. Pour y arriver, on
suppose diffrentes valeurs de yc et calculer le dbit Q et la pente Ic correspondant
respectivement aux quations (3-6a) et (3-8). Le minimum IL correspond une pente
limite requise pour avoir un rgime critique pour le canal donn (Figure 3-4).

Figure 3-4 Variation du dbit critique avec la pente critique pour un Canal donn

4.3 TYPES DE CANAUX (CLASSIFICATION SELON LA PENTE)

On a dj tabli que :
La profondeur critique est indpendante de la pente 1 du canal pour une section
dtermine et un dbit donn ;
0

La pente critique 1, dans ce cas est la pente pour laquelle la profondeur normaie y,,
correspond la profondeur critique yc.

74

Donc il est intressant de connatre, dans le cas gnral, les relations entre la profondeur
normale et la profondeur critique (ou entre la pente 1 et la pente critique 1, dun canal) pour un
dbit donn et une forme de section dtermine.
Si la section est vase vers le haut, la profondeur normale diminue si la pente augmente. En
effet, le dbit constant Q est li la profondeur y et la pente 1par la relation suivante o la
dbitance du canal 96)est une fonction croissante de y :

Donc pour un dbit fix, si la pente 1 est suprieure la pente critique Ic, ~(y,,)sera infrieure
cp(yc) ;do la profondeur nurmale sera infrieure la profondeur critique. Cest le contraire
qui se produit si la pente 1 est suprieure la pente critique L.

Selon la position de la profondeur normale yn par rapport la profondeur critique yc (ou la

pente 1 par rapport la pente critique Ic), on classifie les canaux en leur attribuant une lettre
majuscule latine qui dsigne ce positionnement relatif (Tableau 3-2).
Cette classification (type de canal) est distinguer de lcoulement rel qui y a lieu (Tableau
3-3). Par exemple on peut avoir un coulement graduellement vari dans un canal fond plat
mais pas dcoulement uniforme. On peut galement avoir un coulement graduellement vari
torrentiel dans un canal fluvial (cest lcoulement normal qui est fluvial).
4.4 Applications

Exemple 4 : Dterminer la pente critique et le type de canal dans les exemples 1, 2 et 3 du

paragraphe 3.3 ci-dessus.
Rponses
ExemDle 1 : Ic=6,6%et type de canal=T ;
ExemDle 2 :I,=3,0% et type de canal=F ;
Exemde 3 : IC=4,6%0
et type de c d = T
Exemple 5 : La pente critique dun canal trapzodal, de largeur au radier b=Q,5 m , de
rugosit Ks=70 et de b i t des berges m=l ,est de 4%0.
1 O ) Quelle est sa profondeur critique yc ?
2) Quel est son dbit Q ?

Rponses
1O ) yc=0,50 m
2) Q=894 Vs

75

Tableau 3-2 : Classification d

canaux
Caractristique de
lcoulement
normal

Symbole

Appellation

canal IJorizontal
ou canal fond pla

Lcoulement
normale nexiste pas
(y,=w)

Lcoulement
normale nexiste pas
( y , indtermin)
1< O

Canal Adverse ou
canal contre penh

y, non dfinie

Lcoulement
normale est fluvial

Lcoulement
iormale est torrentiel

Lcoulement
normale est critique

Canal Fluvial ou
canal lent ou canal
pente faible

Canal Torrentiel ou
canal rapide ou
canal pente forte

Canal critique

76

5 COURBES DE REMOUS
La courbe de remous est le profil en long de la surface libre en fonction de labscisse x de la
section considre.
Pour les canaux prismatiques, il sagit dtudier la fonction qui lie la profondeur deau y
labscisse x de la section puis que la pente est constante sur un bief et la section mouille ne
dpend que de la profondeur.

Par contre pour les cours deau naturels ou les canaux non prismatique, la pente et la section
mouille peuvent varier dune section lautre et on cherche plutt dterminer la fonction
qui lie la cote h = zf+ y de la surface libre (par rapport un plan de rfrence horizontal)
labscisse de la section x ;o zf est la cote du fond.

h=G(x)
Ces fonctions ne sont pas connues priori mais on peut tablir lquation diffrentielle qui
donne leur pente (dy/dx ou W d x ) par rapport une droite de rfrence .
A partir des drives de la charge spcifique Hs par rapport labscisse x (quation 32) et par rapport la profondeur y (quation 3-4) on tire lexpression (3-9) donnant la
drive dy/dx qui est pente de la surface libre par rapport au fonds pour les canaux
prismatiques.

1-J
2
I-FR

(3-9)

o J , la pente de la ligne dnergie, est donne par une formule

du type Manning-Strickler daprs les hypothses sur les
coulements graduellement varis :

J=

Q2
4

K~ZS~R~
A partir directement de la drive de la charge hydraulique, on obtient lexpression (310) donnant la drive dh/dx qui est la pente de la surface libre par rapport un plan
de rfrence horizontal pour les chenaux non-prismatiques.
(3-10)

77

o k est un coefficient tenant compte des pertes de charge

singulires sur le tronon dx (k=O si lcoulement est convergent
et O<k<l si lcoulement est divergent).

Pour les canaux prismatiques, ltude de la fonction dyidx donne lallure de la courbe de
remous et lintgration de lquation diffrentielle donnera la courbe si lon connat une
condition initiale appele section de contrle.
5.1 Etude qualitative des courbes de remous dans les canaux prismatiques
Ltude qualitative et la classification des courbes de remous permettent de mieux contrler le
sens des calculs numriques et de fixer les conditions darrt de ces calculs.
5.1.1 Sens et domaine de variation de la profondeur deau

Le signe de dy/dx est le produit des signes de dH&lx et dHddy qui dpendent respectivement
de la position relative de la pente de la ligne dnergie J par rapport la pente du canal 1 dune
part (ou de la profondeur deau y de la courbe de remous par rapport la profondeur normale
y,, du canal) et de la position relative de y par rapport la profondeur critique yc dautre part :
La drive d H ~ / d xest positive si y est suprieure Yn (ou J infrieure 1) et elle est
ngative si y est infrieure yn (ou J suprieure 1) ;
La drive dHs/dy est positive si y est suprieure yc et elle est ngative si y est
infrieure yc (voir courbe Hs(y) de la figure 3-1)

Les domaines de variation des courbes de remous sont dsigns par un chiffre 1 ou 2 ou 3
selon la position de chaque courbe par rapport aux deux profondeurs normale et critique du
canal :

Si la courbe est situe au-dessus des deux valeurs yn et yc, on lui affecte le chiffre 1;
Si la courbe est situe entre les deux valeurs yn et yc, on lui affecte le chinre 2 ;
Si la courbe est situe au-dessous des deux valeurs y,, et yc, on lui affecte le chiffre 3.

Les canaux tant dj classifis selon la position relative de la profondeur normale yn et de la

profondeur critique yc (Tableau 3-2) on peut alors distinguer les diffrents types de courbes
de remous par une lettre latine majuscule (type de canal) suivi du chifie lou 2 ou 3 (domaine
de variation).
On peut ainsi distinguer, comme le montre le tableau 3-3 douze types de courbes de
remous (FI ;FZ;F3 ;TI ;T2 ;T3 ;Cl ;Cz(c0ulement uniforme critique) ;C3 ; A2 ;A3 ;HZ;H3 )
avec les sens de variation des profondeurs y, les types ou rgimes dcoulement (positions de
y par rapport yc) distinguer des types de canal (on peut avoir un coulement torrentiel dans
un canal F cornme la courbe F3) et des exemples doccurrence.

78

5.1.2 Comportement de la profondeur deau aux limites

11 est intressant de c o d t r e le comportement des courbes de remous au voisinage des
limites des 3 domaines dfinies ci-dessus (O, y, yc et 00) :

Les courbes de remous tendent asymptotiquement vers la profondeur normale yn.

En effet, lorsque y tend vers yn, la pente J de la ligne dnergie tend vers la pente 1 du
canal ;do dy/dx tend vers O si le dnominateur de lquation (3-9) reste fini.
A lapproche de la profondeur critique yE, on na plus un coulement
graduellement vari pais plut& une variation brusque de la ligne deau (ressaut
ou chute) que nous tudierons au chapitre 4. En effet si y tend vers yc, cest le
dnominateur de lquation (3-9) qui tend vers O et dy/dx tendrait v m linfini si le
numrateur reste fini. La ligne deau serait alors normale la droite y=yc (Figure 3-5) ;
do une courbure trs prononce qui est contmrire lhypothbse de f ~ b courbure
e
des lignes de coufant faite pour les coulements graduellement varis. On obtient ainsi
un coulement brusquement vari (ressaut ou chute) qui nest plus rgi par les mmes
quations.

Lorsque la profondeur deau y tend vers les grandes valeurs, la surface libre tend
vers lhorizontal. En effet lorsque y tend vers linfini la pente de la ligne dnergie
tend vers O ,donc cette ligne tend vers lhorizontal. Comme la vitesse aussi tendra vers
O, la surface libre est confondue avec la ligne dnergie et tendra vers lhorizontal et
dy/dx tendra vers la pente du canai 1.
Lorsque la profondeur y tend vers les faibles valeurs, lindtermination d a o est
leve par une condition aux limites physique qui dtermine lorigine du dbit.

F i m e 3-5 : comportement des courbes de remous lapproche de la profondeur critique

79

FI

'$3.

di wm

Aq

Tableau 3-3 : Classification des courbes de remous

80

5.2 Mthodes de calcul des courbes de remous

Le calcul de la courbe de remous est bas sur la rsolution de l'quation diffrentielle (3-9) ou
(3-10) qui rgit les coulements graduellement varis et que nous pouvons remettre sous l'une
des formes (3-1 1). Il faut une condition limite c'est dire une profondeur y0 ou cote ho pour
une section xo donne, appele section de rfrence ou section de contrle.
Cette section est localise l'aval pour les courbes de types FI, F, y C , y T, , H, , A, et la
courbe-de remous est calcule de I'aval vers l'amont ; c'est dire quand l'coulement est
fluvial. La section de contrle est localise $ l'amont et la courbe de remous est calcule de
l'amont vers l'aval pour les courbes de remous de types T, , F3, T3, C, , H, et A, ; c'est
dire quand l'coulement est torrentiel (Tableau 3-4).

pour les canaux prismatiques avec I>O

K,=p(y)=$

(3-1 la)

est la dbitance normale ;

"
-

K=p ('y)=K,SRH3
est la dbitance pour la profondeur y ;

Q est le facteur de section critique

z,=-Jg
et Z=

(T)

est le facteur de section pour la profondeur y.

--I
dY

di-

pour les canaux prismatiques avec I>O

(3-llb)

@)

o Q , = K ( v ~est le dbit normal correspondant la profondeur y

et Q, = Z ( y ) & est le dbit critique correspondant la profondeur y.

81

l-@)

--I
dv

+)
2

dx1-r(

pour les canaux prismatiques avec I<O

(3- 11c)

o r=I avec In, tant la pente critique correspondant la profondeur

I"C

normale yn calcule par le dbit rel Q.

-=F(x,y)
dY
a?x

pour les canaux non-prismatiques

(3-1 Id)

On peut distinguer 3 mthodes : la mthode d'intgration graphique, la mthode d'intgration

directe et la mthode des diffrences finies ou par itration.
5.2.1. Mthode d'inttigration graphique
La procdure d'intgration graphique directe est utilise pour le calcul. Considrons 2
sections du canal l'abscisse x1 et x2 auxquelles correspondent respectivement les
profondeurs d'eau y1 et y2.
L'intgration de cette quation variables spares conduit :

On peut donc tracer points par points la courbe l/f(y) et estimer graphique l'aire hachure la
figure 3-6 qui dtermine la distance entre y1 et y2. Cette aire est dtermine par planimtrie ou
en utilisant la mthode des trapzes (remplacer la courbe par des segments de droite entre les
diffrents points) qui est quivalente la mthode numrique des pas directs que nous verrons
plus loin.
Il existe d'autres mthodes graphiques dont celle des diagrammes universels de Silber et celle
de Raytchine et Chatelain qui sont rsumes dans Carlier ou celles de Grimm, d'Ezra et de
Leach qui sont prsentes dans Lencastre. Elles ne sont plus beaucoup utilises avec
l'avnement des ordinateurs de poche et de bureau et la rapidit des moyens de calcul
modernes.

5.2.2 Mthode d'intgration directe

L'intgration analytique de l'quation diffrentielle sous les formes (3-1 1) s'avre
pratiquement impossible dans le cas gnral. Mais plusieurs tentatives de rsolution de
quelques cas spcifiques avec l'introduction d'hypothses simplificatrices ont t faites par
divers auteurs (V. T. CHOW p253. Les toutes premires mthodes (Bresse, BakhmetefX)

82

utilisaient la formule de Chzy pour le calcul de la dbitance K(Y)=CRH% tandis que les
dernires utilisent la formule de Manning-Strickler K(Y)=KsRH~~S.
Dans la mthode de Bakht-neteff modifie par V T CHOW, on introduit gnralement des
exposants hydrauliques N et M pour exprimer les dbitances K et KN et les facteurs de section
Z et ZC en fonction des profondeurs deau y et des constantes Cl et CZ(quation 3-12)

Z2=C2yM
et Zc2=C2Yc
M

(3-12)

La substitution de ces formules (3-12) dans lquation diffrentielle (3-1la) et le changement

de variables U=y/yn conduisent une formulation dont lintgrationdonne lexpression (3-13)
o F(u,N) est une fonction dcoulement graduellement vari, tabuie dans beaucoup
douvrages hydrauliques, mais qui peut se calculer facilement avec la TI92 maintenant.

(3-13a)

et F(u,N)=[--,d.c

1-X

La distance sparant les 2 profondeurs y1 et y2 peut donc sexprimer ainsi :

Pour la mthode de Bresse (hypothse de canal infiniment large avec coefficient de Chzy
constant), on prend N=M=J=3 pour les exposants hydrauliques et la fonction F(u,3) peut se
calculer analytiquement :

u - 1)

83

La procdure de mise en uvre consiste :

0
0

Calculer les profondeurs normale yn et critique y, pour le dbit Q et la pente 1 du canal ;

Derminer les exposants hydrauliques N et M pour une valeur moyenne de la
profondeur deau y estime sur le bief en considration. On suppose que M et N sont
constants.

Calculer0

Calculer u et v aux 2 sections dlimitant le bief

Trouver les valeurs de F(u,N) et F(v,O) partir des tables ou dun calcul dintgrale
avec la TI92
Calculer la longueur du bief par lquation (3-13b)

La procdure ci-dessus peut tre applique pour les c m u x pente adverse @<O) condition
de modifier lquation (3-13a) comme suit :

tcste

(3-13~)

et y,, est calcul avec la valeur absolue de la pente 1.

Quand la pente est nuile ( P O ) , on pose le changement de variable p=/yc et lintgration de
lquation diffrentielle originale donne la formule (3-13d).
(3-13d)

Quand les exposants hydrauliques ne sont plus constants comme on la suppos (canal avec
changement brusque de section ou non vas vers le haut comme un canal circulaire), on
subdivise le canal en biefs courts o N et M sont constants mais peuvent diffrer des
exposants pour les profondeurs normale Nn et critique M,. la formule (3-13a) est modifie en
(3-13e).
(3- 13e)

84

avec u

Y= ;
YnN

F uN/O;
et O=N/(N-M+l).

Pour les conduites circulaires, on peut tablir les relations (3-130 et les substituer A dans
l'quation diffrentielle (3-1 la) .

(3-13f)

et l'intgration de l'quation diffrentielle donnera la formule de l'abscisse de la section de

profondeur y (quation 3-13g) o les fonctions d'coulement Fi et F2 sont tabules et peuvent
maintenant tre calcules facilement avec la TI92.

et 8=2arccos( 1-2x)

85

5.2.3 Mthode itrative ou des diffrences finies

Elle est caractrise par la subdivision du canal en tronons courts et par une progression pas
par pas du calcul. Il ya beaucoup de variantes de cette mthode. Quelques unes sont
prfrables dautres selon les circonstances mais aucune delles ne se trouve tre la
meilleure dans toutes les applications.
La mthode de variation des profondeurs ou pas directs (Ay fix lavance et Ax calcul)
est la plus simple mthode par diffrence finie applicable aux canaux prismatiques.
La mthode des tronons ou mthodes des pas standards (Ax fix A lavance et Ay calcul)
sapplique galement aux canaux prismatique est trs bien adapte aux chenaux naturels o
les lments gomtriques et hydrauliques ne sont connues qu certaines sections o on les a
releves.
Ces 2 mthodes ont chacune des variantes sur la faon dvaluer les valeurs moyennes des
fonctions sur le bief Ax mais galement le nombre de points intermdiaires ou antrieurs pour
valuer ces fonctions.
Dun point de vue numrique, les mthodes que nous prsentons ci-aprs sont au plus dordre
2 (erreur systmatique dordre Ax3 ou Ay3 selon les cas). Il existe des mthodes dordre plus
lev telle que la mthode de Gil-Runge-Kutta dordre 4 et de rang 4 (voir annexe) qui sont
utilises pour les quations diffrentielles du premier ordre condition initiale.
Mais compte tenu de limprcision sur les donnes hydrauliques (rugosit, dbit, forme de la
section pour les chenaux naturels, pente, etc.), les mthodes dordre 2 suffisent avec un bon
choix des pas de calcul (Ax ou Ay) et du sens de progression des calculs.
Le sens de progression du calcul pour toutes ces mthodes (dcrit au tableau 3-4) doit tre
rigoureusement observes, sinon on tend invitablement vers des valeurs divergentes de la
courbe de remous relle.
5.2.3.1 Mthode des pas directs
Cest la mthode la plus simple pour le calcul de la courbe de remous pour un dbit donn et
un canal prismatique. Soit xo une section o la profondeur deau y0 est connue et une section 1
situe une distance inconnue Ax de xo o lon veut que la profondeur soit yi=y~+Ay.
La distance Ax est estime en supposant que la pente de la ligne dnergie est la pente
moyenne sur le bief dans la discrtisation de lquation diffrentielle (3-9). Cette hypothse
donne dassez bons rsultats pour des Ay faibles.

ou

J = JYi+YO
(y)

86

avec pour la formule de Strickler JO)=

Q2

K ,"Soi>'RHW+
La discrtisation de l'quation diffrentielle s'crit alors :

Les valeurs positives de Ax indiquent une progression vers l'aval et celles ngatives une
progression vers l'amont.
Le tableau 3-4 donne les procdures typiques souvent utilise partir d'une section de
contrle du canal. Toutefois on peut procder suivant une direction diffrente de celle
indique au tableau si la profondeur est connue l'autre bout de la courbe avec les risques de
divergence dj mentionns.
La mise en uvre de la procdure de calcul passe par les tapes suivantes :
1") Connaissant yo, calculer HSOet JQo)
2") Pour une valeur de YI= yo+Ay, calculer &let J(y1
3") Calculer7=(J(yo)+ J(y1))/2 ou 7=J((yo+y1)/2)
4") Calculer Ax
5") Passer la section suivante en prenant un nouveau y0 g u au y1 qui vient d
',*
continuer la procdure en 1")

trouv et

Le calcul de AHs=Hsi-Hso est imprcis (gnralement une diffrence faible entre 2 grands
nombres) et la vitesse doit alors tre calcule avec une grande prcision (3 dcimales au
moins). Il est prfrable de calculer AHs par son expression andytiqueAfj,=A&-g2)

-2
o F~

=&g s 3 est value pour la profondeur moyenne Yo+Yl entre les 2 sections.
2

L'usage de pas faibles permet d'avoir une plus grande prcision mais augmente les itrations
et le temps de calcul. Les erreurs se propagent comme dans toutes les mthodes itratives puis
que la solution d'un pas est incorpore dans le pas suivant (le courbe peut osciller dans
certains cas.

87

Pente du
canal

Relations entre
profondeurs
Y d Yn >Yc

Type de
courbe de
remous

Signe
de
dHJdx

FI

> Y+

Yc

F2

Y d Yc

Yn

ci

Yn

I>O

Y d Yc> Yn

Ti

> Yd Yn

T2

T3

H2

Yc

Y d Yc

I=O

-TYc > Y0

Y d Yc

A2

Yc > Y0

A3

I<O

Procdure de calcul
Section de contrle laval et on
progresse vers lamont (Ax <O)
Ay <O
Limite=y,
Section de contrle laval et on
progresse vers lamont (Ax <O)
Ay >O
Limite====,,
Section de contrle lamont et on
progresse vers laval (Ax >O)
Ay >O
Limite==c;ressaut ver5
lcoulement uniforme
Section de contrle laval et on
progresse vers lamont (Ax <O)
Ay <O
Limite-y,= yc
Ecoulement uniforme
y
- Yc-- Yn
Section de contrle lamont et on
progresse vers laval (Ax >O)
Ay >O
LimiteTcTn
Section de contrle laval et on
progresse vers lamont (Ax <O)
Ay <O
LimiteTn ;ressaut
Section de contrle lamont et on
progresse vers laval (Ax >O)
Ay <O
Section de contrle lamont et on
progresse vers laval (Ax >O)
Ay >O
LimiteT,
Section de contrle laval et on
progresse vers lamont (Ax <O)
Ay >O
Limite=valeur fixe sinon pas
Section de contrle lamont et on
progresse vers laval (Ax >O)
Av >O
Lhite-y, :ressaut
Section de contrle 1 aval et on
progresse vers lamont (Ax <O)
AY
Limitmaleur fixe sinon pas
Section de contrle lamont et on
progresse vers laval (Ax >O)
Ay >O
Limite-Y, :ressaut

rableau 3-4 :Procdure de calcul des courbes de remous

88

5.2.3.2 Mthode des pas standards

Les pas Ax sont fixs priori; ce qui permet de collecter les donnes de terrain des
abscisses et biefs choisies avant de procder au calcul de la courbe de remous pour les
chenaux naturels.

II sagit alors de dterminer la cote de la surface libre ces abscisses bien dtermines ey on
procde par approximations successives en tenant compte du sens de progression des calculs
indiqu au tableau 3-4). Lquation aux diffrences finies peut scrire sous la forme
suivante :

h + LV=2h o + a,-O_v 2 -JAx-sign(Ax)

a2g
2g
La prise en compte des singularits se fait travers la diffrence de charge cintique et k qui
est un Coefficient de perte de charge singulire avec :
k=O pour un canal prismatique ;
0
k=O 0,l pour un coulement convergent graduel
0
k=0,2 pour un coulement divergent graduel
0
k=0,5 pour un changement brusque de la section de
lcoulement.
Dautres auteurs expriment la perte de charge singulire par un coefficient k et la charge
cintique la section O seulement
ou bien par la diminution du coefficient de Strickler sur le
bief dans lexpression de J qui se calcule de la mme manire quau paragraphe 5.2.3.1.
Pour la mise en uvre de la mthode, on crira lquation de Bernoulli ci-dessus sous cette
forme :
(3-14a)
(3- 14b)
2

(3- 1 4 ~ )
On procde alors comme suit :
1) Calcul de HOpar lquation (3-14b), connaissant ho
2) Choix initial de hlcomme valeur dessai ;
3) Calcul de Hl par lquation (3-14c) en passant par y, S , v=Q/S et aV2/2g pour la
valeur dessai de hl ;

89

4) Calcul de Hl par lquation (3-14a) en passant par y, RH, RH^'^, J et laperte de charge
singulire ventuellement ;
5 ) Sil ya accorde entre les 2 valeurs de Hl trouves (diffrence en valeur absolue lAH,(
infrieure une petite valeur E), la cote dessai hl set la bonne ; sinon on modifie hl
par une nouvelle valeur et on recommence ltape 3).

Henderson, en supposant que J est proportionnel yJ, suggre de prendre la correction

suivante sur hl pour acclrer la convergence :

AH^ est la diffrence entre les 2 charges trouves par le 3) et

le 4) ;

F RI

a -1 V12

Wx
2g pour un cours deau naturel o I d
__

2
Croley propose une correction sur y du type suivant pour arriver la mme fin :

o lexposant (1) indique la valeur dessai et lexposant (2) la

valeur corrige.
On peut galement tester si la diffrence
est en valeur absolue infrieure E
comme critre darrt et prendre comme nouvelle valeur dessai de h la dernire valeur
trouve hi(+)comme dans la mthode du point fixe. Il faut signaler que hl(+)est donn par
le thorme de Bernoulli :

Lorsque lnergie cintique aV2/2g est faible, on peut conduire les calculs mme dans la
mauvaise direction ( vers lamont pour un coulement torrentiel ou vers laval pour un
coulement fluvial) sans commettre de graves erreurs.
Aussi pour les cours deau naturels, la section de rfrence ou section de contrle nest pas
toujours entirement dfinie pour un dbit donn. Si on se fixe un ho qui est incorrect pour le
dbit donn, la courbe de remous sera peu prs correcte aprs quelques pas de calcul dans la
bonne direction. On peut ainsi commencer une station loin en amont ou en aval du bief qui

90

nous intresse pour le calcul de la courbe de remous sur celui-ci si on ne connat pas la section
de rfrence.
Enfn on peut constater que plus le coeflcient de Strickler est lev, plus la courbe de remous
est longue et plus il est faible, plus la courbe est courte. Do la valeur par excs de Ks
devrait tre choisie dans les calculs si la connaissance de la plus grande influence du remous
est recherche (remous cr par un barrage par exemple). Par contre pour la connaissance de
la plus courte courbe possible, une plus faible valeur de Ks devrait tre utilise (amlioration
de la navigation par exemple puisque la profondeur navigable doit tre suprieure ii une
valeur donne et la plus courte courbe indiquera la plus faible profondeur une section
donne du chenal).

5.2.4 Mthode de la Dnivele

- Dbit pour les cours deau natuml

On peut utiliser cette mthode lorsque lon dispose de profils de la surface libre en
coulement normal (sans leffet de remous) pour un certain nombre de dbits. Elle est simple
et conome en calcul.
Lapplication du thorme de Bernoulli, entre les sections amont 1 et aval 2 spares de la
distance L, donne la pente de la ligne dnergie sur le bief.
J=-

A+S

(3-1%)

o
A est la dnivele (hl - h2) entre les cotes de la surface libre ;

inclut la variation de la charge cintique

et les pertes de charge singulires entres les 2 sections.
1) Si les nergies cintiques sont faibles, on peut ngliger 6

Lquation (3-1Sa), avec la formule de Manning-Strickler, donne

ou

Q = f(h1)
fi

(3 - 15b)

Pour une courbe de remous avec un dbit Q, on peut utiliser la mme quation avec une
dnivele A et trouver la relation (3-15c) suivante :

91

(3-1%)
Al=[%]

Le dbit Q est appel dbit par unit de dnivele. On peut alors utiliser la relation
(3-1%) pour le calcul de la courbe de remous si lon connat la courbe de la relation
(3-15b) d-t

Q en fonction de la cote amont hi. Cette courbe peut tre trace

.Jd

partir de donne de jaugeage (dbit et dnivele entre 2 chelles). Elle est trace comme
une courbe moyenne pour diffrentes conditions de fonctionnement (crue, dcrue,
pousse des herbes,. ..).
Pour une cote hl donne et un dbit Q, elle permet de lire

2 et la formde (3-1%)
$h

permet le calcul de A ; do la valeur de h2 = hl i A selon le sens de progression du

calcul. On peut galement calculer le dbit Q si la dnivele A est lue sur 2 chelles.
Lassimilation de 6 O fait que la mthode est plus utilise pour des problmes o les
vitesses sont bien en de de la vitesse critique et diminuent vers laval (nombre de
Froude faible).
2) Si les nergies cintiques ne sont plus ngligeables

La mme approche peut tre utilise par approximations successives pour la

dtermination du dbit partir de la mesure de la dnivele A = hl - h2. La procdure est
la suivante :
a) calculer K, S R F pour hl ou pour h2 ou par une moyenne (certains auteurs

SI^" RG: (Ks~S Z ~R$:

.

utilisent la moyenne gomtrique (Ks,

b) Initialiser6=0

(J=g)

c) Calcul de Q par la formule de Manning-Strickler. Q = K, S R H f i

d) Calculer VI et V2 avec le dernier dbit trouv Q.
e) Amliorer J par lquation (3-15a) en calculant 6 par VI et V2.
f ) Calculer un nouveau dbit avec la formule de Manning-Strickler et le J

amlior.
g) Tester si la variation du dbit est, en valeur absolue, infrieure une faible
valeur E. Si oui arrter ;sinon aller d).

92

5.2.5 Exemple de calcul de courbe de remous

Un canal trapzodal vhicule un dbit Q=lO m3/s avec une largeur au plafond b=20 m, une
rugosit Ks=50, une pente I=O,1 %O et un fhit des berges m=0,5. la profondeur d'eau une
station x=O est dtermin par un seuil et gale y0=0,5 m.
1 O ) Calculer les profondeurs normale et critique du canal.
2") Dterminer le type de courbe de remous et ses caractristiques.
3") Calculer la courbe de remous.

Rponses
1 ") yn=l ,O2 m ;yc=0,29 m
2") C'est une courbe de remous de type F2 puisque yn>yo>yc. La section de contrle se trouve
l'aval et on fait le calcul vers l'amont (Ax<O). La courbe F2 diminue la profondeur d'eau de
l'amont vers l'aval ;donc avec le sens de progression des calculs, Ay>O. On arrte les calculs
la limite yn=l ,O2 m.
3") Les tableaux 3-5 montrent les rsultats des calculs par diffrentes mthodes avec un
Ay=0,4 m et les diffrences sont minimes.

93

Calcul de la courbe de remous par la mthode des pas directs (delta y fix), systme SI
b
20

m
0,5

Ks
50

alpha

Delta y

0,0001

10

0,04

74,54805576 -169,53692!
99,47492848~-269,01185
131,78@3923-400,79254i
174,3592525675,15180:
-231,888251 -807,04005:
312,4602648-1119,5003:
431,4007049 -1550,901O:
621,6775878-2172,5386
969,54842383142,1270;
1795,2295564937,3565!

Tableau 3-5a : Mhodes des pas directs. Le calcul de ANS peut tre imprcis (gnraiement une diffrence faible entre
-2 grandes valeurs, cest
pour quoi il faudra calculer la vitesse avec prcision (3 dcimales). 11 est plus prcis de calculer A&=AAl-z) o p R2 est la valeur moyenne
entre les 2 sections comme dans le tableau 3-5b.

Calcul de la courbe de remous par la mthode des pas directs (hita y fix), systme SI

b
20

Y
0,5
0.54
0,58
0,62
0,7
0,74
O,78
0,82
01%
0.9
OS94
0,98

i
1

m
O, 5

Ks
50

10,5352
11,3568
12,18
13,0048
13,8312
14,6592
15,4888
16,32
17,1528
17,9872
18,8232
19,6608

21,1627553
21,2521981
21,3416408
21,4310835
21,5205262
21,6099689
21,6994117
21,7888544
21,8782971
21,9677398
22,0571825
22,1466253

I
0,0001

Q
10

-22
RH-^

j7-e

Delta y
0,04

alpha
1

1--

Q'i
g

s3

0,49781797 0,94919888 0,8211125

0,53438237 0,88052973 0,85691762
0,57071526 0,82101806 0,88378656
0,60681953 10,76894685 0,90434016
0,642698W 1 0,72300307 0,92032869 1
0,67835359 0,68216547 1 0,93295137 1
0,71378894 0,64562781 0,94304852 f
0,7490068 0,6127451 0,95122089
0,78400983 0,5829952 0,95790567
0,81880067 0,5559509 0,96342614
0,85338188 1 0,5312593 0,96802487
0,88775602 0,5086263 0,97188625

J-

I-7

2-

K s R H 3

0,000913441
0,000715174
0,00056956
0,00046037
0,000376991
0,000312295
0,000261375
0,000220786
0,000188059
0,000161397
0,00013947
0,000121283

1
1
1

1
1
1
1

-0,000813441
-0,000615174
-0,00046956
-0,00036037
-0,000276991
-0,000212295
-0,000161375
-0,000120786
-8,80592E-05
-6,13969E-O5
-3,94703E-05
-2,12829E-05

40,37725541
-55,71869163
-75,2863676
-100,3792038
-132,9037739
-175,7842025
-233,7527446
-315,0096846
-435,1189178
-627,6705195
-981,O1 50419
-1826,607889

-40,3772554
-96,095947
-171,382315
-271,781518
-404,665292.
1-580,449495*
-814,202239
-1129,21192-1564,33084
-2192,00136
-3173,0164
-4999,62429

'rableau 3-Sb : Mthodes des pas directs. Valeurs moyennes calcules avec l'estimation la profondeur (yi+yi+i)/2.

oc

Ksi

20

0,5

50

Paramtres du canal
t
Q
0,Oooi

10

1
1

apha

DdtaX

variable

hrfanmce

i
i

915

.-

I
e

H-

_
I
_

CG
-0,C

-0.t

-0 c

-OC

-0,c
-0.C
6,C

6t
-0,c
-0 c
4

-0.C
i

Tableau 3-5c : Mthodes des pas standards, Les Ax sont variables et choisis de faon ?
comparer
i
les rsultats avec ceux des tableaux 3-5a et

96

ANNEXES DU CHAPITRE 3
A ECOULEMENTS A DEBIT VARIABLE DANS UN CANAL NON-PRISMATIQUE

Les coulements permanents variables dans lespace ont un dbit non uniforme (coulement
non conservatif) dcoulant dun apport ou dune perte deau le long du chenal dcoulement.
On les rencontre dans les fosss latraux des routes, les vacuateurs de crue coulements
latraux, les canaux dirrigation et de drainage non revtus, les canaux dvacuation autour
des bassins de traitement deau, les biefs avec dversoir latral dans un canal, etc.
Les apports ou les pertes deau causent une modification de lnergie ou de la quantit de
mouvement de leau du chenal. Le comportement hydraulique est ainsi plus compliqu que
dans le cas des coulements conservatifs. En effet lcoulement avec apports latraux se
comporte diffremment dans une certaine mesure de celui avec pertes latrales.
Les coulements avec pertes latrales peuvent tre traits comme une dviation dune partie
de leau qui naffecte pas la charge hydraulique de lcoulement restant et lexprience prouve
cette hypothse. On peut alors utiliser le thorme de Bernoulli pour ltude de ce type
dcoulement tout comme le thorme de quantit de mouvement.
En coulement avec amorts latraux, une part importante de la perte dnergie est due la
turbulence cause par le mlange de leau entrant et de leau scoulant dj dans le canal.
Pour appliquer le thorme de Bernoulli il faut ajouter lapport dnergie cintique latral
Les hypothses dtude sont les mmes que pour les coulements graduellement varis :
lcoulement est unidirectionnel (les courants latraux impulsant un mouvement
hlicodal leau et une surface libre mouvemente dans la section droite sont
ngligs et on tient compte de leurs effets en appliquant le thorme des quantits de
mouvement).
les coefficients de Corriolis a et de Boussinesq p sont constants et gaux lunit.
la rpartition des pressions est hydrostatique dans une section.
la pente du canal est relativement faible.
la pente de la ligne dnergie est donne par une formule de type Manning-Strickler.
ieffet de lentranement dair est ngligeable.
Le thorme des quantits de mouvement appliqu sur la tranche de canal dx donne :
P(Q +kdQ) CV

+ dv)

- P Q = - P g d (SYd + P g S dx 0 - 9

o k = O pour les coulements avec perte (dbit latral dQ ngatif) ;

97

et k=l pour les coulements avec apport (dbit latral dQ positif et

perpendiculaire laxe dcoulement x).

En ngligeant les termes du 2 ordre, cest dire :

dV=dQ_
s gs2
(ldy+g),

dV dQ ngligeable ;
Lquation exprimant le thorme des quantits de mouvement devient, aprs
simplifications :
Q2
3s
1 - J - (1+ k) Qs+ dy-gs2 gs3d x
dx
1
1 - Q2
gS3
~

o q est dbit latral (positif si apport et ngatif si perte).

as = O pour un canal prismatique.

et le terme ax

A. 1 Solutions analytiques

Avec des hypothses simplificatrices, on peut trouver des solutions analytiques cette
quation (A-1) dans des cas simples :
A.I.1 Ecoulemenf dans un canal rectangulaire de faible pente et de hiternent
ngligeable avec une grille BU hnd (figure A-1):
Si fcoulement d travers lu griffeest vertical (grille avec barreaux), on aura :
q = E C b ,/2gHs

o E est le pourcentage douverture de la grille et C est le coefficient de

dbit travers les ouvertures.

Le dbit dans le canal est donn par : Q = by ,/2g(H, - y)

98

On trouve alors la solution:

o la constante C=(HS /

C) (y1 / Hs)

est donne par y=y1 pour

x=o.

La longueur pour laquelle y = O est donne par :

Si lcoulement travers la grille est inclin, alors

q = Ecb

et lintgration de lquation diffrentielle donne

ou

La constante dintgration Ctei ou Cte2 est donne par y=yi x=O ;et la longueur pour
laquelle y = O est donne par :

99

-t

F i m e A-1 : Ecoulement dans un canal rectangulaire de faible pente et de frottement

ngligeable avec une grille au fond
A.1.2 Ecoulement sur un dversoir latral dans un canal rectangulaire
horizontal de frottement ndgligeable (figure A-2)

On suppose que l'coulement travers le dversoir est perpendiculaire (coulement fluvial).

q =C

fi

(y - P)3n o C est le coefficient du dbit du dversoir.

On trouve comme solution de l'quation (A- 1) :

c (d

x=-F

- +Cte

F( L)= 3Hs-3P +-------3arcsin,/----Hs-y

Hs-y
Hs
Hs-P
y-P
Hs-P

Figure A-2 : Ecouiement sur un dversoir latral dans un canal rectangulaire horizontal de
frottement ngligeable

100

A.2.3 Canal dun vacuateur de crue latral de penfe nulle, de frottement

ngligeable et de longueur L (fisure A-3)

q est constant

1 13

au dbouch du canal

La solution de lquation (A-1) est donne, en tenant compte de y = yc pour x=L,par :

l

qmonstants

----.il-

7 c ;ilii i i i i Cc i i t i il i i i + +ICI
ii il i l
!$
l
Y01

I_

_-

t$

-- -

__

-_- -..

- -..

c
I

._

yc

//////,~////////////////////~
-

L
__ __

__

-- - --- ..

__-

--X

Fimue A-3 : Canal dun vacuateur de crue latral de pente nulle, de frottement ngligeable et
de longueur L
A. 1.4 Ruissellement sur une sutface plane &ultant dune pluie

Le dbit par unit de largeur q la so:-tie est donne par :

q=Kym
o y est la profondeu deau la sortie

La valeur de m dpend du type dcoulement

turbulent et 3 pour u n coulement laminaire).

(%

5/3 pour un coulement

101

Le dbit entrant est donn par :

q =. x i
o x est la distance par rapport la crte de la surface
et i est lintensit utile de la pluie (intensit climinue de linfiltration et des
autres pertes).
A lquilibre on aura : q = q ;do le profil applicable si x nest pas trs grand est donn par :

Si un coulement turbulent se fait sur une surface non plane telle quune route ; 1 quation
diffrentiellele long du profil en travers peut scrire de la maaire suivante :

avec k = O si la quantit de mouvement des gouttes de pluie est ngligeable

et k = 1 dans le cas contraire.
et on peut lintgrer numriquement.

Pour un coulement laminaire sur une route o la vitesse et la profondeur deau sont
invariables sur le profil en long, Iwagaki cit par Van Te Chow a 6abli lquation
diffrentielle suivante quon peut intgrer numriquement ou par la mthode des isoclines :

avec

=q=(:)
-K,
3

n-i

isi2xy

3ix

5 u2

Le profil en travers de la route est reprsent par (figure A-4) :

);(

y=-H

Iwagaki en a tir les conclusions suivantes :

102

- le profil de la surface libre est indpendante de la pente du profil en long.

dY = O sauf au sommet de la route ( y = cte).
- dans le cas du ruissellement naturel, -

dx
y = cste en thorie si le profil en travers est parabolique (n = 2) dy/dx > O si n = 1 (la
profondeur deau augmente du sommet la base).

- leffet de la pente longitudinal sur la vitesse et la vitesse de frottement est daugmenter

les 2. Cet effet est plus grand vers la base pour les petites valeurs de n.

- leffet de la pente longitudinal est pratiqueheht ngligeable si celles-ci est petite (<2 A)
dans la condition normale de UL = 0,02.

- pour minimiser lrosion due aux gouttes de pluie sur une route non revtue, la pente
longitudinale devrait tre aussi faible que possible.

pour maintenir une rosion uniforme sur une route non revtue, une section en travers
avec n = 1 est la plus indique.

....* ....... .;
cr

...........,......*,..

m...

U
.
,
.W
*
.

,. ................

. . . ...
...iap, .,.

Fimire A-4 : Ruissellement latral sur une route

A.2 Mthode des points singuliers pour les canaux non prismatiques avec un
dbit constant &=O)

Lquatiqn diffrentielle (A-1) se met sous la, forme suivante avec une dtermination des
sections de contrle qui se fait par la mthode de recherche des points singuliers :

dY

dx

103

Lquation F1=0 (Q=Qn ou y-y, )donnera ou le profil normal fictif ou profil cimi-normal.
Lquation F p O (Q = Qc ou y = yc) donnera le profil critiauG.

Dans les canaux prismatiques Fi O et F2 O sont deux droites paraliltss entre elles et
parallles au fond du canai. Dans les canaux non prismatiques, ces 2 lignes peuvent se
recouper en un point appel point sinrmlier C et lquation diffrentielle prend une forme
indtermine du type /O.
-j

3=

= zc J, on aura y,, = yc qui dame la lime carackristiaue ou lime de

K
Z
transition (Fl/F2 = 1). Cete ligne passe par le point singulier C car Qn= Qc = Q (ou y,, = y, =
y) au pint singulier.

Lorsque Qn= Qc (

Les proprites suivantes peuvent tre notes pour cette ligne caradristique au ligne de
transition ;

w3

un ecoufement traversant ce profil en un autre point que te point singulier doit y 2tre

ecouiement 21 travers t: protiI de transition au paint singulier changera de rgime

(torrentiel fluvial ou vice versa)

horizontal

(2

=I)

sa posirion est dtemine:par feu caractristiques du chenal (@omPtrk, rugositc! et pente)

comme le montre It5quation Q, = Qc. Mais contrairement aux profils quasi-normal et
aiXique, elle est. indpendante du changement du dbit Q.

La pente de Ili ligne deau au p i n i singulier est gaie i la limite, si die existe de la leve de
tindtermination O/ :

i
iSn peut distinguer 4 types de points siriguliers (figure A-5) selon les racines de It?quatioa du

si lquation a deux racine relles de mme signe, on obtient un nud qui est stable
quand elles sont ngatives et instable si dles sont positives ;
0

si les 2 racines sont relles et de signes contraires, on a une & ;

si les 2 racines sont des nombres complexes avec une partie relle non nulle , on a un
fover qui est stable quand cette partie relle est ngative et instable quand elle est
positive
si les 2 racines sont des imaginaires purs, on a un centre

Figure A-5 : Courbes de remous autour dun point singulier. P=point singulier ; F=courbe de
remous passant par P et asymptotique Fl=O ; F=courbe de remous passant par
P mais non asymptotique F1=0 ; f=autres courbes de remous ; FI=ligne de
lcoulement quasi-normal ; et T= ligne caractristique ou ligne de transition.

B METHODES
NUMERIQUES DE RESOLUTION DES EQUATIONS DIFFERENTIELLES DU PREMIER
ORDRE A CONDITION INITIALE

Lquation diffrentielle de la forme (A-2) admet une solution exacte pour toute fonction F
dfinie et continue sur lintervalle [xo, XJ

105

Une mthode simple comme la mthode de la tangente fournit une fonction qui converge vers
cette solution exacte si la condition de Lipschitz (A-3) est satisfaite.

On peut utiliser les mthodes pas spars ou les mthodes pas lis . Dans un algorithme
pas spars, on passe de Yi yi+1 en tenant compte seulement de linformation connue au
point Xi, tandis que dans les algorithmes pas lis, on tient compte des valeurs antrieures au
point xi mais on aura des problmes de dmarrage, dinstabilit et des difficults de
changement de pas h = Ax.

B. 1 Mthode de Runge-Kutter dordm p et de rang q.

Partant dun point initial yi que lon note yi,, on calcule par la tangente un point intermdiaire
dabscisse xi1 = xi0 + 91h
avec xi0 = xi ;do :

On peut calculer la tangente en ce point intermdiaire Yi1 :

On construit un nouveau point intermdiaire yi2 dabscisse XQ

combinaison linaire ay;+Py! des 2 tangentes y,oet y;l :

= Xi0

+ Qz h 1 aide dune

Yi2 = YiO +g2h[.F[Xi0

Yi0)+P F(nil

Yil)]

et cela pour (4-1) points intermdiaires.

La formule qui donne le point yi+l considr comme le qkmepoint intermdiaire est :
a-1

avec Fik = F (x

+ Qk h ,Yi + h

k-1

C A k i Fil) qui est une formule

l=o

de rcurrence.

106

On dtermine les coefficients Akl et $k (k variant de 1 q ;1variant de 0 a 9-1, ce qui donne q

valeurs de 8 k et q2 valeur de Akl ) de sorte que la formule obtenue pour Yi+] concide avec le
dveloppement en srie de Taylor jusqu lordre D le plus lev possible.
6.1.1 Exemples :

a). Mthode dordre 1 et de rang 1

avec A

10

=1

b). Mthodes dordre 2 et de rang 2

Si lon dveloppe la quantit F (xi ,y. ) ,on obtient :

1

F(x. , y i )
1

F[xi+91h,yi+ hA,,F(xi,yi)

avec d x = & h
dY = h Al0 F (Xi, Yi)
En reportant cette valeur dans yi2 on obtient :

Le dveloppement de Taylor donne lui aussi

Par identification de ces 2 dveloppements on trouve :

A20+A21= 1

pour les termes en h F(Xi, Yi)

107

dF
3Y

pour les termes en h F--A

Do

21

a1 = 21-

a-F
pour les termes en h2 aX

91 = Al0 et 3 quations 4 inconnues

1
alors
2
de la tangente amliore.

a).Si 91= Al0 =

A,, = 1 ; A,, = O et on obtient la mthode

p). Si 9, = Al0 = 1 alors

A,, =

; A,, =

1
et on obtient la mthode de Euler
2

Cauchy

Y). 91 = Al0

1
3

cest la mthode de Heun.

c) Mthodes dordre 4 de rang 4

Par la mthode didentification prcdente, on obtiendrait un nouveau systme indtermin,

par exemple quant au choix des Si, do une infint de mthodes. On utilise gnralement
1
9 1 = 9 2 = - et 93=1,cequicorrespond:
2
1
yj+i = Y i + -- [ kl + 2 k 2 + 2 k3 +k4]
6

h
x. +- , y. +1 2 1 2

On rduit le pas h ds que

> 10-, pour minimiser les erreurs (diffrence entre la

solution exacte et la solution approche).

Gill a propos une formule damlioration de la mthode Runge-Kutta dordre 4 qui permet
dobtenir une plus grande prcision par compensation au moins partielle, des erreurs darrondi
faites chaque pas. On prend qo nul au dpart ;sil ny a pas derreur darrondi, q4 serait nul.

108

k3=hF(xO+-,
h
2 y2)

Y2=Y1+ (I-@k2

-ql)

t-

c13=c12+3[(1+g)(k3 - q 2 ) l - [ l - E ) k 3

En pratique, q4 vaut approximativement 3 fois sur 114 due aux erreurs de chute.

Pour compenser, on prendra pour qo au pas suivant cette valeur de 94 :

6.2 Mthodes ii pas lis

109

- Po = O qui donne des algorithmes explicites c'est dire n'utilisant que l'information
antrieur k, on ajoute alors un terme

Pk-r-lFk-r-l
pour avoir 2r+linconnues.

- Po f O

qui donne des algorithmes implicites oh l'on ne peut regrouper les termes en yk
que si F(x,y) est linaire c'est dire de la forme F(x,y) = A(x) . y + B(x).

On dtermine les coefficients ctj et pj par deux formulations possibles :

Formules F rs
On les dtermine de sorte que la relation soit exacte pour toute quation diffrentielle
dont la solution analytique y(x) est un polynme de degr 5 s. Le rang r est le nombre de
points antrieurs utiliss. L'ordre s et le degr maximum du polynme tel que si la
solution est un tel polynme, la solution approche fournie en la solution exacte.
Cela revient dire que la formule est exacte quand on l'applique chaque monme de
degr infrieur ou gal s.
xk = xo + kh ,correspond la solution approche yk

la solution exacte est y&)

degr infrieur ou gal s

=@

(xk) o @ est un polynme de

et on doit avoir yk = y(Xk)

1

On obtient ainsi le systme de s quations 2 r + 1 inconnues pour les monmes x

Pour 1=0 : 1=a ,+a,+....+a

(r+l}'=a]r'+a2(r-ly+a&-2j+. ...+a
,+p,zr'-'+p,Z(r-l~ - l + ~ ~ Z ( r - 2 ~ - ] + . . * + ~ ~
qu'on applique jusqu' au degr s. La plus grande valeur de s est gale 2r.
Exemple : la formule explicite de Adams de rang 2 et de degr 2 .
a2=0 et Po=O ce qui donne :

110

Mthode des olvnmes d 'intemolation

Elle consiste crire la formule exacte suivante :

et par anaiogie l'appliquer la formule approche car yk est diffrent de y(Xk) :

On ne connat pas la primitive de F(x,y) priori mais on connat sa valeur pour tous les
points Xi O 5 i Ik, et on peut remplacer la fonction par son polynme d'interpolation
9m(x)sur m points appartenant l'intervalle [k-r, k+l].
On interpole par exemple avec les m+l points yk-1 ,yk-1.1 , ......., yk-i avec 1-i = m > O
comme le montre la figure A-6 :

Avec Fk-j = F(Xk-j, yk-j) et

En portant ces relations dans l'quation (A-2)' on trouve :

Exemde
Dans le cas o r = 2,1= 3 et i = O ;on obtient la solution suivante :

avec Ap21/8 ;A1=-9/8 ;A2=15/8 et A3=-3/8

111

Toutes les combinaisons r, i, 1 sont possibles avec r 2 O ; i 2 -1 et 1 > i ;toutefois on

prend en gnral :

i =O

avec r=O qui donne la formule d'Adams explicite

avec r=l qui donne la formule de Nystr6m explicite
avec r=O qui donne la formule d'Adams implicite
avec r=l qui donne la formule de Nystrom implicite

-9

integiration

yk+l

Figure A-6 : Mthode des polynmes d'interpolation pour les algorithmes pas lis.

B.2.1 Exemples classiques

Adams exdicites

Adams imdicites

il2

Ces formules explicites ont une erreur dordre 1 + 2 tandis que les formules implicites
ont une erreur dordre 1 + 3

6.2.2 Mise en uvm des mbthodes implicites

J=r

Lquation de dpart est yk+i=yk+yk-,+

...+yk-,+XpiFr-,-,
1-0

1) On dmarre lalgorithme avec une mthode pas spars p u r calculer yo, yi,
...,yr-1 ;do les valeurs de FoyFi, ... Fr-1

2) Pour la mise en uvre, on peut considrer 2 grandes mthodes :

- mthode itraive
Soit yk+l(O) une premire estimation de la solution, on calcule
successivement les termes de la suite :

jusqu la convergence. On peut prendre pour la valeur initialesuivante

pour dmarrer les itrations :
yk+l(O) = yk 4- h F(Xk yk)
Le problme de la convergence est satisfait quand les conditions suivantes
sont respecte (K est une borne donne) :
9

<K

et h+- 1

Po4

- mthode du vrdicteur - correcteur

Elle consiste trouver une premire valeur de yk+l* par une mthode
explicite suffisamment prcise (cest la formule de prdiction).
Puis utiliser une formule de correction p u r trouver yk+l

113

Les deux formules tant du mme ordre.

Par exemples :

114

Chapitre 4

Ecoulements brusquement varis

CHAPITRE 4
ECOULEMENTS BRUSQUEMENT VARIES

! GENERALITES
ET METHODES DETUDE
2 LE RESSAUT HYDRAULIQUE
2.1 DEFINITIONS
2.2 IMPULSION TOTALE

2.3 POSITION DU RESSAUT

2.4 CALCUL DU RESSAUT EN CANAL RECTANGULAIRE

3 LA SECTION DE CONTROLE
3.1 DEFINITION

3.2 APPLICATION
4 LES DEVERSOIRS
4.1 DEFINITIONS
4.2 DEVERSOIRSA SEUIL MlNCE

4.3 DEVERSOIRSA SEUIL EPAIS

4.4 DEVERSOIRSPROFILES
4-5 DEVERSOIRS
NOYES
5 LES VANNES

5.1 VANNES DENOYEES

5.2 VANNES NOYEES

6 LE VENTURI ou PARSHALL
7 LES APPLICATIONS DU RESSAUT HYDRAULIQUE
7.1 BASSINSDE DISSIPATION DENERGIE

7.2 OUVRAGES DE GESTION DE LEAU ET COUPURE HYDRAULIQUE

CHAPITRE 4

ECOULEMENTS BRUSQUEMENT VARIES

1 GENERALITES
ET METHODES D'ETUDE

Un coulement brusquement vari est un coulement permanent dans le temps mais les
variables de l'coulement varient trs vite, voire de maniCre discontinue dans l'espace.
f . 1 Caractristiques

Les caractristiques d'un coulernent brusquement vari sont les suivantes :

courbure trs- prononce des lignes de courant , d'o la rpartition_des

pressions de la scction n'est plus hydrostatique
rpartition trs irrglire des vitesses (le coefficient d'kncrgie cintique cc
>>1) et existence r>arfoisde courants de retour
I'eIet du Frottement sur les Darois peut tre nkgliaeable (courte distance)

la surface libre est souveiit instable et irrgulire.

1.2 Principe d'tude des couiement rapidement varis
D6msle cas gnral, on ne recherche pas i tracer Sa surface libre cause de ce que a t au
paridgraphe 1.1. On diermine deux sections yui englobent au plus prs l'coulement et o
les rpartitions de pression et de vitesse y sont rgulires.

Entre ces deux seciions dfinissant l'coulemcnt, on appliqiie :

soit le thorme des uuantits de mouvement lorsqu'ii s'agit d'coulement

divergent avec forte dissipation d'nergie : ressaut par exemple.
0

soit le thorkme de Bernoulli s'il s'agit d'coulement convergent sans

dissipation d'nergig : coulernent sous une vanne ~ aexenide
r

117

2 LE RESSAUT HYDRAULIQUE
2.1 Dfinitions
Le rcssaiit est unc surelvation brusque de la surfacc libre qui sc produit lors du passage du
rgime torrentiel au rkgime fluvial (figure 4- I ). L'nergie cintique de l'coulement torrentiel
est transforme en nergie potentielle avec une forte dissipation et une grande agitation de
l 'kcoulenient.

Lcs profonciirs avant (y!) ct aprbs (yzj lc ressaut sont appclkcs profondeurs conjugues au
sens du ressaut .

O n appeIle hauteur du ressaut la quantite y2 - y1 el la h g u e u r du ressaut, la distance entre

les sections 1 et 2. Cette longueur ne peut eh-e dtermine qu'exprimentalement en fonction
des caractcristiqucs dc la foimc de ia section et des caractristiques de l'coulement.

lx rendement du ressaut est ie rapport de l'augmentation nergie potentiel (y1 -y]) sur 1
diminiiticm d'nergie cintique (

u:

L
T
S
-- ) quand on passe de l'coulement torrentiel h
2g 2,q
-.

'

I'Gcuuleiizcnt fluvid. Elle peut 2tre tire directement de la figure 4-5 ci-dessous.

Rpartltlon des pressions hydrostatlque

Figure 4-1 : Dfinition du ressaut

118

Certains auteurs donnent une dfinition classique du rendement du ressaut comme tant
lnergie de la section 2 rapporte lnergie de la section 1 :
I

On distingue, selon la valeur du nombre de Froude lamont Frl, diffrents types de ressaut
quon retrouve la figure 4-2.

Dans certaines circonstances, le ressaut est submerg ou noy, sil ne peut pas se mouvoir
vers lamont cause dun obstacle (vanne ou dversoir), lorsque la profondeur deau en aval
(contrle par laval) est suprieure la profondeur conjugue y2 de la profondeur deau en
amont y1 contrl en amont figure (4-3).Le degr de submersion est dfini par S

Y- Y2
Y2

o y est la profondeur deau relle en aval.

Figure 4-3 : ressaut noy ou submerg

2.2 Impulsion totale

Y

La projehion dans la direction de lcoulement du thorme des quantits de mouvement

appliqu entre les sections 1 et 2 o la rpartition des pressions est hydrostatique donne
lquation (4-1).

p--p--=F~~
Q Q
-FP +Mg sin 0 -Ff

SI

s
2

(4-1)

- les forces de frottement Ff sur la paroi sont ngligeables cause de la

courte distance

- la

composante du poids Mg sine est nulle (canal horizontal) ou

ngligeable (canal sub-horizontal).

119

Type de
ressaut

1 < F r < 1.3

1.3 < Fr < 1.6

Caractristiques

Ressaut
ondul

Des ondulations
sur une longueur
importante :
hauteur peu
perceptible.

Ressaut faible

Petits tourbillons
superficiels ;
hauteur faible

Mouvement de
balancement
terticai du systme
I,6

Fr < 2

Ressaut
oscillant

Ressaui bien
localis
Ressaut calme

'aques d'eau dans

tous les sens
Ressaut fort

Figure 4-2 : Definitions des types de ressaut (les valeurs du nombre de FROUDE sont donnes pour

120

- La composante des forces de pressions F, dans chaque section est &ale il

pg S yi; o S est laire de la section et yc; ia profondeur de son centre de

gravitk.

Dans ces conditions, lquation de la quantite de mouveinent peut se rduire lgalit des
quantits M dofinies par lequation (4-2) et appelCs impulsion totalc caliilEes pour Ics
sections 1 et 2 (.quation 4-3) :
(4-2 j

Le thorme des quantits dc mouvement apjdiquk au ressaut se traduit donc par la

conservation de limpulsion totafe M d6finie par lquation (4-2).
On peut donc par ldquatioii (4-31, qui est symtrique en y] et yz, calculer la profondeur
onju&e y2 si lon connat y1 et vice versa.

2.2.7. Variation de limpulsion totale dune section en fonction de la profondeur

deau pour un dbit fix.
Ltude de Ia fonction M(y) pour une section vase vers le haut montre deux branches

paraboliques : lune quand y tend vers zro (la fonciion se compte comme p - )

Si

et lautre

quand y tend vers linfini (la fonction se comporte comme pg S y ~ ) .

La fonction tant positive et continue dans cet intervalle atteint iicessairement un minimum
pour dM/dy = O. Ce minimum correspond la profondeur critique (figure 4-4) et, de la
drivation de M :

121

I-

F k g e 4-4 : Variation de l'impulsion totale M en fonction de la profondeur d'eau y pour un

dbit Q fix.

On retrouve l'quation du rgiine critique dilni par le nombre de E'roucie au carr est gal 21
I 'unit6 ;

La relation dS/dy=i a t dj tablie et pour

W Y , 1 S ,il suffit de constater la figure 4-5

que

* Syc est le moment statique de S par rapport h la surface libre.

ie moment statique de S + AS par rapport la surface libre es1 :

en aiculant la diffkrence des deux termes et en faisant tendre Ay vers zro on retrouve
1' galit.

122

&&.-iure-4-5 : Calcul de dS!dy= I et d(Syc,)idy=S

2.2.2 Dtermination graphique des profondeurs conjuguees.

La figure 3-3 montre que :
I ) Pour M suprieure Mc, il existe 2 fondeurs Y I et y1 qui correspondent 2i la rneme
impulsion totale M. Ce sont des profondeurs conjugues au sens du ressaut :
y{,avant le ressaut est torrentiel
0

yz, aprs le ressaut: est fliivial.

Daprs lquation (4-3), la connaissance de lune dentre elles permet de

dterminer graphiquement lautre une fois que la courbe est connue (forme de la
section et dbit connus),
2 ) Pour M

= h4 qui est minimum, on obtienl

la profondeur critique.

3) Pour M < Mc, il ny a pas dcoulement correspondant cette impulsion totale M

et le dbit Q.

2.2.3 Dtermination graphique de la perte de charge du ressaut.

La perte de charge du ressaut est la diffrence entre les charges Iiydrauliqucs la section
avant le ressaut (section) et celle aprrs le ressaut (section 2). Elle correspond galement a la
diffrence entre les charges spcifiques i-1s

=y

+---Q2 si le canal est horizontal ou suppos

gS

tel :

A I I = TI,! - l i s 2 .

En traant limpulsion totale M et la charge spcifique 111 en fonction de la profondeur sur le

mme graphique, on peut donc dterminer la perte de charge du ressaut connaissant une des
profondeurs y1 ou y2 (figure 4-6).

123

Hst

'

Figure 4-6 : Dtermination graphique de la perte de charge du ressaut.

2.3 Position du ressaut

On utilise la notion de profondeur conjugue pour dterminer A priori la position du ressaut si
l'on connait les lignes d'eau de l'coulement torrentiel avant le ressaut et de l'coulement
fluvial aprs le ressaut ; ces lignes d'eau tant contrbles par l'amont et pax l'aval
respectivement

Dans la pratique, on peut procder de la manire suivante (Figure 4-7) :

124

Tracer les lignes deau amont et aval (courbes de remous). Les extrmits du ressaut
doivent tre sur ces courbes leurs profondeurs deau doivent tre conjugues au sens
du ressaut.
rracer la courbe conjugue de Iuiie des courbes de remous et i lintersection avec
1 autre, on a la position du ressaut.
Tenir compte ventuellement de la longueur du ressaut dans la construction graphique.

Figure 4-7 : position du ressaut.

2.4 Calcul du ressaut en canal rectangulaire

Dans le cas dun canai rectangulaire de largeur b, lquation (4-3) peut se simplifier en notant
q ie dbit par unit de largeur avec y(; = y/2 :

en inuitipliant les membres de iquation ci-dessus par 2gylyz et en mettant au facteur (y? YI), la conservation de limpulsion toiale se rduit une quation symtrique du second degr
en Y1 et Y2 (4-4)

v,y: + Y z Y I

zql

---

-0

(44)

Le calcul des profondeurs conjugues peut alors se faire algbriquement :

y?=&(-1+$+8F:)
2
ou
YI

=<(-1+&8F2
Y

(4-5a)

(4-5b)

L.

25

(4-5c)
La perte de charge, compte tenu de l'quation de continuits Lhyi=U~yz=q et de
i 'quation (4-4)peut se simplifier par des Iransformations algbriques :

Le rendement du ressaut galement se simplifie avec les mmes considrations

La longueur du ressaut est determine de Faqon exprimentale et les rsultats, exprims soit
en fonction de y2 soit en fonction de (39 - y,) divergent (difficult de dterminer ia fin du
ressaut :
mLR=6(1.8(~2-yi);
0

Figure 4-8 o la longueur du ressaut L est dom& en fonction de y2 et le nombre de

Froude Fr I ;
LR =: (4.9S -t- 6, y2 ou S est le degr de submersion pour un ressaut submerg.

*LR=5y2 1+4

/y:-]

- pour un canal trapzodal o L est la largeur en gueule.

Exenmie

Dans un canal rectangulaire de 5 in de large, s'coule wi dbit de 10 m3/s avec une

profondeur d'eau yi = 0,20 m.Dterminer la profondeur conjugue y2 et le rendement du
ressaut.
Rvonse yz=l,92m et RR1=34%

126

Figure 4-8 : Longueur du ressaut en canal rectangulaire en fonction du nombre de Froude aval
(D'aprs V.T. CHOW)

3 LA SECTION DE CONTROLE
3.1 Dfinition

On appelle section de contrle toute singularit (fort accroissement de la pente, rtrcissement

important, seuil dversant dnoy, etc..) qui provoque une augmentation suffisante de la
vitesse de l'eau pour faire passer l'coulement du rgime fluvial au rgime torrentiel.
Il existe donc dans l'emprise du "contrle" une section o le rgime est critique; d'o la
valeur de la profondeur critique yc ( et par suite la cote correspondante H du plan d'eau dans la
section de contrle) est entirement dtermine pour chaque dbit Q par le profil en travers de
la dite section.
Le contrle d'un coulement fluvial est toujours situ l'aval de cet coulement.
Celui d'un coulement torrentiel se trouve l'amont de cet coulement.
3.1 Applications

En fait l'coulement dans un canal quelconque est essentiellement dtermin par les
mcanismes de ((contrles )) qu'on y trouve. La notion de contrle (tout appareil qui

127

dtermine une relation biunivoque hauteur-dbit) est par consquent capital dans ltude des
coulements surface libre et il sera approfondie dans ltude des singularits (chapitre 5 )
Une des applications les plus courantes est la mesure de dbit dans les cours deau. On
sefforce de placer des stations de mesure de dbit dans une section telle qu une cote du plan
deau h corresponde un seul et mme dbit.
La courbe Q = f(h) est appele courbe de tarage. La courbe de tarage doit dors tre univoque
cela ne peut tre obtenu que dans certaines conditions : Lcoulement en amont de la section
de contrle ne doit pas tre influenc par laval. Dans le cas contraire , on est oblig de
prendre deux chelles.

4 LES DEVERSOIRS
4. i Dfinitions

On peut considrer le dversoir comme un orifice incomplet (figure 4-9). La hauteur deau
(nappe dversante) hl au dessus du dversoir est mesure loin en amont (loin de la zone
dapproche).
Si lpaisseur e de la crte (partie du seuil qui touche leau) est faible par rapport la hauteur
deau hl (hl > 2 e), on parle de dversoir seuil mince. Autrement (hl < 1,5 e), cest un
dversoir seuil pais. Entre les 2 limites, les deux modes dcoulement peuvent se produire.

F i m e 4-9 : Dfinition des termes du dversoir

La pelle P du dversoir est la hauteur qui spare la crte et le radier (fond du canal) et elle
dtermine la vitesse dapproche de leau qui devient ngligeable si P est importante.

128

Les formules de dbit sont dtermines selon les conditions de fonctionnement particulires
du dversoir qui sont :
Si hl est faible, cela empche le passage de l'air en-dessous de la nappe et on parle de
nappe adhrente (Figure 4- 1Oa)
Lorsque h augmente, la nappe tend se sparer de la paroi ; cependant si l'air n'entre
pas en quantit suffisante, une zone instable en dpression se forme sous la nappe que
l'on appelle alors nappe dprime(Figure 4- 1Ob).
0

La nappe est libre lorsque l'air peut circuler facilement sous la veine et que l'air
entrah par l'coulement est continuellementremplac (Figure 4-1Oc).
Lorsque pour une raison quelconque l'air n'est pas remplac, on aura une lvation du
niveau de l'eau sous la nappe dversante, qui s'appelle alors noye en dessous (Figure
4- 1Od)
Lorsqu'on augmente le niveau en aval h2 (h2 > 0,75 0,82 hl pour les dversoirs
seuil pais) ,la nappe devient noye (Figure 4-10e).

4.2 Dversoirs seuil mince

4.2.1 Approche theoffque

La figure 4-11 montre la coupe longitudinale d'un coulement au-dessus d'un dversoir
seuil mince et nappe libre sans contraction latral (coulement bidimensionnel)
Une analyse lmentaire peut tre effectue en supposant que la nappe dversante ne se
contracte pas verticalement quand elle passe au dessus du seuil et la pression travers la
section AB est gale la pression atmosphrique.
En supposant en plus qu'il n'y a pas de perte de charge (fluide parfait), on peut tablir par le
thorme de Bernoulli que la vitesse en tout point M de la section AB est gale
o y
est mesure partir de la ligne d'nergie et non de la surface libre en amont (la diffrence
entre les 2 tant lie la vitesse d'approche VO).
2

Le dbit par unit de largeur q est donn par l'intgrale de -&

2g

h + h (valeurs de y en A et
2g

B respectivement) de cette vitesse thorique :

129

a) nappe adhrente

.-

b) nappe dprime

.
.... .

c) nappe libre

d) nappe noye en dessous

air non ionouvol

e) dversoir noy

Figures 4-1O : Diffrentes configurations de la nappe

--

ligne dnergie

... .

rhpartition des preslons

i-

Figure 4-1 1 : Coupe longitudinal dun coulement au-dessus dun dversoir paroi mince.

Leffet de la contraction verticale de lcoulement peut tre exprim par un coefficient de

contraction C, donnant lquation du dbit par unit de largeur :

-I

Lquation (4-6) peut tre rarrange sous la forme (4-7) suivante o cd est un Coefficient de
dbit (quation 4-8).
-

q=C,& h

[(

(4-7)

J 3

2gh

( >)1]
vo

va
2 )

C=2Cc,.1+---

- -

2gh

(4-8)

On sattend ce que le coefficient de contraction C, aussi bien que le t m e exprimant la

2

vitesse dapproche K d p e n d e n t uniquement de la gomtrie des fiontires de

2gh

lcoulement, en particulier du rapport de la hauteur de la nappe dkversante la pelle h/P. Il

sen suit que le coefficient de dbit dpendrait de h/P seulement.
Plusieurs auteurs ont tudi et donn des formules de ce coefficient de dbit qui sera
directement inclus dans les formules empiriques du paragraphe 32 . 1 .
REHBOCK a donn la formule suivante, valable pour 5 < WP < 10 :

cd=2(0,6
3
1 1+0,0&)
P

131

Dans le cas o P=O ( ou h/Ptend vers de grandes valeurs), on obtient une chute et la hauteur h
serait gale la profondeur critique yc et le dbit par unt6 de largeur serait le dbit critique

&y! . est le cas lorsque h/P>20 selon Henderson o lon peut alors tablir :

Tout autre type de dversoir implique un coulement tridimensionnel ( contractions dans les
plans horizontal et verticai). Cest le cas des dversoirs rectangulaires avec contraction
latrale tudis par FRANCIS, des dversoirs triangulaires de GOURLEY et CRIMP, des
dversoirs trapzodaux de CIPOLLETI) qui sont utiliss pour la mesure de dbit dans des
canaux suffisamment large pour que la contraction soit minimale.

Dans le cas gnral, des formules exprimentales sont disponibles mais elles ont t
talonnes dans des conditions bien particulires de gomtrie et de type dcoulement (nappe
libre et dversoir dnoy) quil faut veiller respecter
4.2.2 Formules empiriques
Les dversoirs ci-dessous servent mesurer le dbit dans les canaux o ils sont installs.
Les formules qui sont donnes ont t donnes dans des conditions dtaionnage qui sont
parfois prcises avec une nappe libre. Il faudra par consquent les respecter en particulier
veiller ce lair soit renouvel sous la nappe pour quelle soit libre.
1O) dversoir rectangulaire sans contraction

REHBOCK (19291
Cest le dversoir le plus utilis pour la mesure du dbit en ngligeant le terme 0,001 1
dans la formule suivante :

Q = (1,782 + 0,249) L (h + 0,001 l)3/2

BAZIN f1898)
Cest la formule la plus prcise car on a plus de donnes sur ce dversoir :

Sil y a une nappe dprime il faut corriger en hausse le dbit en fonction de la

dpression observde, Lencastre Tableau 152

132

La construction du dversoir de Bazin doit satisfaire les conditions suivantes (Figure

4-12a):
a) canai de parois verticales et bien lisses longueur du seuil et largeur du canal
b) la crte ne doit pas tre trop basse et le seuil est en mince paroi
(conformment figure ci-dessous).

c) la longueur du canal en amont doit tre au moins gale 20 h avec a l l e s (en

bois ou brique creuse) pour uniformiser la vitesse d'approche.
d) l'aration de la nappe doit tre complte pour maintenir une nappe libre
et on doit installer un tuyau de ventilation si ncessaire et contrler la pression
sous la nappe au moyen d'un manomtre.
e) la charge doit tre lue une distance au moins gale : 5 10 h.

-..ec

lOcm

Figure 4-12a : Dimensions du dversoir de BAZIN

20)Dversoir rectangulaire avec contraction latrale
O

FRANCIS
Francis a trouv exprimentalement que la contraction latrale chaque extrmit du
dversoir rectangulaire est gale 1/10 de la hauteur h ,pourvu que la longueur du
dversoir L soit au moins gale 3h.

la hauteur d'eau h doit tre mesure 2 m au moins l'amont du dversoir.

La sur-largeur Q du canal doit tre au moins gale 3h (Figure 4-12b)

133

Figure 4-12b : Dimensions du dversoir de FRANCIS.

3 y Dversoir triangulaire

GOURLEY et CRIMP
Cest lune des formules les plus employes
Q=1,32(tgT)h
a 2,47
La figure 4-12c dtermine les dimensions respecter pour la prcision de cette
formule.

3 L ).
Figure 4-12c : Dimensions du dversoir de GOURLEY et CRIMP (a 2 4
THOMPSON
Pour le cas particulier
Thompson

des=-,
2

on emploie aussi couramment la formule de

134

HEGLY
n linfluence de la vitesse dapproche ou moyen de
Hegy a m i s en vidence pour a= 2
la formule suivante, valable pour 0,l O m 5 h -< 0,50 m :

S = h2 = crise du dversoir limit au niveau correspondant la

charge h ;

S = section de lcoulement dans le canal

4 y Dversoir trapzodai

GOURLEY et CRIMP
Pour une valeur quelconque de langle a dinclinaison des joues du dversoir sur la
verticale (Figure 4-12d), le dbit set la somme des 2 termes suivants :

CIPOLETTI
Le dversoir de Cipoletti correspond au cas particulier o tga=1/4 (a=) et la formule
empirique est :
Q = 1,86.L.h3/2

Les dimensions doivent obir aux conditions suivantes conformment la figure 412d ci-dessous :

* 0,08m< h

0,6 m

*P>3h
*e> 2h
* L> 3h

135

Figure 4- 12e :Dversoir circulaire

Les dversoirs servent galement rguler les canaux en particulier en contrlant la hauteur
deau dans un bief du canal. Dans ces cas, on peut trouver des cas particuliers par rapport la
disposition du dversoir sur laxe de lcoulement ou par rapport sa construction :

I O) Dversoir proportionnel
Lintrt de ce type de dversoir est de raliser une sensibilit constante pou tous les
dbits (lerreur absolue AQ tant proportionnelle Ah) et toutes les charges.
Le Dbit est donn par la formule suivante :

Q=Kh
Le profil de la joue du dversoir est donn par lquation hyperbolique x = &te/
constante dpend du coefficient de dbit K du dversoir (Figure 4-120

& o la

Figure 4-12f :Forme du dversoir proportionnel

137

Figure 4- 12d : Dimensions du dversoir trapzodal

50) Dversoir circulaire
Il prsente les avantages de la simplicit, de la facilit dexcution et de mise en place et
la facilit de ventilation figure (4-12e).

HEGLY
Le dbit est donn par la formule ci-dessous o S est la section mouille et S est la
section du cercle.

60) Choix du type de dversoir de mesure (daprs Lencastre)

Les principes gnraux suivants doivent tre observs dans le choix du type :
a) pour avoir une nappe libre h 1 6 cm pour les dversoirs triangulaires, h 1 2 cm
pour les dversoirs rectangulaires
b) On doit viter les fortes charges h 5 60 cm
c) La longueur du dversoir rectangulaire doit tre 2 3h

Pour la prcision de certains types de dversoirs est donne au tableau 4- 1

Si Q < 30 Us ,on choisira un dversoir triangulaire
Si Q> 300 Vs, on choisira un dversoir Bazin
Si 40 -< Q 1.300 V s ,les deux donnent la mme prcision.

136

20)Dversoir oblique
La crte est dispose obliquement par rapport laxe de lcoulement dans lequel le
dversoirs est plac figure 4-12g.
Ce type de disposition permet de :
0

augmenter la longueur de dversement

diminuer la charge pour le mme dbit
permettre un rglage plus prcis du niveau deau

F i w e 4-12g : Vue en plan dun dversoir oblique

AICHEL

La formule de dbit est donne ci-aprs o Qn est le dbit dun dversoir de mme
largeur et de mme type dispos perpendiculairement laxe de lcoulement, et K est
un coefficient fonction de langle aigu de la crte avec laxe dcoulement.

@(l-25(+& h l

La formule est valable pour :

a) h/P<0,62 et B>30
ou
b) h/P<0,46 et B<30

20 25

30 35

47 59 70 84
5
5
0
0

40

45

50

55

60

65

70

102 125 151 179 227 298 412

O
5
5
5
O
O
5

Tableau 4-2 Valeurs du coefficient K en fonction de lorientation (daprs Carlier)

139

37 Dversoir inclin
Sa crte est perpendiculaire laxe de lcoulement mais le plan qui le contient est inclin
par rapport la verticale (figure 4- 12h).
Boussinesq donne un coefficient de correction KI (sur le dbit du mme dversoir dispos
verticalement) en fonction de linclinaison i en degr compte positivement si le
dversoir est inclin vers lamont et ngativement dans le cas contraire.

K I = 1 - 0,3902

1
180

+
~- .-

>v

coulement

Figure 4-12h : Dversoir inclin

4O) Dversoir latral

Amnag dans la paroi latrale dun canal, il sert dverser le trop plein et au rglage
prcis du niveau en en-tte du canal.
La forme de la ligne deau le long du dversoir a t dj tudie dans les coulements
dbit variable et le dbit est fonction du rgime dcoulement (torrentiel ou fluvial) dans
le canal (Figures 4-12i)
a

DOMINGUEZ
a) Rgime fluvial dans le canal (les conditions aval imposent la hauteur hl de la
figure 4-12i)
Le dbit est donn en fonction de hl par la formule ci-dessous o les coefficients
Cd et pl sont donns respectivement aux tableaux 4-4 et 4-5

b) Rgime torrentiel dans le canal (les conditions amont imposent la hauteur ho de

la figure 4- 129
Le dbit est donn en fonction de hl par la formule ci-dessous o les coefficients
Cd et cpl sont donns respectivement aux tableaux 4-4 et 4-5

140

1I

h moyen (m)
(hl+ho)/2

0,lO

Crte mince,
nappe libre

0,370

Crte paisse
et arrondie

0,315

Crte paisse
arrtes vives

0,270

Tableau 4-4 : Valeurs moyennes de

dversoir latral

Rgime
fluvial
ho/hl<l

Rgime
torrentiel
homl'l

0,360

0,355

0,350

0,320

0,320

0,325

0,270

0,273

0,275

c d

---H
0,276

0,280

dans les formules de DOMINGUEZ pour un

0,s

0,9

125 1,11 1

Valeurs 0,400
de
cpo ou cpl

3,856 0,924 1

Tableau 4-5 : Valeurs des coefficients cpo et <pl dans les formules de DOMINGUEZ pour les
dversoirs latraux en fonction de
O

ENGELS (pour un rgime fluvial dans le canal)

Le dbit est donn par la formule ci-dessous o b est la largeur du canal. :

La formule est valable si le parement aval du dversoir est inclin 45" et pour les
dimensions suivantes :

P>4 hl et 2P < L < 40P

O

BALMACEDA et GONZALES (pour un rgime fluvial dans le canai)

Le dbit est donn ci-dessous en fonction de la nature du seuil.

141

Q=133Lh,+O,O0304(

i)

Q=l ,43Lhi+,0017 1(

pour un seuil mince

pour un seuil pais

Figure 4-12i : Dversoir latral

4.3 Dversoirs seuil pais

4.3.1 Approche theOri9ue

Si un seuil est assez large (0,08IWeI0,5) pour maintenir une rpartition de pression
hydrostatique de lcoulement dessus , celui-ci sera critique et avec une vitesse uniforme, le
dbit est gal au dbit critique donn par lquation suivante o H est la charge au dessus du
dversoir tenant compte de la vitesse dapproche ( figure 4-13a).

Dans le cas o la vitesse dapproche nest pas ngligeable, des itrations permettront de
calculer aussi bien H que le dbit Q.

142

Figure 4-13a : Dfinition des termes dun dversoir seuil pais

Mais cette approche est une simplification de la ralit cause des hypothses faites
(rpartition de pression hydrostatique, vitesse uniforme, etc.). En effet si le seuil est court
(H/e2 0,5), lcoulement rapidement vari occupera toute lpaisseur e ; et sil est trop pais
(H/e I008) le frottement devient apprciable.
Afin de traduire ces effets secondaires, on introduit un Coefficient de dbit c d qui tient compte
des conditions damene, de la forme du seuil, etc. Les valeurs de cd sont donnes la figure
4-13b pour un dversoir rectangulaire seuil pais et au tableau 4-6 pour un dversoir
triangulaire plac dans un canal rectangulaire.
Dautre part on remplace H (la charge en amont) par la hauteur deau h correspondante en
introduisant un coefficient de vitesse C,,/(H-yJ/(h-yJ qui toujours suprieur 1. Les
valeurs de CV sont donnes la figure 4-13c pour diverses sections en fonction du rapport
cds*/so S est laire mouille de la section de mesure de h et S* serait laire mouille dans la
section de contrle si lon y constatait une hauteur deau gale h (figure 4-13a).
Enfin, la profondeur critique yc est exprime en fonction de h par la relation yc=Kh. La valeur
de K est donne par la rsolution de lquation h=yc+ymJ2o ymcest la profondeur moyenne
pour yyC(K=2/3 pour une forme rectangulaire, K=4/5 pour une forme triangulaire, voir aussi
tableau 1-2 du chapitre 1).

Dans ces conditions, le dbit est donn par la formule (4-9) ci-dessous.
(4-9)
Dans certaines formules on regroupe parfois les coefficients de dbit Cd et de vitesse CV et
ventuellement des constantes de la formule (4-9) dans un coefficient unique p appel
galement coefficient de dbit. Il faudra donc bien faire attention la signification de tous les
termes y figurant.

143

W
a0
0.7

Fimm 4-13b :Valeur de Cd pour un dversoir rectangulaire sans contraction latrale :a est la
pelle et b, l'paisseur du seuil (D'aprs Lencastre).

144

..
c
E

II

e
s.

[
c3

Plusieurs formules empiriques ont t proposes pour diverses formes et profils de dversoirs.
Elles peuvent tre utilises pour mesurer un dbit condition de vrifier les conditions
dtalonnage des coefficients empiriques.
Lencastre donne par exemple des formules assez spcifiques pour des profils trapzodaux
(Pelle, parements amont et aval fixs au tableau 200 page 599) et triangulaires (parements
amont et aval fixs au tableau 201 page 600)
4.3.2 Formules empiriques

BELANGER
Cest un dversoir rectangulaire sans contraction latrale avec un profil
rectangulaire (figure 4-14). Le dbit est donn en fonction de la charge H en amont
par lexpression suivante avec comme valeur approche p=0,385.

La hauteur deau sur le dversoir est telle que b=2H/3. On mesure donc h, et on
en dduit la charge H quon entre dans la formule pour calculer le dbit
Le tableau 4-7 donne des valeurs de p plus prcises en fonction de la charge H et
de lpaisseur du seuil e. On peut le rsumer ainsi :
1. Pour H/e 2 1 2, p est gal 0,414 (dversoir seuil mince) ;
2. Pour H = 0,45m, la nappe laval nest plus adhrente et p 0 , 3 2 8 ;
3. Pour H 4 O, 15m, les valeurs de p sont des valeurs moyennes

Tableau 4-7 : Valeurs de p dans la formule de Blanger.

146

BAZIN

Pour le mme type de dversoir, Bazin donne la formule suivanteo p est le

coefficient de dbit dun seuil mince :

avec K=(0,70+0,18~)

0,6O5+E)[i+O,5~(&#0,41
h

Pour ces 2 formules, le dbit augmente denviron 10% dans les cas suivants :
1. larte en amont est chanfieine 45 jusqu la hauteur 0,12 e ;
2. larte en amont est arrondie par un quart de cercle de rayon
r=0,05e

F i m e 4-14 : Dversoir de Blanger (P m 0,5m) et de Bazin ( h, est mesure et on

prend H=3hJ2

4.4 Dversoirs profils

Certains profils sont utiliss dans les ouvrages hydrauliques pour amliorer le coefficient de
dbit du dversoir : profil CREAGER utilis dans certains vacuateurs de crues des barrages,
profil NEYRTEC utilis dans le modules masque).

Profd CREAGER
Le seuil pouse (avec une lgre surlvation A) la forme de la nappe infrieure
dun dversoir seuil mince fonctionnant librement sous une charge de conception
Ho ou Hd (H est mesure par rapport la crte du Creager et tient compte de la
vitesse dapproche). Sous cette charge, il ny aura pratiquement ni dpression ni
surpression sur la face aval

147

La forme est ainsi dfinie par lquation suivante o les diffrents termes sont
dfinis la figure 4-15.

Y= 1

K Hot-1

Les diffrents coefficients K, t de lquation et les valeurs de raccordement au

parement amont qui peut tre inclin a, b, Ri et R2 sont donns au tableau 4-8 en
fonction du fruit du parement amont.
Parfois pour des raisons de stabilit surtout, le parement aval donn par lquation
ci-dessus se prolonge par une surface plane tangente. La localisation du point de
tangence peut influer sur la capacit dcoulement.

4
aut4

__

VY

Figure 4-15 : Dfintion du profil Creager

Tableau 4-7 : Valeurs des paramtres de forme du Profil Creager

Le dbit est donn par la formule ci-dessus o le coefficient de dbit est influenc
par plusieurs facteurs dont la hauteur de la pelle P,le rapport de la charge relle H

148

sur la charge de conception Ho ou Hd, linclinaison du parement amont et le

prolongement du parement aval.

1. De la figure 4-16a, on peut tirer la valeur de cdO=o,1246Co en fonction de

P/Ho pour un parement amont vertical et pour une charge gale la charge
de conception Ho.
2. Linfluence de la charge de fonctionnement H est traduite par la figure 416b qui donne la rapport c d ( pourla charge H) sur CdO (calcul en 1. pour
la charge Ho) en fonction de HMo. On peut y constater que

si H<Ho, le coefficient de dbit diminue cause de la surpression

sur la face aval du seuil
si H>Ho, le coeficient de dbit augmente cause de la dpression
sur cette face. Ch recherche cette augmentation mais la dpression
ne doit pas dpasser une certaine valeur pour viter le dcollement
de la nappe et les nuisances conscutives (phnomnes de
battement, cavitation, etc.). La valeur limite de H est H,,=l ,60H0.

3. leffet de linclinaison du parement amont est traduit par la figure 4-16c qui
donne le rapport CmICh (coefficient de dbit du parement inclin sur celui
du parement vertical calcul en 2.) en fonction de P/Ho.

Lencastre donne la figure 4-16d illustrant linfiuence dun prolongement du

parement aval donn par lquation parabolique par une surface plane
tangente. Les courbes sont paramtrs en fonction du point de tangence et
de Q/Qo o QOest l dbit pour une charge correspondant Ho.
4. Si lcoulement est noy, la formule de dbit peut tre conserve avec une
rduction en % du coeficient de dbit (cas dun dversoir parement
amont vertical) donne la figure 4-16e (courbes en trait plein). Les
branches verticales des courbes correspondent aux variations provoques
essentiellement par la cote du fond en aval. Dans les branches horizontales,
la diminution du coefficient de dbit est li surtout au degr dennoiement

hdm.

Le graphique est subdivis par des lignes en traits discontinus en 4 zones

correspondant aux types dcoulement indiqu en mdaillon :

Type 1 : la diminution de c d est Uniquement fonction de cote du

radier en aval o le rgime dcoulement est torrentiel
Type II : il se produit un ressaut et il ya une zone de transition
(type 1ou II)
Type III : le ressaut devient ondul, hd/H diminue mais la veine
liquide suit encore la forme du seuil

149

Type IV : cest un coulement nettement noy ; le phnomne est

trs instable ;la veine liquide ne suit plus la forme du seuil pour des
valeurs plus faibles de h&i. Cest dans cette zone que la
diminution de Cd est la plus marque.

Profil NEYRTEC
Le seuil normalis Neyrtec permet deffectuer des mesures de dbit dans des
canaux et il est galement utilis, dans les modules masque de ce mme
fabricant, comme dversoir.

Il assure un rgime dnoy pour une large gamme de dbit (talus avai) et pour des
faibles dniveles disponibles tout en vitant le dpt des matires ea suspensions
dans leau (taius amont).
Son profil est triangulaire avec le sommet suprieur arrondi suivant un rayon de
courbure r fonction de la pelle P (figure 4-17). Les valeurs du rayon de courbure
au sommet suprieur, des angles des sommets la base et de la largeur de celle-ci
sont donnes dans la figure.
Le dbit par unit de largeur q est donn par la formule classique du dversoir avec
les valeurs du coeficient de dbit c d qui sont exprimes en fonction du degr
dennoiement hz/hI dans le tableau 4-8.

h2h

cd

O
0,500

O,167
0,490

0,375
0,470

0,600
0,450

0,800
0,420

0,875
0,405

Tableau 4-8 : Coefficient de dbit cd du dversoir Neyrtec en fonction du degr dennoiement

150

Firrure 4-16a :Valeurs de CdO = 0,1246 CO(COlu sur la courbe). (d'aprs USBR)

Rapport charge de fonctionnement sur charge de conception H/Ho+

F i m 4-16b : Correction du coefficient de dbit pour une charge de fonctionnement H.
(d'aprs USBR)

151

Valeurs de P/Ho
Fime 4-16c : Correction du coefficient de dbit pour un parement amont hclh. (d'apds
USBR)

Fiam 4- 16d : Influence du prolongement de la face aval par un plan tangent au poht d sur le
coefficient de dbit (d'aprs Lencastre)

152

Position du radier aval (y+h)/H

Figure 4-16e : Rduction du coefficient de dbit c d lie aux contions aval (degr
dennoiement et cote du radier). (daprs USBR)

B=5,5P

Fimue 4- 17 : Profil Neyrtec et ses dimensions en fonction de la pelle P.

153

4.5 Dversoirs noys

Lorsque le niveau en amont (ou le dbit) est influenc par le niveau aval, on dit que le
dversoir est noy. La formule dun dversoir noy peut tre mise sous la forme classique
amis avec un coefficient de dbit Cd qui est le coefficient de dbit c d applicable un
dversoir de mme type et nappe libre par un facteur de rduction K qui est fonction du
rapport h2h1 (figure 4- 18) :

On peut constater sur le figure 4-18 que le niveau aval commence influencer notablement le
dbit (rduction de 5%) pour :
0

--=0,80
hz

4h0 , 4 0 pour un dversoir seuil mince.

hi

pour un dversoir seuil pais ;

hi

Les donnes sur le profil Creager et le profil Neyrtec confirment galement ce constat. Les
variations de K en fonction de h ~ hdpendent
l
donc du type de dversoir.

5 LES VANNES
Les vannes sont utilises pour contrler le dbit la sortie des rservoirs ou dans les canaux
dirrigation. Dans son fonctionnement la vanne peut tre dnoye ou noye (figure 4-19).
En coulement dnoy, le jet (de profondeur y2) est torrentiel et ouvert latmosphre. Il nest
pas submerg par la profondeur deau laval y3 car il se produira un ressaut qui ralise une
coupure hydraulique entre les 2 coulements (torrentiel et fluvial).
En coulement noy, la profondeur deau en aval y3, qui est fluviale et donc contrle par
laval, est plus grande que la profondeur conjugue (au sens du ressaut) de la profondeur du
jet YZ.

154

Fimue 4- 18 :Facteur de comtion K de la formule du dversoir noy pour divers profils (1 :

seuil pais ;2 :profil trapzodal,3 : seuil mince,).

T---=l

Rourut

b
/~///////////////////////////,

Vannodhno-

Vrnno no-

Fimue 4- 19 : Dfinition de la vanne dnoye ou noye.

155

5.1 vannes dnoyes

Lexprience montre que :
la profondeur deau y2 au niveau du jet est gale Cc*b o Cc est un coefficient de
contraction (proche de 0,61) et b est la leve de la vanne ;
la rpartition de pression est hydrostatique au niveau de la contraction et les vitesses y
sont uniformes.
Lcoulement entre les sections 1 et 2 tant convergent sans dissipation dnergie, on peut
crire le thorme de Bernoulli suivant o q est le dbit par unit de largeur :
2

4
: %Y,

2=y2+=

En tenant compte de la relation y2=CCb; une transformation algbrique permet de calculer le

dbit unitaire (quations 4-10).
(4-1Oa)
(4-1Ob)

Comme le coefficient de contraction C, dpend de la gomtrie de la hntire (b/yi), le

coefficient de dbit dpendra uniquement de b/yi.
Les diffrentes expriences sur les vannes verticales ont montr que C, est de lordre de 060
a 061 et si le rapport b/yl est ngligeable, le coefficient de dbit c d sera de lordre de 061
daprs lquation (4-lob).
Le coefficient de contraction sera plus grand que 061 si la vanne est incline comme dans le
cas dune vanne segment (Figure 4-2Oa). Pour ce type de vanne, le coefficient de dbit est
donn en fonction des rapports yl/r, ah, b/r figures 4-20b.
5.2 Vannes noyes

La profondeur y2 o lon observe une vitesse uniforme est produite par la vanne et la
profondeur y3 est due un contrle situ laval. Si y3 est plus grand que la profondeur
conjugue (au sens du ressaut) de y2, la vanne sera noye.

Le jet issu de la vanne sera surmont par une masse deau, fortement turbulente parfois, sans
mouvement net dans une direction donne. On peut donc considrer que cest une masse
deau stagnante.

a
YS

//,////,////
a) Vanne secteur

b) coefficient de dbit
Firnues 4-20: Vanne secteur et valeurs du coefficient de dbit de dbit en fonction des
rapports yllr, d r et blr

En premire approximation on suppose que les pertes de charge se produisent uniquement

entre les section 2 et 3. Dans ces conditions, lapplication du thorme de Bernoulli entre les
sections 1 et 2 donne :
2

Il faut bien noter y la place de y2 (la pression est due toute la colmne deau y et non le jet
lieu que sur y2) dans le cas de la vanne noye.

y2 mais la vitesse na

Entre les sections 2 et 3, comme lcoulement est divergent, on applique le thorme des
quantits de mouvement :

157

Dans la pratique, YI, y2=Ccb et y3 sont connus davance et on doit calculer le dbit par unit de
largeur q et ventuellement la profondeur y qui doit tre la solution de lquation du second
degr (4-11) obtenue en liminant q2/2g dans le syseme constitu par les quations de
Bernoulli et des quantits de mouvement. Seule la solution suprieure la profondeur
conjugue de y2 sera retenue.
Y2Y3
4 ( y 2 - y 3 )y

Y172
+Yi - Y2

YiY2 +
Y -y (; -y ; 4(Y* - Y,)

(4-1 1)

avec

Connaissant y et YI, le thorme de Bernoulli donne :

On peur galement mettre le dbit sous la fonne classique (4-1Oa) condition de prendre la
Valeur de c d donne par lquation (4-1Oc) ci-dessous.

(4-1Oa)

(4- 1Oc)

Les valeurs exprimentales de c d (Figure 4-21) ne sont pas loignes de celles issues de la
thbrie ci-dessus (trait discontinu sur la figure) pour les valeurs de y3/b=5 ,Cc=0,6 et y calcul
par lquation (4-11) pour y1 et y3 connus (erreur relative de 5%).

Finute 4-21 : Coefficient de dbit d'une vanne verticale.

5.3 Formules approches

PONCELET donne la formuie approche ci-dessous pour une vanne dnoye sans
contraction. La hauteur d'eau en amont H est mesure partir de l'axe du pertuis et non du
radier (figure 4-22).

avec
C=O,70 pour une vanne verticale ;
C=0,74 pour une vanne inclin 1W2V (60") ;
C=0,80 pour une vanne inclin 2W2V (45").
Si la vanne est noye, on remplace H par AH la diffrence des niveaux d'eau entre l'amont et
l'aval.

159

Figure 4-22 : Formule approche de Poncelet.

6 LE VENTURIou PARSHALL
Le canal Venturi est constitu par une rduction progressive de la section transversale
jusquau col (section rtrcie) qui est suivi dun largissement progressif (divergent) pour
retrouver la section initiale (figure 4-23) ,

Pour un radier horizontal, en supposant un rgime critique au niveau du col avec une
rpartition de pression hydrostatique, on retrouve le dbit critique :

Compte tenu des hypothses faites, on affecte des coefficients de dbit c d et de vitesse CVqui
transforme cette 6quation comme dans le cas de la thorie des dversoirs seuil pais :

avec K=2/3 pour un col recangulaire (voir les dversoirs pais).

Les valeurs exprimentales de c d sont donnes au tableau 4-9 pour des sections de contrle
rectangulaires et trapzodales et les valeurs de Cv la figure.

Tableau 4-9 : Valeurs des coefficients de dbit c d et de vitesse des canaux Venturi longs

160

Parshall a apport une amlioration cette conception au niveau du profil en long , de la

normalisation des dimensions et de la transition afin de dnoyer le Venturi mme aux forts
dbits (Figure 4-23).
La formule de dbit dans le systme d'units SI pour un coulement dnoye est donne cidessous o X est un exposant fonction de la largeur du col 12 (Tableau 4-10) et la hauteur
d'eau hl est mesure aux 2/3 ge la distance entre le seuil d'entre du convergent et l'entre du
col :

Les dimensions principales du Parshall sont, dans le systme S.I., donnes par les quations
SLliVanteS :

A=l, 1961,+0,479
B=0,491,+1,194

F i m 4-23 : Dimensions principales du canal Parshall : a) vue en plan ;b) profil longitudinal
du radier.

161

Largeur
du col
12 (m)
Exposant

02
1,506

06

1,548

08
1,560

10
1,569

20

26

1,598

1
,609

Tableau 4-10 : Valeurs de lexposant x dans la forIm.de de dbit du canal Parshail (daprs
Carlier)
Pour le dtails des dimensions et des abaques des a Parshall )) commerciaux, on peut se
rfrer lannexe du chapitre 4 tir de Lencastre.
En choisissant convenablement le Parshall (limite de lennoiement), on peut limiter la perte de
charge (surlvation du niveau amont aprs introduction de lappareil dans le canal) 0,20hi.
En effet, il introduit des pertes de charge de lordre de 0,08 0,25 hl et il prsente lavantage
de ne pas provoquer des dpts de matires en suspension dans leau comme dans le cas dun
dversoir.

7 LES APPLICATIONS DU RESSAUT HYDRAULIQUE

Les applications du ressaut sont nombreuses :
0

dissipation dnergie la sortie des vacuateurs de crues ou des vannes de fond pour
viter lrosion du talus aval ;
augmenter ou maintenir le dbit dune vanne ou dun dversoir en ralisant une
coupure hydraulique avec le niveau deau aval qui ninfluence plus le dbit ;
mixer des produits chimiques utiliss dans le traitement des eaux ou arer leau si
ncessaire ;
etc.

7.1 Bassins de dissipation dnergie

Dans quelques types de ressaut, la dissipation dnergie est apprciable. Les vitesses laval
(rgime fluvial) sont rduites par rapport celles lamont (r6gime torrentiel). Ces proprits
font utiliser le ressaut comme dissipateur dnergie pour viter les rosions laval dun
(t coursier de lvacuateur de crue dun barrage )) ou dune vanne de fond.
Pour stabiliser le ressaut au voisinage du pied de louvrage pour diffrents dbits possibles, on
approfondit habituellement la radier du canal ralisant ainsi un bassin de dissipation bien
stabilis.

162

On peut placer la fin du bassin un seuil form par une pente adverse avec ventuellement
des blocs de chute et des blocs chicanes ; formes qui intensifient la dissipation dnergie et
aident stabiliser le ressaut lintrieur de bassin (Figures 4-24).
Le ressaut devra tre lgrement submerg dans le calcul du bassin de dissipation avec un
degr de submersion S de lordre de 0,06 012.
Dans les figures 4-25, les profondeurs conjugues sont dsignes par hl et hz ; et la longueur
du ressaut en canai rectangulaire horizontal est note C et donne par la figure 4-8.
Bassin type II USBR (figure 4-24a)
Cest un basin ,avec blocs de chute et seuil dent, quon utilise pour des chutes
infrieures 65 m, des dbits par unit de largeur infrieurs 45 m3/s/ml et des
nombres de Froude lamont Fr1 infrieurs 43.
Bassin type III USBR (figure 4-24b)
Le seuil nest plus dent et on trouve des blocs damortissement ou blocs chicanes
entre la range des blocs de chute et le seuil. Ce dispositif rduit considrablement la
longueur du bassin. La figure 4-24c permet de dimensionner les blocs.
On utilise ce type de bassin pour U 1 4 8 m/s ,q 4 8 m3/s/ml (pour viter la cavitation
sur les blocs damortissement) et Fr1>4,5.
Le niveau minimal de leau laval est celui qui correspond 0,83h2 ou 0,83y2.
Bassin type IV USBR (figure 4-24d)
Cest un bassin avec dflecteurs et un seuil terminal continu, qui est indiqu dans le
cas du ressaut oscillant (2,5<Fr1<4,5). Le rle des dflecteurs est dattnuer les
ondulations de manire significative.

7.2 Ouvrages de gestion de leau et coupure hydraulique

On utilise le ressaut pour raliser une coupure hydraulique dans le fonctionnement des
appareils de gestion de leau sur un primtre irrigu: modules masques , partiteurs
proportionnels. Cest ainsi que ces appareils ne seront pas influencs par laval.
7.2.1 Module masque

Il maintient un dbit plus ou moins constant pour une large gamme de hauteur deau. La
constance du dbit est obtenue par des moyens entirement statiques : association dun seuil
profil ( Q W ~et~dun
~ ) masque fixe plac au dessus qui corrige leffet sur le dbit (Q-hln)
dune lvation du niveau amont.

163

Figure 4-24a : Bassin type II USBR

Fimire 4-24b : Bassin type III USBR

164

Figure 4-24c : Dimensionnement des blocs de chute pour bassin type IIi USBR.

Figure 4-24d : Bassin type IV USBR.

165

Le dbit ne varie pas non plus en fonction du niveau aval cause de la coupure hydraulique
ralise par le ressaut.
Les modules se prsentent sous la forme densembles mtalliques monoblocs quil suffit de
sceller bonne cote dans des maonneries en attente pour raliser des ouvrages de prise
perfectionns.
Ils se fabriquent dans 4 types de dimensions en profil diffrentes, caractriss par le dbit
nominal passant par unie de largeur :
Srie X :
SrieXX:
Srie L :
Srie C :

10 Y d d m
2Ol/ddm
50 Y d d m
100 Ysldm

La dsignation de la srie est assortie dun indice 1 ou 2 selon que les appareils sont du type
1 ou 2 masques (Figures 4-25).
Le tout est suivi dun nombre qui dsigne le dbit nominal de lappareil, toutes vannes
ouvertes.

Le dbit est fkactionnable par chelon de :

5 l/s dans la srie X qui comprend des vannettes de 5,10,15,30 Ys ;
1O l/s dans la srie XX ( vannettes de 1O, 20,30,60, et 90 Ys) ;
50 Ys dans la srie L ( vannettes de 50,100,200, et 400 Ys) ;
1O0 l/s dans la srie C (vannettes de 100,200,400,600 et 1O00 Ys)
7.2.2 Parteur proportionnel
Le partiteur est une prise spciale partageant le dbit arrivant en 2 (ou plusieurs) fractions
constantes. Les volets mobiles permettent de faire varier la proportion des dbits.
Un bon partiteur doit prsenter une section critique suivi dun ressaut afin que le partage ne
soit pas influenc par les conditions aval (figure 4-26).

166

Fiaire 4-25a : Schma de principe et courbe de fonctionnement d'un module 1 masque

F i m 4-25b : Schma de principe et courbe de fonctionnement d'un module 2 maques

167

Yigure 4-2% : Courbe de lonctionnement aes moames a 1 masque

si

I
s

i
I

j
1

Figure 4-25d : Courbe de fonctionnement des modules 2 masques

169

F i m 4-26 :Partiteu' proportiomel.

170

Chapitre 5
Etude des singularits

CHAPITRE 5
ETUDE DES SINGULARITES

1 POSITION DU PROBLEME

1.1 OBJECTIFS

1.2 DON~PIEES
1.3 TYPE DE PROBLEMES
1.4 METHODEDETUDE

2 CHANGEMENTS DE PENTE
2.1 AUGMENTATION DE LA PENTE
2.2 DIMINUTION DE LA PENTE

3 CHANGEMENTS DE NIVEAU DU RADIER

3.1 EXHAUSSEMENT DU RADIER

3.2 ABAI~SEMENT
DU RADIER
4 CHANGEMENTS DE SECTION

4.1 RETRECISSEMENTDE LA SECTION

4.2 ELARGISSEMENTDE LA SECTION
4.3 RETRECISSEMENTLOCALISE

5 PHENOMENES LOCAUX
5.1 PILES DE PONT
5.2 GRILLES
5.3 COUDES
O

6 EXTREMITESDES CANAUX
6.1 PASSAGE DUN RESERVOIR A UN CANAL

6.2 PASSAGE DUN CANAL A UN RESERVOIR

172

CHAPITRE 5
ETUDE DES SINGULARI'TES

1 POSITION DU PROBLEME

Les ouvrages hydrauliques sont constitus essentiellement de canaux prismatiques prvus

pour acheminer le dbit en rgime uniforme. Mais les contraintes dues au
fonctionnement des ouvrages empchent de respecter partout les conditions d'coulement
uniforme.
Ces conditions ne sont pas respectes autour des "singularits" de l'coulement (changement
de pente, changement de niveau du radier changement de section, extrmits, etc.).

Dans les coulements en char%, une singularit provoque un abaissement brusctue de la ligne
de charge ; au contraire, dans un coulement surface libre dans un canal sufisarnment
long, une singularit provoque une perturbation qui se ressent le long d'un tronon plus
ou moins long. Mais en dehors de ce tronon, le rgime demeure avec des caractristiques
nergtiques similaires celles qui se produiraient s'il n'y avait pas de singularit.
Ainsi l'expression de la perte de charge localise a une signification un peu diffrente de
celle qui correspondait une conduite en charge et sa dtermination est gnralement plus
dificile.
1.1 Objectifs
L'tude des singularits de l'coulement permet :
soit de prvoir le dimensionnement local des ouvrages de telle manire qu'ils
contiennent les perturbations mentionnes ci - haut
soit de modifier les dispositions des singularits de telle manire que les irrgularits
de l'coulement soient rduites et que le dimensionnement courant des ouvrages les
contiennent.

Parfois des cas apparemment similaires provoquent des conditions d'coulement et des
pertes de charge bien diffrentes. C'est pourquoi toutes les fois que la nature du problme
le justifie, il est recommand de faire des essais sur modle rduit.
1.2 Donnes

Les donnes thoriques de l'tude des singularits ont t dj vues dans les chapitres
prcdents :

173

'

L'tude qualitative des coulements graduellement varis qui permet de dfinir l'allure
des courbes de remous proximit des singularits.
L'tude quantitative des coulements graduellement varis qui permet de tracer les
liynes d'eau correspondantes. connaissant la section de contrle correspondante.
L'tude des coulements brusquement varis qui permet de dfinir auantitativement
certaines section de contrle.
L'tude des coulements uniformes qui permet de dfinir quantitativement certaines
autres sections de contrle.
1.3 Type de problmes
a

Les problmes rencontrs sont de deux types :

Soit le dbit est connu : Il s'agit alors de dterminer la ou les sections de contrle et
de dfinir la ligne d'eau dans la zone d'influence de la singularit.
Soit le dbit n'est pas connu : Dans ce cas, une condition de niveau (Rservoir par
exemple) permet de calculer le dbit et ensuite de tracer la ligne d'eau.

1.4 Mthode d'tude

L'tude des singularits est base dans la plupart des cas sur l'interprtation des courbes Hs (y)
correspondant la singularit.
Le tableau 5.1 synthtise les conclusions de l'tude qualitative des coulements graduellement
varis. On peut en tirer les remarques suivantes :
1) Connaissant la charge spcifique Hs dans une section donne, la courbe Hs(y) dfinit deux
profondeurs y1 et y2 correspondant au .rgime torrentiel et au rgime fluvial
respectivement.

II faut une charge spcifique minimale Hsc pour tablir un coulement dans une section
donne.
Certains types de courbe de remous permettent d'augmenter la charge spcifiaue Hs et
l'coulement devra s'adapter suivant de telles courbes pour avoir la charge spcifique
minimale Hsc au niveau de la singularit.

dHs > O
Ce sont les courbes F 1, T 1, T2,C 1 pour lesquelles dx
8

Le passane du raime fluvial au rgime torrentiel ncessite le passage par une section
critique (chute)

174

Type Profondeur!
de
d'eau
Courbc Y
Y'Yn>Yc

Yn>Y>Yc

F1

F2

F3

Allures

Courbes Hs

w.

Yn>Yc>Y

-t

T1

Y>Yc>Yn

----i;

T2

.a
El

T3

c3

H2

Y>Yc

H3

Yc>Y

A2

A3

Tableau 5-1 : les divers types de courbes de remous

175

5 ) Le passage du rgime torrentiel au rgime fluvial ncessite un ressaut. L'quation du

ressat permet de calculer y2 (profondeur cdnjugue de YI)connaissant YI.

6) Il n'y a pas de courbe de remous permettant de se raccorder une profondeur normale
torrentielle l'amont. Cette profondeur doit tre atteinte l'amont immdiat de la

singularit ou par l'intermdiaire d'un ressaut. L'coulement torrentiel est influenc par
l'amont jusqu' la singularit ou au ressaut.
7) Il n'y a pas de courbe de remous permettant de se raccorder une profondeur normale
fluviale l'aval. Cette profondeur doit donc tre atteinte l'aval immdiat de la sinaularit
ou par l'intermdiaire d'un ressaut. L'coulement fluvial est influenc par l'aval jusqu' la
singularit ou au ressaut.

La combinaison de ces diverses remarques permet de dfinir les singularits de l'coulement.

Les singularits suivantes seront tudies en fonction du caractre rapide ou lent des biefs
amont et aval de la singularit :
I

changement de pente
changement de radier
changement de section
extrmits de canaux
Ces divers types de singularits sont tudis en fonction du caractre rapide ou lent des
tronons de canaux amont et aval de la singularit.
Dans ces divers cas le critre de l'nergie spcifique ncessaire pou franchir la singularit
cre des variantes supplmentaires.

2 CHANGEMENTS DE PENTE

Considrbns un canal de section constante dont fa pente longitudinale varie brusquement.

Les biefs situs l'amont et l'aval du changement de pente seront supposs uniformes et de
grande longueur ; aux extrmits amont et aval du canal le rgime est donc uniforme.
Proposons de dterminer la forme mise par la ligne d'eau de part et d'autre du changement de
pente, le dbit du canal tant connu.
On commencera par prciser la classe de la ligne d'eau (F, 1, C, H,A) dans chaque portion du
canal en comparant les profondeurs normale yn et critique yc (yc ne dpend que du dbit Q et
yn varie en sens inverse de la pente 1.
Les conditions aux limites permettront ensuite de dterminer les tvpes de courbes de remous
et leurs contrles.

176

2.1 Augmentation de la pente

On distingue trois cas possibles comme le montre le tableau 5-2a :

l mcas :Le rgirne normal est fluvial dans les 2 biefs (NFpuis NF)

Le rnime normal fluvial aval vient iusau' la singularit ( puisqu'il n'y aura pas de ressaut et
pas de courbe de remous permettant de se raccorder un coulement normal fluvial aval).
L'action de la singularit se fait ressentir l'amont par une courbe de remous F?- dans le bief 1
dont le contrle est Y.2.
-

>'

22 Cas :Le rgime normal est torrentiel dans les 2 biefs (NTpuis NT)
Le rgime normal torrentiel amont vient jusau' la singularit ( puisqu'il n'y aura pas de
ressaut et pas de courbe de remous permettant de se raccorder un coulement normal
torrentiel amont). L'action de la singularit se fait ressentir l'aval par une courbe de remous
T2
- dans le bief 2 dont le contrle est va
32 Cas :Le rgime normal estjluvial 1 'amont et torrentiel l'aval (nFpuis N J
L'action de la singularit se fait sentir l'amont et l'aval de la singularit. Une courbe de
remous F2 permet le raccordement yn1 (contrle yc) dans le bief 1 et une autre courbe de
remous 7'2 permet le raccordement y& (contrle yc) dans le bief 2
I

2.2 Diminution de la pente

On distingue 4 cas possibles comme le montre le tableau 5-2b :

l mCas Le rgime normal est fluvial sur les 2 bit@

Le rgime normal fluvial aval s'tablit jusqu' la singularit (puisqu'il n'y aura pas de ressaut
et pas de courbe de remous permettant de se raccorder un coulement normal fluvial aval).
La singularit fera sentir son effet l'amont (bief 1) par une courbe de remous Fi dont le
contrle est yn2.
2"'

Cas :Le rgime normal est torrentiel sur les 2 biefs

Le rgimJe normal torrentiel amont s'tablit jusgu'a la singularit (puisqu'il n'y aura pas de
ressaut et pas de courbe de remous permettant de se raccorder un coulement normal
torrentiel amont). La singularit fera sentir son effet l'aval (bief 2) par une courbe de remous
T3 dont le contrle est Yn1.

177

Cas

Interprtation sur les courbes Hs(y)

Lignes d'eau

Ha

/ 0

NF

NT
4

/ I

-.- \..-

<

...,

Tableau 5-2a: Augmentation de la pente. yn1 et yn1 sont les profondeurs normales
respectives en amont et en aval ; yc ne varie pas avec la pente.

178

Cas

Interprtation sur les courbes

HsW

Lignes d'eau

$a

4
A
..

!i

NF

v,

A
..

crl

1 4c

NF

L
&
..

E
4

d
.y

4)
O

6
.O

NF

Tak eau 5-2b : Diminution de la pente. ynl et ynl sont les profondeurs normales respectives en
amont et en aval ;yc ne varie pas avec la pente.

179

3h et P

cas : Le rgime normal est torrentiel lamont et fluvial laval

La ligne deau doit obligatoirement passer du raime

torrentiel au rgime fluvial ; il y aura
donc un ressaut. La singularit fera sentir ses effets une trs faible distance lamont o il
est torrentiel et laval o il est fluvial. Le rgime ne sera plus graduellement varie mais
brusquement vari, ce que traduit le ressaut.
La position du ressaut dpendra des grandeurs relatives des profondeurs normales vgLyn2L
leurs con-iugues Vnl&
- --dans les deux biefs.
3&

ca~: Si ynl> yn2 ou siyn2> S>nl ,.Le ressaut se produira dans le bief aval (2)

Le rgime normal torrentiel amont arrive jusqu la singularit ; une courbe de

remous F3 ayant comme contrle yni augmentera la profondeur deau y jusqu
I
atteindre ya puis un ressaut ( ~ ~y2i ) se forme.

Cas : Si ym1<yn2ou siyn2<

ynl

,.Le ressaut se produira dans le bief amont (1)

Le rgime normal fluvial aval arrive jusqu la singularit ;une courbe de remous T1
ayant comme contrle y d diminuera la profondeur deau y vers lamont dans le bief
1jusqu atteindre yn1, puis le ressaut se forme pour diminuer Yn1 y,l.

CHANGEMENTS DE NIVEAU DU RADIER

On suppose que la section du canal et sa pente restent constantes. On supposera galement

que les Dertes de charge sinaulires sont ngligeables seuil suffisamment progressif et
convenauement profil).
!
,
Les biefs situs lamont et laval du changement de niveau seront supposs uniformes et
sufisamment longs ; aux extrmits amont et aval du canal, le rgime est uniforme. Le
problme est de dterminer la forme de la 1igne.deau de part et dautre du changement de
niveau du radier.
On commencera par dterminer la classe de la ligne deau ; le dbit tant connu. Les
conditions aux limites permettront ensuite de dterminer les types de courbes de remous.

3.1 Exhaussement du radier

La courbe Hs(y ) du canal aval est translat de + h vers les ordonnes Hs (figure 5-1). Les
profondeurs critique et normale sont les mmes dans les 2 biefs; elles sont seulement
translateh de +h vers le haut dans le bief 2. onkiistingue 3 cas possibles comme le montre le
tableau 5-3a.

180

AV

yc

Figure 5-1 : Exhaussement du radier de +h : la courbe du bief aval est translat vers le haut de
la quantit h.
i)

1 Cas : Le rgime normal est fluvial

Le rgime fluvial aval stablitjusqu la singularit (puisquil ny aura pas de ressaut et pas
de courbe de remous permettant de se raccorder4 un coulement normal fluvial aval). Une
courbe de remous F1 de contrle y& diminue la profondeur jusqu y,1 dans bief 1.

*cas

9
et

: Le rgime normal est torrentiel

Cas lnergie spcifiilue du rgime normal amont est suffisante pour

franchir la singularit :Hsl(y,J H s C yHscl+h
Le rgime torrentiel amont est conserv jusqu la singularit quil peut fianchir
(puisquil ny aura pas de ressaut et pas de courbe de remous permettant de se
raccorder un coulement normal fluvial aval). Une courbe de remous T2 de contrle
y2
tel que Hsl(y,J=Hs2(y$ diminuera la
profondeur jusqu y& dans le bief 2.
r
>*
,

3& Cas : 1nergie spcijique du rgime normal amont nest pas suffisante pour
franchir la singularit :H s & J< HsC2=Hscl+h
Il faut une courbe de remous pour augmenter lnergie spcifique Hsl(y,J ; mais il
ny a pas de courbe de remous permettant de se raccorder un coulement normal
torrentiel amont.
Un ressaut transforme y,, en Yni; une courbe de remous T1 de contrle y3 tel
Hsl(yJ)=HsC2 augmente la charge spcifique jusqu la valeur requise HsC2pour
franchir la singularit ; tout ceci dans le bief 1.
La singularit est franchie avec HsCz et une courbe de remous T2 de contrle yc2
diminue la profondeur jusqu y& dans le bief 2.
J

181

Lignes d'eau

Interprtation sur les courbes

wu
Hs

/"Ys1

Hscz

Hsci

Hs

Hs2

Hsc2

HS1(yJ> Hkc2

H~ci

Ha

HSCl

Tableau 5-3a : Exhaussement du radier de +h. Les profondeurs y , et

~ yc2 l'aval sont les
mmes qu'en amont ;elles sont seulement exhausses de +h.

182

3.2 Abaissement du radier

La courbe Hs(y) du canal aval est translat de -h vers les ordonnes Hs (figure 5-2). Les
profondeurs critique et normale sont les mmes dans les 2 biefs; elles sont seulement
translates de h vers le bas dans le bief 2.

Figure 5-2 : Abaissement du radier de +h : la courbe du bief aval est translat vers le bas de la
quantit h.
1" et 2""cas :Le rgime normal estfiuviai

Les 2 cas sont similaires et ils se distinguent seulement par le contrle de la courbe F2 dans le
bief 1. Le rgime fluvial aval s'tend jusqu' la singularit dans le bief 2 (puisqu'il n'y aura
pas de ressaut et pas de courbe de remous permettant de se raccorder un coulement normal
fluvial aval).
le' Cas :l'nergie spctjque du rgime normal aval est suprieure l'nergie
critique de la s e c t h amont :Hs2GyiS =HsI@,,) h >H&

Le contrle y2 de la courbe de remous F2 dans le bief 1 est tel que Hs2&J =HS&,,)
h=H&> calcul sur la courbe amont

2"""

Cas :l'nergie spcwque du rgime normal aval est infrieure C'nergie

critique de la s e c t h amont :Hs2U>ns =Hs~@,,)h < HsC1

Le contrle de la courbe de remous F2 dans le bief 1 est yc et une chute se produit

la singularit pour atteindre le rgime normal fluvial l'aval.

183

Lignes d'eau

Interprtation sur les courbes

&(Y)
Hsi

Ha

HSCI
Htk2

Y=

iisa(y4.c

Y2

Yn

'Y

HSCI

Hsi

Ha

Hscl
Hsu2
/

1 I I.,
I V !
)"l
I
I

Y2Y" Y,

+Y

Tableau 5-3b : Abaissement du radier de +h. Les profondeurs y& et yc2 a l'aval sont les
mmes qu'en amont ; elles sont seulement abaisses de +h.
Q

184

3"'

cas :Le rgime normal est torrentiel

Le rgime torrentiel amont dtend jusqu' la singularit (puisqu'il n'y aura pas de ressaut et
pas de courbe de remous permettant de se raccorder un coulement normal torrentiel
amont). Une courbe de remous T3 de contrle y2 s'tablit dans le bief 2. le contrle y2 est tel
que HsI(VJ=H s ~ ( V Hk(Vj-h
~=
U

CHANGEMENTS DE SECTION

Le canal a la mme pente avant et aprs le changement de section. Les pertes de charge
localises seront ngliges (examiner leurs importances aprs)
Les courbes Hs(y) sont dcales en Hs et y et seront traces en fonction des gomtries. Dans
les exemples qui suivent, la forme du canal ne change pas (seule la largeur b pour un canal
rectangulaire ou trapzodal change.
4.1 Rtrcissement de la section

La courbe H s est
~ au dessus de Hs1 et elle est dcale en y vers la droite (figure 5-3).

Figure 5-3 : Rtrcissement de la section.

185
a

On peut distinguer les cas typiques suivants (tableau 5-4a) :

le*Cas :Le rgime normal estfluvial avant et aprs le rtrcissement

Le rgime normal fluvial aval s'exerce jusqu' la singularit (puisqu'il n'y a de courbe de
remous permettant de se raccorder un coulement normal fluvial aval et il n'y a pas de
ressaut).
Une courbe de remous F1 de contrle y2 tel que Hsl(y2)
dans le bief Ide ynl y2.

= Hsz(yd) remontera

la ligne d'eau

2"" cas Z e rgime normal estpuviai avant et torrentiel aprs le rtrcissement

Une courbe de remous F1 de contrle y2 tel que Hsl(y2)
le bief Ide ynl y2.

= Hsc2 remontera

la ligne d'eau dans

Une courbe de remous T2 de contrle yc2 diminuera la profondeur d'eau de yc2 yd dans le
bief 2.
3h cas ; Le rgime normal est torrentiel avant etfluvial aprs ie rtrcissement.

L'coulement normal fluvial vient jusqu' la singularit (puisqu'il n'y a de courbe de remous
permettant de se raccorder un coulement normal fluvial aval) dans le bief 2.
Un ressaut raccorde l'coulement normal torrentiel amont une courbe de remous T1 dont le
contrle y2 est tel que Hsl(y2) = Hsz(yd) dans le bief 1.

4""" et SL" cas ; Le rgime normal est torrentiel avant et aprs ie rtrciwement.

dmCas :L'nergie spcifique de l'coulement

normal amont n'est pas sufliante pour

franchir la singularit :HsICynS < Hscl

Il faut qu'une courbe de remous augmente Hsl(Yn1) avant le franchissement ; ce qui ne peut
se faire que par l'intermdiaire d'un ressaut qui fera passer y , ~ ,torrentiel y'nl fluvial. On aura
donc un ressaut suivi d'une courbe de remous T1 dont le contrle y2 est tel que Hsi(y2)= Hsc2
dans le bief 1.
Une courbe de remous T2 de contrle yc2 diminuera la profondeur d'eau de yc2 y a dans le
bief 2.
C'est le contrle de la courbe T1 dans le bief 1 et la courbe T2 la place de l'coulement
normal qui distingue le qemecas du le 3mecas.

Cas :L'nergie spcrque de l'coulement normal amont est suJjfibante pour franchir
la singularit :HsICynl) > Hsd

.Le rgime torrentiel amont s'tablit jusqu' la singularit (puisqu'il n'y a de courbe de remous
permettant de se raccorder un coulement normal torrentiel amont et il n'y a pas de ressaut).

186

Interprtation sur les courbes Hs(y)

Lignes d'eau

Hs2

Hs

'A
-- Y,

h o 2

Hsol

Iiclcl

141c2

Ycl, ycp

y*

1 H62

Hs
k c 2

Hscl

'Y

Y",I

Y"2

y2

Tableau 5-4a : Rtrcissement de la section. La courbe de la section aval (2) est au-dessus de
celle de la section amont (1).
d

187

Cas

Lignes d'eau

Interprtation sur les courbes Hs(y)

h

Hs'

//

Hs
h C 2

)btl

..
y"l

rcp

Yc2

Y"2
k

Y2

1 HB2

kb
WC2
H8cl

,,
'

Y",

,
1 Ycl Yc2

.Y

Y",

Tableau 5-4a (suite) : Rtrcissement de la section. La courbe de la section aval (2) est audessus de celle de la section amont (1).

188

4.2 Elargissement de la section

O

'

La courbe &2 est en dessous de H,s~

et elle est dcale en y vers la gauche (figure 5-4).
l

At

Figure 5-4 : Elargissement de la section.

On peut distinguer les cas typiques suivants (tableau 5-4b) :
1" cas et 2"e cas : Le rgime normal estfluvial avant et aprs I'iargissement
1" Cas iL 'nergie spcifique de 1 'couleqenf normalfluvial aval dpasse 1 'nergie
spcifique critique de l'amont: Hs&,S>HsCr

Une courbe de remous F2 dont le contrle y2 est tel que Hsi(y2) = &2(yn~) diminue la
profondeur d'eau et l'nergie spcifique dans le bief 1 pour franchir la singularit.
Le niveau normal fluvial aval se conserve jusqu' la singularit dans le bief 2 (puisqu'il n'y a
pas de courbe de remous permettant de se raccorder un coulement normal fluvial aval et il
n'y aura pas de ressaut).
cas : L'nergie spcflque critique de l'amont dpasse l'nergie spcflque de
1'coulement normalfluvial aval :Hsz(ynz)<HsC~
Il faut une courbe de remous et ventuellement un ressaut pour diminuer l'nergie spcifique
critique atteinte au niveau de la singularit dans le bief 1.

Une coulabe de remous de type F2 de contrle!yclz diminue la profondeur d'eau et l'nergie

spcifique dans le bief 1 comme dans le l mcas.

189

-a1 et cc3 sont les coefficients dnergie cintique aux sections

respectives 1 et 3.

- a i et a ; sont les coefficients pour exprimer que la rpartition

des pressions nest pas hydrostatique dans les sections 1 et 3
T72
avec lhypothse que a +i et cz ;=l+k,- Y 3
%Y3

On peut dduire de cette expression la vitesse V3 :

Le dbit sera alors donn par Q = C, S3 V3 o Cc est le coefficient de contraction de la

veine. Les diffrents coefficients sont regroups pour donner lexpression suivante du
dbit o le coefficients C doit tre dtermin empiriquement :

Ah=y1-y3est la chute de la surface libre entre les sections

1 et 3.

LUSGS a donn les diffrents facteurs qui influencent le coefficient C sous la forme
suivante (figures 5-7) :

C est une valeur de base correspondant des conditions bien

particulires chaque type de rtrcissement.

Les coefficients K sont des facteurs correctifs pour des

conditions scartant de celles de C :
0

KFtient compte de leffet du nombre de Froude aval

la section 3. Ce nombre de Froude doit tre Srieur
0.8. A priori il nest pas connue mais on procde par
approximations successives en prenant une valeur

197

initiale du dbit. On peut admettre en premire

7

approximation

F,=F,,=V,/Jgy,

Kr tient compte de leffet darrondi des cules

verticales.
KW tient compte du chanfrein W (longueur du mur en
aile dans une direction normale la section 3) et de
langle douverture 9 (angle entre le mur en aile et le
plan du rtrcissement ) dans le cas des cules
verticales.

$ tient compte de leffet du biais du pont. (4 est

langle entre laxe du pont et la section droite de
1coulement).

Ky Tient compte de leffet de la profondeur deau

moyenne rapporte la largeur au droit du pont:
(yb+ya)/2b avec ya et Yb les profondeurs deau aux
pieds des cule, calcules sans le remous.

& tient compte de lexcentricit du rtrcissement

par rapport au champ dcoulement : e = Kfib o Ka
et Kb sont les dbitances (ou coefficients de transfert)
des sections latrales barres par le pont.

K, tient compte du h i t E (horizontavvertical) des

remblais et de la pente des cules reprsente par le
rapport x/b o x est la distance horizontale entre
lintersection du remblai et de la cule, et un point du
remblai amont ayant la mme cote que le plan deau
amont (section 1).
Kt tient compte de leffet de la submersion ventuelle
du tablier
Kj Tient compte de leffet des piles de pont ; fonction
de leurs formes, du rapport j = Aj/A3 o Aj est la
section transversale mouille des piles et A3 la section
mouille aval.

Dans ces calculs :

L est la largeur moyenne de la cule dans la direction de lcoulement excluant la
partie divergente.

198

m = 1 - G est le coefficient douverture o CT = KJKB , avec Kc et KB les

dbitances (ou coefficients de transfert) respectives de la section rtrcie et de la
section totale pour une profondeuq d?eau sans le pont, est le coefficient de
b
rduction. (CT
= - dans le cas dun canal rectangulaire).
B
Dans le cas o les diffrents facteurs correctifs conduisent une valeur de C suprieure a
1, les auteurs de cette mthode prconisent de prendre C = 1.
Lorsque seules les piles de pont constituent lobstacle lcoulement, on peut prendre
tous les facteurs gaux 1 sauf Kj quon prendra sur les courbes (figure 5-7F). Les flches
en traits discontinus indiquent la mthode de dtermination de Kj. Si j est plus grand que
O, 1O, on prendra Kj = 1.
La perte de charge par frottement hf est donne par la formule suivante :
U

- K1 et K3 sont les dbitances (ou coefficients de Transfert

respectives des sections 1 et 3.

- La

est la longueur sparant la section 1 de lentre du

rtrcissement. Sa valeur est gale b la largeur du
rtrcissement.

Si la section 1 est compose de sections S l i ayant des dbitances Kli, i variant de 1 N, le

coefficient a1 est donn par :
B

i=l

1=1

Le remous cr par le pont h; est rapport la valeur Ah qui est la diffrence des cotes
des plans deau aux sections 1 et 3. La figure 5-8a donne la valeur de en hf/Ah en
fonction du coefficient de Manning n et du coefficient douverture m pour le type 1 des
figures 5-7. La rugosit a une influence mineure sur h;/Ah sur cette courbe.
Pour les autres types, on apporte un facteur correctif Ka fonction du coeGcient
douverture m et du rapport C/C type 1 o C et C fype 1 sont respectivement les coefficients
ci-deqsus de dbit du cas tudi et du cas type 1,(figure 5-8b).

199

Figure 5-7A : Coefficients de dbit du type 1 : cdes verticales sans murs en aile;
remblais vertical (E=O) (daprs Chow).

200

Fimire 5-7B : Coefficients de dbit du type II : cuies verticales sans murs en aile ;
remblais talut 111 (E=l) (daprs Chow).

202

02

Fipure 5-7C : Coefficients de dbit du type III : cules et remblais taluts 111 (E=l) avec
raccordement en quart de cne (daprs Chow).

204

1.10

ib

-Gel 1

-}

..,+___t__.

F i m e 5-7C (suite) : Coeflcients de dbit du type III : cdes et remblais taluts 1/2
(E=2) avec raccordement en quart de cne (d'aprs Chow).

205

B.jLnaisdanowu

(4
Figure 5-7D : Coefficients de dbit du type IV : cules verticales avec murs en aile,
remblais taluts 1/1 (E=l) avec raccordement en quart de cne (d'aprs
Chow).

206

Fipure 5-7D (suite) : Coefficients de dbit du type IV : cules verticales avec murs en aile,
remblais taluts 1/2 (E=2) avec raccordement en quart de cne
(d'aprs Chow).

207

. .

Figure 5-7F : Influences de l'excentricit, de la submersion et des piles de pont pour les
types 1, II, III et IV (d'aprs Chow).

208

F i m e 5-Sa : Calcul du remous h*l partir de la perte de charge Ah pour les ponts type 1
(daprs Chow).

W%aiLc

Figure 5-Sb : Facteur de correction de h*1 pour les autres types (daprs Chow).

209

2O)

Formule approche de Carlier

Une formule approche pour les ponceaux ou petits ouvrages transversaux dune seule
trave et comportant seulement 2 pidroits est donne par Carlier :

o Ah est le remous d au paysage du pont

U, la vitesse moyennesans le pont
Cl, le coefficient de perte de charge lentre

C2, le coefficient de perte de charge sous le pont

C3, le coefficient de perte de charge la sortie

Les valeurs de Cl sont donnes ci-dessous :

0
0

C1=0,5 pour des pidroits rectangulaires ;

C1=0,2 pour des pidroits arrondis.

La valeur de C2 est donne par une fornhe de perte de charge linaire (Chzy ou
Strickler) :

o L est la longueur di pont dans le sens di courant ;

C, le coefficient de CHZY la traverse de louvrage ;
et RH,le rayon hydraulique moyen sous le pont.

Le coefficient C3 est donn par :

o s est la section mouille soug le,pont

et S, celle laval du pont.

Pour les ouvrages ordinaires on peut prendre les valeurs de p =1/,/(li+(12+(-:3

suivantes:

p = 0,8 pour de trs petits ponceaux.

p = 0,7 pour des ponceaux dont la longueur mesure dans la direction

du courant excde 20 30 m.

210

Lorsque le frottement (C2) est ngligeable (cas des petits ouvrages), on peut prendre les
volumes de p suivantes :

- Si le radier de louvrage est la mme cote que le fond du canal

pidroits rectangulaires
pidroit arrondis

p = 0,80
p = 0,90

Si le radier de louvrage est plus lev que le fond dun canal :

JJ= 0,72
seuil et pidroit rectangulaires
0
seuil arrondi et pidroit rectangulaires
p = 0.76
p = 0,85
0
seuil et pidroits arrondis

5 PHENOMENES LOCAUX
5.1 Piles de pont

Les piles de pont ressemblent lanalyse effectue pour les rtrcissements localises si lon
considre une tranche longitudinale du courant limite par des plans verticaux de symtrie de
2 piles conscutives (figure 5-9)

Figures 5-9 : Schmatisation de ltude dune range de piles de pont.

Outre uqe diminution de la largeur de Icgulement, la tte de la pile provoque une
contraction de la veine liquide dautant plus marque que sa forme scarte dun profil
arodynamique.
On dtermine la surlvation provoque par les piles (remous) en rgime noy par la formule
de Rehbock (quation 5-4). Cette surlvation est du mme ordre de grandeur que la perte de
charge la borda dans le cas du rgime fluvial. En effet cette perte de charge est donne par la

211

formule donne par Escande :

mme ordre de grandeur.

AH=K(Ui-u3)*/2g+Ah,
avec UI et U3 trs faibles et du

Ah=[6 -0(6-1)](0,4+0~+0~)(1+F,&
2 u2
2g

(5-4)

O=-=-

B-b est le coeffrcient de rduction de la section

B

(rapport de la section rduite par les piles la section totale) ;

F R 3 = U 3 / d Gle, nombre de Froude laval (coulement

non perturb) ;
C,=b/b, le coefficient de contraction (rapport de la section
contracte la section entre les piles) ;
et E = L/e, lallongement de la pile (rapport de sa longueur son
paisseur).
La valeur de 6 es donne en fonction de B et de la forme de la pile par labaque de la figure 51O. Cet abaque est valable dans les conditions suivantes :
Allongement E peu diffrent de 4 ;
Des piles disposes paralllement au courant ;
Ecoulement dans la classe (1) de labaque.

O, 30
tf

Q.%b
1,

a) Zone de validit

b) Valeur de 6 :

- - -Tte rectangulaire ;

-Tte circulaire

Figure 5-10 : Calcui du remous provoqu par des piles par la formule de Rehbock , zone de
validit et valeur de 6 en fonction de o. Valeurs de E =: 4 et 2~ =FZ3 (daprs Lencasre).

212

5.2 Grilles

Une grille constitue dune range de barreaux parallles provoque un remous Ah donn par
la formule suivante :

o U est la vitesse moyenne du courant ;

s, lpaisseur des barreaux dans le sens transversal ;

b, la largeur libre entre les barreaux conscutifs ;

0, langle du plan de la grille avec lhorizontal (sa valeur est
proche de 80 pour faciliter le nettoyage manuelle de la grille);
et fi,un coefficient qui est fonction de la forme de la section
transversale des barreaux (figure 5-1 1).

Figure 5-1 1 : Grilles constitues dune range de barreaux parallles et valeurs de p.

a

5.3 Coudes
Les coudes introduisent diverses perturbations du rgime normal dans un canal. Ces
perturbations se traduisent par la formation de courants secondaires impulsant un mouvement
hlicodal aux particules fluides, des pertes de charge singulires et la modification de la
surface libre avec des ondes transversales dans le cas dun rgime torrentiel.

213
P

On trouve dans la littrature ( Lencastre ou Chow) des mthodes destimation des pertes de
charge et des ondes transversales introduites par un coude. Ces pertes de charges sont
toutefois ngligeables pour des angles 0 infrieurs 45 et pour des rayons de courbures R,
suprieurs 2 fois la largeur au plafond Z (figure 5-12).
On observe galement un dvers Ay d aux effets centrifuges avec une surlvation de la
surface libre du cot extrieur du coude et une diminution du cot intrieur. Un valeur
approche de ce dvers est donn par la formule (5-6) en supposant que les filets liquides ont
un mme rayon de courbure & et que la vitesse est uniforme dans la section et gale la
vitesse moyenne U.

Dans le cas des fonds mobiles, les courants secondaires entranent une excavation du matriau
prs des talus extrieurs des coudes et un dpt prs des talus intrieurs.

Coupe A-A

Vue en plan

Figure 5-12 : Dvers dans les coudes.

8

6 EXTREMITES
DES CANAUX
6. i Passage dun rservoir un canal

Le dbit et la ligne deau ne sont pas connus priori dans le canal. On suppose que la surface
libre dans le rservoir a un niveau constant et que le canal est suffisamment long pour que le
rgime uniforme stablisse laval.
On utilise les courbes dbits profondeur pour une charge spcifique donne qui est gale a la
diffrence de cotes h entre le plan deau dans le rservoir et la crte de la jonction canalrservoir.
On peut distinguer les 2 cas suivants (Tableau 5-6a)
8

214

1 Cas :le rgime normal dans le canal estfluvial

Le rgime normal fluvial aval stablit jusqu la singularit est on observera un abaissement
de la ligne deau dans le rservoir lapproche du canal.
Le dbit du canal et la profondeur normale y smt donns par lintersection (figure 5-13a) de
la courbe Q(y) pour Ii, = h (quation 5-7) et de lquation (5-8) de Manning Strkkler .
(5-7)

La rsolution pour un canal de forme trapzodale donne :

0

une formule implicite en y quon calcule par approximations successives :

un dbit Q qui peut tre calcul par lune des quations, connaissant y par exemple

A,+-

Manning-trickler

Q(y) pour Hr fixlo ii h

Q Q

Q max

Q(Y)

Figure 5-13a : Passage dun rservoir un canal. Mthode de rsolution graphique dans le cas
dun rgime uniforme fluvial

215

2""" Cas :Le rgime normal dans le canal est torrentiel.

Le rgime normal torrentiel sera atteint avec une courbe de remous T2 de contrle yc,
profondeur critique du canal. L'coulement est similaire l'coulement par un dversoir et le
dbit correspondant au dbit maximum est donn par :
(5-1 1)

o yc est la profondeur critique

et Sc = S(yc ) est la section critique.
d

La profondeur normale torrentielle est donne par l'ordonne y correspondant Q = Qmaxsur

la courbe reprsentant l'quation de Manning-Strickler (figure 5-13b) ou par la rsolution de
1' quation (5-12).

Q ~ , = K ~ X Y ):(R
.y@
~

(5-12)

Pour un canal trapzodal, on obtient les formules analytiques suivantes suivantes :

Q,=y,(b+myJ,/2g(h-yc)

pour le calcul du dbit

pour le calcul de y,, par itrations

'9

(5-13)

(5-14)

t'
Q(y) pour Hs f i x b ih

Y'
I
1
l

Figure 5-13b : Passage d'un rservoir un canal. Mthode de rsolution graphique dans le cas
d'un rgime uniforme torrentiel.

216

Interprtation sur les corbes

&(Y)

Lignes d'eau

Hs

dans le rservoir ,,/

,,'

,,,

H&

Hs

Hsc

Tableau 3-6a : Passage d'un rservoir un canal:"

'

217

On peut observer que la diffrence entre les 2 cas est la situation de l'abscisse Q'
correspondant la profondeur critique y, sur la courbe correspondant l'quation de
Manning-Strickler.
Dans le premier cas, Q' Qma et dans le second Q' > Qmm On peut donc tablir sur cette
base l'algorithme de calcul suivant pour un canal de forme trapzodale :
1"). Calcul de y,

b+2myc
3b+5my ,

par itrations si m#O

ou

b-4hm)'+4Ohb
Y,= 4hm-3b+,/(31 Om

par

rsolution

de

l'quation du second
degr si m#O

2"). Calcul de Qmaxpar l'quation (5-13)

3"). Calcul de Q'

4"). Si Q' < Qma, l'coulement normal est fluvial. La profondeur normale et le
dbit Q sont donns respectivement par les quations (5-9) et (5-10).
5"). Si Q' > Qma, le dbit Q est gal au dbit maximal donn par l'quation (5-13)
et la profondeur normal par l'quation (5-14).

6") Pour tenir de la perte de charge l'entre, on multiplie le dbit ainsi trouv
par un coefficient K variant de 0,8 0.9 suivant que l'entre est bords
fiancs ou profils.
6.2 Passage d'un canal B un r6servoir

Le dbit dans le canal est connu ou peut tre dduit de la mthode vue au paragraphe 6.1.
ler,2"" et Pm cas :Le rgime normal dans le canal estfluvial.

1" cas :La cote du plan d'eau est suprieure au niveau normal la sortie du canal
CVI'YnS
Il y aura dans le canal une courbe de remous FI de contrle y. tel que Hs(y,) - AH = y1.

218

2e

cas :Le niveau du plan deau est compris entre le niveau critique et le niveau normal

la sortie du canal (y.<y~<yJ.

On observera dans le canal une courbe de remous F 2 de contrle y. tel que Hs(yo) - AH
3me

= y1 .

cas :Le niveau du plan deau est infrieur au niveau critique la sortie du canal

OI<Y~
<yn>
Il y aura galement une courbe de remous F2 mais la section la sortie est critique. Le
raccordement au rservoir se fera par une chute. Le contrle de la courbe F2 est donc yc.

p.s et dme cas :Le rgime normal dans le canal est torrentiel
d-

cas et

#me

cas :La

au niveau critique la sortie du

cas :La cote du plan deau estI suprieure
,
canal (y1 >yJ.

cote du plan deau est suprieure au niveau conjugu du niveau normal la

sortie du canal (ypy,J.

On observera un ressaut se raccordant lcoulement normal torrentiel amont, suivi dune

courbe de remous Ti de contrle y. tel que Hs(yo) - AH = y1 dans le canal.

cas :La cote du plan deau est infrieure au niveau conjugu du niveau normal la
sortie du canal (yl<y,J.

Lcoulement normal torrentiel viendra jusquau rservoir o lon observera un faux ressaut.
On observe le mme phnomne lorsque la cote du plan deau est comprise entre les niveaux
normal et torrentiel a la sortie du canal.
B

cas :La cote du plan deau est infrieure au niveau normal la sortie du canal
O.(l<Y,J.
On observera une chute brusque la sortie du rservoir.

219

Lignes deau

Interprtation sur les courbes

HSb)

Hs

Hsc

!
!

I
I
I

I
I

Ha
3

Hsc

Hs

H&
l
I

YI Yc

*Y

Tableau 5-6b : Passage dun canal un rservoir.

220

Cas

Lignes deau

Interprtation sur les courbes

Hdv\
Hs

a
h

..h

Hsc

&.a

Faux ressaut
dans le rservoir

Hs

Hsc

Hs

Hsc

Tableau 5-6b(suite) : Passage dun canal un rservoir

221

Chapitre 6
Notions sur les coulements non-permanents

CHAPITRE 6

NOTIONS SUR LES ECOULEMENTS NON-PERMANENTS

1 GENERALITES

2 ECOULEMENTSNON-PERMANENTS GRADUELLEMENT VARIES

2.1 EQUATIONS DE SAINT-VENANT
2.2 CELERITE DE L'ONDE

2.3 MET~ODES
DE RESOLUTION

'

2.4 APPLICATIONS EN HYDROLOOIE

3 ECOULEMENTS
NON-PERMANENTS BRUSQUEMENT VARIES (ONDES DE TRANSLATIONS)
3.1 NOTIONS ET EXEMPLES
3.2 MISE EN EQUATION ET CELERITE DES ONDES
3.3 CONSIDERATIONS SUR LES ONDES NEGATIVES

230

CHAPITRE 6
NOTIONS SUR LES ECOULEMENTS NON-PERMANENTS

1 GENERALITES
Bien que les ouvrages hydrauliques soient dimensionns par les formules des coulements
permanents, lingnieur peut se trouver confront des problmes de grandes importances
ncessitant des notions sur les coulements non-permanents.
De tels coulements se rencontrent dans les cas suivants :
Propagations des mares dans les fleuves et estuaires ;
Ecoulements transitoires louverture ou la fermenture dune vanne place dans un
canal ;
Onde de crue conscutive la rupture dun barrage ( pas souhaitable mais on ltudie
dans le cadre des tudes dimpacts parfois) ;
Propagation des crues dans les cours deau.
Les quations du mouvement ne peuvent tre rsolues que numriquement dans le cas gnral
mais avec certaines simplifications assez proches de la ralit, on peut trouver des solutions
explicites pouvant dgrossir le problme.

2 ECOULEMENTS
NON-PERMANENTS GRADUELLEMENT VARIES

On considre les coulements o le dbit Q et la profondeur deau y varient au cours du temps

et le long de laxe dcoulement de faon progressive. Seuls les coulements filaires, ayant
une direction privilgie (axe du chenal x) et qui sont des coulements unidirectionnels, seont
traits dans le cours.

2.1 Equations de Saint-Venant

2.1.1 Equations compltes de Saint-Venant
Les hypothses sont les mmes que celles faites lors de ltude des coulements permanents
graduellement varies du chapitre 3 :
Variations lentes et progressives des variables de lcoulement (Q et y sont des
fonctions continues et drivables) ;
Filets liquides presque parallles (rpartitions des pressions hydrostatiques dans une
section droite) ;
Vitesse uniforme dans la section (coefficient dnergie cintique CI N 1) ;

231

Vitesse uniforme dans la section (coefficient dnergie cintique a N 1) ;

Pas dapports latraux ;
Le fiottement sur les berges est donn par une formule du type Manning-Strickler.
Lquation de continuit (volume entrant moins volume sortant gale la variation de stock
pendant un intervalle de temps dt) crite pour une tranche de chenal dx donne (stock=Sdx
avec dS / y = 1,largeur au plafond) :

Les termes de lquation dynamique pour le mme volume projets dans la direction x de
lcoulement se dcomposent ainsi :

(ar 3;

Variation de quantit de mouvement : pSdx -+ U-

as,
Forces de pression : - ~ S & ;
dx
4
0

Frottement laparoi : -pgsdxJ avec J = U2/(Ks

2~

-~

;3

Poids de leau pour une pente du chenal, 1, faible : a S d x I .

Le thorme des quantits de mouvement donne alors lquation dynamique qui peut se
mettre sous lune des formes suivantes en tenant compte de lquation de continuit :

au +u-au = - g ay
at

dx

dx

+g ( I - J )

(6-2a)

(6-2b)

aQ + gS(1- F~2) %
QaQ
+ 2 -at

dx

sdx

= gS(I - J)+

Pochat donne une formulation plus complte en partant des quations de Navier-Stokes, en
tenant compte du fiottement de lair qui a une vitesse w par rapport la surface libre et des
apports latraux possibles q [L2T-]. 11 a tablit les 2 quations suivantes pour un coulement
filaire :
1is/ +aQ
- =q
at

dx

232

Pu
aQ + a(Q21s)+ gS(-9)- 1) = 2Cd---W2

at

dX

ax

- gsJ +

Q
x

+ kq-Q
S

o :

Cd est un coefficient de fiottement air-eau et pa la masse

volumique de lair ;
est un cofficient de dispersion tenant compte de la
rpartition relle des vitesses dans la section. Ce
coefficient masque la diffusion visqueuse v, la diffusion
turbulente Kt et la dispersion verticale & ;
&S

k reprsente le rapport de la vitesse dapport et de la

vitesse moyenne dans la section dune part, et langle
entre la direction de lapport et la direction x de
lcoulement avec :
- k=O si dbit latral est perpendiculaire laxe x,
- k=O si dbit latral vitesse trs faible par rapport la
vitesse moyenne axiale U,
- k=l si dbit latral ngatif (perte).
Dans la suite du cours, nous ne retiendrons ni les pertes latrales q, ni la vitesse du vent w ni
la dispersion &s.
2.1.2 Equations simplifies
Il arrive quun ou plusieurs termes de lquation dynamique de Saint-Venant (6-2) puissent
tre ngliges devant les autres :
1 Onde dynamique quasi-permanent

Le terme dU / t de lquation (6-2a) est ngligeable devant les autres.

20)Onde diffusive
Le terme dinertie dU / t + UdU 1dx est ngligeable dans lquation dynamique. Cest le cas
des coulement faible nombre de Froude (Fr<0,3) correspondant aus rivires en rgime
fluvial.
Pour traiter la propagation du dbit, on drive lquation de continuit par rapport x et
lquation dynamique simplifie par rapport t pour liminer le terme a2y/(dxt)et lon
obtient :

233

avec a=-

1
qui est un terme dattnuation et
dJ

e-aQ

a + (1- J ) at

1ay
c = --P aJ

+-021 a x

dx qui est terme de clrit.

dJ

30) Onde cinmatique

On ne retient que les 2 termes 1 et J dans lquation dynamique. Pour un observatuer fixe, la
pente de ligne dnergie J est gale la pente du canal.
2

Lquation

de

Q = KsS ~~3

Manning-Strickler
-2

(@/dx)+d(s~~3)/&
Ks

fi

permet

de

calculer

en fonction de la gomtrie de la

section , de la rugosit et de la pente J=I quon substitue dans lquation de continuit pour
avoir lquation de londe.
Londe cinmatique est une onde de translation qui ne diffuse pas, do il ny a pas
dattnuation de londe. On fait cette hypothse dans le cas des canaux dassainissement :
relative courte longueur et section bien dfinie (circulaire, parois en bton, etc.)
2.2 Clrft de londe
Pour mettre en vidence la nature des quations de Saint-Venant, en peut driver les formules
(6-1) et (6-2c) par rapport x et t et les combiner linairement pour obtenir les quations (6-3)
et (6-4) qui sont les quations classiques des ondes ; en faisant les hypothses skiplifkatices
suivantes :

canal rectangulaire ! . =constante ;

pente du cand 1= pente de la ligne dnergie J ;
nombre de Froude faible devant lunit ;
petites perturbations avec les termes de second ordre qui sont ngligeables.

La clrit de ces petites perturbations c est donc canne par les quations (6-3) et 6-4 et les
hypothses faites ne changent pas sa valeur comme on le verra dans la mthode des
caractristiques du paragraphe 3.3.1 :

234

On constate que la clrit des ondes augmente avec la profondeur deau, do une onde
positive aura tendance dferler avec les hypothses faites. Ce pendant dans certaines
conditions (nombre de Froude initial infrieur 2, fioiement la paroi, etc.) ce dferlement
peut tre retard.

2.3 Mthodes de rsolution

2.3.1 Mthode des caractristiques

Elle est un peu plus complexe que dans le cas des coulements en charge ( paragraphe 3-2 du
chapitre 6 du tome 1) parce que la clrit de londe c est du mme ordre de grandeur que la
vitesse moyenne de leau U et que les caractristiques ne sont plus des droites.
En ajoutant lquation dynamique (6-2c) lquation de continuit (6-1) multiplie par IL, puis
en posant dc / dt = il + 2Q / S = gS(1- F;) /(Al), on obtient lquation diffrentielle suivante

o
.!A2

R = - U * , / g m est lune des solutions de lquation du second degr

+ (2Ql / S ) A - gS(1- ~ f =)O
-Al*-g+gS(I-J)+dt dt

Cette forme condense est quivalente aux 2 systmes dquations diffrentielles (6-5) et (6-6)
appeles quations de compatibilit :

(U-,/m)tdy-dQ+

(6-5)

(6)

&

-= -U - JgsTe
dt

(4
(6-6)

(U+,/~)edy-dQ+

(6)

235

Les (6-5a) et (6-6a) sont respectivement les caractristiques positive C' et ngative C- qui ne
sont plus ncessairement des droites dans le plan (x,t).
Pour un coulement fluvial, la pente de la caractristique positive C' est positive et celle de
C-est ngative, d'o il faut une condition la limite amont de la forme adQ+pdyydt et une
condition la limite aval du mme type.
Pour un scoulement torrentiel, les pentes de 0
' et C- sont toutes les deux ngatives et il faut
2 conditions la limite amont: Q=Q(t) d'une part, et y=y(t) ou Q=Q(y) d'autre part.
Aucune condition n'est ncessaire la limite aval.
Le rgime critique correspond une rsonnance'entre la vitesse de l'coulement et la clrit
des ondes, il est instable.

Dans les deux cas, on doit connatre les conditions initiales.

On utilise .la mthode des diffrences finies explicites pour rsoudre le systme des quations
(6-5) et (6-6) qui dtermine a la fois XM, tM, YM,QM, du poht M de rencontre des 2
caractrisque partant des points Mi et Mz o l'on connat ces valeurs (Figure 6-1).
La mthode des caractristiques reste complexe cause des conditions de stabilit et de
convergence qui impliquent la condition de Courant sur le choix de la discrtisation :
J

Fimre 6- 1 : Schma de discrtisation de la mthode des caractristiques.

236

2.3.2 Mthodes des dinerences finies implicites

Elles consistent partir des quations (6-1) et (6-2c) directement et estimer les drives par
des diffrences finies comme cest expos au paragraphe 3.1 du chapitre 6 du tome 1 du cours.
La mthode implicite est la mthode la plus utilise pour s&chir
courant.

de la condition de

2.4 Applications en hydrologie

2.4.1 Courbe de tarage des cours deau

Si lon utile la formule de Manning-Strckler avec lquation dynamique (6-2a), le dbit est
donn par la formule suivante :

Dans le cas dun coulement permanent uniforme ,la relation entre le dbit Q et la profondeur
deau y travers S et RHest univoque.

Dans le cas contraire, cette relation est non univoqu, en boucle dont la lageur indique
UdU
1 au et de pression @ (figure 6-2a)
limportance des termes dinertie -- +-g h gat
ax

La courbe de tarage (figure 6-2b) donne les renseignements suivants :

Dans un coulement non permanent, le dbit a 2 valeurs diffrentes pour la mme
profondeur deau, y, suivant que le niveau deau monte ou descend.
A une section donne fixe, on observe tout dabord le maximum de vitesse moyenne
Uma, puis le maximum de dbit Qm(uL, puis le maximum de hauteur deau y-; tous

prcds de la pente de la surface libre maximale

I
4

F i m 6 - 2 : Reprsentation schmatique de la relation Q=f(y) et y=f(t) pour un coulement

non permanent graduellement vari.

237

2.4.2 Mthodes simples de propagation des crues

La rsolution des quations compltes de Saint-Venant demande un temps et des moyens de

calcul significatifs. Cest pourquoi on lui prfre, pour ltude de la propagation des crues dans
les fleuves, des mthodes hydrologiques beaucoup plus simples telles que la mthode de
Muskingum, base seulement sur lquation de conservation de la masse applique sur un bief:

dv
-=QE-QS

dt

Les hypothses la base de la mthode de Muskigum sont rsumes par la formule suivante :

Y=

K b Q ; + (1 - X ) Q j

o K, x et X sont des constantes caractristiques de la section

mouille lentre et la sortie du bief. On prend gnralement
une valeur de X comprise entre O et 1 pour pomdrer linfluence
des 2 sections. Dans un canal prismatique on prend X=0,5 et
pour le laminage dune retenue deau, on prend X=O car le
volume nest li quau dbit la sortie. Dans la propagation des
crues, on prend gnralement x=l .

La combinaison de ces 2 quations donne lquation (6-7) o le terme de gauche,

correspondant lhydrogramme entrant, est connu.

( 1 - X ) K ~ + Q S = Q , - k X -d Q E
dt
dt

(6-7)

Une mthode explicite des diffrences finies permet de dterminer lhydrogramme sortant du
bief et on calcule ainsi de proche en proche la propagation de la crue.

3 ECOUL~MENTSNON-PERMANENTS BRUSQUE~IENT VARIES (ONDES DE TRANSLATIONS)

3.1 Notions et exemples

Le front donde prsente un changement rapide de la profondeur deau, y, avec une courbure
prononce de la surface libre. Pour un N observateur )) se dplaant avec la clrit nde londe,
1coulement est permanent brusquement vari.
Cet effet peut avoir pour cause une ouverture instantane dune vanne par exemple, donc une
variation rapide du dbit AQ. On peut distinguer 4 types dcoulements de ce type (figure 6-3) :

238

a) Onde positive daval

Ce phnomne peut tre observ lamont dune vanne place dans un canal et qui est ferme
instantanment. Une onde positive (Ay>O), cest--dire une surlevation de la surface libre par
rapport lcoulement initial se propage vers lamont avec une clrit q.

6) Onde ngative damont

On lobserve laval de la vanne du cas prcdent. Une dpression de la surface libre par
rapport lcoulement initial (Ay<O) se propage vers laval avec une clrit -Q oppose au
cas prcdent.
c) Onde positive damont
On lobserve laval dune vanne place dans un canal et qui est ouverte instantanment. Une
onde positive (Ay>O) se propage vers laval avec une clrit q.
d) onde ngative daval
On lobserve lamont de la vanne du cas prcdent. Une onde ngative se propage vers
lamont avec une clerit -q oppose au cas prcdent.

b) Onde ngative damont

intirrrie,rcrnce

Cconlenrctzt
initiai

Q+AQ

AQc O

d ) Onde ngative dma1

c ) Onde positive damont

Ay> O

l
t i r # ,

ititittricormce

Q+AQ

.+,
1

L
r

4cuiclciiirrit
initial

AQ> O

intumcsccncc

cnicletnctrt

itiitiaf

Q+AQ

iQ> O

F i m e 6-3 : Intumescescences dues une variation brusque, AQ, du dbit dun canal.

239

3.2 Mise en quation et clrit des ondes

a

3.2.1 Mise en quation

Considrons par exemple une onde positive damont (figure 6-4). Un << observateur )) ,li a
londe et qui se dplace avec la clrit G, verra un coulement brusquement vari permanent
qui peut tre tudi de la faon classique.
On encadre lcoulement par 2 sections o la rpartition des pressions est hydrostatique et les
vitesses rgulires, et on applique le thorme des quantits de mouvement en ngligeant le
frottement sur les parois et le poids de leau (courte distance).
Notre observateur )) verra respectivement dans les sections 1 et 2 de la figure 6-4 :
0
Des vitesses relatives de [Ji - ct et 1J2 - ct ;
0
Des dvits relatifs Q,, = Si(IIi - ct) et Qr2= Sz(U2- ct) entrant et sortant ;
0

Des forces de pression de

FP 1 = P g3siY , 1 et F PZ = P 8 szYGZ-

Sil applique le thorme des quantits de mouvement projet dans la dierction de

lcoulement, il trouve :

et lquation de continuit scrit :

(6-9a)
do :

(6-9b)

La combinaison des quations (6-8) et (6-9a) donne, aprs simplification par pS2, la relation
(6- 1O) suivante :
(6-10)
De lquation de continuit (6-2a), on peut tirer lexpression U 2= (VISi + ctSz - crSi)/Sz
quon substitue dans la formule (6- 1O) pour tirer la clrit absolue des ondes ct :
(6-1 la)

(6- 11b)

240

Les relations (6-1 lb) et (6-9b) permettent de calculer par approximations successives q et
Ay=y2-y1 lorsque U1, y1 et AQ sont connus.
Par exemple, la coupure instantanne du dbit dans un canal rectangulaire donne une onde

[ ).

positive daval de hauteur Ay = c

&!

2 ~ 1

Yi+Y2

oii c = i/p-(yl
2YI

+ y2) ; comparer avec la

CU

formule de Michaud vue en ccoulement en charge Ah = - dans le cas yi = y 2 .

g

3.2.2 Crrit des ondes

Dans lquation (6-1 lb) qui exprime la clrit relative de londe de translation, le signe +
correspond aux ondes damont et le signe - aux ondes daval.
Dans le cas dune section rectangulaire (S=by et yG=y/2), si la hauteur de londe A n 2 - y l est
ngligeable devant y1 et y2, la clrit peut sexprimer ainsi :

La clrit relative devient donc ct = uif

fidans le cas des petites ondes

(n

<< 1 .

La clrit absolue de londe de translation est donne par le terme de la racine carre de
lquation (6-1 lb) en considrant une vitesse initiale nulle (U1=0)qui correspond leau au
repos.
On a admis dans la mise en quation que le fiont donde reprsente une discontinuit de la
surface libre qui est stable. En ralit, le constat est le suivant :
Si londe est positive (ypyl), le front est raide et reste stable.
0

Si londe est ngative (y2<y1), londe devient instable et elle tendra vers une courbe
continue.

En effet, londe de translation peut tre considre comme tant la superposition dondelettes.
Chaque ondelette se dplaant avec une clrit absolue c = &. Les ondelettes du haut
auront de ce fait une clrit plus importante que celles du bas. On constate par consquent :
Pour londe positive, les ondelettes du haut vont, en dferlant, absorber celles du bas et
il en rsulte un front raide.
0

Pour londe ngative, les ondelettes du haut se dplacent plus vite que celles du bas et
il en rsulte un front de moins en moins raide.

241

Le frottement la paroi aussi joue un rle dans la modification du front de londe (voir
Henderson).

3.3 Considrations sur /esondes ngatives

Londe ngative, raide au dpart, se dforme plus rapidement que londe positive. On nobteint
pas rellement un vritable front raide et la hauteur Ay a mois de sens physique que pour
londe positive.
Lcoulement devient vite un coulement graduellement vari quon peut traiter par les
quations de Saint-Venant du paragraphe 1.
3

Pour les petites ondes et dans le cas dun canal rectangulaire ou infiniment large, Grafmontre
que la vitesse moyenne de leau U et la clrit relative de londe sont donnes par les
quations (6-12) et (6-13) si on nglige la pente 1 et le frottement J.
(6- 12)
(6-13)

Le signe du haut correspond aux ondes lmentaires qui se propagent ves le sens de
lcoulement initial.
3.3.1 Onde ngative damont
d

On peut lobserver la fermeture dune vanne dans un canal (figure 6-4). La position de londe
est donne par x=c;T en prenant lorigine du temps T linstant de la fermeture (cration de
londe).
Llexpression de la courbe de la surface libre de londe x(y) , et la vitesse moyenne de leau
U(y) deviennent alors :

3.3.2 Rupture de barrage (onde ngative daval)

3

Ce problme a t rsolu par Ritter, cit par Goutx, qui a considr un canl fond horizontal
et sans fiottemnt. Leau dans le barrage a une vitesse nulle en amont et laval est sec linstant
initial (figure 6-5).

242

ly
Figure 6-4 : Fermeture dune vanne ;onde ngative damont.

Ct

= Ui + 3&

-2 f i

La vitesse de leau une profondeuy et linstant T est donne par lquation 6-12 avec le
signe du bas:
6

W Y )= - 2 6 + 2 f i
la clrit de londe (vers lamont et laval), selon iquation (6- 13) devient :
Cr

( Y ) = - 3 G+ 2 f i

Cette formule montre que pour le front de londe positive daval (y=O), la clrit est de
cf(0)= 2&
et pour le creusement de londe ngative damont (y-yl), elle est de
Cf(Y,)=

JE.

La courbe de la surface libre un instant T est donne par lquation x(y) = T ct(y).
Au droit du barrage rompu [(x=O) donc ~(y)=0],la profondeur deau Y b est constante en
4
fonction du temps et gale ya= - y i theoriquement, et la vitesse de leau est de
9

&

2
U(y,) = . Donc le dbit stationnaire correspondant mais aussi la hauteur maximale de
3
londe positive daval peuvent tre estims si la section de la rivire varie peu.
En ralit, la hauteur G pivot )) thorique Y b doit tre corrig a cause de la prsence deau a
laval du barrage qui raidit le front donde sur une hauteur verticale de quelques dizaines de
9

pour-cent de la hauteur pivot. La correction de la hauteur (< pivot )) conduit plutt y , = -yi
16
daprs Stocker cit par Goutx.

243

Cette approche thorique nempche pas de simuler les ondes de rupture de barrag par les
quations compltes de Saint-Venant. Mais dans tous les cas, des hypothses sont ncessaires
sur le mode de rupture et la gense de londe de crue.

Barrage

--+u,=o

Figure 6-5 : Rupture instantane dun barrage ; onde ngative daval

244

BIBLIOGRAPHIE SOMMAIRE

Auteurs

Edition

Titre

A. LENCASTRE

Hydraulique Gnrale

Eyrolles - SAFEGE, Paris,

1996

Ven Te CHOW

Open-channel hydraulics

McGraw-Hill, New Yoork,

1959

M. CARLIER

Hydraulique gnrale et applique

HENDERSON

Open Channel flow

Mac Millan, New York

1966

W.H. GRAF

Hydraulique fluviale
(2Tomes)

Presses Polytechniques
Romandes, EPFL,
Lausanne, 1993

IDELCIK

Mmento de pertes de charge

H. VARLET

Usines de drivaion (2 Tomes)

Eyrolles, Paris, 1959

USBR

Design of small dams

USG Pinting Office,

Washington, 1977

NG. VAN TU

Hydraulique routire

SOGREAH
R.V. GILES
Ch. OBLED
M.POCHAT
D. GOUTX

ALSTHOM
FLUIDES
J. DEBOISSEZON

Eyrolles, Park, 1972

Eyrolles, Paris, 1969

BCEOM, Ministre franais

de la Coopration et du
Dveloppement, Paris, 1981
Ministre franais des
Manuel du technicien du gnie rural
Relations Extrieures, Paris,
1985
Srie SCHAUM,
Mcanique des fluides et hydraulique
Ediscience, Paris,1975.
INPG,Grenoble, 1989
Analyse numrique
ENGREF, Paris, 1979
Ecoulements surface libre
C E W F , Minist&e
Hydraulique des cours deau
franais de lEquipement,
des Transports et du
Logement, Paris, 2000
Catalogues Douvrage de gestion et de - ALSTHOM FLUIDES,
rgulation des coulements dans les
Services Techniques et
canaux. Procds NEYRTEC
commerciaux, La
Courneuve.
Ecoulements surface libre
1 EIER,Ougadougou, 1986