******

Unité thématique d’enseignement et de recherche (UTER)

Gestion Valorisation de l’Eau et Assainissement (GVEA)
******

COURS D’HYDRAULIQUE A SURFACE

LIBRE

Angelbert Chabi BIAOU

Enseignant-Chercheur, Hydraulique Générale/Hydrogéologie

1
Version Octobre 2010
Sommaire
Avant propos ................................................................................................................................................. 6
Notation ........................................................................................................................................................ 7
1. Introduction ........................................................................................................................................... 8
2. Définitions ............................................................................................................................................. 9
2.1 Paramètres géométriques ............................................................................................................. 9
2.2 Notion de vitesse ......................................................................................................................... 10
2.2.1 Vitesse instantanée ............................................................................................................. 10
2.2.2 Vitesse moyenne ................................................................................................................. 10
2.2.3 Profil de vitesse ................................................................................................................... 10
2.3 Le courant liquide ou veine liquide ............................................................................................. 11
2.3.1 Ligne de courant .................................................................................................................. 11
2.3.2 Surface de courant, tube de courant, filet liquide, veine liquide ........................................ 11
2.4 Notion de charge hydraulique, théorème de Bernoulli .............................................................. 12
2.5 Différents types d’écoulement .................................................................................................... 13
1. Définition ............................................................................................................................................. 14
2. Etude des écoulements uniformes ...................................................................................................... 14
2.1 Principe fondamental .................................................................................................................. 14
2.2 Calcul de la pente d’énergie. ....................................................................................................... 14
2.3 Equations de l’écoulement uniforme .......................................................................................... 14
2.4 Estimation des valeurs de ..................................................................................................... 15
2.5 Estimation des valeurs de dans le cas des des sections composées..................................... 15
3. Types de problème dans un écoulement normale ............................................................................. 17
3.1 Problème de type 1 ..................................................................................................................... 17
3.2 Problème de type 2 calcul de la profondeur normale ................................................................ 17
3.2.1 Notion de profondeur normale ........................................................................................... 17
3.2.2 Méthode des abaques de recherche de la profondeur normale ........................................ 19
3.2.3 Méthodes graphiques de calcul de la profondeur normale ................................................ 20
3.2.4 Méthode numérique de calcul de la profondeur normale. ................................................ 22
3.2.5 Exemple d’application : ....................................................................................................... 22
3.3 Problème de type 3 forme optimale ........................................................................................... 28

2
 3.3.1 Section hydrauliquement favorable de forme circulaire .................................................... 28
3.3.2 Section hydrauliquement favorable de forme trapézoïdale ............................................... 28
3.3.3 Section hydrauliquement favorable de forme trapézoïdale : construction graphique ...... 29
1. Définition ............................................................................................................................................. 31
2. Propriétés des Ecoulement graduellement variés (EGV) .................................................................... 31
3. Charge spécifique ................................................................................................................................ 31
3.1 Définition ..................................................................................................................................... 31
3.2 Etude de la charge spécifique ..................................................................................................... 32
3.2.1 Etude de la variation de la charge spécifique en fonction de la profondeur ...................... 32
3.2.2 Classification des régimes d’écoulement : le nombre de Froude ....................................... 34
3.2.3 Variation du tirant d’eau pour un débit donné, lorsque est fixé : ................................ 35
4. Etude du régime critique ..................................................................................................................... 36
4.1 Méthode des abaques. ................................................................................................................ 37
4.2 Méthode graphique..................................................................................................................... 37
4.3 Variante de la méthode graphique ....................................................... Erreur ! Signet non défini.
4.4 Méthode numérique ................................................................................................................... 39
4.5 Exercice d’application .................................................................................................................. 40
4.5.1 Cas d’une section trapézoïdale pure ................................................................................... 40
4.5.2 Cas d’une section circulaire ................................................................................................. 42
4.5.3 Cas d’une section composée ............................................................................................... 45
4.6 Notion de pente critique ............................................................................................................. 46
5. Les courbes de remous ........................................................................................................................ 46
5.1 Equation général des écoulements graduellement variés .......................................................... 46
5.2 Etude analytique des courbes de remous ................................................................................... 47
5.2.1 Position relative de I et J, par rapport aux positions relative de et de : ..................... 47
5.2.2 Signe de =1 2 3 en fonction des positions relatives de y et de ............... 48
5.2.3 Etude qualitative des courbes de remous ........................................................................... 48
5.2.4 Quelques illustrations.......................................................................................................... 50
5.2.5 Notion de section de contrôle ............................................................................................. 51
5.3 Calcul des lignes d’eau ou encore axe hydraulique .................................................................... 52
5.3.1 Méthode d’intégration graphique ....................................................................................... 52
5.3.2 Méthodes de Bakhmeteff.................................................................................................... 55

3
 5.3.3 Méthode de Bresse ............................................................................................................. 58
1. Généralité ............................................................................................................................................ 62
2. Le ressaut hydraulique ........................................................................................................................ 62
2.1 Définition ..................................................................................................................................... 62
2.2 Classification des ressauts hydrauliques ..................................................................................... 63
2.3 Théorème de quantité de mouvement (canal rectangulaire) ..................................................... 63
2.4 Notion d’impulsion totale ........................................................................................................... 64
2.4.1 Définition ............................................................................................................................. 64
2.4.2 Etude des variations de l’impulsion totale .......................................................................... 64
2.5 Calcul de ressaut dans le cas d’un canal rectangulaire ............................................................... 66
2.5.1 Position du ressaut .............................................................................................................. 66
2.5.2 Estimation graphique de la perte de charge du ressaut ..................................................... 67
2.5.3 Cas d’un canal rectangulaire ............................................................................................... 68
Annexe A : Rappel de calculs géométriques élémentaires ........................................................................ 72
A1‐ Rappels sur les Systèmes de coordonnées ........................................................................................... 72
A1.1 Système de coordonnées cartésiennes .......................................................................................... 72
A1.2 Système de coordonnées polaires .................................................................................................. 72
A1.3 Système de coordonnées cylindriques ........................................................................................... 73
A1.4 Système de coordonnées sphériques ............................................................................................. 74
A2‐ Elément de calcul de périmètre ............................................................................................................ 74
A2.1 Périmètre de cercle ......................................................................................................................... 74
Cercle plein .......................................................................................................................................... 74
Cercle non fermé ................................................................................................................................. 75
A2.2 Périmètre d’une parabole ............................................................................................................... 75
A2.3 Périmètre d’un rectangle arrondi sur les bords de sa base ............................................................ 77
A2.3 Périmètre d’un triangle à fond arrondi ........................................................................................... 78
A3 Elément de calcul de surface ................................................................................................................. 79
A3.1 Calcul de l’aire de rectangle ............................................................................................................ 79
A3.2 Calcul de l’aire d’un triangle connaissant la base et la hauteur. .................................................... 79
A3.3 Calcul de surface d’un disque ......................................................................................................... 80
A3.2 1 Cas d’un disque plein ............................................................................................................... 80
A3.2 2 Cas d’un disque en partie plein ................................................................................................ 80

4
 A3.4 Calcul de surface d’une sphère ....................................................................................................... 81
A3.4 Aire d’une parabole......................................................................................................................... 83
A3.5 Aire d’un rectangle arrondi sur les bords de la base ...................................................................... 83
Annexe B : Table de la fonction de Bakhmeteff ..................................................... 84
Annexe B : Table de la fonction de Bakhmeteff ..................................................... 85
Annexe B : Table de la fonction de Bakhmeteff ..................................................... 86
Annexe B : Table de la fonction de Bakhmeteff ..................................................... 87

5
 Avant propos
Ce cours d’hydraulique est destiné aux étudiants de divers niveaux, régulièrement
inscrits au sein 2iE, tant en formation présentielle qu’en formation à distance. Sans être
exhaustif, ce cours présente les notions essentielles qu’il faut connaitre pour aborder un
problème d’hydraulique à surface libre et pour préparer les étudiants à bien d'autres cours
enseignés, dans le cadre des formations dispensées dans notre institution académique à savoir
le cours d'hydraulique fluviale, le cours d'hydraulique routière, le cours de barrage, le cours
d'assainissement pluvial, le cours d'irrigation gravitaire, le cours d'évacuation des eaux usées, le
cours d’AEP, pour ne citer que ces quelques uns. Son caractère sommaire fait que seules
quelques démonstrations sont présentées (notamment celles qui ne sont pas présentées dans
les ouvrages courants). Ceux qui désirent un approfondissement peuvent se référer aux
ouvrages spécialiser présentés à la fin de ce document, ou s'appuyer également sur le cours de
Lamine Mar, également disponible au 2iE, afin d’étudier certaines notions plus en profondeur.
Par ailleurs, les exercices d'applications et de synthèses présentés dans ce cours lui donnent une
bonne compréhension des notions présentées ce qui devrait faciliter un self‐apprentissage à
ceux qui désirent s'initier par eux mêmes au cours d'hydraulique à Surface Libre. L'une des
choses essentielles qu'il faut savoir pour un cours d'HSL étant les calculs de section mouillée, de
périmètre mouillé, une partie des annexes est consacrée à ces notions surtout sur les sections
définies analytiquement.
Il n’est pas inutile d’attirer l’attention sur le fait que tous les calculs effectués dans cette
discipline nécessitent des notions sommaires de calculs trigonométriques. Les méthodes de
résolution des équations différentielles telles que la méthode des différences finies ou encore
de Runge‐Kutter sont nécessaire. Par ailleurs, la connaissance de quelques méthodes de
résolution des équations implicites telles que la méthode de newton‐Raphson serait d’une
grande utilité

6
 Notation
SYMBOLE UNITE DEFINITION
C m1/2s‐1 coefficient de Chézy
L m largeur au miroir d’un chenal
m longueur d'un seuil déversant
Lr m longueur d'un ressaut
M kg.m/s quantité de mouvement
P m périmètre mouillé
Q m3/s débit liquide
F ‐ nombre de Froude
H m charge moyenne dans une section
HP m charge en un point P
Hs m charge spécifique Hs = y + V2/2g
Hsc m charge spécifique critique
K coefficient de rugosité (ou de Strickler) global du lit mineur
R m Rayon hydraulique R= S / P
S m2 section mouillée
U m/s vitesse moyenne U= Q/S
Vc m/s vitesse critique
cc m/s célérité de propagation d'une onde cinématique
cd m/s célérité de propagation d'une onde diffusive
p Pa m pression hydrostatique en un point hauteur du seuil coté amont (pelle)
g m/s² accélération de la pesanteur g =9,8m/s²
i ‐ pente du fond
j ‐ perte de charge linéaire
l m distance entre 2 sections
t s temps
v m/s vitesse en un point dans une section d'écoulement
x m abscisse d'un point
y m tirant d'eau
yc m tirant d'eau critique

yn m tirant d'eau normal

zf m cote du fond du chenal
zP m cote d'un point P dans une section d'écoulement
rw kg/m3 masse volumique de l'eau ( rw= 1000 kg/m3)
3
gw kN/m poids volumique de l'eau (gw »9,8kN/m3)
s ‐ coefficient d'atténuation d'une onde diffusive
α degrés angle avec l'horizontale du fond du chenal
β ‐ coefficient relatif à la répartition des vitesses dans une section
m kg.m‐2.s‐1 coefficient de débit d'un seuil viscosité dynamique de l'eau m = [Link]

7
 I‐ Généralités
1. Introduction
L’hydraulique à surface libre est la partie de l’hydraulique générale qui étudie les canaux
à ciel ouvert, les cours d’eau naturels aménagés, les canalisations ne fonctionnant pas à pleine
section. L’hydraulique générale utilise les principes fondamentaux de la mécanique comme la
Relation fondamentale de la dynamique, le théorème de l’énergie cinétique, le théorème des
quantités de mouvement, le théorème de la conservation de l’énergie totale…
Les équations générales de l’hydraulique sont établies par la mécanique des fluides qui
est une branche de la mécanique des milieux continus ou déformables qui elle aussi est une
branche de la mécanique au sens large, définie comme « la science qui étudie le mouvement des
corps matériels, de ses causes et de ses effets ». Les simplifications de ces équations obéissent à
des hypothèses qui dépendent du type d’écoulement. Ainsi, pour chacun des trois types
d’écoulements définis ci‐dessus, des simplifications des équations générales permettent d’avoir
des équations qui permettent de les étudier.
Ecoulements à surface libre.
Ce sont les écoulements pour lesquels la surface libre est en contact avec l’atmosphère.
Pour ces écoulements, la surface libre est soumise à la pression atmosphérique. On rencontre
ce type d’écoulement dans les rivières naturelles ou aménagés, dans les canaux d’irrigation, de
drainage, dans les conduites dont le liquide n’occupe pas toute la section, dans les réseaux
d’assainissement urbain (eaux pluviales, eaux usées).
Ecoulements en charge
Pour ces types d’écoulement, le liquide ne possède pas de surface en contact avec
l’atmosphère. Ces types d’écoulements sont généralement rencontrés dans les conduites des
réseaux d’alimentation en eau potable, les conduites d’irrigation, dans les conduites forcées
alimentant des turbines hydrauliques.
Ecoulements souterrain
Ces écoulements ont lieux dans les milieux dits poreux, constitués de matériaux
granulaires. Ces écoulements ont lieux de façon générale dans les sols ou dans certains
ouvrages artificiels comme les barrages, les digues, remblais. Selon les perméabilités des
différentes couches rencontrées, ces écoulements peuvent être à surface libre (comme cela
peut être le cas dans certaines nappes superficielles dans des milieux grossiers) ou en charge.
En fonction de certaines caractéristiques des milieux poreux, une nappe, siège d’un tel
type d’écoulement peut être captives, libre ou confinée.
Pour aborder l’hydraulique à surface libre qui fait l’objet de ce cours, nous admettons les
quelques hypothèses ci‐après :
‐ Les chenaux qui seront étudiés dans ce cours seront supposés suffisamment
rectilignes, de ce fait, les problèmes que nous aborderons, tant en écoulement
permanents qu’en écoulement transitoire, sont des problèmes à une dimension.
‐ L’eau véhiculée dans les différents chenaux (rivières, canaux…) est supposée être une
eau claire
‐ les parois et fond des chenaux sont supposés fixes
8
 2. Définitions
