formule de Joukowsky (sur- & dpression) : h =

au
p
aQ
= 0 =
g
g
gS

Chapitre 7

ETUDE DU
COUP DE BELIER [ Marteau dEau ]
[ Ondes lastiques ]
Si lun quelconque des appareils en service subit une variation dans le
temps : ouverture ou fermeture de vannes, ralentissement ou acclration
de la vitesse angulaire dune pompe, changement douverture dune
turbine, dtente dun gaz comprim, variation du plan leau amont ou aval
dun rservoir, etc., la pression p et le dbit Q varient aussi dans le
temps en tous points du systme [on dit : rgime hydraulique vari ] :
exemple du rservoir niveau constant
Ps = pression la sortie
h0
A
R
Caractristique

OA = h 0
te

Ps = h 0 C Q 2
O

La caractristique est obtenue en retranchant les pertes de charge de la

dnivele h 0 .
HISTORIQUE :
Historiquement au moyen ge pour enfoncer les portails des forteresses,
les romains et les vikings utilisaient des troncs darbre dont lavant avait la
forme de tte de blier, par ce moyen les envahisseurs parvenaient ouvrir
en enfonant et en dtruisant celles ci par des coups rpts :
disons quils donnaient des coups de blier aux portails.
Par analogie les surpressions dpressions qui se produisent dans un
coulement en charge par une variation quelconque du dbit et la
propagation de londe qui en rsulte est btit : COUP DE BELIER
(Pour plus dinformation sur la formulation mathmatique < par la mthode
des caractristiques > de ce phnomne ainsi que les diffrents cas
pratiques consultez louvrage de Louis. BERGERON :
DU COUP DE BLIER EN HYDRAULIQUE
AU COUP DE FOUDRE EN ENECTRICITE
Mthode graphique gnrale
DUNOD Paris 1950 pp 336). Dans la suite on prsente une formulation
simplifie pour comprendre ce phnomne et ses mcanismes.
E.H.T.P.
Zorkani Mohammed
7-1
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

1) Gnralit :
Considrons une conduite, en charge, alimente par un rservoir deau, le
rgime permanent tant tabli : voir figure suivante

v0
Rservoirr

x0

L
Si a un instant t, choisi gal zro comme origine des temps, on diminue la
section de passage de la vanne A une surpression prend immdiatement
naissance proximit de la vanne A car le fluide avant larrt total continu
couler en saccumulant proche de A.
Cette surpression de A va se propager vers B, o elle a beaucoup de
chance de se rflchir en une dpression se propageant vers A, et ainsi de
suite : trains dondes qui samortissent la longue.
Limportance de la surpression sera fonction de la manuvre de la vanne
A, et du rgime prcdemment tabli dans la conduite.
2) Mise en quation : On dfinit le coefficient de compressibilit par
1 V
1
1 V
1
-1
-1

=
et S enthropie
=
T temprature
V p T p T
V p S p S
2 1 ) Hypothses :
A) La conduite est rectiligne uniforme, et de plus, mince et
parfaitement lastique, ceci mous permettra de relier variation de
section et variation de pression.
B) Le fluide, bien que peu compressible, a un coefficient de
dp
compression constant [ T = ]. Ce module, pour les
d T

liquides, est sensiblement gal au coefficient de compressibilit

adiabatique S = T = et sa valeur peu fonction de la temprature.
C) Le fluide est suppos non visqueux. Les vitesses u , les pressions

motrices p = p + g z sont constantes dans une section droite

[dbitante].
E.H.T.P.
Zorkani Mohammed
7-2
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

D) Nous ne tiendrons pas compte des efforts longitudinaux sur la

conduite qui seront absorbs par les dispositifs dancrage.
D
E) Pour calculer les drives particulaires
nous admettrons que
Dt
laspect ondulatoire du phnomne nous permettra dcrire :
D

=
: c d ladvection ngligeable (thorie linaire)
Dt t
Cette dernire hypothse devra tre vrifier la fin du calcul et sur
des exemples.
2 2 ) Mise en quation :
Nous allons immdiatement traduire ces hypothses sous forme dquation
a) Elasticit de la conduite : daprs les relations modifies de la
rsistance des matriaux : La contrainte n e vaut :
[ Voir : Notes sur Influence ETC. p42 pour la justification ]
formule des chaudires
D
D
p D = 2 e n e n e = contrainte = p dn e
dp
2e
2e
D
En introduisant le module dYOUNG E de la conduite : dD dn e
E
[ Caractrise llasticit des parois de la conduite ]
dD
dS
D2
D
car S =
alors
Sachant que :
=2
D
S
4
dS 1 D
=
dp (1)
e
D
S
E
e
ne
10
2
E acier 2 10 Kg/m
dD
b) Equation dtat du fluide :
1 1
d 1

= dp (2 ) : eau 2,07 108 Kg/m 2

p T ou S

c) Equation de continuit :

x + dx

u (x )

u (x + dx )

La quantit de matire qui entre en x au temps t est ( Su )x et celle qui

sort en x + dx au temps t est ( Su )x + dx . Pendant le temps dt la
quantit de matire qui saccumule vaut : d 2 m = [( Su )x ( Su )x + dx ] dt

E.H.T.P.
Zorkani Mohammed
7-3
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Dautre part, linstant t, la masse contenue entre x et x + dx vaut

(dm )t = (Sdx )t . A linstant t + dt cette masse vaut (dm )t + dt = (Sdx )t + dt
la diffrence d 2 m scrit alors

d 2 m = [( Sdx )t + dt ( Sdx )t ]

( S u ) = ( S) (3)
Soit

x
t
d) Equation du mouvement : dynamique
1 p
Du
(4 ) Euler
=
x
Dt
D

c d on
=
Dt t
nglige ladvection (linarisation) alors les quations [(1) (4)] deviennent :
dS 1 D
(1)
= dp
S E e
d 1
(2 )
= dp

En confondant la drive particulaire avec celle local

( S) = S u
t
x

(3)

1 p
u
(4 )
=
t
x
Dautre part, on admet que :
S / variation de S, faible devant S0
S = S0 + S / o
S0 section de passage en rgime normal
De mme pour linariser les quations on pose :
= 0 + /

/ , p / , u / Les effets du coup de blier

/
=
+
avec
p
p
p

0
{ 0 , p 0 , u 0 } Les valeurs en rgime normal

/
=
+
u
u
u

0
La prime / signifie les variations par rapport aux valeurs normales.
Ces 4 relations deviennent :

E.H.T.P.
Zorkani Mohammed
7-4
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

dS / 1 D 0 /
=
dp
S0 E e

(1)

d /

(2 )

dp /

/
u/
/
/ /
S 0 + 0S + S = 0S0
t
x

(3)

1 p/
u/
=
0 x
t

(4 )

En liminant S/ et / de ces 4 quations, nous obtenons :

D0
u/
1 p/

+
=
x
Ee t

(a1)

1 p/
u/
(b1)

=
0 x
t
En prenant pour origine des x, le distributeur (point A) (voir figure), ces 2
quations deviennent :
D0
1 p/ u /

(a2)
+
=
x
Ee t
1 p/ u /
=
0 x
t

(b2 )

En liminant u / ou p / de lune de ces 2 quations, et en reportant dans

lautre, et en posant :
p / = 0 gh

dsigne la clrit des ondes planes (Allievi 1904)

1
a:
D
se propageant dans le milieu fini
1+
e E
c = 0 : dsigne la clrit des ondes planes se propageant dans le
c
tend vers
mme milieu mais infini : D + ainsi : a =
D
1+
e E
Ces 2 relations deviennent :
a=

E.H.T.P.
Zorkani Mohammed
7-5
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Postulat
1 h 2h
(A )

=
a2 t2 x2
2

1 2u / 2u /
( )

=
a2 t2
x2

une variation d' tat ou de rgime provoque en un lieu d' une

conduite , se propage de part et d' autre sous forme d' onde plane

Equations dONDES de pression (prorogatives)

x
x
te

F t + f t + = h + C
a
a

Lintgration de ces 2 relations donne :

x
x
te

G t + g t + = u / + C
a
a

Compte tenu des relations diffrentielles existantes entre p / (ou h ) et u / on

x
x

Aval
onde d' amont
h = F t + f t +
f
a
a

peut crire:
g
x
x
x
u = u 0 F t f t +
F
onde
d'
aval
a
a
a
variation de pression = h = F + f
on a
c - - d : En x et l' instant t donns
a a
(
)
variation
de
vitesse

= u u =Ff

g
g 0

Ces 2 relations sont dites dAllievi : quations fondamentales du

rgime vari dans les conduites. [L. Allievi (1895 1942) : Toria generale del

moto pertubado dellacqua nei tubi in pressione , Milan 1903 , Teoria del colpo dariete,
Milan 1913 { traduit en franais par Allievi Gaden : Thorie du coup de blier, Dunod,
Paris 1921} ].

Remarque : leau 150 C : a = 1425m / s [ 2,07 10 4 [Link] -2 2,15 10 9 m -2 ]

Pour une conduite mince en acier :
D
= 100
E 2,15 10 6 Kgf/cm 2 Et e
On obtient pour a une valeur : a = 1017m/s

3
3
= 10 Kg/m
Remarque : Il faut maintenant vrifier hypothse de linarisation :
u g /

= 2 F +f/
x a
u u
g

u
et comme u est de la forme [F f ] alors

g /
u
x t
a
/
= F f

a
t

E.H.T.P.
Zorkani Mohammed
7-6
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

g u
u u
u
+
= 1 F / 1 + f /
a a
a
x t

u
sont trs petits [u tant de lordre de quelque m/s au
a
u
u u
g
u
F/ f / =
+
maximum] : 1 . Dans ce cas : u
a
a
t
x t
Lhypothse (e) est bien vrifie.
dS
Remarque : Si la conduite nest pas mince, paroi paisse, la loi
= f (dp )
S
nest plus celle nonce. Dans ce cas, on obtient pour a, lexpression :
D : diamtre intrieur
c
a=
avec 1
de la conduite
D
:
diamtre
extrieur
2

D 22 + D12

1+ 2

D 22 D12
Les termes en

(
(

)
)

x
x

2 3 ) Relation entre les fonctions donde F t et f t + :

a
a

Nous admettrons ce qui parat intuitif, mais que nous justifierons dans la
suite qu lextrmit amont B de la conduite, de longueur L la
surpression h est nulle tout instant. Soit :
L
L

h (L, t ) = F t + f t + = 0
a
a

L
m
=
t
+
En
posant
:
alors
D
a
u0
Rservoirr
B
2L

