Chap 5 Esti Ecoulement Surface Libre COURSfinal

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

RAPPELS HYDRAULIQUE

• Équations de base en hydraulique


– Équation de continuité
– Équation de quantité de mouvement
– Équation du moment de quantité de mouvement (moment
cinétique)
– Équation de Bernoulli
– Équation d’énergie
• Problèmes en hydraulique
– Calcul du débit
– Calcul de pertes de charge
– Calcul de diamètre (dimensionnement)
• Formules en hydraulique
1
• Formules de Darcy weisbach
– Formule de pertes de charge
– Facteur f
• Cas du régime laminaire
• Cas du régime turbulent
– Formule de colebrooke White
– Formules approchées
– Abaque de Moody
– Formule quadratique des pertes de charge
– Formule de débit
– Formule de vitesse
2
• Formules de Manning Strickler
–Formule de pertes de charge
• Facteur de rugosité n
• Rayon hydraulique Rh
–Formule quadratique des pertes de
charge
–Formule de Débit
–Formule de vitesse
3
• Formule Hazen Williams
– Formule de pertes de charge
– Cœfficient de Hazen williams
– Formule quadratique des pertes de
charge

– 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

¾ Les paramètres suivants sont à retenir :


¾ Largeur du canal à sa surface libre B
¾ La section d’écoulement ou section mouillée A
¾ Le périmètre mouillé P
¾ Le rayon hydraulique R=A/P
¾ La profondeur hydraulique D = A/B
¾ Profondeur d’eau h: profondeur maximale (sauf
précision).
Canal trapézoïdal

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)

Section Section pleine


partiellement pleine

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
  2y  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

A bymy  P2y 1m ymy  cP


2 2 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 b2 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 QA
• et la profondeur critique peut être obtenue par
dérivation comme suit
dE s Q 2 dA
 0 1
dy gA 3 dy

Avec T = largeur de la surface libre, nous avons D =


profondeur hydraulique
A = T*Y et D = A/T Ainsi dA = Tdy
48
2
Q Ac
3
Tc  1 Es  yc 
gA c 2Tc

• Cette équation montre que Es minimum est


obtenue avec l’énergie cinétique est moitié de
A/T ou D
• L’obtention de la profondeur critique pour Q
max est donnée par la résolution de l’équation
2
Q T
f y  3
gA
49
2
Q
Es  y  2 Q  A 2g Es  y
2gA
2
dE s Q dA
 0  1  3
dy gA dy
2
Q T
f y   3
gA
50
51
52
12

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.1983  y 
f y     1 .0

9.806 3y  y / 2
2
 3y  0.5y 
3 2 3

•On trouve Yc = 0.984 m 55


q:=10;g:=9.81;T:=3+y;
A:=3*y+1/2*y^2;
q:=q^2*T/(g*A^3)=1;
solve(eq,y);

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

A bymy  P2y 1m ymy  cP


2 2 2 25

• On a 3
m 
3
58
• Solutions pour le canal trapézoidal

3
P  2 3y b2 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

a 0,121 0,147 0,171 0,193 0,214 0,233 0,250 0,265 0,278

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

Remarque : les valeurs de m et p données pour celles


de r qui dépassent 80 % n’ont qu’un intérêt théorique :
en pratique, on ne rencontre pas de conduite dans
lesquelles r présent de façon stable des valeurs voisines
de 100,
Q est donc le débit de la conduite supposée entièrement
pleine pour lequel la perte de charge est égale à sa
67
pente.
q et v s’obtiennent en fonction
de Q et V, par relations
q = m Q et V = p.V
dans lequel m et p sont des
coefficients donnés par le
tableau suivant en fonction r
r= Y/D
68
69
A 
1
  sin  d 0
2 1
8 P  d 0
2
1  sin 
R  1   d
 
0
4
T  2 y (do  y )
70
 
1    sin  
D  do
8  sin 1  
 2 
1 .5
 
2    sin  
Z    do 2 .5

8  32 sin 1  
 2 
71

Vous aimerez peut-être aussi