Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica y Electrónica
Escuela Profesional de Ingeniería Eléctrica
Curso: EE354
Análisis de Sistemas de Potencia II
Expositor : Manfred F. Bedriñana Aronés
Doctor en Ingeniería Eléctrica
Ingeniero Electricista C.I.P. Nº 95644
2023
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
Capítulo 1
Elementos de Estabilidad de Sistemas
de Potencia
09/04/2023 [Link]@[Link] 2
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
Ecuación de Oscilación
09/04/2023 [Link]@[Link] 3
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
A. Introducción (1)
Operación bajo perturbaciones
◼ Bajo condiciones de operación normal:
◼ La posición relativa del eje del rotor y el eje de campo
magnético resultante es fija.
◼ El ángulo entre ellos es conocido como ángulo de
potencia o ángulo de torque.
◼ Durante perturbaciones:
◼ El rotor desacelera o acelera con respecto a la fmm
rotatoria síncrona del entrehierro, lo que origina
movimiento.
◼ La ecuación que describe este movimiento relativo se
denomina “Ecuación de Oscilación”.
09/04/2023 [Link]@[Link] 4
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
A. Introducción (2)
Estabilidad y nuevo punto de operación
◼ Si después de un periodo oscilatorio, el
rotor regresa a la velocidad síncrona
(frecuencia 60 Hz), el generador
mantendrá su estabilidad.
◼ Si la perturbación no envuelve un cambio Nuevo valor de
neto en la potencia de la máquina, el estado estacionario
rotor retorna a su posición original.
◼ Si existe un cambio en generación, carga
o en las condiciones de la red, el ángulo
de potencia alcanza un nuevo valor.
09/04/2023 [Link]@[Link] 5
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
A. Introducción (3)
Definición del ángulo de potencia
◼ Usando un generador de rotor cilíndrico (dos polos), el
ángulo de potencia 𝛿𝑟 viene a ser el ángulo entre la
fmm del rotor 𝐹𝑟 y la fmm resultante del entrehierro
𝐹𝑠𝑟 , ambos girando a la velocidad síncrona.
◼ El ángulo de potencia también es igual al ángulo entre la
tensión en vacío 𝐸 y la tensión del estator resultante 𝐸𝑠𝑟 .
◼ Si se incluyen la resistencia de armadura y el flujo de
dispersión:
◼ El ángulo de potencia 𝛿 se define como el ángulo entre
la tensión en vacío 𝐸 y la tensión terminal 𝑉.
09/04/2023 [Link]@[Link] 6
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
B. Desarrollo de la ecuación de oscilación (1)
Torque electromagnético y mecánico, ley de rotación
◼ Considere un generador síncrono desarrollando un torque
electromagnético 𝑇𝑒 girando a una velocidad síncrona.
◼ Si 𝑇𝑚 es el torque mecánico desarrollado, entonces para la
operación en estado estacionario (despreciando las pérdidas):
𝑇𝑒 = 𝑇𝑚 (1)
Acción del torque en
◼ Desde una condición en estado estacionario, una un generador síncrono
perturbación resulta en un torque acelerante (𝑇𝑚 > 𝑇𝑒 ) o
desacelerante (𝑇𝑚 < 𝑇𝑒 ) en el rotor:
𝑇𝑎 = 𝑇𝑚 − 𝑇𝑒 (2) Turbinas del
motor primo
◼ Si J es el momento de inercia combinado del motor primo
y generador (un único cuerpo rígido), despreciando los
torques de amortiguamiento y de fricción, de la ley de
rotación se tiene:
𝑑2 𝜃𝑚 Rotores acoplados en el eje
𝐽 = 𝑇𝑚 − 𝑇𝑒 (3) (unidad turbo vapor)
𝑑𝑡 2
HP: Turbina de alta presión
donde: 𝜃𝑚 es el desplazamiento angular del rotor con LP-A, LP-B: Turbinas de baja presión
respecto a un eje estacionario en el estator. Gen: Generador
Exc: Excitattiz
09/04/2023 [Link]@[Link] 7
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
B. Desarrollo de la ecuación de oscilación (2)
Ángulo del rotor y ley de rotación
◼ Desde que estamos interesados en la velocidad del rotor relativa a la velocidad síncrona,
la referencia angular es elegida relativa a la referencia rotatoria síncrona moviéndose a
velocidad angular constante 𝜔𝑠𝑚 , esto es:
𝜃𝑚 = 𝜔𝑠𝑚 𝑡 + 𝛿𝑚 (4)
donde: 𝛿𝑚 es el ángulo del rotor antes de la perturbación en el tiempo 𝑡 = 0, medido
desde el marco de referencia rotatoria síncrona.
◼ Derivando la ecuación (4) se obtiene la velocidad angular del rotor:
𝑑𝜃𝑚 𝑑𝛿𝑚
𝜔𝑚 = 𝑑𝑡
= 𝜔𝑠𝑚 + 𝑑𝑡
(5)
y la aceleración del rotor será:
𝑑2 𝜃𝑚 𝑑2 𝛿𝑚
𝑑𝑡 2
= 𝑑𝑡 2
(6)
◼ Sustituyendo (6) en (3) se tiene:
𝑑2 𝛿𝑚
𝐽 𝑑𝑡 2
= 𝑇𝑚 − 𝑇𝑒 (7)
09/04/2023 [Link]@[Link] 8
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
B. Desarrollo de la ecuación de oscilación (3)
Constante de inercia
◼ Multiplicando la ecuación (7) por 𝜔𝑚 :
𝑑2 𝛿𝑚
𝐽𝜔𝑚 𝑑𝑡 2 = 𝜔𝑚 𝑇𝑚 − 𝜔𝑚 𝑇𝑒 (8)
◼ Desde que la velocidad angular multiplicado por el torque
es igual a la potencia, la ecuación anterior se coloca
en término de potencia:
𝑑2 𝛿𝑚
𝐽𝜔𝑚 2 = 𝑃𝑚 − 𝑃𝑒 (9)
𝑑𝑡
◼ La cantidad 𝐽𝜔𝑚 es llamada de constante de inercia y es denotada por M. Incluyendo la
energía cinética de las masas rotantes 𝑊𝑘 se tiene:
1 1
𝑊𝑘 = 𝐽𝜔𝑚 2 = 𝑀𝜔𝑚
2 2
(10)
2𝑊
o 𝑀= 𝜔 𝑘 (11)
𝑚
Comentario
◼ Mientras M es llamado de constante de inercia, este no es realmente una constante cuando la
velocidad del rotor varía desde la velocidad síncrona.
◼ Sin embargo, cuando 𝜔𝑚 no cambia en cantidades grandes antes de la pérdida de
sincronismo, M es evaluado en la velocidad síncrona y permanece constante:
2𝑊
𝑀=𝜔 𝑘 (12)
𝑠𝑚
09/04/2023 [Link]@[Link] 9
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
B. Desarrollo de la ecuación de oscilación (4)
Constante de inercia, ecuación de oscilación
◼ La ecuación de oscilación en términos de la constante de inercia será:
𝑑2 𝛿𝑚
𝑀 𝑑𝑡 2 = 𝑃𝑚 − 𝑃𝑒 (13)
◼ Es más conveniente escribir la ecuación de oscilación en términos del ángulo de potencia
eléctrica 𝛿. Si p es el número de polos del generador síncrono, el ángulo de potencia
eléctrica se define como:
𝑝
𝛿 = 2 𝛿𝑚 (14)
𝑝
también: 𝜔 = 𝜔𝑚 (15)
2
◼ La ecuación de oscilación en términos del ángulo de potencia eléctrica será:
2 𝑑2 𝛿
𝑝
𝑀 𝑑𝑡 2
= 𝑃𝑚 − 𝑃𝑒 (16)
◼ Desde que el análisis de sistemas de potencia se realiza en valores por unidad, se divide la
ecuación (16) por la potencia base 𝑆𝐵 y sustituyendo M de (12):
2 2𝑊𝑘 𝑑2 𝛿 𝑃𝑚 𝑃
𝑝 𝜔𝑠𝑚 𝑆𝐵 𝑑𝑡 2
= 𝑆𝐵
− 𝑆𝑒 (17)
𝐵
09/04/2023 [Link]@[Link] 10
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
B. Desarrollo de la ecuación de oscilación (5)
Constante de inercia por unidad y ecuación de oscilación
◼ Se define la constante de inercia por unidad H como:
𝑒𝑛𝑒𝑟𝑔í𝑎 𝑐𝑖𝑛é𝑡𝑖𝑐𝑎 𝑒𝑛 𝑀𝐽 𝑎 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑛𝑜𝑚. 𝑊𝑘
𝐻 = 𝑝𝑜𝑡𝑒𝑛𝑐𝑖𝑎 𝑛𝑜𝑚𝑖𝑛𝑎𝑙 𝑑𝑒 𝑙𝑎 𝑚á𝑞𝑢𝑖𝑛𝑎 𝑒𝑛 𝑀𝑉𝐴 = 𝑆𝐵
(18)
◼ H tiene unidades en segundos. El valor de H está entre 1 – 10 segundos, dependiendo
del tamaño y tipo de máquina.
◼ Valores típicos para la constante de inercia H (en segundos), respecto a la potencia MVA de
la máquina.
Comentarios
◼ Significado físico de la constante de inercia H:
◼ Cuanto mayor es el valor de H, mayor resistencia presenta la máquina a cambiar el ángulo de
potencia 𝛿.
09/04/2023 [Link]@[Link] 11
Escuela Profesional de Ingeniería Eléctrica, FIEE-UNI
Curso EE354 – Análisis de Sistemas de Potencia II
B. Desarrollo de la ecuación de oscilación (6)
Ecuación de oscilación
◼ Sustituyendo en (17) en términos de H:
2 2𝐻 𝑑2 𝛿
𝑝 𝜔𝑠𝑚 𝑑𝑡 2
= 𝑃𝑚(𝑝.𝑢.) − 𝑃𝑒(𝑝.𝑢.) (19)
donde 𝑃𝑚(𝑝.𝑢.) y 𝑃𝑒(𝑝.𝑢.) están en valores por unidad.
◼ La velocidad angular eléctrica se relaciona con la velocidad angular mecánica por
2
𝜔𝑠𝑚 = 𝑝 𝜔𝑠 . Sustituyendo en (19):
2𝐻 𝑑2 𝛿
= 𝑃𝑚(𝑝.𝑢.) − 𝑃𝑒(𝑝.𝑢.) (20)
𝜔𝑠 𝑑𝑡 2
◼ La ecuación anterior se expresa también en términos de frecuencia 𝑓0 , 𝑃𝑚 y 𝑃𝑒 están
en por unidad:
𝐻 𝑑2 𝛿
= 𝑃𝑚 − 𝑃𝑒 (21)
𝜋𝑓0 𝑑𝑡 2
09/04/2023 [Link]@[Link] 12