Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
TAREA N° 4: Sistemas Multimáquina
SIMULACIÓN DE SISTEMAS ELÉCTRICOS DE POTENCIA
Estudiante: Benjamín Godoy Ayala
Profesor Fernando Torres Leal
RESUMEN
El objetivo de este trabajo es estudiar las principales variables y fenómenos eléctricos asociados a la
implementación de un grupo de generadores que incluyen un sistema de control que permite mejorar la
respuesta transitoria ante distintos tipos de contingencia. Estos estudios incluyen análisis estático
mediante flujo de carga, análisis modal y estudios dinámicos.
Como resultados de estos estudios, se obtienen gráficas y curvas que ilustran la importancia del
estudio de los ángulos rotóricos de los generadores, su influencia sobre las variables eléctricas de interés y
como estos determinan si el sistema converge ante una determinada falla. Además, es posible desprender
que los estudios de flujo de carga pueden, en algunos casos, ser predictivos en cuanto a los valores de
convergencia, a pesar de no considerar sistemas de control. Por último, es posible comprobar que un
correcto balance entre demanda y generación, y sistemas de control en generadores, contribuyen a una
estabilidad en el sistema, lo que puede comprobarse mediante análisis modal y convergencia en curvas
transitorias.
INTRODUCCIÓN
Este informe responde a las actividades de la tarea N°4 del módulo “Simulación de Sistemas Eléctricos
de Potencia” impartido por el profesor Fernando Torres Leal, Universidad de Talca. La tarea es abordada
mediante el software DIgSILENT PowerFactory. Cada una de las actividades se aborda de la siguiente
forma:
i. Se realizan estudios estáticos de flujo de carga pre y post contingencia mediante el uso de los
datos flexibles. Estos son insertados en un archivo Excel, a partir del cual se produce una
gráfica comparativa pre y post contingencia. Las variables analizadas son el voltaje en barras y
la potencia activa producida por los generadores.
ii. Se lleva a cabo un análisis modal pre y post contingencia
iii. Se realizan estudios dinámicos mediante simulación RMS de las variables eléctricas de interés.
En este caso, corresponden a Voltaje y frecuencia en barras, y ángulos rotóricos de
generadores. Todas las gráficas incluyen comentarios.
1
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Teniendo en cuenta la forma en que se realiza cada actividad, a continuación, se presenta la temática
de cada sección:
- En la sección I. Se aborda la respuesta a la actividad a): Aumento de demanda, extensión del
sistema en estudio y conexión de más de 4 generadores.
- En la sección II. Se aborda la respuesta a la actividad b): Impacto de carga de 30 MW con
generadores en servicio.
- En la sección III. Se aborda la respuesta a la actividad c): Rechazo de carga de 30 MW con
generadores en servicio.
- En la sección IV. Se aborda la respuesta a la actividad d): Con generadores en servicio, pérdida
de unidad con mayor inyección de potencia al sistema.
- En la sección V. Se aborda la respuesta a la actividad e): Con generadores en servicio, pérdida de
una línea de transmisión relevante.
DESARROLLO
El desarrollo de este informe está contextualizado en la implementación de nuevas centrales que
contribuyan a la generación y al control de frecuencia además de la central Ralco implementada en el
trabajo 3. En dicho estudio se presentan respuestas indeseables que es importante tener en cuenta para
contrastar con las respuestas obtenidas en este trabajo. Estas pueden resumirse en ciclos límite y
divergencia. Por último, en este estudio se busca establecer una buena respuesta transitoria de señales sin
la contribución de elementos controladores de voltaje y frecuencia (sin considerar los sistemas incluidos en
las máquinas generadoras implementadas), por lo cual se establecen todas las medidas necesarias para
lograr este objetivo.
I. ACTIVIDAD A: Aumento de demanda y generación.
En primera instancia, para esta actividad se incorpora un nuevo subsistema a partir de la barra 5. El
criterio para la elección de esta barra como origen de esta nueva red es su ubicación, cercana a la central
Ralco y a las generadoras hidroeléctricas; su voltaje de 220 kV el cual puede otorgar robustez al nuevo
subsistema; y su característica, pues es una barra sin muchos cubículos y conexiones.
La nueva red alimenta 5 nuevos consumos a niveles de tensión variados (ver Fig.1) y consta de líneas
y transformadores obtenidos directamente de la base de datos del Sistema Eléctrico Nacional. Los niveles
de tensión del nuevo sistema son 220 kV, 110 kV y 66 kV y los colores que se visualizan obedecen al
estándar chileno.
2
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Las cargas suman en su totalidad 190 MVA y los factores de potencia de cada una son asignadas según
la norma técnica de seguridad y calidad de servicio. Las características de estas cargas se muestran en la
tabla 1.
Fig. 1: Nuevo sistema a partir de la barra 5.
Tabla 1: Nuevas demandas asociadas al nuevo subsistema.
Carga Consumo en MVA Factor de Potencia
Carga A.1 – 220 kV 70 0,98 i
Carga A.2 – 110 kV 45 0,98 i
Carga A.3 – 66 kV 25 0,96 i
3
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Carga A.4 – 66 kV 25 0,96 i
Carga A.5 – 66 kV 25 0,96 i
Ante estos nuevos consumos, se implementan 8 nuevas generadoras: 6 hidroeléctricas de 35 MW y
13,8 kV y 2 termoeléctricas, una de 30 MW y otra de 37 MW, ambas a 13,8 kV (ver tabla 2). El criterio para
elegir la potencia de estas centrales generadoras toma en consideración el hecho de que se busca mejorar
la respuesta de control de frecuencia ante contingencias, por lo tanto, todos los equipos mencionados
anteriormente poseen un control respondiendo a esta necesidad. Además, la central Ralco funciona muy
cerca del límite de su capacidad, de manera que es conveniente sobredimensionar la generación,
considerando posibles contingencias de impactos de carga y las mayores pérdidas de transmisión
provocadas por la extensión del sistema.
Tabla 2: Potencias nominales de generadores.
Generador Potencia Activa en MW Potencia Reactiva en MVar
GS Ralco 381,9 125,5
GS Hidroeléctrica 1 – Barra 25 35 10
GS Hidroeléctrica 2 – Barra 26 35 5
GS Hidroeléctrica 3 – Barra 29 35 7
GS Hidroeléctrica 4 – Barra 32 35 7
GS Hidroeléctrica 5 – Barra 30 35 7
GS Hidroeléctrica 6 – Barra 31 35 7
GS Termoeléctrica 1 – Barra 27 30 5
GS Termoeléctrica 2 – Barra 28 37 5
II. ACTIVIDAD B: Impacto de Carga.
i. Flujos de carga pre y post contingencia
Las figuras 2 y 3 ilustran los niveles de tensión y potencia activa en barras del sistema y en
generadores respectivamente. Se observa que, ante un aumento de carga, los niveles de tensión no se ven
contundentemente alterados. Sin embargo, es posible apreciar una ligera caída en los valores de voltaje en
las barras. Esto es esperable, pues una modificación en la demanda de potencia activa altera el factor de
potencia de los elementos del sistema, provocando una variación en la potencia reactiva y por lo tanto en el
voltaje. Esta misma tendencia se presenta ante un sistema con solo una central generadora. Por otro lado,
se aprecia que la central que aporta la potencia activa demandada en la contingencia es suministrada por
la central Ralco, mientras que las otras generadoras permanecen a niveles constantes de operación. Esto
es esperable pues la central Ralco está conectada a la barra de oscilación necesaria para resolver el
problema de flujo de carga, actuando para suplir potencia ante nuevos consumos, desbalances y pérdidas.
4
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Voltajes pre y post contingencia
1.02
1
0.98
0.96
0.94
Voltaje en p.u.
0.92
0.9
0.88
0.86
0.84
0.82
01 03 05 07 09 11 13 15 17 19 21 23 25 27 29 31
ra ra ra ra ra ra ra ra ra ra ra ra ra ra ra ra
ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar
B
Voltaje pre contingencia Voltaje post contingencia
Fig.2: Voltajes pre y post contingencia de impacto de carga.
Potencia activa pre y post contingencia
350
300
250
Potencia en MW
200
150
100
50
0
GS Ralco Hidro 1 Hidro 2 Hidro 3 Hidro 4 Hidro 5 Hidro 6 Termo 2 Termo1
MW pre contingencia MW post contingencia
Fig.3: Potencias activas pre y post contingencia suministrados por generadores.
5
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
ii. Análisis modal pre y post contingencia
El análisis modal para los casos de pre y post contingencia es mostrado en las figuras 4 y 5. Se
aprecia que, aunque los valores propios del sistema permanecen en el semiplano izquierdo, una
significativa cantidad de estos se encuentran justo al límite del eje imaginario. Aunque se cumple el criterio
que determina la estabilidad ante pequeñas perturbaciones, es deseable que un sistema posea sus valores
propios lo más a la izquierda posible en el plano mostrado en las figuras. Además, es posible observar que
existen valores propios con un valor no nulo de componente imaginaria, lo que implica oscilaciones en el
sistema.
Fig.4: Análisis modal pre contingencia.
6
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.5: Análisis modal post contingencia.
iii. Análisis dinámico
En la figura 6 se muestra la respuesta transitoria de voltaje y frecuencia en barras ante la contingencia.
En primera instancia, se aprecia que el impacto de carga provoca una drástica disminución en la frecuencia
del sistema. Sin embargo, esta contingencia logra ser mitigada en un porcentaje considerable gracias al
control presente en las generadoras. Cabe destacar que esto corresponde al control primario de frecuencia
y que, en la realidad, se establece un control jerárquico que reestablece la frecuencia a su nivel nominal.
Por otro lado, en la figura 7 se ilustra la respuesta de ángulos rotóricos ante impacto de carga. Se
observa que, ante esta contingencia, los ángulos rotóricos disminuyen su valor en una primera instancia.
Sin embargo, el control presente en las generadoras logra que los ángulos converjan a un determinado
valor. Cabe destacar que una de las principales aristas de la estabilidad transitoria está dada por la
convergencia de estos ángulos. Si los ángulos convergen, entonces existe una muy alta probabilidad de
que las magnitudes del sistema también lo hagan. Tomando en cuenta lo anterior, todas las dinámicas
obtenidas son consistentes con la teoría asociada a la estabilidad transitoria.
7
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.6: Respuesta dinámica de Voltajes y frecuencia.
Fig.7: Respuesta dinámica de ángulos rotóricos. La curva de ordenada cero es la referencia (Ralco).
8
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
III. ACTIVIDAD C: Rechazo de Carga.
i. Flujos de carga pre y post contingencia
En las figuras 8 y 9 se muestran, respectivamente, las gráficas de voltajes y potencias activas pre y
post contingencia. El efecto sobre el sistema es contrario a lo mostrado en la contingencia de impacto de
carga, lo cual es esperable por las mismas razones. Es decir, la modificación de factor de potencia y
consecuentemente, los valores de potencia. Además, nuevamente la tendencia de los generadores es a
mantenerse entregando una potencia activa constante, con la excepción de la central Ralco, que constituye
el generador conectado a la barra de oscilación del sistema.
Voltajes pre y post contingencia
1.02
0.98
0.96
0.94
Voltaje en p.u.
0.92
0.9
0.88
0.86
0.84
01 03 05 07 09 11 13 15 17 19 21 23 25 27 29 31
ra ra ra ra ra ra ra ra ra ra ra ra ra ra ra ra
ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar
B
Voltaje pre contingencia Voltaje post contingencia
Fig.8: Voltajes pre y post contingencias.
9
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Potencia activa pre y post contingencia en Ge-
neradores
350
300
250
Potencia activa en MW
200
150
100
50
0
GS Hidro 1 Hidro 2 Hidro 3 Hidro 4 Hidro 5 Hidro 6 Termo 2 Termo1
Ralco
MW pre contingencia MW post contingencia
Fig.9: Potencias activas en generadores pre y post contingencia.
ii. Análisis modal pre y post contingencia
En las figuras 10 y 11 se muestran los resultados del análisis modal realizado para esta contingencia.
Al igual que en el caso de la contingencia anterior, el sistema logra ser estable ante pequeñas
perturbaciones para los dos casos analizados.
10
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.10: Análisis modal pre contingencia.
11
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.11: Análisis modal post contingencia.
iii. Análisis dinámico
La respuesta transitoria de voltajes y frecuencias en barras del sistema es mostrada en la figura 12. Se
puede apreciar que, para esta contingencia, el control primario en frecuencia incluido en las máquinas
generadoras actúa de forma efectiva. Además, la frecuencia no alcanza valores prohibitivos en estado
estacionario. Dado que la frecuencia no presenta una variación contundente, los niveles de potencia y por
lo tanto, el nivel de las demás variables eléctricas no se ven afectados con una gran variación en sus
valores.
La gráfica de los ángulos rotóricos, mostrada en la figura 13, es consistente con la tendencia de la
frecuencia y el voltaje, pues es apreciable claramente que los ángulos convergen. Dado que ocurre esto, la
estabilidad transitoria para esta contingencia está asegurada.
En este punto, es importante mencionar que el control primario de frecuencia no posee la capacidad
alcanzar un valor nulo de error en estado estacionario. Por lo tanto, el valor de convergencia de la
frecuencia mostrada en la gráfica es esperable.
Fig.12: Respuesta dinámica de voltaje y frecuencia en barras.
12
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.13: Respuesta dinámica de ángulos rotóricos. La curva de ordenada cero es la referencia (Ralco).
IV. ACTIVIDAD D: Pérdida de Generación.
La central con mayor capacidad del sistema en estudio es la central Ralco. Por este motivo, la pérdida
de generación está constituida por la pérdida de esta central.
i. Flujos de carga pre y post contingencia
En la figura 14 se muestra el perfil de tensión pre y post contingencia. La pérdida de generación
constituye una disminución de tensión en las barras. Esto es esperable no solo por la relación inherente
entre potencia activa y voltaje, si no también por la estrecha relación entre los reactivos inyectados al
sistema y el voltaje a nivel sistémico. Dado que la central Ralco es la generadora asociada a la barra de
oscilación, es este elemento el que, en gran medida, suministra los reactivos necesarios para regular los
niveles de tensión. Por lo tanto, es esperable que, ante esta pérdida de fuente de reactivos, los niveles de
voltaje desciendan en su valor, llegando a niveles prohibitivos en una gran parte de las barras.
Además, los niveles de tensión en barras conectadas a los nuevos generadores se mantienen en 1 em
por unidad debido al control de tensión presente en los generadores. Todas estas barras son tipo PV.
En la figura 15 se muestran las potencias activas suministradas por los generadores. En este caso, la
potencia post contingencia que no puede ser suministrada por la central Ralco (ahora fuera de operación)
es suministrada por la central Termoeléctrica 2. Este estado de operación es absurdo, pues la cantidad de
potencia otorgada por esta central es mucho mayor que la potencia nominal de la generadora. En
13
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
particular, este equipo queda sobrecargado sobre un 500% de su capacidad. Es importante tener en
cuenta que el flujo de carga no considera elementos de control, y que, en la realidad, para un sistema
similar a este, no hay más alternativa que dejar parcial o totalmente sin suministro a uno o más consumos.
Voltajes pre y post contingencia
1.05
0.95
Voltaje en p.u.
0.9
0.85
0.8
01 03 05 07 09 11 13 15 17 19 21 23 25 27 29 31
ra ra ra ra ra ra ra ra ra ra ra ra ra ra ra ra
ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar
B
Voltaje pre contingencia Voltaje post contingencia
Fig.14: Voltajes pre y post contingencia.
Potencia activa pre y post contingencia
350
300
Potencia activa en MW
250
200
150
100
50
0
GS Ralco Hidro 1 Hidro 2 Hidro 3 Hidro 4 Hidro 5 Hidro 6 Termo 2 Termo1
MW MW
Fig.15: Potencia activa pre y post contingencia en generadores.
14
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
ii. Análisis modal pre y post contingencia
En las figuras 16 y 17 se muestra el análisis modal pre y post contingencia respectivamente. Se aprecia
que el sistema permanece estable ante pequeñas perturbaciones. Este resultado es inusual, pues ante una
contingencia de pérdida de principal fuente de potencia, lo esperable es que el sistema presente
inestabilidad. Sin embargo, se debe tener en cuenta que, para el caso en estudio, se decide conectar una
gran cantidad de equipos con control de frecuencia y tensión en terminales de barra. Además, no existe un
mayor desbalance entre la generación, y las demandas y las pérdidas. Estas características son deseables
en un sistema, y pueden ayudar a soportar contingencias.
Fig.16: Análisis modal pre contingencia.
15
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.17: Análisis modal post contingencia.
iii. Análisis dinámico
La respuesta transitoria de los voltajes y frecuencias en barras del sistema se aprecia en la figura 18.
Se observa que la frecuencia y el voltaje descienden a niveles claramente prohibitivos, pero logran
converger a un determinado valor. Esto es esperable, pues la gran cantidad de generadores con control de
frecuencia logran establecer estos valores en una determinada banda. Sin embargo, esta contingencia
posee un alto nivel de contundencia, y no es posible que el sistema se estabilice después de que esta falla
ocurra.
La figura 19 muestra la respuesta transitoria de los ángulos rotóricos. Esta respuesta es esperable,
pues el hecho de que la frecuencia y el voltaje converjan puede ser tomado como una consecuencia de la
convergencia de los ángulos rotóricos ante la falla. Cabe destaca que la curva que llega a la ordenada igual
a cero (curva celeste) corresponde a la central termoeléctrica 2. Con esto se comprueba que el estudio
estático fue, en cierta medida, predictivo, logrando anteceder ciertos niveles de voltaje y la actuación de la
central termoeléctrica 2 como la nueva fuente principal de potencia (su ángulo ahora es el mismo valor que
la referencia).
16
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.18: Respuesta dinámica de voltaje y frecuencia en barras.
17
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.19: Respuesta dinámica de ángulos rotóricos. La curva de ordenada cero es la referencia (Ralco).
V. ACTIVIDAD E: Pérdida de una línea.
Con el objetivo de comprobar la robustez del sistema, las líneas sometidas a esta contingencia son las
líneas 1.1 y 1.2, que constituyen un doble circuito y su voltaje nominal es de 500 kV.
i. Flujos de carga pre y post contingencia
En la figura 20 se muestran los niveles de voltaje pre y post contingencia. Se observa una disminución
en las caídas de tensión, la cual puede estar dada por la alteración de los flujos de potencia reactiva que
en condiciones normales de operación circulan por las líneas.
Por otro lado, en la figura 21 se muestran los niveles de potencia activa suministrados por la
generación. Se observa una disminución de potencia activa por parte de la central Ralco. Esta disminución
se debe a la redistribución de flujos de potencia provocada por la pérdida de esta línea, la cual posee una
alta capacidad de transmisión debido a su alto voltaje. Cabe destacar que, en cierta medida, la pérdida de
una línea puede comportarse como un rechazo de carga de potencia activa debido a la resistencia
distribuida de este tipo de elemento. No obstante, es importante tener en cuenta los flujos de potencia
reactiva asociada a la línea.
18
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Voltajes pre y post contingencia
1.2
0.8
Voltaje en p.u.
0.6
0.4
0.2
0
01 03 05 07 09 11 13 15 17 19 21 23 25 27 29 31
ra ra ra ra ra ra ra ra ra ra ra ra ra ra ra ra
ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar B ar
B
Voltaje pre contingencia Voltaje post contingencia
Fig.20: Voltajes en barras pre y post contingencia.
Potencia activa pre y post contingencia
350
300
250
Potencia activa en MW
200
150
100
50
0
GS Hidro 1 Hidro 2 Hidro 3 Hidro 4 Hidro 5 Hidro 6 Termo 2 Termo1
Ralco
MW pre contingencia MW post contingencia
Fig.21: Potencia activa suministrada por generadores pre y post contingencia.
19
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
ii. Análisis modal pre y post contingencia
En las figuras 22 y 23 se muestra el análisis modal pre y post contingencia respectivamente.
Nuevamente, se obtiene el resultado obtenido a lo largo de este trabajo. Esto es, que el sistema logra ser
estable ante pequeñas perturbaciones.
Fig.22: Análisis modal pre contingencia.
20
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.23: Análisis modal post contingencia.
iii. Análisis dinámico
En la figura 24 se muestra la respuesta transitoria de frecuencia y voltaje en las barras del sistema.
Ante esta contingencia, nuevamente el control establecido por los sistemas de control de frecuencia logra
mitigar la falla, siempre considerando el error en estado estacionario inherente a este tipo de estrategias.
En la figura 25 se observa la consistencia con los resultados obtenidos, pues nuevamente los ángulos
rotóricos convergen, condicionando a las variables eléctricas de interés a poseer esta misma tendencia a la
convergencia.
21
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
Fig.24: Respuesta dinámica de frecuencia y voltaje en barras.
Fig.25: Respuesta dinámica de ángulos rotóricos. La curva de ordenada cero es la referencia (Ralco
22
Universidad de Talca
Facultad de Ingeniería
Escuela de Ingeniería Civil Eléctrica
CONCLUSIONES
- La implementación de un mayor número de generadoras con control de frecuencia mejora
significativamente la respuesta transitoria ante contingencias.
- La convergencia del ángulo rotórico está directamente relacionada con la convergencia de
variables eléctricas de interés, y, por lo tanto, también concierne a la estabilidad transitoria.
- El control primario de frecuencia posee un error en estado estacionario.
- Un mayor enmallamiento y una buena relación de balance de energía contribuye (Generación,
pérdidas y demandas) a la estabilidad ante pequeñas perturbaciones.
- El ángulo rotórico está relacionado directamente con la frecuencia, y, por lo tanto, condiciona las
variables eléctricas del sistema en mayor o menor medida.
- Los estudios estáticos mediante flujo de carga pueden no ser predictivos para ciertas ventanas de
tiempo en donde la respuesta de un sistema depende críticamente de las herramientas de control.
REFERENCIAS
[1] Comisión nacional de energía. (diciembre de 2019). Norma técnica de calidad y seguridad de
servicio. [Link]
[2] Torres, F. (noviembre de 2020). Trabajo 4. Simulación de Sistemas Eléctricos de Potencia.
[Link]
[3] Godoy. B. (octubre de 2020). Tarea 3 de Simulación de sistemas Eléctricos de Potencia.
[4] Torres, F. (2020). Apunte Sistemas Eléctricos de Potencia.
23
