Una aplicacin de Matlab al Motor de Induccin

Autor: M. Sc. Ing. Ignacio Romero Rueda

Profesor Asistente, Dpto. Ingeniera Elctrica, Instituto superior Minero Metalrgico de Moa.
Email: iromero@ismm.edu.cu
Resumen:
En el trabajo Una aplicacin de Matlab al Motor de Induccin se ofrece la teora general de un
motor de induccin de rotor en jaula en rgimen estable, las ecuaciones fundamentales de
momento, corriente y potencia a partir del circuito equivalente como va para el trazado de las
curvas caractersticas de corriente, mecnica y de funcionamiento, aspecto de gran inters en
el proceso de enseanza y prctica de la temtica de mquinas asincrnicas.
En la segunda parte se expone la realizacin a travs del Matlab de las caractersticas de
corriente, mecnica y de funcionamiento, con varias posibilidades a partir de un men de
seleccin y de interaccin con fines tcnicos y de aprendizaje con salida grfica y en forma de
datos numricos, que pueden ser obtenido en una base de dato creado por la aplicacin, dando
la posibilidad de realizar la corrida, con cambio de tensin y frecuencia cuantas veces se
precise.
La aplicacin permite tener en cuenta o no el efecto pelicular y asimetra de tensin. en
motores de diseo bsico en las caractersticas de funcionamiento y mecnica.
Palabras claves: motor de induccin, comportamiento asimtrico, efecto pelicular.
Generalidades del comportamiento energtico del Motor de Induccin.
Uno de los problemas principales en el funcionamiento del motor asincrnico es la
transformacin de la energa elctrica, consumida de la red, en energa mecnica. La
transformacin de energa de una forma a otra est inevitablemente enlazada con prdidas en
los diferentes rganos de la mquina. Por esta razn la energa, el flujo de potencia, las
prdidas y los momentos de rotacin representan un conjunto de problemas que tienen gran
importancia al estudiar las propiedades de trabajo del motor asincrnico.
Momentos de rotacin del Motor de Induccin.
El momento electromagntico M
em
desarrollado por el motor asincrnico en el rgimen
permanente de funcionamiento (n=const) vence el momento de carga en el eje, M
car
compuesto por el momento de marcha en vaco, M
0
y el momento de freno til M
2
.
;
2 0
M M M
em
+ = ;
2
0
n
p p p p
M
ad mec ad mec

+
=

+
= .
2
2 2
2
n
P P
M

= (1)
Donde P
2
es la potencia mecnica entregada por el motor en el eje y n es la velocidad de
rotacin del rotor del motor.
Momento electromagntico y corriente del Motor de Induccin.
La potencia mecnica total, correspondiente al momento electromagntico M
em
que se
desarrolla en el rotor del motor y se determina por la relacin:
, n 2 M M P
em em mec
= = de donde .
n 2
P P
M
mec mec
em

= = (2)
El momento M
em
surge como resultado de la interaccin del flujo magntico giratorio
m
y la
corriente I
2
en el rotor. Pero el flujo
m
gira en el espacio con la velocidad angular
1
=2n
1
donde n
1
=f
1
/p. La potencia desarrollada en este caso es la potencia electromagntica del
motor P
em.
Por lo tanto,
, M P
1 em em
= de donde .
n 2
P P
M
1
em
1
em
em

= = (3)
La expresin para la magnitud del momento electromagntico M
em
se puede obtener a travs
del circuito equivalente, en funcin de la potencia electromagntica transmitida del estator al
rotor P
em.
.
2
cos
2 2 2 2
p
f
M I E m P
em em

= = De aqu ( ) , cos 2 1
2 2 2 2 2
I k p m M
m dev em
=
si se tiene en cuenta que en el circuito en T equivalente esto corresponde a la potencia
liberada en la resistencia, s / r R
'
2 2
= y la corriente en el circuito en ngulo equivalente (L) es:
( )
,
X X
s
R
R
U
I
2
2 1
2
2
1
1 ''
2
+ + |
.
|

\
|
+
= (4)
La variacin de la corriente del rotor de acuerdo con la frmula (4) viene representada en la
fig.1,a).
Fig.1.Variacin de la corriente I

