Modelluntersuchungen betreffend die Stabilität - Institut für ...

Anhang
Q2mges=3*(-A2m +(Bm*cos(thetalr+a12m)));
ukleff_m=[ukleff_m ukleff];
Epeff_m=[Epeff_m Epeff];
P1tges_m=[P1tges_m P1tges];
P2tges_m=[P2tges_m P2tges];
Q1mges_m=[Q1mges_m Q1mges];
Q2mges_m=[Q2mges_m Q2mges];
thetagesg_m=[thetagesg_m thetagesg];
%Ende erste While Schleife
end
% Begin der zweiten while Schleife
% Ab EP = 288 V kann nicht mehr von konstanter Generatorklemmenspannung ausgegengen werden. Es wird nun, solange
%der Polradwinkel kleiner 180 Grad ist
% von Ep auf ik zurückgerechnet, um dann auf die tatsächliche Generatorklemmenspannung rückzurechnen.
% Danach folgt wieder die Berechnung von Wirk und Blindleistung
while thetagesg < 180
thetagesg=thetagesg + 0.1;
thetagesr=thetagesg*pi/180;
Epeff=288;
u2eff=127;
%Netzspannung komplex
u2k=u2eff;
Epk=(Epeff*cos(thetagesr))+(i*Epeff*sin(thetagesr));
I1k=y11*Epk+y12*u2k;
duk=I1k*zd;
uklk=Epk-duk;
ukleff=abs(uklk);
%Berechnung der Faktoren A1t, A2t und Bt zur Ermittlung der übertragbaren Wirkleistungen
A1t=(Epeff*Epeff*sin(a11)) / z11abs;
A2t=(u2eff*u2eff*sin(a22)) / z22abs;
Bt=(Epeff*u2eff) / z12abs;
%Berechnung der Faktoren A1m, A2m und Bm zur Ermittlung der übertragbaren Blindleistung (nur Leitung und
Transformator!)
A1m=(ukleff*ukleff*cos(a11m)) / z11mabs;
A2m=(u2eff*u2eff*cos(a22m)) / z22mabs;
Bm=(ukleff*u2eff) / z12mabs;
%Berechnung der übetragbaren Wirk- und Blindleistung
P1tges=3*(A1t + (Bt*sin(thetagesr-a12)));
P2tges=3*(-A2t +(Bt*sin(thetagesr+a12)));
Q1mges=3*(A1m - (Bm*cos(thetalr-a12m)));
Q2mges=3*(-A2m +(Bm*cos(thetalr+a12m)));
ukleff_m=[ukleff_m ukleff];
Epeff_m=[Epeff_m Epeff];
P1tges_m=[P1tges_m P1tges];
P2tges_m=[P2tges_m P2tges];
Q1mges_m=[Q1mges_m Q1mges];
Q2mges_m=[Q2mges_m Q2mges];
thetagesg_m=[thetagesg_m thetagesg];
%Ende zweite While Schleife
end
% Maximalwert der übertagbaren Leistung,der Polradspannung und der Generatorklemmenspannung
p2maxt=max(P2tges_m);
epmax=max(Epeff_m);
uklmin=min(ukleff_m);
% Achsenbeschriftung und Plotkonfiguration
hold on;
plot(thetagesg_m,P1tges_m,'-',thetagesg_m,P2tges_m,'--')
plot(thetagesg_m,Q1mges_m,'-.m',thetagesg_m,Q2mges_m,':k')
plot(thetagesg_m,Epeff_m,'-c')
plot(thetagesg_m,ukleff_m,':b')
grid on
legend('In die Leitung hineingeschickte Wirkleistung P1','Aus der Leitung herauskommende Wirkleistung P2','In die Leitung
hineingeschickte Blindleistung Q1','Aus der Leitung herauskommende Blindleistung Q2','Polradspannung
Ep','Generatorklemmenspannung Ukl')
xlabel('Theta [°] (Winkel zwischen Polradspannung Ep und Netzspannung U2)')
ylabel('Wirkleistung [W],Blindleistung [Var],Polradspannung [V],Generatorklemmenspannung [V]')
title('Übertragene Leistung als Funktion des Polradwinkels (mit kompensierter Längsreaktanz der Leitung)')
% Anzeigen der berechneten Werte am Bildschirm (Epmax)
Teilung=p2maxt/10;
Diplomarbeit
94