Maquinas Eléctricas-Chapman-5ta-edición

6.6 Variaciones en las características par-velocidad del motor de inducción 265
% Primero, inicializar los valores que requiere el programa.
r1 5 0.641;
% Resistencia del estator
x1 5 0.750;
% Reactancia del estator
r2 5 0.300;
% Resistencia del rotor con jaula sencilla
r2i 5 0.400;
% Resistencia de la jaula interna del rotor
% con motor de doble jaula
r2o 5 3.200;
% Resistencia de la jaula externa del rotor
% con motor de doble jaula
x2 5 0.500;
% Reactancia del rotor del motor de jaula sencilla
X2i 5 3.300;
% Reactancia de la jaula interna del rotor
% del motor de doble jaula
X2o 5 0.500;
% Reactancia de la jaula externa del rotor
% del motor de doble jaula
xm 5 26.3;
% Reactancia de la rama de magnetización
v_phase 5 460 / sqrt(3);
% Voltaje de fase
n_sync 5 1800;
% Velocidad síncrona (r/min)
w_sync 5 188.5;
% Velocidad síncrona (rad/s)
% Calcular el voltaje y la impedancia de Thevenin con las
% ecuaciones 6-41a y 6-43.
v_th 5 v_phase * ( xm / sqrt(r1^2 1 (x1 1 xm)^2) );
z_th 5 ((j*xm) * (r1 1 j*x1)) / (r1 1 j*(x1 1 xm));
r_th 5 real(z_th);
x_th 5 imag(z_th);
% Ahora, calcular la velocidad del motor para varios
% deslizamientos entre 0 y 1. Nótese que el primer valor
% de deslizamiento es 0.001 en lugar de 0 exactamente para
% evitar la división entre cero.
s 5 (0:1:50) / 50;
% Deslizamiento
s(1) 5 0.001;
% Evitar división entre cero
nm 5 (1 - s) * n_sync;
% Velocidad mecánica
% Calcular el par del rotor de jaula sencilla.
for ii 5 1:51
t_ind1(ii) 5 (3 * v_th^2 * r2 / s (ii)) / ...
(w_sync * ((r_th 1 r2/s(ii))^2 1 (x_th 1 x2)^2) );
end
% Calcular la resistencia y reactancia del rotor de doble
% jaula con este deslizamiento y luego utilizar estos valores
% para calcular el par inducido.
for ii 5 1:51
y_r 5 1/(r2i 1 j*s(ii)*x2i) 1 1/(r2o 1 j*s(ii)*x2o);
z_r 5 1/y_r;
% Impedancia efectiva del rotor
r2eff 5 real(z_r);
% Resistencia efectiva del rotor
x2eff 5 imag(z_r);
% Reactancia efectiva del rotor
% Calcular el par inducido del rotor de doble jaula.
t_ind2(ii) 5 (3 * v_th^2 * r2eff / s(ii)) / ...
(w_sync * ((r_th 1 r2eff/s(ii) )^2 1 (x_th 1 x2eff)^2) );
end
% Hacer la gráfica de las curvas par-velocidad
plot (nm,t_ind1,'b','LineWidth',2. 0);
hold on;
plot(nm,t_ind2,'k-','LineWidth',2.0);
xlabel('bf\itn_{m}'¡
ylabel('\bf\tau_{ind}';
title ('\bfCaracterísticas par-velocidad de un motor de inducción';
legend ('Diseño de jaula sencilla', 'Diseño de doble jaula');
grid on;
hold off;