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