SIMULACION DE LA MAQUINA SINCRONICA
Alfredo Rifaldi
1 - Motivación
Uno de los temas de interés es el funcionamiento en cortocircuito de la maquina sincrónica, esta se encuentra en marcha a la velocidad sincrónica y ocurre una falla en su proximidad.
Se pretende encarar este tema con la teoría unificada de las maquinas rotantes, utilizando el esquema de la maquina universal, que permite representar correctamente a cualquier maquina, y en particular la que nos interesa en este trabajo
2 - La maquina sincrónica
Los estudios de Kron y Park sobre la maquina sincrónica fueron los primeros ejemplos de tratamiento de las maquinas trifásicas con la teoría de los dos ejes.
La maquina se describe con tres matrices, resistencias
| R |
0 |
0 |
0 |
0 |
0 |
| |
Directo |
| 0 |
R |
| |
Cuadratura |
||||
| 0 |
R |
| |
Homopolar |
||||
| 0 |
Re |
| |
Excitación |
||||
| 0 |
Rs |
| |
Amortiguador directo |
||||
| 0 |
Rt |
| |
Amortiguador cuadratura |
Inductancias de transformación
| Ld |
0 |
0 |
Mde |
Mds |
0 |
| |
| 0 |
Lq |
Mqt |
| |
|||
| 0 |
Lgamma |
| |
||||
| Mde |
Le |
Mes |
| |
|||
| Mds |
Mes |
Ls |
| |
|||
| 0 |
Mqt |
Lt |
| |
Inductancias rotacionales
| 0 |
-Lq |
0 |
0 |
0 |
-Mqt |
| |
| Ld |
Mde |
Mds |
| |
|||
| 0 |
| |
|||||
| 0 |
| |
|||||
| 0 |
| |
|||||
| 0 |
| |
Las ecuaciones de funcionamiento eléctrico y mecánico son las siguientes
| RLG |
| = |
| R |
| + |
| L |
| * p + omegar |
| G |
| |
| Vd |
| = |
| R LG |
| * |
| id |
| |
| Vq |
| |
| |
| |
| iq |
| |
| Vgamma |
| |
| |
| |
| igamma |
| |
| Ve |
| |
| |
| |
| ie |
| |
| 0 |
| |
| |
| |
| is |
| |
| 0 |
| |
| |
| |
| it |
| |
Cupla eléctrica
Ce / pares de polos |
= |
[ id iq igamma ie is it ] * |
| G |
| * |
| id |
| |
| |
| |
| iq |
| |
|||
| |
| |
| igamma |
| |
|||
| |
| |
| ie |
| |
|||
| |
| |
| is |
| |
|||
| |
| |
| it |
| |
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 * |
| - flujoq |
| |
| flujod |
| |
|
| 0 |
| |
|
| 0 |
| |
|
| 0 |
| |
|
| 0 |
| |
Ce / pares de polos = flujod * iq - flujoq * id
3 – El cortocircuito
Ya hemos visto como se plantean las ecuaciones de funcionamiento de la maquina sincrónica trifásica, la maquina se encuentra en marcha. Las ecuaciones diferenciales que describen y rigen su comportamiento son:
[ V ] = [ R ] * [ I ] + [ L ] * delta [ I ] / deltaT + omega * [ G ] * [ I ]
No olvidemos que se trata de un sistema de varias ecuaciones diferenciales, del tamaño de las matrices, 4, 5 o 6, de estas se logra despejar delta [ I ]
.delta [ I ] = inversa [ L ] * ( [ V ] – [ R + omega * G ] * [ I ] ) * deltaT
El sistema de ecuaciones diferenciales se puede resolver por métodos numéricos, el programa SINPET resuelve por Runge Kutta.
Para el lote de datos que se detalla a continuación
POTENCIA NOMINAL |
20 |
KVA |
TENSION DE FASE |
200 |
V |
CORRIENTE NOMINAL |
58 |
A |
TENSION DE CAMPO |
110 |
V |
VELOCIDAD NOMINAL |
1000 |
VPM |
FRECUENCIA NOMINAL |
50 |
HZ |
RESISTENCIA ESTATOR |
0.1 |
OHM |
RESISTENCIA ROTOR |
32.4 |
OHM |
INDUCTANCIA ARMADURA D |
4.777 |
MHY |
INDUCTANCIA ARMADURA Q |
2.866 |
MHY |
INDUCTANCIA DE CAMPO |
9300 |
MHY |
MUTUA CAMPO ARMADURA |
187 |
MHY |
RESISTENCIA AMORTIGUADOR D |
0.186 |
OHM |
RESISTENCIA AMORTIGUADOR Q |
0.605 |
OHM |
INDUCTANCIA AMORTIGUADOR D |
4.777 |
MHY |
INDUCTANCIA AMORTIGUADOR Q |
2.866 |
MHY |
MUTUA AMORTIGUADOR D-ARMADURA d |
4.29 |
MHY |
MUTUA AMORTIGUADOR Q-ARMADURA q |
1.952 |
MHY |
MUTUA AMORTIGUADOR D-CAMPO |
195 |
MHY |
INDUCTANCIA SUBTRANSITORIA D |
0.819 |
MHY |
INDUCTANCIA SUBTRANSITORIA Q |
1.537 |
MHY |
INDUCTANCIA TRANSITORIA D |
1.017 |
MHY |
La ejecución del programa arroja los resultados con los que se han desarrollado los siguientes gráficos:
Figura 4 – corriente de las tres fases Figura 3 – corriente de una fase Figura 2 – variación de la corriente de excitación, desde que inicia el cortocircuito hasta que vuelve el valor inicial (la tensión de excitación se mantiene constante). Figura 1 – variación de la corriente en el circuito amortiguador (según el eje directo) Figura 5 – variación del par durante el cortocircuito (con velocidad constante)Se dispone de un modelo de la maquina sincrónica que permite observar la corriente de cortocircuito que se presenta en el estator, y también en el devanado de excitación, y en las jaulas amortiguadores.
4 – Notas sobre el programa SINPET
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 SINPET.TXT (ver archivo) se trata de ocho registros de datos, con mas o menos valores en cada registro, entre cada par de valores un espacio en blanco de separación.Entre estos registros se pueden tener registros de comentarios (renglones que inician con C o con *) utiles para la preparacion de datos y que en la ejecucion del programa este descartan.
Tratemos de aclarar un poco mas el significado de los datos que requieren estos cálculos, los datos deben permitir armar las matrices arriba citadas, las resistencias en ohm son fáciles de comprender,
R resistencia estator,
Re resistencia rotor,
Rs y Rt resistencias de amortiguador directa y en cuadratura,
cuya evaluación requiere conocimientos de cómo esta construida la maquina.
Las inductancias en milihenry,
Ld y Lq inductancia de armadura directa y en cuadratura,
Le inductancia de campo,
Ls y Lt inductancias de amortiguador directa y en cuadratura, valores mas difíciles de evaluar.
Por ultimo las inductancias mutuas:
Mde campo – armadura,
Mds amortiguador d - armadura d,
Mqt amortiguador q - armadura q,
Mes amortiguador d – campo.
Estos valores permiten determinar las inductancias (reactancias) subtransitorias y transitorias, que son los datos clásicos utilizados para la determinación de las corrientes de cortocircuito.
Ld’’ = Ld – (Le * Mds^2 + Ls * Mde^2 - 2. * Mde * Mds * Mes) / (Le * Ls – Mes^2)
Lq’’ = Lq – Mqt^2 / Lt
Ld’ = Ld – Mde^2 / Le
Esto es lo que el programa hace antes de iniciar los cálculos de resolución de las ecuaciones diferenciales, los valores de reactancias fuera de rangos lógicos, pueden descalificar los datos introducidos.
El programa genera dos archivos de resultados, el primero separado por punto y coma ";" (con extensión CSV) 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.
El segundo archivo muestra las matrices de datos, útiles para control u otros cálculos.
- Bibliografía consultada