A
f (m ) = F m

x0

amont

aval

Cette identit est valable, quel que soit m, en particulier lorsque m = t +

L
.
a

L
2L x

Soit : f t + = F t
. Ceci a pour consquence, en remplaant
a
a

E.H.T.P.
Zorkani Mohammed
7-7
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

L
L

f t + par son expression en fonction de F t dans les 2 quations

a
a

p/
x
2L x

h=
= F t F t

g
a
a

dAllievi :
g
x
2L x

u / = u 0 F t + F t

a
a
a

En posant = 2L a : cest le temps mis par une onde pour effectuer un

aller et retour dans la conduite et [en x = 0 {distributeur (vanne)} ] on a
alors :
(5)
h = F(t ) F(t )
g
u / = u 0 [F(t ) + F(t )] (6 )
a
On constate que dans ces 2 groupes dquations une seule fonction F est
inconnue.
2 4) Etude dune fermeture totale instantane :
N. B. : Le mode opratoire, qui est dcrit est gnral : il est valable lors
dune fermeture, ou dune ouverture, lente ou rapide, dun distributeur :
On tudie ce qui se passe au distributeur pendant le temps t compris entre
0 et , puis et 2 LL La conduite tant initialement parcourue par
Q
courant uniforme : u 0 = 0
, on ferme compltement le distributeur au
s0
temps : t = 0 . Nous nous proposons tout dabord dtudier comment varie la
surpression h lextrmit aval de la conduite.
Daprs les quations 5 et 6 dAllievi, on a pour des instants t 0 :
h = F(t ) F(t )

au 0
g = [F(t ) + F(t )]

Remarquons tout dabord que londe de retour [F(t )] napparat au

2L
.
distributeur quau bout du temps : =
a
Si t1 dsigne un temps quelconque de lintervalle 0 , daprs la seconde
des quations ci dessus : F(t1 ) = au0 g
Connaissant ainsi F(t ) dans le premier intervalle de temps [0, ] , celle ci
est dduite nimporte quel instant. On a en effet, successivement :
E.H.T.P.
Zorkani Mohammed
7-8
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

au 0
: t 2 = t1 +
F(t1 )
g
au
F(t 3 ) = 0 F(t 2 )
: t 3 = t 2 + = t1 + 2
g
LLLLLLLLLLLLLLLLLLLL
au
F(t n +1 ) = 0 F(t n ) : t n +1 = t n + = t1 + n L L
g
F(t1 + ) = 0
F(t 2 ) =

F(t1 + 2 ) =

Soit :

au 0
g

1 + ( 1)n au 0

L ETC.
L
F(t1 + n ) =
2
g
Do le graphe de la fonction F(t ) , et par suite celui de la fonction F(t - ) :

F(t )

au 0
t
F(t - )

au 0
g
u distributeur : x = 0

h x = 0 = h (x = 0, t )
T0 = 2

2L
=
a

Fig1 : Fermeture instantane dun distributeur (vanne)

E.H.T.P.
Zorkani Mohammed
7-9
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

 x
F t -
a

x
a

2L - x
a

au 0
g

t
2L - x
F t
a

x
a

x
a
au 0
g

h (x , t )

la distance : x 0

t
x
a
u/

x
a
+ u0
t
u0

Fig2 : Fermeture instantane dun distributeur (vanne) [suite : Fig1]

On a dduit immdiatement, la courbe reprsentant les variations de
surpression : h x = 0 = F(t ) F(t ) .
On observe ainsi une srie dintervalles successifs , pendant lesquels h
a u0
a u0
passe alternativement de la valeur
la valeur
; la priode
g
g
4L
doscillation du phnomne est :
T0 = 2 =
a

E.H.T.P.
Zorkani Mohammed
7-10
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Exprime en hauteur deau, la pression au distributeur varie donc entre les

a u0
a u0
limites : p 0
et p 0 +
( les variations de pression sajoutent
g
g
algbriquement la pression rgnante normalement en chaque point M considr )
a u0
Le terme
a habituellement une valeur trs grande, pour
g
a = 1000 m/s a u 0
= 500m d' eau 50 Atmosphre

u 0 = 5 m/s
g
Il est donc craindre, dune part, que la conduite et son quipement ne
rsiste pas leffet de telles surpressions et, dautre part, que le
phnomne de cavitation apparaisse du fait des dpressions. Cherchons
maintenant dterminer comment voluent en fonction du temps, la
vitesse u et la surpression h, en un point quelconque M, dabscisse x. On a
vu que :
x
2L x

h (x, t ) = F t F t

a
a

g
x
2L x

u / (x, t ) = u 0 F t + F t

a
a
a

La fonction F(m ) est connue.

x
2L x
, on obtient alors les
En posant successivement m = t et m = t
a
a
x
2L x

variations de F t et de F t
. On en dduit aprs les
a
a

x
2L x
translations respectives et
de laxe des abscisses, les
a
a

x
2L x

courbes donnant les variations de F t et F t

, en fonction de
a
a

t , en un point M dabscisse x(voir figure 2) : do le graphe de h (t ) et u (t )

au point M. On constate que la priode du phnomne est toujours 2 et
a u0
a u0
que la surpression h varie toujours entre
et +
. Par contre, le
g
g
temps pendant lequel la conduite est en surpression, ou en dpression,
2(L x )
2x
nest pas , comme au distributeur , il se rduit
=
.
a
a
E.H.T.P.
Zorkani Mohammed
7-11
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Ainsi toutes les sections de la conduite, lextrmit amont excepte, sont

a u0
.
donc soumise au coup de blier maximum, damplitude +
g

Nous avons admis jusquici que la fermeture du distributeur (vanne) tait

instantane. Ce qui est physiquement inacceptable. Les rsultats que nous
avons obtenus sont dailleurs entachs dune incertitude puisque nous
observons une onde de choc et que, dans ces conditions, lhypothse de
continuit, admis pour tablir les quations dAllievi, nest plus valable.
Il convient donc de reprendre ltude prcdente en considrant que la
fermeture du distributeur seffectue pendant un temps fini. Deux cas
2L
seront considrer suivant que est infrieur ou suprieur =
.
a
En effet :
Dans le premier cas [ ], londe F atteint sa valeur maxima au
distributeur avant que londe rflchie f , de signe contraire,
natteigne celui ci : des lors, on peut prvoir que le coup de
blier sera particulirement brutal, du moins au distributeur.
Dans le deuxime cas [ ], londe F naura pas atteint sa
valeur maxima au moment o londe rflchie f se manifeste au
distributeur : le coup de blier sera certainement moins violent
que prcdemment.

Q m 3 / s = S.v = Se .v e

section : S

section de sortie Se (vitesse v e )

Pour turbine Pelton

pointeau rglable
DISTRIBUTEUR : injecteur

2 5 ) Etude dune fermeture brusque et complte :

On dit quune manuvre de fermeture est brusque lorsquelle se produit en
un temps , infrieur . La fermeture est complte si elle correspond un
arrt total du dbit, dans le cas contraire on dit quelle est partielle.
Supposons connue la loi donnant, durant le temps , la vitesse u au
distributeur ; daprs les hypothses faites on a donc :
E.H.T.P.
Zorkani Mohammed
7-12
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

u = u 0 1 (t )

pour 0 t

(t ) tant une fonction connue , variant de 0 1 lorsque t varie de 0 ;

u = 0 pour t . Si t1 dsigne un temps infrieur , londe F(t ) nest
pas encore au distributeur, de sorte que lon a, daprs lquation (6) :
au
F(t1 ) = 0 (t1 ) pour t1
g
au
F(t1 ) = 0
pour t1
g
Soit maintenant, un instant t 2 = t1 + de lintervalle , 2 . Comme u est nul
au
F(t 2 ) = 0 [1 (t 2 )] pour t 2 +
g
au distributeur, on a simplement :
F(t 2 ) = 0
pour + t 2 2
Do le graphe suivant de la fonction F(t ) dans lintervalle 0 2 :
F(t )

F(t - )
t

h x = 0 = h (x = 0, t )

2L
=
a

au 0
g

2
t

Fig3a Fermeture brusque et complte du distributeur : x x c

Pour t 2 la courbe reprsentative de la fonction F(t ) sobtient

immdiatement en remarquant que cette fonction est priodique, de
priode 2 ; on a en effet, pour t :
a u0
a u0
= F(t ) + F(t ) ,
= F(t + ) + F(t )
g
g
et par suite :

E.H.T.P.
Zorkani Mohammed
7-13
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

F(t + ) = F(t )
La fonction F(t ) tant ainsi connue, quel que soit t, on en dduit
immdiatement les fonctions F(t ) et h x = 0 = F(t ) F(t ) .
au
au
On constate que la surpression au distributeur varie entre 0 et + 0
g
g
comme dans le cas dune fermeture instantane. Par contre, le temps
pendant lequel la contrainte est maxima se trouve rduite de .
Pour dterminer la surpression h dans une section dabscisse x, la
mthode est identique celle que nous avons utilise dans le paragraphe
prcdent, nous tracerons, au pralable, les courbes reprsentant les
te
2L x
x

variations de F t et F t
en fonction de t x = C .
a
a

2[L x]
Deux cas sont considrer selon la valeur du dphasage
de
a

ces courbes.
2[L x]
a
1. Supposons que lon ait :
soit : x L
a
2
Les courbes reprsentes sur les figures (3a) au dessus,
montrent que lamplitude maximum h max du coup de blier est
a u0
; le temps pendant lequel sexerce cette
encore gale
g
2x
contrainte maxima est gal :
.
a
a
2. Supposons maintenant que lon ait : x L
2
Dans ce cas, lamplitude maximum du coup de blier est
a u0
infrieure
; h max diminue progressivement depuis la valeur
g
a u0
a
, jusqu zro , lorsque x augmente de x c = L L
g
2
(voir figure 3b au dessous ) . La loi de variation h max (x ) , lie
la loi de fermeture du distributeur (vanne), pourra tre tablie
dans chaque cas particulier.

E.H.T.P.
Zorkani Mohammed
7-14
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

 x
F t -
a

x
a

F t - -
a

h (x )

au 0
g
t

x
a

2
t

2L
=
a

