“Ecuación de oscilación para en

análisis de estabilidad”
Esteban Zambrano esteban.zambrano5216@utc.edu.ec

(Universidad Técnica de Cotopaxi UTC) Ingeniería Eléctrica

Latacunga - Ecuador

Resumen. En el presente artículo se mostrará la genera inestabilidad en el sistema. Con base

ecuación de oscilación para el análisis de en el estudio de severidad, se pueden
estabilidad, así también como definiciones de determinar algunas medidas correctivas, que
cada una y su aplicación dentro de los sistemas incluyen: tiempo crítico de falla, diseño de
eléctricos de potencia. esquemas de reparto de carga y diseño de
planes de apagado escalonado.
Abstract. This article will show the oscillation
equation for stability analysis, as well as II. DESARROLLO
definitions of each one and their application in
electrical power systems. A. ESTABILIDAD TRANSITORIA

Palabras claves. Síncrona, Parámetros, Se refiere a la capacidad del sistema

Operación, Estabilidad, Falla, Perturbación eléctrico para mantener el sincronismo
cuando es sometido a una perturbación
I. INTRODUCCIÓN fuerte, por ejemplo, una falta en la red de
transporte, perdida de generación o perdida
La estabilidad en los sistemas eléctricos de de una cantidad importante de carga.
potencia es una propiedad que les permite a
los sistemas mantenerse en un estado de
operación en equilibrio bajo condiciones
operativas normales y recuperar otro estado
de equilibrio después de haber sido
sometidos a una perturbación.
Al referirse a máquinas sincrónicas
interconectadas entre si se debe tener
presente que pueden perder su estabilidad
debido al: tipo de falla, localización y
tiempo de despeje. Figura 1. Estabilidad transitoria

La estabilidad angular está relacionada con B. ESTABILIDAD ESTACIONARIA

el comportamiento de los ángulos del rotor Engloba a la capacidad de un sistema de
bajo perturbaciones grandes o pequeñas.
energía eléctrica para mantener su
Durante una falla, el ángulo del rotor del
condición inicial después de una pequeña
generador comienza a oscilar, si la falla no
interrupción o para alcanzar una condición
es tan severa, el sistema se estabiliza hasta
que el ángulo del rotor vuelve al estado muy cercana a la inicial cuando la
inicial o a un nuevo estado estable, mientras perturbación aún está presente.
que la falla es demasiado grande, el ángulo
del rotor oscila sin amortiguamiento, lo que
 Figura 2. Estado estacionario
Figura 4. Sistema inestable de un grupo de
C. ESTABILIDAD ANGULAR generadores.

Esta es la capacidad requerida para que una D. RELACIÓN ENTRE ESTABILIDAD

máquina síncrona reanude su ANGULAR Y ESTABILIDAD TRANSITORIA
funcionamiento normal después de una falla
menor o mayor para que pueda mantener o La estabilidad transitoria se usa cuando el
restablecer el equilibrio entre la energía SEP experimenta grandes perturbaciones.
eléctrica y la energía mecánica. Si no se donde el tiempo de investigación de
alcanza este equilibrio, el sistema se vuelve estabilidad es de 3-5 segundos después de
inestable, aumenta las oscilaciones que ocurre la falla, puede extenderse hasta
angulares y pierde sincronización [6][14]. 10 segundos dependiendo de la severidad y
probabilidad de ocurrencia de la falla, y es
La figura 3 muestra el sistema en estado principalmente para el desempeño en
estable después de una perturbación. La ángulos del generador del sistema
acción del generador es oscilante donde los [3][10][16].
ángulos no se proyectan demasiado.
E. CURVA DE POTENCIA MÁXIMA

Para este caso la potencia máxima viene

dada por la siguiente relación, en donde el
valor de la tensión nominal y la máxima
corriente de armadura participan de la
manera que se muestra:
S = V x Ia max (1)

De una manera gráfica se vería de la

siguiente manera, donde la capacidad del
generador, donde se puede determinar
Figura 3.Sistema estable de un grupo de trazando una recta paralela al eje Q distante
generadores.
en una magnitud igual a la potencia activa
En la siguiente figura se podrá observar el máxima que impone el rotor.
mismo sistema de generación anterior, pero
los ángulos proyectan una separación
considerable, indicando que el sistema esta
perdiendo sincronismo por lo que se
encuentra en inestabilidad.

Figura 5. Curva de potencia máxima.

La corriente de excitación máxima trae

