164

164 Rozdział 13

Oznacza to. że znając wartości własne macierzy stanu, można wyznaczyć współczynniki tłumienia kołysań wirnika ze wzoru

(13.31)

__ Re(/ł) abs(/l)

W celu analizy stanu nieustalonego elektromechanicznego opracowano trzy funkcje:

•    stabdx - definiuje równania różniczkowe ruchu wirnika generatora,

•    stab - bada stabilność lokalną, a następnie po wciśnięciu ENTER symuluje zwarcie na generatorze i wyznacza przebieg czasowy kąta wirnika,

•    stabsym - tworzy okno GUI, w którym użytkownik może zmieniać współczynnik tłumienia.

funcŁion dxdt. = s tabdx (t, x) global Eprim Us X Pm Tm ws global tpzw tzw utratastab

global wm psi Dpusuw vslider vedit rys Dpu

dk=x(1); dw=x(2);

ddk=dw;

Pe = Eprim*Us/X*sin(dk) ; if t>=tpzw & t<=(tpzw+tzw)

Pe=0; % moc el. w czasie zwarcia

end

ddw=ws/Tm*(-Dpu*dw+Pm-Pe); dxdt=[ddk; ddw] ; i f abs(dw)>1000

fprintf('\nutrata Stabilności - przerwano obliczenia'); dxdt =[0;0] ; utratastab=l;

end

return

function stab(k,Eprim,Us,X,Pm,Tm,Dpu, tpzw, tzw) global Eprim Us X Pm Tm ws

global wm psi Dpusuw vslider vedit rys Dpu global tpzw tzw utratastab

%Opcjonalne wywołania z 7 parametrami, wg poniższych warunków; if nargin<9

tzw=0.2; % czas trwania zwarcia if r,argin<8

tpzw=l; % chwila wystąpienia zwarcia if nargin<7

Dpu=5; %wsp. tłumienia Dpu if nargin<6

Tm=12; %stala czasowa mechaniczna if nargin<5

Pm = 1; %moc mechaniczna turbiny if nargin<4

X = 0.5; % reaktancja układu GENERATOR - SYSTEM if nargin<3

Us=l; % napięcie SYSTEMU