RTH SUR UNE ADDICTION GRAVITAIRE ALIMENTEE PAR DEUX RESERVOIRS
RANDRIAMIRADO Misaela
INTRODUCTION
En généralement, 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
sns…cette étude visé pour observer et apprécier la réaction de l’ensemble de l’installation, déterminer
les zones à risque et calculer les sollicitation maximales pour dimensionner correctement, prévoir et
dimensionner des dispositifs de protection. Mais les but de notre étude à base sur la détermination de
valeur maximale de la pression pour dimensionner correctement de la conduite c’est-à-dire pour qu’elle
résiste à la contrainte maximale qui se produit lors du Régime Transitoire Hydraulique.
1
1. Description du système physique
1.1 Schéma du système physique
Figure 1:Installation d’embranchement à double réservoirs
1.2 Explication du principe de fonctionnement
Le système représente sur la figure au-dessus est une installation d’embranchement à double réservoir.
La conduite forcée à caractéristique unique (D, A, f constant tout le long de la conduite) de longueur L
alimentée par un réservoir à niveau constant sous une hauteur HD. La conduite soumis une vanne
partiellement ouverte à l’extrémité aval.
1.3 Données
Voici les données pour faire la simulation de système.
Tableau 1: Données de l'installation
Dimension Valeur Unité
D 0,15 m
f 0,02 -
Hd 30 m
CdSv0 0,012 m2
tmax 30 s
A 1200 m/s
g 9,81 N/kg
tf 3 s
Dx 50 m
2. Mise en équation
2
Equation de mouvement :
V V h fV V
+V +g + =0
t x x 2D (1)
Equation de continuité :
h h A2 h
+V + =0 (2)
t x g x
2.1 Conditions initiales
Pour trouver Vi et hi
On applique le théorème de Bernoulli entre (A)surface libre dans le réservoir et (J) dans le jet libre
issu de vanne.
PA V 2A P V2
+ zA + − h = J + zC + J (3)
2g 2g
2 gH D
D’après la démonstration de l’équation de Bernoulli on a trouvé : V ( 0, x ) = V0 =
S2 fL
+
( CD SV )0 D
2
fx V 20
et h ( 0, x ) = H D − . , alors
D 2g
A t =0 pour i=1 à 18
Vi = V0 (4)
hi = H D − f
( i − 1) x . V 20 (5)
D 2g
2.2 Condition limite
Pour entrer on a une méthode des caractéristiques :
A l’entrée du conduite 1 :
Pour i=1, h1 = H D ;
h1 − y N
V1 = avec yN = h2− − BV2− + RV2− V2− et Q1 = V1.S (6)
B
Sections internes :
Pour i=2 à 5
3
yM = hi−−1 + BVi −−1 + RVi −−1 Vi −−1 et (7)
yN = hi−+1 − BVi +−1 + RVi +−1 Vi −−1 (8)
yM + y N y − yN
Avec hi = et Vi = M (9)
2 2.B
A la sortie du conduite 1 avec l’embranchement :
Pour i=6 ; i=10 ; et i=11
yM 1 = h5− + BV5− − RV5− V5− (10)
yM 2 = h9− + BV9− − RV9− V9− (11)
yN 3 = h12− + BV12− − RV12− V12−
La hauteur piézométrique commune au niveau de l’embranchement
Alors : hi1 = hi 2 et hi1 = hi 3 donc hi1 + hi 2 + hi 3 = 3hi1 et on a trouvé
yM 1 + yM 2 + y N 3
h11 = h10 = h6 = (12)
3
yM 1 − h6
V6 = (13)
B
A l'entrée du conduite 2 :
Pour i=7 h7 = H D ;
h7 − y N
V7 = avec (14)
B
yN = h8− − BV8− + RV8− V8− et Q7 = V7 .S (15)
Sections internes :
Pour i=8 à 9
yM = hi−−1 + BVi −−1 − RVi −−1 Vi −−1 et (16)
yN = hi−+1 − BVi +−1 + RVi +−1 Vi −−1 (17)
yM + y N y − yN
Avec hi = et Vi = M (18)
2 2.B
A la sortie du conduite 2 :
yM 2 = h9− + BV9− − RV9− V9− (19)
yM − h10
h10 = h6 et V10 = (20)
B
A l'entrée du conduite 3 :
4
yM 3 = h12− − BV12− + RV12− V12− (21)
h11 − y N
h11 = h6 et V11 = (22)
B
Sections internes :
Pour i=12 à 17
yM = hi−−1 + BVi −−1 − RVi −−1 Vi −−1 et (23)
yN = hi−+1 − BVi +−1 + RVi +−1 Vi −−1 (24)
yM + y N y − yN
Avec hi = et Vi = M (25)
2 2.B
A la sortie du conduite 3 :
yM = h17− + BV17− − RV17− V17− (26)
Fermeture de la vanne : si t < tf
1.5
t
CD SV = ( CD SV )0 . 1 − (27)
tf
2
S 1
C1 = . , C2 = B, C3 = − yM (28)
CD SV 2 g
= C2 2 − ( 4.C1C3 ) (29)
−C2
V18 = ; (30)
2.C1
Si non
CD SV = 0;V18 = 0 (31)
h18 = yM − ( BV
. 18 ) ; Q18 = V18 .S (32)
3. Résolution numérique des équations
5
3.1 Procédure de calcul
Données :
H D = 30; ( CD SV )0 = 0.012; tmax = 30; A = 1200; g = 9.81; tf = 3
x = 50; D = 0.15; f = 0.02
Constantes :
D2 L f x A x
S= ; x = ; R= ; B = ; t = ;
4 n 2 gD g A
L1 = 5.x; L2 = 3.x; L3 = 7.x
f .L1 f .L2 f .L3
K1 = 2
; K2 = 2
; K3 =
2.g.D.S 2.g.D.S 2.g.D.S 2
a = K1 − K 2 ; b = 2.K 2 ; c0 = − K 2 ; = b 2 − ( 4.a.c0 ) ;
−b S 1
C= ; K4 = . ;
( CD SV )0 2.g.S 2
2 −1
2.a
HD
Q30 = ; Q10 = C.Q30 ; Q20 = (1 − C ) .Q30 ;
K1.C + K 3 .K 4
2
Q10 Q Q
V10 = ;V20 = 20 ;V30 = 30 ;
S S S
Régime permanant initial : t = 0; k = 1;
i =1 à 6
Q10
V (i ) =
S
f ( i − 1) xV
. 210
h (i ) = H D .
2 gD
Q ( i ) = V ( i ) .S
i = 7 à 10
2
1
6
1
2
Q20
V (i ) =
S
f ( i − 7 ) x.V 2 20
h (i ) = H D .
2 gD
Q ( i ) = V ( i ) .S
i = 11 à 18
Q30
V (i ) =
S
f ( i − 11) xV
. 230
h (i ) = H D .
2 gD
Q ( i ) = V ( i ) .S
Régime variable :
t<tmax ; t= t+dt
k = 1+k
i = 1à18
h i− = h i
Vi − = Vi
3 4
7
5
A l'entrée du conduite 1 :
h1 = H D
y N = h2− − BV2− + RV2− V2−
h1 − YN
V1 =
B
Q1 = V1.S
i = 2à5
Sections internes 1 :
yM = h(−i −1) + BV(i−−1) − RV(i−−1) V(i−−1)
y N = h(−i +1) − RV( i−+1) + RV(i−+1) V(i−+1)
yM + y N y − yN
hi = ; Vi = M ;
2 2B
A la sortie du conduite 1 avec l’embranchement :
yM 1 = h5− + BV 5− − RV5− V5−
yM 2 = h9− + BV9− − RV9− V9−
y N 3 = h12− − BV12− + RV12− V12−
yM 1 + yM 2 + y N 3
h6 =
3
y −h
V6 = M 1 6
B
8
6
A l'entrée du conduite 2 :
h7 = H D
y N = h8− − BV8− + RV8− V8−
h7 − y N
V7 =
B
Q7 = V7 .S
Sections internes 2 :
i = 8à9
yM = h(−i −1) + BV(i−−1) − RV(i−−1) V(i−−1)
y N = h(−i +1) − RV( i−+1) + RV(i−+1) V(i−+1)
yM + y N y − yN
hi = ; Vi = M ;
2 2B
A la sortie du conduite 2 :
h10 = h6
yM = h9− + BV9− − RV9− V9−
yM − h6
V10 = ;
B
A l'entrée du conduite 3 :
h11 = h6
yM = h12− + BV12− − RV12− V12−
h11 − y N
V11 = ;
B
Sections internes 3 :
i = 12à17
8
7
9
7 8 4
yM = h(−i −1) + BV(i−−1) − RV(i−−1) V(i−−1)
y N = h(−i +1) − RV( i−+1) + RV(i−+1) V(i−+1)
yM + y N y − yN
hi = ; Vi = M ;
2 2B
A la sortie du conduite 3 :
yM = h17− + BV17− − RV17− V17−
Fermeture de la vanne :
t<tf
1.5
t
CD SV = ( CD SV )0 . 1 −
tf
2
S 1
C1 = . , C2 = B, C3 = − yM
CD SV 2 g
= C2 2 − ( 4.C1C3 )
−C 2
V18 =
2.C1
CD SV = 0;V18 = 0
h18 = yM − ( B.V18 )
Q18 = V18 .S
Enregistrement :
tps ( k ) = t ; h18 ( k ) = h18 ; h11 ( k ) = h11 ; h9 ( k ) = h9 ; h3 ( k ) = h3 ;
Q18 ( k ) = Q18 ; Q7 ( k ) = Q7 ; Q1 ( k ) = Q1 ;
10
4. Résultats et discussion
Voici la simulation a fait sur le MATLAB, on a trouvé les résultats suivants donc il varie en fonction de
temps alors que ce temps qu’on parle ici est temps d’ouverture et temps du fermeture de la vanne. Le
voilà les résultats des grandeurs hydraulique (hauteurs et débits) qu’on a trouvé sur le MATLAB.
Figure 2 : Variation des hauteurs
➢ Interprétation de résultat figure 2 :
D’après le résultat que nous avons vu sur la figure ci-dessus, les quatre hauteurs piézométries ne sont
pas des mêmes allures, en fonction de temps tf = 6s. La variation de hauteur à la sortie de h18 est le plus
instable par rapport à l’autre suivant Hd=30m jusqu’au temps maximale de tmax=30s. 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.
11
Figure 3: Variation des débits
➢ Interprétation de résultat figure 3 :
D’après le résultat que nous avons vu sur la figure ci-dessus, les trois débits ne sont pas des mêmes
allures, en fonction de temps tf = 6s aussi. La variation de débit à la sortie Q18 est constante suivant
Hd=30m jusqu’au temps maximale de tmax=30s. parce que à temps du zéro jusqu’à 6s le débit fait le
transitoire après qu’il fait le régime permanent final jusqu’à tmax=30s et puis l’autre débit après la chute
de transitoire qui fait le mouvement d’oscillation.
12
Figure 4 : Variation des hauteurs
➢ Interprétation de résultat figure 4 :
Sur le résultat que on a trouvé au-dessus, les trois hauteurs piézométries ne sont pas des mêmes allures
aussi, en fonction de temps tf = 15s. La variation de hauteur à la sortie de h18 est le plus instable suivant
Hd=30m jusqu’au temps maximale de tmax=30s. car le régime de transitoire de zéro à 15s donc la
hauteur est plus de deux fois de la hauteur maximale après il fait le mouvement d’oscillation entre le
temps de tf =15s jusqu’à 30s.
13
Figure 5 : Variation des débits
➢ Interprétation de résultat figure 5 :
D’après le résultat que nous avons vu sur la figure ci-dessus, les trois débits ne sont pas des mêmes
allures, en fonction de temps tf = 15s aussi. La variation de débit à la sortie Q18 est constante suivant
Hd=30m jusqu’au temps maximale de tmax=30s. parce que à temps du zéro jusqu’à 15s le débit fait le
transitoire après qu’il fait le grandeur temps jusqu’à tmax=30s et puis l’autre débit après la chute qui
fait le mouvement d’oscillation.
14
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êmes
hauteurs ou de débit. Car le temps de la fermeture de la vanne, il y a deux temps différents tf = 6s et tf
=15s et puis l’installation d’embranchement à double réservoir est la hauteur piézométrique commune
au niveau de l’embranchement, la continuité donc la somme des débits sont égale zéro et le même
diamètre alors que la somme des vitesses est égale zéro aussi.
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.
15
ANNEXE
%% Régime Transistoire Hydraulique à deux réservoirs%
% Données :
D=0.15; f=0.02; Hd=30; CdSvo=0.012; tmax=30; A=1200; g=9.81; tf=15; Dx=50;
% Constants:
L1=5*Dx; L2=3*Dx; L3=7*Dx; Dt=Dx/A; S=pi*D^2/4; B=A/g; R=f*Dx/(2*g*D);
K1=f*L1/(2*g*D*S^2); K2=f*L2/(2*g*D*S^2); K3=f*L3/(2*g*D*S^2);
a=K1-K2; b=2*K2; co=-K2; delta=b^2-(4*a*co);
c=(-b+sqrt(delta))/(2*a); K4=((S/CdSvo)^2-1)*(1/(2*g*S^2));
Q30=sqrt(Hd/((K1*c^2)+K3+K4))
Q10=c*Q30
Q20=(1-c)*Q30
V10=Q10/S;
V20=Q20/S;
V30=Q30/S;
% Régime permanant initiale
t=0;k=1;
for i=1:6
V(i)=V10;
h(i)=Hd-((f*(i-1)*Dx*V10^2)/(D*2*g));
Q(i)=V(i)*S;
end
for i=7:10
V(i)=V20;
h(i)=Hd-((f*(i-7)*Dx*V20^2)/(D*2*g));
Q(i)=V(i)*S;
end
for i=11:18
V(i)=V30;
h(i)=h(6)-((f*(i-11)*Dx*V30^2)/(D*2*g));
Q(i)=V(i)*S;
end
% Enregistrement du RPI:
Tps(k)=t;h18(k)=h(18);h11(k)=h(11);h9(k)=h(9);h3(k)=h(3);
Q18(k)=Q(18);Q1(k)=Q(1);Q7(k)=Q(7);
while t<tmax
t=t+Dt; k=1+k;
% Affectation:
for i=1:18
Vm(i)=V(i);
hm(i)=h(i);
end
% A l'entrée du conduite 1:
h(1)=Hd;
yN=hm(2)-(B*Vm(2))+(R*Vm(2)*abs(Vm(2)));
V(1)=(h(1)-yN)/B;
Q(1)=V(1)*S;
a
% Sections internes:
for i=2:5
yM=hm(i-1)+(B*Vm(i-1))-(R*Vm(i-1)*abs(Vm(i-1)));
YN=hm(i+1)-(B*Vm(i+1))+(R*Vm(i+1)*abs(Vm(i+1)));
h(i)=(yM+YN)/2;
V(i)=(yM-YN)/(2*B);
end
% A la sortie du conduite 1 avec l’embranchement :
yM1=hm(5)+(B*Vm(5))-(R*Vm(5)*abs(Vm(5)));
yM2=hm(9)+(B*Vm(9))-(R*Vm(9)*abs(Vm(9)));
yN3=hm(12)-(B*Vm(12))+(R*Vm(12)*abs(Vm(12)));
h(6)=(yM1+yM2+yN3)/3;
V(6)=(yM1-h(6))/B;
% A l'entrée du conduite 2:
h(7)=Hd;
yN=hm(8)-(B*Vm(8))+(R*Vm(8)*abs(Vm(8)));
V(7)=(h(7)-yN)/B;
Q(7)=V(7)*S;
% Sections internes:
for i=8:9
yM=hm(i-1)+(B*Vm(i-1))-(R*Vm(i-1)*abs(Vm(i-1)));
YN=hm(i+1)-(B*Vm(i+1))+(R*Vm(i+1)*abs(Vm(i+1)));
h(i)=(yM+YN)/2;
V(i)=(yM-YN)/(2*B);
end
% A la sortie du conduite 2 :
yM=hm(9)+(B*Vm(9))-(R*Vm(9)*abs(Vm(9)));
h(10)=h(6);
V(10)=(yM-h(10))/B;
% A l'entrée du conduite 3 :
h(11)=h(6);
yN=hm(12)-(B*Vm(12))+(R*Vm(12)*abs(Vm(12)));
V(11)=(h(11)-yN)/B;
% Sections internes:
for i=12:17
yM=hm(i-1)+(B*Vm(i-1))-(R*Vm(i-1)*abs(Vm(i-1)));
YN=hm(i+1)-(B*Vm(i+1))+(R*Vm(i+1)*abs(Vm(i+1)));
h(i)=(yM+YN)/2;
V(i)=(yM-YN)/(2*B);
end
% A la sortie du conduite 3 :
yM=hm(17)+(B*Vm(17))-(R*Vm(17)*abs(Vm(17)));
% Fermeture de la vanne :
if t<tf
CdSv=CdSvo*(1-(t/tf))^1.5;
C1=((S/CdSv)^2-1)*(1/(2*g)); C2=B; C3=-yM;
Delta=C2^2-(4*C1*C3);
V(18)=(-C2+sqrt(Delta))/(2*C1);
b
else
CdSv=0; V(18)=0;
end
h(18)=yM-(B*V(18));
Q(18)=V(18)*S;
% Enregistrement:
Tps(k)=t;h18(k)=h(18);h11(k)=h(11);h9(k)=h(9);h3(k)=h(3);
Q18(k)=Q(18);Q1(k)=Q(1);Q7(k)=Q(7);
end
figure
plot(Tps,h18,'r',Tps,h11,'b',Tps,h9,'y',Tps,h3,'g')
legend('h18','h11','h9','h3')
xlabel('t')
ylabel('h')
figure
plot(Tps,Q18,'b',Tps,Q1,'g',Tps,Q7,'r')
legend('Q18','Q1','Q7')
xlabel('t')
ylabel('Q')