x
a
Fig3b Fermeture brusque et complte du Distributeur : x x c
Nous avons donc mis en vidence lexistence dune section particulire
(S c ) dabscisse : x c = L a
2
Pour x x c , le coup de blier est aussi important , bien que plus faible
dure , que dans le cas dune fermeture instantane. Par contre pour x x c
on observe une rduction de lamplitude de ce coup de blier.
Ces rsultats auraient pu tre rtablie immdiatement. En effet, dans une
section dabscisse x, la premire onde de surpression atteint sa valeur
x
tandis que londe de dpression correspondante
maxima linstant +
a
2L x
.
se manifeste, dans cette mme section, linstant
a
Dans ces conditions, le coup de blier ne pourra tre induit que si lon a :
2L x
x
a

+
x x c avec x c = L

a
a
2

Sur les figures (4a 4b ) nous dessinons les courbes correspondant au cas
o la loi de fermeture du distributeur (vanne) serait de la forme :
E.H.T.P.
Zorkani Mohammed
7-15
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

u
t

t
u = u 0 1 pour t linaire

Les premires courbes correspondant aux sections pour lesquelles x x c et
les secondes aux sections pour lesquelles x x c .
x

x
F t -
a
a

F t - -
a

h (x )

x
a
x
F t -
a

2L
a

x
a

Fig4a Fermeture linaire d' un Distributeur : x x c

F t - a

h(x )

L - x au 0
L - xc

t
x
a

x
a

Fig4b Fermeture linaire d' un Distributeur : x x c

E.H.T.P.
Zorkani Mohammed
7-16
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

hmax
au0
g

Fig4c

xc

Dans une section dabscisse x suprieure x c lamplitude h max (x ) du coup

de blier a pour expression :
2 (L x ) a u 0
1
varie en
h max =

a

g
Soit :
L x a u0
a

o x c = L
h max =

L xc g
2

Depuis la section (Sc ) jusqu lextrmit amont de la conduite, la

rpartition du coup de blier est donc linaire (voir Figure 4c).
2 6 ) Etude dune fermeture lente et complte :
Nous venons de monter que si , toutes les sections dabscisse x x c
a u0
sont le sige de surpression et de dpression dont lamplitude
a
g
gnralement une valeur inadmissible.
En pratique, on est donc amen nenvisager que des fermetures lentes
cest dire dont la dure est suprieure .
Bornons-nous, ici, tudier le cas particulirement simple pour lequel la loi
de fermeture du distributeur (vanne) est donne par les relations :
t

u = u 0 1 pour t

u=0
pour t

tant un multiple de , que nous crirons : = k k tant un entier 1

Nous cherchons, tout dabord, dterminer la fonction F(t ) pour t .
E.H.T.P.
Zorkani Mohammed
7-17
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

En dsignons par t1 un instant quelconque de lintervalle [0 , ] on a

successivement :
t au
F(t1 ) = 1 0
g

t + t1 au 0 au 0
F(t1 + ) = 1

=
g g

t + 2 au 0 t1 au 0
F(t1 + 2 ) = 1

= +
g g

t + 3 t1 au 0 2 au 0
F(t1 + 3 ) = 1

=

g
g

LLLLLLLLLLLLLLLLLLLL
n au 0
F[ t1 + (2n 1) ] = LLLLLLLL =

g
n au 0
t
F[ t1 + 2n ] = LLLLLLLLL = 1 +

LLLLLLLLLLLLLLLLLLETC.
Cette suite doit tre arrte lune des relations suivantes selon la parit
de k =

p au 0 1 au 0
F[ t1 + (2n - 1) ] =
pour k = 2p
=

2 g
g
t p au 0 1 t1 au 0
F[ t1 + 2n ] = 1 +
= +
pour k = 2p + 1

g 2 2 g

Pour t on a :
au
F(t ) + F(t ) = 0
g
Et par suite selon la parit de k :

E.H.T.P.
Zorkani Mohammed
7-18
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

(a ) F[ t1 ] =

1 au 0

2 g

pour k =

= 2p

ou
1 au 0 t1 au 0

2 g
2 g
au
1 au
t
F[ + t1 + ] = 0 + 1 0
2 g
2 g
LLLLLLLLLLLLLLL

(b ) F[ + t1 ] =

au

1 au 0
t

( 1) p 1 0 pour k = = 2p + 1

2 g
2 g
Nous somme maintenant en mesure de tracer la courbe reprsentant F(t )
et den dduire comment varie la surpression h x = 0 au distributeur, et la
surpression h x dans toute section dabscisse x.
Les deux courbes reprsentes sur les figures (5a) et (5b) correspondant
respectivement aux 2 cas : = 2p et = (2p + 1)
F[ + t1 + p ] =

E.H.T.P.
Zorkani Mohammed
7-19
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

1 au 0

= 2k

F(t )

k N

2 g

au 0

1 au 0
2 g

F(t - )

au 0

g
2

h x = 0 (t )
au 0

1 au 0
2 g

x
F t -
a

x
a

au 0

F t

-x
2L

(2L - x)
a

h (x )

t
Figure (5a) : = 2 k
E.H.T.P.
Zorkani Mohammed
7-20
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

au 0 +
g 2

= (2k + 1)

F(t )

k N

au 0

F(t - )

au 0

g
2

h x = 0 (t )
+

au 0

au 0
g

au 0
2
g

F t -

x
a

au 0

F t -

2L - x
a

(2L - x)
a

h (x )

t
Figure (5b) : = (2 k + 1)

E.H.T.P.
Zorkani Mohammed
7-21
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

On constate que, quelle que soit la parit de k, lamplitude du coup de

blier est donne par les relations suivantes :
Au distributeur : x = 0

a u 0 1 au 0 2L u 0
=
=
g
k
g
g
Dans une section d' abxcisse ' ' x' ' :
h max =

avec : u 0 =

Q0

p max (x )
x 2L u 0
= 1
g
L g
Ce sont les formules de Michaud, qui demeurent dailleurs valables
h max (x ) =

mme si k =

nest pas un nombre entier. Elles prsentent lavantage de

ne pas faire intervenir la clrit a sur laquelle rgne une certaine

incertitude du fait des hypothses admises. En particulier, celles qui sont
relatives la dformation de la conduite. [ J. Michaud : coup de blier dans
les conduites, Bull. de la Soc. Vaudoise des Ing. et Arch. , 1878].
Ainsi ces formules sont elles trs utilises en pratique, dans tous les cas
o la fermeture du distributeur (vanne) est progressive, bien quen toute
rigueur, elles ne doivent tre appliques que si la loi de fermeture est de la
t
forme linaire :
u = u 0 1

2 7 ) Dtermination exprimentale de la clrit a par la mthode
de la dpression brusque :
Les hypothses admises pour calculer a ne sont gnralement pas
valables dans le cas dune installation industrielle. Il est dont prudent de
dterminer exprimentalement la valeur de a.
Monsieur Camichel a propos la mthode suivante :
Le distributeur tant initialement ferm, et si le fluide tat au repos, on cre
une petite dpression en ouvrant et en refermant un robinet de purge R
dispos lextrmit aval de la conduite ; la manuvre dure un temps trs
court (voir figure 6).
R Robonet

Figure(6 )

capteur
enregistreur

E.H.T.P.
Zorkani Mohammed
7-22
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Le passage des ondes successives est enregistr grce un capteur de

pression dispos, lextrmit aval de la conduite. On obtient ainsi un
enregistrement dont les caractristiques sont indiques sur la figure (7) :

h x =0
2a
95

Q max
a
95

a
95

2a
95

Q max

Q max

Q max

Figure (7 ) : Enregistrement

( Q max dsigne la valeur maxima du dbit vacu durant la manuvre du

robinet R . La mesure se dduit immdiatement de celle de celle de , ou
mieux, de lun de ses multiples).
2 8 ) Phnomne de rsonance dans les conduites :
Il est possible de faire entrer en rsonance une conduite, dont (la vanne) le
distributeur est ferm, on disposant lextrmit aval de celui ci un
robinet que lon fait tourner vitesse variable.
La rsonance stablit lorsque la priode T des ouvertures et des
4L
fermetures du robinet devient gale la priode propre 2 =
de la
a
conduite. On dmontre sans difficult qu la fin du rgime transitoire, la
pression au distributeur varie de 0 2p 0 o p 0 est la pression statique, et
que le dbit du robinet est constamment nul.
Lexprience confirme ces rsultats ; cependant le dbit nest pas
rigoureusement nul cause des frottements.
La courbe reprsente sur la figure (8a) en dessous donne les variations
de lamplitude de la surpression le long de la conduite ; celle ci vibre en
quart donde. A fin quil ait rsonance la conduite il faut que lorgane
dobturation ouvre et se ferme de manire rythmique.

E.H.T.P.
Zorkani Mohammed
7-23
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

h max

1 = 4L fondamental

h max

Fig(8a)

Fig(8b)
L

h max

2 = 2L

L
h max

4
3

3 = L

harmoniques

4 = L
Fig(8d)

Fig(8c)

Les figures(8) : Modes de mises en rsonance

En augmentant la vitesse de rotation du robinet, de manire atteindre des
valeurs multiples de la prcdente, on peut mettre galement en rsonance
les harmoniques du fondamental. Pour lharmonique 2, la conduite vibre en
demi onde ; pour lharmonique 3, elle vibre en trois quart donde
(figures 8a 8d).
2 9 ) Clapet automatique, moteur hydraulique synchrone :
Remplaons le robinet tournant, dispositif laval de la conduite, par un
clapet tar de telle sorte quil soit appliqu sur son sige sous laction de la
pression statique p 0 (voir figure 9) :

Figure (9): Clapet

Appuyons sur lextrmit de la tige du clapet de manire dcoller de son
sige et laissons stablir un rgime permanent pour lequel la vitesse dans
la conduite est u 0 . Au moment o on lche le clapet, celui ci est
brusquement rabattu sur son sige. Londe de surpression qui en rsulte
E.H.T.P.
Zorkani Mohammed
7-24
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

au 0
) , maintient le clapet ferm pendant le temps . Dans
g
au
lintervalle de temps 2 londe de dpression ( 0 ) se manifeste au
g
distributeur. Si le clapet est convenablement tar, cest dire si, pour la
au
pression : p 0 0 [ ici ( g) est le poids spcifique du fluide qui
g
remplit la conduite ] il nest plus appliqu sur son sige, il se produit, dans
cet intervalle de temps un coulement qui amplifie dailleurs la dpression ;
le clapet se referme alors brusquement, et ainsi de suite.
Cette proprit a t mise profit, par Monsieur Camichel, pour raliser un
moteur hydraulique synchrone dont le principe est le suivant :
Au voisinage du clapet, que vient de dcrire, on monte, sur la conduite un
cylindre lintrieur duquel coulisse un piston reli une roue par un
systme de bielle manivelle (voir figure 10 au dessous).
dintensit (

roue
piston
ressort

clapet

Figure(10): Moteur Hydraulique Synchrone

Le piston est mis en mouvement du fait des surpressions et des
dpressions provoques par le mouvement du clapet.
2 10 ) Mthode graphique de Bergeron : Les caractristiques
2 10 1 ) Principe de la mthode :
Les quations de dpart tant :
E.H.T.P.
Zorkani Mohammed
7-25
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

