EJERCICIOS DE MAQUINAS ELECTRICAS CON MATLAB

1.      Mediante MATLAB, grafique la corriente de campo que se requiere para alcanzar un factor de potencia unitario de operación para el motor sincrónico de 500 kW, 4160 V, trifásico, tetra polar, 1800 rpm que se opera a velocidad proporcional:

 

 codigo fuente de matlab:

clc

clear

%Parámetros del Sistema

P=5000;

Vta = 4160/sqrt(3);

polos=4;

ns=1800;

Rs=0.011;

Xs=0.12;

If=1:120;

%Se resuelve

for n = 1:101

    delta(n) = (pi/2.)*(n-1)/100;

    Ia(n)= (Eaf *exp(j*delta(n)) -  Veq)/(j*(Xs + Xeq));

    Va(n)= abs(Veq + j*Xeq*Ia(n));

    degrees(n)= 180*delta(n)/pi;

end

Vterm = 1.0;

for n = 1:101

    P(n) = (n-1)/100;

    deltat(n) = asin(P(n)*Xeq/(Vterm*Veq));

    Ia(n) = (Vterm *exp(j*deltat(n)) -  Veq)/(j*Xeq);

    Eaf(n) = abs(Vterm + j*(Xs+Xeq)*Ia(n));

end

%Graficar resultados

plot(If,P)

xlabel('Potencia [kW]')

ylabel('Corriente de campo [A]')

title('Corriente de campo vs. Potencia')

RESULTADO:

 

 

a)

clc

clear

%Parámetros del Sistema

S=145000000;

Vta = 13800/sqrt(3);

polos=72;

f=60;

%Circuito abierto

If=[100 200 300 400 500 600 700 800 900];

V=[2.27 4.44    6.68    8.67    10.4    11.9    13.4    14.3    14.5];

plot(If,V)

xlabel('Corriente de campo [A]')

ylabel('Voltaje [V]')

title('Voltaje vs Corriente de campo')

 

 

%Cortocircuito

If2=[0 710];

In2=[0 6070];

plot(If2,In2)

xlabel('Corriente de campo [A]')

ylabel('Corriente de cortocircuito [A]')

title('Corriente de cortocircuito vs Corriente de campo')

polyfit(If,V)

polyfit(If2,In2)