UNI FIEE SECCIÓN POSGRADO

EE55 MODELOS DE MÁQUINAS ELECTRICAS EN SEP

PRIMER TRABAJO 2022-I

Roberto Ramírez A.

INTEGRANTES
Nombre Mijhael P. Tuesta Herrera
Fecha 13-06-2022
 TABLA DE CONTENIDO

1. INTRODUCCIÓN ................................................................................................... 5
2. OBJETIVO ............................................................................................................. 5
3. PARAMETROS DE LA MÁQUINA SINCRONA A CONSIDERAR .......................... 5
4. PROBLEMAS A RESOLVER ................................................................................. 5
4.1 CASO1: MÁQUINA SINCRONA TRIFASICA EN VACÍO ........................................ 6
4.2 CASO2: MÁQUINA SINCRONA TRIFASICA EN VACIO Y EL ESTATOR EN
CORTOCIRCUITO .......................................................................................................... 12
5. CONCLUSIONES ................................................................................................. 18
 ÍNDICE DE TABLAS

Tabla 1. Parámetros eléctricos ...................................................................... 5

Tabla 2. Ecuaciones de la MS caso 1 ........................................................... 7
Tabla 3. Ecuaciones de la MS caso 2 ......................................................... 13
 ÍNDICE DE FIGURAS

Figura 1. Transformada de Clarke y Parck ........................................................... 5

Figura 2. Máquina sincrona trifasica .................................................................... 6
Figura 3. Máquina sincrona trifasica en vacío ...................................................... 6
Figura 4.Ecuaciones MS en los ejes d – q ........................................................... 7
Figura 5.Modelado en SIMULINK ........................................................................ 7
Figura 6.Solución de ecuación (3) y (4) en Live Editor ......................................... 8
Figura 7.Transformada a sistemas 𝛼𝛽 y d-q ........................................................ 8
Figura 8.Diagrama de bloques Problema N° 1 ..................................................... 8
Figura 9.Corriente de Campo If(t) y del devanado D ............................................ 9
Figura 10.Tensión Vd y Vq ................................................................................ 10
Figura 11.Tensión 𝑽𝜶 y 𝑽𝜷................................................................................ 10
Figura 12.Tensión 𝑽𝒂𝒃𝒄 .................................................................................... 11
Figura 13.Diagrama de Bloques con bobina “Q” estuviera en circuito abierto .... 12
Figura 14. IQ con bobina “Q” estuviera en circuito abierto ................................. 12
Figura 15. Máquina sincrona trifasica en vacío y estator en cortocircuito ........... 12
Figura 16. Diagrama de bloques Problema N° 2 ................................................ 13
Figura 17. Corriente de Campo If y devanado ID-IQ .......................................... 14
Figura 18. Corriente de Id y Iq ........................................................................... 15
Figura 19. Corrientes Iα - Iβ ............................................................................... 16
Figura 20. Corriente de abc ............................................................................... 17
Figura 21. Corriente de If(t) – Id con y sin devanado Amortiguador ................... 18
1. INTRODUCCIÓN
La Máquina Sincrona son máquinas capaces de operar sólo a la velocidad sincrónica
(Wo) y son usadas en centrales para la generación de energía eléctrica. Para el estudio
dentro del sistemas de potencia se usan modelos para estudiar el comportamiento
estático y dinámico. Debido a u naturaleza las ecuaciones de máquina síncrona tiene
características no lineales y hacen muy complejo su estudio. Por tanto, se usan
transformaciones que permiten convetir estas ecuaciones trifásicas en ecuaciones mas
simples para su tratamiento en el modelamiento de la máquina sincrona.
Las relaciones importantes son la transformada de Clarke y Parck considerando que la
potencia es invariente, defina como:

Figura 1. Transformada de Clarke y Parck

2. OBJETIVO
Determinar y graficar la evolución en el tiempo de la corriente y tensión de la máquina
síncrona (MS) trifásica usando las ecuaciones eléctricas usando el modelo d-q y la
transformada de Clarke y Park para obtener valores trifásicos.

3. PARAMETROS DE LA MÁQUINA SINCRONA A CONSIDERAR

Se indica la siguiente información para resolver y determinar la evolución en el tiempo
de la corriente y tensión de la máquina síncrona propuesta.
Tabla 1. Parámetros eléctricos

Lad Laq RD RQ Wr(*)

Vf Ra (Ω) Ld (H) Lq (H) R f (Ω) L f (H) LD (H) LQ (H)
(H) (H) (Ω) (Ω) (rad/s)
50 4.70 2.70 2.40 1.60 1.30 0.56 2.80 10.00 2.59 2.40 1.40 157.08
(*) Wr = 1500 RPM

