TP N° 1 : Ecoulement de puissance
But de TP :
Calcules d’un programme MATLAB et traçage les courbes de écoulement de
puissance.
Part 1:
Théorique :
L’exemple ci-dessous concerne l’écoulement de puissance entre deux sources
de tension.
Dans ce premier cas nous allons fixer l’angle δ de la source 1 à 5°. (δ=5 ° fix ¿
a) Calcule la puissance complexe pour chaque source :
R=1 Ω , X=7 Ω donc : Z=R+jX = 1+7j
V1= 120. [cos(5π/180)+[Link](5π/180)]= 119.54 + 10.46 j
V2= 100. [cos(0)+[Link](0)]= 100 + 0 j = 100
(V1-V2)= Z × I12 donc : I12 = (V1-V2)/ Z
(V2-V1)= Z × I21 donc : I21 = (V2-V1)/ Z = - (V1-V2)/ Z = - I12
I12 = (V1-V2)/ Z = [(119.54 + 10.46 j)-(100+0j)] / (1+7j) = 1.86 - 2.53j
I21 = - I12 = -(1.86 - 2.53j) = -1.86 + 2.53j
La puissance apparente d’un la source 1 :
* *
S1 = V1 × I12 = (119.54 + 10.46j) × (1.86-2.53j) = (119.54 + 10.46j) × (1.86+2.53j)
S1 = 195.33 + 321.48 j
La puissance active et réactive d’un la source 1 :
S1 = P1 + Q1 . j donc : P1 est la partie real et Q1 est la partie imaginer
(P1 : La puissance active d’un la source 1)
(Q1 : La puissance réactive d’un la source 1)
Donc :
P1 = 195.33 w
Q1 = 321.48 var
S1 = √ ¿ ¿
S1 = 376.169 VA
----------------------------------------------------------------------------------------------------------
La puissance apparente d’un la source 2 :
* *
S2 = V2 × I21 = (100 + 0j) × (-1.86 + 2.53j) = (100 + 0j) × (-1.86 - 2.53j)
S2 = - 186- 253j
La puissance active et réactive d’un la source 2 :
S2 = P2 + Q2 . j donc : P2 est la partie real et Q2 est la partie imaginer
(P2 : La puissance active d’un la source 2)
(Q2 : La puissance réactive d’un la source 2)
Donc :
P2 = 186 w
Q2 = 253 var
S2 = √ ¿ ¿
S2 = 314 VA
b) Calcule les pertes :
Pa= R × (|I 12|)2 = R × (|I 21|)2 = 1 × (√ ¿ ¿)2 =¿
Pa = 9,8605 w (perte active)
Qr= X × (|I 12|)2 = X × (|I 21|)2 = 7 × (√ ¿ ¿)2 =7 . ¿]
Qr = 69,0235 var (perte réactive)
Calcule dans le programme MATLAB:
Le programme script :
clc;
clear all ;
%les donnée
R=1; X=7; Z=R+j*X;
V1=120*(cos(5*pi/180)+j*sin(5*pi/180));
V2=100+j*0;
I12=(V1-V2)/Z;
I21=-I12;
%la puissance
S1=V1*conj(I12)
P1=real(S1)
Q1=imag(S1)
S2=V2*conj(I21)
P2=real(S2)
Q2=imag(S2)
%les perts
Pa=R*abs(I12)^2
Qr=X*abs(I12)^2
Exécution
S1 = 195.33 + 321.48i
P1 = 195.33
Q1 = 321.48
S2 = -185.51 - 252.69i
P2 = -185.51
Q2 = -252.69
Pa = 9.8265
Qr = 68.786
Conclusion
Rapprochement entre le calcul théorique et l'arithmétique appliquée.
Part 2:
Dans cette deuxième partie on fixer l’angle de phase de la source 2 à 0° et on
fait varie l’angle de phase de la source 1 de ± 30° par un pas de 5 ° les
amplitudes des tensions des deux sources sont maintenus constants.
δ =[-30° : 5° : 30°] varie
R=1 Ω,
X=7 Ω
La même calcul dans la partie 1 la seule différence on varie l’angle de phase de
la source 1 de ± 30° par un pas de 5 °.
On calcul (la puissance active, réactive et apparente et les pertes de la ligne) et
traçages les courbes (P1,P2 et PL ; Q1,Q2 et QL en fonction de l’angle δ ) sur le
programme MATLAB.
Calcule et traçage dans le programme MATLAB:
Le programme script :
clc
clear
clear all
R=input('Resistance de la ligne=');
X=input('Riactance de la ligne=');
E1=input('Amplitude de la tenstion de phase2=');
E2=input('Amplitude de la tenstion de phase1=');
a2=input('Angle de phase de la source2=');
a1=(-[Link])';
ar1=a1*(pi/180);
Z=R+j*X;
V1=E1.*(cos(ar1)+j*sin(ar1));
V2=E2.*(cos(a2*(pi/180))+j*sin(a2*(pi/180)));
I12=(V1-V2)/Z; I21=-I12;
S1=V1.*conj(I12);
P1=real(S1); Q1=imag(S1);
S2=V2.*conj(I21);
P2=real(S2); Q2=imag(S2);
Pl=R.*(abs(I12)).^2;
Ql=X.*(abs(I12)).^2;
%tracage puissances actives
Resultats01=[a1,P1,P2,Pl];
disp(' a1 P1 P2 Pl');
disp(' ');
disp(Resultats01)
plot(a1,P1,'r')
grid; hold;
plot(a1,P2,'b')
plot(a1,Pl,'g')
Xlabel('angle de phase en degree')
Ylabel('puissance en watt')
gtext('P1')
gtext('P2')
gtext('Pl')
%tracage puissances réactives
%Resultats02=[a1,Q1,Q2,Ql];
%disp(' a1 Q1 Q2 Ql');
%disp(' ');
%disp(Resultats02)
%plot(a1,Q1,'r')
%grid; hold;
%plot(a1,Q2,'b')
%plot(a1,Ql,'g')
%Xlabel('angle de phase en degree')
%Ylabel('puissance en var')
%gtext('Q1')
%gtext('Q2')
%gtext('Ql')
Exécution
a1 P1 P2 Pl
-30.00000 -759.84610 832.15390 72.30781
-25.00000 -639.51255 692.48481 52.97226
-20.00000 -512.12007 549.06761 36.94754
-15.00000 -378.63819 402.99380 24.35560
-10.00000 -240.08280 255.37508 15.29228
-5.00000 -97.50838 107.33492 9.82654
0.00000 48.00000 -40.00000 8.00000
5.00000 195.33492 -185.50838 9.82654 résultat dans la partie (1) δ=5 °
10.00000 343.37508 -328.08280 15.29228
15.00000 490.99380 -466.63819 24.35560
20.00000 637.06761 -600.12007 36.94754
25.00000 780.48481 -727.51255 52.97226
30.0000 920.1539 -847.8461 72.3078
a1 Q1 Q2 Ql
-30.0000 681.0773 -174.9227 506.1546
-25.0000 594.8313 -224.0255 370.8058
-20.0000 519.4012 -260.7684 258.6328
-15.00000 455.36118 -284.87196 170.48922
-10.00000 403.19854 -296.15259 107.04595
-5.00000 363.31029 -294.52447 68.78581
0.00000 336.00000 -280.00000 56.00000
5.00000 321.47553 -252.68971 68.78581 résultat dans la partie (1) δ=5 °
10.00000 319.84741 -212.80146 107.04595
15.00000 331.12804 -160.63882 170.48922
20.00000 355.23156 -96.59877 258.63279
25.00000 391.97453 -21.16870 370.80584
30.00000 441.07732 65.07732 506.15464
P1
P2
PL
Q1
Q2
QL
Interprétation des résultats
A noter que par suite de la première correspondance avec le résultat
de la seconde partie dans le cas où ¿ 5 ° , Cela signifie que les résultats
exacts on variation de l’angle δ , et puis par le traçage.
