152

152 Rozdział 12

Numeryczne rozwiązanie układu równań różniczkowych wykorzystano do symulacji rozruchu silnika. W tym celu opracowano 2 funkcje:

sasdx() - funkcja definiująca układ równań różniczkowych, sas() - funkcja sterująca obliczeniami.

Analizie poddano rozruch silnika asynchronicznego o następujących parametrach zastępczych w jednostkach względnych:

•    rezystancja stojana R, = 0.04656, rezystancja wirnika R, = 0.01628,

•    reaktancja stojana X_s = 0.08346, reaktancja wirnika X_r = 0.08346,

•    reaktancja magnesowania X,„ - 3.755,

•    stała czasowa mechaniczna T_m = 2 s,

•    moment obciążenia m_obc = 1,

•    wymuszone przez sieć napięcie stojana U_s = 1.

function dxdt-sasxdt(t,x)

% funkcja definuje równania różniczkowe % sianu nieustalonego silnika asynchronicznego global ws Ki R2 Ll L2 Lm Lip L2p Tm Us mobc filx=x(l); f i 1y=x(2) ; fi2x=x(3); fi2y=x(4); s=x(5);

tsek=t/ws;

% pochodne zmiennych stanu u1x- Us; uly-0;

Df i l.x=-R1/Ll p* ( f ilx-Lm/L2*f i2x) + fily + ulx;

Dfily=-Rl/Llp*(fily-Lm/L2*fi2y) - filx + uly;

Dfi2x--R2/L2p*(fi2x-Lm/Ll*filx) + S*fi2y;

Df i 2y----R2/L2p* ( f i2y-Lm/Ll* f ily) - S*fi2x;

% moment elektromagnetyczny

me = Lm*I.2/I,lp* (£ily*fi2x-filx*fi2y) ;

Ds=(- me + mobc)/(ws*Tm);

dxdt=[Dfilx; Dfily; Dfi2x; Dfi2y; Ds];

return

function sas

% funkcja steruje symulacja rozruchu silnika asynchronicznego global ws R1 R2 Ll L2 Lm Lip L2p Tm Us mobc j =sqrt(-1) ;

% m-plik wyniki analizy [fd,msg]= fopen('sasout.m' , 'wt' ) ;

% dane silnika asynchronicznego ws = 100 *pi ;

Rs=0.04656; Rr=0.01628;

Xs=0.08346; Xr=0.08346; Xm=3.755;

Tm=2 ; mobc=l;

Us = 1; % napięcie stojana sztywne spocz=l; % rozruch silnika % parametry zastępcze silnika Rl=Rs; Ll=Xs +Xm; R2 =Rr; L2=Xr+Xm; Lm=Xm;

Llp=L 1 -Lrrb2 /L2 ; L2p=L2-Lm"2/Ll ;

% warunki początkowe strumieni stojana i wirnika