consigo un valor máximo de la fuerza
electromotriz inducida. El lugar geométrico
de este límite se representa como una velocidad en función de la frecuencia
semicircunferencia en el semiplano positivo Ѡ=2πf, siendo f la frecuencia en ciclos (Hz).
P de centro en el punto (- 𝒗^𝟐/𝑿𝒔 , 0) el cual
1 Ѡ
produce un ángulo en este caso representado 𝐸𝑐 = GH = IѠ2 = M (2)
2 2
por σ de acuerdo a la ecuación anterior.
Se supone ahora:

Figura 6. Curva de potencia máxima. Figura 8.Máquina síncrona conectada a una barra
de potencia infinita.
F. MÉTODO DE LAS ÁREAS IGUALES

En general, la estabilidad o inestabilidad de Ahora, durante su funcionamiento esta

un sistema depende en gran medida de si el puede entregar una potencia que viene
error persiste (permanentemente) o se dada por:
elimina con el tiempo. Hay un tiempo que se
mide relativo al inicio de la falla. Esto se
llama el "tiempo crítico de liberación". La
única determinación del tiempo de despeje
crítico utilizando la curva de vibración es un
método de aproximación iterativa, que no es
eficiente. El método consiste en calcular la
curva de vibración para falla sostenida y
luego para varios tiempos de liberación
estimados. Los tiempos críticos de
liberación se determinan a partir del análisis
de todas las diferentes curvas.
Figura 9. Máquina síncrona conectada a una barra
de potencia infinita.

Si no se considera:

- El par originado por el

rozamiento mecánico.

- El rozamiento del aire.

- Pérdidas en el núcleo.
Figura 7.Método de las áreas iguales.
- Pérdidas por corrientes de
G. ECUACIÓN DE OSCILACIÓN Focault en los arrollados
amortiguadores.
Si se considera que la parte giratoria,
corresponde al rotor de una maquina que Entonces cualquier diferencia entre la
gira a una velocidad Ѡ en grados eléctricos potencia mecánica (𝑃𝑚𝑒𝑐 ) y la eléctrica
por segundo, entonces, se puede tomar la (𝑃𝑒𝑙𝑒𝑐 ) debe actuar sobre la máquina
como:
𝑇𝑎𝑐𝑒𝑙 = 𝑇𝑚𝑒𝑐 - 𝑇𝑒𝑙𝑒𝑐 (3)

- Esta potencia acelerante es causada

por una diferencia entre el torque
mecánico y el torque
electromagnético
Figura 11. Torque girando a una velocidad síncrona.
La potencia acelerante queda expresada
por: 𝜃𝑚𝑒𝑐 = 𝛿 + Ѡ𝑠 t (7)
𝑑𝜃𝑚𝑒𝑐 𝑑𝛿
𝑃𝑎𝑐𝑒𝑙 = 𝑇𝑎𝑐𝑒𝑙 Ѡ (4) = + 𝜔𝑠 (8)
𝑑𝑡 𝑑𝑡
Relación entre el torque de aceleración y el 𝑑 2 𝜃𝑚𝑒𝑐 𝑑2 𝛿
= (9)
momento de inercia: 𝑑𝑡 2 𝑑𝑡 2

𝑇𝑎𝑐𝑒𝑙 = Iα (5) Se deduce que la aceleración absoluta es

igual a la relativa por lo que tenemos:
Aplicando la segunda Ley de Newton
𝑑2
aplicada a los torques quedaría: I = 𝛴 − 𝑇 = 𝑇𝑚𝑒𝑐 - 𝑇𝑒𝑙𝑒𝑐 = 𝑇𝑎𝑐𝑒𝑙 (10)
𝜕𝑡 2
𝑑 2 𝜃𝑚𝑒𝑐
ƩT = I = 𝑇𝑚𝑒𝑐 = 𝑇𝑒𝑙𝑒𝑐 - 𝑇𝑎𝑐𝑒𝑙 (6) Multiplicando en ambos lados por la
𝑑𝑡 2
velocidad angular tenemos:
𝜕2𝛿
𝜔𝐼 = 𝜔𝑇𝑚𝑒𝑐 - 𝜔𝑇𝑒𝑙𝑒𝑐 = 𝜔𝑇𝑎𝑐𝑒𝑙 (11)
𝑑𝑡 2

Resultado de la definición del momento

angular:
𝑑2 𝛿
𝑀′ = 𝑃𝑚𝑒𝑐 - 𝑃𝑒𝑙𝑒𝑐 = 𝑃𝑎𝑐𝑒𝑙 (12)
𝑑𝑡 2