4. PROBLEMAS A RESOLVER
Se propone una máquina síncrona que se muestra en la Figura 2, considerando lo
siguiente:
• Caso1: Considerar la máquina síncrona girando en vacio.
• Caso2: Considerar la máquina síncrona con el estator en cortocircuito.
 Figura 2. Máquina sincrona trifasica

4.1 CASO1: MÁQUINA SINCRONA TRIFASICA EN VACÍO

Se considera la máquina síncrona de la Figura 3, utilizando el modelo d-q, para
resolver el caso 1 se tiene en cuenta las siguientes condiciones:
• Máquina síncrona se impulsa en vacio
• En el tiempo t=0 se aplica al devanado “f” una tensión continua Vf.
• El devanado de compensación en cortocircuito ( 𝑉𝐷 = 𝑉𝑄 = 0)

Figura 3. Máquina sincrona trifasica en vacío

Sobre las condiciones iniciales, se puede deducir que la corriente en el estator es cero
y usando las relaciones de Clarke y Park se determina que la corriente del estator en
los ejes 𝛼 y 𝛽 como en d – q es igual a cero.
𝑰𝒂 = 𝑰𝒃 = 𝑰𝒄 = 𝑰𝜶 = 𝑰𝜷 = 𝑰𝒅 = 𝑰𝒒 = 𝟎 (1)
Asimismo, considerano las ecuaciones de la máquina síncrona reflejado en los ejes d-
q, se puede dedudir lo siguiente:
𝑰𝑸 = 𝟎 (2)
 Figura 4.Ecuaciones MS en los ejes d – q
Sobre las consideraciones iniciales tenemos las siguientes ecuaciones, las mismas
que se puede implementar en SIMULINK que es un toolbox de MATLAB.
Tabla 2. Ecuaciones de la MS caso 1

4.1.1 Implementación en SIMULINK del toolbox MATLAB

Antes de mostrar la implementación completa, se explica la estrategia usada para
la implementación de las ecuaciones en SIMULINK de todo este informe.
Considerando Fixed Step con un paso de 0.001.

𝑉𝑓 = 𝑟𝑓 ∗ 𝐼𝑓 + 𝐿𝑓 ∗ 𝐼𝑓′ + 𝐿𝑎𝑑 ∗ 𝐼𝐷′ (3)

𝑉𝑓 − 𝑟𝑓 ∗ 𝐼𝑓 − 𝐿𝑎𝑑 ∗ 𝐼𝐷 ′
𝐼𝑓 ′ =
𝐿𝑓
Para la implementación de las ecuaciones se considera solucionar cada ecuación
diferencial como semuestra en el diagrama de bloques de la Figura 5 y relacionando
con las otras varibales (p.e. corrientes del vanado) usando bloques de integración,
sumadores, ganancia y una fuente se representa con un escalón (Vf). Toda la
implementación que se muestra en este informe se realiza siguiendo el mismo
criterio.

Figura 5.Modelado en SIMULINK

 Para la determinación y solución de las ecuaciones diferencial se usa la herramienta
Live Editor de MATLAB, a continuación, se muestra un ejemplo de solución para las
ecuaciones (3) y (4), ver Figura 6.
0 = 𝐿𝑎𝑑 ∗ 𝐼𝑓 ′ + 𝑅𝐷 ∗ 𝐼𝐷 + 𝐿𝐷 ∗ 𝐼𝐷′ (4)

Figura 6.Solución de ecuación (3) y (4) en Live Editor

Para la transformación de señales se usa el bloque de la Figura 7, se considera que
el sistema d-q gira a una velocidad Wr, representando por una rampa.

Figura 7.Transformada a sistemas 𝜶𝜷 y d-q

4.1.2 Solución al problema N° 1 en SIMULINK del toolbox MATLAB

Considerando las condiciones iniciales indicado en la ecuación (1) - (2) y
considerando las ecuaciones que muestra la Tabla 2 se construye en SIMULINK
los diagramas de bloques para determinar la evoluación en el tiempo de las
corrientes del estator y de campo.

Figura 8.Diagrama de bloques Problema N° 1

a) Las corrientes en el devanado de excitación y en el devanado D
Después de aplicar una tensión Vf=50 Volts en t=0 se obtiene las corrientes
de campo que alcance un valor máximo de 89.23 A en el tiempo infinito y la
corriente en el devanado amortigiador D (en cortocircuito) en el tiempo
infinito regresa a su valor de cero (después de la perturbación) como se
esperaba, como se muestra en la Figura 9.

Figura 9.Corriente de Campo If(t) y del devanado D

b) Las tensiones inducidas en el estator d-q

