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