2
a) y mecnica b) en funcin del deslizamiento s
Caracterstica mecnica del motor de induccin
De acuerdo con el circuito equivalente en L la expresin para el momento electromagntico se
puede escribir en la forma:
( )
.
2
2
2 1
2
2
1 1
2 2
1 1
(
(

+ + |
.
|

\
|
+
=
X X
s
R
R f
s
R
U pm
M
em

(5)
As pues, para U
1
=const, la curva M
em
=f(s) tiene la forma representada en la fig.1,b). Cuando
los valores del deslizamiento son grandes, M
em
disminuye a pesar de que aumenta I
2
ya que
cos
2
en este caso disminuye rpidamente.
A partir de la ecuacin (5) se pueden deducir los momentos caractersticos del motor de
induccin:
Momento electromagntico mximo.
Es preciso determinar el deslizamiento para el cual el se obtiene el mismo, partiendo de los
datos del circuito en L :
.
X R
R
s
2
c . c
2
1
2
m
+
= (6)
En virtud de (5) y (6), para el momento electromagntico mximo hallamos:
| |
2
.
2
1 1 1
2
1 1
2 2
c c
mx
X R R f
U pm
M
+ +
=

(7).
La relacin
n mx m
M / M k = se llama factor de capacidad de sobrecarga del motor
asincrnico.
Momento de arranque del motor asincrnico.
A la par con el momento mximo, el momento de arranque del motor M
arr
es una de las
caractersticas de explotacin ms importantes del motor. El valor del momento M
arr
se obtiene
de la frmula general para el momento (5) siendo s=1:
( ) ( ) | |
.
X X R R f 2
R U pm
M
2
2 1
2
2 1 1
2
2
1 1
arr
+ + +
=

(8)
El momento M
arr
se expresa habitualmente por la relacin , M / M K
n arr arr
= la cual se llama
multiplicidad del momento de arranque.
Frmula de Kloos para el momento relativo.
Las relacin (5) determinan el valor del momento cualquiera que sea el deslizamiento y el
valor del momento mximo. En las mquinas asincrnicas, si no existe la resistencia adicional
en el circuito del rotor, habitualmente R
1
R
2
. Por eso, cuando la relacin R
1
/R
2
no se puede
establecer ms exactamente se puede aceptar que
( )
.
2
1 2
m
m
m
m
mx
s
s
s
s
s
s
M
M
+ +
+
= (9)
Conociendo por los datos de catlogos M
n
, k
m
y s
n
por la frmula (9) se puede determinar s
m
y
entonces por esta frmula se puede hallar la magnitud relativa del momento para cualquier
deslizamiento s.
Para pequeos valores del deslizamiento, en la zona de cargas normales, en la frmula (10) se
puede tener en cuenta en el denominador slo el trmino s
m
/s entonces
Si en la frmula (9) despreciamos s
m
en el numerador y el trmino 2s
m
en el denominador,
obtendremos la relacin aproximada, empleada a veces en los clculos:
.
s
s
s
s
2
M
M
m
m
mx
+
(10)
Caractersticas de trabajo del Motor de Induccin.
Por caracterstica de trabajo de un motor asincrnico se comprende:
( ) . const f const U siendo P f cos y , M , n
1 1 2
= = =
Adems, pertenecen tambin a los ndices de importancia el coeficiente de capacidad de
sobrecarga k
s
y, para los motores en cortocircuito, tambin el orden de la corriente de
arranque y del momento de arranque.
A. Velocidad de rotacin del motor ( ).
2
P f n = De la frmula para el deslizamiento se
desprende que ). 1 (
1
s n n = Por otro lado
m . e 2 cob
P / p s =
As pues, el deslizamiento del motor asincrnico numricamente es igual a la relacin de
las prdidas en el devanado del rotor p
cob2
a la potencia electromagntica P
em
; por eso aqu,
s0 y nn
1
. A medida que aumenta la carga crece la relacin (12), pero, teniendo en cuenta
el alto rendimiento del motor, est limitada entre lmites bastante estrechos.
Habitualmente, a carga nominal, %. 5 5 , 1 / /
2 2 . 2
= P p P p
cob m e cob
La cifra menor se refiere
a los motores de mayor potencia, del orden de 3 a 10 Kw. La dependencia
) (
2
P f n = representa una curva poco inclinada con respecto al eje de abscisas. El motor
posee una caracterstica de velocidad del tipo de la caracterstica del motor de corriente
continua en derivacin.
B. Dependencia ). (
2
P f M = en el rgimen permanente de funcionamiento
,
0 2
M M M + = donde M
2
es el momento de frenado del motor, y M
0
es el momento en
vaco. Puesto que al variar la carga en los lmites desde la marcha en vaco hasta la
nominal la velocidad de rotacin de los motores asincrnicos permanece casi constante, la
dependencia ) (
2
P f M = del motor asincrnico es casi rectilnea.
C. Prdidas y rendimiento del motor. En los motores asincrnicos tienen lugar las mismas
prdidas que en otras mquinas elctricas: prdidas mecnicas, prdidas en el acero,
prdidas en el cobre y prdidas adicionales. Las prdidas totales en el motor son
. p p p p p p
ad mec ac 2 cob 1 cob
+ + + + = (13)
Al aumentar la carga la suma de las prdidas p
ac
+p
mec
disminuye algo a causa de la
disminucin del flujo principal, as como de la disminucin de la velocidad de rotacin,
generalmente estas prdidas se incluyen en las prdidas constantes del motor. En
contradiccin con las prdidas en el acero, las prdidas en el cobre varan
proporcionalmente al cuadrado de la corriente.
D. Factor de potencia ). ( cos
2
P f = El motor asincrnico, consume de la red corriente
retrasada que casi no depende de la carga. Por eso su cos es siempre menor de la unidad.
En marcha en vaco cos habitualmente no sobrepasa de 0,2, pero luego, bajo carga, crece
bastante rpidamente y alcanza el mximo cuando la potencia es prxima a la nominal. Al
seguir aumentando la carga disminuye la velocidad del motor, correspondientemente a lo
cual aumenta el ngulo
2 2 2
/ r sx arctg = y disminuyen el cos
2
y el cos.
E. Capacidad de sobrecarga del motor. Se llama capacidad de sobrecarga del motor
asincrnico o, de otra manera, momento de vuelco, a la relacin del momento mximo del
motor a su momento nominal, es decir,
. M / M k
n mx m
=
Habitualmente en los motores de pequea y media potencia k
m
=1.6-1.8. En los motores de
media y gran potencia k
m
=1.8-2.5. en los motores especiales k
m
alcanza 2,8-3,0 y ms.
Influencia de la asimetra
La diferencia de las corrientes en las fases, cuando se trata de un devanado con simetra de
impedancia en el estator y rotor, se manifiesta debido a la asimetra de tensiones primarias en
la red. Si el punto neutro del devanado del estator, conectado en estrella ests aislado,
entonces las componentes de la secuencia cero en las corrientes y las tensiones no surgen. Las
corrientes en el devanado del estator se representa como la suma de las corrientes de
secuencias directas e inversas.
La influencia de la asimetra de tensin se tiene en cuenta en la caracterstica mecnica y de
funcionamiento a travs de las componentes simtricas, analizado la secuencia negativa y
positiva por separado y luego se unen para dar el resultado final.
Debido a que las impedancias de secuencia positiva son mayores que para secuencia negativa,
la distorsin de la simetra de corriente supera en muchas veces la distorsin de asimetra de
tensin.
En relacin con esto se plantean exigencias altas con respecto a la simetra de las tensiones de
la redes que alimentan las mquinas asincrnicas.
Influencia del efecto pelicular.
Al circular corriente por los devanados es necesario tener en cuenta que debido al efecto
pelicular las densidades de corrientes se distribuyen irregularmente por la ranura. Debido a
esto las componentes activa y reactiva de la impedancia varan con un coeficiente kr y kx
respectivamente.
En la bibliografa se destaca que ambos coeficientes dependen de la altura reducida del
conductor elemental ,
* h =
definida como la relacin de la altura del conductor (h) a la profundidad de penetracin del
campo electromagntico ().
r t
1 o
b
b f

=
donde b
r
y b
1
es la ancho de la ranura y ancho de los conductores en la ranura,
respectivamente, f frecuencia del circuito y
t
la resistencia especfica a la temperatura t del
material que se trate.
Para tener en cuenta lo anterior se dedujeron ecuaciones de la altura del conductor elemental
en funcin de parmetros del motor segn clase de diseo, los cuales sirven para determinar
los coeficientes para la resistencia y reactancias a travs de expresiones empricas existente
en las bibliografas analizadas.
Aplicacin FUNCMOTOR
La nomenclatura FUNCMOTOR ha sido designada para identificar al motor de induccin de
rotor en jaula en la aplicacin. Como todo algoritmo de interaccin se divide en tres partes
principales: entrada de datos y seleccin de tareas; clculo y por ltimo resultados, estos
pueden ser en forma numrica o en forma grfica. Se ha realizado con ayuda de la Interfase
Grfica del Usuario (GUI) del Matlab con visualizacin de datos a travs de ventanas. A
continuacin se ofrece una panormica de cmo se estructura.
Entrada de datos.
Lo primero que se introduce son los datos nominales del motor objeto de anlisis, los cuales
servirn para el trabajo futuro en determinada tareas, constituido por: Tensin, Potencia,
Rendimiento, Velocidad, Frecuencia, Factor de potencia y Momento de inercia, adems de
seleccionar la conexin y clase de diseo del motor.
Como se observa en la figura 2, los datos se solicitan a travs de una ventana, las cuales
tienen valores por defecto que sirven para una demostracin.

Fig. 2. Muestra de la forma de entrada de datos nominales y opciones de clculo.
Al seleccionar una de las vas fundamentales aparecer la ventana de entrada de datos
especfico, como se muestra en la figura 3 para la cuarta opcin.
Existe una interrelacin entre las distintas opciones, con el objetivo de visualizar los datos
solicitados y decidirse por una de ellas, se puede pasar de una seleccin a otra y solo se hace
efectiva al pulsar el botn Siguiente.
Luego aparece la figura 4, sin las grficas que se muestran, solo se observan las opciones
generales de Fun y Perd y Mec y Corr, la primera permite visualizar las caractersticas de
funcionamiento y prdidas y la segunda visualiza las caractersticas de corriente y mecnica.
Fig. 3. Ventana para la introduccin de datos segn la opcin seleccionada.

Fig. 4. Entrada de frecuencia y tensin diferentes a los nominales. Resultados grficos
Al ejecutar unos de los botones anteriores aparecen las opciones de Cerrar e Iniciar con lo
cual se puede hacer valer su significado. En el men de opciones se puede tener en cuenta o
no el efecto pelicular en las caractersticas.
Los valores reales de operacin de tensin y frecuencia, se toman por defecto igual a los
nominales, figura. 4, si se quiere ver el comportamiento en rgimen asimtrico, basta con
cambiar los valores por defecto para otros valores diferentes.
Otras opciones interna de la aplicacin.
Durante la ejecucin de FUNCMOTOR se salvan datos en cuatro ficheros que sirven para
disponer de un conjunto de datos en forma numrica. Se pueden ejecutar con el comando
load del MATLAB.
1. datosnominales : Guarda en su interior los datos nominales fundamentales, como se
relacionan a continuacin y en el mismo orden.
Magnitud Descripcin de la magnitud Magnitud Descripcin de la magnitud
Un Tensin en V Us1 Tensin de fase en V
Pn Potencia en W Sn Deslizamiento en v.r
Rn Rendimiento en v.r p Nmero de pares de polos
fn Frecuencia en Hz n1 Velocidad sincrnica en rev/min.
Fpn Factor de potencia en v.r In Corriente nominal en A
Sumperd Prdidas nominales en W nr Velocidad del rotor en rev/min.
Mn Momento en N.m Qn Potencia reactiva en VAr
Ta Temperatura ambiente
o
C
2. datoscircuito : Guarda los datos del circuito equivalente en L.
Mag. Descripcin de la magnitud Mag. Descripcin de la magnitud
Rcs Resistencia del estator xm Reactancia de magnetizacin
Rcr Resistencia del rotor Sm Deslizamiento crtico en v.r
Xs Reactancia del estator rcm Resistencia de magnetizacin
Xr Reactancia del rotor Todos a temperatura ambiente
3. servicio : Almacena datos para el trazado de las caractersticas de servicio.
Mag. Descripcin de la magnitud Mag. Descripcin de la magnitud
P11s Potencia de entrada del motor en W S Deslizamiento
Memt Momento electromagntico en N.m CosFi Factor de potencia de entrada v.r
Iestator Corriente del estator (Valor absoluto) A Re Rendimiento en v.r.
Nr Velocidad del rotor rev/min Pss Potencia de salida en W
ptotal Suma de prdidas totales
4. mecanica : Almacena los datos fundamentales para el trazado de las caractersticas
mecnicas y de corriente. Tiene en cuenta el efecto pelicular y la asimetra de tensiones y
corriente.
Mag. Descripcin de la magnitud Mag Descripcin de la magnitud
Sk Deslizamiento en v.r. I12n Corriente de secuencia negativa A.
Memp Momento elctromagntico sec (+) N.m. Rrp Resistencia del rotor secuencia (+)
Memn Momento electromagntico sec (-) N.m Rrn Resistencia del rotor secuencia (-)
Memr Momento eletromagntico total N.m. Xrp Reactancia del rotor secuencia (+)
I11p Corriente de secuencia positiva A Xrn Reactancia del rotor secuencia (-)
Iestm Corriente media del estator en A
Resultados generales que se pueden obtener del uso de FUNCMOTOR
1. Obtener los datos del circuito equivalente por diferentes vas.
2. Trazar las caractersticas de funcionamiento o de servicio.
3. Trazar las caractersticas de prdidas.
4. Trazar las caractersticas de corriente y mecnica.
5. Obtener los datos teniendo en cuenta el efecto pelicular.
6. Obtener los datos teniendo en cuenta la asimetra de tensin y corriente.
7. Obtener los datos teniendo en cuenta la asimetra y efecto pelicular.
Conclusiones
Todo esto resuelve la problemtica de las ecuaciones ms generales que aparecen en la
mayor parte de la literatura, que se cumplen con gran exactitud en los valores cercanos al
deslizamiento nominal y que carecen de la misma cuando se acercan al punto de arranque
y freno electromagntico.
Esta aplicacin sirve como un medio de anlisis del motor de induccin sujeto a estas
particularidades.
Tiene uso de carcter docente para demostrar y actualizar las teoras que se imparte al
estudiantado.
Tiene carcter cientfico en cuanto a la tcnica utilizada en el algoritmo y los resultados que
arroja, que pueden ser objeto de discusin.
Tiene uso prctico, por cuanto se refiere a una problemtica de los talleres de reparaciones
y pruebas a motores de induccin.
Se observa con claridad como afectan la asimetra y efecto pelicular en la caracterstica
fundamental de un Motor de induccin.
Bibliografa
(1). Mc. Pherson Robert D Laramore. A introduction to electrical machines and transformers.
Second edition. 1992.
(2). Kostenko Piotrosvki, Mquinas Elctricas. Tomo II Editorial Mir. Mosc. (20.1-20.10 y
20.18), ejemplos resueltos.
(3). A. V. Ivanov-Smolenski. Mquinas Elctricas.Tomos I y II. Editorial Mir. Mosc. 1984.
(4). Matlab6.0 GUI
Anexos
0 0.2 0.4 0.6 0.8 1
0
50
100
150
200
Normal y E. Pelicular
S
0 0.5 1 1.5 2
-50
0
50
100
150
200
Normal y asimetra
S
N.m
N.m
S
S
1
2
3
4
5
Fig.A1. Caractersticas mecnicas para diferentes situaciones: 1 Normal, 2 considerando el efecto pelicular, 3
Momento resultante con asimetra leve sin efecto pelicular, 4 Momento de sec +, 5 Momento de sec ( -).
0 0.5 1 1.5 2
-100
-50
0
50
100
E. Pelicular y Asimetra
S
0 0.5 1 1.5 2
0
50
100
150
200
250
Normal, E. Pelicular y asimetra
S
N.m
(A)
S
S
1
2
3
1
2
3
4
Fig.A2. Caractersticas mecnicas: 1 Momento de sec +, 2 Momento de sec (-), 3 Momento resultante con alto nivel
de asimetra y con efecto pelicular; Caracterstica de corriente: 1 corriente normal, 2 corriente con efecto
pelicular, 3 con asimetra baja y 4 corriente con alta asimetra.
0 0.2 0.4 0.6 0.8 1 1.2
0
0.2
0.4
0.6
0.8
1
0 0.2 0.4 0.6 0.8 1 1.2
0
0.5
1
1.5
2
2.5
3
P2/Pn
Re Spr/Spn
P2/Pn
1
2
3
4
Fig.A3. Rendimiento: 1 normal, 2 con asimetra; prdidas totales: 3 normal, 4 con asimetra.