Modelluntersuchungen betreffend die Stabilität - Institut für ...
Modelluntersuchungen betreffend die Stabilität - Institut für ...
Modelluntersuchungen betreffend die Stabilität - Institut für ...
Erfolgreiche ePaper selbst erstellen
Machen Sie aus Ihren PDF Publikationen ein blätterbares Flipbook mit unserer einzigartigen Google optimierten e-Paper Software.
Anhang<br />
% Berechnung der Leitungswinkel in Bogenmaß<br />
psi11m=angle(z11m);<br />
psi12m=angle(z12m);<br />
psi21m=angle(z21m);<br />
psi22m=angle(z22m);<br />
% Berechnung der Komplementärwinkel in Bogenmaß<br />
a11m=(pi/2) - psi11m;<br />
a12m=(pi/2) - psi12m;<br />
a21m=(pi/2) - psi21m;<br />
a22m=(pi/2) - psi22m;<br />
% Betragsberechnung der Kuppel- und Eigenimpedanzen<br />
z11mabs=abs(z11m)<br />
z22mabs=abs(z22m)<br />
z12mabs=abs(z12m)<br />
% Array-Initialisierung von Epeff,P1tges,P2tges,Q1mges,Q2mges,thetagesg,u1eff<br />
Epeff_m=[];<br />
P1tges_m=[];<br />
P2tges_m=[];<br />
Q1mges_m=[];<br />
Q2mges_m=[];<br />
thetagesg_m=[];<br />
ukleff_m=[];<br />
% Beginn der ersten While Schleife zur Berechnung der Leistungen bei konstanter Generatorklemmenspannung Ukl<br />
%Initialisierung des Polradwinkels und der Polradspannung<br />
thetalg=0;<br />
Epeff=0;<br />
%Solange <strong>die</strong> Polradspannung kleiner 288 V ist wird in der Schleife<br />
%in Abhängigkeit des Polradwinkels <strong>die</strong> Polradspannung bei konstanter Generatorklemmenspannung berechnet.<br />
while Epeff < 288<br />
thetalg = thetalg +0.1;<br />
thetalr=thetalg*pi/180;<br />
%Solange <strong>die</strong> maximale Polradspannung nicht erreicht ist ,soll <strong>die</strong> Klemmenspannung des Generators auf konstante 127 V<br />
gehalten werden.<br />
ukleff=127;<br />
%Die Netzspannung ist betragsmäßig und winkelmäßig Bezugsgröße (127 Veff,Winkel=0 Grad)<br />
u2eff=127;<br />
%Generatorklemmenspannung komplex (Index k ... komplex)<br />
uklk=(ukleff*cos(thetalr))+(i*ukleff*sin(thetalr));<br />
%Netzspannung komplex<br />
u2k=u2eff;<br />
% Berechnung des in <strong>die</strong> Leitung fließenden (komplexen) Stromes I1<br />
I1k=y11m*uklk+y12m*u2k;<br />
%Komplexer Spannungsabfall an der Synchronreaktanz<br />
duk=I1k*zd;<br />
%Rückgerechnete komplexe Polradspannung<br />
Epk=duk+uklk;<br />
%Effektivwert der Polradspannung<br />
Epeff=abs(Epk);<br />
%Berechnung des gesamten Polradwinkels (zwischen Ep und U2=Unetz)<br />
thetagesr=atan2(imag(Epk),real(Epk));<br />
thetagesg=(thetagesr*180)/pi;<br />
% Berechnung der ges. Wirkleistungen P1 und P2 <strong>für</strong> <strong>die</strong> gesamten Übetragungsstrecke (a12,a11,a22 und z11abs usw.)<br />
% Berechnung der ges. Blindleistungen Q1 und Q2 nur <strong>für</strong> <strong>die</strong> Leitung und Transformator (a12m,a11m,a22m und<br />
z11mabs usw.)<br />
%Solange Polradspannung nicht ganz ausgeregelt ist (kleiner 288 V) , ist u1eff auf konstante 127 V<br />
ukleff=127;<br />
%Berechnung der Faktoren A1t, A2t und Bt zur Ermittlung der übertragbaren Wirkleistungen<br />
A1t=(Epeff*Epeff*sin(a11)) / z11abs;<br />
A2t=(u2eff*u2eff*sin(a22)) / z22abs;<br />
Bt=(Epeff*u2eff) / z12abs;<br />
%Berechnung der Faktoren A1m, A2m und Bm zur Ermittlung der übertragbaren Blindleistung (nur Leitung und<br />
Transformator!)<br />
A1m=(ukleff*ukleff*cos(a11m)) / z11mabs;<br />
A2m=(u2eff*u2eff*cos(a22m)) / z22mabs;<br />
Bm=(ukleff*u2eff) / z12mabs;<br />
%Berechnung der übetragbaren Wirk und Blindleistung<br />
P1tges=3*(A1t + (Bt*sin(thetagesr-a12)));<br />
P2tges=3*(-A2t +(Bt*sin(thetagesr+a12)));<br />
Q1mges=3*(A1m - (Bm*cos(thetalr-a12m)));<br />
Diplomarbeit<br />
93