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

Anhang
% Induktivitätsbelag der Freileitung
XLstrich=0.3;
% Berechnung der Betriebsreaktanz XL und der Induktivität L
XL=(XLstrich*X);
L=(XL/w);
% Halber wert der Betriebskapazität, laut Ifea Modell
cbh=L/(2*274*274)
% Berechnung des Betriebswiderstandes RL
RL=Rstrich*X;
% Berechnung der Suszeptanzen der Freileitung, wobei der Index p für Pi-Ersatzschaltung steht
y1p=w*cbh*i;
y2p=y1p
% Berechnung der Betriebsimpedanz der Freileitung bzw. der Betriebsadmittanz
z3p=RL+XL*i
y3p=inv(z3p)
% Umrechnung Admittanz-Impedanz
z1p=inv(y1p)
z2p=z1p
% Impedanzkorrektur
z3pk=((sinh(sqrt(2*z3p/z2p)))/(sqrt(2*z3p/z2p)))*z3p
z2pk=((sqrt(z3p/(2*z2p)))/(tanh(sqrt(z3p/(2*z2p)))))*z2p
z1pk=z2pk
% Umrechnung Impedanz-Admittanz
y1pk=inv(z1pk)
y2pk=inv(z2pk)
y3pk=inv(z3pk)
% Vereinfachung-Faktorisierung
k1=y1pk*y2pk
k2=y2pk*y3pk
k3=y1pk*y3pk
k=k1+k2+k3
% Berechnung der Ersatzimpedanzen der T Ersatzschaltung
z1t=y1pk/k
z2t=y2pk/k
z3t=y3pk/k
% Parameter der Trafo und Synchronreaktanz
zd=118*i;
zt=5.74*i;
% Zusammenfassung der Impedanzen
z1g=z1t+zd+zt
z2g=z2t
z3g=z3t
% Vereinfachung-Faktorisierung
k4=z1g*z2g
k5=z2g*z3g
k6=z1g*z3g
k7=k4+k5+k6
% Berechnung der Ersatzadmittanzen der Pi Ersatzschaltung
y1pg=z1g/k7
y2pg=z2g/k7
y3pg=z3g/k7
% Anschreiben in Matrixform
y11=y3pg+y2pg
y12=-(y3pg)
y21=-(y3pg)
y22=y1pg+y3pg
Y=[y11 y12;y21 y22]
% Umrechnung auf Kuppel und Eigenimpedanzen der Pi Ersatzschaltung
z11=inv(y11)
z12=-(inv(y12))
z21=-(inv(y21))
z22=inv(y22)
% Anschreiben in Matrixform
z=[z11 z12;z21 z22]
%Initialisierung der Spannungen
u1=127
u2=127
% Berechnung des Leitungswinkels in Bogenmaß
psi11=angle(z11)
Diplomarbeit
89