x
x

h = F t + f t +
a
a

x
g
x
u = u 0 F t f t +
a
a
a
En y introduisant les dbits volumiques que lon a diffrents temps. Soit :
Q = S u et Q 0 = S0 u 0 , La deuxime quation devient :
a
(Q Q 0 ) = F f

g S
La premire quation devient
a
(Q Q 0 ) + 2F ou h = a (Q Q 0 ) + 2f
h=
g S
g S
Ces 2 quations sont susceptibles dune interprtation graphique.
En effet, soit h et Q la surpression et le dbit linstant t labscisse x, h
est lintersection des 2 droites en fonction de Q.
caractristique

h (x, t )

F=C

Les droites caractristiques :

te

f =C

(Q Q 0 ) + 2F
h
=

g S

h = a (Q Q ) + 2f
0

g S

te

caractristique

Q
Q(x, t )
On remarque que le lieu des points correspondant F = Cte est une droite
a
de pente tg =
, le lieu des points correspondant f = Cte est
g S
a
galement une droite de pente tg =
, et dans les deux cas , pour
g S
diffrentes valeurs de F et f , on a affaire des faisceaux de droites de
a
pentes
[ dite pure de Bergeron ].
g S
Remarquons ensuite que si la chambre de mise en charge (rservoir) est
de grande dimension, le lieu des points reprsentant ltat de la section
amont de la conduite est laxe des abscisses (h = 0 ) .
E.H.T.P.
Zorkani Mohammed
7-26
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

En outre, le lieu des points reprsentatifs de lextrmit lorsque du

distributeur (vanne) est ferm, est Q = 0 (laxe des ordonnes).
La mthode consiste donc suivre pas pas lvolution des fonctions F et f
dans le temps et dans lespace. Pour cela, il faut se rappeler que :
Si une onde daval prend la valeur F linstant t, labscisse x,
cette mme onde garde la mme valeur linstant (t + t)
x
labscisse (x + x) avec t =
.
a
Si une onde damont prend une valeur f linstant t , labscisse
x
x , cette mme onde garde la mme valeur au temps t = t

a
labscisse (x x) .
3- COMPLEMENT SUR
Rflexion et transmission dondes
Quand il se produit un changement de section dune conduite ou une
variation de ses caractristiques physiques (paisseur, matriau
responsable du changement de la clrit des ondes dans le rseau
industriel) : il en rsulte une rflexion partielle des ondes principalement au
niveau de ces discontinuits [une mthode consiste dcomposer la
conduite en plusieurs tronons ayant des caractristiques constantes
(application pratique de la mthode de Bergeron)]. A titre dexemples on va
prsenter les cas de changement de section et celui dun trs grand
rservoir de mise en charge :
3 1 Cas dun changement de caractristique de la conduite :

S1

S2

S0

On peut crire les quations dAllievi de la manire suivante :

h1 = F1 + f1
h 2 = F2 + f 2
et
u1 = u 01

Or : u =

g
(F1 f1 )
a1

et

u 2 = u 02

g
(F2 f 2 )
a2

Q
et par consquent :
S

E.H.T.P.
Zorkani Mohammed
7-27
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

g S1
(F1 f1 )
a1

g S2
(F2 f 2 )
a2
Q = Q 2 conservation de la masse
Comme dans la section S0 on a : 1
h1 = h 2 continuit de la pression
On en dduit alors que :
F2 + f 2 = F1 + f1
Q1 = Q 0

et

Q2 = Q0

S2
(F2 f 2 ) = S1 (F1 f1 )
a2
a1

Exemple dapplication :
Soit le temps compris entre 0 et t 2 pour lequel londe de retour f 2 ne sest
pas encore manifeste alors :
F2 = F1 + f1
S2
S
F2 = 1 (F1 f1 )
a2
a1
On en dduit, en liminant F2 ou f1 des 2 quations :
a S
1 1 2
S1 a 2
2
F2 =
F1 et f1 =
F
a1 a 2
a1 S2 1
1+
1+
S1 S2
S1 a 2
b
a
a
En posant :
b1 = 1 , b 2 = 2 et m = 2 alors
S1
S2
b1
2m
m 1
F2 =
F1 et f1 =
F1
m +1
m +1
Considrons maintenant le temps t t 2 : londe de retour F2 est connue. Il
suffit de respecter les conditions aux limites et de procder de proche en
proche.

E.H.T.P.
Zorkani Mohammed
7-28
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

3 2 Rservoir de grande dimension :

v
v p
=0
+ v
+
t
x x

L v
1
v 2x = 0 = (p 0 p) dx
2
0 t

Rservoir

F1

f1

f
1

D
A

amont : p
L

F1

aval: p 0 , v x = 0

x=0

On peut assimiler le rservoir comme une conduite dans la section est trs
S2
grande, cest dire que le rapport
.
a2
Dans ce cas : F2 = 0 , f1 = F1
Le coup de blier ne se transmet pratiquement pas dans un rservoir. Il est
donc lgitime dadmettre que la surpression dans la section B est nulle.
Dautre part, dans le cas o la section S2 nest pas infinie, ce qui est
gnralement le cas, le rsultat acquis reste valable, si toutefois les
manuvres sont effectues lentement.
Soit F2 , dans ce cas , une onde damplitude trs faible , issus au temps t de
la section B et qui se rflchit la surface libre du rservoir pour donner
t
f 2 = F2 au temps t + .
2
Apparat alors une onde f 2 (t + t ) = F2 (t ) en B.
La surpression en B devient : h B (t + t ) = F2 (t + t ) = F2 (t )
h
La cote h tant faible, t = trs faible est que linfluence en B des
c
dimensions finies du rservoir est nulle.
Remarque :
Dans le cas dune chemine dquilibre une tude dtaille devient
indispensable. Cest ce quon va tudier maintenant {$4}.

E.H.T.P.
Zorkani Mohammed
7-29
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

[Link] dquilibre
[ Oscillations en masse (en bloc) ]

Une chemine dquilibre est un rservoir de section S finie.

Un systme hydraulique fondamental est schmatis par :

la retenue

chemine d' quilibre

La galerie en charge
S

organe d' obturation

la conduite force

Sytme Hydraulique fondamental

La galerie en charge est sensible aux oscillations de pression. Le rle de la

chemine dquilibre est de la protger, en interceptant les ondes de
pression provenant de la conduite force. Grce linstallation de cette
chemine on verra dans la suite quelle rduit les oscillations de pression
dans la galerie en charge, on signale que la conduite force doit, dans la
plupart des cas, tre construite de manire quelle puisse supporter des
coups de blier considrables.
4 1 ) Gnralits de mise en quation :

Un rservoir, gnralement appel chambre dquilibre ou chemin

dquilibre, judicieusement plac dans la conduite permettra dviter ou au
moins damortir les coups de blier.
Ltude des chemines dquilibre repose sur les phnomnes doscillation
en masse qui prennent naissance par suite de la prsence dune surface
dformable : la surface libre.
La figure suivante prsente les lments dune chemine dquilibre
simple.
E.H.T.P.
Zorkani Mohammed
7-30
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

pa

pa

origine des cotes

S
chemine
conduite

Rservoir
A

V=

dz
dt

V 0 si haut
bas

v 0 de A vers C
C
s

Q
L
Mise en quation :
Egalit des dbits au nud : quation de continuit
dz
s v =SV + Q =S + Q
dt
Thorme de Bernoulli entre A et B : quation dynamique
r
p
v2 1 v 1 r r
r r

+
+

+
+

rot
v
v
lap
v
z
grad
grad

dl = 0
g

g

2g
g
t
g

r
v2 1
v
d l + PCH
0=z+
+
2g g (A B ) t
Formule dans laquelle PCH reprsente les pertes de charge rgulires et
singulires entre A et B [exprime en mtre du fluide].

r
v et d l : sont toujours colinaires.
La vitesse V (de C B ) tant faible (S : chemine de grande section ) .
Les oscillations en masse, devant avoir une priode grande, telle que la
vitesse v soit, chaque instant, la mme dans toute la conduite (ceci devra
tre vrifi par la suite).
Lquation scrit alors
L dv
+ z + PCH = 0
g dt
Il est dusage courant dcrire que PCH , les pertes de charge sont de la

forme quadratique : PCH = K v 2

Lorsque v 0 , leau allant vers lutilisateur [ c--d leau coule de la retenue
(rservoir) vers la chemine] , lquation de Bernoulli est alors
E.H.T.P.
Zorkani Mohammed
7-31
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

L dv
+ z + Kv 2 = 0
g dt
Lorsque v 0 , leau refoulant dans la conduite[ c--d leau coule de la
chemine vers la retenue (rservoir)], lquation de Bernoulli est alors
L dv
+ z Kv 2 = 0
g dt
Souvent on pose Q = S U afin de donner lquation de Bernoulli une
forme plus homogne.
En rsum : les 2 quations notre disposition sont
s v = S (V + U)

L dv
avec = 1 si v 0
+
z
+

P
=
0

CH
g dt
= 1 si v 0

4 2 ) Rpartition des surpressions le long de la conduite L :

Ecrivons Bernoulli entre la surface libre du rservoir et une section s
dabscisse l , en utilisant les mmes hypothses que prcdemment et
avec v 0 et Pl_CH tant les pertes de charge entre la cote A et le point l .

v2 l d v

+ z +
+
+ Pl_CH = 0
g
l 2 g g d t
La surpression est
p

, soit :
+ z
+ z
+ z
g
g
l, v = 0
l g
l _ lorsque v est nul

ainsi
l d v v2

p = p l p l_v = 0 = g
+
+ PCH
g d t 2g

En tirant

L dv
dv
de lquation gnrale
+ z + PCH = 0 on obtient :
g dt
dt
l

dv
g
v2
l
+ PCH
= (z + PCH ) p = g z PCH +
dt
L
2g
L
L

p devient :

v2
l
p = g z
1 PCH
2g L
L

E.H.T.P.
Zorkani Mohammed
7-32
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

La figure suivante nous renseigne sur PCH qui comprend :

