RTH SUR UNE CHEMINE D’EQUILIBRE
RANDRIAMIRADO Misaela
INTRODUCTION
En généralement, les régimes variables ou instationnaires lorsque les grandeurs hydraulique (Pression
P, Hauteur piézométrique H, Débit Q et Vitesse V) varient non seulement dans l’espace mais aussi dans
le temps. Le Régime Transitoire Hydraulique c’est une méthode de simulation dans les domaines
hydrauliques par exemples : centrales hydroélectrique, réseaux de distribution d’eau, réseaux
hydraulique industrielle, rivières et canaux à surfaces libres, sante : circulation sanguine artérielles etc…
Ce devoir consiste à faire l’étude de la centrales hydroélectrique, ce système qu’on étude l’installation
à chemine d’équilibre alors elle a trois section (galerie, chemine d’équilibre et conduite force) pour
observer et apprécier la réaction de l’ensemble de l’installation, déterminer les zones à risque et calculer
les sollicitations maximales pour dimensionner correctement, prévoir et dimensionner des dispositifs de
protection.
Ce travail est divisé en quatre parties : la première partie va parler la description du système physique,
le second parlera mise en équation, le troisième va être dédier la résolution numérique des équations et
pour finir la quatrième partie qui parlera le résultat et discussion.
1
I. Description du système physique
I.1 Schéma du système physique
Figure I-1 : système physique
I.2 Les données
• Pour la galerie
Tableau I-1 : donnée de galerie
Dimension Valeur Unité
LG 100 m
DG 0.75 m
fG 0.03 _-
AG 1000 m/s
HD 150 m
CDSVo 0.009 m2
nG 20 -
• Pour la cheminée
Tableau I-2 : donner de la cheminée
Dimension Valeur Unité
DE 3 m
2
• Pour la conduite forcée
Tableau I-3 : donnée de la conduite force
Dimension Valeur Unité
LF 600 m
DF 0.5 m
fF 0.018 -
AF 1200 m/s
nF 10 -
II. Présentation des équations
II.1 Régime permanent initial
Pour calculée Vi et hi
On applique Bernoulli entre (A) surface libre dans le barrage et (B) dans le jet libre issu de vanne.
PA V P V
+z A + A -Δh G -Δh F = B +z B + B (1)
γ 2g γ 2g
VB2
On a : H D -Δh G -Δh F = avec (2)
2g
2
f L VG
Δh G = G G × 0 et (3)
DG 2g
2
f L VF
Δh F = F F × 0 (4)
DF 2g
VGo et VFo sont les vitesses au régime permanent initial avec Qo sont le débit
Qo Q
Avec VGo = et VFo = o (5)
SG SF
2gHD
Donc Qo = (6)
1 f L f L
2
+ G G2 + F F2
(CdSv) o DGSG D FSF
A t = 0 pour i=1 à n1
2
Q f (i-1)Δx G VGo
Vi = o et h i =H D - G × (7)
SG DG 2g
Pour i=n2 à n3
Qo f (i-i )Δx F Vi 2
Vi = et h i =h i1 - F 2 × (8)
SF DF 2g
3
Avec :
LG L
Δx G = et Δx G = G (9)
nG nG
II.2 Régime variable
On applique les formules établies pour la méthode des caractéristiques :
II.2.1 A l’entrée
Pour i=1à n1, hi=HD
h1 -y N
V1 = avec y N =h -2 -BG V2- +R G V2- V2- (10)
BG
II.2.2 Sur le nœud interne de la galerie
Pour i=2 à n1-1,
y M +y N yM +h i
hi = et Vi = (11)
2 BG
- - - -
Avec yM =h i-1 +BG Vi-1 +R G Vi-1 Vi-1 et (12)
- - - -
y N =h i+1 +BG Vi+1 +R G Vi+1 Vi-1 (13)
II.2.3 A la base de la cheminée
La continuité :
QG -QF =QE (14)
SG Vi1 -SFVi2 =SE ×VE (15)
VE vitesse de déplacement de la surface libre dans la cheminée
y E -y -E
VE = avec yE = hn1 ou yE = hn2 (16)
Δt
Avec :
h n1 =yM -BG Vn1 ; (17)
h n2 =yN -BFVn2 ; (18)
hn2 = hn1 (19)
yM -h n1 h n -y N
SG -SE 2 =SE VE , (20)
BG BF
yM -h n1 h n2 -y N h n1 -h n1 -
on a : Vn1 = ; Vn 2 = et VE = (21)
BG BF BF
4
Posons :
SG S S
=CG , G =CG , F =C E (22)
BG BG Δt
SG S
(yM -h n1 )- F (h n1 -y N )=CF (h n1 -h n1 - ) (24)
BG BF
CG y M +CF y N +CE h n1 -
h n 2 =h n1 = (25)
CG +CF +CE
II.2.4 Sur le nœud interne du conduit forcée
Pour i=n1+1 à n3-1
y M +y N y -h
hi = ; Vn 2 = M i (26)
2 BF
II.2.5 A la sortie (juste avant la vanne)
On à la même condition aux limites que pour le cas simple traitée précédemment.
On va résoudre l’équation en second dégrée
C1Vn32 +C2 Vn3 +C3 =0 (27)
1 S2
Avec : C1 = 2
-1 , C2 = BF , C3 =y M (28)
2g (CdSv)
On prend le racine positive
-C2 + Δ
Vn3 = (29)
2C1
Lorsque la vanne est complètement fermée t t F
Temps de fermeture Vn3 = 0
Dans tous les cas :
h n3 =yM -BFVn3 (30)
yM =h -n3 -1 -BF Vn-3 -1 +R F Vn-3 -1 Vn-3 -1 (31)
t 1.5
CdSv=(CdSv)o (1- ) (32)
tf
5
III. Organigramme pour la procédure calcule
Données :
H D ; ( CDSV )0 ;t max ;A G ;A F ;g;
DG ;D F ;D E ;L G ;L F ;f G ;f F ;n G ;n F ;
Constantes :
πDG 2 πD F 2 L L f Δx
SG = ;SF = ;Δx G = G ;Δx F = F ; R G = G G
4 4 nG nF 2gDG
f FΔx F A A Δx
;R F = ; BG = G ;BF = F ; Δt= G ;
2gD F g g AG
f G .LG f .L 1
KG = ;K F = F F2 ;K V =
( CDSV )
2 2
DG .SG D F .SF
0
2gH D S S
Q0 = ;CG = G ;CF = F ;
K G +K F +K V BG BF
t=0;k=1;tg=0
Régime permanant initial t f =0;t=max(tg,t f )
Δt=Δt f
Pour la galerie :
i = 1 à n1
Q0
V (i ) =
SG
f ( i-1) Δx G .V 2 (i)
h ( i ) =H D .
2gDG
Chemine :
h (n2) =h (n1)
V(n2) =V(n1)
6
1
Pour le conduit force
i = n2à n 3
Q0
V (i ) =
SF
f ( i-n 2 ) .Δx F .V 2 (i)
h ( i ) =h(n1 ).
2gD F
Régime variable :
t<tmax ;
k = 1+k ;
i = 1 à n1
h i- = h i
Vi- =Vi
i = n2à n 3
h i- = h i
Vi- =Vi
Entre de la galerie :
2 4
7
2
h1 =H D
y N =h -2 -BV2- +RV2- V2-
h1 -YN
V1 =
BG
Section interne dans la galerie :
i = 2 à n1−1
y M =h (-i-1) +BG V(-i-1) -R G V(-i-1) V(-i-1)
y N =h (-i+1) -BG V(-i+1) +R G V(-i+1) V(-i+1)
y M +y N y -y
hi = ; Vi = M N ;
2 2BG
A la base de la cheminée :
- - - -
y M =h (n1-1) +BG V(n1-1) -R G V(n1-1) V(n1-1)
- - - -
y N =h (n2+1) -BF V(n2+1) +R F V(n2+1) V(n2+1)
h n1 =y M -BG Vn1;h n2 =y N +BF Vn2
Δt(CE y M +C F y N )+SE h - n1
h n1 =h n2 =
Δt(CG +CF )+SF
y M -h n1 h -y
Vn1 = ;Vn2 = n2 N
BG BF
Section interne dans la conduite forcée :
8
3
4
i = n2+1 à n 3−1
y M =h (-i-1) +BF V(-i-1) -R FV(-i-1) V(-i-1)
y N =h (-i+1) -BF V(-i+1) +R F V(-i+1) V(-i+1)
y M +y N y -y
hi = ; Vi = M N ;
2 2BF
A la sortie de la conduite :
yM =h -(n3-1) +BF V(n3-1)
- -
-R F V(n3-1) -
V(n3-1)
Fermeture de la vanne :
t tF
1.5
t
CDSV = ( CDSV )0 . 1- CDSV =0;Vn3 =0
tf
C1Vn32 +C2 Vn2 +C3 =0
S 2 1
C1 = F -1 . ,C2 =B,C3 =-y M
CDSV 2g
- - - -
y M =h (n3-1) +BF V(n3-1) -R F V(n3-1) V(n3-1)
Δ=C2 2 - ( 4.C1.C3 )
-C2 Δ
Vn3 = ;h n3 =y M -BF Vn3
2.C1
Enregistrement :
tps ( k ) =t;Pn3 ( k ) =Pn3 ;Vn2 ( k ) =Vn2 ;Vn3 ( k ) =Vn3 ;V1 ( k ) =V1;
Vn1 ( k ) =Vn1;;y E ( k ) =y E ;
9
Après cet organigramme, on va simuler dans le logiciel MATLAB.
IV. Représenter graphiquement et analyse de résultats
Le voilà les résultats des grandeurs hydraulique (Pression et Vitesse) qu’on a trouvé sur le MATLAB.
IV.1 Pression dans la section située juste avant la vanne pour t = 0 à 10 s
Figure IV-1 : courbe de pression juste avant la vanne pour t entre 0 à 10s avec et sans cheminée
✓ Interprétation :
D’après le résultat que nous avons vu sur la figure ci-dessus, l’allures de la pression à la sortie est variée
jusqu’au temps maximale de tmax=10s. car il y a un phénomène d’oscillatoire dans un conduite ça veut
dire l’oscillation en masse c’est un mouvement d’ensemble d’oscillatoire amortie de la masse d’un fluide
contenue dans la conduite. Donc si en ferme le vanne, l’eau retour fait le mouvement d’oscillation
pendule.
10
IV.2 Vitesse à l'entrée et à la sortie de la conduite forcée pour t = 0 à 10 s
Figure IV-2 : vitesse à l’entrée et à la sortie de la conduite forcée pour l’intervalle de temps de 0 à
10s
✓ Interprétation :
D’après le résultat que nous avons vu sur la figure ci-dessus, les deux vitesses ne sont pas des mêmes
allures. La variation du vitesse Vn2 à la sortie de la conduite forcée n’est pas constante suivant la
hauteur (Ye) jusqu’au temps maximale de tmax=10s. parce que à temps du zéro jusqu’à 10s la vitesse
fait le régime transitoire jusqu’au final et puis la vitesse Vn3 à l’entrée de la vanne est constante
jusqu’au temps maximale de tmax=10s.
11
IV.3 Vitesse du niveau libre dans la cheminée pour t = 0 à 600 s
Figure IV-3 : vitesse à l’entrée et à la sortie de la galerie pour l’intervalle de temps de 0 à 600s
Sur le résultat que on a trouvé au-dessus, les deux vitesses (V1 (bleu), Vn1 (vert)) sont des mêmes
allures, jusqu’au temps maximale de tmax=600s, car le chemine d’équilibre a équilibré les deux vitesses
jusqu’au temps final.
12
IV.4 Niveau d’eau dans le chemine de l’équilibre.
Figure IV-4 : chemine d’équilibre
D’après le résultat que nous avons vu sur la figure ci-dessus, l’eau fait le mouvement d’oscillation
jusqu’à qu’elle est arrivé au niveau statique.
13
Conclusion
Pour conclure, l’objectif de ce devoir et aussi le plus important doit être savoir dans ce cours RTH les
grandeurs hydraulique (pression P ; hauteur piézométrique H ; débit Q et vitesse V) entant que
l’hydraulicien. Pour faire la simulation du système on a utilisé la logiciel MATLAB pour savoir l’état
final, le but de cette simulation pour éviter le risque des problèmes après l’installation du système.
D’après le résultat que nous avons trouvé sur les figures ne sont pas des mêmes allures, ni des même
temps fermeture. Car pour faire la simulation les deux temps différents tmax = 10s et tmax =600s pour
avoir quelque résultat à bien réalisé le projet et pour améliore le rendement dans la centrale
hydroélectricité.
Cette devoir qui apprend toute l’étude théorique et la simulation réel sur le MATLAB, donc
l’opportunité pour faire de ce cours sont pour m’aider quand fait l’installation Hydraulique et ce cours
est l’outil de base de l’addiction d’eau dans le domaine hydraulique.
14
Annexe
%Données:
%Galérie
Lg=1000;Dg=0.75;
Ag=1000;fg=0.03;ng=20;
%Conduite forcé
Lf=600;Df=0.5;
Af=1200;ff=0.018;nf=10;
%Cheminée d'équilibre
De=3;
%Réservoir
Hd=150;
%Vanne
CdSvo=0.009;tmax=10;g=9.81;
%Constantes:
%Galérie d'ameneé
Sg=(pi*Dg^2)/4;Bg=Ag/g;
dxg=Lg/ng;Rg=(fg*dxg)/(2*g*Dg);
dtg=dxg/Ag;Cg=Sg/Bg;
%Conduite forcé
Sf=(pi*Df^2)/4;Bf=Af/g;
dxf=Lf/nf;Rf=(ff*dxf)/(2*g*Df);
dtf=dxf/Af;Cf=Sf/Bf;
%Cheminée d'équilibre
dt=dtf;Se=(pi*De^2)/4;Ce=Se/dt;
%Débit
Qo=sqrt((2*g*Hd)/((1/CdSvo^2)+((fg*Lg)/(Dg*Sg^2))+((ff*Lf)/(Df*Sf^2))));
tF=2,1;
i1=ng+1;
i2=ng+2;
i3=i2+nf;
z=0;
gamma=9180;
%initialisation
t=0;k=1;
tg=0;
tf=0;
t=max(tg,tf);
dt=dtf;
for i=1:i1
V(i)=Qo/Sg;
h(i)=Hd-fg*dxg*(i-1)*V(i)^2/(Dg*2*g);
end
%galérie
h(i2)=h(i1);
V(i2)=V(i1);
for i=i2:i3
V(i)=Qo/Sf;
h(i)=h(i1)-((ff*dxf*(i-i2)*(V(i))^2)/(Df*2*g));
end
%enregistrement RPI
k=k+1;
tps(k)=t;
Vi1(k)=V(i1);
Vi2(k)=V(i2);
Vi3(k)=V(i3);
A
ye(k)=h(i1);
Pi3(k)=(h(i3)-z)*gamma*10^-5;
V1(k)=V(1);
while (t<tmax)
k=k+1;
tg=tg+dtg;
tf=tf+dtf;
t=max(tf,tg)
for i=1:i3
hn(i)=h(i);
Vn(i)=V(i);
end
%Entrée
h(1)=Hd;
yN=hn(2)-Bg*Vn(2)+Rg*Vn(2)*abs(Vn(2));
V(1)=(h(1)-yN)/Bg;
%Noeuds interne dans le galerie 1 à 11
for i=2:i1-1
ym=hn(i-1)+Bg*Vn(i-1)-Rg*Vn(i-1)*abs(Vn(i-1));
yn=hn(i+1)-Bg*Vn(i+1)+Rg*Vn(i+1)*abs(Vn(i+1));
h(i)=(ym+yn)/2;
V(i)=(ym-h(i))/Bg;
end
%Base de la cheminée
ym=hn(i1-1)+Bg*Vn(i1-1)-Rg*Vn(i1-1)*abs(Vn(i1-1));
yn=hn(i2+1)-Bf*Vn(i2+1)+Rf*Vn(i2+1)*abs(Vn(i2+1));
h(i1)=(Cg*ym+Cf*yn+Ce*hn(i1))/(Cg+Cf+Ce);
h(i2)=h(i1);
V(i1)=(ym-h(i1))/Bg;
V(i2)=(h(i1)-yn)/Bf;
%Noeuds interne de la conduite forcée
for i=i2+1:i3-1
ym=hn(i-1)+Bf*Vn(i-1)-Rf*Vn(i-1)*abs(Vn(i-1));
yn=hn(i+1)-Bf*Vn(i+1)+Rf*Vn(i+1)*abs(Vn(i+1));
h(i)=(ym+yn)/2;
V(i)=(ym-h(i))/Bf;
end
ym=hn(i3-1)+Bf*Vn(i3-1)-Rf*Vn(i3-1)*abs(Vn(i3-1));
if(t<tF)
CdSv=CdSvo*(1-t/tF)^1.5;
c1=(1/(2*g))*((Sf^2/(CdSv)^2)-1);
c2=Bf;
c3=-ym;
delta=(c2)^2-4*c1*c3;
V(i3)=(-c2+sqrt(delta))/(2*c1);
else
CdSv=0;
V(i3)=0;
end
h(i3)=ym-Bf*V(i3);
%Enregistrement RPI
tps(k)=t;
Vi1(k)=V(i1);
Vi2(k)=V(i2);
Vi3(k)=V(i3);
ye(k)=h(i1);
Pi3(k)=(h(i3)-z)*gamma*10^-5;
V1(k)=V(1);
end
figure
plot(tps,Pi3,'r')
B
legend('Pi3')
xlabel('Temps t en (s)')
ylabel('Pression en (bar)')
grid on
figure
plot(tps,Vi2,'magenta',tps,Vi3,'cyan')
legend('Vi2','Vi3')
xlabel('Temps t en (s)')
ylabel('Vitesse en (m/s)')
grid on
hold on
%Données:
%Galérie
Lg=1000;Dg=0.75;
Ag=1000;fg=0.03;ng=20;
%Conduite forcé
Lf=600;Df=0.5;Af=1200;
ff=0.018;nf=10;
%Cheminée d'équilibre
De=3;
%Réservoir
Hd=150;
%Vanne
CdSvo=0.009;tmax=600;g=9.81;
%Constantes:
%Galérie d'ameneé
Sg=(pi*Dg^2)/4;Bg=Ag/g;
dxg=Lg/ng;Rg=(fg*dxg)/(2*g*Dg);
dtg=dxg/Ag;Cg=Sg/Bg;
%Conduite forcé
Sf=(pi*Df^2)/4;Bf=Af/g;
dxf=Lf/nf;Rf=(ff*dxf)/(2*g*Df);
dtf=dxf/Af;Cf=Sf/Bf;
%Cheminée d'équilibre
dt=dtf;Se=(pi*De^2)/4;Ce=Se/dt;
%Débit
Qo=sqrt((2*g*Hd)/((1/CdSvo^2)+((fg*Lg)/(Dg*Sg^2))+((ff*Lf)/(Df*Sf^2))));
tF=2,1;
i1=ng+1;
i2=ng+2;
i3=i2+nf;
z=0;
gamma=9180;
%initialisation
t=0;k=1;
tg=0;
tf=0;
t=max(tg,tf);
dt=dtf;
for i=1:i1
V(i)=Qo/Sg;
h(i)=Hd-fg*dxg*(i-1)*V(i)^2/(Dg*2*g);
end
%galérie
h(i2)=h(i1);
V(i2)=V(i1);
for i=i2:i3
V(i)=Qo/Sf;
h(i)=h(i1)-((ff*dxf*(i-i2)*(V(i))^2)/(Df*2*g));
end
%enregistrement RPI
C
k=k+1;
tps(k)=t;
Vi1(k)=V(i1);
Vi2(k)=V(i2);
Vi3(k)=V(i3);
ye(k)=h(i1);
Pi3(k)=(h(i3)-z)*gamma;
V1(k)=V(1);
while (t<tmax)
k=k+1;
tg=tg+dtg;
tf=tf+dtf;
t=max(tf,tg)
for i=1:i3
hn(i)=h(i);
Vn(i)=V(i);
end
%Entrée
h(1)=Hd;
yN=hn(2)-Bg*Vn(2)+Rg*Vn(2)*abs(Vn(2));
V(1)=(h(1)-yN)/Bg;
%Noeuds interne dans le galerie 1 à 11
for i=2:i1-1
ym=hn(i-1)+Bg*Vn(i-1)-Rg*Vn(i-1)*abs(Vn(i-1));
yn=hn(i+1)-Bg*Vn(i+1)+Rg*Vn(i+1)*abs(Vn(i+1));
h(i)=(ym+yn)/2;
V(i)=(ym-h(i))/Bg;
end
%Base de la cheminée
ym=hn(i1-1)+Bg*Vn(i1-1)-Rg*Vn(i1-1)*abs(Vn(i1-1));
yn=hn(i2+1)-Bf*Vn(i2+1)+Rf*Vn(i2+1)*abs(Vn(i2+1));
h(i1)=(Cg*ym+Cf*yn+Ce*hn(i1))/(Cg+Cf+Ce);
h(i2)=h(i1);
V(i1)=(ym-h(i1))/Bg;
V(i2)=(h(i1)-yn)/Bf;
%Noeuds interne de la conduite forcée
for i=i2+1:i3-1
ym=hn(i-1)+Bf*Vn(i-1)-Rf*Vn(i-1)*abs(Vn(i-1));
yn=hn(i+1)-Bf*Vn(i+1)+Rf*Vn(i+1)*abs(Vn(i+1));
h(i)=(ym+yn)/2;
V(i)=(ym-h(i))/Bf;
end
ym=hn(i3-1)+Bf*Vn(i3-1)-Rf*Vn(i3-1)*abs(Vn(i3-1));
if(t<tF)
CdSv=CdSvo*(1-t/tF)^1.5;
c1=(1/(2*g))*((Sf^2/(CdSv)^2)-1);
c2=Bf;
c3=-ym;
delta=(c2)^2-4*c1*c3;
V(i3)=(-c2+sqrt(delta))/(2*c1);
else
CdSv=0;
V(i3)=0;
end
h(i3)=ym-Bf*V(i3);
%Enregistrement RPI
tps(k)=t;
Vi1(k)=V(i1);
Vi2(k)=V(i2);
Vi3(k)=V(i3);
ye(k)=h(i1);
D
Pi3(k)=(h(i3)-z)*gamma;
V1(k)=V(1);
end
figure
plot(tps,V1,'b',tps,Vi1,'g')
xlabel('Temps t en (s)')
ylabel('Vitesse en (m/s)')
grid on
figure
plot(tps,ye,'b')
legend('ye')
xlabel('Temps t en (s)')
ylabel('ye en (m)')
grid on
