PARTIDA CON AUTOTRANSFORMADOR CONECTADO A BARRA INFINITA

...

Los siguientes datos se obtuvieron mediante simulaciones:

Magnitud

Valor

[pic 24]

0.3474[pu]

[pic 25]

5.8[pu]

[pic 26]

2.17 [pu]

Tiempo de Reaceleración

1.85[s]

ANEXOS

I. Partida reducida con barra infinita

%% Partida con Tap

clear all

close all

clc

%Datos de Placa

P_rtd = 55e3; V_rtd = 400; N_rtd = 982 ; N_sinc = 1000;

T_rtd = P_rtd/(N_rtd*pi/30); I_rtd = 99; FP_rtd = 0.86;

Jm = 1.5;

%Inercia de la carga e Inercia Equivalente

JL = 14*Jm; J = Jm + JL;

%Características de Carga.

wm_pu = [0:0.05:0.9 0.92 0.95 1];

Te_pu = [2.60 2.50 2.45 2.40 2.35 2.28 2.2 2.17 2.13 2.13 2.13 2.1 2.13 2.13 2.13 2.17 2.25 2.35 2.38 2.4 2.25 0];

Ia_pu = [6.2 6.1 6.1 6.1 6.02 6.00 5.98 5.8 5.7 5.5 5.3 5.15 5.00 4.85 4.55 4.3 4.00 3.5 2.8 2.35 1.78 0.6];

PF_pu = [0.31 0.31 0.31 0.31 0.33 0.33 0.34 0.375 0.38 0.39 0.41 0.43 0.45 0.47 0.51 0.54 0.6 0.66 0.8 0.86 0.87 0.5];

% Interpolación por medio de Spline

Nptos = 1000+1;

N_spl = linspace(0,1.00,Nptos);

wm = N_spl;

Te = spline(wm_pu, Te_pu, N_spl);

Ia = spline(wm_pu, Ia_pu, N_spl);

PF = spline(wm_pu, PF_pu, N_spl);

%Torque de Carga

TL = 0.9*N_spl/(N_rtd/N_sinc);

%Variables en función del tiempo

for k=1:Nptos

wm_t(k) = wm(k);

Valim_t(k) = 1;

if 0.8^2*Ia(k)>1.5

Vmot_t(k) = 0.8*Valim_t(k);

Te_t(k) = 0.8^2*Te(k);

Ia_t(k) = 0.8^2*Ia(k);

Tacel_t(k) = 0.8^2*Te(k)-TL(k);

else

Vmot_t(k) = Valim_t(k);

Te_t(k) = Te(k);

Ia_t(k) = Ia(k);

Tacel_t(k) = Te(k)-TL(k);

end

Pmot_t(k) = sqrt(3)*Vmot_t(k).*Ia_t(k).*PF(k);

Palim_t(k) = sqrt(3)*Valim_t(k).*Ia(k).*PF(k);

Ialim_t(k) = Ia(k);

PF_t(k) = PF(k);

TL_t(k) = TL(k);

if wm(k) >= 0.98

t_lenght = k;

break

end

end

%Tiempo en que se evaluará la operación

t(1) = 0;

for k=2:t_lenght

delta_wm = wm(k) - wm(k-1);

delta_t(k)=J*delta_wm/Tacel_t(k)*(N_sinc*(pi/30)/T_rtd);

t(k) = t(k-1)+delta_t(k);

if t(k)>4.298 & t(k)

wcambioetapa = wm(k);

end

end

%Gráficos

f1=figure;clf;

set(f1, 'name', 'Motor & Load Curves[pu]')

subplot(3,1,1), pl = plot(wm, Te, '-b', wm_pu, Te_pu, 'ob'...

, wm, TL,'r', 'linewidth',2);

title('Motor & load curves [pu]')

ylabel('Torque[pu]')

legend('Torque Eléctrico','Puntos T_e [pu]', 'Torque de carga')

grid on

axis([wm(1) wm(Nptos) -.5 4.0])

subplot(3,1,2), plot(wm, Ia, '-b', wm_pu, Ia_pu, 'ob', 'linewidth',2);

ylabel('I_a[pu]')

grid on

axis([wm(1) wm(Nptos) 0.0 9])

subplot(3,1,3), plot(wm, PF, '-b', wm_pu, PF_pu, 'ob', 'linewidth',2);

ylabel('PF[pu]')

axis([wm(1) wm(Nptos) 0.0 1])

set(gca, 'fontsize',9)

xlabel('Speed [pu]')

grid on

f2=figure;

set(f2, 'name', 'Variables as Function of time')

subplot(3,1,1), plot(t, wm_t*100,