SIMULACION DE LA MAQUINA ASINCRONICA
Alfredo Rifaldi
1 - Motivación
Interesa observar el funcionamiento del motor asincrónico, consideramos que este arranca, o se encuentra en marcha y ocurren variaciones de la tensión de alimentación.
Se pretende encarar este tema con la teoría unificada de las maquinas rotantes, utilizando el esquema de la maquina universal, que permite representar cualquier maquina, y en particular la que nos interesa en este trabajo
2 - La maquina asincrónica
La maquina se describe con tres matrices, resistencias
| R1 |
0 |
0 |
0 |
0 |
0 |
| |
Directo |
| 0 |
R1 |
| |
Cuadratura |
||||
| 0 |
Ra |
Rc |
| |
Jaula externa directa |
|||
| 0 |
Ra |
Rc |
| |
Cuadratura |
|||
| 0 |
Rc |
Rb |
| |
Jaula interna directa |
|||
| 0 |
Rc |
Rb |
| |
Cuadratura |
Inductancias de transformación
| L1 |
0 |
M1a |
0 |
M1b |
0 |
| |
| 0 |
L1 |
M1a |
M1b |
| |
||
| M1a |
La |
Mab |
| |
|||
| 0 |
M1a |
La |
Mab |
| |
||
| M1b |
Mab |
Lb |
| |
|||
| 0 |
M1b |
Mab |
Lb |
| |
Inductancias rotacionales
| 0 |
0 |
0 |
0 |
0 |
0 |
| |
| 0 |
| |
|||||
| 0 |
-M1a |
-La |
-Mab |
| |
||
| M1a |
La |
Mab |
| |
|||
| 0 |
-M1b |
-Mab |
-Lb |
| |
||
| M1b |
Mab |
Lb |
| |
Las ecuaciones de funcionamiento eléctrico y mecánico son las siguientes
| RLG |
| = |
| R |
| + |
| L |
| * p + omegar |
| G |
| |
| Vd |
| = |
| R LG |
| * |
| id |
| |
| Vq |
| |
| |
| |
| iq |
| |
| 0 |
| |
| |
| |
| iad |
| |
| 0 |
| |
| |
| |
| iaq |
| |
| 0 |
| |
| |
| |
| ibd |
| |
| 0 |
| |
| |
| |
| ibq |
| |
Cupla eléctrica
Ce / pares de polos |
= |
[ id iq iad iaq ibd ibq ] * |
| G |
| * |
| id |
| |
| |
| |
| iq |
| |
|||
| |
| |
| iad |
| |
|||
| |
| |
| iaq |
| |
|||
| |
| |
| ibd |
| |
|||
| |
| |
| ibq |
| |
Ecuación de equilibrio mecánico
Ce + Cr = J p omegam + D omegam + De (omegam - omega / pares de polos)
Durante el funcionamiento como motor Ce cupla (par) electromagnética es positiva, durante el funcionamiento como generador Cr será la cupla motriz, positiva.
Los flujos concatenados están relacionados con inductancias y corrientes
[ flujos ] = [ L ] * [ i ]
Las tensiones pueden relacionarse con los flujos
[ V ] = [ R ] * [ i ] + p * [ flujos ] + omegar * [ G ] * [ i ] |
3 - Programa de calculo
Con esta metodología, se desarrollo el programa MOTPET2 que permite resolver el problema que plantea un motor asincrónico de jaula simple o doble. en funcionamiento transitorio, se preparan los datos en un archivo respetando algunas reglas.
El lote de datos (ver MOTPET2.DAT) incluye comentarios (renglones que inician con * o C), y que ayudan a prepara los datos pero se descartan al ejecutar el programa.
De la ejecución se obtiene un archivo de resultados aptos para la graficación de distintos valores, los gráficos pueden verse con EXCEL.
Es interesante observar los picos de la cupla eléctrica en los primeros instantes del arranque, fenómeno que pasa desapercibido cuando se resuelve el problema con métodos mas simplificados.
Mientras se presentan los picos de cupla, se observan también variaciones de velocidad, oscilaciones alrededor del valor creciente que caracteriza la aceleración inicial.
Variando los valores (amplitudes) de vd y vq se pueden simular distintas condiciones, por ejemplo:
Los ejemplos incluidos corresponden a un motor simple jaula, la representación de un motor de doble jaula es totalmente análoga, los arrollamientos deben considerarse representando en el rotor una y otra jaula (jaulas superficial y profunda respectivamente - también llamadas de arranque y marcha).
La matriz que representa los arrollamientos es de 6 x 6 (4x4 para jaula simple) observamos que trabajando con matrices la mayor complicación del problema prácticamente no se nota.
El motor de jaula profunda se puede representar en forma análoga con mas arrollamientos, la matriz será de tamaño n x n.
Otra forma de encarar el tema de estas representaciones es mantener el tamaño de la matriz en 4 x 4 y variar los parámetros que correspondan en función de la velocidad, así también puede resolverse la inserción de una resistencia de arranque en un motor con rotor bobinado.
Por la simetría que presentan los motores asincrónicos los parámetros directos y en cuadratura son iguales, para el motor de simple jaula la
figura 21 muestra la relación entre los parámetros del modelo de la teoría clásica, y el que aquí presentamos.Análogamente se relacionan los parámetros del motor doble jaula, obsérvese la
figura 22.Hasta este punto se ha considerado que las tensiones vd y vq son las que alimentan el motor, en rigor la fuente real tendrá también su impedancia, en la que se presenta una caída de tensión que depende de la corriente que el motor absorbe.
4 - ejemplos de resultados
De la bibliografía (Petrecca) fueron extraídos los ejemplos que se reproducen y sirven como comprobación del buen desempaño del programa. La tabla muestra los datos que corresponden a los ejemplos que se presentan:
SIMULACION DE MOTOR ASINCRONICO |
||
PASO DE TIEMPO |
0.001 |
SEGUNDOS |
CANTIDAD DE PUNTOS |
500 |
|
TIEMPO DE SIMULACION |
0.5 |
SEGUNDOS |
ANGULO DE INSERCION |
0 |
GRADOS |
POTENCIA NOMINAL |
4 |
KVA |
TENSION DE RED |
660 |
V |
CORRIENTE NOMINAL |
0 |
A |
TENSION DE FASE |
380 |
V |
VELOCIDAD NOMINAL |
1440 |
VPM |
FRECUENCIA NOMINAL |
50 |
HZ |
PARES DE POLOS |
2 |
- |
TRANSITORIO ELECTOMAGNETICO DE UN MOTOR ASINCRONICO |
||
motor representado como maquina universal |
||
se trata de rotor simple jaula o bobinado |
||
R1 resistencia estator |
1.45 |
|
RA resistencia rotor |
1.18 |
|
HL1 inductancia estator |
0.15088 |
|
HLA inductancia rotor |
0.15088 |
|
HM1A mutua rotor estator |
0.14324 |
|
WJ momento de inercia: |
0.1 |
|
Cupla antagonica = AW |
0 |
|
velocidad (radianes/seg) * BW |
0 |
|
velocidad ^ 2 * CW |
0 |
|
velocidad ^ 3 * DW |
0 |
|
variaciones de tensión V1 T1 |
1 |
0 |
V2 T2 |
0 |
************ |
V3 T3 |
0 |
************ |
Se observa la evolución de la velocidad
figura 1 el par motor figura 2 siendo nulo el par antagónico y la corriente de la fase A figura 3.Si el momento de inercia es menor, WJ = 0.02 el tiempo de aceleración es menor, obsérvese la velocidad
figura 4 el par motor figura 5 siendo también nulo el par antagónico y la corriente de la fase A figura 6.Se vuelve a fijar el momento de inercia WJ = 0.1, y par antagónico
Cupla antagonica = AW |
0.7 |
velocidad (radianes/seg) * BW |
0 |
velocidad ^ 2 * CW |
0.001 |
velocidad ^ 3 * DW |
0 |
Se observa la velocidad
figura 7 los pares motor y antagónico figura 8 y la corriente de fase A figura 9El ultimo ejemplo se desarrolla nuevamente con el motor en vacío, se lo arranca a tensión nominal, luego se reduce la tensión al 50% 0.3 segundos después se vuelve a aplicar tensión 100%
Variaciones de tensión V1 T1 |
1 |
0 |
V2 T2 |
0.5 |
1.1 |
V3 T3 |
1 |
1.4 |
obsérvese la velocidad
figura 10, el par motor figura 11 (el antagónico es nulo) y por ultimo la corriente de la fase A figura 125 - Notas sobre el programa MOTPET2
Los datos se preparan en un archivo de texto en formato libre, que incluye cierta cantidad de registros, si todos los registros son nulos el programa contiene un ejemplo que se ejecuta.
Las reglas de preparación de datos son las incluidas en el
archivo MOTPET2.TXT (ver archivo) se trata de nueve registros de datos, con mas o menos valores en cada registro, entre cada par de valores un espacio en blanco de separación.Tratemos de aclarar un poco mas el significado de los datos que requieren estos cálculos,
motor representado como maquina universal |
se trata de rotor doble jaula (o profunda) |
R1 resistencia estator |
RA resistencia jaula interna |
RB resistencia jaula externa |
RC resistencia común jaulas |
HL1 inductancia estator |
HLA inductancia jaula interna |
HM1A mutua interna estator |
HLB inductancia jaula externa |
HM1B mutua externa estator |
HMAB mutua jaulas |
WJ momento de inercia: |
Cupla antagónica = AW |
velocidad (radianes/seg) * BW |
velocidad ^ 2 * CW |
velocidad ^ 3 * DW |
Si RB es igual o mayor de 999 se considera que el motor es de simple jaula, y ciertos datos no se utilizan.
El programa genera un archivo de resultados separados por punto y coma ";" (con extensión CSV) que permite ser levantado con Excel, separar los resultados en columnas, y generar los gráficos temporales.
Las figuras que se observan en el punto anterior, se han generado en esta forma.
Un segundo archivo muestra las matrices de datos, útiles para control u otros cálculos.
6 - Datos - relación entre modelos
Veamos nuevamente la relación entre modelos tradicional y de la maquina generalizada, observemos nuevamente la
figura 21, consideramos la situación de funcionamiento en régimen, las ecuaciones que corresponden al modelo que incluye el transformador son:V1 = (R1 + j omega * L1) * I1 + j omega * M * I2
0 = j omega * M * I1 + (R2 / s + j omega * L2) * I2
Las ecuaciones corresponden a la figura a), si se modifican como sigue corresponderán a la figura b).
V1 = (R1 + j omega * (L1 - M)) * I1 + j omega * M * (I2 + I1)
0 = j omega * M * (I1 + I2) + (R2 / s + j omega * (L2 - M)) * I2
La
figura 23 (obtenida a partir de la figura 21) muestra el modelo tradicional (en régimen permanente) esta representado por la misma figura b), cambiando las variables figura c), y el sentido de las corrientes.Ip2 = - I2 * (N2 / N1)
los parámetros indicados en la figura b) corresponden al modelo de la maquina generalizada, la figura c) al modelo tradicional, consideramos la situación de funcionamiento en régimen, las ecuaciones que corresponden al modelo que incluye el transformador
V1 = (R1 + j X1) * I1 + j Xpm * (I1 - Ip2)
0 = j Xpm * (I1 - Ip2) - (Rp2 / s + j Xp2) * Ip2
nótese que para estas ecuaciones la corriente cambio de sentido:
En el calculo de los motores asincronicos, o en el ensayo de los mismos, los valores que se obtienen son estos últimos, entonces a partir de estos valores (todos referidos al estator) se debe llegar a los parámetros de la maquina generalizada:
Xpm = omega * M * (N1 / N2)
Xp2 = omega * L2 * (N1 / N2)^2 - Xpm
Rp2 = R2 * (N1 / N2)^2
X1 = omega * L1 - Xpm
Se conocen o se han calculado X1, Xp2, Xpm, R1, Rp2, la relación N1 / N2 es conocida, o arbitrariamente se hace igual a 1, y se pueden despejar los valores que son de interés para realizar el modelo de la maquina generalizada:
R1 = R1
L1 = (X1 + Xpm) / omega
M = Xpm * (N2 / N1) / omega
R2 = Rp2 * (N2 / N1)^2 / omega
L2 = (Xp2 + Xpm) * (N2 / N1)^2 / omega
I2 = - Ip2 * (N1 / N2)
Con estos valores obtenidos a partir del modelo tradicional se dispone de los datos para trabajar con el modelo de la maquina generalizada, con la que se representa un motor trifasico asincronico de simple jaula o de rotor bobinado.
Si: N1 / N2 = 1, se tiene:
R1 = R1
L1 = (X1 + Xpm) / omega
M = Xpm / omega
R2 = Rp2 / omega
L2 = (Xp2 + Xpm) / omega
I2 = - Ip2
Análogamente a partir de la
figura 24 que corresponde al modelo de una maquina doble jaula (u otra representada en dicha forma) se preparan los datos del modelo generalizado (figura superior) partiendo del modelo tradicional de una maquina asincrónica (doble jaula - figura inferior).En esta caso hacemos N1 / N2 = 1 por simplicidad, es inmediato que las resistencias de ambos modelos son iguales, con las reactancias se tiene:
X1 = omega * (L1 - M)
Xm = omega * M
Xa = omega * (La - Mab)
Xb = omega * (Mab – M)
0 = omega * (Lb – Mab)
De estas surge
L1 = (X1 + Xm) / omega
M = Xm / omega
Mab = (Xb + Xm) / omega
La = (Xa + Xm) / omega
Lb = Mab
Asi se puede pasar de un modelo a otro, generalmente se conocen datos (de calculo o de ensayo) que corresponden al modelo tradicional, y en esta forma se obtienen los valores con los que se carga el modelo de maquina universal (generalizada) que representa el motor asincronico.
7 - Bibliografía consultada