Ein Beitrag zur rechnerischen Bestimmung von ...

Diplomarbeit Simulationsmodell
end;
Ruwrechts = input( ' Betrieberdungswiderstand der Trafostation 2 [ in Ohm ]: ' );
while isempty(Ruwrechts) == 1, %Abfrage ob Eingabe getätigt wurde,
Ruwrechts = input (' ' ); %wenn nein..dann erneute Abfrage !!
end;
end; % if eingabe = 0
ausgabe = 0;
%**********************Berechnung***************************************
%ro = durchschnittlicher spezifischer Bodenwiderstand in [Ohm - m ]
%dm = mittleren Abstandes Erdseil-Phasenseil in [m]
%d = Durchmesser des Erdseiles in [m]
%mue = Permeabilität des Seilmaterials
%R1,R2 = ohmsche Widerstandswerte für das Phasenseil bzw. Erdseil in [Ohm pro Spannfeld]
%L1,L2 = Längsinduktivität des Phasenseils bzw. Erdseils in [H pro Spannfeld]
%Lm = Gegeninduktivität Phasenseil/Erdseil in [H pro Spannfeld]
%Rm = mag. Widerstand der Gegeninduktivität in [Ohm pro Spannfeld]
mue0 = 1.26e-8; %Permeabilität von Luft in [H / cm]
mue1 = 1; %Permeabilität des Phasenseil-Leitermaterials in [H / cm]
mue2 = 1; %Permeabilität des Erdseil-Leitermaterials in [H / cm]
k = 1/(ro * 1e2); %k=spezifische Leitfähigkeit in [1e-6 Siemens / cm]
h = 0.0185 / sqrt( mue0 * k * 2 * pi * 50); %mittlere Eindringtiefe in [m] nach
POLLACZEK
L1 = 0.2 * ( log(2*h/d1) + 0.25 * mue1 )* 1e-3; %Berechnung der Längsinduktivität des
%Phasenseils in [H / km]
%Log ist in MATLAB der natürliche
%Logarithmus ln
L2 = 0.2 * ( log(2*h/d2) + 0.25 * mue2 )* 1e-3; %Berechnung der Längsinduktivität des
%Erdseils in [H / km]
Lm = 0.2 * log(h/dm) * 1e-3; %Berechnung der Gegeninduktivität
%Phasenseil/Erdseil in [H / km]
% Benötigte Werte für die Matlabsimulation
freq = 50; %Frequenz
Rm = 1e-5; %mag. Widerstand der Gegeninduktivität in [Ohm pro Spannfeld]
Lm = Lm*l; %Gegeninduktivität Phasenseil/Erdseil in [H pro Spannfeld]
L1 = L1*l; %Längsinduktivität des Phasenseils in [H pro Spannfeld]
L2 = L2*l; %Längsinduktivität des Erdseils in [H pro Spannfeld]
R1 = R1*l; %Widerstandswerte des Phasenseils in [Ohm pro Spannfeld]
Gabbauer Anton Seite 85