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

Anhang
8 Anhang
8.1 Matlab – M.Files
8.1.1 Vergleich der Maximalwerte der übertragbaren Leistung
% Vergleich der Maximalwerte der Übertragbaren Leistung
hold on;
%Array von EUOR, EUMR und EUMRCE initialisieren
%Array der Längenwerte initialisieren
L=[100 200 300 400 500 600];
EUOR=[310.44 260.28 226.71 203.23 186.38 174.20];
EUMR=[711.13 585.10 463.57 359.86 296.51 265.78];
EUMRCE=[877.41 796.37 743.79 711.53 667.11 608.16];
%Plotbefehl
plot(L,EUOR,'-',L,EUMR,'--',L,EUMRCE,':')
grid on
legend('Energieübertragung ohne Spannungsregler','Energieübertragung mit Spannungsregler','Energieübertragung mit
Spannungsregler und TCSC ')
xlabel('Leitungslänge [km]')
ylabel('max. übertragbare Leistung [Watt]')
title('Maximal übertragbare Wirkleistung in Abh. der Leitungslänge')
8.1.2 Kuppel- und Eigenimpedanz als Funktion der Leitungslänge
%Programm zur Berechung der Kuppel- und Eigenimpedanzen als Fkt. der Leistungslänge
format long
z12abs_m=[];
z11abs_m=[];
z22abs_m=[];
X_m=[];
for X= 0 : 10 : 2000
% Berechnung von Omega
w=2*pi*50;
% Widerstandsbelag der Freileitung
Rstrich=0.03;
% 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=1/(z3p);
% Umrechnung Admittanz-Impedanz
z1p=1/(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);
Diplomarbeit
87