2.1 Paramètres géométriques
La position de toute section du chenal dans le plan perpendiculaire à son axe sera définie par son
abscisse x. Les paramètres essentiels sont (voir figures 1 et 2) :
¾ le tirant d’eau (y) qui est la hauteur (mesurée par rapport à la partie la
plus profonde de la section considérée) d’eau entre la surface libre du
liquide et le fond du chenal ;
¾ la largeur en miroir (L) ou largeur en gueule ou encore largeur du plan
d’eau qui est largeur de la surface mouillée au contact de la surface libre ;
¾ la section mouillée (S) qui est la surface de la section du chenal depuis la
surface libre jusqu’au fond du chenal ;
¾ le périmètre mouillé(P) qui est la longueur de la ligne de contact, pour la
section donnée, de l’eau avec les parois et le fond du chenal ;
¾ le rayon ( ) hydraulique qui est le rapport entre section mouillée et
périmètre mouillé.
¾ le Diamètre hydraulique ( 4

L
y

Figure 1: Paramètres géométrique d'une section de chenal naturel

L
dz y Ligne d’eau
Rev
Z1 fond
Z2
y 1
m dx

b
Figure 2: Paramètres géométrique d'une section de canal artificiel (section trapézoïdale).

Dans le cas particulier d’un canal artificiel de forme trapézoïdale, on peut compléter ces
définitions de paramètres géométriques par les suivantes :
¾ la largeur en radier ou base(b)

9
 ¾ les fruits des berges (m) permettent de définir l’inclinaison du canal par
rapport à la verticale (ou horizontal)

Enfin, un des paramètres géométriques indispensable à l’étude des chenaux, artificiels

comme naturels, siège d’un écoulement à surface libre est la pente du fond du canal,
généralement notée . Elle est définie par , z étant la cote du fond du chenal à la section
d’abscisse x (figure 2).
Il est important d’attirer l’attention sur le fait que l’appellation « géométriques » des
paramètres ci‐dessus définis est liée au fait que ces paramètres caractérisent la géométrie de la
section du canal étudié. Certain de ces paramètres (L, y, S, P, R) dépendent du débit et ne
peuvent pas être considérés comme des constantes géométriques simples. La pente quant à elle
est indépendante du débit, mais peut dépendre de l’abscisse du chenal. Pour bien aborder un
cours d’hydraulique à surface libre, il est important de savoir calculer ces paramètres.
Cependant, Mar, 2004 a proposé les résultats de calculs pour quelques géométries usuelles, et
quelques démonstrations sont proposées en annexe de ce document.
2.2 Notion de vitesse
2.2.1 Vitesse en u point
La vitesse en un point de l’écoulement, généralement notée ( minuscule), est celle de
la particule qui passe en ce point au moment considéré. Elle est fonction donc et de l’espace et
du temps : , . On peut adopter la vitesse moyenne (notée ) en un point comme étant
la moyenne dans le temps de plusieurs valeurs de vitesse en ce point :
1

Puisqu’il est pratiquement impossible de suivre chacune des particules de l’écoulement,

on fait recours, pour un certains nombre de calcul, à ce qu’on appelle vitesse moyenne.

2.2.2 Vitesse moyenne

C’est la moyenne des vitesses ponctuele sur une section élémentaire ds. Elle peut
s’exprimer comme U (U majuscule) dans ce cours. Son expression est

On peut l’exprimer à partir de la connaissance du débit qui passe dans la section par
Profil de vitesse

Lorsqu’on mesure la vitesse de l’écoulement, en quelques points, sur une même

verticale, depuis la surface jusqu’au fond, le profil des vitesses qu’on observe est présenté sur
la fFigure 3

10
 δ

Figure 3: Profile des vitesses sur la verticale d'un point donné

La vitesse croit depuis la surface jusqu’à une vitesse maximale, puis décroit vers le fond
pour s’y annuler. On peut y distinguer deux zones différentiées, l’une ayant un faible gradient
transversal dans laquelle les forces de viscosité sont négligeables par rapport aux forces
d’inertie et de turbulence, et l’autre ayant un fort gradient transversal appelée encore couche
limite, au voisinage de la paroi, dont l’épaisseur est matérialisée par δ. Se référer à Carlier, 1972
pour plus de précision sur la notion de la couche limite
2.3 Le courant liquide ou veine liquide
2.3.1 Ligne de courant
Une ligne de courant est une courbe tangente en chacun de ses points au vecteur vitesse
en ce point. En écoulement non permanent, la vitesse v en tout point évolue dans le temps et
les lignes de courant se déforment avec le temps. En écoulement permanent, les lignes de
courant ne se déforment pas et constituent des trajectoires de particules d’eau. Le profil de la
surface libre est une ligne de courant particulière.

2.3.2 Surface de courant, tube de courant, filet liquide, veine liquide

Lorsqu’une infinité de lignes de courant s’appuient sur une courbe C donnée, on obtient
ce qu’on appelle une surface de courant. Si de plus la courbe C est fermée, la surface de courant
devient un tube de courant. Un filet de courant ou filet liquide est le liquide contenu dans un
tube de courant dont la section tend vers zéro.

V
ds

Figure 4 : Un filet liquide

11
 Le filet liquide est l’élément infinitésimal à partir duquel sont établies les équations que
nous étudions en hydraulique, ainsi que leur simplification. Lorsque plusieurs filet liquides sont
juxtaposés, on obtient ce qu’on appelle le courant liquide ou encore veine liquide

2.4 Notion de charge hydraulique, théorème de Bernoulli

La charge hydraulique dans une section est l’énergie totale (au sens mécanique du
terme), composée d’une énergie potentielle de pesanteur, d’une énergie de pression (pression
hydrostatique) et d’une énergie cinétique. Elle est définie à constante près car la référence des
énergies potentielles est toujours arbitraire. La charge hydraulique (comme d’ailleurs chacune
de ses composantes) est généralement rapportée à l’unité de poids et s’exprime, de ce faite en
longueur, l’unité est donc le mètre (m). Géométriquement, la charge hydraulique est comptée
verticalement (Figure 5) et a pour expression :
.
2
1
=I est la pente de fond du canal. Dans le cadre de ce cours, cette pente est supposée
faible, ce qui permet d’écrire l’expression de la charge moyenne en un point donné :

2 é

é é
Figure 5 : Charge hydraulique en un point d'une section donnée

On peut définir une charge moyenne dans une section donnée. Dans ce cas, on introduit
un coefficient α, appelé coefficient de Coriolis ou encore coefficient d’énergie cinétique, qui
permet de tenir compte de la répartition des vitesses dans la section considérée. L’expression
de la charge moyenne est
2-1

1
En réalité, , étant l’angle que fait le fond avec un plan horizontal. Mais puisque cet angle est supposé
faible, on a

12
 Le coefficient de Coriolis est défini comme suit :
1

La somme des deux premiers termes de la charge hydraulique ( ) donne ce qu’on

appelle cote piézométrique.
D’une manière générale, la charges entre totale varie entre deux section distincte à
cause des frottements sur les parois (perte d’énergie). La différence entre les deux charge est
appelée perte de charge. La pente de la ligne d’énergie n’est rien d’autre que la perte de charge
par unité de longueur ou encore perte de charge unitaire. Elle est notée J et a pour expression
2‐2

Cette définition traduit en même temps le théorème de Bernoulli, à savoir que la

variation de charge entre deux sections distinctes est égale à la perte de charge entre ces deux
sections :
∆ 2‐3

Plusieurs formules, que nous présenterons dans le chapitre suivant, permettent de

calculer la valeur de J.
2.5 Différents types d’écoulement
De façon générale, le débit et le tirant d’eau sont des fonctions du temps et de la
position à laquelle on les observe. Si les mesures montrent que ces deux grandeurs dépendent
du temps, on dira qu’on a affaire à un écoulement transitoire ou non permanent. Ce type
d’écoulement ne fait pas l’objet de ce cours. En revanche, lorsque ces grandeurs ne dépendent
pas du temps, on dira que l’écoulement est permanent. On distingue des écoulements
permanents uniformes, les écoulements permanents graduellement variés et les écoulements
permanents brusquement variés. Nous étudierons chacun de ces écoulements dans les
chapitres suivants.

13
 II‐ Ecoulement uniforme
1. Définition
Un écoulement permanent est dit uniforme lorsque les paramètres géométriques
( , , et hydrauliques , ne changent pas d’une section à l’autre. Pour un écoulement
uniforme, les filets liquides sont rectilignes, parallèles entre eux et parallèles aux parois, les
trajectoires ne se distinguent pas des lignes de courant. Les propriétés d’un type d’écoulement
sont résumées dans Mar, 2004. Attention, dans une section donnée, dans les conditions d’un
écoulement uniforme, la vitesse n’est pas constante dans toute la section, mais les vecteurs
vitesses sont parallèles entre eux.

2. Etude des écoulements uniformes

2.1 Propriétés fondamentales
1‐ De la définition de l’écoulement uniforme, on arrive à une définition synthétique qui
résume les définitions et hypothèse d’un écoulement uniforme : c’est la première
propriété fondamentale qui s’énonce comme suit : « Dans un écoulement uniforme à
surface libre, la perte de charge unitaire est égale à la pente du fond du canal elle‐
même parallèle à la surface libre : ». Tous les calculs, dans le cas d’un tel
type d’écoulement, se feront à partir de ce principe fondamental. Voir Carlier, 1972
pour La justification de cette propriété fondamentale
2‐ La deuxième propriété fondamentale peut s’énoncer comme suit : « ce type
d’écoulement ne s’observe que loin des singularités et loin des extrémités ».
2.2 Calcul de la pente d’énergie.
Dans le principe fondamental énoncé ci‐dessus, on voit que la connaissance de la perte
de charge suffit à connaitre la pente du fond ou celle de la surface libre. Mar, 2004 (1) a montré
que la relation entre le frottement contre la paroi s’exprime en fonction de la pente de la ligne
d’énergie par , ou est le rayon hydraulique dans la section considérée. En outre,
il a montré (voir également Carlier 1972) que est fonction de la vitesse moyenne U de
l’écoulement, mesurée dans cette section : . Différents auteurs ont donné une
expression à cette fonction en définissant un coefficient de rugosité, appelé souvent par abus
de langage la rugosité. Certains de ces coefficients (que nous allons définir ci‐dessous) varient
comme la rugosité, alors que d’autre varient inversement à la rugosité. Il faudra donc faire
attention dans la compréhension de cette notion de rugosité et de coefficient de rugosité.
2.3 Equations de l’écoulement uniforme
La relation devient dans le cas des écoulements uniformes en
vertu du principe fondamental énoncé ci‐dessus. La connaissance de permet d’obtenir
une relation qui permet de résoudre les problèmes qui peuvent se poser dans un écoulement
uniforme (voir plus loin).

14
 Chézy a montré que la fonction peut s’écrire sous la forme ou C est
appelé coefficient de Chézy qui tient compte de la rugosité du canal, de sa forme, mais aussi des
conditions d’écoulement.
On obtient ainsi, ce qui s’appelle la formule de Chézy qui s’écrit :

Il existe d’autres formules de ce type pour le calcul de l’écoulement, les ouvrages tels
que Lencastre, 2005 et Carlier, 1972 renseignent sur toutes les formules existantes, nous ne
citerons que celles de Manning et Strickler qui sera celles utilisées dans ce document (voir
également Mar, 2004).

En effet Manning et Strickler ont reformulé la relation de Chézy sous la forme :

/
√ , pour Strickler, étant le coefficient de rugosité de Strickler. Il est
fonction de la nature de la paroi ; il est d’ autant plus grand que celle‐ci est lisse.
/
√ , pour Manning n étant le coefficient de rugosité de Manning
On voit que ces deux formules sont identiques, avec la relation . On peut
également déduire une relation entre l’un de ces coefficients et le coefficient de Chézy :
/ /
ou encore .

2.4 Estimation des valeurs de

Le coefficient de Strickler est le paramètre le plus difficile à connaitre et la plus grande
source d’erreur. Les méthodes d’estimation du coefficient de rugosité de Strickler , et donc de
Manning , ainsi que des tables permettant de choisir des valeurs de ces coefficients en
fonctions des types de parois, pour des matériaux cohérents, sont présentées dans Mar 2004.
Dans les cas des matériaux non cohérents, on peut estimer ce coefficient à l’aide de la formule
(il y a bien d’autres, tout dépend des auteurs) :
/
26 /
Où d35 est le diamètre des grains tel que 35% en poids des particules ont un diamètre
inférieur à d35.
2.5 Estimation des valeurs de dans le cas des des sections composées
La formule de Manning‐Strickler peut s’écrire également en fonction du débit :

/
√

S étant la section mouillée. En réarrangeant cette formule en tenant compte de

l’expression du périmètre mouillée, on peut encore écrire :
/

/ √
Considérons une section mouillée composée de n autres sections ayant des rugosités
propres. Deux cas de figure peuvent se présenter :

15
 A’ H
A H
1
3
2

B C F
G

D E

Figure 6: Rugosité dans le cas d'une section composée. Par exemple, La section mouillée 1 vaut l’aire ABSA’ alors que le périmètre
mouillé 1 a pour valeur la somme des longueurs de segment AB+BC.

¾ Dans le cas des cours d’eau en crue par exemple, il est très peu probable d’avoir
les mêmes vitesses moyennes dans toutes les n sous‐sections. Dans ce cas, le
débit s’écrira tout simplement sous la forme.

/
√ /

¾ Les vitesses moyennes dans les sous‐sections sont identiques entre elles et
identiques à la vitesse moyenne de toute la section (hypothèse d’Einstein). Dans
ce cas, on calcule d’abord la rugosité équivalente puis le débit avec cette rugosité
équivalente . Pour calculer la rugosité équivalente, on écrit que la somme des
sections élémentaires est égale à la section mouillée totale:
Partant de

/
/
√ / √

On peut extraire
/

√
En faisant de la somme on a
/ /

√ √
Soit en remarquant U=U1=U2=…=Un et que

16
 On arrive à
/ /

∑
∑ / ∑ /

Le débit correspondant à cette section s’écrira donc

/
√

3. Types de problème dans un écoulement normale

3.1 Problème de type 1
C’est le cas le plus facile à étudier. Lorsqu’on connait la géométrie de notre canal et
qu’on désire calculer un des paramètres, connaissant les deux autres (comme présenté dans le
tableau ci‐dessous), il suffit d’appliquer directement la formule de Manning‐Strickler.

On connait On cherche
I Q
Q I
Q I

Tableau 1 : Problème de type 1

Le calcul induit par ce type de problème est facile, il suffit simplement de poser
l’équation de Manning‐Strickler et d’extraire l’inconnue cherchée. La plus part des calculatrices
ayant des fonctions puissance peuvent faire ce type de calcul.

3.2 Problème de type 2 calcul de la profondeur normale

Ce type de problème se pose dans les deux cas suivants :
¾ lorsque que La valeur limite de la hauteur normale est choisie, on cherche la
forme à donner.
¾ Lorsque la forme du canal est connue et qu’on cherche la position de la ligne
d’eau, on dit qu’on cherche la profondeur normale, ou encore hauteur normale
Dans les deux cas, le calcul est itératif sauf dans des cas où la géométrie est simplifiée.
L’étude de ce type de cas se résume au calcul de la profondeur normale. En effet, la
méthode utilisée pour calculer cette profondeur normale peut s’appliquer au calcul de la forme.

3.2.1 Notion de profondeur normale

La profondeur normale pour un débit Q donné et une pente I donnée est la
profondeur de l’écoulement en régime uniforme. Dans ce cas, la surface libre représente le
niveau normal. Le rayon hydraulique, la section mouillée et le périmètre mouillé, correspondant
à la profondeur normal, sont appelés les éléments normaux.

17
Figure 7 : Abaque de détermination de la profondeur normale cas des sections circulaires et trapézoïdales
18
 Pour deux pentes différentes les courbes de la profondeur normale en fonction du débit
présentent les caractéristiques de la figure :

Figure 8 : Variation de la profondeur normale pour une pente donnée, en fonction du débit. Influence de la pente

Trois méthodes peuvent être appliquées pour la détermination de la profondeur

normale.

3.2.2 Détermination de la profondeur normale

[Link] Méthode des abaques de recherche de la profondeur normale
Les canaux les plus courants ont souvent des formes trapézoïdales (de préférence isocèle
s’il n’y a pas de contrainte particulière) ou circulaires, ou encore rectangulaire.
Pour des canaux ayant des sections trapézoïdales (les sections rectangulaires et
triangulaires sont des sections trapézoïdales particulières) ou circulaires, des abaques ont été
établies afin de lire directement la profondeur normale à partir des caractéristiques de la
section Figure 7.
¾ Pour les sections circulaires, il suffit de calculer la quantité

De la projeter sur la courbe (graduation logarithmique) correspondant à la section

circulaire pour obtenir la valeur de . Connaissant le diamètre du canal, on en déduit la
profondeur normale
¾ Pour les sections trapézoïdales, on calcule la quantité

On projette le résultat sur la courbe correspondant au fruit de berge de la section

trapézoïdale considérée pour obtenir la valeur de on en déduit la profondeur normale.

19
 NB : Cette méthode n’est pas applicables aux sections composées (qui ne sont pas
purement circulaires ni purement trapézoïdales).
[Link] Méthodes de la débitance

[Link].1 Méthode graphique

Cette méthode s’applique à n’importe quelle section dont la géométrie est bien connue.

√
Figure 9 : Méthode graphique utilisant la débitance pour la détermination de la profondeur normale.

La première méthode dite méthode de la débitance, part de l’égalité suivante, issue de la

formule de Manning‐Strickler :
/
√

On peut remarquer que la partie gauche de l’égalité est fixe, puisque le débit et la
√
/
pente sont entièrement connus. La deuxième partie, , appelée débitance, en
revanche est variable. Le but de la méthode est de construire la courbe de la débitance en
fonction du tirant d’eau . La profondeur normale est obtenue en projetant la débitance
particulière sur la courbe obtenue (Figure 9).
√

[Link].2 Méthode d’exploration itérative

La méthode graphique est très approximative, en ce sens qu’on trace la courbe en
quelques point puis on utilise des projections sur cette courbe qui au départ était imprécise. Il y
a possibilité alors d’être plus précis, en faisant le calcul en quelques points, mais de façon
optimale par exploration. Il s’agit dans un premier temps de calculer le rapport , puis de partir
√
d’une valeur de y, d’incrémenter par pas assez grands au départ tout en essayant d’aller dans le

20
sens de la convergence de la débitance vers la valeur du rapport . Il faudra réduire ce pas au
√
fur et à mesure qu’on s’approche de la solution pour assurer la précision. L’exemple (Tableau 3
à Tableau 6) des exercices d’application illustre parfaitement cette méthode qui peut être
facilement mise en œuvre sous Excel.
[Link] Méthode d’identification des pentes
La deuxième méthode part de l’égalité, elle aussi obtenue à partir de l’équation de
manning‐Strickler :

D étant la débitance définit ci‐dessus. La partie de gauche de l’égalité ne change pas,

alors que la partie de droite varie en fonction de la profondeur alors que la pente du canal étant
fixée une fois pour toute. On trace la courbe pour différentes valeurs de y, puis on
cherche le point d’intersection de cette courbe avec la droite horizontale . Cette
intersection n’est rien d’autre que la profondeur normale (Figure 10).

Figure 10:Autre méthode graphique de détermination de la profondeur normale

21
[Link] Méthode numérique de calcul de la profondeur normale.
Aujourd’hui, il existe des calculatrices capables de résoudre numériquement les
équations implicites par itération. Avec l’avènement des ordinateurs, il est tout simplement
recommandé de d’écrire de petits programmes pour la résolution de ce type d’équations.

[Link] Exemple d’application :

[Link].1 Cas d’une section trapézoïdale pure

On considère un canal de section trapézoïdale ayant pour largeur au radier b=20 m, avec
un fruit des berges m=1. La rugosité des parois de ce canal est de 70 m1/3/s et la pente de fond
est de 0.1%. Calculer la profondeur normale de ce canal, respectivement par la méthode
graphique et par utilisation d’abaque, lorsqu’on écoule un débit Q=170m3/s dans ce canal.
Solution
• Méthode graphique

/
L’expression de la débitance est .

Tableau 2 : Méthode graphique de détermination de Yn Canal trapézoïdal

/
S P
0 0 20 0 0
0.5 10.25 21.41 0.48 439.04
1 21.00 22.83 0.92 1390.42
1.5 32.25 24.24 1.33 2730.62
2 44.00 25.66 1.71 4412.83
2.25 50.06 26.36 1.90 5373.76
2.5 56.25 27.07 2.08 6411.58
3 69.00 28.49 2.42 8711.63
3.5 82.25 29.90 2.75 11303.56
4 96.00 31.31 3.07 14181.66
4.25 103.06 32.02 3.22 15726.94

0
0 5000 10000 15000

Figure 11 : Méthode graphique de détermination de Yn Canal trapézoïdal

22
 170
5375.87202
√ √0.001
En prenant différentes valeurs de y, on obtient le tableau ci‐dessus. La courbe obtenue () permet
d’estimer, pour la valeur de 5375.87202, 2.25
√

• Méthode d’exploration itérative

Tableau 3 : Méthode d’exploration itérative de détermination de Yn Canal trapézoïdal

y S P R D √
1 21 22.8284271 0.91990569 1390.42064

2 44.00 25.66 1.71 4412.83

3 69.00 28.49 2.42 8711.63 5375.76
2.5 56.25 27.07 2.08 6411.58
2.25 50.06 26.36 1.90 5373.76
2.251 50.09 26.37 1.90 5377.76
2.2505 50.07 26.37 1.90 5375.76
On converge avec une précision meilleure que ci‐dessus, Le tout reste joué sur le
nombre de chiffres après la virgule que l’on s’autorise
• Méthode des abaques
•
On calcule la quantité
170
0.026
√ 70 √0.001 20

En projetant cette valeur sur la courbe correspondant à m=1, on obtient environ

0.2

On en déduit que

20 0.2 2.4

• Méthode numérique
Avec un petit logiciel développé par le 2iE (Biaou, 2008), o a les résultats présentés sur
l’interface de la Figure 12.

23
 Figure 12 : Yn Canal trapézoïdal Méthode numérique, utilisant le logiciel Hydraulique hydrau_Gene du 2iE (Biaou, 2008)

[Link].2 Cas d’une section circulaire

Calculer, par diverses méthodes, la profondeur normale du canal circulaire de diamètre D= 1m, de pente
I=0.1%, de fruit des berges KS=70 qui écoule un débit Q=128l/s.

Méthode graphique

/
. Avec θ sin , , avec 1 cos 2

Ici, c’est la valeur de θ qu’on fait varier et non la valeur de y. En effet, il existe une
relation entre y et θ. En fixant la valeur de θ, on déduit la valeur de y par cette relation. Le
principe de détermination de la profondeur normale reste le même que celui du ca s du canal
de section trapézoïdale.
Le tableau de calcul ainsi que la figure ci‐dessous permettent d’estimer à 0.29 la
profondeur normale et d’en déduire l’angle normal à 2.82 radian

24
 Tableau 4 Méthode d’exploration itérative de détermination de Yn Canal circulaire

θ y S P R D
0.1 0.00062487 2.08229E‐05 0.05 0.00041646 8.1287E‐06
0.5 0.01554379 0.002571808 0.25 0.01028723 0.00851534
0.9 0.04977645 0.014584136 0.45 0.03240919 0.10377456
1.3 0.1019581 0.042055227 0.65 0.06470035 0.47444853
1.7 0.17000843 0.088541899 0.85 0.10416694 1.37214542
2.1 0.25121448 0.154598829 1.05 0.14723698 3.01751086
2.282 0.29165717 0.190552909 1.141 0.16700518 4.04514022
2.9 0.43974862 0.332593834 1.45 0.22937506 8.72388022
3.3 0.53956044 0.432218212 1.65 0.26195043 12.3864593
3.7 0.63779512 0.528729518 1.85 0.28579974 16.0585434

0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
‐5 0 5 10 15 20

Figure 13 : Application du calcul de yn par la débitance pour un canal circulaire

• Méthode d’exploration itérative

On utilise les expressions de sections et périmètre mouillé ci‐dessus.

Tableau 5 : Méthode d’exploration itérative de détermination de Yn Canal circulaire

θ y S P R D √
1 0.06120872 0.019816127 0.5 0.03963225 0.16124367
2 0.22984885 0.136337822 1 0.13633782 2.52808667
3 0.4646314 0.357359999 1.5 0.23824 9.61347566
2.5 0.34233882 0.237690982 1.25 0.19015279 5.50189806 4.04771541
2.3 0.29575628 0.194286848 1.15 0.16894509 4.1562834
2.28 0.29120275 0.190139911 1.14 0.1667894 4.03289534
2.282 0.29165717 0.190552909 1.141 0.16700518 4.04514022
Puis on déduit que 0.29 et que 2.82

25
 • Méthode des abaques
On calcule la quantité
0.128
0.058
√ 70 √0.001 1

En projetant cette valeur sur la courbe correspondant au canal circulaire, on obtient environ

0.3

On en déduit que 1 0.3 0.3

• Méthode numérique
Ici aussi, c’est l’application du 2iE (Biaou, 2008 qui est utilisée). Le résultat des calculs est
présenté sur la figure suivante :

Figure 14 : Canal circulaire, méthode numérique de calcul de la profondeur normale

26
[Link].3 Cas d’une section composée
On considère un canal de pente 0.1%, avec un coefficient de Manning n=0.02. la partie
trapézoïdale isocèle a un fruit de berge
berge m=0.5. Ce canal véhicule un débit 8 / . On
donne 1.49 et 1.84 Calculer la profondeur normale de ce canal.

r=1m r=1m

Figure 15

Il faut au départ remarquer et insister sur le fait que pour les sections composées, on ne
peut plus utiliser la méthode des abaques. La méthode numérique reste valable si les
expressions littérales des ections et périmètres mouillés sont bien définies (par exemple par le
solveur de la TI, le solveur intégré sous Excel) L’application du 2iE dans ses versions supérieures
devrait permettre d’introduire ces expressions et d’en déduire cette profondeur. La méthode la
plus rapide (sous Excel ou tout autre tableur) est la méthode d’exploration itérative.

• Méthode d’exploration itérative

Donnons dans un premier temps les expressions de caractéristiques géométriques en
remarquant que la figure peut se décomposer en deux sous figures :
Section mouillée
2
Commentaire de ce calcul : La première expression représente l’expression de la section
de la partie trapézoïdale. La deuxième est celle de la partie rectangulaire du canal. Il est facile de
remarquer que cette partie a pour dimensions : hauteur et largeur 2
Périmètre mouillé
2 1 2 2
Commentaire du calcul du périmètre mouillé :
La première partie de cette expression représente l’expression du périmètre mouillé de
la partie trapézoïdale. La deuxième (2 ) est l’expression des deux décrochements en fin de
section trapézoïdale. La troisième est la somme des deux hauteurs de la partie rectangulaire.
NB : Il faut remarquer que c’est juste les parties physiques effectivement en contact avec
canal qui on été prises en compte dans l’expression du périmètre mouillé en conformité avec la
définition du périmètre mouillé.
Les résultats de calcul, renseignés sur le tableau ci‐dessous donnent 2.1

27
 Tableau 6 : Méthode d’exploration itérative de détermination de Yn Canal composé

y S P R D √
1 1,97005 6,19174129 0,31817382 45,908195
2 5,81005 8,19174129 0,70925702 231,039181
2,5 7,73005 9,19174129 0,84097776 344,356788 252,982213
2,2 6,57805 8,59174129 0,76562478 275,26088
2,1 6,19405 8,39174129 0,7381126 252,945209
3.3 Problème de type 3 forme optimale : section hydrauliquement favorable
Ce problème se pose lorsque l’on dispose d’un débit fixé et que l’on désire dimensionner
le canal (largeur et hauteur d’eau correspondant à ce débit) de manière à évacuer ce débit, pour
une pente I et une rugosité données. Ce problème est donc à deux degrés de liberté dont la
résolution consiste à résoudre un problème d’optimisation. La section qu’on obtient par la
résolution de ce problème s’appelle la « Section hydrauliquement favorable » ou encore section
économiquement favorable (Mar, 2004). Cet optimum sera obtenu pour périmètre mouillé
minimisé avec une section mouillée. Se référer à Mar, 2004 pour les propriétés de ce type de
problème.

3.3.1 Section hydrauliquement favorable de forme circulaire

De tous les canaux, les sections circulaire et demi‐circulaire, de part leur géométrie, sont
celles qui offrent pour une section mouillé donnée, le plus petit périmètre. De ce fait, ces
sections sont naturellement hydrauliquement favorables. Ceci dit, ces formes ne sont pas faciles
à mettre en œuvre. Pour les canaux en terre par exemple, les risques d’éboulement et la
dégradation des parois sont à craindre pour ce type de canal. C’est pour cela qu’on fait recours
généralement aux sections trapézoïdales de préférence isocèle.

3.3.2 Section hydrauliquement favorable de forme trapézoïdale

Pour déterminer les relations existant entre les dimensions (largeur et tirant d’eau)
permettant de réaliser un tel type de problème, on écrira que le périmètre mouillée est
minimum, comme la section mouillée également. Cela revient à écrire que les dérivées
partielles de ces deux grandeurs sont nulles.
Pour une section trapézoïdale, on a et 2 √1 , on a ainsi

0 2 0
0 2 1 0

Ce système n’a de solution différente de la solution triviale (0,0) que si le déterminant du

système est nul. Ainsi, on a
2 1 2 0
Soit

28
 2 1
A partir de cette relation entre la largeur au radier et le tirant d’eau, on peut déduire
celles de la surface mouillée, du périmètre mouillé et du rayon hydraulique :
2 1 2 1
En posant 2√1
S s’écrit tout simplement

Par ailleurs,
2 1 2 1 2 1 2 2 1
Soit finalement
2
Le rayon hydraulique s’en déduit immédiatement :

2
Remarque : En adoptant la même démarche pour une section rectangulaire et pour une
section demi‐circulaire, on arrive au même résultat, à savoir que .
/
/ / /
√ √ 2 √
2
Soit
/ /
2
√
3.3.3 Section hydrauliquement favorable de forme trapézoïdale : construction graphique

Figure 16 : Section hydrauliquement favorable

En considérant la section trapézoïdale ci dessus, dans le cas de la section hydraulique

ment favorable, on peut écrire que

2 1

Remarquons d’une part que

29
 1
tan

D’où on tire
D’autre part, on a

D’où l’on tire

1
Enfin on a
En revenant à l’expression de b ci‐dessus :
2 1
On peut écrire que
1
2
Soit

Ou encore

Ce qui veut dire tout simplement que le triangle OBC est isocèle en B, ce qui fait que les
hauteurs CF et OK sont égales.
Pour les propriétés de sections hydrauliquement favorables, se référer à Mar, 2004

Il existe d’autres formes de problèmes de type 3. Il s’agit des problèmes de vitesse imposé (une
trop faible vitesse peut entraîner des problèmes de dépôt et une trop forte vitesse peut
entraîner des problèmes d’érosion, surtout dans le cas des canaux en terre) et de pente limite.
Pour comprendre la résolution de ces types de problèmes, se référer à Mar, 2004

30
 III‐ Ecoulement graduellement varié
1. Définition
Nous avons vu qu’un écoulement uniforme est caractérisée par une constance de tous
les paramètres qui le définissent. Ceci dit, lorsque nous avons par exemple un changement de
pente, partant d’une pente (Figure 8) les hauteurs normales sont différentes (plus faibles
dans le cas de la pente forte que dans le cas de la pente faible, pour un débit donné). Il est
important de préciser que ces hauteurs normales ne sont observées que loin en amont et loin
en aval du changement de pente (appelée encore singularité). Il est important de savoir ce qui
se passe au voisinage de la singularité. Dans ce voisinage, nous obtenons des lignes d’eau
particulières (que nous allons étudier plus loin), en fonction du type de singularité. Cette ligne
d’eau ne respecte plus les hypothèses d’un écoulement uniforme, on arrive dans le cas des
écoulements variés. Pour ces écoulements, les pentes de la surface libre, de la ligne d’énergie, la
hauteur d’eau et donc section mouillée varie d’une section à l’autre. Selon l’allure du
changement observé, on peut avoir les écoulements graduellement variés et les écoulements
brusquement variés. Dans ce chapitre, nous n’étudions que les premiers cas.

2. Propriétés des Ecoulement graduellement variés (EGV)

Comme introduit ci‐dessus, dans un EGV, on assiste à une variation (lente) et continue
de la section mouillée, des pentes d’énergie et de surface libre, de la profondeur d’eau et donc
de la section mouillée. Les frottements ne sont donc plus négligeables et on assiste ainsi à ce
que l’on appelle la perte de charge (ou encore perte d’énergie) dans le canal. Les vitesses sont
considérées en première approximation comme parallèles entre elles. La répartition des
pressions restent, dans ce cas, hydrostatique dans la section, comme c’est le cas dans un
écoulement uniforme, ce qui entraîne que La ligne piézométrique se confond à la surface libre.
Ici, la ligne d’eau, la ligne d’énergie et le fond du canal ne sont plus parallèles. Le mouvement
est permanent (constance des paramètres dans le temps) mais on assiste à une variation des
paramètres dans l’espace.

3. Charge spécifique
3.1 Définition
Nous avions défini plus haut la notion de charge hydraulique moyenne notée H :

2
Tout au long de ce chapitre, nous allons donner une valeur 1 au coefficient de Coriolis .
Il faut rappeler que la charge moyenne est définie par rapport à une référence arbitraire.
Ainsi, pour deux références quelconques, les charges calculées ne sont pas les mêmes. C’est
pour cela qu’on introduit une charge particulière qui est définie par rapport à une référence

31
particulière, c'est‐à‐dire le fond du canal. Cette charge moyenne est appelée la charge
spécifique :

2 2
3.2 Etude de la charge spécifique
La charge spécifique s’écrit de deux manières différentes, la première

Et la conséquence de cette dernière :

2
Ce chapitre sera consacré à l’étude de cette grandeur sous toutes ses formes. Dans un
premier temps, l’étude de cette charge spécifique en fonction de la profondeur critique, pour
un débit donné. Cette profondeur critique avec la profondeur normale, permettra d’étudier les
lignes d’eau plus loin. Nous allons par ailleurs étudier la variation du débit en fonction de la
profondeur critique pour une charge spécifique donnée.

3.2.1 Etude de la variation de la charge spécifique en fonction de la profondeur

Pour étudier la variation de la charge spécifique, nous allons utiliser son expression

Dont la dérivée par rapport à s’écrit

1
1 1
2

Figure 17 : Estimation de la variation de la section mouillée en fonction du tirant d'eau

Cette dérivée s’annule en un point donné qui a des propriétés particulières. La charge
spécifique admet un optimum en ce point.

32
 Lorsque y tends vers l’infini, tend vers l’infini et donc tend vers 0, ce qui fait
que
lim ∞
De même, lorsque y tend vers 0, S tend vers 0 et
lim ∞
2
ce qui fait que
lim ∞

Par ailleurs,
lim lim 1 1
2
Suggère que cette fonction admet pour asymptote, à l’infini, la première bissectrice.

lim lim 1 ∞
2
Donc la fonction admet pour asymptote l’axe des charges spécifiques lorsque le tirant
d’eau tend vers 0.

Q 0

45°

Figure 18: Courbe en fonction du tirant d’eau pour un débit donné, constant, dans une section donnée

33
 En résumé, et au vu de ce qui précède, la charge spécifique est une fonction qui admet
une profondeur minimum qu’on appelle profondeur critique au point d’abscisse telle
que
1 0
L’équation ci‐dessus définit ce qu’on appelle la « condition critique », celle qui
caractérise la profondeur critique de l’écoulement.
La fonction pour un débit donné tend vers l’axe des ordonnées lorsque y
tend vers 0 et tend vers la première bissectrice lorsque y tend vers l’infini. On obtient alors la
représentation de la Figure 16. Cette courbe permet de montrer que pour un débit donné et
pour une charge spécifique , supérieure à la charge critique, on a deux tirants d’eau.

¾ Pour les valeurs de y inférieures à la profondeur critique, la courbe est

décroissante et donc la dérivée 1 de la charge spécifique est négative.
Cette plage de valeur du tirant d’eau, on dira que le régime torrentiel ;
¾ Pour les valeurs de y supérieures à la profondeur critique la courbe est
croissante, ce qui fait que la dérivée est positive. Pour toutes ces valeurs de tirant
d’eau, on dira que le régime fluvial.
Les deux profondeurs sous les quels peut s’écouler ce débit unique, fixé pour une charge
donnée sont appelées des profondeurs conjuguées.

3.2.2 Classification des régimes d’écoulement : le nombre de Froude

Le nombre de Froude (Noté Fr) est un nombre sans dimension dont le carré représente
le rapport de l’énergie cinétique du liquide en mouvement à l’énergie potentielle de pesanteur.
Son expression est.

Dans cette expression, U est la vitesse moyenne de l’écoulement et , la profondeur

moyenne, c'est‐à‐dire
Expérimentalement, le nombre de Froude est le rapport entre la vitesse moyenne de
l’écoulement U et la célérité des ondes infinitésimales
En remplaçant cette profondeur moyenne ainsi que l’expression de la vitesse en fonction
du débit dans le l’expression du nombre de Froude, on obtient

Soit

On conclut donc que la condition critique correspond également à

34
 A partir de cette transformation, on peut classifier les régimes d’écoulement avec le
nombre de Froude :
¾ Le régime est dit torrentiel ou encre supercritique lorsque
1 0 1
Pour ce type de régime, les hauteurs d'eau sont généralement faibles et les
vitesses fortes
¾ Le régime est dit fluvial ou encore subcritique lorsque
1 0 1
Pour ce type de régime, les hauteurs d'eau sont généralement fortes et les
vitesses faibles

¾ Le régime est dit critique lorsque

1 0 1

3.2.3 Variation du tirant d’eau en fonction du débit, lorsque est fixé :

En reconsidérant l’expression de la charge spécifique

2
On peut écrire
2
La dérivée de cette fonction s’écrit :
1 2
2 . .
2 2
On avait vu que

Puisque
2
Alors

Le maximum du débit est obtenu lorsque cette dérivée est nulle , c'est‐à‐dire lorsque
0
Ou encore
1 0
Ce qui correspond bien à la fonction critique.

35
 En somme, la courbe
2

a les propriétés suivantes :

¾ Pour un tirant d’eau nul, on a un débit nul, puisque la surface est nulle ;
¾ Le débit est nul également lorsque le tirant d’eau prend la valeur de la charge
spécifique. En effet, signifie que l’énergie cinétique est nulle donc que le
débit est nul ;

Fluvial

Torrentie

Figure 19 : Variation du débit en fonction de la profondeur

¾ Il existe un tirant d’eau pour lequel le débit devient maximum. Ce tirant d’eau
n’est rien d’autre que le tirant d’eau critique défini ci‐dessus et le débit
correspondant est le débit critique ;
¾ Pour une charge spécifique donnée, on peut écouler un débit Q donné sous deux
tirants d’eau, appelés profondeurs conjuguées. La première, 1 correspond
au régime torrentiel et la deuxième, 2 correspond au régime fluvial ;
¾ Pour une charge constante donnée, H, une augmentation de débit se traduit par
une élévation de la hauteur d’eau dans en régime torrentiel, et par un
abaissement de la hauteur d’eau en régime fluvial.

4. Etude du régime critique

Le régime critique est par définition le régime d’écoulement pour lequel le tirant d’eau
est critique. Il correspond à un débit maximum pour une charge spécifique donnée et à une
charge spécifique minimum pour un débit donné. Il est caractérisé par un nombre de Froude
égal à l’unité. L’étude du régime critique consiste principalement à déterminer les éléments
critiques de l’écoulement, ce qui se résume généralement à la détermination de la profondeur
critique. Comme dans le cas de la détermination de la profondeur normale, ici aussi, trois
méthodes peuvent être utilisées pour la détermination de la profondeur critique.

36
4.1 Méthode des abaques.
Pour les sections trapézoïdales (dont un cas particulier est la section rectangulaire) et
circulaires, il est possible de déterminer la profondeur critique à partir d’un abaque
préalablement établi (Figure 18)
¾ Pour les sections circulaires, calculer la quantité

et la projeter sur la courbe (en lisant sur l’axe du bas en graduation logarithmique), de la
section circulaire pour lire (en graduation logarithmique) la valeur de . Connaissant le
diamètre du canal, on en déduit la profondeur normale ;
¾ Pour les sections trapézoïdales, calculer la quantité

et la projeter sur la courbe (en lisant sur l’axe de dessus de l’abaque en graduation
logarithmique) qui correspond au fruit des berges du canal pour lire(en graduation
logarithmique) .
Il est à remarquer que la section rectangulaire est un cas particulier de section
trapézoïdale (fruit de berge nulle) et donc que la profondeur critique peut être déterminée à
l’aide de l’abaque, pour m=0. En revanche, à cause du paramètre b au dénominateur (à cause
de la méthode de construction de l’abaque), l’autre cas particulier, c'est‐à‐dire la section
triangulaire (qui s’obtient pour b=0) peut utiliser cette méthode pour déterminer la profondeur
critique.

4.2 Méthode graphique

La condition critique étant donnée par
1
On peut écrire que

D’une part, est une constante, puisque le débit est fixé et que g est une constante.
√
D’autre part, est une grandeur qui est fonction de y. La profondeur critique est donc le
l’intersection de la partie fixe avec la partie variable compte tenue de l’égalité ci‐dessus..
Ainsi, la méthode graphique consiste à tracer la courbe , puis à chercher la
profondeur (en abscisse) pour laquelle prend la valeur (Figure 19)
√

37
 1
m b

Figure 20 : Abaque de détermination de la profondeur critique pour les sections circulaires et trapézoïdales
2iE/2008/BIAOU

38
 Figure 21 : Méthode graphique pour la détermination de la profondeur critique

4.3 Méthode d’exploration itérative

Elle est identique à celle développée dans la partie [Link].2
Ainsi, on calcule dans un premier temps le rapport , puis on fait varier y avec un pas
√
assez élevé au départ, en cherchant à aller dans le sens de la convergence vers la valeur du
rapport , puis on ressert le pas au fur et à mesure qu’on s’approche de la solution pour aller
√
chercher la précision qu’il faut pour la solution.
Cette variante est expliquée dans l’exercice d’application ci‐dessous
4.4 Méthode numérique
La condition critique
1
Est une équation implicite qui ne peut être résolu directement, de façon analytique que
pour certaine section simple. Sa résolution requiert donc un solveur, ce qui est possible
aujourd’hui avec la généralisation de l’informatique. La méthode numérique consiste donc à
utiliser ces solveurs ou des programmes implicites.

39
4.5 Exercice d’application
4.5.1 Cas d’une section trapézoïdale pure
En utilisant les diverses méthodes, calculer la profondeur critique du canal de la partie
3.2. du chapitre 2
Solution
• Méthode graphique
On calcule la quantité

170
54.28
√9.81

Les résultats des calculs sont présentés dans le tableau ci‐dessous

S L
0.001 0.02 20.00 0.00 0.00
1 21.00 22.00 0.95 20.52
1.75 38.06 23.50 1.62 48.44
1.88 41.13 23.76 1.73 54.12
2.5 56.25 25.00 2.25 84.38
2.75 62.56 25.50 2.45 97.99
3 69.00 26.00 2.65 112.41
3.25 75.56 26.50 2.85 127.60
3.5 82.25 27.00 3.05 143.56
4 96.00 28.00 3.43 177.76
4.25 103.06 28.50 3.62 195.99

200,00
160,00
120,00
80,00
40,00
0,00
0 1 2 3 4

Figure 22 Exemple de profondeur critique par la méthode de débitance

On trouve 1.88

40
 • Variante de la méthode graphique

Tableau 7 : Variante de la méthode graphique

S L
1.000 21.000 22.000 0.955 20.517
2.000 44.000 24.000 1.833 59.576
1.500 32.250 23.000 1.402 38.188 54.277
1.800 39.240 23.600 1.663 50.598
1.850 40.423 23.700 1.706 52.791
1.882 41.182 23.764 1.733 54.213
Ce tableau montre que la profondeur critique est estimée à 1.882 m
• Méthode des abaques

On calcule

170
0.03
√9.81 20

En projetant cette valeur (attention ! par le haut, car section trapézoïdale) sur la courbe qui correspond à
m=1 (attention, ici, il s’agit d’une section trapézoïdale, alors la lecture de U se fait sur la graduation
de dessus), on obtient alors,

0.095 0.095 20 1.90

• Méthode numérique
L’application utilisée pour cette méthode est celle développée par le 2iE (Biaou, 2008). La solution du
problème est présentée sur l’interface ci‐dessous.

On peut remarquer méthode d’exploration itérative est aussi précise que la méthode numérique. Elle
requiert cependant un temps de calcul plus élevé. Les méthodes graphiques et des abaques sont assez
rapides, mais moins précises. Elles permettent cependant d’avoir une idée assez rapide de la profondeur
critique, comme c’est le cas d’ailleurs pour la profondeur normale.

41
 Figure 23 : Méthode numérique de calcul de profondeur critique par l’application Hydrau‐Géné (Biaou, 2008)

4.5.2 Cas d’une section circulaire

Calculer, par diverses méthodes, la profondeur normale du canal circulaire de diamètre D= 1m, de pente
I=1%, de fruit des berges KS=70 qui écoule un débit Q=200l/s.

Les expressions des différents paramètres géométriques sont

θ sin
8

2
2

A partir de ces différentes expressions, on peut appliquer les différentes méthodes :

42
 • Méthode numérique

S L
0.1000 0.0006 0.0000 0.0500 0.0004 0.0000
0.5000 0.0155 0.0026 0.2474 0.0104 0.0003
1.0000 0.0612 0.0198 0.4794 0.0413 0.0040
1.5000 0.1342 0.0628 0.6816 0.0922 0.0191
2.0000 0.2298 0.1363 0.8415 0.1620 0.0549
2.0900 0.2490 0.1527 0.8649 0.1766 0.0642
3.0000 0.4646 0.3574 0.9975 0.3583 0.2139
4.0000 0.7081 0.5946 0.9093 0.6539 0.4808
4.5000 0.8141 0.6847 0.7781 0.8800 0.6423
5.0000 0.9006 0.7449 0.5985 1.2446 0.8310
5.5000 0.9622 0.7757 0.3817 2.0324 1.1058

0,5000

0,4000

0,3000

0,2000

0,1000

0,0000
‐0,1000 0,0000 0,1000 0,2000 0,3000 0,4000 0,5000 0,6000

On en déduit la valeur de la profondeur critique : 0.25

• Méthode d’exploration itérative

Ici également, on s’appuie sur les expressions ci‐dessus

S L
1.0000 0.0612 0.0198 0.4794 0.0413 0.0040
2.0000 0.2298 0.1363 0.8415 0.1620 0.0549 0.0642
2.1000 0.2512 0.1546 0.8674 0.1782 0.0653
2.0900 0.2490 0.1527 0.8649 0.1766 0.0642
On en déduit la valeur de 0.249 pour un angle critique estimé à 2.09

43
 • Variante des abaques
0.2
0.0.064
√9.81 1

En lisant sur l’abaque (attention, ici il faut lire sur la graduation du bas.), on en déduit

0.25 Soit 0.25

• Variante numérique

44
4.5.3 Cas d’une section composée

On considère un canal de pente 0.1%, avec un coefficient de Manning n=0.02. la partie

trapézoïdale a un fruit de berge m=0.5. Ce canal véhicule un débit 8 / . On donne
1.49 et 1.84 Calculer la profondeur normale de ce canal.

r=1m r=1m

Ici aussi, il faut insister sur le fait que dans le cas de la profondeur critique pour les
sections composées, on ne peut plus utiliser la méthode des abaques. La méthode d’exploration
itérative est la plus indiquée (Sous Excel).
Section mouillée
2
Commentaire de ce calcul : La première expression représente l’expression de la section
de la partie trapézoïdale. La deuxième est celle de la partie rectangulaire du canal. Il est facile de
remarquer que cette partie a pour dimensions : hauteur et largeur 2

Largeur en miroir
2 2
Commentaire de ce calcul : La largeur en miroir est obtenue en additionnant celle de la
section trapézoïdale 2 aux deux revanches2 .

S L
1 1.24283612 2.35152568 0.5285233 0.90353699
2 6.57253514 4.35152568 1.51039788 8.07753038 2.57895697
1.5 3.90768563 3.35152568 1.16594232 4.21947241
1.2 2.30877593 2.75152568 0.83908936 2.11488079
1.269 2.67652516 2.88952568 0.9262853 2.57598729
On arrive à 1.269

45
4.6 Notion de pente critique
La pente critique est par définition une pente pour laquelle un débit donné s’écoule à
une profondeur normale critique, de façon uniforme. Sa détermination obéit donc
simultanément à la condition critique et à la formule de Manning Strickler avec les paramètres
du régime critique ( :, ):

1
/ /

Se référer à Mar, 2004 Carlier, 1972 pour plus d’information sur la pente critique

5. Les courbes de remous

5.1 Equation général des écoulements graduellement variés
C’est la courbe décrite par la surface libre (par rapport au fond du canal) lorsqu’on
parcourt le canal (ou le cours d’eau) sur sa longueur, par son axe par exemple. C’est la courbe
donnée par , x étant la position du canal à partir d’une référence sur son cours.
La fonction n’est pas connue explicitement, mais est définie par une équation
différentielle :

La dérivée du dénominateur est connue (cf III‐3.2.1) :

Ligne d’eau

Figure 24 : Définition de la courbe de remous ou encore ligne d’eau

1
Pour le numérateur, on part de la définition de :

46
 Alors pour les canaux prismatiques
(2)
Au final, l’équation différentielle recherchée s’écrit

1
J étant donné par la formule de Manning‐Strickler
En regardant de près cette formule, on peut faire des observations suivantes :
‐ Il serait difficile d’attribuer une expression analytique à la ligne d’eau, sauf pour
certaines sections suffisamment simple qui permettent d’avoir des primitives
ordinaires pour l’équation différentielle ci‐dessus, mais la ligne d’eau est
entièrement bien définie à partir de sa dérivée et d’une section de contrôle (voir plus
loin).
‐ La fonction caractéristique de la ligne d’eau n’est pas dérivable au point critique,
puisque cette dérivée s’annule en ce point. C’est d’ailleurs ce qui rend ce point
critique !
Pour les canaux non prismatiques, se référer à Mar, 2004 pour plus d’information
5.2 Etude analytique des courbes de remous
Cette étude se fait en s’intéressant aux positions relatives du tirant d’eau par rapport
aux tirants d’eau particuliers et ou en étudiant la pente du canal par rapport à la pente
d’énergie et celle critique .

5.2.1 Position relative de I et J, par rapport aux positions relative de et de :

Considérons une section dans laquelle est la pente (dite positive si elle descend dans le
sens de l’écoulement et négative ou adverse dans le cas contraire) et la perte de charge
unitaire, c'est‐à‐dire la pente d’énergie, toujours positive ou nulle.
Dans la méthode graphique de détermination de la profondeur normale, on peut se
rendre compte (Figure 9) que la débitance

est une fonction croissante de .

Nous savons par ailleurs que pour la profondeur normale de l’écoulement est donnée
par la formule de Manning‐Strickler pour la pente de fond :

√
Et que pour une profondeur d’eau quelconque, cette formule s’applique avec la pente
d’énergie J :

Ainsi, pour un débit donné, on peut écrire le rapport

2
Remarquer que l’énergie étant décroissante sur le trajet. Il en est de même pour la cote du fond

47
 √
Puisque la débitance, donc , est une fonction croissante, ce qui conduit au résultat
du Tableau 8

Tableau 8: Etude comparée de I et de J en fonction de y et yn

profondeur
Débitance
pente

On peut en déduire la position relative de I et de Ic résumé dans le Tableau 9: Etude

comparée de I et de Ic en fonction de yn et ycTableau 9

Tableau 9: Etude comparée de I et de Ic en fonction de yn et yc

profondeur
Débitance
pente
Ces tableaux donne en même temps le sens de variation du numérateur de , c'est‐à‐
dire de
I J

5.2.2 Signe de =1 en fonction des positions relatives de y et de

2

Tableau 10 : Signe de 1 3
de la position relative de y et yc.

Profondeur
décroissant croissant
négatif positif
1

Nous avons étudié ci‐dessus les variations de (voir la Figure 16). Cette étude a
montré que cette fonction est décroissante pour les valeurs de y inférieures à la valeur critique
et croissante pour les valeurs de y supérieures à la valeur critique. C’est ce constat qui
permettra de donner le signe de =1 que nous résumons sur le Tableau 10

5.2.3 Etude qualitative des courbes de remous

Le tableau récapitulatif des différentes courbes de remous, avec les exemples est résumé
dans Mar, 2004.
Nous détaillons ici, en fonction des résultats qui précèdent, l’évolution de ces courbes
de remous. En fonction des positions relatives de , et , on l’indice de la courbe sera
désigné par un numéro dans une numérotation qui commence par les fortes valeurs de . Le

48
type de la courbe sera caractérisé par la première lettre de la classe à laquelle appartient le
remous. Nous adoptons la nomenclature présentée dans Mar, 2004.

1er cas : pente douce, soit Classe : Fluvial, notation : F

‐ ‐ +
1 ‐ + +
+ ‐ +
Sens de

Régime d’écoulement Torrentiel Fluvial Fluvial

Type

On peut remarquer que, quand bien même la classe générale de ce type d’écoulement
est fluvial, le régime d’écoulement peut être Torrentiel ( ). Se référer à Mar, 2004 pour les
exemples de ces type d’écoulement et dans quelles conditions on les rencontre
2ième cas : canal horizontal, , rejeté à l’infini Classe : Horizontale, notation : H

‐ ‐
1 ‐ +
+ ‐
Sens de

Régime d’écoulement Torrentiel Fluvial

Type
3ième cas : pente critique, , , Classe : Critique, notation : C

‐ +
1 ‐ +
+ +
Sens de

49
Régime d’écoulement Torrentiel Uniforme critique Fluvial
Type
4ième cas : pente forte, soit Classe : Torrentiel, notation : T

‐ + +
1 ‐ ‐ +
+ ‐ +
Sens de

Régime d’écoulement Torrentiel Torrentiel Fluvial

Type

5ième cas : pente contraire, 0, rejeté à l’infini Classe : Adverse, notation : A

‐ ‐
1 ‐ +
+ ‐
Sens de

Régime d’écoulement Torrentiel Fluvial

Type

5.2.4 Quelques illustrations

Ressaut

Figure 25 : Courbe de remous F1 (remous d'exhaussement) Figure 26 : Courbe de remous de type F3

50
 Figure 27 : Courbes de remous F2 et T2 Figure 28 : Courbe de remous de type F2 (remous d'abaissement)

Les figures ci‐dessus illustrent les conditions dans lesquelles surviennent les remous dans
le cas de la classe F et du remous de type1. Les remous dans les autres classes surviennent
pratiquement dans les mêmes conditions avec une atténuation ou amplification du phénomène
observé dans chaque type (représenté par les numéros 1 à 3). Se référer à Mar, 2004 ou à
Carlier, 1972 pour les passer en revue

5.2.5 Notion de section de contrôle

Le calcul de la courbe de remous consiste à la résolution de l’équation différentielle

1
Q, I, Ks, connus et J donné par

De façon générale, la résolution de cette équation différentielle est très difficile à cause
de sa complexité. Dans tous les cas, l’intégration est obtenue à une constante près. Il faut donc
connaître obligatoirement l’un des points de cette ligne d’eau. Ce point est l’intersection de la
ligne d’eau avec une section particulière appelée section de contrôle ou tout simplement
contrôle. La section de contrôle est matérialisée sur chacune des figures 19 et 20 et est
généralement définie à partir des propriétés hydrauliques de la singularité qui provoque le type
de remous. On peut citer comme exemple de singularité, les barrages, les déversoirs, les vannes
de fond, les changements de pentes etc….
Cette section est localisée à l’aval pour les remous de types , , , , , , pour
les quels le régime d’écoulement est de type fluvial et pour ces types d’écoulements, le calcul
de de la ligne d’eau se fera de l’aval vers l’amont. En revanche, pour les remous de types , ,

51
 , , , , pour lesquels le régime d’écoulement est torrentiel, cette section de contrôle se
trouve en amont et le calcul des lignes d’eau se fait de l’amont vers l’aval.

5.3 Calcul des lignes d’eau ou encore axes hydrauliques

La courbe de remous est le profil de la ligne d’eau dans le cas des mouvements
graduellement variés. C’est le lieu des points , solutions de l’équation différentielle

1
à partir d’un point de référence situé à l’intersection de la section de contrôle et de la
surface libre.

On peut écrire l’équation différentielle ci‐dessus sous la forme d’une équation

différentielle à variables séparées :
1

Le second membre n’étant fonction que de la variable , à condition bien sûre que la
section droite ne se modifie pas le long du canal.
En choisissant deux sections d’abscisses et , la solution de cette équation
différentielle entre ces deux section s’écrit :

Plusieurs méthodes existent pour le calcul des axes hydraulique, voir Mar, 2004, Carlier,
1972, Lencastre, 2005.

5.3.1 Méthode d’intégration graphique

Elle permet de calculer, pour un abaissement fixée à partir d’une profondeur
connue (le premier abaissement est mesuré à partir de la section de contrôle), la distance à
laquelle on observe cet abaissement, par la méthode graphique. On utilise la solution utilisant le
principe de variable séparée :

En traçant point par point la fonction dans un graphique, on obtient la distance

en estimant graphiquement l’aire entre et , soit par planimétrie ou par la méthode
des trapèzes. On n’a donc pas besoin de déterminer des points intermédiaires. Cette méthode
n’a pas besoins d’hypothèse simplificatrice rendant peu fiables les résultats.
Application
Calculer la ligne d’eau pour le canal due la partie 3.2 du chapitre 2, sachant qu’en un
point d’bscisse 0, on mesure une hauteur d’eau 3.1
On a calculé la profondeur normale ainsi que la profondeur critique de ce type de canal,
dans pour le débit considéré dans les exemples précédents. Puisque 2.25 et
52
1.88 , on en déduit que nous avons un écoulement de classe F, la section de contrôle se
trouve donc en aval et nous déroulons le calcul de la ligne d’eau de l’aval vers l’amont. La
section de contrôle étant le point 0,3.1 , on en déduit que nous avons un remous de type .
Le calcul démarre depuis l’ordonnée du M et on arrête dès qu’on s’approche de la profondeur
normale (voir le tableau ci‐dessous).
Quelques observations sont nécessaires pour la compréhension de la démarche
• La pente d’énergie est donnée par

Avec
/

Théoriquement, la valeur de l’aire estimée par la méthode des trapèzes est positive. Il
faut cependant tenir compte du fait sens de comptage des abscisses positives (d’amont en aval)
est contraire au sens d’évolution de la ligne d’eau (d’aval en amont, donc une variation
négative) en mettant un signe négatif dans le cas présent. Si le calcul se faisait pour un
écoulement torrentiel, ce signe n’est plus nécessaire.
Les résultats présentés dans représentent les abscisses relatives par rapport à
l’abscisse initiale (c'est‐à‐dire la section de contrôle).
Les résultats et les étapes intermédiaires des calculs sont présentés dans le tableau ci‐
dessous
3,5
3
2,5
2
1,5
1
0,5
0
‐2400 ‐2000 ‐1600 ‐1200 ‐800 ‐400 0

Figure 29 : Allure de la ligne d’eau ( ) déterminée par la méthode graphique à partir du tableau ci‐dessous

53
 Tableau 11 : Ligne d'eau calculée à partir de la méthode graphique

1 /

3.1 71.61 28.7681241 2.4892134 26.2 0.78981143 9206.89381 0.00065907 1198.38167 0

3.05 70.30 28.63 2.46 26.1 0.77871242 8957.80841 0.00063984 1217.04057 60.3855559
3 69.00 28.49 2.42 26 0.76683954 8711.63188 0.0006192 1238.43923 121.772551
2.95 67.70 28.34 2.39 25.9 0.75412494 8468.37121 0.00059701 1263.17715 184.31296
2.9 66.41 28.20 2.35 25.8 0.74049358 8228.03376 0.00057312 1292.03912 248.193367
2.85 65.12 28.06 2.32 25.7 0.7258623 7990.62727 0.00054738 1326.07248 313.646157
2.8 63.84 27.92 2.29 25.6 0.71013875 7756.15985 0.0005196 1366.70665 380.965636
2.75 62.56 27.78 2.25 25.5 0.69322022 7524.64004 0.00048958 1415.9444 450.531912
2.7 61.29 27.64 2.22 25.4 0.67499227 7296.07677 0.0004571 1476.68085 522.847543
2.65 60.02 27.50 2.18 25.3 0.65532712 7070.47943 0.0004219 1553.26196 598.596113
2.6 58.76 27.35 2.15 25.2 0.6340818 6847.85786 0.00038371 1652.52251 678.740725
2.55 57.50 27.21 2.11 25.1 0.61109603 6628.22239 0.00034219 1785.86388 764.700385
2.5 56.25 27.07 2.08 25 0.58618973 6411.58386 0.00029698 1973.83043 858.692742
2.45 55.00 26.93 2.04 24.9 0.5591601 6197.95361 0.00024768 2257.56862 964.477719
2.4 53.76 26.79 2.01 24.8 0.52977823 5987.34358 0.00019382 2733.28526 1089.24907
2.35 52.52 26.65 1.97 24.7 0.49778507 5779.76624 0.00013488 3690.6287 1249.84691
2.3 51.29 26.51 1.94 24.6 0.46288678 5575.23472 7.0239E‐05 6590.2089 1506.86785
2.26 50.31 26.39 1.91 24.52 0.43265334 5413.81169 1.3967E‐05 30977.3359 2258.21875

54
5.3.2 Méthodes de Bakhmeteff
[Link] Principe de la méthode
L’équation différentielle de la ligne d’eau (pour les canaux prismatiques) peut s’écrire
sous la forme

/
1 √
1 1 / 1

1 1 1
1 √

est la débitance correspondant à l’écoulement normal, donc à

√
/
est la débitance correspondant à la profondeur quelconque
Considérons , la pente pour laquelle le régime d’écoulement est en même temps
normal et critique, qui obéit à
1
1
Et à
/ 1

D’où
1

1
Il est à remarquer que ces transformations peuvent s’effectuer en utilisant la formule de
Chézy de calcul de J, c’est d’ailleurs à partir de la formule de Chézy que Bakhmeteff à établi
cette expression, on écrivant la débitance correspondante.
Ces formules sont établies pour n’importe quel type de section. Le rapport varie
peu avec la profondeur. Dans la quasi‐majorité des cas, il a été montré qu’on peu poser
, ce qui conduit à
1 1

1 1
En posant

On peut écrire

55
 1
1 1
1
En séparant les variables, on arrive à
1
1
1 1
Comme dit plus haut, varie très peut, donc sera considéré comme une constante. Mais
rigoureusement, il faut avoir à l’esprit que dépend de , et donc de . Cette hypothèse
acceptée, on peut écrire la solution de l’équation différentielle à variables séparées ci‐dessus :
1 1
1 1
1 1
Soit
1 1
1
1 1
Bachmeteff a construit une table pour une fonction portant son nom, définie par
l’intégrale
1
,
1
La solution de l’équation différentielle de la ligne d’eau peut donc se mettre sous la
forme :
1 , ,
Qui la longueur d’un bief dont les extrémités ont des profondeurs respectives et
[Link] Utilisation de la méthode.
Les étapes ci‐dessous sont à respecter pour l’utilisation de la méthode de Bakhmeteff :
• Calculer les profondeurs normales , critique et mesurer la profondeur de
la section de contrôle. A partir de ces trois valeurs, préciser le type d’écoulement.
• Calculer l’exposant hydraulique en représentant dans un diagramme log‐log la
/
fonction en fonction de . En effet, de la relation admise ci‐
dessus :
log
5
y = 1,669x + 3,144
4 R² = 1
3
2
1
0 log
‐0,8 ‐0,3 0,2 0,7

Figure 30 : Détermination de l'exposant hydraulique n de la méthode de Bakhmeteff. Ici 1.7, soit n=3.4

56
 On peut écrire que
1
log
log log
2 2
Ainsi, la pente de la courbe obtenue (assimilable à une droite) n’est rien d’autre
que
• Tracer, en diagramme naturelle, la courbe . Cette courbe permet
de récupérer les valeurs de , pour en fonction des profondeurs donnée

0,7

0,6

0,5

0,4

0,3
0 1 2 3 4 5

Figure 31 : Courbe en fonction de y

La section de contrôle étant identifiée, on connait ,

1. on estime, à partir des 2 courbes, et donc .On en déduit en lisant sur le
table de Bakhmettef la valeur de .
2. On se donne une profondeur voisine de celle et on calcule la moyenne
arythmique de on fait la moyenne et de de pour obtenir . On estime par la suite
, puis on en déduire en lisant sur la table de Backhmeteff la valeur de
3. On en déduit la valeur de .
4. A partir de la section 1, on se donne une profondeur puis on recommence
[Link] Cas simplifié ou la variation de l’énergie cinétique est négligée
Si la variation de l’énergie cinétique est négligée, alors on aura

0 0
2 2

Il s’ensuit que l’équation différentielle de la ligne d’eau devient

1
1 1 1 1 1

On en déduit
57
 1 1 1 1
1
1 1 1
1
On peut en déduire la solution sous la forme
1
1
Ou encore sous la forme
1 1
1 1
On peut écrire finalement
, ,
En posant
,
La solution cherchée est alors :

Il existe des tables pour l’application de cette méthode, voir en annexe.

5.3.3 Méthode de Bresse

[Link] Principe de la méthode
Cette méthode s’applique dans le cas d’un canal infiniment large. Cette méthode utilise
la formule de Chézy pour exprimer les pentes, avec pour hypothèse que le coefficient de Chézy
est constant.

La section étant de largeur infiniment grande, le rayon hydraulique est confondu à la profondeur
normale . En effet, en prenant le cas de section trapézoïdale par exemple, on a

1
2 √1 1 2 √1

De la même manière, le canal étant infiniment large, la section mouillée peut être approximée comme
suit

On a donc

58
On en déduit une relation identique pour la condition normale (en tenant compte du fait que C est
supposé constant) :

On a alors

Par ailleurs, la condition critique peut s’écrire

1 0 1

Alors

1 1 1

1 1 1

Finalement, on a

1
1
1

En posant une fois encore, comme dans le cas de la méthode de Bakhmeteff :

On a

1
1

La solution de cette équation différentielle entre deux sections de profondeurs et et d’abscisses

et conduit à :

1 1
1 1 1

En posant comme dans le cas de la méthode de Bakhmeteff :

59
 , la solution de l’équation différentielle s’écrit alors

Il existe de table pour la fonction de Bresse, mais on peut déterminer cette fonction intégrale en utilisant
la fraction en une somme de fractions élémentaires :

1 1
1 1 1

Soit

1
1 1 1

Pour déterminer les coefficients a, b, c, on procède de la façon suivante :

• En multipliant les deux membre par 1 et en faisant 1, on déduit .

• En injectant 0 dans les deux membres et , on déduit que
• Enfin, pour 1, on obtient, avec les valeurs de a et de c :

Ce qui fait que la fraction peut s’écrire

1 1 1 2
1 3 1 1

Alors

1 1 1 2
1 3 1 1

C'est‐à‐dire que

1 1 1 1 2 1 3 1
1 3 1 2 1 2 1

Ou encore que

1 1 1 1
ln 1 ln 1
1 6 2 1

1 1 1 1 1
ln
1 6 1 2 1

Soit I, l’intégrale

60
 1 1 1 1 2 1
4 4 2. .
4 4 4 2 1 √3 √3 2 √3

On arrive finalement à

1 1 1 1 2 1
ln
1 6 1 √3 √3

La constante est déterminée en remarquant qu’on l’intégrale que nous venons de

calculer représente l’abaissement de la ligne d’eau sur une longueur . L’intégrale doit donc
être nulle pour le point de la courbe de remous situé sur le plan horizontal qui sert de référence
pour la mesure de cet abaissement. Or l’étude qualitative des courbes de remous (voir carlier) à
montré que ce cas de figure est obtenu lorsque ∞ en somme lorsque ∞.
En remplaçant cette remarque dans l’intégrale, on déduit la valeur de la constante C :

2√3
1
Au final, l’intégrale 3 1
donne
1 1 12 1 2 1
ln
3 1 6 2 1 √3 2 √3
Cette expression permet d’utiliser la table mise en place par Bresse pour estimer la
valeur de l’intégrale. Si on dispose d’une calculatrice, on peut utiliser directement l’expression
littérale de la fonction pour calculer la ligne d’eau.
Son expression de est la suivante :
ln .
√ √

[Link] Mise en œuvre de la méthode

• On calcule bien évidement les éléments particuliers , pour déduire le type de
remous, à partir de la connaissance de afin de donner le sens de progression des
calculs.
• On estime, à partir de , le paramètre à partir duquel on démarre les calculs. On
utilise, soit la table de Bresse, soit une calculatrice pour calculer .
• On se donne une valeur , voisine de celle de (plus petite ou plu grande en fonction
du type de remous) et on calcule . On peut en déduire la variation de longeur

• On se donne une autre valeur et on recommence le processus jusqu’à construire la

courbe de remous.
D’autres méthodes comme la méthode des différences finies, des approximations successives,
de Silber, la méthode graphique de Raytchine et Chatelain… permettent de calculer la ligne
d’eau. Elles sont présentées dans les ouvrages cités ci‐dessus

61
I. IV‐Ecoulement rapidement varié
1. Généralité
Les écoulements rapidement variés ou encoure écoulements brusquement variés sont des
écoulements permanents, qui se rencontrent en cas de changements brutaux de géométrie, par
exemple des convergents ou des divergents. En conséquence, les variables de l’écoulement
varient très rapidement dans l’espace. Pour ces écoulements, les lignes de courant deviennent
très courbes de sorte que l répartition des pressions cesse d’être hydrostatique, contrairement
à ce qui a été admis dans le cas des écoulements graduellement variés. Ces écoulement se
déroulent souvent sur de courtes distance, ce qui fait que l’effet des frottement peut être
négligé. L’étude de ce type d’écoulement utilise généralement le principe de la quantité de
mouvement. Dans ce document, l’étude de ces écoulements se résume à l’étude des réssaut
hydraulique.

2. Le ressaut hydraulique
2.1 Définition
Le ressaut hydraulique est une surélévation brusque du niveau d’eau lors du passage
d’un écoulement de type torrentiel à un écoulement de type fluvial. Ce passage se fait avec une
forte discontinuité du tirant d’eau, entrainant une forte agitation qui provoque une dissipation
d’une grande partie de l’énergie acquise dans la partie torrentielle du tronçon. On peut
observer un ou plusieurs rouleaux qui se produisent avec déferlement et turbulence. On
observe de grands tourbillons, des remous et de grands bulles d’air en mouvement.

Ligne de charge

2
2

Figure 32 : Définition du ressaut hydraulique

Pour ne pas être obligés de tenir compte des frottements internes au fluide entre deux sections
données, encadrant le ressaut, nous allons utiliser deux notions, conduisant au même résultat,
pour étudier les ressauts hydraulique : Il s’agit 1‐ du théorème de quantité de mouvement : la
variation de la quantité de mouvement entre les deux sections est égale à la somme des forces

62
appliquées (poids, le frottement sur les parois et les forces de pression) et 2‐ de la notion
d’impulsion totale. Mais avant, nous passons à la classification des ressauts hydrauliques.
2.2 Classification des ressauts hydrauliques
La classification des réssauts hydraulique que nous proposons ici est faite à l’aide du
nombre de Froude de la section amont du ressaut ( ).
√

Nbre de Froude Type de ressaut Forme du ressaut Caractéristique

Ondulé Seules quelques
légères rides sont
1 1.7 observées en surface

Faible Naissance en
1.7 2.5 surface de petits
tourbillons ou
rouleaux
Oscillant Production de
2.5 4.5 turbulences fortes
(surface et font) de
manière régulière
qui se propagent
loin à l’aval.
Etabli Ressaut bien
4.8 9 Localisé et efficace
en terme de
dissipation de
l’énergie
Fort Paquets d'eau
9 projetés par
intermitence
Tableau 12 : Typologie des ressauts Hydrauliques. Les figures proviennent de Mar, 2004

2.3 Théorème de quantité de mouvement (canal rectangulaire)

La variation de quantité de mouvement entre la section 1 et 2 s’écrit :
1 1

Par, ailleurs, en supposant que les deux sections sont assez proches, on peut négliger la
force de frottement sur les parois ainsi que le poids. Il reste alors la force de pression. Si le canal
est supposé uniforme, les forces exercées par les parois peuvent être considérées comme
perpendiculaire à l’axe ou s’annules tout simplement par symétrie. Il reste alors les forces
pressantes sur les sections d’entrée et de sortie et .
Or l’expression de la force pressante dans une section s’écrit

G étant le centre de gravité et S la section considérée.

63
 Finalement, le théorème de quantité de mouvement conduit à l’expression suivante :
1 1

Figure 33 : Application du théorème de quantités de mouvement

2.4 Notion d’impulsion totale

2.4.1 Définition
Le théorème de quantité de mouvement étudié ci‐dessus nous permet d’écrire
1 1

Soit

Cette égalité montre que l’application du théorème des quantités de mouvement conduit à la
conservation d’une grandeur M appelée impulsion totale qui se définit comme suit :

2.4.2 Etude des variations de l’impulsion totale

Nous allons étudier dans cette partie, les variations de la fonction impulsion totale.
Commençons par écrire sa dérivé

∆ ∆
∆

Figure 34 : Estimation du moment statique ∆

64
Le moment statique de S par rapport à la surface libre est
En cherchant le moment statique de ∆ par rapport à la surface libre, il faudra remarquer
∆
que le nouveau bras de levier de S est ∆ et que le bras de levier de ∆ est
Par ailleurs, la variation de surface due à la variation de y de ∆ s’écrit ∆ ∆
Au final, le moment statique de ∆ sera :

∆
∆ ∆ ∆
2

∆
∆ ∆ ∆
2
∆
∆ ∆ ∆
lim lim 2
∆ ∆ ∆ ∆
∆
∆ ∆ ∆
lim lim 2
∆ ∆ ∆ ∆
Par ailleurs, on avait montré, en étudiant les écoulements graduellement variés ci‐dessus, que

Au total, on a

La dérivée de l’impulsion totale s’annule en un point tel que

0 1

Au final,
0 1
Au final, l’impulsion totale s’annule au point critique tout comme la charge spécifique. Ainsi,
pour une impulsion totale donnée, on peut écouler un même débit à deux profondeurs
différentes et appelées profondeurs conjuguées. L’une, la plus petite ( ) correspond
à un écoulement en régime torrentiel et l’autre, la plus grande ( ), correspond à un
écoulement en régime fluvial.

Les variations de l’impulsion totale sont représentées sur la

Figure 33 ci‐dessous :

65
 M

Figure 35 : Variation de l’impulsion totale, pour un débit Q fixé, en fonction de la profondeur

2.5 Calcul de ressaut dans le cas d’un canal rectangulaire

2.5.1 Position du ressaut
Le ressaut est un phénomène qui dissipe énormément d’énergie. Il est donc important
de bien positionner car la zone dans laquelle ce phénomène a lieu peut être le siège d’une
importante activité d’érosion des berges et du fond. Le but de l’estimation de la position du
ressaut est donc de pouvoir mettre en place des ouvrages de protection des berges et du fond,
dans cette zone de dissipation d’énergie.
En supposant que la longueur du ressaut est nulle, la position du ressaut est
théoriquement le point où le tirant d’eau torrentiel en amont est égal au tirant d’eau fluvial en
aval. Dans la réalité, cette longueur n’est pas nulle et s’estime comme :
4.5 à 6

Fluvial

Conjugué

Torrentiel

Lr

66
 Figure 36 : Longueur du ressaut

2.5.2 Estimation graphique de la perte de charge du ressaut

Elle est synthétisée dans la figure ci‐dessous.

M

Figure 37 : Détermination graphique de la perte du ressaut hydraulique

67
2.5.3 Cas d’un canal rectangulaire
[Link] Calcul des profondeurs conjuguées
Si nous considérons un canal rectangulaire, alors on peut écrire la conservation de
l’impulsion totale :
1 1

Pour un canal rectangulaire,

;
2
Alors
. .
2 2 2
D’où
1 1 2
2
La solution est une solution triviale et n’est donc pas pris en compte, alors on a
2

En posant que

On a
2
0
Pour résoudre cette équation du second degré en , on calcule le discriminant :
2 8 8
∆ 4 1 1 1 8
La seule solution physiquement acceptable est
1 8
2
Soit finalement
1 1 8
2
De la même manière, on arrive à la solution
1 1 8
2
On calcul ainsi les tirants d’eau conjugués en fixant l’un comme variable et en le calculant en
prenant l’autre comme paramètre. Dans tous les cas, il suffit de fixer l’une des profondeurs
(amont ou aval) comme variable et de la calculer en fonction de la deuxième. La symétrie du
problème fait qu’on peut déduire la deuxième profondeur en fonction de la première, cela
dépendra de l’usage qu’on veu en faire.

68
[Link] Calcul de la perte de charge dans le Ressaut
Comme nous l’avions signalé plus haut, le théorème de Bernoulli n’est pas applicable dans le cas
d’un ressaut hydraulique. C’est d’ailleurs la raison pour laquelle nous avons fait recours au
théorème de quantité de mouvement pour estimer les tirants d’eau conjuguées. Or, il peut
s’avérer utile de connaitre la quantité d’énergie perdue dans un tel phénomène, raison pour
laquelle on va chercher à calculer la perte de charge (perte d’énergie par unité de poids)
générée dans ce phénomène. On mettra à profit les expressions définies ci‐dessus pour la
calculer. La perte de charge du ressaut par définition s’écrit :
∆
2 2
∆
2 2
1 1
∆
2 2 2
On sait que
2

on en déduit que
1 1
∆
4 4
Soit
3
∆
4 4 4 4 4
3 3 3 3
∆
4 4 4 4
Finalement
∆
4
[Link] Rendement du ressaut
Le rendement du ressaut hydraulique se définit comme étant le rapport du gain d’énergie
potentielle (énergie utile) par la variation observée de l’énergie cinétique(énergie fournie).

2 2
En développant, on a

1 1 1 1
2 4 4 4 4 4
Soit
69
 4

4 4
3 3 4

4 4
4 2

70
II. Bibliographie
Biaou, A, 2008, Logiciel de petits calculs en hydraulique générale, 2iE, Ouagadougou,
Burkina Faso
Carlier, M., 1972, Hydraulique générale et appliquée, Eyrolles, Paris, France
Lencastre, A., 2005, Hydraulique générale, Eyrolles, Paris, France
Mar, A. L., 2004, Cours d’hydraulique. Tome 2 : Ecoulement à surface libre, Polycopié de
cours de l’Institut International d’Ingénieur de l’Eau et de l’Environnement, Edité au 2iE,
Ouagadougou, Burkina Faso
Mar, A. L., 2004, Cours d’hydraulique. Tome 2 : Ecoulement en charge, Polycopié de
cours de l’Institut International d’Ingénieur de l’Eau et de l’Environnement, Edité au 2iE,
Ouagadougou, Burkina Faso

71
 Annexe A : Rappel de calculs géométriques élémentaires
Dans les problèmes d’hydraulique générale, surtout dans le cas de l’hydraulique à
surface libre, nous seront amenés à estimer des surfaces et des périmètres. En effet, la notion
de rayon hydraulique qui est le rapport de ces deux grandeurs, intervient assez constamment
dans les problèmes posés. Dans les cas simples, on connait les formules d’estimation d’un
périmètre, d’une surface, ou d’un volume. Mais lorsque les figures deviennent complexes, il
devient nécessaire de décomposer les figures en éléments simples dont on connait les formules.
Cet annexe est introduit pour montrer comment, à partir d’un infiniment petit, on établit
les formules des figures courantes et comment on peut décomposer une figure complexe en
plusieurs éléments simples facilement manipulable. Cette démarche permettra à l’étudiant de
s’inspirer des cas traités pour établir bien d’autres formules !

A1­ Rappels sur les Systèmes de coordonnées

Les formulations des problèmes rencontrés, en hydraulique et dans bien d’autres
disciplines, conduisent à l’établissement d’équations. Il est important de préciser le système de
coordonnées dans lequel ces équations ont été établies. Quatre systèmes d’équations sont
couramment utilisés :
¾ Le système de coordonnées cartésiennes
¾ Le système de coordonnées polaires
¾ Le système de coordonnées cylindriques
¾ Le système de coordonnées sphériques
Il est parfois nécessaire de connaitre le mode de passage d’un système donné à l’autre et
inversement. Dans ce qui suit, nous montrons sommairement comment passer de l’un des trois
derniers systèmes au système cartésien.
A1.1 Système de coordonnées cartésiennes
C’est de loin le système de coordonnées le plus utilisé dans les problèmes courants. Dans
ce système, les positions (x,y) d’un point sont représentées dans un repère orthonormé plan.
Une troisième composante (z) permet de représenter tout point de l’espace dans un repère
orthonormé à trois dimensions.

A1.2 Système de coordonnées polaires

Les coordonnées polaires sont pratiques dans les problèmes du le plan qui font
intervenir une rotation. Dans ce système, chaque point du plan est déterminé par une
coordonnée radiale (souvent notée r, et appelée rayon) qui exprime la distance du point à un
point central appelé pôle (équivalent à l’origine des coordonnées cartésiennes) et une
coordonnée angulaire (également appelée angle polaire ou azimut, et souvent notée θ), qui
exprime la mesure, dans le sens trigonométrique, de l’angle entre le point et la demi‐droite
d’angle 0°, appelé axe polaire (équivalent à l’axe des abscisses en coordonnées cartésiennes).

72
 Le passage du système de coordonnées polaires (r, ) au système de coordonnées
cartésiennes (x,y) se fait de la manière suivante :
cos
sin

y θ
x
Figure 38: coordonnées polaires

A1.3 Système de coordonnées cylindriques

En ajoutant une troisième dimension (z) au système de coordonnées polaires, on obtient
le système de coordonnées cylindriques. Cette troisième dimension permet de mesurer la
hauteur d’un point par rapport au plan repéré par les coordonnées polaires. Si les coordonnées
polaires permettent de repérer un point dans un plan, les coordonnées cylindriques permettent
de repérer tout point dans l’espace.
Le passage d’un système de coordonnées cylindriques au système de coordonnées
cartésiennes se fait de la manière suivante :

cos
sin

y
x θ

Figure 39: Cordonnées cylindriques

73
 A1.4 Système de coordonnées sphériques
C’est un système de coordonnées qui généralise le système (plan) de coordonnées polaires. Un point de
l’espace y est déterminé par une distance par rapport à l’origine et par deux angles, le premier par
rapport à un axe vertical et le deuxième par rapport à un axe horizontal (

Figure 38). Les coordonnées sphériques sont pratiques lorsque l’on a à faire à des
problèmes à symétrie sphérique comme c’est souvent le cas en géographie (latitude, longitude)
et en astrophysique.

[Link]θ
θ
x

ϕ
Figure 40:Cordonnées sphériques

Un point de coordonnées (x,y,z) dans un système cartésien est repéré dans un système
de coordonnées sphériques par les coordonnées , θ, , reliées aux coordonnées cartésiennes
par les équations suivantes :

cos θ
sin θ cos
sin sin

A2­ Elément de calcul de périmètre

A2.1 Périmètre de cercle
On considère un arc infinitésimal (voir Figure 39). Sa longueur est :

Cercle plein

La longueur total du cercle (son périmètre donc) s’obtient en sommant les arcs
infinitésimaux de ce type (intégration) sur tour le pourtour du cercle. Puisque seul l’angle

74
azimutal varie, il suffira donc d’intégrer la longueur élémentaire pour θ allant de 0 à 2 (faire
un tour complet):

2 .

θ dl=R dθ
R

dθ
x

Figure 41: Calcul du périmètre d'un cercle

Cercle non fermé

2R

y
θ

Figure 42: Calcul du périmètre d'un cercle non fermé

Dans le cas d’un cercle non fermé (ce cas sera rencontré dans les calculs de périmètre
mouillé de conduite à section circulaire), délimité par l’angle , il suffira d’intégrer la longueur
élémentaire pour l’angle allant de allant 0 à θ :

A2.2 Périmètre d’une parabole

Quand bien même on ne rencontre pas très souvent des canalisations à section
paraboloïdale, nous montrons comment exprimer le périmètre d’une figure complexe. L’objectif
étant simplement de montrer comment procéder..
L’équation de la parabole est du type .

75
 Pour calculer le périmètre de la partie de la parabole définie par L et y, on considère
l’élément de longueur dP tel que

Soit

1 1

Figure 43:Calcul du périmètre d'une paraboloïde

Comme le problème est symétrique, on a

2 1 2 2 2 1 dx

Pour k=2,

2 1 4

Faisons le changement de variable

2 2

Alors l’intégrale exprimant le périmètre cherché P peut s’écrire encore sous la

forme (sans les bornes):
1
2 1 4 1

On montre que

76
 1
1 2 1 2 ln 1 2
2 2
Alors
1
2 1 4 1 2 ln 1 2
2
1
√1 4 ln 2 √1 4
2

En introduisant les bornes 0 et de l’intégral de P, on obtient

1
4 1 4 ln 4 1 4
2

A2.3 Périmètre d’un rectangle arrondi sur les bords de sa base

r
r
b
Figure 44: Calcul du périmètre d'un rectangle dont les bords sont arrondis à la base.

A présent, nous nous proposons de déterminer le périmètre du rectangle arrondi ci‐

dessous ( , sur sa partie en trait épais (Figure 42). Cette détermination a pour but de
montrer comment on peut décomposer un problème compliqué en des problèmes simples.
Ainsi, pour déterminer le périmètre d’une telle figure, une décomposition s’impose.

Le rectangle peut être décomposé en 2 rectangles de dimensions respectives (y‐r,b+r) et

(r,b) et deux quarts de cercle identiques dont le rayon commun est r.
Pour le premier rectangle, juste deux côtés nous intéressent, puisque les parties
supérieure et inferieure ne font pas partie du périmètre total. Donc la partie du périmètre qui
nous intéresse dans ce rectangle reste les deux côtés de longueur y‐r, soit 2 .
Le deuxième rectangle ne contribue au périmètre initial que par sa partie inférieure,
c'est‐à‐dire b. La longueur d’un quart de cercle peut se calculer comme ci‐dessus :

Au total, le périmètre cherché est donc

77
 A2.3 Périmètre d’un triangle à fond arrondi
L’approche de résolution de ce type de problème est semblable à celle utilisée dans le cas ci‐
dessus. Le fond du triangle arrondi est un demi‐cercle de rayon .

Le triangle initial (A’,B’, C’) est tangent au cercle (de centre A et de rayon r) en ses points de contact B et
C avec celui‐ci. En conséquence, l’angle au sommet (A) de la partie circulaire est égal à 2 . Par ailleurs,

alors ou encore .

Or,

L
A’C’ A’B’ cos
2

Il s’en suit immédiatement que

L
CC’ BB’ cos . 1
2 2

L
B’ C’

y’
1
ϕ
y
A
r
D
B O ϕ C

A’

Figure 45 : triangle à fond arrondi

Le périmètre cherché est donc égal à

2 1 2
2
2 2
1 1 1

78
 A3 Elément de calcul de surface
A3.1 Calcul de l’aire de rectangle
On considère, dans un système de coordonnées cartésiennes, une surface infiniment petite

. . .

l
d2s
dx
dy
L x
Figure 46: Calcul de la surface d'un rectangle

A3.2 Calcul de l’aire d’un triangle connaissant la base et la hauteur.

b

S1
S2 h

b1 b2
Figure 47:calcul de l'aire d'un triangle

L’aire S du triangle rouge peut être calculée à partir de la connaissance de l’aire de

chacun des rectangles S1 et S2 (on sait déjà calculer l’aire d’un rectangle). Le triangle dont nous
voulons estimer la surface peut être décomposé en deux triangles droits (le côté commun de
ces deux triangles est la hauteur du triangle initial) dont les aires sont les moitiées des
rectangles S1 et S2. Ainsi, l’aire recherchée se calcul comme suit :

2 2 2 2 2 2
Ce calcul confirme ce qui a été enseigné depuis le cours primaire :

Aire trianle
2

79
 A3.3 Calcul de surface d’un disque
A3.2 1 Cas d’un disque plein
On se place cette fois en coordonnées polaire, puis on considère de ce système, un élément infinitésimal
dont on calcul la surface :

1
. .2 .
2

dr r dθ
θ d2S
dθ x

Figure 48: Calcul de la surface d'un cercle

A3.2 2 Cas d’un disque en partie plein

La partie pleine du disque est délimitée par la ligne épaisse. Ce type de problème se
rencontre lorsqu’on cherche les sections mouillées des conduites circulaires.
La surface cherchée est la somme de la partie inferieure non hachurée et de la partie
triangulaire hachurée.
Pour la partie non hachurée on considère l’élément de surface en coordonnées polaire comme le cas du
disque plein :
.

En intégrant, on obtient

80
 1
.
2

2R
y
θ

Figure 49: Calcul de l'air d'un disque en partie plein

Pour la partie triangulaire, nous appliquons la formule établie ci‐dessus :

S
2
2
base 2 2
2 2
2
hauteur
2 2

Alors

1
S
2 2 2

La surface cherchée s’écrit donc

1 2
2

A3.4 Calcul de surface d’une sphère

Nous nous plaçons cette fois‐ci dans un système de coordonnées sphériques. Comme
dans les cas précédents, considérons un élément de surface infinitésimal
. ϕ

. . 0 .2 4

81
 z
R dθ
Rsinθ.dϕ
θ
ds
y

ϕ
x

Figure 50: Calcul de la surface d'une sphère

Profitons au passage pour calculer le volume d’une sphère

. . ϕ
Le volume total est obtenu en prenant l’intégrale triple du volume élémentaire ci‐dessous

1 4
. . . 0 .2
3 3

z
r dθ
dr rsinθ.dϕ
θ
y

x ϕ

Figure 51:Calcul de volume d'une sphère

82
 A3.4 Aire d’une parabole
/ /

Figure 52:Calcul de l’aire d'une paraboloïde

Rappelons que l’équation de la parable (cas n=2 traité ici) s’écrit :

La figure présentant une symétrie par rapport à l’axe horizontal, il suffit donc de calculer
l’aire pour les x positif et de la multiplier par 2 pour obtenir l’aire totale :

y 2
x y dy dy
a 3

4 2
2 1 , 2
3 3

A3.5 Aire d’un rectangle arrondi sur les bords de la base

r r

b
Figure 53: Calcul de l’aire d'un rectangle dont les bords sont arrondis à la base.

En se basant sur la décomposition faite dans la partie de calcul du périmètre, on a

2 . 2 2 2

83
 Annexe B : Table de la fonction de Bakhmeteff

n=2.8 n=3.0 n=3.2 n=3.4 n=3.6 n=3.8 n=4.0 n=4.2 n=4.6 n=5.0 n=5.4
0.00 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.02 0.020 0.020 0.020 0.020 0.020 0.020 0.020 0.020 0.020 0.020 0.020
0.04 0.040 0.040 0.040 0.040 0.040 0.040 0.040 0.040 0.040 0.040 0.040
0.06 0.060 0.060 0.060 0.060 0.060 0.060 0.060 0.060 0.060 0.060 0.060
0.08 0.080 0.080 0.080 0.080 0.080 0.080 0.080 0.080 0.080 0.080 0.080
0.1 0.100 0.100 0.100 0.100 0.100 0.100 0.100 0.100 0.100 0.100 0.100
0.12 0.120 0.120 0.120 0.120 0.120 0.120 0.120 0.120 0.120 0.120 0.120
0.14 0.140 0.140 0.140 0.140 0.140 0.140 0.140 0.140 0.140 0.140 0.140
0.16 0.160 0.160 0.160 0.160 0.160 0.160 0.160 0.160 0.160 0.160 0.160
0.18 0.180 0.180 0.180 0.180 0.180 0.180 0.180 0.180 0.180 0.180 0.180
0.20 0.201 0.200 0.200 0.200 0.200 0.200 0.200 0.200 0.200 0.200 0.200
0.22 0.221 0.221 0.220 0.220 0.220 0.220 0.220 0.220 0.220 0.220 0.220
0.24 0.241 0.241 0.241 0.240 0.240 0.240 0.240 0.240 0.240 0.240 0.240
0.26 0.262 0.261 0.261 0.261 0.260 0.260 0.260 0.260 0.260 0.260 0.260
0.28 0.282 0.282 0.281 0.281 0.281 0.280 0.280 0.280 0.280 0.280 0.280
0.30 0.303 0.302 0.302 0.301 0.301 0.301 0.300 0.300 0.300 0.300 0.300
0.32 0.324 0.323 0.322 0.322 0.321 0.321 0.321 0.321 0.320 0.320 0.320
0.34 0.344 0.343 0.343 0.342 0.342 0.341 0.341 0.341 0.340 0.340 0.340
0.36 0.366 0.364 0.363 0.363 0.362 0.362 0.361 0.361 0.361 0.360 0.360
0.38 0.387 0.385 0.384 0.383 0.383 0.382 0.382 0.381 0.381 0.381 0.380
0.40 0.408 0.407 0.405 0.404 0.403 0.403 0.402 0.402 0.401 0.401 0.400
0.42 0.430 0.428 0.426 0.425 0.424 0.423 0.423 0.422 0.421 0.421 0.421
0.44 0.452 0.450 0.448 0.446 0.445 0.444 0.443 0.443 0.442 0.441 0.441
0.46 0.475 0.472 0.470 0.468 0.466 0.465 0.464 0.463 0.462 0.462 0.461
0.48 0.497 0.494 0.492 0.489 0.488 0.486 0.485 0.484 0.483 0.482 0.481
0.50 0.521 0.518 0.514 0.511 0.509 0.508 0.506 0.505 0.504 0.503 0.502
0.52 0.544 0.540 0.536 0.534 0.531 0.529 0.528 0.527 0.525 0.523 0.522
0.54 0.568 0.563 0.559 0.556 0.554 0.551 0.550 0.548 0.546 0.544 0.543
0.56 0.593 0.587 0.583 0.579 0.576 0.574 0.572 0.570 0.567 0.565 0.564
0.58 0.618 0.612 0.607 0.603 0.599 0.596 0.594 0.592 0.589 0.587 0.585
0.60 0.644 0.637 0.631 0.627 0.623 0.620 0.617 0.614 0.611 0.608 0.606
0.61 0.657 0.650 0.644 0.639 0.635 0.631 0.628 0.626 0.622 0.619 0.617
0.62 0.671 0.663 0.6567 0.651 0.647 0.643 0.640 0.637 0.633 0.630 0.628
0.63 0.684 0.676 0.669 0.664 0.659 0.655 0.652 0.649 0.644 0.641 0.638
0.64 0.689 0.690 0.683 0.677 0.672 0.667 0.664 0.661 0.656 0.652 0.649
0.65 0.712 0.703 0.696 0.689 0.684 0.680 0.676 0.673 0.667 0.663 0.660
0.66 0.727 0.617 0.709 0.703 0.697 0.692 0.688 0.685 0.679 0.675 0.672
0.67 0.742 0.731 0.723 0.716 0.710 0.705 0.701 0.697 0.691 0.686 0.683

84
 Annexe B : Table de la fonction de Bakhmeteff

n=2.8 n=3.0 n=3.2 n=3.4 n=3.6 n=3.8 n=4.0 n=4.2 n=4.6 n=5.0 n=5.4
0.68 0.757 0.746 0.737 0.729 0.723 0.718 0.713 0.709 0.703 0.698 0.694
0.69 0.772 0.761 0.751 0.743 0.737 0.731 0.726 0.722 0.715 0.720 0.706
0.70 0.787 0.776 0.766 0.757 0.750 0.744 0.739 0.735 0.727 0.722 0.717
0.71 0.804 0.791 0.781 0.772 0.764 0.758 0.752 0.748 0.740 0.734 0.729
0.72 0.820 0.807 0.796 0.786 0.779 0.772 0.766 0.761 0.752 0.746 0.741
0.73 0.837 0.823 0.811 0.802 0.793 0.786 0.780 0.774 0.765 0.759 0.753
0.74 0.854 0.840 0.827 0.817 0.808 0.800 0.794 0.788 0.779 0.771 0.766
0.75 0.872 0.857 0.844 0.833 0.823 0.815 0.808 0.802 0.792 0.784 0.778
0.76 0.890 0.874 0.861 0.849 0.839 0.830 0.823 0.817 0.792 0.784 0.778
0.77 0.909 0.892 0.878 0.866 0.855 0.846 0.838 0.831 0.820 0.811 0.804
0.78 0.929 0.911 0.896 0.883 0.872 0.862 0.854 0.847 0.834 0.825 0.817
0.79 0.949 0.930 0.914 0.901 0.889 0.879 0.870 0.862 0.849 0.839 0.831
0.80 0.970 0.950 0.934 0.919 0.907 0.896 0.887 0.878 0.865 0.854 0.845
0.81 0.992 0.971 0.954 0.938 0.925 0.914 0.904 0.895 0.881 0.869 0.860
0.82 1.015 0.993 0.974 0.958 0.945 0.932 0.922 0.913 0.897 0.885 0.875
0.83 1.039 1.016 0.996 0.979 0.965 0.952 0.940 0.931 0.914 0.901 0.890
0.84 1.064 1.040 1.019 1.001 0.985 0.972 0.960 0.949 0.932 0.918 0.906
0.85 1.091 1.065 1.043 1.024 1.007 0.993 0.980 0.969 0.950 0.935 0.923
0.86 1.119 1.092 1.068 1.048 1.031 1.015 1.002 0.990 0.970 0.954 0.940
0.87 1.149 1.120 1.095 1.074 1.055 1.039 1.025 1.012 0.990 0.973 0.959
0.88 1.181 1.151 1.124 1.101 1.081 1.064 1.049 1.035 1.012 0.994 0.978
0.89 1.216 1.183 1.155 1.131 1.110 1.091 1.075 1.060 1.035 1.015 0.999
0.90 1.253 1.218 1.189 1.163 1.140 1.120 1.103 1.087 1.060 1.039 1.021
0.91 1.294 1.257 1.225 1.197 1.173 1.152 1.133 1.116 1.088 1.064 1.045
0.92 1.340 1.300 1.226 1.236 1.210 1.187 1.166 1.148 1.117 1.092 1.072
0.93 1.391 1.348 1.311 1.279 1.251 1.226 1.204 1.184 1.151 1.123 1.101
0.94 1.449 1.403 1.363 1.328 1.297 1.270 1.246 1.225 1.188 1.158 1.134
0.95 1.518 1.467 1.423 1.385 1.352 1.322 1.296 1.272 1.232 1.199 1.172
0.96 1.601 1.545 1.497 1.454 1.417 1.385 1.355 1.329 1.285 1.248 1.217
0.97 1.707 1.644 1.590 1.543 1.501 1.464 1.431 1.402 1.351 1.310 1.275
0.975 1.773 1.707 1.649 1.598 1.554 1.514 1.479 1.447 1.393 1.348 1.311
0.980 1.855 1.783 1.720 1.666 1.617 1.575 1.536 1.502 1.443 1.395 1.354
0.985 1.959 1.880 1.812 1.752 1.699 1.652 1.610 1.573 1.508 1.454 1.409
0.990 2.106 2.017 1.940 1.873 1.814 1.761 1.714 1.671 1.598 1.537 1.487
0.995 2.355 2.250 2.159 2.079 2.008 1.945 1.889 1.838 1.751 1.678 1.617
0.999 2.931 2.788 2.663 2.554 2.457 2.370 2.293 2.223 2.102 2.002 1.917
1.001 2.399 2.184 2.008 1.856 1.725 1.610 1.508 1.417 1.264 1.138 1.033
1.005 1.818 1.649 1.506 1.384 1.279 1.188 1.107 1.036 0.915 0.817 0.737
1.005 1.818 1.649 1.506 1.384 1.279 1.188 1.107 1.036 0.915 0.817 0.737

85
 Annexe B : Table de la fonction de Bakhmeteff

n=2.8 n=3.0 n=3.2 n=3.4 n=3.6 n=3.8 n=4.0 n=4.2 n=4.6 n=5.0 n=5.4
1.010 1.572 1.419 1.291 1.182 1.089 1.007 0.936 0.873 0.766 0.681 0.610
1.015 1.428 1.286 1.166 1.065 0.978 0.902 0.836 0.778 0.680 0.602 0.537
1.02 1.327 1.191 1.078 0.982 0.900 0.828 0.766 0.711 0.620 0.546 0.486
1.03 1.186 1.060 0.955 0.866 0.790 0.725 0.668 0.618 0.535 0.469 0.415
1.04 1.086 0.967 0.868 0.785 0.714 0.653 0.600 0.554 0.477 0.415 0.365
1.05 1.010 0.896 0.802 0.723 0.656 0.598 0.548 0.504 0.432 0.374 0.328
1.06 0.948 0.838 0.748 0.672 0.608 0.553 0.506 0.464 0.396 0.342 0.298
1.07 0.896 0.790 0.703 0.630 0.569 0.516 0.471 0.431 0.366 0.315 0.273
1.08 0.851 0.749 0.665 0.595 0.535 0.485 0.441 0.403 0.341 0.292 0.252
1.09 0.812 0.713 0.631 0.563 0.506 0.457 0.415 0.379 0.319 0.272 0.234
1.10 0.777 0.681 0.601 0.536 0.480 0.433 0.392 0.357 0.299 0.254 0.218
1.11 0.746 0.652 0.575 0.511 0.457 0.411 0.372 0.338 0.282 0.239 0.204
1.12 0.718 0.626 0.551 0.488 0.436 0.392 0.354 0.321 0.267 0.225 0.192
1.13 0.692 0.602 0.529 0.468 0.417 0.374 0.337 0.305 0.253 0.212 0.181
1.14 0.669 0.581 0.509 0.450 0.400 0.358 0.322 0.291 0.240 0.201 0.170
1.15 0.647 0.561 0.490 0.432 0.384 0.343 0.308 0.278 0.229 0.191 0.161
1.16 0.627 0.542 0.473 0.417 0.369 0.329 0.295 0.266 0.218 0.181 0.153
1.17 0.608 0.525 0.458 0.402 0.356 0.317 0.283 0.255 0.208 0.173 0.145
1.18 0.591 0.509 0.443 0.338 0.343 0.305 0.272 0.244 0.199 0.165 0.138
1.19 0.594 O.494 0.429 0.375 0.331 0.294 0.262 0.235 0.191 0.157 0.131
1.2 0.559 0.480 0.416 0.363 0.320 0.283 0.252 0.226 0.183 0.150 0.125
1.22 0.531 0.454 0.392 0.341 0.299 0.264 0.235 0.209 0.168 0.138 0.114
1.24 0.505 0.431 0.371 0.322 0.281 0.248 0.219 0.195 0.156 0.127 0.104
1.26 0.482 0.410 0.351 0.304 0.265 0.233 0.205 0.182 0.145 0.117 0.095
1.28 0.461 0.391 0.334 0.288 0.250 0.219 0.193 0.170 0.135 0.108 0.088
1.30 0.442 0.373 0.318 0.274 0.237 0.207 0.181 0.160 0.126 0.100 0.081
1.32 0.424 0.357 0.304 0.260 0.225 0.196 0.171 0.150 0.118 0.093 0.075
1.34 0.408 0.342 0.290 0.248 0.214 0.185 0.162 0.142 0.110 0.087 0.069
1.36 0.393 0.329 0.278 0.237 0.204 0.176 0.153 0.134 0.103 0.081 0.064
1.38 0.378 0.316 0.266 0.226 0.194 0.167 0.145 0 .127 0.097 0.076 0.060
1.40 0.365 0.304 0.256 0.217 0.185 0.159 0.138 0.120 0.092 0.071 0.056
1.42 0.353 0.293 0.246 0.208 0.177 0.152 0.131 0.114 0.087 0.067 0.052
1.44 0.341 0.282 0.236 0.199 0.169 0.145 0.125 0.108 0.082 0.063 0.049
1.46 0.330 0.273 0.227 0.191 0.162 0.139 0.119 0.103 0.077 0.059 0.046
1.48 0.320 0.263 0.219 0.184 0.156 0.133 0.113 0.098 0.073 0.056 0.043
1.50 0.310 0.255 0.211 0.177 0.149 0.127 0.108 0.093 0.069 0.053 0.040
1.55 0.288 0.235 0.194 0.161 0.135 0.114 0.097 0.083 0.061 0.046 0.035
1.60 0.269 0.218 0.179 0.148 0.123 0.103 0.087 0.074 0.054 0.040 0.030

86
 Annexe B : Table de la fonction de Bakhmeteff

n=2.8 n=3.0 n=3.2 n=3.4 n=3.6 n=3.8 n=4.0 n=4.2 n=4.6 n=5.0 n=5.4
1.60 0.269 0.218 0.179 0.148 0.123 0.103 0.087 0.074 0.054 0.040 0.030
1.65 0.251 0.203 0.165 0.136 0.113 0.094 0.079 0.067 0.048 0.035 0.026
1.70 0.236 0.189 0.153 0.125 0.103 0.086 0.072 0.060 0.043 0.031 0.023
1.75 0.212 0.177 0.143 0.116 0.095 0.079 0.065 0.054 0.038 0.027 0.020
1.80 0.209 0.166 0.133 0.108 0.088 0.072 0.060 0.049 0.034 0.024 0.017
1.85 0.198 0.156 0.125 0.100 0.082 0.067 0.055 0.045 0.031 0.022 0.015
1.90 0.188 0.156 0.117 0.094 0.076 0.062 0.050 0.041 0.028 0.020 0.014
1.95 0.178 0.139 0.110 0.088 0.070 0.057 0.046 0.038 0.026 0.018 0.012
2.00 0.169 0.132 0.104 0.082 0.066 0.053 0.043 0.035 0.023 0.016 0.011
2.1 0.154 0.119 0.092 0.073 0.058 0.046 0.037 0.030 0.019 0.013 0.009
2.2 0.141 0.107 0.083 0.065 0.051 0.040 0.032 0.025 0.016 0.011 0.007
2.3 0.129 0.098 0.075 0.058 0.045 0.035 0.028 0.022 0.014 0.009 0.006
2.4 0.119 0.089 0.068 0.052 0.040 0.031 0.024 0.019 0.012 0.008 0.005
2.5 0.110 0.082 0.062 0.047 0.036 0.028 0.022 0.017 0.010 0.006 0.004
2.6 0.102 0.076 0.057 0.043 0.033 0.025 0.019 0.015 0.009 0.005 0.003
2.7 0.095 0.070 0.052 0.039 0.029 0.022 0.017 0.013 0.008 0.005 0.003
2.8 0.089 0.065 0.048 0.036 0.027 0.020 0.015 0.012 0.007 0.004 0.002
2.9 0.083 0.060 0.044 0.033 0.024 0.018 0.014 0.010 0.006 0.004 0.002
3.0 0.078 0.056 0.041 0.030 0.022 0.017 0.012 0.009 0.005 0.003 0.002
3.5 0.059 0.041 0.029 0.021 0.015 0.011 0.008 0.006 0.003 0.002 0.001
4.0 0.046 0.031 0.022 0.015 0.010 0.007 0.005 0.004 0.002 0.001 0.000
4.5 0.037 0.025 0.017 0.011 0.008 0.005 0.004 0.003 0.001 0.001 0.000
5.0 0.031 0.020 0.013 0.009 0.006 0.004 0.003 0.002 0.001 0.000 0.000
6.0 0.022 0.014 0.009 0.006 0.004 0.002 0.002 0.001 0.000 0.000 0.000
7.0 0.017 0.010 0.006 0.004 0.002 0.002 0.001 0.001
8.0 0.013 0.008 0.005 0.003 0.002 0.001 0.001 0.000
9.0 0.011 0.006 0.004 0.002 0.001 0.001 0.000 0.000
10.0 0.009 0.005 0.003 0.002 0.001 0.001 0.000 0.000
20.0 0.003 0.002 0.001 0.001 0.000 0.000 0.000 0.000

87