H A : perte de charge lentre de la conduite.
H RAC : perte de charge linaire (rgulire) proportionnelle l .
H C : perte de charge singulire la drivation en C, pied de la
conduite.
Perte de charge ( dbit constant ) avec : v 0

H C

H RAC

Pl_CH

H A

PCH

A
= H A + H RAC + H C

Pl_CH = H A +

l
H RAC
L

v2
l
p = g (z + H C )
1 H A
2g L
L

La loi p(l ) est une loi linaire, la surpression la plus dangereuse se situe
en C, pour l = L :

(p )l = L = g z + H C v

2g

Remarque : Habituellement le canal damene dbouche directement dans

r r
v2 v v
la chemine dquilibre, de sorte que : H C =
=
. Dans ce cas on
2g 2g
peut crire : (p )l = L = g z .
4 3 ) Emploi de variables adimensionnelles :
Il sagit dexposer la mthode de Calame et Gaden, qui consiste donner
aux quations gnrales un caractre universel, et dont la justification sera
conduite, lors des exemples tudis ultrieurement.
Soit U 0 vitesse correspondant au dbit nominal Q 0 dans la section S.
Soit v 0 vitesse correspondant au dbit nominal Q 0 dans la section s.
T0 = 2

L S
2
correspondant 0 =
=
g s
T0

g s

L S

E.H.T.P.
Zorkani Mohammed
7-33
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

 z = v 0

L s
= U0
g S

L S
, Q 0 = SU 0 = sv 0 ( rgime permanent )
g s

z est la dnivellation maximum la fermeture brusque de la vanne.

2

T0 z
L S T02 z
Entre U 0 , 0 et z existent les relations suivantes : = 2 = 2 =

g s 4
U 0 2 U 0
Les 2 quations rappeles ici, deviennent :
s v = S (V + U)

L dv
avec = 1 si v 0
+
z
+

P
=
0

CH
g dt
= 1 si v 0

dv S dV S dU S dV
S dU
En drivant la premire quation :
=
+
=
V+
dt s dt s dt s dz
s dt
L S dV L S dU
Et en reportant dans la seconde :
V
+
+ z + PCH = 0
g s
dz g s dt

En remplaant

z
LS
par 2 dans le premier terme.
g s
U0
T0 z
LS

par
dans le deuxime terme.
g s
2 U 0

V
U

d
d
V U 0 1 U 0
Lquation devient :
z + z + PCH = 0
z +
2 t
U0 z

d
d
T
0
z
En divisant par z , il apparat les termes { sans dimension } :
PCH
t
/ z
/
/
z = z , PCH = z , = T

V/ = V , U/ = U , v/ = v

U0
U0
v0
La prime ' / ' signifie sans dimension (adimensionnel).
Les 2 relations scrivent :

E.H.T.P.
Zorkani Mohammed
7-34
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

 V/

d V/
d z/

1 d U/
avec = 1 si v 0
/
/

+
z
+

P
=
0

CH
2 d T /
= 1 si v 0

v/ = U/ + V/
Sous cette forme rduite, les quations ne renferment plus quun seul
/
paramtre PCH
, on conoit alors que cest celle ci que lon a recours
pour dresser des tables et abaques pouvant tre utilises pour nimporte
quelle installation.
Remarque :
dz
1 d z/
/
Sachant que : V =
et par suite V =
on trouve que :
2 d T /
dt
d V/
1 d2 z/
V
=
d z / 4 2 d T / 2
/

Enfin les pertes de charge de la forme quadratique Kv 2 imposant que les

0
pertes de charge en rgime permanent sont de la forme : PCH
= Kv 02 et par
suite : PCH = Kv =
2

Kv 02

v
0
= PCH
v / 2 . En divisant par z , il vient alors :
v0

/0 /2
/
PCH
= PCH
v

/0
avec PCH
=

0
PCH

z
Ainsi les quations gnrales peuvent aussi avoir la forme :
1 d 2z /
1 d U/
/0 /2
v =0

+
+ z / + PCH
2
/2
/
2 d T
4 d t
v/ = U/ + V/
Ces 2 quations, sous cette forme, permettront la rsolution dune
fermeture totale instantane en tenant compte des pertes de charge(voir
ltude analytique dans la suite).
4 4 ) Etude dune fermeture instantane :
Hypothse de travail : pertes de charge ngliges
Le rgime permanent correspond au dbit Q a tabli pour t 0 .
A linstant initial t = 0 , on modifie le dbit Q en lui donnant la valeur Q b . Les
quations du mouvement et les conditions initiales sont :
L dv
+ z = 0 avec t 0 : Q = Q a

g
d
t

t = 0 : Q = Q b
sv = SV + Q

E.H.T.P.
Zorkani Mohammed
7-35
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

sv (t = 0 ) = Q a

sv (t = 0 ) = SV(t = 0 ) + Q b
En liminant v entre les 2 quations du mouvement nous obtenons :
dv
dV
d 2z
=S
=S 2
s
dt
dt
dt

2
L S d z + z = 0
g s d t2

La rsolution gnrale donne : z = A sin ( 0 t + ) avec

t = 0 : z = 0 = 0
2
g s

0 =
=
comme :
1
T0
LS
(
)
(Q a Q b )
=
=
=

t
0
:
v
A
sin
t
0
0
0

s
Q
Q
Ls
En dsignant par : v b = b , v a = a ; A = (v a v b )

s
s
gS
Nous tirons :

z = 0 pour t 0 et

z = (v a v b )

Ls
sin
gS

Ls
Ls
2
t = (v a v b )
sin
t
gS
gS
T0

Remarque :
Lamplitude des oscillations, non amorties, car les pertes de charge ont t
ngliges, autour du niveau statique, est uniquement fonction de (Q a Q b )
z
z

Q b Q a
T0 = 2

t
Qb Qa
t

LS
g s

vec pertes de charge ngliges

manuvre brutale d' un distributeur :

oscillations du niveau de la chemine d' quilibre

E.H.T.P.
Zorkani Mohammed
7-36
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

 Calcul de la section minimale dune chemine simple :

Soit p max la surpression maximale tolrable lors de la fermeture complte
et instantane du distributeur (la vanne).
Ls
soit :
gS

p max = g z max = g v a

L v a2
S=s
g z 2max

S
z max

On constate que la section de la chemine (chambre) dquilibre est

inversement proportionnel au carr de la surpression.
L v a2
Quant au volume : = 2 z maxS = 2 s
, il nest inversement
g z max
p max
proportionnel qu la surpression : z max =
.
g
4 5 ) Etude dune fermeture linaire dun distributeur :
Hypothses de travail :
Q
Q = Q 0 pour t = 0
Q0

Q = Q 0 1 pour 0 t

Q = 0 pour t
Posons

avec T0 = 2

T0
4 5 1 ) Etude pour 0 t :

LS
gs

Lquation de continuit est : s v = S V + Q 0 1

Q
t

s v = S V + U 0 1 avec U 0 = 0
ou
S

dv S d 2 z U 0
=

Drivons par rapport au temps t :

dt s dt 2

L S d 2z U0
Reportons dans lquation de Bernoulli alors :

+z=0
g s dt 2

E.H.T.P.
Zorkani Mohammed
7-37
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

LS d 2z
L S U0
+
z
=
g s dt 2
g s
LS d 2z
La solution gnrale de
+ z = 0 est : z = A sin ( 0 t + ) laquelle il
g s dt 2
L S U0
faut ajouter la solution particulire : z p =
.
g s
Soit

U T
t
En posant z = 0 0 et =
lquation devient :
2
T0

2
z
T0 d z
+
=
z

2
2 dt 2
2

t
z
Dont la solution est : z = A sin 2
+ +
T
0
2

Les conditions initiales nous donnent :

t = 0 : z = 0 (pas de perte de charge )

ce qui entrane
dz

t
=
0
:
V
=
0
soit
=
0
0

dt

A = z

Ainsi

z
2 t

z
1
cos

2
T
valable pour

0 t

z
2
t

V = T sin T

0
0

z
z
(1 cos 2 ) = sin 2
z =

2
Donc, linstant t = on a :

2z
z
V =
sin 2 = sin cos
T0
T0

4 5 2 ) Etude pour t :
s v = S V

Les quations sont :

L dv
g dt + z = 0

Avec comme conditions initiales les valeurs particulires de z et v pour z

et v respectivement au temps .

E.H.T.P.
Zorkani Mohammed
7-38
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

2
T0 d z
Il en rsulte lquation diffrentielle :
+ z = 0 dont la solution est :
2 dt 2

t
z = A / sin 2
+ /

T0
/
/
Les 2 constantes dintgration A et sont dtermines par les conditions
initiales : il en rsulte alors :
/
z
/
sin 2
A sin 2 + =

2 A cos 2 + / = 2 z sin cos

T
T0
0

En divisant membre membre les 2 galits, il vient : tg 2 + / = tg

sin
Do : / = et / = z

Alors

z = z

sin

valable pour
sin 2 t
t
2

Les 2 formules donnant z(t ) pour t et t vont permettre de tirer z max

cote maximale atteinte par le plan deau dans la chemine :
T
z
Pour 0 t : z max =
et aura lieu pour t1 = 0 cette valeur ne

2
T
sera effectivement atteinte que si 0 .
2
T
sin
Pour t : z max = z
et aura lieu pour t 2 = 0 + , cette

2 2
T
valeur ne sera atteinte que si 0 .
2
Ainsi le maximum z max atteint par le plan deau dans la chemine aura pour
expression :
z
T
si 0

L s
LS

z max =
& =
o z = v 0
avec T0 = 2
g S
T0
gs
z sin si T0

2
E.H.T.P.
Zorkani Mohammed
7-39
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Remarque :
Le tableau ci dessous permet de se rendre compte de linfluence du
temps relatif de la fermeture linaire selon Calame et Gaden

0,05

0,1

0,2

0,3

0,996

0,984

0,935

0,858

0,4

0,5

0,6

0,8

z max

0,757 0,637 0,531 0,398

0,318

z
Dautre part, le temps habituellement est de lordre de quelques
secondes, alors que T0 est trs grand de lordre de quelques minutes.
Do , trs petit, ce qui implique gnralement
surlvation voisine de z .
z max

1 et alors une

Q0
U 0 T0

=
=
z
avec
:
U
0

2
S

L
S

T0 = 2 g s

1.00

0.80

0.60

0.40

0.20
0.00

0.20

0.40

0.60

0.80

1.00

Influence du temps de fermeture sur z max