Ante la perturbación de Vf = 50 Volts se obtiene los voltajes Vd y Vq
mostrado en la Figura 10.
 Figura 10.Tensión Vd y Vq

c) Determine y grafique la tensión inducida en los devanados α y β.

Considerando las ecuaciones de la máquina síncrona bifásica con el rotor
reflejado al estator de obtiene los volatjes en α y β.
𝑉𝛼(𝑡) = 𝐿𝑎𝑑 ∗ cos(𝑊𝑟 ∗ 𝑡) ∗ (𝐼𝑓 ′ + 𝐼𝐷′ ) − 𝐿𝑎𝑑 ∗ 𝑊𝑟 ∗ 𝑠𝑒𝑛(𝑊𝑟 ∗ 𝑡) ∗ (𝐼𝑓 + 𝐼𝐷)
𝑉𝛽(𝑡) = 𝐿𝑎𝑑 ∗ 𝑊𝑟 ∗ cos(𝑊𝑟 ∗ 𝑡) ∗ (𝐼𝑓 + 𝐼𝐷) + 𝐿𝑎𝑑 ∗ 𝑠𝑒𝑛(𝑊𝑟 ∗ 𝑡) ∗ (𝐼𝑓 ′ + 𝐼𝐷′ )
Resolviendo las ecuaciones diferenciales (3) y (4), se obtiene:

Se obtiene las tensiones considerando un tiempo infinito o muy grande:

𝑉𝛼(𝑡 → ∞) = 𝑉𝑚𝑎𝑥𝑆𝑒𝑛(𝜃) = −3.3654 ∗ 104 𝑠𝑒𝑛(𝜃)
𝑉𝛽(𝑡 → ∞) = 𝑉𝑚𝑎𝑥𝐶𝑜𝑠(𝜃) = 3.3654 ∗ 104 𝑐𝑜𝑠(𝜃)

Figura 11.Tensión 𝑽𝜶 y 𝑽𝜷
d) Determine y grafique las tensiones inducidas en los devanados a, b y c
de la máquina sincrona trifásica.
Para determinar tensiones en el sistema trifasico se usa la transformada de
Clarke, considerano Vo=0.
𝟐
𝑽𝒂 = √ ∗ 𝑽𝜶
𝟑

𝟐 𝑽𝜶 𝑽𝜷
𝑽𝒃 = √ ∗ (− + )
𝟑 𝟐 𝟐

𝟐 𝑽𝜶 √𝟑𝑽𝜷
𝑽𝒄 = √ ∗ (− − )
𝟑 𝟐 𝟐

Figura 12.Tensión 𝑽𝒂𝒃𝒄

e) Comentar los resultados si la bobina “Q” estuviera en circuito abierto.

Las ecuaciones considerando abierto el devanado Q se transforman como se
indica a continuación (color rojo) y las mismas se agregan al diagrama de
Bloques de la Figura 8.
𝑉𝑄 = 𝑟𝑄 ∗ 𝐼𝑄 + 𝐿𝑄 ∗ 𝐼𝑄′ (5)
𝑉𝑑 = 𝐿𝑎𝑑 ∗ (𝐼𝑓 ′ + 𝐼𝐷′ ) − 𝐿𝑎𝑞 ∗ 𝑊𝑟 ∗ 𝐼𝑄 (6)
𝑉𝑞 = 𝐿𝑎𝑑 ∗ 𝑊𝑟 ∗ (𝐼𝑓 + 𝐼𝐷) + 𝐿𝑎𝑞 ∗ 𝐼𝑄′ (7)
 Figura 13.Diagrama de Bloques con bobina “Q” estuviera en circuito abierto

Figura 14. IQ con bobina “Q” estuviera en circuito abierto

Se demuestra que al abrir la bobina Q no afecta los resultados, lo cual es

correcto ya que en este devanado no hay flujo.

4.2 CASO2: MÁQUINA SINCRONA TRIFASICA EN VACIO Y EL ESTATOR EN

CORTOCIRCUITO
Se considera la máquina síncrona de la Figura 15, utilizando el modelo d-q, para
resolver el caso 2.

Figura 15. Máquina sincrona trifasica en vacío y estator en cortocircuito

 Considerando la Figura 15 se puede tener las condiciones iniciales:
𝑽𝒂 = 𝑽𝒃 = 𝑽𝒄 = 𝑽𝜶 = 𝑽𝜷 = 𝑽𝒅 = 𝑽𝒒 = 𝑽𝑫 = 𝑽𝑸 = 𝟎 (8)

Teniendo las condiciones iniciales y las ecuaciones de la MS en los ejes d-q de la