Donde 𝑀′ :
Figura 10. Segunda Ley de Newton aplicada a los
torques. 𝜔𝐼𝜔𝑠
𝑀′ = 𝐼𝜔 = (13)
𝜔𝑠

𝜔𝑀
𝑀′ = (14)
𝜔𝑠
El torque acelerante, será positivo cuando
el torque mecánico supere al La cantidad es conocido como el momento
electromagnético, con lo que la máquina se angular, reescribiendo la ecuación resulta:
acelera, de lo contrario pierde aceleración.
𝜔 𝑑2 𝛿
𝑀 2 = 𝑃𝑚𝑒𝑐 - 𝑃𝑒𝑙𝑒𝑐 = 𝑃𝑎𝑐𝑒𝑙 (15)
𝑇𝑎𝑐𝑒𝑙 > 0 𝜔𝑠 𝑑𝑡

Todos los siguientes datos giran a una Se puede aproximar a 1:

velocidad síncrona 𝜃𝑚𝑒𝑐 , el ángulo de rotor 𝜔
=1 (16)
𝜔𝑠
medido respecto a una referencia Ѡ𝑠 ,
velocidad síncrona de la máquina y el El momento angular de M de una maquina
desplazamiento angular del rotor respecto a no es constante, pero como varia la
un eje 𝛿. velocidad angular se puede tomarle en
cuenta como una variable.
𝜔 𝜔
𝑀′ = M; ≈ 1; 𝑀′ ≈ 𝑀 (17)
𝜔𝑠 𝜔𝑠
𝑑2 𝛿 alimentar cargas aisladas o para
= 𝑃𝑚𝑒𝑐 - 𝑃𝑒𝑙𝑒𝑐 = 𝑃𝑎𝑐𝑒𝑙 (18)
𝑑𝑡 2
entregar potencia a una red eléctrica
La ecuación anterior recibe el nombre de
compleja.
ecuación de oscilación y se caracteriza la
reposición del rotor de la máquina síncrona IV. BIBLIOGRAFIA
durante la perturbación. [1] N. Zhou, D. Trudnowski, J. W.
La ecuación de oscilación es una ecuación Pierre, S. Sarawgi, and N. Bhatt, “An
diferencial transcendental de segundo
algorithm for removing trends from
orden y su solución da origen a una integral
elíptica power-system oscillation data,” in

La ecuación de oscilación se puede tratar 2008 IEEE Power and Energy Society
en cantidades por unidad: General Meeting - Conversion and

M= 𝐽 =
2𝐺𝐻
(19) Delivery of Electrical Energy in the
𝜔𝑠
21st Century, 2008, pp. 1–7.
2𝐺𝐻 𝑑 2 𝛿
= 𝑃𝑚𝑒𝑐 - 𝑃𝑒𝑙𝑒𝑐 = 𝑃𝑎𝑐𝑒𝑙 (20) [2] J. Fang, C. Su, Z. Chen, H. Sun, and P.
𝜔𝑠 𝑑𝑡 2
Lund, “Power System Structural
Se divide ambos lados de la expresión por Vulnerability Assessment Based on
la potencia base 𝑆𝑏𝑎𝑠𝑒 : an Improved Maximum Flow
2𝐻𝑏 𝑑 2 𝛿 𝑃𝑚𝑒𝑐− 𝑃𝑒𝑙𝑒𝑐 𝑃𝑒𝑙𝑒𝑐 Approach,” IEEE Trans. Smart Grid,
= = (21) vol. 9, no. 2, pp. 777–785, Mar.
𝜔𝑠 𝑑𝑡 2 𝑆𝑏𝑎𝑠𝑒 𝑆𝑏𝑎𝑠𝑒
2018.
2𝐻𝑏 𝑑 2 𝛿
= 𝑃𝑎𝑐𝑒𝑙 [ P.U.] (22) [3] P. Kundur ; J. Paserba ; V. Ajjarapu ;
𝜔𝑠 𝑑𝑡 2
G. Andersson ; A. Bose ; C.
Cañizares ; N. Hatziargyriou ; D.
Hill ; A. Stankovic ; C. Taylor ; T. Van
III. CONCLUSIONES
Cutsem ; V. Vittal, “Definition and
• Se pueden usar varios métodos Classification of Power System
numéricos para resolver la ecuación Stability IEEE/CIGRE Joint Task
de oscilación. Force on Stability Terms and
Definitions,” IEEE Trans. Power
• El método Runge-Kutta permite un
Syst., vol. 19, no. 3, pp. 1387–1401,
mejor análisis de las curvas de
Aug. 2004.
oscilación que surgen en el
[4] M.-E. Hamedani-golshan, A.
problema de estabilidad con mayor
Davoudi, and A. Bidram, “Loading
precisión.
constraints for first swing stability
• Mediante este trabajo investigativo, margin enhancement of distributed
se logró entender de una manera generation,” IET Gener. Transm.
clara y concisa como es el Distrib., vol. 6, no. 12, pp. 1292–
funcionamiento de la maquina 1300, Dec. 2012.
síncrona, así también como su [5] Xiaoyu Wang and W. Freitas,
ecuación de oscilación para el “Influence of Voltage Positive
análisis de estabilidad. Feedback Anti-Islanding Scheme on
• Para que un SEP sea confiable debe Inverter-Based Distributed
ser seguro y para ellos debe ser Generator Stability,” IEEE Trans.
estable. Power Deliv., vol. 24, no. 2, pp. 972–
973, Apr. 2009
• La máquina sincrónica se puede
[6] R.H. PARK, "TWO-REACTION
utilizar como generador, tanto para THEORY OF SYNCHRONOUS
 MACHINES", AIEE TRANS,PARTE PROC. IEE (LONDON), VOL. 124,