NOTES SUR
Influence du choix du matriau dune conduite
de refoulement sur le coup de blier
- Dterminer : le choix dune conduite [diamtre, paisseur, matriau]
- Connatre : la pression de service {donne par le constructeur}
- Connatre : la nature et temprature du liquide
- Calculer : la surpression pour le dimensionnement etc.
E.H.T.P.
Zorkani Mohammed
7-40
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

 Matriaux compars :

(Mpa )
Matriau
200000
0,3
Acier
0,38
2758
PVC
90000
0,25
Fonte
170000
0,28
F.D.
1035
0,45
PE
Tendance dans le monde municipal :
Pour une conduite de refoulement, est-il prfrable de spcifier un
matriau plastique plutt que mtallique ?
aV
H =
= surpression en m
g
a : clrit de londe (en m/s) ; dpend du matriau et de son paisseur
V : vitesse de lcoulement (en m/s) en rgime permanent qui dpend du
diamtre.
La Clrit de londe est donne par :
K

a=

K D
1 + c1
E e

Paramtres en jeu :
K = : module de compressibilit du fluide transport
: masse spcifique de du fluide transport
E : module de Young du matriau de la conduite
D : diamtre de la conduite
e : paisseur de la paroi de la conduite
c1 : coefficient dencastrement [ancrage] (1 2)
: coefficient de Poisson du matriau (voir tableau au-dessus)
Coefficient dencastrement : Ancrage
- Conduite rigide : c1 = 0
-Conduite lastique paroi mince :
Ancre contre tout mouvement longitudinal sur toute sa
c1 = 1 2
longueur :
Encastre l amont seulement : c1 = 1.25
Avec de nombreux joints d expansion : c1 = 1
E.H.T.P.
Zorkani Mohammed
7-41
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Exemples Pratiques :
Surpression
Remarque : il s agit ici
dune des conduites
forces d une centrale
Hydro-lectrique.

Sous-pression !
Remarque: il s agit ici d une
section de la mme conduite
force, situe en amont de
la section rupture par
clatement.
DR : dimension ratio DR
: = D0
e

F = force

D0

e
p

pD

F
ne

e
On travaille par mtre linaire

pD

F = 2 = n e e

n = p D
e
2e

or

Ainsi
F
p D = 2F
Fx = 0
n e = contrainte
e
Contrainte circonfrencielle :
p
n e = (DR 1) Equation ISO
2

E.H.T.P.
Zorkani Mohammed
7-42
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

4 6 ) Etude de la fermeture instantane en tenant compte

des pertes de charge :
Immdiatement aprs la manuvre de fermeture, le dbit appel Q est nul
U / = 0 le niveau deau slve dans la chemine v / = V / . Tant que dure
cette lvation la vitesse v dans le canal reste positive, , indiquant le
signe des pertes de charge est gal + 1 .
Inversement lorsque le niveau dans la chemine dcrot : = 1 .
4 6 1 ) Examinons ce qui se passe durant la 1re monte
du plan deau :
On a tablit au paragraphe (4 3 : pp.34) que
d V/
1 d U/
avec = 1 si v 0
/
/
z
P
0
V/
+
+
+

CH
/
/
2

dz
dT
= 1 si v 0

v/ = U/ + V/
Pour la monte deau dans la chemine on a alors
/
1 d U/
/0 /2
/ dV
V
v =0
+
+ z / + PCH
/
/
2 d T
dz

U/ = 0 , v/ = V/
Lquation devient

d V/
/0 /2
V
+ z / + PCH
V =0
/
dz
/

En posant : V / 2 =

d
dz

/0
+ 2 PCH
= 2z /

2 P

/0

z/

CH
Lintgrale gnrale de lquation est : = C e
Il reste trouver une solution particulire que nous cherchons sous la
forme : = a z / + b .
1
1
Par identification, on obtient : a = / 0 et b =
/0 2
PCH
2 PCH
Lintgrale a donc pour solution finale :
/0
1
2 PCH z /
/0 /
/2
=V =Ce
+
1 2 PCH
z
2
/0
2 PCH

( )

( )

E.H.T.P.
Zorkani Mohammed
7-43
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Pour trouver la constante dintgration C, nous ferons appel aux conditions

initiales dabord avant deffectuer la manuvre, le rgime tait tablit, et
dv
lon avait :
=0
dt
0
Lquation de Bernoulli nous donne : z = PCH
Immdiatement aprs la fermeture, la cote du plan deau na pas chang
/0
/0
do z / (t = 0 ) = PCH
v 0 / 2 = PCH
.
Dautre part, immdiatement aprs la fermeture, le dbit dans le canal
damen na pas chang, do :
Q 0 = sv (t = 0 ) = SU 0

V / (t = 0 ) =

V
=1
U0

/0
En tenant compte de z / (t = 0 ) = PCH
et que V / (t = 0 ) = 1 lquation prend la
forme :
/0
/0

2 PCH PCH
+ z /
1
/
0
/2
/

1 2 PCH z e
V =
2
/0

2 PCH

/
Cette quation nous permet de dterminer la cote z max = z max
z c d la

( )

/
plus grande dnivellation ; la valeur de z max
sera obtenue en crivons que
V = 0 . Soit il faut rsoudre lquation :
/0
/0

/
+ z max
2 PCH PCH

/0 /
=0

1 2 PCH z max e
Rsolution de lquation du maximum :
On peut la rsoudre numriquement par approximation successive. On
peut galement faire appel une formule approche, dont voici la
justification :
/0 /
/0
/0
Log 1 2 PCH
z = 2 PCH
PCH
+ z/

/0 /
Si lon a : 2 PCH
z max 1 on peut dvelopper le Log :

/0
/0
/0
/
/
/
4 PCH
z max
8 PCH
z max
16 PCH
z max
/0 /
/0 /
Log 1 2 PCH z = 2 PCH z max

L
2
3
4
Notre quation se simplifie ainsi
2
4
4 /0 / 3
/0 2 /
/
1 z max
PCH
z max 2PCH
z max LL = 0
3

E.H.T.P.
Zorkani Mohammed
7-44
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

/0
/
Pour PCH
= 0 , on retrouve la solution z max
=1 .
Ainsi, il apparat naturel de rechercher une solution de la forme :

/0
/0
/
z max
= 1 + a PCH
+ b PCH
+ LL
2
1
En se limitant au 2me ordre, on obtient : a = et b =
3
9
/

Do la formule approche due D. E ydoux :

z
2 /0 1 /0 2
SI
/0 /
/
max

2 PCH
z max 1
PCH + PCH
z max
1

9
3
z

( )

/0
Constatation : si PCH
0,7 cette formule scarte peu de la solution prcise
(tableau) , qui est obtenue numriquement, soit :
/0
0
0,1
0,2
0,4
0,6
0,8
0,9
1
PCH
=
/
1
0,933 0,875
0,76
0,65
0,555
0,512 0,475
z max
=
D. Eydoux Directeur des Etudes lcole Polytechnique de Paris vers
les annes 20 du sicle dernier [Ouvrages : Les mouvements de leau et
les coups de blier dans les chemines dquilibre {Toulouse : 1119} ,
Hydraulique gnrale et applique {Paris 1921}, Hydraulique industrielle :
Paris : 1921, Turbines Hydrauliques 1225 en collaboration avec A. Rateau
et M. Gabriel Etc.].
4 6 2 ) Examinons ce qui se passe durant la 1re descente
du plan deau :
Lquation de Bernoulli est
d V/
/0 /2
V/
+ z / PCH
V =0
/
dz
/0
/0
Cette quation est identique la prcdente si lon change PCH
en PCH
,
do
/0
1
2 PCH z /
/0 /
/2
V =Ke
+
1
+
2
P
z
CH
/0 2
2 PCH

( )

/
Les nouvelles conditions initiales sont : V / = 0 pour z / = z max
valeur trouve
dans la premire partie de ltude.
Par un calcul identique celui du $4 6 1, on trouvera :
/0 /
2 PCH z max
1
/0 /
1 + 2 PCH z max e
K=
/0 2
2 PCH

( )

E.H.T.P.
Zorkani Mohammed
7-45
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

/2

( )

/0
2 PCH

/0
/

2 PCH zmax
z/
/0 /
/0 /
1
2
P
z
1
2
P
z
e
+

CH
CH max

) (

La cote minimale atteinte par le plan deau sera obtenue pour V / = 0 ainsi
lquation annulant V / est :
/0
/
/
zmin1
2 PCH zmax1
/0 /
/0 /
1 + 2 PCH z min1 1 + 2 PCH z max1 e
=0
L encore on peut rsoudre cette quation numriquement par
/0
/
/
approximation successive pour dterminer : z min1
= f PCH
, z max1
.

) (

z max1

z max2

t
/0
/
z min1
1 + 0,13 PCH

selon CALAME et GADEN

z min1

variation du niveau du plan d' eau dans la chemine d' quilibre

4 6 3 ) Etude des oscillations suivantes :
Lors de la remonte du niveau de leau dans la chemine, lquation
gnrale donnant V / 2 reste valable, savoir :
/0
1
2 PCH z /
/0 /
/2
V =Ke
+
1 2 PCH
z
2
/0
2 PCH
La dtermination de la constante dintgration sobtient en crivant que
/
pour z / = z min1
on a : V / 2 = 0 ce qui donne :
/0
/

2 PCH z / zmin1
1
/0 /
/0 /
/2
V =
1
2
P
z
1
2
P
z
e

CH
CH min1
/0 2
2 PCH

( )

