058

58 Rozdział 5

tau =0.025; % - stała czasowa zanikania

%    nieokresowych prądów zwarcia

taubisd =0.028; % - stała czasowa podprzejscowa podłużna

tauprimd =0.85; % - stała czasowa przejściowa podłużna

xq=1.85; % - reaktancja sychroniczna poprzeczna

xd=2; % - reaktancja synchroniczna podłużna

xprimd=0.39; % - raktancja przejściowa podłużna

xbisq=0.32; % - reaktancja podprzejsciowa poprzeczna

xbisd=0.28; % - reaktancja podprzejsciowa podłużna

Ug=1.05; % - napięcie generatora przed zwarciem

psi=0; % - początkowy kat napięcia

Um=sqrt(2)*Ug; % amplituda napięcia

t0=0; tmax=0.5; dt=0.001;

t=t0;dt;tmax;

omega=2*pi*50;

%

ct=xbisd.*(l/xbisd-l/xprimd).*exp(-t./taubisd)...

+(l/xprimd-l/xd).*exp(-t/tauprimd)+l/xd;

%

%skladowe prądy prądu zwarcia ia %skladowa 50 Hz

i50Hz = Um/xbisd*{ ct.*cos(omega.*t + psi )    );

% składowa nieokresowa

inok = Um/xbisd*( -0.5*exp(-t./tau).* ...

(    (l/xbisd+l/xbisq )*cos(psi) )    );

% składowa 100 Hz

ilOOHz = Um/xbisd*( -0.5*exp(-t./tau).*...

( -(l/xbisd-l/xbisq).*cos(2*omega.*t+psi)    )    );

% prąd zawrciowy w fazie A iA = i50Hz + ilOOHz + inok;

^poszczególne wykresy subplot(2,2,1),plot(t,iA,'k—¹); title('prad zwarciowy¹),grid on; subplot(2,2,2),plot(t,i50Hz, 'k--'); title('składowa okresowa 50Hz'),grid on; subplot(2,2,3),plot(t,ilOOHz,'k-.'); title(¹ składowa okresowa 100Hz'),gnid on; subplot(2,2,4),plot(t,inok,¹k:'); title('składowa nieokresowa'),grid on; return

Wynik działania funkcji rysz3fgs() pokazano na rysunku 5.6.