I, VOL. 48, PP. 716-730, 1929; NO. 12, PP. 1170-1178, 1977.
PARTE II, VOL. 52, PP. 352-355, [16] F.P. DEMELLO Y J.R.
1933. RIBEIRO , "DERIVATION OF
[7] G. SHACKSHAFT Y A.T. PORAY SYNCHRONOUS MACHINE
"ESTABILIDAD EN LOS SEP". PARAMETERS FROM TESTS"
SHARE AND DISCOVER (DERIVACIÓN DE LOS
KNOWLEDGE ON SLIDE SHARE , PARÁMETROS DE LAS MÁQUINAS
CENTENO ASTUDILLO 2020. SÍNCRONAS A PARTIR DE LAS
[8] C. CONCORDIA , SYNCHRONOUS PRUEBAS ), IEEE T RANS., VOL.
MACHINES - THEORY AND PAS096, PP. 1211-1218,
PERFORMANCE , JOHN WILEY & JULIO /AGOSTO DE 1977
SONS, 1952.
[9] B. ADKINS, THE GENERAL V. AUTOR
THEORY OF ELECTRIC MACHINES,
CHAPMAN AND HALL, 1964. ESTEBAN FERNANDO ZAMBRANO
[10] G. SHACKSHAFT , "NEW CAJIAS.
APPROACH TO THE
DETERMINATION OF
SYNCHRONOUS MACHINE
PARAMETERS FROM TESTS",
PROC. IEE (LONDON), VOL. 121,
NO. 11, PP. 1385- 1392, 1974.
[11] M. CANAY, "CAUSAS DE
LAS DISCREPANCIAS EN EL
CÁLCULO DE LAS CANTIDADES
DEL ROTOR Y EL DIAGRAMA
Nacido en Latacunga el 26 de diciembre
EQUIVALENTE EXACTO DE LAS
MÁQUINAS SÍNCRONAS ", IEEE
de 1997, estudio el bachillerato en el
TRANS., VOL. PAS- 88, PP. 1114- Colegio Americano Gypsi, México -
1120, JULIO DE 1969. Campeche, graduándose en el año 2015
[12] INFORME EPRI EL-1424, obteniendo el titulo de Bachiller en
"DETERMINATION ciencias, en la actualidad se encuentra
OFS YNCHRONOUS MACHINE cursando el último año en la Facultad
STABILITY CONSTANTS", VOL. 2, (CIYA) en la Carrera de Ingeniería
PREPARADO POR ONTARIO Eléctrica en la Universidad Técnica de
HYDRO, DICIEMBRE DE 1980. Cotopaxi, su línea de investigación está
[13] NORMA ANSI/IEEE 115- ligado a la automatización.
1983, GUÍA IEEE:
PROCEDIMIENTOS DE PRUEBA
PARA MÁQUINAS SÍNCRONAS .
[14] Y. TAKEDA Y B. ADKINS,
"DETERMINATION OF
SYNCHRONOUS MACHINE
PARAMETERS ALLOWING FOR
UNEQUAL MUTUAL
INDUCTANCES", PROC. IEE
(LONDON), VOL. 121, NO. 12, PP.
1501-1504, DECEMBER 1974.
[15] G. SHACKSHAFT Y A.T.
PORAY, "IMPLEMENTATION OF
NEW APPROACH TO
DETERMINATION OF
SYNCHRONOUS MACHINE
PARAMETERS FROM TESTS",