( )

) (

(1

/
z max2

) ( 1

/0
2 PCH

Lquation du 2me maximum est :

/0
2 PCH

/
z min1

/0
/
/
zmin1
2 PCH zmax2
e
=0

A partir de cette formule, lquation de ce maximum et de lquation du 1r

minimum on peut gnraliser ainsi :
E.H.T.P.
Zorkani Mohammed
7-46
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

/
Pour le ime minimum : z min
I

(1+

/0
2 PCH

/
z min
I

) ( 1 +

/0
2 PCH

/
z max
I

/
Pour le ( + 1)ime maximum : z max
I +1

(1

/0
2 PCH

/
z max
I +1

) ( 1

/0
2 PCH

/
z min
I

/0
/
/
zmin
2 PCH zmax
I
I =0
e

)e

/0
/
/
zmin
2 PCH zmax
I +1
I

)=0

A titre dexemple pour PCH = 0,7 :

= n de loscillation
/0

/
/
z max
z min

0,60
1
0,38
0,28
0,22
2
0,18
0,15
3
4 6 4 ) Manuvre douverture instantane en tenant compte
des pertes de charge :
Avant toute manuvre, on suppose que le rgime est bien statique.
A louverture instantane, le dbit Q 0 prend immdiatement naissance, et

/ U

d U/

.
est fourni par la chemine, le terme est
=
=
1
nul
U

U
d T/
0

Rappel : on a tablit en variables adimensionnelles que

PCH
t
/ z
/
=
=
, / =
,
P
z
CH

T0
z
z

variables
adimension
nelles

V/ = V , U/ = U , v/ = v

v0
U0
U0
/
1 d U/
avec = 1 si v 0
/ dV
/
/

V
+
+
z
+

P
=
0

CH
d z / 2 d T /
= 1 si v 0
v/ = U/ + V/
V
Ainsi nos quations se rduisent :

d V/
d z/

/0
+ z / + PCH
v/ = 0
2

v/ = U/ + V/ v/ = 1 + V/
d V/

=0
d z/
Par voie analytique, on ne sait pas rsoudre cette quation.
Soit :

V/

/0
+ z / + PCH
1 + V/

E.H.T.P.
Zorkani Mohammed
7-47
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Les mthodes sont soit numriques ou graphique.

/
Nanmoins le tableau le tableau ci dessous, consigne la cote z min1
, cote
z
/0
/
/
du premier minimum en fonction de PCH
: z min1
= min1
1 + 0,13 PCH

z
/0
PCH

0,1

0,2

0,4

0,6

0,8

0,9

/
z min1

1,0125

1.025

1,05

1,075

1,103

1,12

1,15

1.16

/
z min1
=

z min1
z

z = v 0

L s

g S

1.12

1.08

1.04

/0
PCH
1.00
0.00

0.20

0.40

0.60

0.80

1.00

/
/
Cote du premier z min1
minimum en fonction de PCH

On constate plus les pertes de charge sont importantes, plus la

dnivellation est importante, ce qui tait tout fait prvisible.
En rsum :
Pour dimensionner une chemine dquilibre, il faut considrer 2
/0
valeurs de PCH
:
/0
Soit PCH
: valeur sous estime, utilise pour calculer z max1 , lors
dune fermeture.
/0
Soit PCH
: valeur sur estime, utilise pour calculer z min1 , lors
dune ouverture.

E.H.T.P.
Zorkani Mohammed
7-48
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

4 7 ) Chemine dquilibre tranglement [ rtrcissement ] :

Etranglement d' une chemine

Les surpressions les plus importantes ont lieu lextrmit aval de la

conduite et ont pour valeur : p = g z .
Leffort de freinage, fonction de p natteint une valeur maximale qu la fin
de la monte du plan deau dans la chemine. Le freinage serait plus
efficace si, ds le dbut de la manuvre de fermeture p pouvait avoir une
valeur notable. Cest la raison pour laquelle la base de la chemine, on
dispose un tranglement C. Soit R V la perte de charge provoque par
ltranglement, lorsque la vitesse dans la chemine est V .
La surpression a pour expression : p = g (z + R V ) . Ainsi, comme R V
est fonction de V 2 , on peut prvoir le fonctionnement de la chemine sur
les graphiques suivant :

p
g

p
= z (t ) + R V (t )
g

0
PCH

z (t )

R V (t )
effet optimal

=C

te

Effet de freinage

Lors dune fermeture progressive, R V est importante au dpart, puis

devient faible lorsque la cote z est grande.

E.H.T.P.
Zorkani Mohammed
7-49
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Remarquons enfin, que ltranglement nintroduit pratiquement pas de perte

de charge en rgime permanent

0 : Lim R V (t ) 0 .
t
t
4 7 1 ) Equations gnrales :

P = 1 si v 0

L
d
v
avec

P = 1 si v 0
+
z
+

P
+

R
=
0

P
CH
R
V
= 1 si v 0
g dt
R

R = 1 si v 0

s v = S (V + U)
De mme que nous avons pos [$(4 3)] :
P
/0 /2
/
/
PCH
= CH
PCH = PCH
z = PCH
z V

z
R
/0 /2
/
/
Par analogie nous posons : R V
= V R V = R V
z = R V
z V
z
Les quations dans dimension deviennent :
/
1 d U/
/0
/0
/
/2
/ dV
V
+
+
z
+

V
+

R
V/ 2 = 0

P
R
V
CH
/
/
2 d T
dz

v/ = U/ + V/
4 7 2 ) Rpartition des surpressions le long de la conduite L :
Le calcul est identique celui du $ 4 2.
En crivant le thorme de Bernoulli entre le point A et un point B la cote
l de la conduite (voir figure au dessous), la perte de charge en C valant
l

l
v2
1 H A
R V , nous obtenons : p = g (z + R V )
2g L
L

origine de z
air

B
V
C

Rservoir

A
Q
E.H.T.P.
Zorkani Mohammed
7-50
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

l
v2
1 H A
Car on a tablit que : p = g (z + H C )
2g L
L

l
Cest donc pour = 1(cest dire en C) que la surpression est la plus
L
v2
importante. savoir : au terme de
: p = g (z + R V )
2g
4 7 2 ) Etude dune fermeture instantane :
Proposons-nous simplement de rechercher sil est possible dobtenir
te

p = C pendant la priode de monte deau.

Lquation simplifie (U = 0 ) est donc :
d V/

V/

/0
/0
+ z / + PCH
+ RV
V/ 2 = 0

dz
Cette quation est identique celle obtenue lors de la fermeture
instantane avec une chemine dquilibre sans rtrcissement la base,
savoir :
/
/0
/
/ dV
V
+
z
+
P
V/ 2 = 0
CH
/
dz
Soit :
/0
/0
+R V z/
2 PCH
1
/0
/0

V/ 2 = K e
+
1 2 PCH
+ RV
z/
2
/0
/0
2 PCH
+ RV

( )

) ]

/0
A linstant initial, z / = PCH
et V / = 1 ; ce qui impose :

K=

/0
/0
2 PCH
+ RV

et par suite
2

V/ =
2

/0
PCH

/0
/0
/0
1 2RV
PCH
+ RV

/0 2
+ RV

) e 2P

/0 /0
/0
P +R V

CH CH

) [

)]

/0
/0
/0 /0
/0
/
1 2 PCH + R V z 1 2R V PCH + R V

/0
/ 0 / 0
2 PCH
+ R V PCH + z /

E.H.T.P.
Zorkani Mohammed
7-51
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

/ 0 /2
La surpression : p = g (z + R V ) = g z z / + R V
V do :

gz
/0
/0
/0
/0
/
p =
R V + 2 PCH PCH + R V z
2
/0
/0

+ RV
2 PCH

)]

/0
/0 /0
/0
e
RV
1 2R V
PCH + R V

/0
/ 0 / 0
2 PCH
+ R V PCH + z /

Cas idal : il est obtenu que si le freinage est constant.

Cette expression ne peut tre constante quel que soit z / que si
/0
PCH
=0
/0
/0
/0
2 PCH PCH + R V = 0
1
/0
/0 /0
/0
1 2R V PCH + R V = 0 R V =
2

1
Ainsi dans le cas idal on a :
p idal =
g z
2
Ce cas idal ne sera donc obtenu que si les pertes de charge dans la
conduite sont ngligeables, et que si ltranglement a t correctement
1
/0
construit afin que R V
=
.
2
Comparaisons des chemines avec et sans tranglement :
Supposons ces conditions ralises et comparons les sections S1 et S2 des
2 chemines dquilibre, le premier tant une chemine ordinaire alors que
1
/0
le second est tranglement [rtrcissement] avec R V
=
et qui, toutes
2
les deux, donnent la mme surpression p M dans le canal damen.
Pour le premier : chemin sans rtrcissement

p M =

g z1

avec

z1

= v0

L s
L g

on en dduit S1 = v 02 s
g S1
g p M

Pour le second : chemin avec rtrcissement

2

v 02 L g
S
L s
1

= 1
on en dduit S2 =
s
p M =
g z 2 avec z 2 = v 0
g S2
2 g p M
2
2
Do lintrt indniable (conomie de 50% sur le volume) des chemines
munies dun rtrcissement [tranglement].
E.H.T.P.
Zorkani Mohammed
7-52
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Mthodes de Protection
(5) Rservoir air
N.B. : Dimensionnement dune conduite & dune chemine
Le coup de blier en dehors de rupture spectaculaire de conduites et de
destructions dappareils de pompage peut avoir pour consquence des
destructions de joints, des dbotements de conduites, causes des pertes
deau importantes, et des dtriorations de robinetterie ou dappareils de
comptage. On a tabli que la surpression [une conduite] est donne par :
Q

0
u 0 = s le courant uniforme dans la conduite [rgime permanent normal]

K = [ ] module de compressibilit d' eau

E module d' lasticitde' la matre' paroi

a u 0
1420(m/s )
h=
la clrit ' onde':
: a
D diamtre de la conduite
g
K D
1 +

E e
e paisseur de la conduite

g = 9,81m/s 2 champ de pesanteur

: 100
acier
fonte
: 37 50

K
fonte centrifuge : 41

On donne : = pour tube en amiate ciment : 12

a=

K
K D
plomb

1 +
:9
E e

bton
:
10

bois

:5
DATA : calcul des clrits donde de pression
Pour une galerie creuse dans une roche seine on calcule la clrit des
ondes de pression par :

2b

g
a=

2
+
E2

tel que =

roche

g
Galerie dans roche saine

E.H.T.P.
Zorkani Mohammed
7-53
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

O E 2 est le module dlasticit du Rocher, et = 0,001g/cm 3 sont

respectivement le module dlasticit et le poids spcifique de leau.
Alors que pour une galerie blinde on utilise la formule :

b2
Ee

g
a=

2b
(1 3 )
+
Ee
1

avec 3 =

b2
+
.e

