Chap 5 Esti Ecoulement Surface Libre COURSfinal
Chap 5 Esti Ecoulement Surface Libre COURSfinal
Chap 5 Esti Ecoulement Surface Libre COURSfinal
– Formule de débit
– Formule de vitesse
4
REYNOLDS Re
VD VD 4Q
Re
D
visc.dynamiqueN.S m
2
m s2
7
• Canal: système de transport dans lequel l’eau
s’écoule et dont la surface libre est à la
pression atmosphérique
Canaux
Naturels Artificiels
Géométrie des canaux
B
B = b + 2mh
h
A = (b + mh)h 1
m
b
P = b + 2h(m + 1)
2 1/ 2
R = A/ P
Conduite circulaire
(diamètre D)
h
θ = 2 cos (1 − 2 )
−1
D
θ p = 2π
A=
D2
(θ − sin θ ) πD 2
Ap =
8 4 θ h
θD Pp = πD
P=
2
A D sin θ Rp =
D
R = = (1 − )
P 4 θ 4
Rayon Hydraulique Rh
Rayon hydraulique Rh
Rh A P
A section mouillée
P Périmètre mouillé
Conduite circulaire pleine
2
A D P D R h D 4
4
12
Rayon Hydraulique
Amouillée
Rh
Pmouillée
D2
Amouillée *1 2
4
Pmouillée D *1 2
Pm b 2h
b * h
h
b Am
13
Manning Strickler
1
V R 2 3
S 1 2
n h
Q V .A
S Hf L
14
Manning Strickler
0 , 397
V D 2 3
S 1 2
n
1 2
0 , 397 hf
V D 2 3
L
n
1
n k strickler
15
Manning Strickler
0 , 312
Q D 8 3
S 1 2
n
1 2
0 , 312 hf
Q D 8 3
L
n
S gradient Hydrauliqu e
pente
S 10,29n 2
Q 2 16 3
D
16
2
Q
S 10 , 29 * n * 16 / 3
2
D
0 , 397
V *D *S2/3 1/ 2
n
0 , 312
Q VA 8/3
*D *S 1/ 2
n
R D / 4 (plein) A D / 4 2
17
Manning Strickler
2
Hf Q
S 10,29 * n * 16 3
2
L D
Hf 10,29 * n * L * Q D
2 2 16 3
RQ n
18
2 2
10.29* L * n *Q
Hf 16/ 3
RQ n
D
2
10.29* L * n
Résistance R 16/ 3
D
19
20
HAZEN - WILLIAMS
V 0,849 * Chw * R 23
h *S 12
23 12
A Hf
V 0,849 * Chw * *
P L
A
R h Rayon Hydrauliqu e
P
21
Écoulement à surface libre
• Optimisation des sections des
canaux
– Sections rectangulaires
– Sections trapézoïdales
• Énergie spécifique Es
• Profondeur critique Yc
– Sections rectangulaires
– Sections quelconques
28
• Formule de débit selon Manning strickler
• Equation de débit
• Ecoulement uniforme
Q CA R h
S
C = Coefficient de CHEZY
A = Section d’écoulement (mouillée)
Rh = Rayon hydraulique = A/Pm Pm=
Périmètre mouillé
S = Pente du canal
Pm = Périmètre Mouillé
C prendra une expression selon l’auteur
29
• Pour la formule de Manning
• C= coefficient de Chézy prendra la forme
suivante
1 16 Cm
C Rh Q 23 12
ARh S
n n
Cm = 1 dans le système SI
A= section d’écoulement
S = pente
Rh = rayon hydraulique 30
Q CA S
h
1 16
C Rh Manning
n
87
C Bazin
1 / Rh
31
• optimisation
2 5
P= périmètre mouillé
C= constante
Selon cette formule A minimum
entraîne P minimum
32
• Cas d’une section rectangulaire
•
• P = b+2y et b = P-2y A = by
A P 2 y y cP 25
dP 2 3 5 dP
2y P 2y cP
dy 5 dy
33
2
34
• On trouve pour l’extremum b 2y
1
y
m
b
P b 2y 1 m 2 A by my cP
2 25
35
• Procédure :
• Dériver par rapport à y en supposant m la pente
constante puis par rapport à la pente m.
P 4y 1 m 2
2 my
2m 3
1 m
1 m 2 3
Pente m = H/V 36
avec la valeur m (pente) constante nous avons
P 2
P 4y 1 m 2 my 0
y m cte
P 4y 1 m 2 2 my
avec la valeur y constante nous avons
P 2m
0 .0
m y cte 2
1 m
3
m
3 37
• Solutions pour le canal
3
P 2 3y b2 y
3
A Y
A 3y 2 R
P 2
3
m
3 38
• Solutions pour le canal (pour m donné)
P 4 y 1 m² 2my
2
A by my
A Y
R
P 2 39
•Remarque :la section
semi-circulaire
(demicercle) est la
section optimale de
toutes les sections de
canaux
40
• Exemple: Déterminer les dimensions optimales d’un
canal trapézoïdal de débit Q = 200 m3/s et de pente
S= 0.0004 et n= 0.016
1
Q ARh S
23 12 23
1.00 2 Y
n 200 3y 0.0004
0.016 2
y 83
146.64 y 6.492m
3 3
b 2 y m
3 3
41
23 12
Q ARh S
n
23
1.00 2 Y
200 3y 0.0004
0.016 2
42
3
b 2 y
3
3
m A Y
3 R
P 2
43
• Énergie spécifique et Profondeur critique
d’un canal
• C’est l’énergie rapportée à un repère
horizontal (datum = fond du canal)
2 2
Es y
V Q 1
Es y
2g A 2g
Cas d’un canal rectangulaire
Surface A = b*y avec b= 1 q est un débit unitaire
(par unité de largeur) nous aurons
q2
Es y 2
44
2 gy
• On peut s’intéresser à la variation de Es en
fonction de la profondeur y avec q = constant (
exemple q = 5.5 m3/s) courbe Es en abscisse et y
en ordonnées
• Selon cette courbe (une hyperbole) Es prend une
valeur minimale en dessous de laquelle Es
n’existe pas. On l’obtient par l’équation suivante :
q 2 dE s q2
Es y 0 1 3
2 gy 2 dy gy
45
2
q dEs q 2
Es y 2 0 1 3
2gy dy gy
1 3
q
2 3
y c E s min yc
g
2
c c 46
• Ici yc est la profondeur critique et la vitesse
associée Vc ou vitesse critique
13
q2
y c
3
g E s min yc
2
Vc gy c
On réalise que q est maximum pour une profondeur
égale à la profondeur critique Yc
47
• Dans le cas général nous avons
2
2 g E s y
Q
Es y
2gA 2 QA
• et la profondeur critique peut être obtenue par
dérivation comme suit
dE s Q 2 dA
0 1
dy gA 3 dy
10
0
0 1 2 3 4 5 6 7 8
53
12,00
10,00
8,00
6,00
4,00
2,00
0,00
-0,20 0,00 0,20 0,40 0,60 0,80 1,00 1,20 1,40 1,60
-2,00
54
• Exemple : Déterminer la profondeur critique Yc
d’un canal trapézoïdal de largeur b= 3m et de
pente m = 1/2 (1H/2V) Q= 10 m3/s
2
y
A 3y T 3 y
2
10 3 y
2
10.1983 y
f y 1 .0
9.806 3y y / 2
2
3y 0.5y
3 2 3
56
• Récapitulons
• Section optimale formule de base
A cP 2 5
• Section rectangulaire b, y
• P = b+2y et b = P-2y A = by
• On a
• b = 2*y et P= 4*y et A = 2*Y 2
57
• Récapitulons
• Section optimale formule de base
A cP 25
• Section trapézoidale b m et Y
P b 2y 1 m 2
• On a 3
m
3
58
• Solutions pour le canal trapézoidal
3
P 2 3y b2 y
3
A Y
A 3y 2 R
P 2
59
Sections critiques
2
Q
E s y 2
2 gA
Q A 2 g Es y
2
f y
Q T
3
1
gA 60
61
62
Ecoulement partiel en conduite
• Abaque de Manning
• En ordonnées Y/D
• En abscisses paramètres
• Débit
• Rayon hydraulique
• Vitesse
• Section d’écoulement
63
Rayon hydraulique des collecteurs à
section circulaire, non pleins en fonction
du remplissage
Supposons la conduite remplie à r% (rapport de
la hauteur du liquide dans la conduite au
diamètre de celle-ci = Y/D) et soit a le rapport
entre le rayon hydraulique de la conduite et son
diamètre (a = Rh/D).
Les valeurs de a en fonction de r ont données
par le tableau suivant :
64
r% 20 25 30 35 40 45 50 55 60
r% 65 70 75 80 85 90 95 100
a 0,288 0,296 0,302 0,304 0,303 0,298 0,286 0,250
65
Débits et vitesses dans des collecteurs à section circulaire,
non pleins, en fonction du remplissage.
Considérons une conduite à section circulaire de pente
uniforme, transportant un liquide et partiellement
remplie soient :
Q le débit de cette conduite, supposée remplie à r
% (rapport de la hauteur du liquide dans la conduite au
diamètre de celle-ci), et
V la vitesse correspondante à q
Q son débit si elle était entièrement remplie sans
que la veine supérieure du liquide soit soumise à
aucune pression
V sa vitesse correspondante Q, c’est-à-dire à pleine
66
section ;
r% 10 15 20 25 30 35 40 45 50 55
m 0,21 0,049 0,088 0,137 0,196 0,263 0,337 0,416 0,500 0,586
p 0,401 0,517 0,615 0,700 0,776 0,843 0,902 0,954 1,000 1,039
r% 60 65 70 75 80 85 90 95 100
m 0,672 0,756 0,837 0,912 0,977 1,030 1,066 1,075 1,000
p 1,072 1,099 1,120 1,133 1,140 1,137 1,124 1,095 1,000
8 32 sin 1
2
71