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 9

El 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 12

5 - 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