(c

b2
2cE1
2

(m 2 + 1)b
m 2 2

rocher
c
b
E1 : le module dlasticit du bton
E 2 : le module dlasticit de la roche
bton
E : le module dlasticit de la tle dacier
tle
e : lpaisseur de la tle dacier
Galerie blinde
c : le rayon extrieur du revtement en bton
b : le rayon intrieur du revtement en bton entre la tle et le rocher
a u0
permet en premire approximation
g
de choisir (le matriau) et de dimensionner une conduite contre le coup de
blier dfaut de calcul (graphique ou numrique) prcis. On signale que
si H est la pression existante dans la conduite en un point de celle ci
avant le coup de blier, la pression au mme point prend la valeur : H h .
Pour la protection dune conduite damene on a prsent ltude dune
chemine dquilibre et son dimensionnement contre le coup de blier.
[ Cest Joukowsky qui a le premier publi les quations fondamentales du
coup de blier en 1898 (en russe) {Nanmoins La thorie du coup de
blier reste attach au nom dAllievi depuis 1913}, ce travail a t traduit
en anglais en 1904 par O. Simin dans Proc. Ann. Water Works Assoc. ].

La formule de Joukowsky : h =

Etude du Rservoir dAir :

Dans le cas o il sagit dentranement par un moteur lectrique, o larrt
de fourniture de courant peut tre inopin et o les parties en rotation nont
quune inertie faible, il est indispensable de prvoir des dispositifs de
protection.
Le dispositif le plus courant est celui du rservoir dair (voir figures) . Lors
de larrt des groupes lectriques, lair du rservoir se dtend et refoule un
E.H.T.P.
Zorkani Mohammed
7-54
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

dbit deau dans la conduite, le ddit qui se substitue celui des pompes.
Au retour de londe de pression (phase de surpression), lair du rservoir se
comprime et leau pntre dans le rservoir au lieu de scraser sur
lextrmit de la conduite avec tous les effets nfastes quil sagit dviter.

eau

air
rservoir d' air
eau
By - passe

clapet

pompe
Lamortissement du processus est amlior en interposant entre le
rservoir dair anti blier et la conduite un clapet trou qui laisse passer
librement leau dans le sens anti blier conduite (phase de dpression)
et qui ralentit, en le laliment, le courant deau en sens inverse (phase de
compression). Il est utile galement de munir le dispositif dun by passe
permettant dapporter lusage quelques corrections de finition par
manuvre de la vanne.
Il existe une formule permettant le calcul du volume du rservoir dair :

Hd Hd
Hd
v 02 Vc
1

=
log
H0 H0
H 0 2gH 0 Vr
O

H0
Hd
Vc
Vr
v0

: pression initiale (en mtre deau) de lair du rservoir.

: pression en fin de dtente.
: volume de la conduite.
: volume du rservoir.
: vitesse de leau dans la conduite en rgime normal.

E.H.T.P.
Zorkani Mohammed
7-55
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Quand les hauteurs de refoulement sont relativement faibles, le rservoir

dair peut tre avantageusement remplac par une chemine dquilibre ;
cette disposition est reprsente la figure en dessous :

chemine

pompe

N.B. : mode de fonctionnement thermodynamique dun rservoir dair

Soit un rservoir air branch en un lieu M dune conduite (voir figure
suivante) dont le volume dair est V0 la pression P0 du rgime permanent
initial ; si H 0 est la hauteur deau correspondante cette pression, celle ci
est donne par : P0 = g H 0 + H a o H a est, en hauteur deau, la pression

atmosphrique.
air

a / / S/

Q moyen =

Qt

A
Q /t

Q /t = Q t + Q //t

Q //t

Q t + Q t + dt
2

a // // S //

En rgime vari ce volume V0 se dtendra et se comprimera selon

respectivement la descente [ Q t 0 ] ou la monte du plan deau dans la
chemine [ Q t 0 ], et lon aura entre la pression absolue P = g H + H a et

le volume V la relation :
P0 V0 = P V si le phnomne (compression dtente ) est isotherme

Ou P0 V0 = P V sil est adiabatique. Lexposant , qui est le rapport

des chaleurs spcifiques volume et pression constants est gal 1,4
pour lair. Si le rgime est lent cest la relation isotherme qui est valable. Si
elle est trs rapide cest lquation adiabatique qui est valable. Dans la
pratique cest une loi intermdiaire entre les 2 qui est acceptable.
Par exemple choisissons lquation isotherme, et si le volume du rservoir
dair est Vt linstant t on aura : Vt (H t + a ) = V0 (H 0 + a )
E.H.T.P.
Zorkani Mohammed
7-56
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Alors qu linstant (t + dt) on aura :

Vt + dt (H t + dt + a ) = V0 (H 0 + a )
Do lon tire : Vt + dt (H t + dt + a ) = Vt (H t + a ) o Vt + dt = Vt Q t dt
V (H + a ) Vt + dt a
V
H t + dt = t t
= (H t + a ) t a
Vt + dt
Vt + dt
Vt
(H t + a ) a = ( (H 0 + ) a ) a
H t + dt =
H0 + a
Vt Q t dt
Q dt
(H t + a ) t
Or H t est connu (tude du rgime transitoire), on trouve ainsi H t + dt et on
V (H + H a )
en dduit alors Vt + dt = t t
. Do lon calcule un dbit moyen sur
H t + dt + H a
V
Vt Q t + Q t + dt
lintervalle [t , t + dt ] soit Q moyen = t + dt

do lon dduit
dt
2
alors Q t + dt connaissant Q t & ainsi de proche en proche
H + Ha
dV 2V (H + H a ) d H t
N.B. :
Vt = V0 0
Qt = t = 0 0
H t + Ha
dt
(H t + H a )2 d t

Une autre approche :

Le rservoir air est une capacit qui se remplit ou se vide. La fonction qui
lie le volume dair V la pression H est : Vt (H t + a ) = V0 (H 0 + a ) alors
H + Ha
V = Vt = V0 0
= (H )
H + Ha
On pose souvent : z H + H a z 0 H 0 + H a .
Quand la pression H en M augmente le volume dair dans le rservoir
diminu, et vis versa, car la drive de / (H ) est toujours ngative :
H0 + Ha
d d V
/ (H ) =
=
= V0
0
dH dH
(H + H a )2
On peut crire :

V(H + dH ) V(H ) +

H0 + Ha
d
dH = V(H ) V0
dH
dH
(H + H a )2

Ainsi en remplaant V(H ) par son expression on obtient :

H + Ha
V(H + dH ) V0 0
H + Ha

dH
1
H + Ha

E.H.T.P.
Zorkani Mohammed
7-57
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

remarque :
Dans le cas adiabatique, on aurait :
1

(H ) = V0

(H 0 + H a )
1

(H + H a )

et / (H ) =

(H + H a )
0

V0

(H + H a )

+1

avec (air ) 1,41

FONDEMENT DE LA
THEORIE DES CONDUITES SOUS PRESSION
AVEC RESERVOIR DAIR
Nous devons Rateau et Camichel les premiers travaux dans ce
domaine. Une tude dtaille de A. Foch [contribution ltude des coups
de blier dans les conduites munies dun rservoir dair, Toulouse 1920] et
un mmoire dAllievi doivent tre nots tout particulirement.

r
v 0 si v
r r
v 0 si v//g

air

r
g

S v
S2

S1

v2

v1

Soit v la vitesse dbitante de leau qui entre dans le rservoir dair, v1 et v 2

les vitesses de leau respectivement dans les conduites darrive et de
dpart. On a conservation du dbit volumique (incompressible) :
S v = S1 v1 S2 v 2
[ Le plus souvent on a S = S1 = S2 de sorte quon a : v = v1 v 2 ]. Soient :
P
V Le volume et P La pression de lair dans le rservoir [
= z = H + Ha ]
g
r
v 0 si v
O H a est la pression atmosphrique, on a : dV = S vdt :
r
v 0 si v
De plus on a :

z V = [H + H a ] V = C

Soit z V = [H + H a ] V = C

te

te

(compression isotherme)

(compression adiabatique avec = 1,41 ).

1

Foch crit dune manire gnrale :

V = avec 1 1,41 .

E.H.T.P.
Zorkani Mohammed
7-58
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

Il en rsulte :
1

1
V d z
d

v1 v 2 = z = 0
S dt
S d t z 0

Ou pour petites surpressions :

V0 d z
V0 d H
v1 v 2 =
=
S z0 d t S z0 d t

[A] Thorie de Foch pour grands rservoirs dair

(sans tenir compte de llasticit de leau)
Foch considre des conduites courtes avec des rservoirs dair
relativement grands. Supposons que le rservoir dair se trouve
lextrmit infrieure de la conduite et cre une pression z cet endroit. A
lextrmit suprieure de la conduite rgne la pression z 0 .
On peut crire pour cette conduite : ( z 0 = H 0 + a )
dv
L + g(z z 0 ) = 0
dt
z
Do en utilisant [ dV = S vdt ] et en posant
= , aprs intgration par
z0

rapport au temps [ 0t dV = Vt V0 = 0t Svdt Vt = V0 S0t vdt ] { si 1 }

on obtient (il ne faut pas oublier de tenir compte du signe de v) :
2

LV0 d
1 1
te
+

= C

2z 0Sg d t
1
Cette quation diffrentielle convient la discussion et la comprhension
thorique du phnomne physique malgr quelle prsente des difficults
mathmatique (une tude numrique est ncessaire).

E.H.T.P.
Zorkani Mohammed
7-59
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent

[B] Thorie dAllievi pour rservoirs dair lastiques {petits}

[ L. Allievi, Comere daria nelle turbazionni prementi, Elettrotecnica,
vol. 23, n0 20, 25 octobre 1936 ]
Par contre, si le rservoir dair nest pas assez grand pour que la
compressibilit de leau puisse tre nglige. Il faut de nouveau faire appel
aux quations fondamentales du problme du coup de blier :
a
H t - t + t 2H 0 = ( v t - t m v t )
g
[le signe suprieur sapplique aux turbines et le signe infrieur aux pompes]
Allievi considre le cas dune pompe situe au pied dune conduite avec
rservoir dair. Au lieu de lquation [ dV = Svdt ] il crit son intgrale au
temps t :
r
le rservoir d' air
t
v 0 si v

se remplie
Vt = V0 S v dt avec :
r
le rservoir d' air
se vide
0
v 0 si v
Lquation isotherme : z 0 V0 = (H + H a ) Vt = z t Vt apparat ici comme
auxiliaire. Si on pose de plus :
av 0
av 0
=
=
2g (H 0 + a ) 2g z 0
Allievi crit pour les pompes :
v0
(z t + 2z t + t + z t + 2t 4z 0 )
v t + v t + 2 t =
2 z 0
Nous poserons ensuite avec Allievi :
v t - t + v t
-1
t
[] = T s -1
t - t vdt
2
Pour une fermeture brusque on a pendant la premire phase :

2V0 z 0
2V z

1 v 0 = 0 0 1
v 0 + v t =
Sv 0 z t
S z t

Pour les phases ultrieures nous crivons :

1
2V0
2
1

v t - t + v t + t =
z 0
+
S
z
z
z
t
t + t
t - t
Cette quation rsout le problme.

( )

E.H.T.P.
Zorkani Mohammed
7-60
Dpartement dHydraulique
CHAPITRE 7 : ETUDE DU COUP DE BELIER coulements en charge non permanent