Figura 4, se obtiene las ecuaciones indicadas en la Tabla 3.
Tabla 3. Ecuaciones de la MS caso 2

Se considera el mismo criterio incado en el numeral 4.1.1 para la elaboración del

diagrama debloques de las ecuaciones de la Tabla 3 en SIMULINK.

Figura 16. Diagrama de bloques Problema N° 2

a) Las corrientes en el devanado de excitación, asi como las corrientes las

bobinas D y Q.
Después de aplicar una tensión Vf=50 se obtiene la evoluación en el tiempo de la
corriente de campo, esta corriente en tiempo infinito tiene un valor de 90 A y las
corrientes del devanado amortiguador en el transitorio tratan de inhibir la
perturbación y en el tiempo infinito (estato estacionario) volver a un valor de cero.
 Figura 17. Corriente de Campo If y devanado ID-IQ

b) Las corrientes i d e i q en el estator d-q.

Aplicando la tensión de campo Vf que esta ubicado en el eje d, se nosta que tiene
mucha influencia en la corriente Id, esto por la aparición del flujo en dicho eje d.
 Figura 18. Corriente de Id y Iq

c) Determine y grafique las corrientes en los devanados α y β

Considerando la ecuación de transformación y los resultados de las corrientes Id –

Iq se determina que el valor máximo de corriente en un tiempo infinito tiene una
magnitud similar a la Id, es decir, la 𝑰𝒑𝒊𝒄𝒐 = 𝐼𝑑 = 𝟕𝟗. 𝟑 𝑽. La evoluación en el tiempo
se muestra en la Figura 19.

𝐼𝛼 = 𝐶𝑜𝑠(𝜃) ∗ 𝐼𝑑 − 𝑆𝑒𝑛(𝜃) ∗ 𝐼𝑞
𝐼𝛽 = 𝑆𝑒𝑛(𝜃) ∗ 𝐼𝑑 + 𝑐𝑜𝑠(𝜃) ∗ 𝐼𝑞
 Figura 19. Corrientes Iα - Iβ

d) Determine y grafique las corrientes en los devanados a, b y c de la máquina

sincrona trifásica.
Para determinar las corrientes en el sistema trifasico se usa la transformada de
Clarke, considerano Vo=0.
𝟐
𝑰𝒂 = √ ∗ 𝑽𝜶
𝟑

𝟐 𝑽𝜶 𝑽𝜷
𝑰𝒃 = √ ∗ (− + )
𝟑 𝟐 𝟐

𝟐 𝑽𝜶 √𝟑𝑽𝜷
𝑰𝒄 = √ ∗ (− − )
𝟑 𝟐 𝟐
 Figura 20. Corriente de abc

e) Comentar los resultados si la máquina no tuviera devanado amortiguador.

Los resultados se muestran en la Figura 21, se puede notar que considerando el
devanando amortiguador la corriente antes del segundo se muestra una curva
suave el incremento de la curva. Sin embargo, cuando no se considera el
devanando amortiguador, la evolución de las corrientes en el tiempo muestra una
onda o una pequeña oscilacoón y se amortigua después del seguno de la inyección
de la tensión en el devanado de campo Vf.
 Figura 21. Corriente de If(t) – Id con y sin devanado Amortiguador

5. CONCLUSIONES

• La transformada de Clarke convierte las componentes del dominio del tiempo de un

sistema de tres fases (de un marco abc) en dos componentes de un marco
estacionario ortogonal (αβ). La transformada de Park convierte las dos
componentes del marco αβ a un marco de referencia rotatorio ortogonal (dq) que
gira a Wr del rotor, lo cual es muy util para el estudio de la máquina sincrona.
• Con lo parametros entrgados, se determinó que la frecuencia es 25 Hz a la que gira
la máquina sincrona.
• Para el problema N°1, sin considerar el devanado amortiguador (eje D )se nota que
la corriente de campo (If) se atrza en la respuesta.

• Se verifica que el devanando de amortiguados en ambos casos, actua ante la

perturbacicón ( inyección de tensión en el devanado de campo Vf) y luego de varios
segundos en el estado estacionario regresa a un valor de cero en la corriente de
devanado.
• En el caso1 al abrir solo el devando de amortiguamiento Q, se observa que no tiene
una influencia significatitava en la evolución en el tiempo de las corrientes Iq, Id.
• Los devanado amortiguador tiene una influencia en los primeros instantes de la
perturbación y cuando no se considera estos, la corrientes en los ejes d-q muestran
pequeña ondulaciones que toman hasta 1 segundo para amortiguar la señal de
salida y cuando estos estan presentes la señal de corriente en los ejes d-q tiene
una respuest suave en todo momento de su evoluación en el tiempo.