UNIVERSIDAD DE LA FRONTERA

FACULTAD DE INGENIERIA CIENCIAS Y

ADMINISTRACION
DEPARTAMENTO DE INGENIERIA ELECTRICA
TEMUCO

SISTEMAS ELECTRICOS DE POTENCIA II

Profesor: Manuel Villarroel Moreno

2003
 i
INDICE

CAPITULO 1 LA MAQUINA SINCRONA 1

1.1. Generalidades 1
a. Generador conectado a una barra infinita 1
b. Condiciones de puesta en paralelo 1
c. Generador alimentando una carga individual 2
1.2. Parámetros de la máquina síncrona 2
a. De régimen permanente 2
b. De régimen transitorio (sin considerar enrollados amortiguadores) 2
1.3. Circuitos equivalentes de la máquina síncrona 2
1.4. Ecuaciones de Potencia 3
a. Generador de rotor cilíndrico 3
b. Motor de rotor cilíndrico 5
c. Generador de polos salientes 6
Problemas propuestos 9

CAPITULO 2 REGULACION DE TENSION O CONTROL DE VOLTAJE 10

2.1. Introducción 10
2.2. Clasificación de las variaciones de tensión 10
2.3. Formas de regular las variaciones lentas de tensión 10
2.4. Formas de actuar de los medios de regulación de tensión 11
2.5. Regulación de tensión por inyección de potencia reactiva 11
a. Línea corta de transmisión 11
b. Línea de longitud media o larga 14
c. Sistema enmallado 15
2.6. Regulación de tensión mediante transformadores con cambio de TAP (TCT) 17
2.6.1 Sistema de transmisión radial con un TCT en el extremo transmisor 17
2.6.2 Sistema de transmisión radial con un TCT en ambos extremos de la línea 19
2.7. Regulación de tensión mediante el uso combinado de TCT e inyección de 21
Potencia Reactiva
Problemas propuestos 23

CAPITULO 3 CALCULO DE FLUJOS DE POTENCIA 25

3.1. Introducción 25
3.2. Planteamiento del problema básico 25
3.3. Modelo de representación del SEP 26
3.3.1. Tipos de barras 26
3.3.2. Representación de los elementos del SEP 27
a. Líneas 27
b. Transformadores 27
c. Generadores 27
3.4. Planteamiento matemático para un SEP de “n” barras 27
3.4.1. Ecuaciones de barras 27
3.4.2. Ecuaciones del flujo de potencias 27
3.4.3. Potencia perdida en la transmisión 28
3.4.4. Cálculo de las tensiones de barras 28
3.5. Técnicas de solución para el problema del flujo de potencias 29
3.5.1. Método de Gauss 29
3.5.2. Método de Gauss-Seidel 30
a. Caso general 30
b. Aplicación del método de Gauss-Seidel YB al cálculo de flujos de potencias 30
b.1 Sistemas con barras de carga y flotante solamente 30
 b.2 Sistemas con barras de carga, tensión controlada y flotante 31
3.5.3. Método de Newton-Raphson 35
a. Formulación general 35
b. Aplicación al cálculo de flujos de potencias 37
3.5.4. Método de Newton-Raphson desacoplado (MNRD) 41
3.5.5. Método de Newton-Raphson desacoplado rápido (MNRDR) 42
3.6. Flujo de Carga “DC” o Flujo de Potencia de Corriente Continua 45
Problemas propuestos 48

CAPITULO 4 ESTUDIO Y ANALISIS DE FALLAS 50

4.1. Introducción 50
4.2. Tipos de fallas 50
4.3. Cálculo de Cortocircuitos 51
4.3.1. Objetivos 51
4.3.2. Aproximaciones 51
4.4. Cortocircuitos trifásicos simétricos 52
4.4.1. Comportamiento de un generador en condiciones de cortocircuito trifásico 52
simétrico
a. El generador en vacío antes de producirse la falla 52
b. El generador en carga antes de producirse la falla 53
c. Concepto de potencia de cortocircuito 54
d. Algunos antecedentes relativos a la selección de interruptores 55
4.4.2. Cortocircuitos trifásicos simétricos en un SEP 57
a. Método tradicional 57
b. Cálculo sistemático (Método general) 57
4.5. Cortocircuitos Asimétricos 61
4.5.1. Componentes simétricas 61
a. Consideraciones generales 61
b. Relación entre voltajes (corrientes) de secuencia y de fase 61
c. Potencia en función de las componentes simétricas 63
4.5.2. Circuitos equivalentes de secuencia de los elementos componentes de un 63
SEP
a. Líneas 64
b. Generadores 64
c. Transformadores 66
d. Redes de secuencia 69
4.5.3. Análisis de algunos tipos de corctocircuitos 70
4.5.3.1.Cortocircuito monofásico a tierra a través de una impedancia de falla ZF 70
4.5.3.2.Cortocircuito bifásico a tierra a través de una impedancia de falla Z F 72
4.5.3.3.Cortocircuito entre dos fases a través de una impedancia de falla 75
4.5.3.4.Observaciones finales respecto a los cortocircuitos asimétricos 76
4.6. Fases abiertas 77
4.6.1. Introducción 77
4.6.2. Una fase abierta 77
4.6.3. Dos fases abiertas 81
4.6.4. Impedancias serie desequilibradas 82
Problemas propuestos 83

CAPITULO 5 ESTABILIDAD 88
5.1. Introducción 88
5.2. Estabilidad transitoria 88
5.2.1. Objetivos 88
5.2.2. Suposiciones básicas para los estudios simplificados de estabilidad transitoria 88
5.2.3. Ejemplo ilustrativo 89
5.2.4. Ecuación de oscilación de una máquina síncrona 91
5.2.5. Ecuaciones de oscilación en un sistema multigenerador 94
5.2.6. Estudio de la ecuación de oscilación de un generador conectado a una barra 95
infinita
5.2.7. Algunos métodos de solución de la ecuación de oscilación 99
a. Método punto a punto 99
b. Método de Euler 102
c. Método de Euler Modificado 103
d. Método de Runge Kutta de 4º orden 103
5.2.8. Criterio de Areas Iguales 104
5.2.9. Reducción a un sistema simple 106
5.2.10. Ecuación de oscilación de dos máquinas finitas 107
5.2.11. Factores que condicionan la estabilidad transitoria 108
5.3. Estabilidad estática o permanente 109
5.3.1. Introducción 109
5.3.2. Generador conectado a una barra infinita 110
a. Despreciando el efecto de amortiguación 110
b. Con efecto de amortiguación 111
5.3.3. Estabilidad de un sistema de “n” máquinas 111
a. Despreciando el efecto de los reguladores y del roce 111
b. Considerando el efecto del amortiguamiento 112
Problemas propuestos 112

BIBLIOGRAFIA 117
 CAPITULO 1

LA MÁQUINA SÍNCRONA

1.1. Generalidades

Desde el punto de vista del comportamiento de un generador sincrónico en un SEP se pueden identificar
las variables mostradas en la Figura 1.1. En general, el cambio de cualquiera de las variables de control, afecta a
las cuatro variables de salidas.

En un sistema de potencia, la actuación sobre el torque de la máquina motriz hace posible mantener un
balance exacto entre la potencia activa generada en el sistema y la potencia activa consumida por las cargas más
las pérdidas. Este balance permite que el sistema trabaje a una frecuencia constante.

Por otra parte, actuando sobre la corriente de excitación de cada generador, es posible mantener un
balance exacto entre las potencias reactivas generadas y las del consumo y pérdidas. Este balance permite
mantener una tensión constante (en módulo) en las barras del sistema.

a. Generador conectado a una barra infinita: La tendencia actual es trabajar con un sistema eléctrico de
transmisión y distribución que interconecte las diferentes centrales generadoras y los puntos de consumo, de
modo que una sola central y más aún, una sola máquina representa un pequeño porcentaje de la potencia total del
sistema. Evidentemente esta máquina no estará capacitada para alterar ni el voltaje ni la frecuencia del sistema
eléctrico. En el límite se puede considerar que el sistema mantiene el voltaje y la frecuencia invariables, lo que se
puede asimilar a una máquina que tiene cero impedancia interna e inercia rotacional infinita. Se habla entonces
de una “barra infinita” (Figura 1.2).

I exc
Variables P Variables Barra
Máquina
de Q
de
Tm Síncrona V
Control F Salida
/V/=cte, f=cte
Figura 1.1.- Representación esquemática de una Máquina Síncrona
Figura 1.2.- Barra infinita

Un gran SEP puede considerarse como una barra infinita para muchos fines prácticos. En este caso
entonces, dado que el voltaje y la frecuencia se mantienen constantes, las variables de salida se reducen a las
potencias activa P y reactiva Q donde la corriente de excitación I exc influye fundamentalmente en Q y el torque
motriz Tm influye básicamente en P; es decir, hay un acoplamiento débil entre I exc y P y entre Tm y Q.

b. Condiciones de puesta en paralelo: Para operar dentro de un sistema eléctrico, un generador síncrono
debe trabajar normalmente en paralelo con otros generadores, unidos a una barra de generación común, si
pertenecen a una misma central o separados por alguna impedancia (líneas, transformadores) si pertenecen a
centrales diferentes. Tratándose de alternadores polifásicos, las condiciones de puesta en paralelo son las
siguientes:

 Módulo de la tensión en bornes igual al existente en la barra. Se actúa sobre la corriente de excitación

 Frecuencia igual a la de la barra. Se actúa sobre el torque motriz

 Igual secuencia de fases

 Igualdad en la forma de onda

 2
c. Generador alimentando una carga individual: Supongamos por simplicidad que la carga es de tipo
impedancia estática y se aumenta T m. En este caso, al aumentar T m, aumenta la velocidad y por lo tanto varían
tanto la frecuencia como el voltaje y en consecuencia, varían también P y Q.

1.2. Parámetros de la máquina síncrona

a. De régimen permanente

Rotor cilíndrico: Impedancia síncrona Zs=Rs+jXs, habitualmente Zs=jXs

Rotor de polos salientes: Reactancia de eje directo: X d
Reactancia de eje en cuadratura: Xq

b. De régimen transitorio (sin considerar enrollados amortiguadores)

Reactancia transiente de eje directo: X 'd
Constante transiente de tiempo: 'd

Adicionalmente, para el estudio del comportamiento de la máquina síncrona en régimen desbalanceado

es necesario introducir los parámetros impedancia de secuencia cero (Z0), positiva (Z1) y negativa (Z2) que se
estudiarán posteriormente, así como la reactancia subtransiente ( Xd" )

1.3. Circuitos equivalentes de la máquina síncrona

Los circuitos equivalentes dependen del tipo de máquina (rotor cilíndrico o de polos salientes), del tipo
de estudio que se desea realizar y del grado de precisión deseada. Se pueden clasificar en dos grupos: aquellos
utilizados en los estudios de régimen permanente (balanceado o desbalanceado) y los utilizados en los estudios
de régimen transitorio, particularmente cortocircuitos y estabilidad. Las Figuras 1.3. a) y b) muestran los
circuitos equivalentes de una máquina síncrona de rotor cilíndrico, usuales en los estudios de régimen
permanente. En a) se considera la resistencia síncrona y en b) se desprecia. Para una máquina de polos salientes
no es posible establecer un circuito equivalente simple y su comportamiento se estudia en base a su diagrama
fasorial d-q.
Rs jXs jXs

+ + + +
E V E V
- -
- -

a) b)
Figura 1.3.- Circuitos equivalentes de una máquina síncrona para estudios de régimen permanente

Los circuitos de la Figura 1.4. a) y b) se usan en los estudios de cortocircuitos simétricos: a) se emplea
en aquellos casos en que se supone la máquina en vacío antes de producirse el cortocircuito. La tensión E, en
este caso, es la de vacío en los terminales de la máquina. Si la corriente antes de producirse el cortocircuito es
significativa, se emplea el circuito de la Figura 1.4.b). La elección de E' y E'' así como X 'd y X "d en ambos
casos, depende del instante de tiempo en que se evalúa la corriente después de ocurrida la falla.

Es conveniente hacer notar que los circuitos de la Figura 1.4. se emplean tanto para máquinas de rotor
cilíndrico como de polos salientes. Esto se debe a que en condiciones de cortocircuito, las corrientes son
prácticamente reactivas y en consecuencia la componente en cuadratura de la corriente (I q) en las máquinas de
polos salientes se puede despreciar sin cometer gran error.

En los estudios simplificados de estabilidad transitoria, una máquina síncrona de rotor cilíndrico o de
polos salientes, se representa por el circuito equivalente mostrado en la Figura 1.4 b), considerando E' y X ' d.
 jX'd jXd'

+ jX'd' + + jXd' +
E'
E V
V E ''
- - - -

a) b)

Figura 1.4.- Circuitos equivalentes de una máquina síncrona para estudios de fallas simétricas

1.4. Ecuaciones de Potencia

a. Generador de rotor cilíndrico: Consideremos un generador suministrando potencia directamente a un

consumo. El circuito equivalente por fase y el diagrama fasorial se muestran en la Figura 1.5.

jXs
I E
+
+ Con-
E V  jI Xs
sumo

-
I V
-
a) b)

Figura 1.5.- Generador síncrono de rotor cilíndrico: a) Circuito equivalente; b) Diagrama fasorial

La potencia suministrada por el generador al consumo es:

S& g  Pg  jQ g  V& I  (1.1)

Si V&  V0º y E&  E; del circuito se tiene que:

&
&I  &E V  E
V j (1.2)
jXs Xs

Expresando la corriente en la forma rectangular, reemplazando en (1.1) y separando en sus partes real e
imaginaria se obtiene:

Pg  E V EV V2
sin Qg  cos   (1.3)
 Xs Xs
Xs
que corresponden a las potencias activa y reactiva, respectivamente.

 Potencia activa

La primera ecuación de (1.3) se denomina Ecuación Potencia Activa-Angulo y se puede visualizar en la

Figura 1.6, para E y V constantes. Como en estas condiciones, P g es sólo función de , significa que si se supone
velocidad constante y Tm aumenta, el rotor avanza un cierto ángulo respecto a su posición original, lo que
implica que el fasor E se aparta del fasor V, ya que su fase depende de la posición del rotor. En consecuencia 
aumenta y Pg aumenta. En el gráfico se puede apreciar también que en el caso en que la máquina funcione como
motor, el ángulo  es negativo, es decir la potencia activa es negativa (según la referencia considerada), esto es,
la potencia activa llega a la máquina.
 Pg

Motor
Generador 

Figura 1.6.- Gráfico Potencia Activa-Angulo para una máquina síncrona

Por otra parte, para E y V constantes, la máxima potencia activa que el generador puede entregar se tiene
cuando el ángulo  es de 90º, lo que representa el límite de estabilidad permanente de la máquina. Si el ángulo
sigue aumentando, la potencia activa disminuye y el generador pierde el sincronismo.

Ejemplo 1.1. Un generador de rotor cilíndrico que tiene una reactancia síncrona del 50%, está conectado a una
barra infinita a través de una línea corta de transmisión cuya reactancia es del 25%. La tensión en bornes de la
máquina es la nominal, cuando la potencia recibida por la barra infinita es de (90+j40 ) MVA. Datos en base
común 100 MVA. Determinar:
a) La potencia activa máxima que el generador le puede entregar a la barra infinita.
b) La potencia reactiva que el generador entrega en las condiciones anteriores.

Solución:

a) Potencia activa máxima

En el circuito de la Figura 1.7 se tiene: j0,5 j0,25

+ I +
V&   V 0º ; V&  1 +
E V V
La corriente queda: -
* - -
&  S  0,9  j0,4
I 
V0º Figura 1.7
 V 
Para calcular Pg se debe conocer E y V. A partir del circuito se puede plantear la siguiente ecuación:

(0,9  j0,4)
V&  j0,25   V&  ; reemplazando valores y amplificando por V&  se puede escribir:
V& 

1  V  0,1  j0,225  V 2

 separando en parte real e imaginaria quedan las dos ecuaciones siguientes:
;

V  cos   0,1  V 2 Elevando al cuadrado ambas ecuaciones y sumando se tiene:


V  sin   0,225

V 4  0,8  V 2  0,0606  0; cuyas soluciones son: V  0,2911 &  0,84570º

  
  V 
 0,8457
 
A partir del valor de V&  y del circuito se puede determinar E&
y luego Pgmáx como sigue:

0,9  j0,4
E&  j0,75  0,84570º  1,441633,62º
 0,84570º

E  V 1,4416  0,8457
P
P   sin 90º   =1,6255 (pu)
gmáx
g máx
XT 0,75

b) Potencia reactiva

EV V 2
Qg  0,8457
X T  cos 90º  Qg=-0,9536 (pu)
XT  
2

0,75

 Potencia Reactiva

La segunda ecuación de (1.3) representa la potencia reactiva que la máquina puede generar (consumir), y
se puede escribir como sigue:

Qg  V (E cos   V) (1.4)
Xs
Respecto a esta ecuación, se puede observar que:

- Si E cos   V, entonces Qg  0. Esto significa que la máquina genera potencia reactiva y actúa desde el
punto de vista de la red, como un condensador. Generalmente, la desigualdad anterior se cumple cuando el
generador trabaja sobreexcitado (Figura 1.5.b); y por lo tanto, la carga servida por el generador es inductiva.

- Si E cos   V, entonces Qg  0. Esto I

E
significa que el generador consume jI X s
potencia reactiva "desde la red" y actúa  
como una reactancia shunt vista desde ésta. V
Esta condición se cumple cuando la
máquina trabaja subexcitada y tal como se Figura 1.8.- Diagrama fasorial de un generador subexcitado
muestra en el diagrama fasorial de la Figura
1.8, la corriente adelanta al voltaje, lo que
significa
que la carga servida por el generador es de tipo capacitivo.

- Si =0, significa que la potencia activa es cero y en la expresión (1.4) se tiene:

Qg  V (E  V) (1.5)
Xs
Ecuación que en el caso en que la tensión V del sistema es constante, depende sólo de la corriente de
excitación, lo que da origen a los “condensadores síncronos” (EV) o “reactores síncronos” (EV).

b. Motor de rotor cilíndrico: Tal como se indicó anteriormente, en el caso en que 0, la máquina absorbe
potencia activa. La potencia reactiva en cambio, es independiente del signo de  y por lo tanto, para ella se
cumplen las mismas condiciones dadas para el generador, esto es:

- El motor entrega potencia reactiva si está sobreexcitado, pero consume corriente en adelanto (Figura 1.9.
b); esto es, se comporta como un consumo capacitivo.
- El motor absorbe potencia reactiva si está subexcitado, pero consume corriente en atraso (Figura 1.9.c); o
sea se comporta como un consumo inductivo.

jXs IM I
IM V
+ I +  V
I j Xs I M  
E V 
jI Xs
- E
- E IM

a) b) c)

Figura 1.9.- Motor síncrono de rotor cilíndrico: a) Circuito equivalente, b) Diagrama Fasorial sobreexcitado, c)
Diagrama Fasorial subexcitado

La Figura 1.10 muestra un diagrama con la posición de la corriente para las cuatro condiciones de
operación descritas anteriormente. En todas ellas se supone que la corriente sale de la máquina.

De acuerdo con las convenciones más usadas, cuando la corriente tiene una proyección positiva sobre el
voltaje, la máquina entrega energía (generador) y viceversa. Si se conectara un fasímetro, éste indicaría los
ángulos mostrados en la figura, suponiendo siempre que la corriente sale de la máquina. Sin embargo, las
designaciones mostradas para el trabajo como motor (atraso o adelanto) corresponden más exactamente a estas
corriente giradas 180º.

P
Generador, P > 0 Generador, P > 0
Subexcitado, Q < 0 V
Corriente en adelanto Sobreexcitado, Q > 0
Corriente en atraso
I V
I  

-Q Q I
 

I
I
Motor, P < 0 Motor, P < 0
Subexcitado, Q < Sobreexcitado, Q > 0
0 Corriente en Corriente en adelanto
atraso
-P
a) b)

Figura 1.10.- a) Diagrama que muestra las cuatro condiciones de operación de una máquina síncrona, b)
Diagrama unilineal con el sentido adoptado para la corriente

Otra forma de representar las condiciones de operación de una máquina síncrona, considera al voltaje V
positivo en el eje de abcisas, lo que es equivalente a girar en 90º hacia la derecha el diagrama mostrado en la
Figura 1.10.

c. Generador de polos salientes: Como ya se dijo, no es posible representar el generador de polos

salientes a través de un circuito equivalente. Por ello, se trabaja con su diagrama fasorial. La Figura 1.11 muestra
el diagrama fasorial de un generador de polos salientes, a partir del cual se pueden obtener las ecuaciones de
potencias activa y reactiva.
 d

Iq E
 q
 j Xq I q

Id
V j Xd I d
I

Figura 1.11.- Diagrama fasorial de un generador de polos salientes

 Potencia activa

La potencia activa es de la forma Pg=V I cos  es decir V multiplicado por la proyección de I sobre V.
Por lo tanto, proyectando las componentes de la corriente, se tiene:

I cos   Iq cos   Id sin  (1.6)

Por otra parte, del mismo diagrama se puede escribir:

E  V cos 
I d V sin 
Xd I q (1.7)
Xq

Por lo que la potencia activa Pg queda de la forma:

P VE X d  Xq
g  sin   V 2 sin 2
Xd (1.8)
2 Xd X q

La Figura 1.12 muestra la curva de P g=f() dada por la ecuación (1.8). En líneas de segmentos están
dibujadas las curvas correspondientes al primer término (la fundamental) y al segundo término (la segunda
armónica).

El segundo término de (1.8) se denomina componente de reluctancia o de saliencia y es pequeño

comparado con el primero (10 a 20% usualmente). No depende de la excitación E y por ello existe aunque la
corriente de excitación sea nula. Para corrientes de excitación grandes no se comete un error importante al
despreciarlo.

Si la máquina está conectada a un sistema relativamente grande (barra infinita), como ocurre en la
mayoría de los casos, V será constante o tendrá un rango de variaciones posibles bastante estrecho; E, que
depende de la corriente de excitación, también se mantendrá constante. Por lo tanto, para fines prácticos se puede
afirmar que Pg depende sólo de . Si   0, Pg es positiva, es decir, la máquina genera potencia activa si E
adelanta a V. A un aumento de  corresponde un aumento no proporcional de la potencia, ya que el denominado
“coeficiente de dP d no es constante.
sincronización” g

Existe un , para el que, dPg d  0, o sea para el cual se obtiene la máxima potencia activa compatible
con los valores de V y E adoptados. Si  sigue aumentando, Pg se reducirá y la máquina perderá el sincronismo.
 Pg

Motor


Generador

Figura 1.12.- Curva Potencia-ángulo para una máquina de polos salientes

 Potencia Reactiva

Recordemos que Qg=V I sin , es decir es el producto del voltaje por la componente de la corriente, en
cuadratura con él. Por lo tanto, a partir del diagrama fasorial y considerando las ecuaciones (1.7), la potencia
reactiva queda de la forma:
VE X d  Xq
Q  cos   Xd  Xq
V 2 cos 2  V2 (1.9)
g
X 2 X d Xq 2 Xd X q
d

Ejemplo 1.2. Un motor síncrono de polos salientes con Xd=80% y Xq=50%, está conectado a una barra infinita
cuyo voltaje es el nominal. Los datos están en la base propia de la máquina. Despreciando las pérdidas calcular:
a) El valor mínimo de E para que la máquina permanezca en sincronismo a carga nominal.
b) La potencia reactiva (consumida o generada).

Solución:

a1) Límite de estabilidad:  dPg d  0 ; es decir:

dPg V E X d  Xq 1 E 0,8  0,5 2
 cos   V 2 cos 2   cos   1  cos (2  )  0
d Xd Xd X q 0,8 0,8  0,5

despejando E de la ecuación anterior se tiene: cos (2  (a)

E  0,6 
)
cos 
a2) Carga nominal:  Pg  1 en ecuación (1.8), de donde se obtiene:
1 E
1  sin 
0,8  0,5
12  sin (2  )  0 ; es decir:
 2  0,8 
0,8 0,5

 1  1,25  E  sin   0,375  sin (2 (b)

 )

Introduciendo (a) en (b) y ordenando se puede escribir:

1
F()  0,75  tg   0,375  tg (2  ) 0
 cos (2  )
Empleando un método iterativo (Newton-Raphson, por ejemplo), se obtiene que con   60,5º F()  0,008 .

; Reemplazando  en la ecuación (a) se obtiene finalmente: E=0,6276 (pu)

b) Potencia reactiva: Reemplazando valores en la ecuación (1.9), se tiene:

1 0,6276
Q   cos (60,5º )
0,8  0,5 0,8  0,5
12  cos (2  60,5º ) 12  QM= 1,4318 (pu)
 2  0,8   2  0,8 
M
0,8 0,5 0,5

Problemas propuestos

1.1. Un generador síncrono tiene una impedancia síncrona de 10  y una resistencia síncrona de 1,2  (Valores
por fase). La tensión generada es 1,1 kV (por fase) y se conecta a una barra infinita de 1 kV (por fase).
Determinar la corriente y factor de potencia cuando se tiene la máxima potencia activa generada.

1.2. Un generador síncrono de polos salientes con X d=1,5 (pu) y Xq=1 (pu) y resistencia despreciable, se conecta
a una barra infinita cuyo voltaje es el nominal, a través de una línea que tiene una reactancia de 0,3 (pu).
Determinar la potencia activa de salida para un ángulo de torque de 30º, si E es 1,4 veces el voltaje nominal
terminal del generador. Calcular además, el coeficiente de sincronización en estas condiciones. Todos los valores
en (pu) están en base 75 MVA. El generador es de 75 MVA, 11 kV.

1.3. En el problema 1.2 determinar el valor máximo de la potencia activa que el generador le pueda entregar a la
barra infinita. En estas condiciones, ¿cuál es el valor de la potencia reactiva?

1.4. ¿Qué porcentaje de su potencia nominal podrá desarrollar sin pérdida de sincronismo un motor síncrono de
polos salientes cuando se le aplica su tensión nominal siendo la corriente de excitación igual a cero?. X d=0,8;
Xq=0,5, ambos en por unidad en base propia. Determine también, la corriente en el inducido a la potencia
máxima.

1.5. Un motor síncrono, de polos salientes, tiene Xd=0,8; Xq=0,5; en tanto por unidad base propia y está
conectado a una red de potencia infinita cuyo voltaje es el nominal. Depreciando las pérdidas, determine la
excitación mínima, para que la máquina se mantenga en sincronismo con el torque de plena carga.

1.6. Un motor shunt de corriente continua está acoplado mecánicamente a un generador trifásico síncrono de
rotor cilíndrico. El motor se conecta a una red de tensión continua constante de 230 V y el generador a una red
trifásica de frecuencia y tensión constantes, siendo ésta de 230 V entre fases. La máquina síncrona, de cuatro
polos, conectada en Y, es de 25 kVA a 230 V, y su reactancia síncrona es de 1,6 ohm/fase. La máquina de
continua, también de cuatro polos, tiene los valores nominales de 25 kW y 230 V. Se desprecian todas las
pérdidas.
a. Si el conjunto de ambas máquinas trabaja como un grupo motor-generador recibiendo potencia de la red de
continua y suministrándola a la de alterna, ¿cuál es la tensión simple (entre fases y neutro) inducida por el
campo inductor de la máquina de alterna cuando trabaja a sus kVA nominales y factor de potencia unitario?
b. Dejando la excitación de la máquina de alterna como en el apartado anterior, ¿qué ajustes deberán realizarse
para reducir a cero la transferencia de potencia entre los sistemas de alterna y de continua? En estas
condiciones de transferencia nula, cuál es la corriente en el inducido de la máquina de continua y en el
inducido de la de alterna?
c. Dejando la excitación de la máquina de alterna como en los dos apartados anteriores, ¿cómo podrá regularse
para transferir 25 kW de la red alterna a la de continua? En estas condiciones ¿cuál es la corriente en el
inducido de la máquina de continua? Cuáles son la magnitud y fase de la corriente en la máquina alterna?
 CAPITULO 2

REGULACIÓN DE TENSION O CONTROL DE VOLTAJE

2.1. Introducción

Un sistema de potencia bien diseñado debe ser capaz de entregar un servicio confiable y de calidad.
Entre los aspectos que caracterizan una buena calidad de servicio están, la adecuada regulación del voltaje así
como de la frecuencia. El Control de Voltaje tiene como objetivo mantener los niveles de tensión dentro de
límite razonables. El problema; sin embargo, es diferente según se trate de una red de distribución o una de
transmisión.

En una red de transmisión se pueden admitir variaciones de tensión mayores que en una red de
distribución, ya que no existen aparatos de utilización directamente conectados a ella. Por lo tanto, dentro de
ciertas limitaciones, no hay mayores inconvenientes en que la tensión en un punto dado de la red de
transmisión varíe dentro de límites relativamente amplios, generalmente  10% como máximo de un valor que
habitualmente es diferente del nominal.

En las redes de distribución, las variaciones de tensión están limitadas por las características de los
consumos. Estos sólo funcionan adecuadamente con tensiones cercanas a la nominal y admiten variaciones
lentas que no sobrepasen un  5% en aplicaciones térmicas (cocinas, lámparas, calentadores) y un  8% en el
caso de motores, lavadoras, receptores de radio y televisión, etc. Con tensiones muy altas habrá calentamiento
y menor vida útil. Si son muy bajas habrá mal rendimiento, malas características de torque (motores), etc.

2.2. Clasificación de las variaciones de tensión

Según sus características, las variaciones de tensión se pueden clasificar en:

– Variaciones lentas: Tanto previsibles (periódicas), originadas en los cambios periódicos de los
consumos que presentan máximos a ciertas horas del día y mínimos en otras; como aleatorias, debidas
a las conexiones y desconexiones de los consumos, que pueden ocurrir en cualquier momento.

– Variaciones bruscas: Tanto regulares como aleatorias (“pestañeos”), debidas a los “golpes de
corriente” causados por el funcionamiento intermitente de equipos tales como refrigeradores,
ascensores, soldadoras, etc.

– Caídas de tensión: De breve duración (desde fracciones de segundo hasta algunos segundos) y de
amplitudes muy variables (hasta un 100% de la tensión). Su efecto es casi equivalente al de una
interrupción de servicio

En este capítulo, se analizará exclusivamente el control de las variaciones lentas de tensión en las redes
de transmisión y distribución.

2.3. Formas de regular las variaciones lentas de tensión

La regulación lenta de tensión tiene por fin mantener el módulo de la tensión en todo el sistema en el
mejor valor posible. Los métodos más empleados son:

– Inyección (absorción) de potencia reactiva: permite modificar la potencia reactiva circulante en el

sistema, que es una importante causa de variación de la tensión. Se consigue con el empleo de
condensadores estáticos, compensadores síncronos, reactores y los generadores de las centrales.

– Inserción de una tensión serie adicional: para compensar la caída que se desea regular. Se consigue,
por ejemplo, con los transformadores con derivaciones, que permiten variar discontinuamente la razón
de transformación. Las derivaciones pueden ser operables en vacío (más baratas) o bajo carga.
 11
– Modificación de la reactancia: para mantener constante la caída longitudinal ZI. Se consigue por
ejemplo, usando conductores fasciculados, empleando condensadores serie, colocando líneas en
paralelo o disminuyendo el largo de las líneas, acercando los transformadores de distribución a los
consumos.

2.4. Formas de actuar de los medios de regulación de tensión

– Regulación continua: Reguladores de inducción y compensadores síncronos.

– Regulación cuasi-contínua: Cambiadores de derivación bajo carga de los transformadores.

– Regulación intermitente: Condensadores estáticos.

– Regulación fija: Condensadores serie y cambiadores de derivación en vacío.

2.5. Regulación de tensión por inyección de potencia reactiva

a. Línea corta de transmisión: Con el objeto de establecer algunos conceptos básicos consideraremos
en primer lugar, el caso de una línea corta de transmisión y su diagrama fasorial correspondiente, tal como se
muestra en la Figura 2.1. PR y QR son las potencias activa y reactiva que llegan a la carga a través de la línea.

PT + jQT PR + jQR VT E
R + jX
I

VT VR Carga
A D jXI
 B
R VR C
RI
I V
a) b)

Figura 2.1.- Línea corta de transmisión: a) Circuito equivalente; b) Diagrama fasorial

De la Figura 2.1 y considerando que  es un ángulo pequeño, se puede escribir:

V& T   VT  VR  V (2.1)
V& R

pero, V=AD+DC, con: AD=RI cos R y DC=XI sin R, por lo que se obtiene finalmente:

RPR  XQ R
V  VR (2.2)

De la misma forma:

XPR 
EC  (2.3)
RQR
VR

Por lo que:
1 
EC 
  tg  V  V  (2.4)
 R 

Las ecuaciones (2.2) y (2.4) muestran claramente que el transporte de P R y QR desde el extremo
transmisor T al receptor R, va acompañado de una caída de tensión V y de un desfase  entre las tensiones de
ambos extremos de la línea. Debido a que usualmente, en los sistemas de transmisión y generación R X, el
término dominante en la ecuación (2.2) es XQ R, o sea la potencia reactiva transferida desde T a R es la
causante principal de la caída de tensión V. Por la misma razón anterior, el término dominante en (2.4) es
XPR, de donde se concluye que el desfase  es decisivo en el valor de la potencia activa suministrada al
consumo.

Según lo planteado, para reducir la caída de tensión V es necesario evitar (disminuir) el transporte de
potencia reactiva por la línea, es decir, la potencia reactiva debería ser suministrada en lo posible, en el mismo
punto donde será consumida.

Supongamos que se desea mantener constante V para cualquier valor de PR. Entonces, a partir de la
ecuación (2.2), QR debe variar según la ecuación (2.5), para cumplir con la condición señalada.
R
Q K P
(2.5)
R R
X
Donde el término K=(V/X)VR Q
es constante, puesto que para VT y V R
constantes, entonces, VR es constante y K
por lo tanto, la potencia reactiva Q R
enviada al consumo a través de la línea
debe variar tal como se muestra en la
Figura 2.2.
Q
R0
Controlando la potencia reactiva
QR en la forma mostrada en la Figura QR1
2.2, es posible mantener la tensión VR
constante para diversos valores de PR
potencia activa PR. Este control se hace, P P XK/R
R R
evidentemente en el extremo receptor y 0 1
el procedimiento se conoce con el Figura 2.2.- Variación de la potencia reactiva QR
nombre de “regulación de tensión por entregada por la línea en función de la potencia activa del
inyección de potencia reactiva”. consumo, para mantener VR constante

Conviene insistir en que la potencia QR correspondiente a la ecuación (2.5) representa solamente la

potencia reactiva que llega al consumo a través de la línea, cuyo valor, en caso de producirse una variación de
PR, será diferente al correspondiente al consumo ya que su ley de variación para V constante está definida por
la ecuación (2.5)

Con el objeto de precisar lo expuesto, consideremos:

QRO: Potencia reactiva requerida por el consumo

PRO: Potencia activa inicial del consumo
QR1: Nuevo valor de potencia reactiva que se debe suministrar al consumo a través de la línea.
PR1: Nuevo valor de potencia activa requerida por el consumo.

A partir de (2.5) se puede escribir:

R
QR 0  K  P
R0
X
(2.6)
R
Q K P
R1 R1
X
Según sea el nuevo valor PR1, se pueden presentar dos casos: a
2
a1) PR1  PR0; entonces QR1 QR0 )
 PR1  PR0; entonces QR1 QR0

La potencia que se debe inyectar en el consumo será:

Q  QR 0  QR1

(2.7)

El valor de Q será positivo o negativo, según corresponda al caso a1) o a2) y se puede considerar en
general como una potencia inyectada según se muestra en la Figura 2.3.
PR + jQR PR+j(QR+ Q)
V VR
T
Línea
CONSUMO
Q

Figura 2.3.- Representación de Q como fuente de potencia reactiva

Ejemplo 2.1. El sistema de la Figura 2.4, alimenta una carga que en demanda máxima (D máx) es de 45 MVA,
cos  =1. Durante la demanda mínima (D mín) la carga es de (10+j2) MVA. Si la tensión en el extremo
transmisor puede tomar como valores extremos 110 y 100% respectivamente, determinar la compensación de
potencia reactiva necesaria en el extremo receptor durante D máx y Dmín, si se desea mantener en dicho extremo
una tensión de 100 y 95% respectivamente. SB=100 MVA

Solución:

a) Para demanda máxima: S& L  0,45  j0,0 G

T R
0,38+j0,6 (pu) S R SC
(pu) ; VT=1,1 (pu); VR= 1,0 (pu)

Aplicando la aproximación de una línea corta (Ecuación

2.2) se tiene: Figura 2.4 S L Comp.

V 
VT  V  RPR  XQR  1,1  1  0,38  0,45  0,6  QR
R V R 1

De donde se obtiene: QR=0,1183 (pu)  QC=QLQR=0(0,1183)=0,1183 (pu)  QC=11,83 MVAR

(inyectar).

b) Para demanda mínima: S& L  0,1  j0,02 (pu) ; VT=1,0 (pu); VR= 0,95 (pu)

Aplicando la ecuación (2.2) se tiene:

Q R  0,95  0,38  0,1  0,6   Q =0,01583 Q =Q -Q =0,02-0,01583=0,00417 (pu)

1,0 C L R
0,95
R

 QC=0,417 MVAR (inyectar).

b. Línea de longitud media o larga: El estudio de la regulación de tensión en una línea media o larga,
mediante inyección de potencia reactiva, se realiza usualmente, empleando las ecuaciones de potencia activa y
reactiva en función de las constantes generalizadas A-B-C-D ya vistas, que se resumen a continuación.
D 2 1
P V cos (  )  VV
B  T cos (  )
T   RT
D  2  1B 
Q V sin (  )  V V sin (  )
T B  T   RT
 1    B
A 2 (2.8)
P V V cos (  )  V cos (  )
R  B R T  R
 1 AB 2
Q  V V sin (  )  V sin (  )
R  B R T  R
  B

Ejemplo 2.2. Repetir el Ejemplo 2.1, empleando parámetros ABCD y las ecuaciones de potencia del extremo
receptor.

Solución: Los parámetros ABCD de la línea son:

A=D=10º; C=0; B=Z=0,38+j0,6=0,710257,65º

a) Para demanda máxima: VT=1,1 (pu), VR= 1 (pu); PL=0,45 (pu) y QL=0

Las ecuaciones correspondiente al extremo receptor se obtienen de (2.8), es decir:

1 A 2
P V V cos (  )  V cos (  )
R  B R T  R
 1 AB 2
Q  V V sin (  )  V sin (  )
R  B R T  R
  B

A partir de PR despejamos cos (), donde:

PR B  AVR2 cos (  0,45  0,7102  11 2  cos(57,65º0º  0,77698

cos (  )   =18,64º
) VR VT ) 11,1


reemplazando valores en la ecuación de QR:

 1 
Q  1 1,1 sin (57,65  18,64º ) 1  2
 1 sin (57,65º0º )  0,2146 (pu)
R  0,7102   0,7102 
   

De donde se obtiene que: QC=QL-QR=0-(-0,2146)=0,21146 (pu) QC=21,46 MVAR (inyectar).

b) Para demanda mínima: VT=1,0 (pu), VR= 0,95 (pu); PL=0,1 (pu) y QL=0,02

A partir de PR despejamos cos (), de donde:

cos (  ) PR B  AVR2 cos (  0,1 0,7102  1 0,952  cos (57,65º0º  0,5831  =3,32º
  ) 0,95 1,0
) V V R T
reemplazando valores en la ecuación de QR:
 1 
Q  0,95 1,0  sin (57,65  3,32º ) 1 
 0,952  sin (57,65º0º )  0,01315 (pu)
R  0,7102   0,7102 
   

De donde se obtiene que: QC=QLQR=0,020,01315=0,00685 (pu) QC =0,685 MVAR (inyectar).

c. Sistema enmallado: En un caso general, interesa estudiar la influencia de la inyección o extracción de

potencia reactiva en una barra cualquiera del sistema, sobre la tensión de dicha barra. El estudio parte de la
premisa que existe un relación definida entre las potencias P, Q y la tensión V de la barra, relación que puede
escribirse como V=V(P,Q), de donde:
V V
dV  dP  dQ
(2.9)
P Q
o bien:
dP dQ
dV   (2.10)
(P V) (Q V)

Es decir, la variación de tensión en un punto cualquiera p de un sistema, debido a cambios dP y dQ en

las potencias activa y reactiva allí entregadas, está completamente determinada si se conocen los coeficientes
(P/V) y (Q/V). De los dos, el (Q/V) es el más importante ya que orienta sobre la amplitud de la
variación de potencia reactiva que es necesario producir para provocar una variación determinada de tensión en
el punto considerado.

Para determinar estos coeficientes, consideremos que la carga S = P+jQ existente en el punto p del
sistema de la Figura 2.5 a) se aumenta en un consumo puramente inductivo Q, de valor muy pequeño (y que
en el límite se hace tender a cero). Este hecho modifica la tensión en p, que de V pasa a V + V.

Pi + jQi
P+ j(Q + Q)
Pn+ jQn
ZT
+
P2+ jQ2 Q +
V0 V + V
p
-
jX -
P1 + jQ1 S=P+jQ
a) b)
Figura 2.5.- Determinación del coeficiente (Q/V). a) Diagrama unilineal; b) Circuito equivalente

Según el teorema de Thevenin, el sistema puede ser representado por el circuito de la Figura 2.5 b),
donde V0 es el voltaje que existía en p cuando no había consumo y ZT es la Impedancia equivalente de
Thevenin del sistema vista desde p. En el circuito de la Figura 2.5 b) se puede aplicar lo planteado en las
ecuaciones (2.1) y (2.2); que en este caso quedan:

V   R T P  XT (2.11)
V
Q
V0 V
es decir:

V0 V  V 2  R T P  X T Q 
(2.12)
0
 La ecuación (2.12) es una función implícita y por lo tanto, como en este caso, sólo se produce un
cambio en Q, (P permanece constante) se puede derivar parcialmente respecto a V y se obtiene:

Q V0  2V

(2.13)
V XT

Si hubiera un cambio en P (P) con Q constante, se encuentra que:

P V0  2V
 (2.14)
V RT

En el caso en que no exista consumo previo S en el punto p, V0=V, por lo que (2.13) y (2.14) quedan:

Q V
V   X T (2.15)
P V
V   RT (2.16)

Valores típicos de (Q/V) van desde los 3 a 15 [MVAr/kV]. El signo menos indica que si Q (ó
P) es positivo, lo que significa que la potencia que sale de p aumenta, entonces la tensión en el nudo p
disminuye. Por lo tanto, la conexión de un reactor produce una baja de tensión mientras que la conexión de un
condensador produce un aumento en la tensión.

Es interesante hacer notar que (P/V) y (Q/V) tienen las dimensiones de corriente. Por otra parte; la
corriente de cortocircuito trifásico en p (I CC), (como se verá en el Capítulo 4), se puede calcular directamente
de la Figura 2.5, despreciado RT y la corriente previa a la falla y vale:
V
I 
(2.17)
CC
XT
Las expresiones de (Q/V) e ICC (en módulo) son prácticamente iguales. La diferencia está en que la
corriente de cortocircuito se calcula con las reactancias transitorias de las máquinas; mientras que (Q/V)
considera reactancias permanentes o sólo las reactancias hasta aquellas barras en que la tensión permanece
constante por efecto de la acción de los reguladores de tensión. Para nudos algo alejados de las máquinas, casi
no habrá diferencias entre ambas definiciones. Por ello se suelen considerar como aproximadamente iguales.

Para el cálculo real de (Q/V) no se consideran las impedancias de las cargas y se suponen constantes
las tensiones en bornes de las máquinas síncronas.

Por otra parte, aunque hasta el momento, sólo se ha considerado inyectar potencia reactiva, es evidente
que, dependiendo de las condiciones de carga en el sistema, habrá que absorber potencia reactiva en algunos
casos, para mantener la tensión dentro de los límites prefijados. Generalmente es necesario absorber potencia
reactiva en las horas de poco consumo, lo que se debe fundamentalmente a la potencia reactiva “generada” por
las líneas de transmisión que operan a niveles de tensión elevados o por redes de cables.

Ejemplo 2.3. En el sistema de la Figura 2.6, las barras A A B

y C se mantienen a tensión nominal de 220 kV. La barra X=50  M X=50
B se mantiene a tensión nominal de 154 kV. Suponiendo
X=50 
que una cierta variación de carga en el sistema hace X=10%
disminuir en 5 kV la tensión en la barra M, calcular la
inyección de potencia reactiva necesaria en esta barra X=10%
para restablecer su tensión primitiva. Los valores en % C
están con SB=500 MVA. Figura 2.6
SB=500 MVA; VBL=154 kV  j0,1 (pu) j50 ohm j50 ohm V=Ctte
M
1542
V=Ctte
ZB   47,432 Ohm ; j50 ohm
50050
X   1,05414 (pu) .
L
47,432
j0,1 (pu)
Por otra parte:
5 V=Ctte
V   0,03247 (pu)
154 Figura 2.7
Solución:

La reactancia de Thevenin en M es:

XTH 1
  0,3729 (pu) A partir de la ecuación (2.10):
1 1 1
1,1541  1,1541  1,0541
4 4 4

P Q V Q 1
V   ; con P    y  V   2,6815 se obtiene:
PV Q V V XTH 0,372

V RTH 9

Q  0,03247 *(2,6815)  0,08707 (pu)  Inyectar 43,53 MVAr

2.6. Regulación de tensión mediante transformadores con cambio de TAP (TCT)

El coeficiente (Q/V) de una barra puede en algunos casos alcanzar e incluso superar valores del
orden de 15 [MVAr/kV]. En estas condiciones no resulta adecuado el sistema de inyección de potencia
reactiva debido a la magnitud de las cantidades que habría que poner en juego para compensar las variaciones
de tensión en la barra correspondiente.

Puede suceder que la tensión regulada se desee en la p q

barra q correspondiente al secundario de un transformador
conectado a la barra p como se muestra en la Figura 2.8.

A pesar de la reactancia introducida entre p y q,

(Q/V) puede aún resultar demasiado elevado. En estas
condiciones, lo más adecuado es emplear un TCT que
permita regular la tensión en la barra q “inyectando Figura 2.8.- Transformador con cambio de
una Tap conectado entre las barras p y q de un
tensión adicional”. Esto no genera potencia reactiva, pero sistema
modifica su distribución en el sistema. El empleo de TCT
para regular tensión en una barra determinada, se ilustrará mediante el análisis de los casos siguientes:

2.6.1. Sistema de transmisión radial con un TCT en el extremo transmisor

En el sistema de la Figura 2.9 se tiene: V V2
T
ZL: Impedancia de la línea de transmisión VR Z L
(/fase).
ZT: Impedancia del TCT (/fase), referida al
secundario. ZT P+jQ
V1N, V2N: Tensiones nominales del TCT.
aN=V1N/V2N: Razón nominal del TCT. Figura 2.9.- Sistema radial con un TCT ubicado
en el extremo transmisor
 VT, V2, VR: Tensiones en los puntos que se indican antes del cambio de Tap.
a=V1N/V’2N: razón del transformador después del cambio de Tap (razón no nominal)
V' ; V' ; V' : Tensiones en los puntos que se indican después del cambio de Tap.
T 2 R

La razón de cambio en el Tap t (suponiéndolo en el secundario del transformador), es:

a
t V
 2N (2.18)
a N V2N '

Si se supone que VT=V1N y que ZT, prácticamente no varía al cambiar el TAP, las redes equivalentes
“por fase” antes y después del cambio de TAP son las que se muestran en las Figuras 2.10 a) y b)

aN: 1 ZT
ZL I
+
+
+
+
VT V2N
VR P+jQ
V2
- - -
-
a) ZL
ZT I'
a : 1
+
+ + +
'
VT V2N V2' V 'R P+jQ
- - -
-
b)
Figura 2.10.- Redes equivalentes por fase correspondientes al sistema de la Figura 2.9. a) Antes del cambio de
Tap, b) Después del cambio de Tap
A partir de la Figura 2.10 b) se puede escribir:
'
 (ZT
V2N  ZL ) &I '  (2.19)
&
V&
'
R

Considerando la ecuación (2.2) se puede determinar en forma aproximada, el módulo de la caída de

tensión en las impedancias del transformador y la línea como sigue:

(ZT  ZL ) I '  (R T  R L ) P  (XT  XL ) Q (2.20)

 VR'
VTL

Despejando V ' de (2.18), introduciendo este valor y (2.20) en (2.19) y despejando t, se obtiene:
2
N V2N V'
t R
(2.21)
(R T  R L ) P  (XT  XL ) Q  (V' R) 2

Utilizando esta ecuación se puede determinar el cambio de tap necesario para cumplir determinadas
condiciones. Por ejemplo, si se quiere compensar la caída de tensión en la línea, V '  y por lo tanto:
R
V2

V2N V2 2 (R T 
t
 R L ) P  (XT  XL ) Q  V 2
(2.22)

Si adicionalmente se desprecia la caída en el transformador, se tiene que: R T y XT son iguales a cero y

V2=V2N y por lo tanto:

 2N
V2
t (2.23)
R L P  XL Q  V 2N
2

2.6.2. Sistema de transmisión radial con TCT en ambos extremos de la línea

En el sistema de la Figura 2.11 se tiene: VT TCT "A" V1 V2 TCT "B" VR

ZL
ZL: Impedancia de la línea de transmisión
(/fase)
ZA y ZB: Impedancia de cada TCT (/fase), ZA ZB P+jQ
referida al sector de la línea.
a1N=VAN/VaN: Razón nominal del TCT A Figura 2.11.- Regulación de tensión con TCT en
a2N=VBN/VbN: Razón nominal del TCT B ambos extremos
a1= VAN/V’aN: razón del transformador A después del cambio de Tap
a2= VBN/V’bN: razón del transformador B después del cambio de Tap
VT, V1, V2, VR: Tensiones en los puntos indicados antes del cambio de Tap.
V' ; V' ; V' ; V' : Tensiones en los puntos indicados después del cambio de Tap.
T 1 2 R

La razones de cambio en los Tap t1 y t2 se pueden definir como:

a1 V
t1   aN
' (2.24)
a1N VaN

a 2  VbN
t2 ' (2.25)
a VbN
2N

Las redes equivalentes antes y después del cambio de Tap se muestran en la Figura 2.12.

a : 1 Z
1N A
+ ZL I ZB a :1
+ 2N
+ + +
+
VT VaN
V2 V VR P+jQ
V-1 - - -B
-
- a)
ZL I' ZB
a : 1
a : 1 Z 2
1 A
+
+ + + + +
'
VaN V1' V' '
VB V'
VT P+jQ

- - - - -
-
b)

Figura 2.12.- Redes equivalentes por fase correspondientes al sistema de la Figura 2.11. a) Antes del cambio
de Tap, b) Después del cambio de Tap
 A partir de la Figura 2.12 b) se puede escribir:
'
 (ZA
VaN  ZL  ZB ) &I '  (2.26)
&
V&
'
B

Considerando la ecuación (2.2) se puede determinar en forma aproximada, el módulo de la caída de

tensión en las impedancias de los transformadores y la línea como sigue:

(ZA
 ZL  ZB ) I '  (R A  R L  R B ) P  (XA  XL  XB ) Q (2.27)
 VB'
VALB

Despejando V'
aN de (2.24), introduciendo este valor y (2.27) en (2.26) y despejando t 1, se obtiene:

V V'
t1 aN B (2.28)
(R A  R L  R B ) P  (XA  XL  XB ) Q  (V' B) 2

Por otra parte, de la Figura 2.12.b) y considerando la ecuación (2.25) se tiene que:
'
a2 V
VB VR'  t 2 a 2N R' (2.29)

Introduciendo (2.29) en (2.28), se obtiene finalmente:

VaN t 2 a 2N V'
t1  R (2.30)
(R A  R L  R B ) P  (XA  XL  XB ) Q  (t 2 a 2N V' )R2

Si VT permanece constante, la ecuación (2.30) permite determinar los cambios de TAP necesarios t 1 y
t2 de tal modo que se obtenga un valor deseado V' . En particular, si se desea que VR' sea igual a la tensión
R

nominal del secundario del transformador B, o sea V'

R  , entonces, a partir de (2.30) se obtiene:
VbN
VaN t 2 VBN
t  (2.31)
1
(R  R  R ) P  (X  X  X ) Q  (t V ) 2
A L B A L B BN
2

Si además, se desprecian las caídas internas de los transformadores y se considera t 1 = t2, entonces:

t1  t 2  VaN  R L P  XL Q
(2.32)
VBN V BN
2

Ejemplo 2.4. En el sistema de la Figura 2.13, calcular los cambios de TAP necesarios para mantener 33 kV en
el consumo. Qué potencia reactiva se debería inyectar en la barra de carga para mantener 33 kV, en el caso de
que los transformadores no tuvieran cambiadores de TAP?

Solución:

a) Usando la ecuación (2.30), con t1=t2 se tiene:

 VaN * a 2N * VR'  R T * P  X T * Q V=33 kV
t1  t 2  (12+j40) 
(a 2N * VR' )2
33/230 kV 230/33 kV
en que todos los valores son por fase y por lo tanto: 200 MVA 200 MVA 150 MW
j5,671% b.p. j5,671% b.p. cos =0,8
P=150/3=50 MW/fase
Q=(50/0,8)sin( cos-1 0,8)=37,5 MVAR/fase Figura 2.13
XTR=0,05671ZB2=0,05671(2302/200) 230 230 33
  12  50  70  37,5
XTR=15 /fase 3 33 3
t t   0,90394
RT=12 /fase; 1 2
2
XT=(15+40+15)=70 /fase  230 33 
  
 33 3

b) Usando aproximación de la línea corta:

RPR  XQR 12  50
V  0 R   MVAr / fase
VR R
P
QR 8,57
X
70

Luego: QC=QL-QR=37,5-(-8,57) QC=46,07 MVAR/fase. Es decir; se deberían inyectar 138,21 MVAr

2.7. Regulación de tensión mediante el uso combinado de TCT e inyección de Potencia Reactiva

Los dispositivos de producción o absorción de potencia reactiva se conectan usualmente al enrollado

terciario de transformadores que interconectan redes de transmisión y de distribución. Si el transformador tiene
cambiador de tap, es posible regular independientemente los voltajes primario y secundario, inyectando
potencia reactiva para controlar uno de ellas y modificando los tap para controlar el otro. La Figura 2.14
muestra un esquema típico, en que las variables se suponen “por fase” y las impedancias referidas al primario.

V1 V2 Zp Zs
N
S1 S2 + P+jQ I + +

Zt
V1 VN
V3

V2 CS
CS
- - -
a)
b)
Figura 2.14.- Esquema típico de regulación de tensión combinando TCT e inyección de potencia reactiva. a)
Diagrama unilineal, b) Red equivalente por fase

En general, el problema consiste en determinar los cambios de Tap necesarios en el transformador,

para ciertas condiciones de potencia reactiva del compensador sincrónico y determinada potencia compleja
transferida entre los enrollados primario y secundario, de tal modo que las tensiones V 1 y V2 estén dentro de
ciertos valores especificados. El análisis se efectúa normalmente haciendo las siguientes aproximaciones:

– Se desprecia la potencia activa perdida en el V1

transformador, por lo que Zp, Zs y Zt se
suponen sin resistencia V
VN jXp I
– Las caídas de tensión se calculan por medio 

de expresiones simplificadas
Sea P+jQ, la potencia compleja que llega al I V
nudo N por la rama Zp. El diagrama fasorial de la
Figura 2.15 muestra esta situación, considerando Figura 2.15.- Diagrama fasorial para determinar el
como referencia el fasor de voltaje VN . voltaje en el nudo N

A partir de esta figura se puede escribir:

 V1 Xp Q
 VN  V  (2.33)
VN
Xp P
V 
VN (2.34)

Por otra parte, en la Figura 2.15 se tiene que:

V 2  (VN
1  V) 2  (V) 2 (2.35)

Introduciendo (2.33) y (2.34) en (2.35) y despejando VN se obtiene:

V12  2Xp Q
VN  2
1 V 211(V 2  4Xp Q)  4(X P)
p
2 (2.36)
2 2

En algunos cálculos más simplificados, la tensión V N se determina despreciando V, con lo que resulta
una expresión más sencilla que (2.36).

A partir del valor de VN, se puede determinar el tap en que debe quedar el transformador utilizando,
por ejemplo, la expresión (2.18).

Ejemplo 2.5. El TTTE de la Figura 2.16 tiene conectado en el terciario un compensador síncrono. La potencia
de la carga es de (75+j30) MVA y los voltajes V 1 y V2 deben permanecer en 225 y 70 kV respectivamente. Las
reactancias del TTTE son: X ps=11,5%; Xpt=19,2%; Xst=7,5%, todas en base común 75 MVA. Determinar,
despreciando Xs, el campo de regulación del cambiador de TAP en carga, si se sabe que el compensador
síncrono funciona: a) en vacío; b) Como condensador síncrono (40 MVAR).

Solución:

Parámetros del transformador de 3 enrollados: V1 P T S V2 DATOS TTTE

Sólo interesa el valor de Xp: P/T/S
X 
1
X  X  X  0,116 (pu) ,  como: Y/ /Y
p ps pt st 220/13,2/66 kV
2 2
220 75/45/60 MVA
ZB1   645,33 ohm CS
75
Se tiene que: Xp=74,86 ohm/fase. (XS=-0,001 pu)
Figura 2.16

Usando la ecuación (2.36) y considerando el circuito equivalente de la Figura 2.14 b), se tiene que:

a) Compensador en vacío  P=25 MW/fase; Q=10 MVAR/fase; V1=225/3 kV/fase; V2=70/3 kV/fase.
Reemplazando valores en la ecuación (2.36) se obtiene:

VN1=122,92 kV/fase y VN2=16,4 kV/fase (No válida); por lo tanto; VN=122,92 kV/fase y:

122,92 3 220
 
a'  3,0415 . Además: a   3,3333 . Entonces:
70 N 66

a' 3,0415
ta    0,9125
a N 3,3333
b) Compensador entregando 40 MVAR  Q=3040= 10 MVAr; es decir, Q= 10/3 MVAr/fase.
Reemplazando valores en la ecuación (2.36) se obtiene:

VN1=131,021 kV/fase y VN2=14,41 kV/fase (No válida); por lo tanto; VN=131,021 kV/fase y:

3 220
a" 
131,021  3,2419 . Además: a   3,3333 . Entonces:
70
N 66
a" 3,2419
t b   0,9726 . Es decir, el tap debe variar entre 0,9125 y 0,9726
a 3,3333

Problemas propuestos

2.1. En el sistema de la Figura 2.17, el generador no tiene regulador de voltaje y S B=100 MVA. Antes de
energizar la línea, la tensión en el punto A es del 100%. Determinar:
a. Las tensiones en ambos extremos de la línea después de ser energizada (sin carga).
b. La capacidad del reactor necesario a conectar en “B” para que la tensión en este punto sea del 100%.
c. La tensión en “A” con el reactor conectado.

2.2. En el terciario del Transformador Trifásico de Tres Enrollados (TTTE) de la Figura 2.18, hay un
condensador síncrono que tiene por objeto mantener la tensión de la carga en 100%. Cuando la carga es de
(40+j15) MVA, el condensador síncrono (CS) está en vacío. Determinar la potencia reactiva que entrega el CS
cuando la carga aumenta a (80+j30) MVA. El generador no tiene regulador de voltaje, la potencia base es de
100 MVA y los parámetros del transformador son: X p=13,1%; Xs=-1%; Xt=8%

A B 154 P T S 66
G
kV X=0,05 (pu)
B/2=0,285 (pu)
X=0,05 X=0,067
13,8 kV
GA
Figura 2.17 CS
Figura 2.18

2.3. Una línea corta de transmisión que tiene una impedancia serie de (0,02+j0,1) pu (S B=100 MVA), alimenta
una carga de 200 MVA, Factor de Potencia 0,8 inductivo. En el extremo receptor existe un compensador que
permite que con 220 kV en el extremo transmisor, el voltaje en el extremo receptor sea de 225 kV.
Considerando que el desfase entre ambos voltajes es pequeño y que el voltaje en el extremo transmisor se
mantiene constante, Determinar:
a. La capacidad (MVAr) del compensador síncrono que está conectado en el extremo receptor.
b. El voltaje que aparece en el extremo receptor al desconectar la carga.

2.4. El esquema de la Figura 2.19, muestra un diagrama simplificado de la Central Termoeléctrica Bocamina.
Los coeficientes (P/V) y (Q/V) en la barra de 66 kV son respectivamente: –130 [MW/kV] y 56,5
[MVAr/kV]. La tensión en dicha barra es de 67,3 kV, cuando está la carga de 4 MVA conectada, a Factor de
Potencia 0,9 (inductivo). En estas condiciones se pretende hacer partir un motor de inducción tipo Jaula de
Ardilla, que mueve las bombas de la caldera, sabiendo que la tensión mínima de partida del motor es del 90%
y que durante la partida consume 5 veces la corriente nominal con cos =0,3. Determinar si es posible hacer
partir el motor en estas condiciones.

2.5. Los sistemas de la Figura 2.20, están unidos a través de un autotransformador regulador cuyo rango es
VN15 x 1,25%. Cuando está en su derivación nominal y las tensiones son V1=1,1 (pu); V2=0,9 (pu), el
Sistema 1 entrega (100 + j80) MVA. Determinar la derivación en que se debe ajustar el regulador para que la
potencia reactiva que entregue el Sistema 1 no sea superior a 30 MVAR (V 1 y V2 constantes). Considerar
SB=100 MVA.

4 MVA
66 kV cos =0,9 V1 VN _+ 15 x 1,25% V
2

66/4,16 kV t:1
11 MVA MI
X=35%
X=8% 1,75 MVA
Sistema 1 Sistema 2
cos =0,8
Figura 2.19 Figura 2.20

2.6. En el sistema de la Figura 2.21, los generadores 1 y 2 tienen reguladores de tensión que mantienen
constante la tensión en bornes. El transformador de 154/69 kV tiene cambiador de TAP bajo carga en el lado
de 69 kV con 19 derivaciones en pasos de 1,25% y con el valor nominal en la derivación central. El cambiador
de TAP es comandado por un sistema de control que trata de mantener constante la tensión en barras de 69 kV
y que reacciona cuando la tensión varía en más de 1,5% del valor de ajuste. La tensión en barra de 69 kV es la
nominal y el cambiador se encuentra en la posición (69+4x1,25%). En ese instante se conecta un banco de
condensadores de 6 MVAr a la barra de 69 kV. ¿Cómo reacciona el sistema de regulación? ¿En qué valor
queda la tensión en la barra de 69 kV después que ha reaccionado el sistema de regulación? Todos los valores
en % están en base propia.

(12,5+j37,5) G2
G1 100 MVA (8,33+j30)

50 MVA
50 MVA
Xs=120%
200 MVA X=10%
100 MVA X=12%
Xs=80%
154/69 kV
X=16% 120 MW
F. de P. = 0,96

Figura 2.21

2.7. En el sistema de la Figura 2.22, la potencia en el consumo puede variar del siguiente modo: (0  j3)  S&
C  (30  j20) MVA . Considerando que la tensión en el extremo transmisor VT se mantiene constante
en 12 kV y que la del extremo receptor V R no debe variar en más de  5%; determinar si es necesario inyectar
o absorber potencia reactiva en el extremo receptor y a partir de que valor o valores de la potencia del
consumo. Valores en porcentaje en base 100 MVA.

G1 VT VR
j40%

80 MVA 12/66 kV 66/33 kV

X s =30% X=10% X=12% Sc
Figura 2.22
 CAPITULO 3

CALCULO DE FLUJOS DE POTENCIA

3.1. Introducción

El cálculo y análisis del flujo de potencias en la red de un Sistema Eléctrico de Potencia (SEP) es uno
de los aspectos más importantes de su comportamiento en régimen permanente. Consiste en determinar los
flujos de potencia activa y reactiva en cada línea del sistema y las tensiones en cada una de las barras, para
ciertas condiciones preestablecidas de operación.

Hasta el año 1950, el Cálculo del Flujo de Potencias (CFP) se realizaba utilizando principalmente los
Analizadores de Redes de Corriente Alterna (ARCA) y en algunos casos, los Analizadores de Redes de
Corriente Contínua (ARCC) que corresponden a una simulación a escala del Sistema Real. En la actualidad,
el CFP se realiza fundamentalmente, utilizando los computadores digitales por las grandes ventajas que éstos
presentan respecto a los analizadores de redes.

El análisis del flujo de potencias (AFP) permite:

 Programar las ampliaciones necesarias del SEP y determinar su mejor modo de operación, teniendo
en cuenta posibles nuevos consumos, nuevas líneas o nuevas centrales generadoras.
 Estudiar los efectos sobre la distribución de potencias, cuando se producen pérdidas temporales de
generación o circuitos de transmisión.
 Ayudar a determinar los programas de despacho de carga para obtener un funcionamiento óptimo.

3.2. Planteamiento del problema básico

Considérese el SEP elemental de dos barras de la Figura 3.1 y su circuito equivalente por fase que se
muestra en la Figura 3.2. La línea L 12 se ha representado por su circuito  nominal y donde:

S& 1 y S& 2 : Potencias complejas netas de las barra 1 y 2 respectivamente, representadas como fuentes
de potencia activa y reactiva, que corresponden a la Potencia Generada menos la Potencia
Consumida.

S& 12 y S& 21 : Flujos de potencia compleja que van desde la barra 1 a la barra 2 y viceversa.
1 2
1 2
G1 G2 V1 RL+jXL
V2
L12

S12 S21
S G1 SG2 S1 S2
Y/2 Y/2
SC1 SC2

Figura 3.1.- Sistema elemental de dos barras para Figura 3.2.- Circuito equivalente por fase del sistema
plantear el problema básico de la Figura 3.1

En la Figura 3.1, las potencias complejas netas de las barras 1 y 2 son:

S& 1  S& G1  S& C1  (PG1  PC1 )  j (Q G1  Q C1 )  P1  jQ1

(3.1)
S& 2  S& G 2  S& C2  (PG2  PC2 )  j (Q G 2  Q C2 )  P2  jQ 2
 26
En el circuito de la Figura 3.2 se puede escribir:
S* P  jQ V&  V&
Y
1
 1
 V& 1  1
1 2
V1* V1* 2 R L  jX L
(3.2)
S* P  jQ Y V&  V&
2
 2 2
 V& 2  2 1
V2* V2* 2 R L  jX L

Estas ecuaciones, que relacionan las tensiones con las potencias activas y reactivas, presentan las
siguientes características.

 Son algebraicas y no lineales.

 La frecuencia no aparece en forma explícita porque se la supone constante.
 El sistema de cuatro ecuaciones, tiene 12 variables en total: PG1, PG2, QG1, QG2, PC1, PC2, QC1, QC2, V1,
1, V2, 2, por lo que no es posible obtener una solución para ninguna de ellas a menos que se reduzca
el número de incógnitas, fijando de antemano algunas variables.

En relación a esto último, una forma posible de resolver el problema es la siguiente:

 A partir de los datos del consumo suponer conocidas e independientes del voltaje, las potencias de las
cargas PCi, QCi, con i = 1,2.
 Fijar a priori dos variables de generación P G2 y QG2 por ejemplo. No se pueden fijar las cuatro
variables de generación debido a que las pérdidas en el sistema no son conocidas inicialmente.
 Fijar el módulo y ángulo de la tensión en barra 1; es decir; suponer conocidos V 1, 1. En particular,
puede tomarse esta tensión como referencia, o sea, 1=0

En estas condiciones, el sistema de 4 ecuaciones (3.2) queda con sólo 4 variables: P G1, QG1, V2, 2.

3.3. Modelo de representación del SEP

Teniendo presente el análisis del problema básico y con el objeto de establecer un procedimiento
general para el CFP, es necesario considerar lo siguiente:

3.3.1. Tipos de Barras

Asociados a cada barra p de un SEP existen cuatro variables, P p; Qp; Vp; p. Según las variables
conocidas y desconocidas, las barras se clasifican en los siguientes grupos:

 Barras de Carga (Barras P-Q): Pp y Qp están especificadas; Vp y p son las incógnitas

 Barras de tensión controlada (Barra P-V): Pp y Vp están especificadas; Qp y p son las incógnitas.
En este tipo de barra debe existir alguna fuente controlable de potencia reactiva.
 Barra flotante (Barra V& ): Vp y p están especificados; Pp y Qp constituyen las incógnitas. En
esta barra debe existir por lo menos un generador. La necesidad de definir esta barra nace del hecho
que no es posible especificar a priori, la potencia total que es necesario generar en el sistema debido a
que inicialmente no se conocen las pérdidas en el mismo. La barra flotante debe suministrar la
diferencia entre la potencia compleja inyectada al sistema en el resto de las barras y la carga total más
las pérdidas. Esta barra se conoce también con otros nombres tales como: de referencia, oscilante, de
relajación (slack).
3.3.2. Representación de los elementos del SEP

a. Líneas: Se representan usualmente por su circuito  nominal. Para una línea conectada entre las
barras p y q de un SEP, el circuito equivalente corresponde al mostrado en la Figura 3.3. En algunos casos,
basta representar la línea por su impedancia serie.

b. Transformadores: Cuando funcionan en su razón nominal, se representan por su impedancia de

cortocircuito. Cuando operan con cambio de TAP y razón no nominal, se pueden representar por su circuito
equivalente  que se muestra en la Figura 3.4, cuyos parametros se indican en la ecuación (3.3).

Spq Zpq ( Ypq ) I I2

Sqp 1
p q + B +
I pq Iqp
V C V2
Vp Yp' q /2 Yp' q /2 Vq 1 A

- -

Figura 3.3.- Circuito equivalente  de una línea para Figura 3.4.- Modelación circuital en tanto por
el cálculo de flujos de potencia unidad de un transformador con cambio de TAP

Y1 1 Y Y1 1
A    B C   
      (3.3)
   

Con =1+t1 y =1+t2; y donde t1 y t2, representan el cambio del Tap, en el lado respectivo.

c. Generadores: Se consideran normalmente como fuentes de potencia activa y reactiva.

3.4. Planteamiento matemático del problema para un SEP de “n” barras

3.4.1. Ecuaciones de Barras

Considérese una barra p cualquiera de un sistema p

tal como se muestra en la Figura 3.5. La potencia Resto del Ip
compleja neta, S& y la corriente inyectada en
p
SEP +
Vp
la barra p, &I p están relacionadas por las -
siguientes
ecuaciones, que constituyen las ecuaciones de Figura 3.5.- Representación de una Barra p en un SEP
barras.
S& p  V I *   jQ
P
pp p p

&I  S*p  jQp (3.4)

p 
Vp*
Pp
Vp*

3.4.2. Ecuaciones del flujo de potencias

A partir de la Figura 3.3 se puede escribir:

&I  (V& p  V& q )Y& pq  V& p

pq pq (3.5)
(Y&
'
2)

La potencia compleja que fluye desde la barra p a la q está dada por:

pq pq pq p p q pq
S& pq  V& 
 (Y' 
2) * V 2  V*Y* (3.6)
Y
I* p
*
V&
 Análogamente, la potencia compleja que fluye desde la barra q a la barra p estará dada por:

S& qp  V& 
 (Y' 
2) * V 2  V * V& (3.7)
I*
Y
q
*
qp qp qp Y* q p q qp

Las expresiones (3.6) y (3.7) corresponden a las ecuaciones del flujo de potencia a través de la línea.
Conviene indicar que Y& pq  y
pq
2)  qp 2) . Además, Y& es el inverso de la impedancia entre las
(Y& (Y&
' '
Y& qp
pq

barras p y q, el que no debe confundirse con el valor correspondiente en la matriz Y B, de la ecuación (3.10).

3.4.3. Potencia perdida en la transmisión

De acuerdo con los sentidos adoptados para S& pq y S& qp , la potencia compleja perdida en la línea será:

S& Lpq  S& pq 

(3.8)
S& qp

3.4.4. Cálculo de las tensiones de barras

Las ecuaciones (3.6) y (3.7) indican claramente que para resolver el problema del flujo de potencias
se requiere determinar previamente las tensiones en todas las barras que correspondan. Empleando el método
nodal de resolución de circuitos, en forma matricial, para la red de un SEP de n barras se puede escribir:

IB   YB VB (3.9)


donde IB  esa el vector de corrientes inyectadas a las barras; YB es la matriz admitancia de barras y
VB  es el vector tensiones de barra, definidos como:
 &I1   V& 1   Y& 11 Y& Y& L Y& 1n 
···
12 1p
     Y& 
&
 2  2  21 L Y& L 2n 
I 22 Y&
 M V  YM 2p
M
&M  & L M L
Y
M&
IB   & VB    Y B   Y& L Y& L


(3.10)
  pI   V&  p1 pn
p2 pp Y&
 M  M M
M L M L
 p 
M
&  &  & & & & 
I n  Vn  Yn1 Yn L L Ynn 
2 Ynp

Teniendo presente que según (3.4), las corrientes inyectadas en las barras dependen de las potencias
complejas netas respectivas y considerando (3.9) y (3.10), se puede escribir:
*
S1
 Y& 11V& 1  Y& 12 V& 2  ···  Y& 1p V& p  ···  Y& 1n V& n
V1*
*
S
 2
 Y& 21V& 1  Y& 22 V& 2  ···  Y& 2p V& p  ···  Y& 2n V& n
V2*
M M M M M
* (3.11)
Sp
 p1 V  Y p2  ···  pp  ···  &
pn V
Y& & V& 2 V& p n
Vp* 1
& Y& Y& M
M M
S* M M
n
 Y& n1V& 1  Y& n 2 V& 2  ···  Y& np V& p  ···  Y& nn V& n
Vn*
 Este sistema de ecuaciones es similar al obtenido en el problema elemental de 2 barras; es decir; las
ecuaciones son algebraicas y no lineales, por lo tanto es necesario resolverlo mediante técnicas de
aproximaciones sucesivas.

3.5. Técnicas de solución para el problema del flujo de potencias

Existen actualmente diversos métodos para resolver el problema de cálculo del flujo de potencias, los
que reciben nombres según sea el procedimiento que se aplique para calcular las tensiones. Entre ellos
podemos mencionar el de Gauss, el de Gauss-Seidel, los de Newton-Raphson (Completo, Desacoplado,
Desacoplado rápido), el flujo DC, etc.

Estudiaremos a continuación, cada uno de ellos, considerando en primer lugar el procedimiento

general y luego las aplicaciones al cálculo de flujo de potencias

3.5.1. Método de Gauss

Se emplea para resolver un problema lineal o no lineal. Por simplicidad se considerará un sistema
lineal de ecuaciones, como el indicado en (3.12), para fundamentarlo. Sin embargo, su aplicación a un
sistema no lineal resulta inmediata.

a11x1  a12 x 2  a13 x 3  y1

a 21x1  a 22 x 2  a 23 x 3  y 2 (3.12)
a 31x1  a 32 x 2  a 33 x 3  y 3

Despejando x1 de la primera ecuación, x2 de la segunda y x3 de la tercera se obtiene:

x1  (y1  a12 x 2  a13 x 3)/a11

x 2  (y 2  a 21x1  a 23 x 3)/a 22 (3.13)
x 3  (y3  a 31x1  a 32 x 2 )/a 33

Sean
x 0 , x 0 , x 0 , valores iniciales estimados a priori de la solución del sistema (3.12), entonces,
1 2 3
reemplazando estos valores en (3.13) se tiene:

x1  (y1
1  a12 x 0  a13 x 30 ) / a11
2
x  (y2  a x 0  a 23x 03 ) / a
1
(3.14)
2
21 1 22
1 2
x1  (y3  a x 0  a 32 x 0 ) / a
3 33
31

El procedimiento continua hasta que se satisface algún “criterio de convergencia” tal como, por
ejemplo, el indicado en (3.15), donde  es una cantidad de valor pequeño y positivo. A cada etapa del proceso
se le denomina “iteración”.
x k1  x k  
con i  1, 2, (3.15)
3
i i

Aplicando el método a un sistema de n ecuaciones con n incógnitas; para la incógnita xi, después de k
iteraciones, y con i= 1, 2,.. . .; n; se puede escribir:

 
1 k
x ik1  a   y i  n a ijx (3.16)
j

 j1
ii
 ji 
 Los inconvenientes de este procedimiento son el gran número de pasos que se requiere para llegar al
resultado y la ocurrencia relativamente alta de situaciones en que no hay convergencia, por lo que no se
utiliza para resolver el problema de cálculo de los voltajes de la ecuación (3.11). Sin embargo, constituye la
base para la formulación del Método de Gauss-Seidel, lo que justifica su análisis. Al aplicar la ecuación
(3.16) al problema de cálculo de los voltajes en las barras del sistema de ecuaciones (3.11) se obtiene:

1 P  jQ &Y 
V & k1
  p p
 n
V& k (3.17)
  pq q 
Y& pp
p (Vkp )*
 q1
 
qp 
Donde: p = 1, 2, 3, ….., n; q= 1, 2, 3,......, n y p  s (barra slack).

La ecuación (3.17) se conoce como método de Gauss Y B, porque usa el Método de Gauss y se trabaja
con la matriz admitancia de barras del sistema eléctrico. La expresión es válida sólo para las barras de carga.
En el caso en que el SEP contenga barras de tensión controlada, la ecuación (3.17) debe ser modificada, pues
en este tipo de barras no se conoce el valor de la potencia reactiva Q p. Por lo dicho en el párrafo anterior, la
modificaciones requeridas se estudiarán al considerar el Método de Gauss-Seidel Y B.

3.5.2 Método de Gauss-Seidel

a. Caso general: Corresponde a una modificación del método de Gauss tendiente a acelerar la
convergencia del proceso iterativo. En el método de Gauss se calculan todos los valores de las incógnitas
correspondientes a una iteración y luego se emplean para determinar los nuevos valores de las incógnitas en
la iteración siguiente. En el método de Gauss-Seidel en cambio, los valores calculados en una iteración
determinada, se utilizan inmediatamente para calcular los valores de las incógnitas que restan por calcular en
la misma iteración.

k1 De
k1este modo, si el proceso de cálculo se encuentra en la iteración (k+1) y ya se han determinado
xk1 , x , ....., x entonces, los valores que se utilizan para calcular serán
; xk1
1 2 i1 i
k k
xk k1
,x k1
,...., x k1
, x , x ,....., x
1 2 i1 i2 n
i1

Por tanto, la fórmula iterativa del Método de Gauss-Seidel aplicada a un sistema de n ecuaciones de
la forma dada por (3.12) es:

1
x k1  a y i1 n 
i  i   x k1   a ijx j  (3.18)
a ij k
j
ii  j1 ji1 

b. Aplicación del método de Gauss-Seidel YB al cálculo flujos de potencia: El cálculo de las

tensiones de barras aplicando el procedimiento explicado anteriormente es distinto según sean los tipos de
barras existentes en el SEP. Por ello se considerarán en primer lugar los sistemas con barras de carga y
flotante solamente, por ser el caso más simple. A continuación se analizará la situación de las barras de
tensión controlada.

b1. Sistemas con barras de carga y flotante solamente: Aplicando la ecuación (3.18) al sistema (3.11)
se tiene:

V& k1  1 Pp  jQp p1 & &

k1 n & & k
p
Y & k *
  Y pq   Ypq Vq  (3.19)
Vq
pp  (Vp ) q1 qp1 
Donde: p = 1, 2, 3, ….., n; q= 1, 2, 3,......, n y p  s (barra slack).
 La secuencia de solución según este método es como sigue:

1. Se suponen valores iniciales de tensión para todas las barras a excepción de la flotante, cuya tensión
está especificada, o sea es dato del problema, al igual que P p y Qp en todas las barras de carga y los
términos de la matriz admitancia de barras (YB)

2. Se aplica la fórmula iterativa (3.19) hasta que se cumpla algún criterio de convergencia, por ejemplo:

V k1  Vk  
especificado con p  1, 2, 3,....., n
p p 1
(3.20)
k1  k   especificado con p  1, 2, 3,....., n
p p 2

3. Determinadas las tensiones V& p , se calculan los flujos de potencia S& pq y S& qp aplicando (3.6) y (3.7).

4. Conocidos los valores de S& pq y S& qp se determinan las pérdidas en el sistema, empleando (3.8).

b2. Sistemas con barras de carga, tensión controlada y flotante: Normalmente un SEP incluye
además de las barras de carga y flotante, barras de tensión controlada (BTC) que tienen por objeto permitir
regular la tensión en uno o varios puntos del sistema. En las barras de tensión controlada debe existir una
fuente regulable de potencia reactiva para poder cumplir su cometido.

Debido a que en este tipo de barra sólo se conocen el módulo de la tensión y la potencia activa, es
necesario calcular previamente la potencia reactiva, antes de emplear la ecuación (3.19) para determinar el
voltaje complejo en ella.

A partir de la ecuación para la barra p de la expresión (3.11), se puede escribir:

S*  P
& & & V&  Y&
(YV&
 jQ  Y  Y  ···  Y  ···  V&V *
*
V V V n (3.21)
p p p
p p1 1 p2 2 pp p pn n p  pq q
q1

es decir:

Q *
n Y& V& 
(3.22)
pq q 
p Im ag Vp

 q1 
Cuando se emplea la ecuación (3.19) en una BTC, el valor de Q p, que debe emplearse corresponde al
indicado por (3.22), el que se debe actualizar en cada iteración. Al determinar el voltaje, debe tenerse en
cuenta que su módulo en esta barra está especificado y por lo tanto sólo puede cambiar su ángulo.

Límites de Potencia reactiva en una Barra de tensión Controlada: En el cálculo del flujo de
potencias en un SEP con Barras de tensión controlada es necesario tomar en cuenta los límites de potencia
reactiva de las fuentes de potencia.

Sea p una BTC, entonces el valor de Qp se puede escribir como:

Qp  QGp  QCp
(3.23)

Además:

(QGp)máx : Valor máximo de generación de potencia reactiva de la fuente.

(Qgp)mín : Valor mínimo de generación de potencia reactiva de la fuente.
 QCp : Potencia reactiva de la carga en la barra
Los límites de potencia reactiva para la barra p serán:

(Qp ) mín  Qp  (Qp ) máx (3.24)

donde: (Qp ) mín  (QGp ) mín  QCp y (Qp ) máx  (QGp ) máx  QCp

Si el valor de la potencia reactiva calculado según (3.22) en una iteración cualquiera k, Q pk, excede el
límite máximo o mínimo prefijado, significa que es imposible obtener una solución con la tensión
especificada en esta barra y en consecuencia, ella debe ser considerada como una barra de carga en esa
iteración, en la cual la potencia reactiva es igual al límite superior e inferior según corresponda. En las
iteraciones siguientes, el método intentará mantener el voltaje especificado originalmente en esa barra,
siempre que no se violen los límites de Q p. Esto es posible, porque pueden ocurrir cambios en otros puntos
del sistema, que lo permitan.

Para explicar mejor el procedimiento, 1 2

considérese el sistema de 3 barras de la
Figura 3.6.

Sean:
Barra 1: Flotante
Barra 2: de carga 3
Barra 3: de tensión controlada

La secuencia de cálculo aplicando el

Método de Gauss-Seidel YB es: Figura 3.6.- Sistema de tres barras para explicar el método de
Gauss-Seidel YB

1. Especificar los datos necesarios: V1, 1, P2, Q2, P3, V3 y los parámetros para determinar la matriz YB
2. Suponer los valores iniciales: V 0 , 0 , 0 . Normalmente se usa 1,0 (pu) para los módulos de voltaje y
2 2 3
0º para los ángulos; V3 está especificado

3. Calcular la tensión en la barra 2

V& 
1
 0
1  P2  jQ 2 &
Y& 
V V& (3.25)

Y& 
 21 1 23 3
Y& 22 (V
2 0
2 ) 
 *

4. Calcular la potencia reactiva en la barra 3


Q 0   Im (V 0 ) * (Y& V& 
V&&
1
 V&
0
(3.26)
Y& Y 
3 3 31 1 32 2 33 3

5. Verificar si Q 0
3 está dentro de los límites establecidos

6. Si Q 0
3 está dentro de los límites, determinar la tensión en la barra 3 según (3.27), mantener el valor
especificado para V3 y cambiar el ángulo inicial por el determinado con (3.27)
1 1 P3 & jQ03 & & &1

V& 3  &  V 2   Y 31V 1  Y
32 (3.27)
Y 33  (V30 )* 
7. Si Q 0 0
3 no está dentro de los límites, reemplazar Q 3 en (3.27) por el valor del límite excedido,
tomando el valor de V& calculado en (3.27) completo (módulo y ángulo).
1
3
 8. Verificar que se cumpla el criterio de convergencia tal como se indica en ecuación (3.20), por
ejemplo, en todas las barras. Si se cumple, el proceso de cálculo de las tensiones finaliza y se
determinan los flujos de potencia y las pérdidas en las líneas según ecuaciones (3.6) y (3.7) y (3.8).

9. Si el criterio de convergencia no se cumple, volver al punto 3 y repetir el proceso.

Ejemplo 3.1. Para el sistema de tres barras de la Figura 3.7, los datos en pu, base común, se dan en las Tablas
Nº 1 y Nº 2. Realizar una iteración con el método de Gauss-Seidel Y B, para determinar el voltaje en todas las
barras. Con los valores obtenidos, determinar los flujos de potencia en todas las líneas, las pérdidas del
sistema, la potencia entregada por el generador de la barra slack y la verificación de la potencia entregada a la
carga SC1.

G3 G2 PG2 Tabla Nº 1: Datos de las líneas

3 2 Línea Z (pu) Y’/2 (pu)
1-2 0,04+j0,12 j0,05
1-3 0,02+j0,06 j0,06
2-3 0,06+j0,18 j0,05

Tabla Nº 2: Datos de las barras

1 Barra Nº Tipo V (pu) PG QG PC QC
SC1 1 PQ - - - 0,6 0,25
Figura 3.7 2 PV 1,04 0,2 - 0,0 0,0
3 SL 1,06 - - 0,0 0,0
Límites de generación de Q en la barra 2: 1  QG2  1

Solución:

a) Determinación de la matriz de admitancia de barras YB

Y11 1 1
   j0,05  j0,06  7,5  j22,39  23,6128  71,48º
0,04  j0,12 0,02  j0,06
1 1
Y22    j0,05  j0,05  4,1667  j12,4  13,0813  71,43º
0,04  j0,12 0,06  j0,18
1
  1
Y3  j0,06  j0,05  6,6667  j19,89  20,9775  71,47º
3
0,02  j0,06 0,06  j0,18

Y12 1
 Y21    2,5  j7,5  7,9057108,43º
0,04  j0,12
1
Y2    1,6667  j5  5,2705108,43º
Y32 0,06  j0,18
3 1
  5  j15  15,8114108,43º
 Y13 0,02  j0,06
Y3

b) Valores iniciales y otros datos

V&  10º ; V&  1,040º ; V&  1,060º , especificado (Barra slack)

0 0
1 2 3

S& 1  0  0,6  j(0  0,25)  0,6  j0,25; P2  0,2  0  0,2

Límites de Q en la barra 2 :
 1  Q 2 1
c) Proceso iterativo: Utilizando la ecuación (3.19)

V &1 1   0,6  j0,25 

  7,9057108,43º *1,040º15,8114108,43º *1,060º
 
23,6128  71,48º 
1
10º
V1&
1
 1,0401  1,23º

Antes de determinar V2&

1
se debe calcular la potencia reactiva neta de la barra 2, expresión (3.22)

Q 02   Im7,9057108,43º *1,0401  1,23º13,0813  71,43º *1,040º5,2705108,43º*1,060º*1,040º

Q0  0,2693 (está entre los límites)
2

La tensión V2& se determina usando de nuevo la expresión (3.19), es decir:

1

V &1 1 0,2  j0,2693 

  7,9057108,43º *1,0401  1,23º5,2705108,43º *1,060º
2  
13,0813  71,43º  1,040º
V2&  1,04140,23º V2&  1,040,23º
1 1

d) Cálculo de los flujos de potencia en las líneas: Usando las expresiones (3.6) y (3.7) y considerando que:

Y12
 Y21  1  2,5  j7,5  7,9057  71,57º
0,04  j0,12
Y2 1
  0,06  j0,18  1,6667  j5  5,2705  71,57º
Y32
3
1
 0,02  j0,06  5  j15  15,8114  71,57º
 Y13
Y3

se tiene:

S& 12  7,905771,57º0,05  90º *1,04012  1,0401  1,23º *1,04  0,23º *

7,905771,57º S&  0,2056  j0,0182  0,2064174,94º
12

S& 21  7,905771,57º0,05  90º*1,04 2  1,040,23º *1,04011,23º *

7,905771,57º S&  0,2073  j0,1211  0,2401  30,29º
21

S& 23  5,270571,57º0,05  90º*1,04 2  1,040,23º *1,060º *

5,270571,57º S&  0,0125  0,1654  0,1659  94,33º
23

S& 32  5,270571,57º0,05  90º*1,06 2  1,060º *1,04  0,23º *

5,270571,57º S&  0,0132  j0,0572  0,0587577,0º
32

S& 31  15,811471,57º0,06  90º *1,06 2  1,060º *1,04011,23º

*15,811471,57º S&  0,4617  j0,1345  0,480916,24º
31

S& 13  15,811471,57º0,06  90º*1,04012  1,0401  1,23º *1,060º

*15,811471,57º S& 13  0,4572  0,2533  0,5227  151,01º
e) Pérdidas: Utilizando la ecuación (3.8) se tiene:

S&  S& 12  S& 21  0,2056  j0,0182  0,2073  j0,1211  0,0017 

L12

j0,1029 S& L23  S& 23  S& 32  0,0125  j0,1654  0,0132  j0,0572 

0,0007  j0,1082 S&  S& L31 S&  0,4617  j0,1345  0,4572 
31 13
j0,2533  0,0045  j0,1188

f) Potencia entregada por el generador de la barra slack

S& G3  S& 31  S& 32  0,4617  j0,1345  0,0132  j0,0572  0,4749  j0,1917

g) Verificación de la potencia recibida por la carga SC1

S& C1  S& 12  S& 13  (0,2056  j0,0182)  (0,4572  j0,2533)  0,6628  j0,2351

Observación: No corresponde exactamente al valor especificado para la carga. ¿Porqué?

Factores de Aceleración: La experiencia con el método de Gauss-Seidel Y B para el cálculo de flujos

de potencia ha mostrado que se puede reducir considerablemente el número de iteraciones requeridas si la
corrección en el voltaje de cada barra se multiplica por alguna constante que la incremente, para que el
voltaje sea más cercano al valor al que se está aproximando. El multiplicador que realiza esto, se denomina
factor de
p
aceleración. La diferencia entre el valor de voltaje de la barra p calculado en la iteración actual V& y el
k1

mejor que se obtuvo en la iteración anterior (V& ) se multiplica por un  apropiado para obtener una mejor
k
p
ac
corrección que se agrega a este último. Es decir:

(V&
k1p
) ac  (V&
k
  V& k1
 (V& )
k
(3.28)
)
p ac  p p ac

La elección del valor de  depende del sistema en estudio. Hay valores óptimos de los factores de
aceleración para cualquier sistema. Una mala selección de ellos puede dar como resultado una convergencia
menos rápida o hacerla imposible. Normalmente, en los estudios de flujos de potencia,  varía entre 1,3 y 1,8.
Generalmente, un factor de aceleración de 1,6 para las componentes real e imaginaria es una buena selección.
Sin embargo, es posible que el factor de aceleración utilizado para la componente real pueda diferir del usado
para la componente imaginaria.

3.5.3. Método de Newton Raphson

a. Formulación general: Este método es mas sofisticado que los anteriores y exige un mayor volumen
de cálculos, pero asegura convergencia en un mayor número de veces y además en forma más rápida. El
problema matemático a resolver consiste en n relaciones no lineales del tipo f(x i)=0. Es decir, se trata de un
sistema de n ecuaciones de la forma:

f1 (x1 , x 2 ,..., x n )  0
f 2 (x1 , x 2 ,.. , x n )  0
0
...............................
f n (x1 , x 2 ,.. ., x n ) 
(3.29)

Si se supone una estimación inicial del vector solución:

x  x 0......x 0
0
1 2
n (3.30)
x 0

    
Al que le falta un residuo x 0  x 0 x 0 ...... x 0 , para llegar a la solución correcta; esto es,
1 2 n
tener f (x 0  x 0 )  0 , aunque f (x 0 )  0 , se tiene:
i i i

f1 (x10  x10 , x2 0  x2 0 ,....., nx 0  x

n )0
0

f 2 (x  x , x  x ,....., x  x 0 )  0
0 0 0 0 0
1 1 2 2 n n (3.31)
...............................................................
fn (x 10  x 10 , x 20  x 20 ,....., xn 0  xn0 )  0

Desarrollando cada ecuación en serie de Taylor en torno a los valores x 0

i se tiene:
 f 0  f 0
f1 (x 0  x 0 ,....., x 0  x 0 ) (x 0 ,....., x 0 )  x 0  1   .....  x 0  1   
f
1 1 n n 1 1
n 1
 x1 
n
 x n 1
............................................................... (3.32)
0 0
 f   f 
f n (x 0  x 0 ,....., x 0  x 0 )  (x 0 ,...., x 0 )  x 0  n   .....  x 0  n   
f
1 1 n n n 1
n 1
 x1 
n
 x n  n

donde i es el residuo en la serie de Taylor, que contiene los términos de orden superior
0
 f 
 xj i: representa las correspondientes derivadas parciales, evaluadas en x
0

 i

Como los x 0i son pequeños, se pueden despreciar los términos de orden superior y se obtiene:

 f 0  f 0
f1 (x 0 ,....., x 0 )  x 0  1   .....  x 0  1   0
1 n 1 n
 xn 
 x1 
............................................................... (3.33)
0 0
 f   f 
fn (x ,....., x )  x  n   .....  x  n   0
0 0 0 0

1 n 1 n
 x n 
 x1 
con i, j = 1, 2, …, n

Matricialmente se puede escribir:

 f 0
 f 0
f (x 0 )  1    0 
  L 1 x 0
1 i  x1  x n    
  
  M  M 1  M
 M   L M     (3.34)
 f (x 0 )  
n 0 0
 fn   x 0
 f n 
0

i L
   n
   
 x1   x n  

Es decir: 
 0
f (x ) J  x
0 0 0
 (3.35)

Donde cada vector y matriz está definido según las ecuaciones (3.34) y (3.35), o sea:
 f1 (x i0 ) 
 0 

f (x )   M 
f (x 0 )
0
Vector función evaluada en x i
n
 i 

 f 0
 f 0

 1

 x1  L  xn 1  

 
J0  M
 0
L
M
0
 Matriz Jacobiana evaluada en x 0i

 fn   fn  
L  
 x
 n 
 1

  0
x
   1
x   M 
0
Vector residuo evaluado en x 0i
x 0 
 n

A partir de (3.35), el vector residuo evaluado en x i0 ; x 0   se puede

escribir:
    
x   J
0 0
1
0
f (x )  (3.36)

En general entonces, el residuo en una iteración k es:

x   J  f (x )
k k
1
k
(3.37)
  (vector de valores aproximados de la variable), entonces puede
Suponiendo que se conoce x k

obtenerse una aproximación mejor x  de la k1

forma:
(3.38)
x   x  x   x   J  f (x )
k 1 k k k k
1
k

 
Como se han despreciado los términos de orden superior, x k1 no será la solución correcta y se
debe repetir el proceso en forma iterativa, hasta que se satisfaga algún criterio de convergencia, tal como:

x k1  xk  
(3.39)

b. Aplicación al cálculo de flujos de potencia: En el caso de un Sistema de Potencia, los x i

corresponden a las tensiones de las barras (módulo y ángulo), de manera que la ecuación (3.37) se puede
escribir como:

  k  k 1  P 
 
k
 k 
J  k  (3.40)
V Q 
   

en que:

P   esp calc 

    Pesp  P calc  (3.41)
Q
  Q  Q 
 Donde Pesp y Qesp son los valores de P y Q especificados o programados y P calc y Qcalc son los
respectivos valores que se van calculando en cada iteración. Para mayor comodidad, a los valores calculados,
se les eliminará el superíndice calc, es decir, se designarán simplemente como P y Q.

Según (3.21), los valores de P y Q en la barra p se pueden obtener a partir de:

S&  P 
jQ 
V& n * *
Y V (3.42)
p  pq q
p p p
q1

Expresando los voltajes de barras en forma polar y las admitancias de línea en forma rectangular se
tiene que:

V& p  V& q  Vq y Y& pq  G pq  j (3.43)

Vp  p ;  q B pq

Reemplazando (3.43) en (3.42), con pq  p  q se obtiene:

&
P 
S
 jQ   j sin )
n V V (G  jB 
) (cos
(3.44)
p p p  p q pq pq pq pq
q1

de (3.44) se obtiene finalmente:

n
Pp   Vp Vq (G pq cos pq  Bpq sin pq )
q1
n (3.45)
Qp   Vp Vq (G pq sin pq  B pq cos pq )
q1

A partir de (3.41) y (3.45), P y Q para la barra p se pueden determinar como:

P  Pesp 
n V V (G 
cos B sin  )
p p  p q pq pq pq pq
q1
(3.46)
Q  Qesp  cos  )
n V V (G 
sin  B
p p  p q pq pq pq pq
q1

Por lo tanto se puede escribir a manera de resumen:

– Para las barras PQ y PV

P  Pesp 
n V V (G 
cos B sin  )
(3.47)
p p  p q pq
pq pq pq
q1

– Para las barras PQ

Q  Qesp 
n V V (G 
sin  B cos  )
(3.48)
p p  p q pq
pq pq pq
q1

– Para la barra slack no se requiere ecuaciones.

 En las ecuaciones anteriores, las magnitudes de las tensiones en las barras PV y slack al igual que el
ángulo en la barra slack no son variables, sino que se mantienen en sus valores especificados. Por lo tanto, el
sistema formulado incluye dos ecuaciones para cada barra PQ y una para cada barra PV. Las variables del
problema son V y  para cada barra PQ y  para cada barra PV.
 Por razones prácticas se da a la barra slack el número n y se colocan los primeros números a las
barras PQ. Luego si se tiene m barras PQ, se tendrá (n-m-1) barras con control de voltaje (barras PV).

La Ecuación (3.40) queda entonces:

  
P  J
1
(3.49)
   
V Q
   

Con P y Q calculados según (3.46). Luego, los valores actualizados para  y V son:
 k1   k   k 
 k1     
V k k (3.50)
V V
     

Despejando P y Q de (3.49), considerando (3.34) y arreglando adecuadamente, para hacer más

fácil el manejo de las ecuaciones, se tiene:

P    
 P    H N
P V
 
 
     V      
 V
 (3.51)
  Q Q V 
 V    M L   
Q   V    V 
V 
Si se tiene n nodos, m de carga, 1 libre y n-m-1 de voltaje controlado las dimensiones de las
submatrices que forman el Jacobiano son:

[H] es de (n-1) x (n-1) [N] es de (n-1) x m

[M] es de m x (n-1) [L] es de m x m

Según lo anterior, la matriz Jacobiana completa es cuadrada y de dimensión, [(n-1)+m] x [(n-1)+m].

A partir de (3.51), se obtiene que:

H pq  Pp Pp
N pq  Vq
q Vq
(3.52)
M pq  Qp Qp
q L pq  Vq

Vq
Considerando (3.46), se pueden determinar todos los elementos de la matriz Jacobiana como sigue:

– para p  q

H pq  Vp (G pq sinpq  Bpq cos pq )Vq

N pq  Vp (G pq cos pq  Bpq sinpq )Vq
M pq  Vp (G pq cos pq  Bpq sinpq )Vq (3.53)

Lpq  Vp (G pq sinpq  Bpq cos pq )Vq

Se puede apreciar, que por la forma de la ecuación (3.51): H pq  Lpq

 N pq  M pq

(3.54)
Para p = q

H pp  Vp2  Qp
Bpp
N pp  G V 2  P
pp p p
(3.55)
M pp  G pp Vp  Pp
2

Lpp  Bpp V 2  Q
p p

Las expresiones (3.54) y (3.55) muestran lo importante que fue el haber planteado la matriz Jacobiana
tal como se hizo en (3.51). Utilizando este tipo de coordenadas, el valor de Q en las barras PV puede ser
calculado luego que el proceso haya convergido.

El proceso mediante el método de Newton-Raphson para calcular los voltajes en las barras, se
muestra en la Figura 3.8, luego de lo cual se determinan los flujos de potencia y las pérdidas en las líneas.

Ingreso de Datos
del sistema

de P
Calcular
Elección y Q
Tensiones
Iniciales

Convergió la solución? Salida

Iter=0
no
Iter = Iter máx?
no si

Iter = Iter+1 si

Formación del Jacobiano

Inversión del Jacobiano

Cálculo de
 y V

Corrección de Valores de  y V
Figura 3.8.- Diagrama de Flujo del Algoritmo de Newton-Raphson
 En estricto rigor, la matriz Jacobiana se calcula e invierte en cada iteración. Sin embargo, en la
práctica se recalcula usualmente sólo un determinado número de veces en un rango de iteraciones, con el fin
de aumentar la velocidad al proceso iterativo.

Ejemplo 3.2. En el sistema del ejemplo 3.1 y, considerando los valores de voltajes obtenidos mediante el
método de Gauss-Seidel, determine el Jacobiano completo en la iteración 1

Solución: La matriz Jacobiana queda de la siguiente forma:


P P 
P1 V1 V1 
 1 1 1 H H N 

2  11 12 11
P P
P  
J     2 2
V1 2
 H
V1   21
H N 21

 22
1 2  
  M M12 L11 
Q1  11
Q1 Q1 V1 
 1 2 V1 

Los elementos de la matriz Jacobiana se obtienen utilizando las expresiones (3.53) para pq y (3.55) para
p=q; por lo tanto, para los elementos de la diagonal se deben determinar previamente P y Q, utilizando (3.45)
y considerando los resultados obtenidos en el ejemplo 3.1, que se resumen a continuación:

 7,5 5
 2,5   22,39 7, 15 
5 
G 4,1667 

B  7,5  12,4 5 

2,5 1,6667
  5  6,6667  15 5  19,89
1,6667 

V& 1  1,0401  V& 2  V& 3  1,060º

1,23º ; 1,040,23º

Con (3.45) se obtiene:

P1=0,6628; P2=0,1948; Q1=0,235; Q2=0,2865

Con (3.53) y (3.55), la matriz Jacobiana queda:

 24,4567 8,0412  7,4508 


J 8,1790  2,4967

13,6983
  8,7764  2,91  23,9867 

3.5.4. Método de Newton-Raphson desacoplado (MNRD)

Cuando se resuelven problemas que involucren un número considerable de barras, este método
representa una alternativa para mejorar la eficiencia computacional y reducir los requerimientos de memoria.
Se hace uso de una versión aproximada del método de Newton-Raphson completo visto en 3.5.3., basada en
las siguientes consideraciones:
 Un cambio en el ángulo del voltaje  en una barra afecta principalmente al flujo de potencia activa P
 en las líneas y débilmente a la potencia reactiva Q, lo que significa P Q
que: 
 
 Un cambio en el módulo del voltaje V en una barra afecta principalmente al flujo de potencia reactiva
Q en las líneas y débilmente a la potencia activa P. Por lo tanto: Q P
V  V
 De acuerdo con esto, en la ecuación (3.51) N y M se pueden despreciar frente a L y H
respectivamente y se puede escribir:

P 
P      
 0 H 0   
 
     Q V    V (3.56)
 0 V    0 L   
Q  V     V 
o bien: V 

P H (3.57)

Q  L VV (3.58)

Estas ecuaciones están desacopladas, en el sentido que las correcciones de los ángulos de los voltajes
se calculan usando sólo los cambios en la potencia activa mientras que las correcciones en la magnitud de los
voltajes se determinan sólo con los cambios en la potencia reactiva. Sin embargo, las matrices [H] y [L] son
interdependientes, porque los elementos de [H] dependen de las magnitudes de los voltajes que se están
resolviendo en la ecuación (3.58), mientras que los elementos de [L] dependen de los ángulos de la ecuación
(3.57). El procedimiento natural es ir resolviendo alternadamente los dos sistemas de ecuaciones usando en
uno, las soluciones más recientes del otro. Esta aproximación significa aumentar el número de iteraciones, lo
que en la práctica queda compensado por el menor tiempo en cada iteración, debido a la disminución del
tiempo ocupado en la inversión de las matrices [H] y [L] de menor dimensión que la del Jacobiano completo.

Ejemplo 3.3. A partir del ejemplo 3.2 escriba la matriz jacobiana del método de Newton-Raphson
desacoplado. Eliminando las submatrices N y M se tiene:

 24,4567 8,04124 0 
J 8,1790  13,6983 0


 0 0  23,9867 

3.5.5. Método de Newton-Raphson desacoplado rápido (MNRDR)

El método de Newton Raphson desacoplado, requiere la inversión de dos matrices en cada iteración.
Para evitar estos cálculos, se introducen más simplificaciones que se justifican por las características que
presentan los flujos de potencia en las líneas y los voltajes en las barras de un sistemas de transmisión bien
diseñado y operando apropiadamente. Estas características son:
 Las diferencias angulares pq = p-q entre los voltajes de dos barras típicas son, por lo general, tan
pequeños que:
cos pq = cos (p-q) = 1y sin pq = sin (p-q)  0 (3.59)
 Las susceptancias Bpq de las líneas son mucho más grandes que las conductancias G pq, por lo que:

Gpq sin (p-q)  Bpq cos (p-q) (3.60)

 La potencia reactiva Qp que se inyecta a cualquier barra p del sistema durante su operación normal es
mucho menor que la potencia reactiva que fluiría si todas las líneas de la barra estuvieran en
cortocircuito con la referencia. Esto es:
2
Q p  Bpp
Vp (3.61)
Con estas aproximaciones, y usando (3.45) se tiene:

Pp  Vp n
 (G pq cos pq  Bpqsinpq ) (3.62)
Vq
q=1

Es decir:
P 
G V2V n (G 
cos B sin ) V
(3.63)
p pp p
p  pq
pq pq pq q
q1
qp

Q  V
n B V  B V2V n B V
 (3.64)
p  pq q
pp p p pq q
p q=1
q=1 qp

A partir de (3.63) y (3.64) y considerando (3.46), los elementos de [H] y [L] se obtienen como sigue:

P Qp
H pq  p  Bpq Vp L pq   Bpq Vp Vq (3.65)
Vq
q Vq Vq
Pp
H  B
V2V n B VB V2Q
(3.66)
pp
pp p
p  pq q
pp p p
p q1

L
V  Qp B
2 V2V n BVB V2Q
pp p pp p  pq q pp p p (3.67)
V p
q1
p qp

según (3.61), las ecuaciones (3.66) y (3.67) quedan:

H pp  L pp  BppV 2 (3.68)
p

Considerando (3.57) y (3.58), para el nudo p se puede escribir:

n
Pp  H pp p   H pq q (3.69)
q1
qp

V n Vq
Qp  Lpp p
 L (3.70)
q1 pq Vq
Vp qp

Introduciendo (3.65) y (3.68) en (3.69) y (3.70) se obtiene:

Pp  Vp (Bpp Vp p  n

 Bpq Vq q (3.71)
)
q1
qp
Vp Vp Vp
Vp n
Qp  (Bpp  
q1 qp Vq
Bpq Vq q (3.72)
V )

Dividiendo ambas ecuaciones por Vp y simplificando los voltajes de la ecuación (3.72), las
ecuaciones anterires quedan:
Pp n
 Bpp Vp   Bpq Vq q (3.73)
Vp q1
p qp

Q p
n
Vp  Bpp V   Bpq Vq (3.74)
p q1
qp

Se ha reducido la no linealidad de las ecuaciones, lo que debería permitir mejorar el comportamiento

del proceso iterativo. Si fijamos arbitrariamente los V p y Vq del segundo miembro de la ecuación (3.73) en
1,0 (pu), los coeficientes de ambas ecuaciones quedan iguales, es decir:

Pp
n
Vp  Bpp p   Bpq q (3.75)
q1
qp

Q p n
Vp  Bpp V   Bpq Vq (3.76)
p q1
qp

Matricialmente se tiene:

P 
 B'
 (3.77)
V
Q
 B''V (3.78)

V

donde B’ es el negativo de la parte imaginaria de la matriz Y B, sin la fila y columna correspondiente a

la barra slack y B” es el negativo de la parte imaginaria de la matriz Y B sin las filas y columnas
correspondientes a las barras slack y PV

Una estrategia típica de solución es la siguiente:

1. Calcular los errores iniciales P/V

2. Resolver la ecuación (3.77) para determinar 
3. Actualizar los ángulos  y usarlos para calcular los errores Q/V
4. Resolver (3.78) para calcular V
5. Actualizar las magnitudes del voltaje V y usarlas para calcular los nuevos errores P/V
6. Volver al punto 2. y repetir el proceso hasta que todos los errores estén dentro de la tolerancias
especificadas.

Ejemplo 3.4. Para el ejemplo 3.1 escriba la matriz Jacobiana del método de Newton-Raphson desacoplado
rápido.

Solución: A partir de los valores de la matriz [B] del Ejemplo 3.2 se tiene:

 7,5 y
22,39 12,4 22,39
B'  7,5

B" lo tanto : 7,5  7,5 0 

22,39; por J    12,4 0 
  0 22,39
 0
 3.6. Flujo de Carga “DC” o Flujo de Potencia de Corriente Continua

Recibe este nombre no porque se aplique a sistemas que operan con corriente continua. El método es
aplicable a los Sistemas Eléctricos de Potencia que operan a corriente alterna haciendo ciertas simplificaciones que
permitan analizar la red como si fuera de corriente continua, teniendo la ventaja de que se trabaja con números
reales, simplificando mucho los cálculos. Sin embargo, sólo se obtienen los flujos de potencia activa.

El análisis considera la base Z pq ( Ypq )

S pq
fundamental de los flujos de potencia en p q
corriente alterna ya estudiados. I pq

Vp Y'pq /2 Y'pq /2 Vq
Consideremos la línea conectada
entre los nudos p y q de la Figura 3.9.

La potencia compleja que va desde

la barra p a la barra q, es: Figura 3.9.- Modelo de una línea de transmisión para el
cálculo de flujo de potencia DC
S& pq
 *
 (V  V )Y pq V
'

2) *   jQ pq (3.79)
&
V p V & & & & P
I pq p p
q p (Ypq pq

Sea:

V& p  V& q  Vq Y& pq  G pq  j B pq y '

Ypq jB' (3.80)
pq
Vp  p ;  q ;

R pq X pq
G pq  R 2  X 2 y Bpq   R 2  X 2 (3.81)
pq pq pq pq

Entonces:

S& pq  Vp  p (Vp    p  Vq    q ) (G pq  jB pq )  pjV 2 pq

(B '2) (3.82)

Es decir:

S& pq  V
p 
2

 Vp Vq (cos  pq  j sin  pq ) (G pq  jB pq ) pjV 2pq(B ' (3.83)
2)

Por lo que la potencia activa entre las barras p y q es:

Ppq  (V 2  V V
p p q cos  ) pq B pq
V Vp q sin  (3.84)
G pq pq

Suposiciones simplificatorias

a) Vp = Vq = 1,0 (pu)
b) Rpq  Xpq
c) pq = p - q es muy pequeño

A partir de b) y considerando (3.81) se tiene que:

R pq X pq
G pq  R 2  X 2
 R 2  X 2  G pq  0 (3.85)
B  pq pq

 1
X pq pq 

pq
(3.86)
R 2pq  X 2pq X pq
 Según c)

cos  pq  cos (p  q )  1

sin pq  sin (p  q )  p  q radianes

Por lo tanto (3.84) queda:

p  q
Ppq  Bpq (p  q )  (3.87)
X pq

Para una red eléctrica que tiene n nudos las potencias netas inyectadas en cada nudo se pueden
escribir:
n p  q
P  nP
 (3.88)
p  pqX q1 pq
q1

Donde q incluye todos los nodos directamente conectados al nodo p. Luego para los n nudos se tiene

 P1  B11 L L B1n 1 


   B1p
M M M M M  
M M
     
Pp   L Bpp L Bpn  p  (3.89)
Bp1
M  M M M M   M
     
 P  Bn1 L Bnp L Bnn    
 n   n 

o sea:

 P  (3.90)
B
en que:

P: Vector de potencias activas netas en las barras

B: Matriz de susceptancias del sistema
: Vector de ángulos de las tensiones de las barras
Los elementos de la matriz B se determinan como sigue:

B pp n
1
 (3.91)
q1 X pq

Con p  referencia y donde q incluye todas las líneas conectadas al nudo p

Bpp  0
con p  referencia (3.92)
Bpq 1
 con p y q  referencia (3.93)
X
pq

Bpq  0 con p  referencia ó q  referencia (3.94)

 De acuerdo con lo indicado en las expresiones (3.91) a (3.94), la matrizB, tal como se planteó en
(3.89) es singular ya que la fila y columna correspondientes a la barra de referencia contienen solamente
ceros. Así, para un sistema de n barras sólo hay (n-1) ecuaciones linealmente independientes, por lo que se
debe considerar la submatriz de B con dimensión (n-1)x(n-1), la que puede invertirse para solucionar el
problema del cálculo de los ángulos de fase de las tensiones de las barras. A partir de (3.90) se obtiene
entonces:

  B-1 P (3.95)

Como se aprecia, el procedimiento es muy sencillo, ya que el sistema es lineal y por lo tanto todos los
ángulos de las tensiones de las barras se pueden determinar utilizando (3.95) para luego mediante (3.87),
determinar los flujos de potencia activa en las líneas.

Ejemplo 3.5. Resuelva el Ejemplo 3.1 utilizando el método de Flujo DC.

Solución: A partir de los parámetros dados en el ejemplo 3.1 se tiene:

B1 1 1 1
1
  1  25; B22  
1
 B1  B21    8,3333
0,12 0,06 0,12 0,1 13,8889; 2
0,1
8 2

B33  B31  B13  B23  B32  0

Por otra parte: P1=-0,6; P2=0,2; luego, la ecuación matricial queda, según (3.90):

 P1   25 
 1 
 8,3333
  8,3333 13,8889
  
P2 2
     

Resolviendo la ecuación se obtiene:

    25
 8,3333 -1  P   0,05 0,03  0,6   0,024
1  1      
2   13,8889 P
 2  0,03 0,09   0,2   0 
8,3333 

Con estos valores se determinan los flujos de potencia activa en las líneas, utilizando (3.87)

1  2  0,024  0
P    0,2
0,12
1  3  0,024  0
P    0,4
0,06
2  3 00
P23   0
X 0,18

La potencia que entrega el generador en la barra slack es:

PG3  P32  P31  0  0,4  0,4

La potencia recibida por la carga es:

PC1  P21  P31  0,2  0,4  0,6

 Problemas propuestos

3.1. En el sistema mostrado en la Figura 3.10, empleando el método de Gauss Seidel-Y B, calcular la tensión
en las barras 2 y 3 en la segunda iteración ( V& y V& ) y la potencia compleja entregada por el generador
2 2

en 2 3
estas condiciones. La tensión en la barra 1 se mantiene constante en su valor nominal y los datos en % están
dados considerando SB=100 MVA.

3.2. En el sistema de la Figura 3.11 y empleando el método de Gauss-Seidel-Y B, determinar:

a. La tensión en la barra 3 considerando que k1
V3  Vk3  0,01
b. Las potencias activa y reactiva suministradas por el Generador G y las pérdidas de potencia activa y
reactiva del sistema.

1 2 1
G 3 X = 0,1 2 3
Z=j50% G
. Y/2=j1,5% 0,4 + j0,5
X=10%
SG X= 0,5
1,05 0º
. .
S=(8+j4) S=(25+j18) MVA (20 + j12) MVA
MVA
Figura 3.10 Figura 3.11

3.3. Para el sistema de la Figura 3.12, los datos en pu, base común, se dan en las Tablas Nº 1 y Nº 2.
a. Determinar el voltaje en todas las barras, haciendo dos iteraciones con el método de Gauss-Seidel Y B y
con los valores obtenidos, determinar los flujos de potencia en todas las líneas.
b. Determinar los elementos del Jacobiano desacoplado rápido del sistema.
c. Correr un flujo DC y determinar los ángulos de los voltajes, los flujos de potencia activa en las líneas y
la potencia activa entregada por el generador G3

SC3 Tabla Nº 1: Datos de las líneas

G3 G2
Línea Z (pu) Y’/2 (pu)
3 2 1-2 0,04+j0,10 j0,06
1-3 0,02+j0,08 j0,07
2-3 0,06+j0,20 j0,04

Tabla Nº 2: Datos de las barras

1 Barra Nº Tipo V (pu) PG QG PC QC
SC1 1 PQ - - - 0,8 0,2
2 PQ - 0,2 0,1 0,0 0,0
Figura 3.12 3 SL 1,06 - - 0,4 0,1

3.4. La Figura 3.13 muestra un sistema de tres barras. La barra 1 es la barra slack y su voltaje es 1,05 (pu), la
barra 3 es de tensión controlada (BTC) y el módulo del voltaje en ella está especificado en 1,05 (pu). La
potencia reactiva del condensador síncrono puede variar entre 0 y 0,6 (pu). Determinar:
a. Las tensiones en las barras, haciendo dos iteraciones del método de Gauss-Seidel-YB
b. Los flujos de potencia activa y reactiva en las líneas
c. Las pérdidas de potencia activa y reactiva
d. La potencia compleja entregada por el generador G1
e. La potencia reactiva suministrada por el condensador de la barra 3.
3.5. En el sistema de la Figura 3.13, suponga que los valores de las tensiones en una iteración cualquiera son:
V&  1,05 0º; V&  1,04-2,5º; V&  1,05-1,5º. Determinar:
1 2 3
a. El Jacobiano completo en esa iteración
b. El Jacobiano desacoplado
c. La matrices [B’] y [B”]

G1 G2
S=1,0+j1,4
1 2
j
0,08
S=2,0+j1,0 S=0,4+j0,4
j 0,05 j0,04

3
C
S=1,2+j0,8

Figura 3.13

3.6. En el sistema de dos barras de la Figura 3.14, los datos en pu están en base 100 MVA, la barra 1 es la
barra slack y su voltaje es 1,04 (pu). En la barra 2, hay un compensador síncrono (CS) que mantiene
constante el módulo del voltaje en 1,05 (pu).
a. Determinar el voltaje en la barra 2 haciendo dos iteraciones del método de Gauss-Seidel-Y B y con los
valores anteriores, calcular la potencia reactiva que debe entregar (recibir) el compensador síncrono y las
potencias activa y reactiva que debe entregar (recibir) el generador.
b. Determinar los elementos de la matriz Jacobiana completa, considerando que los voltajes en la barra 1 y
2 son de 1,04 0º y 1,05 -4.5º respectivamente.

PG1 + jQ G1 -0,4  QG2  0,4 QG2

G CS

1 2
j0,1

Sc1 =(40+j20) Sc2 =(100-j40) MVA

MVA
Figura 3.14

3.7. En el sistema de la Figura 3.15, la barra 3 es 1 X = 0,1 2 3

BTC con una tensión de 1 (pu). Empleando el G 0,4 + j0,5
método de Gauss-Seidel-YB, determine:
a. La tensión en la barra 3, considerando que C
X= 0,7 0º
k1  k  0,5 1,05
0  QC  0,4
3 3 (20 + j12) MVA
b. Las potencias activa y reactiva
suministradas por el Generador G y las Figura 3.15
pérdidas de potencia activa y reactiva del
sistema.
c. La potencia reactiva que debe entregar el condensador ubicado en la barra 3. Los datos están en pu base
100 MVA.
 CAPITULO 4

ESTUDIO Y ANALISIS DE FALLAS

4.1. Introducción

Las condiciones anormales de funcionamiento de un Sistema Eléctrico de Potencia (SEP), se deben a

fenómenos transitorios, que se pueden clasificar, según al tiempo de duración en las siguientes categorías:

– Fenómenos transitorios ultrarápidos: Corresponden sustancialmente a descargas atmosféricas sobre

las líneas de transmisión y a los fenómenos producidos por operaciones de conexión y desconexión de
diversos componentes de la red del SEP, tales como, las líneas. Las perturbaciones de este tipo dan
origen a ondas de tensión y corriente que viajan prácticamente a la velocidad de la luz, pero su efecto
dura unos pocos milisegundos después de iniciado. Sin embargo, los procesos de reflexión de las
ondas producen elevadas tensiones que pueden llegar a destruir el equipo asociado a las líneas. La
razón del estudio de estos fenómenos radica en el hecho de que su análisis suministra las bases
necesarias para la selección adecuada del nivel de aislación de los equipos eléctricos asociados a las
líneas y de las líneas mismas.

– Fenómenos transitorios medianamente rápidos: En este grupo se incluyen los fenómenos causados
por cambios abruptos de la estructura del SEP, o sea los cortocircuitos o líneas abiertas. Usualmente,
sólo los 10 primeros ciclos son de importancia práctica y se estudian en el rango de 10 a 100
milisegundos siguientes a la falla.

– Fenómenos transitorios lentos: Cuando ocurre un cortocircuito en una línea de transmisión

importante y no se desconecta oportunamente la sección afectada, puede producirse uno de los
fenómenos más peligrosos de un SEP, esto es, oscilaciones mecánicas de los rotores de los
generadores. Se producen fenómenos transitorios electromecánicos que se estudian bajo el nombre de
estabilidad transitoria. Las oscilaciones mecánicas de los rotores son relativamente lentas, en
consecuencia, los estudios de estabilidad transitoria se realizan en el rango de fracción de segundo
hasta un minuto.

Debido a los fenómenos transitorios se pueden producir en un SEP, diversas alteraciones que reciben
el nombre de fallas. Una falla en un circuito es cualquier evento que interfiere con el flujo normal de
corriente. Sin embargo, dentro de este curso, designaremos como fallas a los cortocircuitos y las fases
abiertas.

4.2. Tipos de fallas

– Cortocircuitos: Trifásico simétrico, aislado o a tierra, bifásico aislado (cortocircuito entre 2 líneas),
bifásico a tierra (entre dos líneas y el conjunto a tierra) y monofásico (una línea conectada a tierra).

– Fases abiertas: Una fase abierta, dos fases abiertas y tres fases abiertas. La última situación significa
que la línea o dispositivo sale completamente de servicio.

Los cortocircuitos trifásicos dan origen a fallas simétricas pues el SEP permanece eléctricamente
balanceado, en cambio los cortocircuitos bifásicos aislado y a tierra y el monofásico, así como 1 ó 2 fases
abiertas corresponden a fallas asimétricas, ya que el sistema queda eléctricamente desbalanceado en el punto
de falla.

En el caso de fallas simétricas, el cálculo se realiza en base a una representación monofásica (por
fase) de la red del SEP y se aplican las técnicas normales de análisis de circuitos. Para el cálculo de las fallas
asimétricas, resulta conveniente utilizar al Método de las Componentes Simétricas.
 51
4.3. Cálculo de Cortocircuitos

En general las corrientes de cortocircuito alcanzan magnitudes mucho mayores que los valores
nominales de los generadores, transformadores y líneas. Si se permite que estas corrientes circulen por un
período prolongado, pueden causar un serio daño térmico al equipo y problemas de estabilidad de
funcionamiento en el SEP.

En este aspecto, el tipo de cortocircuito más severo es el trifásico, el que además de dar valores
elevados de corriente, reduce a cero la capacidad de transmisión de una línea, lo siguen los cortocircuitos
bifásico y finalmente el monofásico. En cambio, el tipo mas frecuente es el monofásico (aproximadamente el
75% de los casos ) y el menos frecuente es el trifásico (aproximadamente el 5% de los casos).

En muchas oportunidades las corrientes de cortocircuito se autoextinguen y se restablece la aislación.

Debido a este hecho, se utilizan en la práctica interruptores que reconectan automáticamente la línea dañada,
una, dos o más veces para probar si la falla se ha eliminado. Sólo en el caso de que la falla persista, el
interruptor desconecta la línea en forma definitiva.

4.3.1. Objetivos

 Definir la capacidad de ruptura de los interruptores necesarios en las diversas partes de un SEP, para
lo que se realiza normalmente un cálculo de cortocircuito trifásico simétrico, debido a que este tipo de
falla produce las corrientes de cortocircuito más elevadas en la mayoría de los casos.
 Ayudar a establecer un sistema adecuado de protección para diversas condiciones de falla, para lo que
se debe realizar un cálculo de distribución de corrientes en la red del SEP tanto para cortocircuitos
simétricos como asimétricos (usualmente el cortocircuito monofásico).
En general, el Cálculo de Cortocircuitos debe proporcionar los siguientes resultados:
 La corriente en el punto de falla

 La potencia de cortocircuito en el punto de falla

 La distribución de corrientes post-falla en todas las líneas del SEP
 Las tensiones post-falla en todas las barras

4.3.2. Aproximaciones

 Las máquinas síncronas se representan por los circuitos equivalentes aproximados, que se muestran
en la Figura 4.1.
jX'd jX'd
+ +
+ jX"d + jX"d
E'
E V V
E"
- -
- -

Figura 4.1.- Circuitos equivalentes para las Máquinas Síncronas

 Las cargas, cuando se estima necesario incluirlas, se suponen independientes de la tensión y se

representan por una impedancia o admitancia equivalente.
 Todas las tensiones internas de los generadores se suponen iguales entre sí e iguales a 1,0 (pu)
 Se desprecian las corrientes de pre-falla
 En muchos casos se desprecian las resistencias de los elementos y sólo se consideran sus reactancias
 Los transformadores con cambio de Tap se consideran en su razón nominal

4.4. Cortocircuitos trifásicos simétricos

4.4.1. Comportamiento de un generador en condiciones de cortocircuito trifásico simétrico

a. El generador en vacío antes de producirse la falla: La corriente que circula por cada fase del
generador en cortocircuito, es similar a la que circula por un circuito R-L serie, alimentado bruscamente por
una fuente de tensión sinusoidal; es decir, la corriente es asimétrica respecto al eje de tiempo y disminuye en
forma exponencial. Sin embargo, existe una diferencia fundamental y ella radica en que la reactancia del
generador no permanece constante durante el fenómeno (Figura 4.1). Las corrientes en las tres fases de un
generador en cortocircuito se muestran en la Figura 4.2.

ia

ib

ic

if

Tiempo
Figura 4.2.- Corrientes de cortocircuito en un Generador Síncrono

Periodo subtransitorio Periodo transitorio Régimen permanente

2 I"

2 I'
2I

Envolvente transitoria
Tiempo
extrapolada

Valor de régimen permanente Envolvente

extrapolado Real

Figura 4.3.- Corriente de cortocircuito en un Generador Síncrono despreciando la

componente unidireccional
 Usualmente la corriente continua no se considera en el análisis y su efecto se incluye posteriormente
en el cálculo de las corrientes instantáneas y de interrupción de los interruptores. Despreciando el efecto de la
componente continua, la corriente de cortocircuito de una fase cualquiera, resulta simétrica, como se muestra
en la Figura 4.3, que corresponde a un generador con enrollados amortiguadores y en vacío antes de
producirse la falla. Directamente de esta figura los valores eficaces de corrientes de cortocircuito quedan:
E
Corriente subtransiente I"  (4.1)
X "d
E
Corriente transiente I'  (4.2)
X 'd
E
Corriente permanente I (4.3)
Xd
b. El generador con carga antes de producirse la falla: En este caso, la fuerza electromotriz (fem)
interna E se va modificando a medida que transcurre el fenómeno y, para determinar las corrientes
subtransiente y transiente de cortocircuito se deben considerar los circuitos mostrados en las Figura 4.4 y 4.5,
respectivamente, donde Ze es una impedancia externa que puede existir entre los terminales del generador y
el punto de Falla F y Zc es la impedancia del consumo.

jXd" Ze I" jXd' Ze

F F I'
+ +
E" Zc Zc
E'
- -

Figura 4.4.- Circuito Subtransiente Figura 4.5.- Circuito Transiente

Para realizar el cálculo se pueden emplear dos procedimientos

– Aplicando el teorema de Thevenin en el punto de Falla: En la Figura 4.4, sean:

V& F (0) : tensión en el punto F antes de producirse la falla

&I" : corriente subtransiente de cortocircuito
ZTH : Impedancia equivalente de Thevenin calculada desde el punto de falla, donde:
(Ze  jX" ) Z
d c (4.4)
ZTH 
Ze  jX d Zc"

Por lo tanto, el circuito de la Figura 4.4, se transforma en el de la Figura 4.6, siguiente

I'' F
ZTH
+
VF(0)

Figura 4.6.- Circuito equivalente de Thevenin en Régimen subtransitorio

Luego: (
0  ZTH
)
&I"  V& F
 d
V&
(Ze  jX" d) Z c(0)

Ze  Zc  jX" F
(4.5)

de la misma forma, para la Figura 4.5 se tiene:

&I'  V& (0)

F Ze  Zc  jX' d &
V (4.6)
 (Ze  jX' )d Z c(0) F
ZTH

– Empleando las tensiones detrás de las reactancias subtransiente o transiente: Cuando circula una
corriente de carga Ic antes de la falla, se pueden visualizar tres tensiones internas posibles, asociadas a
sus correspondientes reactancias, tal como se indicó anteriormente. Las Figuras 4.7 a) y b) muestran
los circuitos equivalentes y los diagramas fasoriales respectivos.
E
jXd"
I" ( I', I)
E'
+
jX' Ic
E" + d E"
E' jXd
V jX"d Ic jX' I jXdIc
E-  dc
V
-
Ic
a) b)

Figura 4.7.- Cálculo de Corriente de Cortocircuito, empleando las tensiones internas. a) Circuito Equivalente
b) Diagrama

Fasorial A partir de la Figura 4.7 se puede escribir:

E& "  V& d &I (4.7)

jX " c (4.8)
E& ''  V&  &I
jX
dc

E&  V&  jX d (4.9)

&I
c

Por lo que, las corrientes de falla son:

&I" 
E& "
jX  Z
" (4.10)
d e

&I'  E& '

jX  Z
' (4.11)
d e
&I  E&
(4.12)
jXd  Ze

c. Concepto de Potencia de Cortocircuito: Durante un cortocircuito trifásico simétrico en un SEP, las

tensiones en las barras no falladas disminuyen. La magnitud de la caída de tensión en las barras es una
indicación de la capacidad de SEP para reaccionar frente al cortocircuito. Es conveniente disponer de una
medida de esta propiedad del sistema como asimismo de la severidad de la falla. Ambos objetivos se pueden
cumplir definiendo una cantidad denominada "Potencia de cortocircuito", "Capacidad de cortocircuito", o
"nivel de falla" de la barra fallada.

Considérese una barra “p” cualquiera del SEP en la cual se ha producido un cortocircuito trifásico
simétrico.

+ Vp (0) p
Ipf
- Ref.
Figura 4.8.- SEP con un cortocircuito trifásico en la barra p

Sean: Vp(0) : tensión en la barra p antes de producirse la falla

Ipf : corriente de cortocircuito o de falla en la barra p.

Entonces, por definición, la potencia de cortocircuito S cc en la barra p será:

Scc  Vp 0 I pf (4.13)

Por otra parte si VB e IB son, respectivamente, el voltaje base y la corriente base en el sector que
corresponde a la barra p, se puede demostrar que:
Vp (0)2 pu 
  I pu  S
Scc  Vp (0) pu  (4.14)
pf B Zpf pu SB

en que Zpf es la impedancia de cortocircuito en la barra “p” y corresponde a la impedancia

equivalente de Thevenin calculada desde la barra p hacia el interior del SEP.

d. Algunos antecedentes relativos a la selección de interruptores: Los valores de corriente

suministrados por un cálculo de cortocircuito, corresponden a corrientes simétricas respecto al eje del tiempo
y por lo tanto no incluyen la componente de corriente continua. En la selección de interruptores debe tenerse
en cuenta la componente de corriente continua; por ello se distinguen dos valores de corriente:

 Corriente instantánea: Es la corriente que debe soportar un interruptor inmediatamente después de

ocurrida la falla. Para determinarla, se calcula en primer lugar la corriente simétrica de cortocircuito
utilizando las reactancias subtransientes de los generadores, motores sincrónicos y de inducción.
Luego, el valor así calculado, se multiplica por un factor que depende de la tensión de operación del
interruptor. Los factores usualmente empleados se indican en la Tabla 4.1.

 Corriente de interrupción: Es la corriente que un interruptor debe ser capaz de interrumpir en el

momento que se abren sus contactos. Para determinar su valor, se procede primero a calcular la
corriente simétrica de cortocircuito y luego se aplica un factor que depende de la velocidad de
operación del interruptor. La Tabla 4.2 muestra algunos valores típicos.

Para el cálculo se recomienda emplear las reactancias subtransientes de los generadores, las
reactancias transientes de los motores y condensadores síncronos. Los motores de inducción no se
consideran.
Tabla 4.1.- Factores de corrección para la corriente Tabla 4.2.- Factores de corrección para
instantánea la corriente de interrupción
Voltaje Factor Interruptores de Factor
 5 kV 1,6 8 ciclos o mas 1,0
 5 kV 1,5 5 ciclos 1,1
 600 V 1,25 3 ciclos 1,2
2 ciclos 1,4

Ejemplo 4.1. En el sistema de la Figura 4.9, los Barra del

datos en pu están en base 2.000 kVA y 480 Volt. Si Sistema m
j0,023 ohm/fase M1 X''=0,8 (pu)
ocurre un cortocircuito trifásico en la barra m,
cuando el voltaje en esta barra es el nominal y los Línea M2 X''=0,25(pu)
motores están trabajando en condiciones nominales,
determinar la corriente subtransitoria de falla VN =480 V 1 MVA, 480 Volt,
(Amperes) en cada uno de los motores y en la línea, SCC3 =9,6 MVA cos  = 0,9 ind.; cada motor
así como la corriente de falla en la barra m,
considerando la corriente de pre-falla y utilizando: Figura 4.9
a. El método de las tensiones internas
b. El método de Thevenin SCC3: Potencia de Cortocircuito trifásico en barra del
Solución: sistema.

a) Cálculos previos:
9,6 jXL
S  2 MVA; V  0,48 kV  S   4,8 (pu) jXS j0,8
m
B B CC3
2 IS I1
V(0) 2 1 I2
XS    0,2083 (pu) + j0,25
S CC3 4,8 E"S Z TH + +
kV 2 0,482 E"
Z B B   0,1152 ; E"M2
M1
MVA B 2
0,023 2 10 6
XL   0,1997 (pu); I B   2405,63 A
0,1152 3  480
Figura 4.10
b) Método de las tensiones internas: Figura 4.10
1
& 4800º & & 1/ 2  cos 0,9
Vm   10º I1  I2   0,5  25,84º (pu)
480 (pu) 10
º
E M1  10ºj0,8  0,5  25,84º  0,9  23,56º (pu);E M2  10ºj0,25  0,5  25,84º  0,9522  6,79º
" "

(pu)
E "S  10º j(0,2083  0,1997) 1,0  25,84º  1,233717,32º (pu

En falla, Vm=0; por lo tanto, del circuito (Figura 4.11): jXS jXL j0,8
m
&I "  1,233717,32º  3,024  72,68º (pu) I"S I"M1
S
j(0,2083  0,1997) + I"M2 j0,25
&I "  0,9  23,56º  1,125  113,56º (pu) E"S I"F +
+
E"M1
M1
j0,8 E"
M2
&I "  0,9522  6,79º
 3,8088  96,79º (pu)
M2
j0,25
&I "  &I "  I "  &I "  7,7  90º (pu)
F S M1 M2
Figura 4.11
Valores en amperes: Multiplicando por la corriente base se obtiene:

&I "  7.274,63  72,68º A;  2.706,33  113,56º A;  9.162,56  91,79º A;

S  18.523,35  90º A
&I " &I " &I "
M M2 F
1

d) Usando el Teorema de Thevenin

A partir de la Figura 4.10, cortocircuitado las fuentes de tensión se obtiene la Impedancia de Thevenin:

ZTH  1
 j0,1299 (pu). Además, VTH=Vm(0)=10º (pu)
1 1 1
 
j0,408 j0,8 j0,25

V& 10º
&I "F 
 j0,1299  7,7  90º
TH
Z& j0,8  j0,25
TH

&I " j(0,8  0,25)

S  7,7   1  25,84  3,024  72,68º (pu)
90º j0,8  j0,25
j0,408 
j(0,8  0,25)
j0,408  j0,25
& "  7,7  j(0,408  0,25)
M  0,5  25,84  1,1257  113,56º (pu)
I 90º
1 j0,408  j0,25
j0,8 
j(0,408 
0,25) j0,408 
& "  7,7  j0,8 j(0,408 
M2  0,5  25,84  3,8083  96,79º (pu)
I 90º 0,8)
j0,408  j0,8
j0,25 
j(0,408  0,8)

4.4.2. Cortocircuitos trifásicos simétricos en un SEP

a. Método tradicional: Como en el caso de un cortocircuito trifásico simétrico, el SEP queda

balanceado, es posible trabajar utilizando el circuito equivalente por fase, con las aproximaciones usuales,
aplicando Thevenin en el punto de falla. El método es cómodo para resolver problemas con pocos nudos; sin
embargo, cuando se trata de sistemas de mayor tamaño, resulta poco práctico. Por otra parte, para calcular un
cortocircuito en otra barra es necesario hacer de nuevo todos los cálculos. Adicionalmente, la determinación
de los voltajes en las otras barras y el cálculo de las corrientes en las líneas significa resolver la red completa
del SEP.

b. Cálculo sistemático (Método general): Cuando se trata de sistemas de gran magnitud, los cálculos
manuales resultan demasiado engorrosos y se debe recurrir al uso de los computadores digitales. El
procedimiento que se sigue, en vez de calcular las corrientes en el punto de falla, para luego repartirlas en
todo el sistema; consiste en calcular directamente las tensiones en los distintos nudos, con ayuda de un
modelo nodal de impedancias. Conocidas las tensiones durante la falla, pueden calcularse a continuación las
corrientes por las diversas ramas. Debido a la rapidez del cálculo digital, la matriz de impedancia puede por
ejemplo, incluir las admitancias paralelo tales como las asociadas a las cargas.

Las tensiones, post-falla se pueden obtener como la superposición de la situación pre-falla (obtenida
normalmente de un cálculo de flujo de potencia) con la situación durante la falla solamente, es decir:
 VB (f )  VB (0) VB (4.15)
(cc)
donde: V (f : Vector de tensiones post-falla
B
) :Vector de tensiones prefalla
VB (0
)
VB (cc) :Vector de tensiones debido sólo a la falla:

V B
(f )  V& ···
&
V& ··· V& ··· V nf (4.16)
V B
···
pf
··· ···
qf
 t
(4.17)
1f V& V&
(0)  V& p0 q0
V& n0 t

10
 t

VB (cc)  &

··· V& ··· V& ··· V ncc (4.18)
1cc
 V& pcc qcc

Aplicando el método de resolución nodal a la red del SEP, después de falla se tiene:

IF   YB VB

 VB (cc)  ZB IF (4.19)
(cc)

En que IF  es el vector de corrientes (de falla) inyectadas en las distintas barras y ZB  es la matriz
de impedancia de barras que corresponde a la inversa de la matriz de admitancia de barras YB  ; definidas
como:

I F   &I1f &I &I &Inf t (4.20)

pf qf
···
··· ···

Z11 ··· Z1p ··· Z1q ··· Z1n 

 
M M M M 
 Zp1··· Zpp ··· Zpq ··· Zpn 
 
Z B  = M M M M (4.21)
 
Zq1··· Zqp ··· Zqq ··· Zqn 
 
M M M M 
 Z nn ··· Z np ··· Z nq ··· Z nn 

En realidad no se inyecta corriente en ninguna de las barras, sino que se extrae corriente
exclusivamente desde la barra fallada (por ejemplo, la barra p); por lo tanto, sólo uno de los elementos del
vector de corrientes inyectadas es distinto de cero y vale  &I pf . Introduciendo (4.19), (4.20) y (4.21) en
(4.15), considerando (4.16), (4.17) y lo dicho recién respecto a la corriente se obtiene:

V& 1f
 V& 10  Z1p &I pf
V& pf
 V& p0  Z pp &I pf
V& qf  V&  Zqp &I (4.22)
q0 pf

V& nf  V& n0  Z np &I pf

Si existe una impedancia de falla ZF entre la barra fallada p y tierra se tiene:

V& pf  Z &I (4.23)

F
pf

Reemplazando (4.23) en la p-ésima ecuación de (4.22) se obtiene finalmente:

V& p0
&I  (4.24)
pf
ZF  Zpp
 expresión que permite calcular la corriente en la barra fallada. Así mismo, el voltaje en esta barra es:

V&  1 
 Zpp (4.25)
 
ZF  Zpp  
p0
pf 

De la misma forma se puede obtener el voltaje en cualquier otra barra y la corriente de falla en una
línea cualquiera conectada entre las barras p y q cuya impedancia es z pq:

V& qf  Zqp
 V& ZF  Zpp V
& (4.26)
q0 p0
V&  V&
&I  (4.27)
qf
pf
pqf
z pq
Ejemplo 4.2. En el sistema de la Figura 4.12, 3 220 kV 220 kV 2
todos los datos en por unidad, están en base j0,05 j0,04
j0,03
común. Determinar las corrientes en cada una de j0,05 j0,08
las líneas y en los generadores, cuando ocurre un G2
3
cortocircuito trifásico simétrico en la barra 2, G G
1
13,2 kV
estando el sistema en vacío y utilizando: 13,8 kV j0,05
j0,05 1 13,2 kV
a. El método tradicional j0,02 j0,03
b. El método general (Matricial) Figura 4.12
Solución:
a) Método tradicional: El circuito equivalente se muestra en la Figura 4.13. Para encontrar la impedancia de
Thevenin en la barra 2 es necesario reducirlo. La Figura 4.14 muestra el circuito anterior donde se ha
realizado una transformación de Delta a Estrella entre los nudos 1, 2 y 3. Los valores de la estrella
equivalente son:
0,05 * 0,08 0,04 * 0,08 0,05 * 0,04
X1   0,02353; X 2  0,01882; X3  0,01176
0,05  0,08  0,05  0,08  0,05  0,08 

0,04 0,04 0,04
A partir del circuito de la Figura 4.14, la impedancia equivalente de Thevenin en la barra 2 queda:
 (0,05  0,05  0,01176)  (0,03  0,02  0,02353) 
  0,01882  0,03  0,05
(0,05  0,05  0,01176)  (0,03  0,02  0,02353)  0,0353
XTH   
 (0,05  0,05  0,01176)  (0,03  0,02  0,02353) 
 0,01882  0,03  0,05
 
(0,05  0,05  0,01176)  (0,03  0,02  0,02353)
 j0,04 3 jX 3 jX2 2
3 2
I 2F I I
X2 2F
j0,05 j0,08 jX1
j0,05 j0,03 j0,05 j0,03
1 V ,Z 1 VTH, Z TH
TH TH
j0,02 j0,02
IG3 IG3
j0,05 I G1 IG2 j0,05 j0,05 IG2 j0,05
j0,03 I G1 j0,03
+ + + + +
+
1 0º 1 0º
1 0º 1 0º 1 0º
1 0º

Figura 4.13 Figura 4.14

El circuito equivalente de Thevenin queda tal como j0,0353
se muestra en la Figura 4.15, donde, debido a que el I 2F +
cortocircuito es directo, se tiene que V 2F=0 y, por lo +
1 0º V2F
tanto:

&I 10º
  28,3286  90º (pu)
2F
j0,0353 Figura 4.15

Considerando la Figura 4.14 se pueden determinar las corrientes y voltajes:

&I 10º0
G   12,5  90º (pu)
2 j0,08
&I  &I  &I  28,3286  90º12,5  90º  15,8286  90º (pu)
X2 F2 G2

&I j(0,05  0,05  0,01176)

G   &I X2  9,5472  90º (pu)
1
&I * j(0,05  0,05  0,01176  0,02  0,03 
X1 0,02353)
&I j(0,02  0,03  0,02353)
G   &I X2  6,2814  90º (pu)
3
&I j(0,05  0,05  0,01176  0,02  0,03 
*
X3 0,02353)
V& 1F  10º j0,05 * I G1  10º j0,05 * 9,5472  90º  0,52260º
(pu) V&  10º j0,1* I  10º j0,1* 6,2814  90º 
3F G3
0,37190º
& (pu) &
&I  V 1F  V 2F  0,5226  6,5325  90º (pu)
L12
j0,0 j0,08
&8 &
&I  V 1F  V 3F  0,5226  0,3719  3,014  90º (pu)
L13
j0,05 j0,05
&I V& 3F  V& 2F 0,3719
   9,2975  90º (pu)
L32
j0,04 j0,04

b) Método general: Considerando el circuito de la Figura 4.13, se determina la matriz Y B, que resulta:

 j12, j20   j0,03 j0,016 j0,0186

j52,5 5 9
 
YB   j12,5  j50 j25 
1 
 ZB   YB    j0,0169 j0,035 j0,0222

 j20 j25   j0,0186 3 j0,035 
j55 j0,0222

Utilizando las expresiones (4.24) a (4.27) con ZF=0 se obtienen:

&I
2F 10º
  28,3286  90º V& 2F  0
j0,0353 (pu);
j0,0169 &I
V&  10º  10º  0,52120º V&
1F (pu); j0,0353
0,52120º0 3F
  6,515  90º (pu);
&I
 10º  j 0,0222 0º0,37110º
10º  0,5212
  3,002  90º (pu)
0,37110º (pu) j0,0353
L1 L13
2 j0,08 j0,05
0,37110º0
&I   9,2775  90º
L3 (pu) j0,04
2

Las corrientes en los generadores se determinan aplicando la Ley de Kirchhoff de corrientes en los nudos 1, 2
y 3 respectivamente:
&I  &I &
G1 L12  I L13  6,515  90º3,002  90º  9,517  90º (pu)
&I  &I  &I &
G2 2F L12  I L32  28,3286  90º6,155  90º3,002  90º  19,1716  90º (pu)
&I  &I  &I  9,2775  90º3,002  90º  6,2775  90º (pu)
G3 L32 L13

4.5. Cortocircuitos Asimétricos

4.5.1. Componentes simétricas

a. Consideraciones generales: El cálculo de cortocircuitos asimétricos en un SEP, se realiza

normalmente empleando el método de las componentes simétricas, por lo que es conveniente iniciar este
estudio resumiendo algunos puntos fundamentales relacionados con su teoría.

El Método de las Componentes Simétricas se basa en el teorema de Fortescue. Se trata de un método

particular de transformación lineal que consiste básicamente en descomponer un conjunto de fasores
desbalanceados en otro conjunto de fasores de características tales que permitan un análisis más sencillo del
problema original. En el caso particular de tensiones y corrientes trifásicas desequilibradas, este método los
transforma en tres sistemas de fasores balanceados. Los conjuntos balanceados de componentes son:

 Componentes de secuencia positiva: formado por tres fasores de igual magnitud, desfasados 120º
entre si y con la misma secuencia de fase que el sistema original.

 Componentes de secuencia negativa: formado por tres fasores de igual módulo, con desfase de 120º
uno de otro y con la secuencia de fases opuesta a la de los fasores originales.

 Componentes de secuencia cero: formada por tres fasores de igual módulo y con desfase nulo.

Cuando se resuelve un problema utilizando componentes simétricas, se acostumbra designar las tres
fases del sistema como a, b y c, de forma que la secuencia de fase de los voltajes y las corrientes en el sistema
es abc. Así, la secuencia de fase de las componentes de secuencia positiva es abc y la secuencia de fase de las
componentes de secuencia negativa es acb. Si los fasores originales de voltaje se designan como
V& a , V& b y los tres conjuntos de componentes simétricas se designan agregando un subíndice (o
V& ,c
superíndice) adicional 1 para las componentes de secuencia positiva, 2 para las de secuencia negativa y 0 para
las de secuencia cero. Una vez obtenidos los resultados en el dominio de las componentes simétricas, los
valores reales en cantidades de fase se calculan haciendo uso de una transformación inversa adecuada.

b. Relación entre voltajes (corrientes) de secuencia y de fase: La Figura 4.16 muestra los tres
sistemas equilibrados de vectores (considerándolos como tensiones) y la suma gráfica de los componentes
para obtener los fasores desbalanceados.

V Va
c1
Va Va
2 Va0
1 Vb2
Vb0
Vc
Vc
Vb 0 Vb
Vc2
1

a) b) c) d)

Figura 4.16.- Componentes de secuencia: a) positiva, b) negativa, c) cero, d) suma gráfica de ellas

Como cada uno de los vectores desequilibrados originales es igual a la suma de sus respectivos
componentes simétricos, se puede escribir:
V&  V&
a  V& a0 a1

 V& a2 V& b  V& b0

(4.28)
 V& b1  V& b2 V& c

 V& c0  V& c1  V&

c2

Si se consideran como referencia los fasores V& a1 , V& a 2 y V& a0 , respectivamente se

tiene:

V& a1  Va10º V& a 2  Va V& 

a 0
V& b1  2 0º Va0 0º V& (4.29)
b0
V 240º V&
a1 c1 V& b2  Va  Va0 0º V&
 Va1120º 2 120º c0  Va0 0º
V&  V
c2 a
2 240º

Designando como a& , al operador complejo que origina un desplazamiento de 120º, es decir:
1 1
a&  1120º   j
(4.30)
3
2 2

e introduciendo las expresiones (4.29) y (4.30) en (4.28), esta última se puede escribir como:

V& a  V& a0  V& a1  V& a 2

V& b  V& a0  a 2 V& a1  a V& a 2 (4.31)
V&  V&  a V&  a 2 V&
c a0 a1 a2

La ecuación (4.31) se puede escribir en forma matricial, tal como se muestra en la expresión (4.32):
 V& a  1 1 1  V& a0 

 &  1a2  & (4.32)
V
 b  a V a1 
 V&  1 a a 2  V& a 2 
c
   

o bien:

a,b,c 0,1,2
 V   T  V  (4.33)
   

donde [T] es la matriz de transformación que permite obtener las componentes de fase abc a partir de
las de secuencia 012, cuyo valor es:

1 1 1
T 1 a 2 a

(4.34)
1 a a2

 La matriz de transformación [T] es no singular y por lo tanto existe su inversa, de manera que es
posible obtener las componentes de secuencia 012 a partir de las de fase abc. Pre-multiplicando (4.33) por la
inversa de T, [T]-1 se obtiene:

0,1,2 a,b,c
 V   T -1  V  (4.35)
   

en que:
1 1 1
1
T 1
 1 a 2
a  (4.36)
3
1 a 2 a


y la ecuación (4.35) queda

V & a0  1   V& a 

1 1 1
V&  1 a 2  & (4.37)
a V
 a1    b & 
 V& 3 1 a 2 a  V c
 a
2
   

Las ecuaciones (4.32) y (4.37) son válidas también para las corrientes, es decir:
 &I a 
1 1 1  &I a0 &I a0  1 1 1  &Ia 
&    &  1 
I  1 a   I  1 a &
2 (4.38)
2 a &I a I
 b  1 a   a1   a1  &b  

&I  a 2  &I a 
&I 3 1 a 2 a  I c
 c 2  a2
       

De la segunda ecuación de (4.38) se puede concluir que si en un sistema trifásico no existen

conductor neutro o conexiones a tierra, o si el sistema está balanceado, la corriente de secuencia cero es nula

c. Potencia en función de los componentes simétricas: Si se conocen las componentes de secuencia

de la corriente y tensión, se puede calcular directamente la potencia suministrada en un circuito trifásico a
partir de dichas componentes.

La potencia compleja total transmitida en un circuito trifásico por 3 líneas; a, b y c viene dada por:

S&  P  jQ  a I *  I *b  V&Icc (4.39)

a b
V& V& *

en que
V& a , V& b son las tensiones respecto al neutro en los terminales e &I , &I , las corrientes
a b

y V& c
&I
c
que entran al circuito por las tres líneas. Puede existir o no neutro.

Matricialmente se
tiene:
t 
I*a  a,b,c a,b,c
* 
S&
V& a V&
b
V
&
c
 I    V
b  
I

(4.40)
  I*     
 c

Introduciendo (4.32) y (4.38) en (4.40) y haciendo las operaciones correspondientes se obtiene:

S&  3
I *a0  V& *  V& a 2 *  (4.41)
 &
V a0 a1 I a1 Ia
2
Es decir, esta transformación no es invariante a la potencia compleja.

4.5.2. Circuitos equivalentes de secuencia de los elementos componentes de un SEP

La aplicación del método de las componentes simétricas al cálculo de cortocircuitos asimétricos

significa que cada componente del SEP se representa por tres circuitos equivalentes monofásicos,
correspondiendo cada uno a una determinada secuencia. En cada uno de estos circuitos equivalentes las
variables tensiones y corrientes corresponden a una misma secuencia y las impedancias asociadas a los
elementos reciben el nombre de impedancia a la secuencia que corresponde. Veremos a continuación, los
circuitos equivalentes de secuencia de los elementos componentes del sistema.

a. Líneas: Las líneas se representan tal como se muestra en los circuitos de la Figura 4.17.
Z1 Z2 Z0

a) b) c)
Figura 4.17.- Circuitos equivalentes de secuencia: a) Positiva; b) Negativa y c) cero, de líneas de transmisión

Generalmente: Z1 = Z2  Z0; ya que en secuencia cero es necesario considerar tanto el efecto del
retorno por tierra, como el de los conductores de guardia, en caso que ellos existan, debido a que la corriente
se reparte por ambos caminos

b. Generadores: Un generador de rotor cilíndrico operando en condiciones de carga balanceada y

despreciando el efecto de la resistencia de sus enrollados, se puede representar según el circuito equivalente
que se muestra en la Figura 4.18. Directamente de esta figura se puede escribir:

V& a  E& a  jX s &I a

V& b  E& b  (4.42)
jX &I s b

V& c  E& c  jX s
&I
c

o bien:

a,b,c a,b,c a,b,c a,b,c a,b,c jXs 0 0


VE ZI con  Z   0 jXs 0 (4.43)

   
       0 jXs 
0
Ea jXs
Ia
- +

Eb jXs Ib
- +

Ec Va
jXs Ic
-+
Vb

Vc
Zn In

Figura 4.18.- Generador de rotor cilíndrico operando en condiciones balanceadas

El análisis de un generador operando en régimen permanente y con carga desbalanceada, es mucho

más complicado que el caso anterior; sin embargo, sus ecuaciones de comportamiento tienen la misma forma,
variando sólo en la matriz de impedancia. Se puede demostrar que en este caso:
  Zs Zm1 Z m2 
a,b.c  
 Z    Z m2 Zs Zm1
 (4.44)
  Zm1 Z m2 Zs 

donde: Zs, Zm1 y Zm2 son funciones complicadas de las inductancias propias y mutuas de todos los
enrollados de la máquina. Esta matriz se puede transformar a una matriz de impedancia de secuencia,
utilizando la siguiente expresión:

0,1,2 1 
a,b,c 
(4.45)
  T
 TZ   Z 
   

Introduciendo (4.44) en (4.45) y haciendo las operaciones respectivas se obtiene:

0,1,2 Z0
 0 0 
 Z  0 Z1 0 (4.46)
 
  
 0 Z 2 
0

donde:

Z0  Zs  Zm1  Zm2
Z1  Zs  a 2 Zm1  a Z m2 (4.47)
Z 2  Zs  a Zm1  a 2 Z m2

Se observa que Z0, Z1 y Z2 son distintas y que no existe impedancia mutua entre las redes de
secuencia, ya que los elementos fuera de la diagonal de la matriz de impedancia de secuencia son todos nulos.
Esto significa, que las redes de secuencia resultan desacopladas.

La expresión (4.43) en componentes simétricas queda:

0,1,2  0,1,2 0,1,2 0,1,2

VEZI (4.48)
       

Es decir:

 V& a0  E& a0 Z0 0 0  &I a0 

  
 &   0 Z  &  (4.49)
V  0 I
E& a1   a1   1   a1 
V&  E&  0 Z 2 &I a 2 
a2  a 0 
2
     

o bien:

V& a0  E& a0  Z0
&I V&  E&
a0 a1 a1 (4.50)
 Z1 & I a1 V& a2 
E& a 2  Z 2 &I a2

Estas ecuaciones permiten representar el generador mediante tres circuitos monofásicos

independientes (uno para cada secuencia). La Figura 4.19 muestra los circuitos equivalentes de secuencia de
un generador síncrono, donde se ha considerado que, como ocurre normalmente, las tensiones generadas son
equilibradas y por lo tanto:
 E& a0  E& a 2  0 (4.51)
Z0 Ia0 Z1 Ia1 Z2 Ia2
+ + +
+
Va0 Va2
Ea1 Va1
-
- - -
a) b) c)

Figura 4.19.- Circuitos equivalentes de secuencia de una generador síncrono: a) Secuencia cero; b) secuencia
positiva; c) secuencia negativa
Algunas observaciones:

 Las impedancias de secuencia Z 0, Z1, Z2, de un generador se pueden calcular en forma analítica a
partir de los parámetros fundamentales de la máquina; sin embargo, usualmente se determinan en
forma experimental.
 La corriente de secuencia cero existirá sólo si el generador está puesto a tierra, directamente o a
través de una impedancia.
 La barra de referencia para las redes de secuencia positiva y negativa es el neutro del generador ya
que por la impedancia Zn sólo circula corriente de secuencia cero (Figura 4.18). La barra de
referencia para la red de secuencia cero es la tierra del generador.
 La corriente In es igual a 3 Ia0, por lo tanto, la caída de tensión de secuencia cero entre una fase
cualquiera y tierra es IaoZ0 3 Ia0 Zn. Como la malla de secuencia cero es un circuito monofásico
(por fase) por el que circula sólo la corriente de secuencia cero de una fase, debe tener una
impedancia total de 3 Zn + Z0. De lo anterior se puede inferir que habrá distintos tipos de mallas de
secuencia cero, dependiendo de la conexión del generador, algunas de las cuales se muestran en la
Figura 4.20.

Z0 Ia0
Z0 Ia0 Z0 Ia0
+ + +
Va0 Va 3Zn Va0
- 0 -
a)
b) - c)

Figura 4.20.- Circuitos equivalentes de secuencia cero de un generador síncrono en conexión: a) Estrella
aislada de tierra; b) Estrella a tierra directa; c) Estrella a tierra a través de una impedancia Z n

c. Transformadores: Consideremos el circuito equivalente de Thevenin de un transformador trifásico

de dos enrollados operando en condiciones balanceadas, que se muestra en la Figura 4.21.
Zeq2 Ia
a
Va' + Va
- n
- -
+ Zeq 2 I b Vb
+
b
Vc' Vb'
Zeq2 Ic Vc
c

Figura 4.21.- Circuito equivalente de Thevenin de un transformador trifásico de dos enrollados

 Zeq2 es la impedancia equivalente referida al secundario.

De la Figura 4.21 se puede escribir:

V& ' a V& a  Z&

&I V& ' V&
eq2 a b b (4.52)
 Z& eq2
&I V& '
b c

V& c  Z& eq2 &I c

o bien:

a',b',c'  a,b,c a,b,c a,b,c

VVZI (4.53)
       

a,b,c
en que,  Z  es la matriz de impedancia del transformador, definida como:
 

 0 
a,b,c Zeq2 0
 Z  0 Zeq2 0  (4.54)

  
 0 0 Zeq2 

Cuando un transformador trifásico de dos enrollados opera con carga desbalanceada, no es posible
emplear directamente la ecuación (4.45) para llevar la impedancia en componentes de fase, a componentes
simétricas y de allí deducir los circuitos equivalentes de secuencia. Se puede verificar sin embargo, que los
circuitos equivalentes de secuencia positiva y negativa son iguales entre sí y corresponden a los ya
estudiados. En cambio, el circuito equivalente de secuencia cero depende del tipo de conexión de los
enrollados primario y secundario y de la existencia o no, de neutros conectados a tierra en los enrollados. La
impedancia de secuencia cero puede tener valores totalmente diferentes según sean los terminales del
transformador que se consideren.

La Figura 4.22, ilustra el diagrama general necesario para determinar experimentalmente la

impedancia de secuencia cero de un transformador trifásico de dos enrollados. Las líneas de segmentos
corresponden al caso que existan neutros conectados a tierra.

A I0 a
I0 Transformador Trifásico b
B
I I0 c
C
+V
Figura 4.22.- Diagrama
- general para determinar la impedancia de secuencia cero

A partir de esta figura, la impedancia de secuencia cero es:

V V
Z  (4.55)
0
I0 I/3
 Tal como se señaló anteriormente, el valor de Z 0 puede ser diferente según se hagan las mediciones
en el primario o secundario. La Figura 4.23 incluye una tabla con las mallas de secuencia cero asociadas a
diferentes conexiones de transformadores trifásicos de dos y tres enrollados.

P Z0 S P Z0 S
P S P S

P Z0 S P Z0 S
P S P S

P Z0 S P Z0 S
P S P S

a)

T T
P S P S

P S P S

T T

T T
P S P S

P S P S

T T

b)

Figura 4.23.- Redes de secuencia cero de transformadores de a) dos y b) tres enrollados

d. Redes de secuencia: Un SEP balanceado se puede representar por tres redes de secuencia
independientes entre si (sin acoplamientos); una red de secuencia positiva, una red de secuencia negativa y
una red de secuencia cero. Cada red de secuencia representa una fase del SEP y todas las impedancias
corresponden a una determinada secuencia.

La red de secuencia positiva es la única que normalmente contendrá fuentes de fem, según lo
expuesto. Por otra parte, teniendo presente las aproximaciones usuales que se realizan en los cálculos de
cortocircuito; es decir, que las fem de todos los generadores se consideran iguales en módulo y ángulo de fase
y que se desprecian las corrientes de prefalla, se concluye que en ausencia de cortocircuitos en el sistema no
existirán corrientes en ninguna de las redes de secuencia. Por lo tanto las redes de secuencia negativa y cero,
son totalmente pasivas antes de falla.

Para los efectos del cálculo de cortocircuitos asimétricos es necesario establecer para cada red de
secuencia, su circuito equivalente de Thevenin mirado desde el punto de falla. Supongamos que se produce
una falla en el punto F de un sistema y que los circuitos equivalentes de Thevenin corresponden a los
indicados en la Figura 4.24, donde la tensión prefalla en F es V a(0). La corriente de falla en dicho punto
tendrá en general componentes de secuencia positiva, negativa y cero. En estas condiciones circularán
corrientes en todas las redes de secuencia y aparecerán tensiones en sus terminales. Esto significa que las 3
redes deben interconectarse en una forma que dependerá del tipo particular de falla y cuyos detalles veremos
luego. La situación post-falla se puede ilustrar esquemáticamente como se muestra en la figura 4.24.

Z1
Ia1 F1
+
+
Va (0)
- Va1 Red
-
Z2 Ia2 F
2
+
Va2 de

-
Z0 Ia0
F0
+ Interconexión
Va0
-

Figura 4.24.- Mallas de secuencia y red de interconexión

Directamente de esta figura se puede escribir:

&
V& a0  0  I a0
Z0
&
V& a1  V& a (0) I (4.56)
 Z1 a1
&
V& a 2  0  I a2
Z2

o bien:

V & a0   0  Z0 0 0  &I a0 

  
V&  V& (0)  0 Z   (4.57)
 a1   a   1 0   &a1I 
 V &   0  0
a2
0 Z 2  &I a 2 
   
 Es decir:

0,1,2   0,1,2  0,1,2 0,1,2

VF   V(0)   Z   I F  (4.58)
      

4.5.3. Análisis de algunos tipos de cortocircuito

4.5.3.1. Cortocircuito monofásico a tierra a través de una impedancia de falla ZF

– Diagrama esquemático: La Figura 4.25 muestra en forma esquemática esta situación

Ia
a
Ib
b Ic
Va ZF
c Vb
Vc

Figura 4.25.- Representación esquemática de un cortocircuito monofásico

– Condiciones impuestas por la falla: A partir de la Figura 4.25 se puede escribir

&I  &I 
b c
& (4.59)
0V  a

Z F &I a

– Componentes simétricas: Las componentes simétricas de las corrientes quedan:

&I a0  1 1 1  &Ia 
&I  1 1 a 2  & (4.60)
a I
 a1  &b  

&I 3 1 a 2 a  I c
 a2
   

de donde se obtiene:

&I 1&
a0   &I a  I (4.61)
a
&I 2 3
a1

Para las componentes simétricas de los voltajes se tiene:

V & a0  1   V& a 
1 1 1
V&  1 a 2  & (4.62)
a V
 a1   b

 V& 3 1 a 2 a    V& c 
 a
2
   

y por lo tanto: o bien, a

partir de
V& a0  V& a1  V& a 2  V& a  Z F &I a (4.61)
V& a0  V& a1  V& a 2  3&I a0 Z F = 3&I a1Z F =
3&I Z
a2 F (4.63)

(4.64)
– Conexión de las mallas: A partir de (4.61) y (4.64), se deduce que las mallas de secuencia quedan
conectadas en serie, tal como se muestra en la Figura 4.26.

Z1
Ia1 F1
+
+
Va (0) Va1
-
-
Z2 Ia2 F
2
+
Va2
3 ZF
-
Z0 Ia0
F0
+
Va
0
-
Figura 4.26.- Interconexión de las mallas de secuencia para una falla monofásica

Del circuito de la Figura 4.26 se tiene:

&I V& a (0)

a0   &I a  (4.65)
&I 2
Z1  Z2  Z0  3ZF
a1

Conocidas las corrientes de secuencia, se pueden determinar las corrientes de las fases, utilizando la
primera ecuación de (4.38) y se obtiene:

&I 3V& a (0)

a  (4.66)
Z1  Z2  Z0  3ZF

&I  &I  0 (4.67)

b c

Para los voltajes de secuencia se puede escribir:

V& a0  &I
a0
Z 0
V& a1  V& a (0) &I (4.68)
 Z1 a1

V&  a2
&I
a2
Z 2

y por lo tanto los voltajes de las fases quedan:

V  3 Z F V& a (0)
(4.69)
& Z1  Z2  Z0  3ZF
a
Z0 (a 2  1)  Z2 (a 2  a)  3 a 2 ZF &
& (4.70)
Vb  Va (0)
Z1  Z2  Z0  3ZF

& Z0 (a  1)  Z2 (a  a 2 )  3 a ZF &
Vc  Va (0) (4.71)
Z1  Z2  Z0  3ZF

Si el cortocircuito es directo a tierra, basta con hacer Z F = 0 en las expresiones (4.65) a (4.71).
4.5.3.2. Cortocircuito bifásico a tierra a través de una impedancia de falla ZF

– Diagrama esquemático: La Figura 4.27 muestra esta situación

Ia
a
Ib
b Ic
c
Va

Vc ZF Vb

Figura 4.27.- Representación esquemática de un cortocircuito bifásico a tierra a través de una impedancia

– Condiciones impuestas por la falla: A partir de la Figura 4.27, se puede escribir:

&I  0
a

V& b  V& c  Z F (&I b  (4.72)

&I )
c

– Ecuaciones en componentes de secuencia: Las componentes simétricas de las corrientes y de los

voltajes quedan:
&I  &I  &I  0
a0 a1 a2
& &
V a0  3 Z F I a0  (4.73)
V& a1

– Conexión de las mallas: A partir de las ecuaciones (4.73), las mallas de secuencia quedan conectadas
en paralelo, tal como se muestra en la Figura 4.28.

Z1
Ia1 F1
+ +
Va (0) Va1
-
-
Z2 Ia2 F
2
+
Va2
3 ZF
-
Z0 Ia0
F0
+
Va
0
-

Figura 4.28.- Conexión de las mallas de secuencia para un cortocircuito bifásico a tierra

Haciendo ZF=0 y ZF= (infinito), en el circuito de la Figura 4.28, es posible modelar el cortocircuito
bifásico a tierra directo y el cortocircuito bifásico aislado de tierra, respectivamente.
Ejemplo 4.3. En el sistema de la Figura 4.29, ocurre una falla bifásica a tierra en la barra 1, a través de una
impedancia de falla ZF=j0,05 (pu). Con las consideraciones usuales del cálculo de cortocircuitos y
considerando que todos los datos en % están en base 100 MVA, determinar:
a. Las corrientes de línea en Amperes, en el punto de falla
b. Los voltajes entre líneas en kV, en el punto de falla
c. Las corrientes de línea en Amperes en bornes de G 1 y de G2
Línea 2-3

G1 2 X1=X2=40%  G2
1 3 4
X0=120%

50 MVA 15 kV
15 kV 75 MVA 75 MVA 50 MVA
X1 =X 15/115 kV 115/15 kV X1=X2=20%
2=20%
X0 =5% X1=X 2=X X1=X2=X0=10% X0=5%
0=10% Y1
Y1 Figura 4.29

Solución
:

Mallas de secuencia: Figuras 4.30 a 4.32

- Secuencia positiva - Secuencia negativa

j0,1 j0,4 j0,1 j0,1 j0,4 j0,1

I'a1 I'a2 I a2 I"a2

I a1 I"a1
j0,2 j0,2
Z1 j0,2 Z2 j0,2
+ +
1 0º 1 0º

Figura 4.31
Figura 4.30
- Interconexión de las mallas
- Secuencia cero
Z1
j0,1 j1,2 j0,1
+ Ia1
+
Ia0 I"a0 Va1
I'a0 Va (0)

j0,05 Z0 j0,05 Z2

Ia2
+
Va2
3ZF
Figura 4.32 Z0

De la Figura 4.33 y a partir de las condiciones de Ia0

+
prefalla se tiene que: Va(0)=10º Va0

Las impedancias Z1, Z2 y Z0 corresponden a las

impedancias de Thevenin de las respectivas mallas de Figura 4.33
secuencia. Sus valores son: Z1=Z2=j0,16 y Z0=j0,05.
a) Corrientes de línea en el punto de falla

a1) Corrientes de secuencia: A partir del circuito de la Figura 4.33 se obtiene:

&I
a1 10º
  4,0179  90º (pu)
j0,16  j0,2
j0,16
 j(0,16  0,2)
&I j0,2
a 2   &I a1  2,232290º
j(0,2 

&I 0,16)
a0  1,785790º
j0,16
 &I a1
j(0,2 

0,16)

a2) Corrientes de línea en (pu): Utilizando la matriz de transformación [T]

 &Ia  1 1 1  1,785790º  &Ia   0 
   
 2    &  
&I  1 a a 4,0179  90º  I  6,0392153,67º
 b   2   & b   
&I
 c 1 a a  2,232290º
  Ic  6,039226,33º 
 
   

a3) Corrientes de línea en Amperes (sector 1):

100 106 & & &

IB1   3.849 A  Ia  0; I b  23.244,88153,67º A; Ic  23.244,8826,33º A
3 15 103

b) Voltajes entre líneas en el punto de falla

b1) Voltajes de secuencia: A partir de la Figura 4.33 se tiene que:

V& a1  10º j0,16  4,0179  90º  0,35710º  V& a 2

V& a0   j0,05 1,785790  0,089290º

b2) Tensiones fase neutro: Utilizando la matriz de transformación [T]

 V& a  1 1
1  0,089290º  V&  0,80350º 
 &      a   
V  1 a2 a 0,35710º  V&  0,2678  180º
 b      b  
 V& 1 a a 2  0,35710º  V & 0,2678  180º 
   
c c
   

b3) Voltajes entre líneas en (pu)

V& ab  V& a  V& b  0,80350º0,2678  180º 

1,07130º V&  V&  V&  0
bc b c

V& ca  V& c  V& a  0,2678  180º0,80350º  1,0713180º

b4) Voltajes entre líneas en kV: Considerando que en ese sector (sector 1), VB1=15 kV, el voltaje base por
fase será de 15 / 3 kV, por lo tanto:

V& ab  1,07130º 3  kV; V& bc  V& ca  9,2777180º kV

*15 / 9,27770º 0;
c) Corrientes de línea en bornes de G1 y G2

c1) Corrientes de secuencia en bornes de G 1: De los circuitos de las Figura 4.30 a 4.32, se obtienen:

&I j0,8
a1G 
 4,0179  90º  3,2143  90º (pu)
1
&I ' a1  j(0,8  0,2)

&I  j0,8
a  2,232290º  1,785890º (pu)
2G1 &I ' a
 j(0,8  0,2)
2
&I
a0G1  a0  1,785790º (pu)
&I '

c2) Corrientes de línea en bornes de G1: Aplicando la matriz de transformación se tiene

 &IaG1  1 1  &I aG1 0,357290
1   1,785790º 

  2    
&  1a a& 3,2143  90º   5150º
I I 
 bG1  2    bG1  
 &I 1 a a  1,785890º  &I  530º 
 cG1  cG1
   

c3) Corrientes de línea en Amperes en bornes de G1 (sector 1):

IB1  100 106

3 15 103
 3.849 A  &I aG1  &I
bG1  19.245150º
&I
cG1  19.24530º A
1.374,8690º ; A;

c4) Corrientes de secuencia en bornes de G 2: De los circuitos de la Figura 4.30 a 4.32, se obtienen:

&I j0,2
a1G
 &I "   60º  4,0179  90º   60º  0,8036  150º (pu)
2
 a1 j(0,8  0,2)

&I j0,2
a 2G  60º  2,232290º 60º  0,4464150º (pu)
2 a
&I " 2  j(0,8  0,2)
&I  &I "'  0
a0G 2 a0

c5) Corrientes de línea en bornes de G2: Aplicando la matriz de transformación se tiene

 &IaG2  1 1  &I aG2
1  0 1,09717 170,63º

     
&  1a2 a& 0,8036 150º    0,357290º
I I 
 &IbG2       &IbG2   
 cG2 1 a a 2  0,4464150º
  cG2  1,09717  9,37º 

   

c6) Corrientes de línea en Amperes en bornes de G2 (sector 3)

100 106 & & &

IB3   3.849 A  IaG2  4.223  170,63º ; I bG2  1.374,8690º A; IcG2  4.223  9,37º A

4.5.3.3. Cortocircuito entre dos fases a través de una impedancia de falla

Este tipo de falla es muy poco frecuente y produce sobrecorrientes inferiores a las de los otros tipos
de corcircuitos, por lo que sólo se calcula en casos excepcionales. Sin embargo, su análisis resulta interesante
ya que es idéntico al de una carga conectada entre dos fases (carga bifásica).
– Diagrama esquemático: La Figura 4.34 muestra esta situación
Ia
a
Ib
b Ic
ZFVa
c
Vc
Vb

Figura 4.34.- Representación esquemática de un cortocircuito entre dos fases

– Condiciones impuestas por la falla: A partir de la Figura 4.34, se puede escribir:

&I  0
a &I  &I 
b c (4.74)
0
V& b  V& c  Z F &I b

– Ecuaciones en componentes de secuencia: Las componentes simétricas de las corrientes y de los

voltajes quedan:
&I  0
a0
&I  &I  0 (4.75)
a1 a2

V a1  V& a 2  Z F
&
&I
a1

– Conexión de las mallas: Las ecuaciones (4.75) nos indican que la malla de secuencia cero no
interviene y que las mallas de secuencia positiva y negativa quedan conectadas en paralelo, según se
muestra en la Figura 4.35
Z1 Ia1 F
1
+ +
Va1
Va (0) ZF
- -
Z2 Ia2 F
2
+
Va2
-

Figura 4.35.- Conexión de las mallas de secuencia para un cortocircuito entre dos fases

4.5.3.4. Observaciones finales respecto a los cortocircuitos asimétricos

 Las corrientes y tensiones de secuencia calculadas corresponden sólo al punto de falla, no a otro. Para
calcular las corrientes y tensiones en otros puntos distintos, se debe resolver los respectivos circuitos.
 Los cortocircuitos asimétricos pueden producir sobre tensiones en las fases no falladas, las que
dependen de la relación entre X0 y X1 y de la existencia o no de impedancia de falla.
 Para limitar los valores de corriente de cortocircuito de las fallas a tierra se utilizan impedancias entre
el neutro y tierra, las que pueden ser de tipo resistiva pura o reactiva pura.
 En aquellas partes del sistema que estén separadas del punto de falla por transformadores Y/  o
viceversa, se deben considerar los desfases de las componentes simétricas de las corrientes y de los
voltajes, introducidos por la conexión del transformador.
4.6. Fases abiertas

4.6.1. Introducción

Las fallas de conductor abierto o las fases abiertas, son los defectos producidos por la interrupción de
una o más fases, sin contacto simultáneo con otras fases o tierra. Aunque no producen corrientes elevadas,
provocan la circulación de corrientes de secuencia (en especial negativa) que son peligrosas para los equipos
por el fuerte calentamiento que pueden originar.

A primera vista, el cálculo empleando componentes simétricos se ve complicado por el hecho de que
las fallas implican una asimetría en las impedancias del sistema, lo que haría necesario considerar los
acoplamientos entre mallas de secuencia.

El problema se resuelve aplicando a las mallas de secuencia, supuestas independientes y sin

impedancias mutuas, las condiciones eléctricas impuestas por la falla. Como las condiciones impuestas a las
tres mallas están relacionadas entre sí, ello equivale a interconectar las mallas en el punto de falla, en una
forma fijada por el tipo de falla.

El fenómeno que sigue a la aparición de la falla es transiente, donde las corrientes máximas se
producen en el instante inicial. Normalmente interesa determinar lo que ocurre al cabo de algunos ciclos de
iniciada la falla (operación de las protecciones, apertura de interruptores, etc.), por lo que en secuencia
positiva, los generadores se representan por la fem E' y la reactancia transitoria X1' . Sólo cuando interesa
verificar los esfuerzos electrodinámicos de los equipos o al especificar interruptores, se considera E'' tras X "1.

Una dificultad preliminar en el estudio de este tipo de fallas será entonces la de calcular las fem E' (o
E"), a partir de las condiciones de operación existentes antes de la falla.

Dada la simetría longitudinal de estas fallas, se acostumbra usar como variables de cálculo, las caídas
longitudinales de tensión V& a , V& b y entre los bornes P y Q de la zona en falla y las corrientes en las
V& c
fases: &I a , &I b e &I c tal como se indica en la Figura 4.36.

Ia
a P+ V -
Q a'
a

Ib b + Vb - b'

Ic c + Vc - c'

Figura 4.36.- Modelación de fallas tipo fases abiertas

Para evitar la aparición de razones de transformación no reales (a, a 2, etc.) en las ecuaciones de
conexión, es preciso mantener en el análisis una simetría respecto a la fase de referencia a , por lo que la falla
monofásica se supone en la fase a y la bifásica en las fases b y c.

4.6.2. Una fase abierta

Esta situación se presenta, por ejemplo, cuando se emplean elementos de apertura que controlan
individualmente cada una de las fases (fusibles o interruptores de accionamiento monopolar). A veces ocurre
también al cortarse un conductor y quedar suspendido de tal forma de no hacer contacto con otra fase o tierra.
– Diagrama esquemático: La Figura 4.37 muestra esta situación

Ia a P+ V -
Q a'
a

Ib b + Vb - b'

Ic c + Vc - c'

Figura 4.37.- Representación esquemática de una falla de una fase abierta

– Condiciones impuestas por la falla: A partir de la Figura 4.37, se puede escribir:

&I  0
a
(4.76)
V&  V&
b c 0

Es decir:

&I  &I '

a a 0
&I  &I ' &I 
b b c c (4.77)
&I '
c
V& b  b   V&  0
'

V&
' V& c

– Ecuaciones en componentes de secuencia: Las componentes simétricas de las corrientes y de las

caídas de voltajes quedan:
&I  &I  &I  0
a0 a1 a2
1 (4.78)
V& a0 V&
 V&  V& 
a
a1 a2 3

– Conexión de las mallas: A partir de (4.78), se puede concluir que las mallas de secuencia quedan
conectadas en paralelo entre los punto P y Q, tal como se indica en la Figura 4.38.

Z1 Ia1 P1 - Q1 Z'1
+  Va
+ 1 +
EA'
- E'B
-
Z2 Ia2
P2 - Q2 Z'2
+ Va2

Z0 Ia0 Q Z'0
P0 0
+  Va0 -
Figura 4.38.- Conexión de las mallas de secuencia para una falla de una fase abierta
 Puesto que las mallas de secuencia negativa y cero son pasivas, su efecto es el de intercalar una
impedancia:
Z2pq Z0pq
Z  Z 2  Z' 2 Z 0  Z'0  (4.79)

20pq
Z2pq  Z0pq Z 2  Z'2    Z'0 
Z 0
entre los bornes P y Q de la malla de secuencia positiva. Por lo tanto, aumenta la impedancia serie de
la malla de secuencia positiva, lo que significa que se reduce la corriente y en consecuencia, la potencia
activa transmitida. En algunos casos particulares y, debido a las conexiones de los transformadores vecinos a
P y Q, puede resultar que Z0pq=, en cuyo caso aumenta aún mas la impedancia serie agregada a la malla de
secuencia positiva, haciendo que la disminución de potencia transmitida sea mayor.

Es conveniente indicar que Z 0pq y Z2pq son las impedancias equivalentes vistas en esas mallas, desde
los bornes P y Q.

Ejemplo 4.4. En el sistema de la Figura 4.39, se abre la fase “a” en la barra 3 cuando el motor M está
recibiendo el 80% de su potencia nominal, con su tensión nominal en bornes, Factor de Potencia 0,8
inductivo. Calcular la potencia recibida por el motor (kVA) y las corrientes en los neutros de los
transformadores en estas condiciones. Datos en % en base común 1.250 kVA.
Línea 2-3

G 2 X1=X2=16% 
1 3 4
X0=20%
M
1.250 kVA 600 V
600 V 1.250 kVA 1.250 kVA 1.250 kVA
X 1' = 600/4.160 V 4.160/600 V X1=X2=10%
X1=X2=10%
X0=4% X1 =X2=X0=5% X1=X2 =X0 X0=4%
Y1 =5%
Figura 4.39 Y11
Solución
:

a) Condiciones de prefalla: El circuito equivalente por fase se muestra en la Figura 4.40

Del Circuito de la Figura 4.14 se tiene: j0,1 j0,05 j0,16 j0,05 j0,1
*
 S  0,8  cos1 0,8 I +
&I  + +
    0,8  36,87º
 MV 10º E'G VM EM
E&  10º j0,36  0,8  36,87º 
'

1,195211,11º
G
E& M  10º j0,1 0,8  36,87º  0,9541 
Figura 4.40
3,85º

b) Condiciones de falla: Como se abre una sola fase las mallas de secuencia quedan en paralelo y se
muestran en la Figura 4.41. A partir de este circuito se tiene:

X 0,46  0,3
 0,46  0,3  0,1816
20
1,195211,11º0,9541    0,5736  36,87º
&I 3,85º
 &I
a1M a1M
j(0,31  0,1816 
0,15)
V&  &I a1  X 20  0,5736  36,87º  0,181690º j0,31 + j0,05 j0,1
V
+ I a1 + I a1M +
E'G
V&  0,104253,13º Va1M E
M

&I V& 0,104253,13º

   0,2265143,13º j0,31 j0,05
a2M
j0,1
X2 j0,46
+I
I a2 a2M
&I Va2M
a0M 0  V& a0M  0

j0,25 j0,05
V& a1M  j0,1 0,5736  36,87º0,9541  I a0M
+
j0,04 I a0
Va0M j0,04
3,85º V& a1M  0,9865  1,06º

Figura 4.41
V& a 2M  j0,1  0,2265143,13º  0,02265 
126,87º

c) Potencia que llega al motor en estas condiciones

S& 3f  3(V&   V&   V&  &I * )

a0 a1 a2
&I* &I*
a a1 a
0 2

S& 3f  3(0  0,9865  1,06º  0,573636,87º0,02265  126,87º  0,2265  143,13º )

S& 3f  1,706736,23º (pu)

S& 3f  1,706736,23º  1.250 

kVA
711,12536,23º 3

S& 3f  (573,63  j420,29) kVA

d) Corrientes en los neutros de los transformadores en pu

&I 0,104253,13º
a0
 &   0,3473  36,87º
V (pu) j0,3
j0,
3

&I &
NT  3  I a0  3 * 0,3473  36,87º  1,0419  36,87º (pu)

d) Corrientes en los neutros de los transformadores en Amperes

IB   173,48 A
1.250 103
3  4.160
 &I
NT 
1,0419

36,87º
173,4
8
180,75

36,87º
A
 4.6.3. Dos fases abiertas

Esta situación se presenta en los mismos casos que una fase abierta, pero con una frecuencia menor.

– Diagrama esquemático: La Figura 4.42 muestra esta situación

Ia a P+ V -
Q a'
a

Ib b + Vb - b'

Ic c + Vc - c'

Figura 4.42.- Representación esquemática de una falla de dos fases abiertas

– Condiciones impuestas por la falla: A partir de la Figura 4.42, se puede escribir:

&I  &I 
b c b  &I ' 
&I ' 0
(4.80)
&I  &I '
a a

V& a  0  V& a 
a V
&' 0

– Ecuaciones en componentes de secuencia: Las componentes simétricas de las corrientes y de las

caídas de voltajes quedan:
&I  &I  &I
a0 a1 a2

V &  V&  V& (4.81)

a0 a1 a
2 0

– Conexión de las mallas: A partir de (4.81), se concluye que las mallas de secuencia quedan
conectadas en serie tal como se indica en la Figura 4.43.

EA' I
Z 1a1
- + P1
+
E'B Va1
- + Z'1 - Q1
P2
+
Z2 Ia2 Va2
-Q2
P0
Z'2 +
 Va0
- Q0
Z0 Ia0

Z'0

Figura 4.43.- Conexión de las mallas de secuencia para una falla de una fase abierta
 De acuerdo con el circuito de la Figura 4.43, este tipo de falla equivale a intercalar una impedancia:

Z20pq  Z2pq
 Z 0pq  Z 2  Z' 2   Z 0  Z' 0  (4.82)

entre estos bornes P1 y Q1 de la malla de secuencia positiva. Con ello se reduce la potencia activa
transmitida en el sistema, en una cantidad mayor que para el caso de una fase abierta, ya que la impedancia es
mas alta. Nótese que la transmisión se interrumpe totalmente si Z 0pq=, es decir, si el sistema no está puesto a
tierra.

4.6.4. Impedancias serie desequilibradas

Un efecto similar, aunque menos grave que el de una o dos fases abiertas, produce la conexión de una
impedancia anormal en una de las fases. Es una situación que se presenta, por ejemplo, en el caso de
reemplazar temporalmente una unidad monofásica defectuosa en un banco de transformadores, por otra de
características diferentes, donde dos de las fases tendrán el mismo valor en su impedancia serie, el que será
distinto al de la tercera. Otra situación de interés práctico se presenta cuando, debido a un cortocircuito
monofásico a tierra en una línea trifásica, se desconecta la fase fallada por acción de los interruptores
(monopolares) que protegen el tramo, que corresponde al caso de una fase abierta en dos puntos.

– Diagrama esquemático: La Figura 4.44 muestra esta situación

Ia a P+ Va - Q a'
ZA
Vb
Ib + - b'
b ZB
Vc
Ic ZB
c + - c'

Figura 4.44.- Impedancias series desequilibradas

– Condiciones impuestas por la falla: A partir de la Figura 4.44, se puede escribir:

V& a = Z A &I a
V& b  Z B &I b (4.83)
V&  Z &I
c B c

– Ecuaciones en componentes de secuencia: Las componentes simétricas de las corrientes y de las

caídas de voltajes quedan:

V& a0 1
 ZB &I  (Z  Z )
B 
&I  &I )
a a2
(&I a0
A
a0 3 1

V& a1 1
 ZB &I  (Z  Z )
B 
&I  &I
a a2 (4.84)
(&I a0
A
a1 3 1 )

V& a 2 1
 ZB &I  (Z  Z )
B  &I  &I )
a a2
a2 3 (&I a0 1
A

– Conexión de las mallas: Considerando las ecuaciones (4.84), las mallas de secuencia quedan
conectadas en paralelo, tal como se indica en la Figura 4.45.
 1/3(ZA - ZB)

Z1 Ia1 P1 ZB Q1 Z'1
-
+ Va1 +
EA'- +
E' B
-
Z2 Ia2 P2 ZB Q2 Z'2
-
+ Va2

Z0 Ia0 P0 ZB Q0 Z'0
+  Va0 -

Figura 4.45.- Conexión de las mallas de secuencia cuando se tiene impedancias serie desequilibrada

Para la transferencia de potencia activa, la conexión de las mallas de secuencia en esta forma,
equivale a intercalar en la malla de secuencia positiva, la combinación de impedancias Z B en serie con el
paralelo de 1/3(ZA - ZB) con (Z2pq + ZB) y con (Z0pq + ZB).

Si ZA= y ZB=0, se tiene el caso de una fase abierta en un punto.

Si ZA=0 y ZB=, se obtiene el caso de dos fases abiertas, pero para llegar a las relaciones ya vistas
hay que calcular primero el equivalente de las impedancias en paralelo, antes de hacer tender Z B a .
Si ZA= y ZB corresponde a las respectivas impedancias de secuencia del tramo, se tiene el caso de
una fase abierta en dos puntos.

Problemas propuestos

4.1. Un generador de 65 MVA, 15,5 kV, conectado en estrella alimenta el primario de un transformador de 65
MVA, 15,5/120 kV, conexión /. La reactancia subtransiente del generador es de 0,12 (pu) y la reactancia
del transformador es 0,1 (pu), ambas en base 65 MVA.
a. Estando el generador en vacío se produce un cortocircuito trifásico simétrico en los terminales del
secundario del transformador. Suponiendo que la tensión interna del generador es 1,0 (pu), determinar la
corriente de falla subtransiente y el máximo valor posible de la componente de corriente continua.
b. Se conecta una carga impedancia trifásica balanceada en los terminales del secundario del transformador,
de valor (0,8+j0,6) (pu) en base 65 MVA, cuando la tensión en bornes del generador es la nominal.
Posteriormente se produce un cortocircuito trifásico simétrico en los terminales de la carga. Determinar
la corriente subtransiente de falla en el generador incluyendo la corriente de prefalla.

4.2. En el sistema de la Figura 4.46, todos los 3 220 kV 220 kV 2

datos en por unidad, están en base común. j0,05 j0,04
1 j0,03
Determinar las corrientes en cada una de las j0,05 j0,08
líneas y en los generadores, cuando ocurre un G3 G2
G1
cortocircuito trifásico simétrico a través de un 13,2 kV
impedancia ZF=j0,04 (pu) en la barra 1, 13,8 kV j0,05
j0,05 220 kV j0,02 13,2 kV
estando el sistema en vacío, utilizando: j0,03
a. El método tradicional (Thevenin) Figura 4.46
b. El método general (Matricial)
4.3. Repetir el Problema 4.2 considerando un cortocircuito trifásico simétrico:
a. En la barra 2.
b. En la barra 3.

4.4. En el sistema de la Figura G1 1 4 3 G

4.47, todos los parámetros están 3 X=0,05 X=0,05
en tanto por unidad, en una base L1 L3
común. Para un cortocircuito X=0,3 X=0,15
X=0,15 L2 X=0,05 X=0,2
trifásico directo en la barra 3 y
utilizando el método general,
X=0,1
(matricial); determinar los voltajes
en todas las barras, las corrientes 2
en las líneas y las corrientes en los
generadores. Suponer que el G
sistema está en vacío antes de 2 X=0,2
producirse la falla.
Figura 4.47
4.5. En el sistema de la Figura 4.48 ocurre 
un cortocircuito monofásico a tierra en la Sistema 220 kV 13,2 kV G
barra de 13,2 kV. Con las aproximaciones
usuales del cálculo de cortocircuitos,
determinar:
100 MVA
SCC3 =800 MVA 100 MVA X80 MVA
1=X2=40%
X1=X2=X0=10% X0=10%
SCC1 =500 MVA3  Y1
a. &I ; &I e &I (A) en el punto de falla.
a b c X1=X2 = X0
b. &I ; &I e &I (A) en barra de 220 kV
a b c Datos en % en base 100 MVA
c. &I ; &I e &I (A) en bornes de G.
a b c Figura 4.48
d. V& a ; V& b y V& c (kV) en bornes de G.

4.6. Repetir el Problema 4.5 considerando un cortocircuito bifásico a tierra en la barra de 13,2 kV

4.7. Repetir el Problema 4.5 considerando un cortocircuito monofásico a tierra en la barra de 220 kV

4.8. Repetir el Problema 4.5 considerando un cortocircuito bifásico a tierra en la barra de 220 kV

4.9. En el sistema de la Figura 4.49, ocurre una falla monofásica a tierra en la barra 2. Con las
consideraciones usuales del cálculo de cortocircuitos y considerando que todos los datos en % están en base
propia, determinar:
a. La corriente de falla en Amperes
b. Los voltajes entre líneas en kV, en el punto de falla
c. Las corrientes en Amperes en bornes de G 1 y G2
d. La corriente en Amperes, en los neutros de los transformadores y de los generadores
e. Las corrientes de secuencia cero en Amperes, en el interior de la delta de los transformadores

 Línea 2-3
G1 1 2
X1=X2=40%  G2
3 4
X0=120%
50 MVA
15 kV
15 kV 75 MVA 75 MVA 50 MVA
X1 15/115 kV 115/15 kV X1=X2=20%
=X2=20%
X0 =5% X1 X1=X2 =X0 X0=5%
=X2=X0=10% =10%
Y1 Figura 4.49 Y1
4.10. Resuelva el problema 4.9, considerando una falla bifásica a tierra.

4.11. En el sistema de la Figura 4.50, los valores en % de los parámetros están en la base propia de cada
equipo. El sistema se encuentra en vacío cuando ocurre un cortocircuito monofásico en el secundario del
transformador conectado al generador G1. Determinar:
a. La corriente de falla en Amperes y la tensión en el punto de falla en Volt.
b. El valor que debería tener una resistencia de falla que limite la corriente de falla máxima a un 40% de
su valor original.

G1 1  4 2 G2
L1 L2
Xn
20 MVA 70 MVA 20 MVA
50 MVA L3 13,2 kV
13,2 kV 13,2/60 kV 60/13,2 kV
X1 50 MVA X1 X1=X2=25%
X1=X2=X0=10%
=X2=30% DY1  60/13,2 kV
X1 =X2 =X0 =10%
=X2=X0=12% X0=20%
X0=20% Xn=10%
YY0
DY1
3
20 MVA
X1L
=X1=L=70
2
2=LOhm
3
13,2 kV
X1=X2=30%
G3 X 0 =100
X0=20%

Figura 4.50 Ohm

4.12. En el sistema de la Figura 4.51 a) se produce una falla en F, del tipo mostrado en la Figura 4.51 b). Los
parámetros están en tanto por unidad, en base común, 20 MVA. Determinar:
a. La interconexión de las mallas de secuencia que permita representar dicha falla
b. Las corrientes de falla en Amperes.
c. Los voltajes entre líneas en el punto de falla en kV.

 a
G F
b
20
11MVA
kV c
P T S MVA
20/10/15
X1 = 15%
Xp1 = Xp2 =Xp0 = 4% ZF
X2 = 10%
X0 =3% X s1 = Xs2 =Xs0 = 6%
X t1 = Xt2 =Xt0 = 9%
Figura 4.51 a)
Figura 4.51 b)

4.13. Un generador alimenta un motor a través de un transformador conectado en Y/. El generador está
conectado al lado en estrella del transformador. Se produce una falla entre los terminales del motor y el
transformador. Las componentes simétricas de la corriente que circulan desde el motor hacia la falla: son &I a1
 0,8  j2,6; &I a2  j2,0 e &I   j3,0 . Desde el transformador hacia la falla, las corrientes
a0
son: &I a1  0,8  j0,4; &I a 2   j1,0 e &I a0  0 . Todos los valores están en por unidad, base común.
Suponer que X1=X2 para el motor y el generador. Determinar:
a. El tipo de falla
b. La corriente previa a la falla en la fase a
c. La corriente de falla.
4.14. Una carga cuya impedancia Zc=(1,1 + j 0,6) (pu) se conecta entre las fases “b” y “c” de un generador de
20 MVA, 13,8 kV; X1=25%; X2=35%; X0=10%. En estas condiciones y considerando que todos los
parámetros están en base 20 MVA, determinar las tensiones entre líneas del generador (en volt).

4.15. Las reactancias de un generador de 100 MVA, 20 kV son: X 1=X2=20%, X0=5%. El generador está
conectado a un transformador Y1 de 100 MVA, 20/230 kV con una reactancia del 10%. El neutro del
transformador está sólidamente aterrizado. El voltaje en terminales del generador es de 20 kV cuando ocurre
una falla bifásica (entre las fases b y c) y a tierra en el lado del alto voltaje del transformador. Calcular:
a. Las corrientes y voltajes de secuencia en el punto de falla (pu)
b. Las corrientes de línea (A) y los voltajes entre líneas en el punto de falla (kV)
c. La corriente en todas las fases del generador (A)
d. La corriente en el neutro del transformador (A)
e. La corriente de secuencia cero en la delta del transformador (A)

4.16. Se tiene un transformador trifásico de

tres enrollados, con carga desequilibrada, tal
como se muestra en la Figura 4.52. Las
tensiones nominales son: 6.600/660/2.200 volt
(P/T/S). Considerando que la corriente en la
carga &I  900º A , calcular la corriente
a
en cada fase del primario y la corriente en el Figura 4.52
terciario.

4.17. En el sistema de la Figura 4.53, se abren las fases “b” y “c” en la barra 4 (bornes del motor) cuando el
motor M está recibiendo el 80% de su potencia nominal, con su tensión nominal en bornes, Factor de
Potencia 0,8 inductivo. Calcular la potencia recibida por el motor (kVA) y las corrientes en los neutros de los
transformadores en estas condiciones. Datos en % en base común 1.250 kVA.

 Línea 2-3
G1 1 2
X1=X2=16%  G2
3
4
1.250 kVA M
X0=50%
600 V
600 V 1.250 kVA 1.250 kVA 1.250 kVA
X 1' = 600/4.160 V 4.160/600 V X1=X2=10%
X1=X2=10%
X0=4% X1 X1=X2=X0 X0=4%
=X2=X0=5% =5%
Y1 Y11
Figura 4.53

4.18. Un generador, está entregando 90 MW (en sus bornes), Factor de Potencia 0,9 inductivo a un motor, a
través de un transformador cuando se abren las fase "b" y "c" en bornes del motor (secundario del
transformador). La tensión en bornes del Generador en ese instante era de 13.860 volt (entre líneas). Los
datos en tanto por unidad, con SB=100 MVA, son los siguientes:
Generador: 13,2 kV; 100 MVA, X1=X'd=0,3; X2=0,1; X0=0,04. Conexión: Estrella con neutro a tierra
Transformador: 13,2/6,6 kV, 100 MVA, X1=X2=X0=0,1. Conexión Delta-Estrella con neutro a tierra: Y1
Motor: 100 MVA; 6,6 kV; X'd=X1=X2=0,1; X0=0,04. Conexión: Estrella con neutro a tierra
Determinar la potencia activa (MW) y reactiva (MVAr) que recibe el motor en estas condiciones.

4.19. El sistema de la Figura 4.54 está trabajando en las condiciones indicadas, esto es, el motor está
recibiendo una potencia de 100 MW, Factor de Potencia 0,95 inductivo a tensión nominal en bornes, cuando
es necesario cambiar el transformador de la fase “a” del banco trifásico. El nuevo transformador monofásico
es de 40 MVA, pero su reactancia es del 7% en base propia. Determine las corrientes que circulan por las
líneas, en estas condiciones. Los datos están en base común 120 MVA.

G
Xn
M
120 MVA
120 MVA
X'd = X1 =
15% X2 = Banco trifásico formado X1 = X 2 =10%
por 3 Unidades
10% X0 = 6%
Monofásicas
X0 = 6%
Cada T/F Monofásico
Xn = 1,5% X1 = X2 = X0 =4%
40 MVA

Figura 4.54

4.20. El sistema de la Figura 4.55, está 

trabajando en las condiciones indicadas, V=100%
G
esto es, la carga, conectada en estrella con
Xn
neutro a tierra, está recibiendo la potencia
de (80+j60) MVA a voltaje nominal en sus 100 MVA
bornes, cuando se abre la fase “a” en este X'd = X1 = 20% 100 MVA (80+j60) MVA
punto. Representando la carga como una X2 = 10% X1 = X2 = X0 =10%
impedancia constante, con Z1C=Z2C=Z0C, X0 = 5% Y11
determine los voltajes V& a ; V& b y V& c Xn = 3%
en los bornes del generador en estas Figura 4.55
condiciones. Los datos están en base 100
MVA.
 CAPITULO 5

ESTABILIDAD

5.1. Introducción

Se dice que un sistema de potencia está en una condición de operación de estado estable si todas las
cantidades físicas que se miden (o se calculan) y que describen la condición de operación del sistema, se
pueden considerar constantes para propósitos de análisis. Si, cuando se está en una condición de estado
estable, ocurre un cambio repentino o una secuencia de cambios en uno o más parámetros del sistema o en
una o más de sus cantidades de operación (variables), se dice que el sistema experimenta un disturbio o una
perturbación de su condición de operación de estado estable. Las perturbaciones pueden ser grandes o
pequeñas de acuerdo con su origen.

Una perturbación grande es aquella para la cual las ecuaciones no lineales que describen la dinámica
del sistema de potencia no se pueden linealizar de forma válida para propósitos de análisis. Las fallas en los
sistemas de potencia, los cambios repentinos y grandes de carga, la pérdida de unidades generadoras y las
maniobras en líneas son ejemplos de perturbaciones grandes y se estudian bajo el nombre de estabilidad
transitoria. Los estudios de estabilidad transitoria, normalmente se hacen en base a la primera oscilación lo
que significa considerar un tiempo de hasta un segundo.

Si el sistema de potencia está operando en una condición estable y experimenta un cambio que puede
ser analizado de manera apropiada a través de versiones linealizadas de sus ecuaciones dinámicas, se dice que
ha ocurrido una perturbación pequeña. Como ejemplos, podemos mencionar, un cambio pequeño y gradual
de carga, un cambio en la ganancia de un regulador automático de voltaje en el sistema de excitación de una
gran unidad generadora, etc.; los que se estudian bajo el nombre de estabilidad permanente. Los estudios de
estabilidad permanente consideran múltiples oscilaciones lo que significa tiempos bastante mayores que los
de la estabilidad transitoria (del orden de los minutos) y por lo tanto, en algunos casos pueden ser importantes
los efectos de los sistemas de control de las unidades generadoras.

En resumen entonces, la estabilidad es la propiedad de un SEP o de sus partes componentes de

mantener un estado de equilibrio (sincronismo), cuando ha sido sometido a acciones perturbadoras. El
concepto puede ser aplicado a una o un grupo de máquinas sincrónicas para señalar la condición de que ellas
permanecen en sincronismo respecto de otras cuando se producen perturbaciones.

5.2. Estabilidad transitoria

5.2.1. Objetivos

Los estudios de estabilidad transitoria suministran la información necesaria para conocer la capacidad
de un SEP de permanecer en sincronismo durante grandes perturbaciones tales como: cambios bruscos
momentáneos o sostenidos de grandes cargas, pérdidas de generación, pérdida de líneas importantes y fallas
tipo cortocircuitos y fases abiertas. Específicamente se obtienen:

 en las máquinas síncronas: Los cambios de tensión, corriente y potencia, velocidad y torque.
 en la red del SEP: Los cambios en las tensiones de las barras y en el flujo de potencia a través de las
líneas.

5.2.2. Suposiciones básicas para los estudios simplificados de estabilidad transitoria

 Usualmente, los generadores se analizan individualmente. Sin embargo, en sistemas multimáquinas

se pueden concentrar grupos de ellos (máquinas equivalentes) para realizar los estudios. El
comportamiento de cada generador se describe mediante una ecuación diferencial, denominada
 89
ecuación de oscilación. Cada generador se representa por una fem constante tras la reactancia
transitoria.
 Como variable fundamental se emplea la posición angular del rotor .
 Debido a la gran inercia de las máquinas, las variaciones de velocidad son pequeñas, por lo que los
elementos estáticos del SEP (líneas, transformadores, etc.) se suponen operando a frecuencia nominal
y las tensiones, corrientes y potencias se calculan mediante ecuaciones algebraicas.
 Los motores síncronos se representan como una máquina más.
 Generalmente no se consideran los efectos de amortiguación.
 En el caso de una perturbación balanceada, las ecuaciones algebraicas de acoplamiento son las
correspondientes al flujo de potencias. Si la perturbación es no balanceada, se deben emplear las
ecuaciones correspondientes a la secuencia positiva, modificando convenientemente la red
equivalente del SEP según el tipo particular de perturbación. Esto significa que en fallas asimétricas,
se utiliza la malla de secuencia positiva, agregándose las de secuencia negativa y cero según
corresponda.

5.2.3. Ejemplo Ilustrativo

Consideremos un generador de rotor cilíndrico suministrando potencia a una barra infinita a través de
una línea de doble circuito, tal como se muestra en la Figura 5.1, donde:

E' : Tensión interna del generador (tensión inducida)

XA : Reactancia total entre la tensión E' del generador y la de la barra infinita V, con ambas
líneas conectadas.
XB :Reactancia total entre las tensión E' del generador y la de la barra infinita V, con una línea
desconectada.
V 0º
G Pg

E' 
X'd
Barra infinita
Figura 5.1.- Sistema generador-linea-barra
infinita

Las ecuaciones potencia ángulo para cada caso son:

PgA V E'
 sin 
XA
V E' (5.1)
PgB  sin 
XB
La Figura 5.2 muestra las curvas potencia ángulo para las ecuaciones (5.1).

Supongamos que inicialmente el generador está entregando una potencia P g0 correspondiente a un

ángulo 0 (curva PgA). Si en estas condiciones se corta bruscamente una línea, entonces el generador pasa a
operar según la curva PgB. Debido a la inercia del rotor, el ángulo  no puede cambiar instantáneamente y en
consecuencia el punto inicial de operación en la curva P gB es a0 y el generador entrega una potencia P1  Pg0.
Como la potencia motriz no ha disminuido, existe un P = Pg0 - P1  0 que acelera el rotor por sobre la
velocidad sincrónica, alcanzándose un punto tal como a 1 en el cual la velocidad es aún mayor que la
sincrónica y en consecuencia este punto es sobrepasado pudiéndose llegar al punto a 2 o bien al punto a3.
 Pg

PgA

M
P a3
3 P
P2 gB
Pmec= Pg0 a2 F
a1
P1
a0

 0 1 2 /2 3  1  

Figura 5.2.- Curvas potencia-ángulo, para el sistema de la Figura 5.1

En estas condiciones el generador suministra mayor potencia que la recibida de la máquina motriz a
expensas de su energía cinética, o sea existe:

P  Pg0  P2  0 par   /2

a
P  Pg0  P3  0 (5.2)
para /2    ( - 1
)

Se produce entonces una desaceleración y el generador vuelve al punto M si había llegado al punto
a3, pero aún P  0 y el rotor continua desacelerándose (ya bajo la velocidad sincrónica), alcanza y sobrepasa
el punto a1 produciéndose ahora una P  0 que acelera nuevamente el rotor. El resultado de este proceso es
una serie de oscilaciones alrededor del punto a1, las que se van amortiguando a medida que transcurre el
tiempo debido a factores tales como: acción electromagnética, fricción, etc.; hasta que se alcanza el punto de
equilibrio (Pg0, 1). En estas condiciones se dice que el sistema generador-línea-barra infinita tiene
"estabilidad transiente".

Supongamos ahora que debido a la falla, el generador se acelera de tal forma que alcanza y sobrepasa
el punto F correspondiente al ángulo -1. En este caso, al sobrepasar el punto F, el generador queda
sometido a un P  0 y continua acelerándose hasta que pierde el sincronismo. De lo anterior se concluye que
el sistema generador-línea-barra infinita tiene "inestabilidad transiente" al sobrepasar el punto F.

Luego, el límite de estabilidad transiente está en  =  - 1

Si   ( - 1)  Operación estable después de un cierto número de oscilaciones

Si   ( - 1)  Operación inestable

Puesto que un ángulo 1 de operación se alcanza después de un tiempo determinado, el análisis

anterior puede visualizarse también a partir de la denominada "curva de oscilación" (Figura 5.3), para lo cual
es necesario resolver la "Ecuación de oscilación" que se estudiará más adelante.


  1 1
 table

I
n
e
s Estable

0

Figura 5.3.- Curvas de oscilación de una máquina síncrona

5.2.4. Ecuación de oscilación de una máquina síncrona

Consideremos el i-ésimo generador de un SEP que contiene m unidades de la Figura 5.4, donde:

Pmi : Potencia mecánica suministrada al generador, vía eje del rotor.

Pgi : Potencia eléctrica generada y suministrada a la red, vía barras del SEP.

 Si Pmi = Pgi, el generador funciona a velocidad constante

 Si Pmi  Pgi, la diferencia de estas potencias se empleará en :

a. Modificar la energía cinética o la velocidad de la unidad motriz-generador.

b. Vencer el torque de amortiguación que se desarrolla fundamentalmente en los enrollados de
amortiguación de la máquina.

Pmi

Generador Pgi

Figura 5.4.- Generador i en un sistema de m unidades

Es decir:
d
P P (EC )  P
(5.3)
mi gi i Di
dt

donde:

ECi : Energía cinética de la máquina i

PDi : Potencia de amortiguación de la máquina i
 Energía cinética ECi: representa la energía cinética total de la unidad i. Se expresa habitualmente en
Megawatt-segundo (MW-s) o MegaJoule (MJ) de acuerdo a la expresión (5.4).
2
ECi  1 Ji i (5.4)
2

Ji : Inercia de la i-ésima máquina

i : Velocidad angular del generador i
Derivando ECi respecto al tiempo y haciendo Mi = Ji i se obtiene
d
(EC )  J 
i i i di d2i
dt  Mi (5.5)
dt dt2
donde:

i : Posición angular del rotor de la máquina i respecto a un eje fijo en grados o radianes
eléctricos.
Mi : Momento angular del Máquina i

Se ha determinado que es más conveniente medir la posición angular del rotor i y su velocidad
angular i respecto a un eje de referencia que gira a velocidad sincrónica, que hacerlo respecto a un eje
estacionario. De este modo, si consideramos la Figura 5.5, donde s es la velocidad sincrónica y i es el
desplazamiento angular del rotor respecto al eje que gira a la velocidad sincrónica, se puede escribir:

i  i  s  i  s t (5.6)
y por lo tanto:
d
(EC )  M d
i 2
i (5.7)
dt i
dt2 s
i
Tmi i
 s = s t

Tgi

Figura 5.5.- Representación esquemática de un generador síncrono

Torque o potencia de amortiguación: A medida que la velocidad del rotor se desvía respecto a la
sincrónica, se inducen corrientes en los enrollados de amortiguación del rotor que desarrollan un torque
resistente. Este torque resistente aumenta con la velocidad relativa (di/dt) y usualmente se supone
proporcionalidad entre el torque o la potencia PDi y la velocidad, a través de un parámetro positivo Di del
generador medido en MW/(rad elect/seg) o sea:

d i
PDi  Di (5.8)
dt
 Por lo tanto, la ecuación de oscilación de la máquina i queda:

Pmi  Pgi  Mi d 2
d i
 Di (5.9)
i dt
dt2

Debido a que el trabajo básico del enrollado de amortiguación es suministrar un torque extra de
estabilización, normalmente se desprecia su efecto en la ecuación de oscilación. De este modo se obtienen
resultados más seguros en el análisis y se simplifican los cálculos. Según esto, la ecuación de oscilación es:
d 2 i
Pai  Pmi  Pgi  Mi (5.10)
dt2

donde Pai  Pmi  Pgi , es la potencia de aceleración.

Constante de Inercia H: Se define como la razón entre la energía cinética almacenada por la
máquina girando a velocidad síncrona y la potencia aparente nominal de la máquina. Es decir:

ECi
H i (5.11)
Gi

Donde Gi es la potencia aparente nominal de la máquina. Si EC i se mide en MegaJoule y Gi en MVA,

H resulta medida en segundos. El valor de H varía muy poco con el tamaño de la máquina (2-8 segundos); en
cambio M depende del tamaño y tipo de máquina, por lo que se prefiere utilizar H en la ecuación (5.10).

A partir de (5.4) y considerando (5.11), se puede escribir:

2
ECi  1 Jii  1 Mii  Hi Gi (5.12)
2 2

de donde se obtiene:

Gi Hi
M i Gi Gi Hi (5.13)
1 i
 H  180 f
f
2 i

Según el ángulo i se mida en radianes eléctricos o grados eléctricos, respectivamente y donde f

representa la frecuencia de operación del sistema. Por lo tanto la ecuación (5.10) queda:

G i Hi d 2 i
Pai  Pmi  Pgi (5.14)
180 f d t 2

Con i medido en grados eléctricos. La expresión (5.14) puede expresarse en tanto por unidad base
propia, dividiéndola por la potencia aparente nominal de generador. De esta forma se tiene finalmente:
H i d 2 i
Pai(pu)  P mi(pu)  P gi (pu)  (5.15)
180 d t 2
f

La potencia eléctrica entregada por el generador durante la perturbación depende de la fem detrás de
la reactancia transitoria E' y de la corriente I, de manera que:

Pgi  Re al E& ' I * (5.16)


 La corriente I depende del comportamiento de toda la red eléctrica. Como la potencia mecánica
permanece prácticamente constante mientras no se modifique el control de velocidad (en el primer segundo)
las variaciones de frecuencia son pequeñas, de manera que se puede suponer sin gran error que todos los
elementos estáticos del sistema (líneas, transformadores, etc.) operan en condiciones cuasi estacionarias,
hecho que simplifica los cálculos. En vez de la corriente I puede usarse también la tensión en bornes V (o la
tensión constante en alguna otra barra cercana).

5.2.5. Ecuaciones de oscilación en un sistema multigenerador

Consideremos un SEP alimentado por m generadores como se indica en la Figura 5.6, donde se
supone que las reactancias transientes de los generadores están incluidas en la red pasiva del sistema.

Empleando el método nodal de resolución de circuitos, para la red de un SEP de m barras (nudos) se
puede escribir:

IB   YB (5.17)

EB 
donde: IB  es el vector de corrientes inyectadas en las barras; EB  vector de tensiones de barras e
YB  es la Matriz de admitancia de barras, que se pueden escribir como:
I B   &I &I &I 

2 i m
&I t (5.18)
1 ··· ···

E B  
E& E& 
' '
1 E& ' (5.19)
E& ' ···
2
···
i
t
m

 Y11 Y1m 
Y21 Y12 ··· Y1i
Y22 ··· ··· Y2m 
 Y2i ··· 
Y B    M M
M M
Yim 
(5.20)
 Yi1 Yi2 Yii
··· ···
 
M M M M 
Ym1 Ym2 Ymi ··· Ymm 
···

I1
E'1
RED PASIVA DEL SEP
I2
E'2
E'i Ii

Im
Em'

Figura 5.6.- Sistema multigenerador

Considerando las ecuaciones (5.17) a (5.20), se puede escribir:

&I  Y E& '
Y &'
&E' 
Y ·····  Y E  ·  &E' E&
m
' Y (5.21)
i i1 1
i2 2 ii i im m  ij j
j1
por lo tanto, la potencia compleja suministrada por el i-ésimo generador es:
S& 
P
E& *I *  E& m
' '
 jQ  Y*E ' (5.22)
gi gi
gi i i i  ij j
j1
 suponiendo que:

E& 'j  E'  j Y& ij  ij (5.23)

j Yij

Se tiene que:
P  jQ  m E' E' Y (   ) (5.24)
gi gi i j ij i j ij
j1

de donde:
m
Pgi   E'i E' jY ijcos( ij   j ) i (5.25)
j1

Es decir, la potencia suministrada a la red por una unidad cualquiera, depende de las posiciones
angulares j de todo el resto de los generadores, además de la propia i, ambas funciones del tiempo. La
admitancia Yij sufre cambios discontinuos debido a los cambios topología de la red (prefalla, en falla, falla
despejada, reconexión, etc.)

En el cálculo de Pgi se ha supuesto que las velocidades de las máquinas permanecen constantes, en
consecuencia se tendrá un sistema de ecuaciones diferenciales simultáneas de la forma:
d 21
Pm1  Pg1 (1 ,  2 ,·····,  i , ,  m )  M1
dt2
d 2 2
Pm2  Pg2 (1 ,  2 ,·····,  i ,· ,  m )  M (5.26)
dt2
2 d 2i
Pmi  Pgi (1 ,  2 ,·····,  i ,· ,  m ) 
Mi dt2
d 2m
Pmm  Pgm (1 ,  2 ,·····,  i ,·,  m )  M
m dt2

En general, la solución del sistema de ecuaciones diferenciales anterior resulta laboriosa, aún para el
caso simplificado de representación de los generadores aquí supuestos.

5.2.6. Estudio de la ecuación de oscilación de un generador conectado a una barra infinita

Esta situación, si bien es cierto, es la más simple posible, tiene importancia desde un punto de vista
didáctico y práctico ya que muchas situaciones reales pueden ser reducidas a esta forma simple equivalente.
Según (5.15) , la ecuación de oscilación en este caso es de la forma:
d 2 H d 2
 (5.27)
Pa  Pm  Pg  M 2 180 f d t 2
dt
donde Pa ,Pm , Pg y M están en “tanto por unidad” y  en grados eléctricos.

– Potencia mecánica suministrada al generador: Los cambios de Pm dependen enteramente de la

acción de los controles potencia-frecuencia. En el instante de la perturbación e inmediatamente
después hasta el lapso de 1 segundo aproximadamente, la potencia de la máquina motriz P m es
sustancialmente constante, o sea:

Pm (0 )  Pm (0 )  Pm
(5.28)
– Potencia suministrada por el generador durante la perturbación al sistema: En régimen
permanente, la potencia activa suministrada por un generador de polos cilíndrico o de polos salientes
está dada respectivamente por:
EV
P sin
g
Xd (5.29)
EV Xd  Xq
P sin  V 2sin2
g
Xd 2X d X q

En los estudios simplificados de estabilidad transiente, para calcular la potencia activa suministrada
por un generador, ya sea de polos cilíndricos o de polos salientes, se utiliza el circuito equivalente
aproximado que se muestra en la Figura 5.7, en cuyo caso:
jX d'

+
+
E' V
-
-

Figura 5.7.- Circuito equivalente de una máquina síncrona para estudios de estabilidad transiente
E' V
P sin
g
X' d (5.30)

– Análisis de la ecuación de oscilación: Supongamos por ejemplo, que la línea que une el generador y
la barra infinita (o sea el resto del sistema) sufre un cortocircuito trifásico. Los relés de protección
harán desconectar la línea y una fracción de segundo después, la línea será reconectada. Suponiendo
que la falla ha terminado se volverá a la situación original o se producirá el fenómeno de inestabilidad
según sea la situación que corresponda y que se analizará luego. Consideraremos dos períodos:

 Durante la falla: En este período Pg=0 y Pm=constante y la ecuación (5.27) queda:

d 2
Pm  M 2
(5.31)
dt
La solución de esta ecuación diferencial, considerando que en t=0,  = 0 y d/d t = 0, queda:
Pm
(t)  0  t2 (5.32)
2M

En la Figura 5.8 se ha dibujado esta solución, donde se aprecia que la máquina aumenta
indefinidamente su ángulo de torque. En consecuencia se producirá una aceleración brusca de la máquina,
perdiéndose rápidamente el sincronismo, si no se elimina la causa; es decir si no actúan las protecciones.

 Reconexión de la línea: Inmediatamente después que se reconecta la línea (y suponiendo que la falla
ha sido eliminada) la potencia Pg pasa bruscamente al punto “a” de la Figura 5.8, en que:
E' V
P sin
g r (5.33)
X
 X corresponde a la reactancia existente entre E' y V es decir, incluye la reactancia de la línea. Se
establece entonces una potencia desacelerante P des que depende de , es decir:

Pdes E' V
  sin (5.34)
Pm X

Pg
a

Pméc =Pg0

0 p r /2  lím

tp

tr

sin reconexión

t
con reconexión

Figura 5.8.- Solución de la Ecuación de oscilación para una máquina síncrona

Por lo tanto, la ecuación de oscilación en este período es:

Pm  E' d 2
sin  M (5.35)
V dt2
X
que no tiene solución analítica simple y en su resolución se deben utilizar métodos numéricos, pero se
comprende que  crecerá a menor velocidad y que el resultado final depende del instante de reconexión. Si la
reconexión ocurre oportunamente  deja de crecer y a pesar de que lím  /2, se conserva el sincronismo, tal
como se muestra en la Figura 5.8.

Es conveniente hacer notar que en muchos casos, durante la falla, el generador sigue suministrando
potencia (cortocircuitos monofásicos, bifásicos, fases abiertas), es decir, P g  0 durante este período y por lo
tanto, la ecuación de oscilación en esta etapa, tendrá la misma forma indicada en (5.35).

Ejemplo 5.1. El sistema de la Figura 5.9 está entregando la potencia indicada a la barra infinita, cuyo voltaje
es el nominal, cuando se produce un cortocircuito monofásico a tierra en el punto P (al medio de la línea). La
falla es despejada simultáneamente por ambos interruptores (en forma monopolar). Determine las ecuaciones
de Pg() antes de falla, en falla y en falla despejada. Todos los datos en % están en base 50 MVA.
Línea
X1=X2=20% 
E'  X0=20%
1 2 (50+j20) MVA V
X'd =X1 =30% P
Barra
X2=15% X1 X1 Infinita
X0=10%
=X2=X0=5% =X2=X0=5%
Figura 5.9
Solución:

a) Antes de falla o en condiciones de prefalla: El circuito equivalente por fase se muestra en la Figura 5.3

En el circuito de la Figura 5.10, con V&  10º a j0,3 j0,05 j0,2 j0,05 b
, se tiene:
I
&  S  1  j0,4 + +
*

I   1,077  21,8º

 V  10º E' V
E& '  10º j0,6 1,077  21,8º  1,377525,82º
E'  V 1,3775 *1
P   sin    sin   2,2958 sin 
gAF
X ab 0,6
Figura 5.10

b) En falla: Las mallas de secuencia cero y negativa se muestran en las Figuras 5.11 y 5.12.

- Secuencia cero - Secuencia negativa

j0,05 j0,1 j0,1 j0,05 j0,1
j0,15 j0,05 j0,1 j0,05

j0,1 X0 X2

Figura 5.11 Figura 5.12

Las impedancias de secuencia cero y negativa - Malla de secuencia positiva incluyendo las
son respectivamente: impedancias de secuencia cero y negativa

0,15 j0,3 j0,05 j0,1 j0,1 j0,05

a b
X  0,075 (pu)
0
2
0,3  0,15 + jX0 +
X  0,1 (pu)
2
0,3  0,15 E' V
jX2
Estas reactancias se agregan a la malla de
secuencia positiva en serie en el punto de falla,
tal como se indica en la Figura 5.13 Figura 5.13
La reactancia Xab se encuentra utilizando la expresión de transformación de estrella a delta en la Figura 5.13;
es decir:

Xab  0,45  0,15  0,15  0,175  0,175  0,45  0,9857

0,175
E'  V 1,3775 1
P   sin  

 sin   1,3975 sin 
gF
X ab 0,9857

c) En falla despejada: Queda una fase abierta en dos puntos. Las mallas de secuencia cero, negativa y
positiva son las indicadas en las Figuras 5.14 a 5.16
 - Secuencia cero - Secuencia negativa
j0,05 P 0 Q j0,05 j0,15 j0,05 P2
0 Q2 j0,05

j0,1

Figura 5.14 Figura 5.15

Las mallas de secuencia quedan conectadas tal - Malla de secuencia positiva
como se muestra en la Figura 5.17, que
j0,3 j0,05 P1 Q1 j0,05
corresponde a la Figura 4.45, donde: a b

E'  E' y E '  V + +

A B
Z1  j0,35; Z  j0,05'
E' V
1
Z 2  j0,2; Z  j0,05
'
2
Z0  Z'  j0,05
0
Figura 5.16
ZA= y ZB corresponde a la impedancia de la línea que está abierta en dos puntos (se abrieron los dos
interruptores simultáneamente); es decir:
ZB  jX1  jX 2  jX 0  j0,2 (pu) 1/3(ZA - ZB )

Por lo tanto, la reactacia Xab entre E’ y V queda: Z1 Ia1 P1 ZB Z'1

a Q1 b
+ +
(X 2  X B  X 2)  (X 0  X B  X 0 )
' '
' + Va1 -
X ab  X1   X B  X1 EA ' EB'
(X 2  XB  X2' )  (X 0  XB  X' )0 - -
Z2 Ia2 P2 ZB Q Z'2
2
(0,2  0,2  0,05)  (0,05  0,2  0,05) + Va2 -
X ab  0,35  (0,2  0,2  0,05)  (0,05  0,2  0,05) 
Z0 Ia0 ZB Q0 Z'0
 0,2  0,05  0,78 P0
+ Va0 -
E'  V 1,3775 1
P   sin    sin   1,766 sin 
gFD
X ab 0,78 Figura 5.17

5.2.7. Algunos métodos de solución de la Ecuación de Oscilación

Existen diferentes métodos para la evaluación numérica de las ecuaciones planteadas en (5.35), tales
como Punto a Punto, Euler, Euler Modificado, Runge-Kutta, etc. La mayoría de ellos son prácticos solamente
cuando se emplean computadores, sobre todo si se estudia sistemas de gran tamaño. Se trata de determinar 
como función de t, graficando la respuesta y de esta forma saber si el sistema es estable o no. Se estudiará a
continuación algunos de ellos, con el fin de ver sus características y posibilidades de aplicación.

a. Método punto a punto (Solución por partes): Es un método simple, que permite realizar cálculos a
mano y por lo tanto es aplicable sólo a sistemas pequeños. Se divide el tiempo total de estudio en n intervalos
de duración  t segundos cada uno, tal como se indica en la Figura 5.18. Generalmente se utiliza un  t =
0,05 segundos y el cálculo se hace bajo las siguientes suposiciones:
– La potencia de aceleración determinada al
comienzo de un intervalo, es constante desde Pa
Pa,n-
la mitad del intervalo anterior hasta la mitad Supuesto
del intervalo considerado. En algunos casos, 2
Pa,n-1 Real
en vez de la potencia de aceleración se emplea
la aceleración angular.
Pa,n

– La velocidad angular es constante en cada

intervalo e igual al valor calculado para la n t
mitad del mismo. n- n-1
' 2
Por supuesto, ninguna de las condiciones  '
n '
anteriores es exacta ya que  está cambiando
continuamente y tanto Pa como ' son -1/2
funciones de . A medida que el intervalo de 'n-3/2
tiempo disminuye, la curva de oscilación Supuesto
Real
calculada de esta forma se hace más exacta.
t
La Figura 5.18 ayuda a visualizar estas n-3/2 n-
suposiciones. La potencia de aceleración se calcula  1/2
en los puntos encerrados en círculos en los
extremos de los intervalos n-2, n-1 y n, que son los
comienzos de los intervalos n-1, n y n+1  n
respectivamente. La velocidad angular '
corresponde a d/dt, es decir al exceso de la n-1
velocidad angular de la máquina, sobre la velocidad
síncrona
un cambiosde
. Entre las ordenadas
velocidad originadon-3/2 y n-1/2
por la hay
potencia t
n-2 n.1 n
de aceleración constante.
Figura 5.18.- Valores supuestos y reales de Pa, ' y 
El cambio de velocidad es igual al producto de la aceleración por el intervalo de tiempo:
d

2
180 f
'n1/ 2 'n3 / 2  t Pa,n1  t (5.36)
2 H
dt
La variación del ángulo  en un intervalo cualquiera es igual al producto de la velocidad ' en el
intervalo por el tiempo. Así, el cambio de  durante el intervalo n-1 es

 n1  n1  n2  'n3 / 2  t (5.37)

y durante el intervalo n

 n  n  n1  'n1/ 2  t (5.38)

Restando (5.37) a (5.38) e introduciendo el resultado en (5.36), se obtiene:

 n 180 f
 k Pa,n1 con k= ( t)2 (5.39)
n 
 1 H
Luego: n  n1 
 n

(5.40)

La ecuación (5.40) permite obtener  como función del tiempo o sea corresponde a la solución paso a
paso de la ecuación de oscilación. Por otra parte, la velocidad ' se puede determinar a partir de (5.39),
dividiendo por t y se obtiene:

k
'n 
n1  Pa,n1 (5.41)
' 
t

 Discontinuidad en la potencia de aceleración: Cuando ocurre una falla, se produce una

discontinuidad en la potencia de aceleración P a que tiene un valor cero antes de la falla y un valor distinto de
cero después de ésta. Esta discontinuidad ocurre al comienzo del fenómeno (cuando t = 0). Lo mismo sucede
cuando se producen aperturas de interruptores, reconexiones, etc. Teniendo en cuenta que este método supone
que la potencia de aceleración calculada al comienzo del intervalo es constante desde la mitad del intervalo
anterior hasta la mitad del intervalo que está siendo considerado y que en este caso se tiene dos valores
distintos para la potencia acelerante, se debe tomar el valor promedio de estos valores como la potencia
acelerante constante. En general, para una discontinuidad en un tiempo t se tiene:
Pa (t  )  Pa (t  )
(5.42)

En el caso en que la discontinuidad ocurra en el punto medio del intervalo no hay necesidad de
emplear (5.42) pues el método contempla una discontinuidad justamente en ese punto. En otro caso, conviene
aproximar al más cercano, esto es, al comienzo o al medio del intervalo, según corresponda.

Ejemplo 5.2. Las ecuaciones de Pg=Pg() de un sistema generador-barra infinita son las siguientes: P gAF=3 sin
 (pu); PgF=0,5 sin  (pu) y PgFD=1,5 sin  (pu). La potencia mecánica de entrada es de 1,2 (pu), la frecuencia
es de 50 Hz y la constante de inercia es de 4 seg. Si la falla se despeja a los 0,2 seg de iniciada y se realiza
una reconexión exitosa 0,15 seg después del despeje:
a. Dibujar las curvas de Pg=Pg(), indicando en forma aproximada, las áreas acelerantes y desacelerantes.
b. Determinar (utilizando exclusivamente la curva de oscilación) si el sistema es estable,.

Solución:
Pg
a) Curvas de Pg=Pg(): Se Acelerante Desacelerante
muestran en la Figura 5.19, 3,0
donde: P gAF

0  sin 1 (1,2 / 3)  23,5782º 2,0 P gFD

lím  180ºsin 1 (1,2 / 3) 1,5 Pm=P g0
1,2
lím  156,4218º P gF
0,5
b) Estabilidad: Se debe 0
resolver la ecuación de 0 d 90 r  lím 180 
oscilación. Usando el método
punto a punto, con t=0,05 seg, Figura 5.19
se tiene:

180 * f * (t)
2
180 * 50 *  5,625
K 
(0,05) 2
H
4
 102

Tabla 5.1. Aplicación del Método punto a punto

T Pm Pg Pa K*Pa   ’
0- 1,2 1,2 0 - 0 23,5782 0
0+ 1,2 0,2 1 - 0 23,5782 0
0p - - 0,5 2,8125 0 23,5782 0
0,05 1,2 0,2222 0,9778 5,5001 2,8125 26,3907 0,9817
0,1 1,2 0,2847 0,9153 5,1486 8,3126 34,7033 2,9016
0,15 1,2 0,3725 0,8275 4,6547 13,4612 48,1645 4,6988
0,2 1,2 0,4578 0,7422 4,1749 18,1159 66,2804 6,3236
0,2- 1,2 0,4578 0,7422 - 18,1159 66,2804 6,3236
0,2+ 1,2 1,3733 -0,1733 - 18,1159 66,2804 6,3236
0,2p - - 0,2845 1,6003 18,1159 66,2804 6,3236
0,25 1,2 1,4963 -0,2963 -1,6667 19,7162 85,9966 6,8823
0,3 1,2 1,4552 -0,2552 -1,4355 18,0495 104,0461 6,3005
0,35 1,2 1,2903 -0,0903 -0,5081 16,6140 120,6601 5,7994
0,35- 1,2 1,2903 -0,0903 - 16,6140 120,6601 5,7994
0,35+ 1,2 2,5806 -1,3806 - 16,6140 120,6601 5,7994
0,35p - - -0,7355 -4,1372 16,6140 120,6601 5,7994
0,4 1,2 2,1892 -0,9890 -5,5631 12,4770 133,1371 4,3553
0,45 1,2 1,9263 -0,7263 -4,0854 6,9139 140,0510 2,4134
0,5 1,2 1,8105 -0,6105 -3,4341 2,8285 142,8795 0,9873
0,55 1,2 1,8357 -0,6357 -3,5758 -0,6056 142,2739 -0,2114

Conclusión: Como el ángulo deja de crecer y empieza a disminuir, el sistema es estable transitoriamente.
b. Método de Euler: Constituye un algoritmo relativamente sencillo, aunque poco preciso para
determinar y en la ecuación diferencial dy/dt = f (y). La fórmula de recurrencia es:
yi1  yi   t f(yi ) (5.43)
Es decir, el valor al final del intervalo se obtiene a partir del que se tiene al comienzo del mismo y de
la derivada evaluada también al inicio de éste. Esto hace que sea poco preciso cuando la derivada cambia con
rapidez dentro del intervalo, lo que obliga a trabajar con intervalos pequeños.
En este problema, el tiempo total t se divide en n intervalos de  t seg y se supone potencia acelerante
constante en cada intervalo e igual al valor que tiene al comienzo. La ecuación (5.35), se puede escribir:
d 2  Pm  Pg ()
 (5.44)
dt2 M
para luego transformarla en dos ecuaciones de primer orden de la forma:
d
 '
dt
d' Pm  Pg () (5.45)
dt  M
cuyas ecuaciones de solución de acuerdo con (5.43), son:
   '  t  180 
i1 i i  

Pm Pg (i )
 (5.46)
'  '  t
i1 i
M
 Como Pg( ) sufre cambios discontinuos es necesario acomodar  t para que cuando ocurra algún
caso de éstos, se quede al comienzo del intervalo, lo que significa dejar flexibilidad en cuanto al intervalo.

c. Método de Euler Modificado: En este método, la ecuación diferencial dy/dt = f (y) se resuelve
utilizando la expresión:
f (y )  f (y )
yi1 
 yi  i1 i
t (5.47)
2
donde y*i+1 es un valor auxiliar que corresponde a una predicción del valor de y i+1, que se obtiene
utilizando la fórmula de Euler, ecuación (5.43). El valor y i+1 de (5.47) es el valor mejorado o corregido.

Aplicado al problema dado por (5.45), la predición se obtienen usando (5.46) y luego se aplica (5.48):
'* ' 180
i1  i  i1 i
t
2  (5.48)

'i1 2Pm  Pg (*i1 )  Pg (i )

 i1 t
'  2M

d. Método de Runge-Kutta de 4º orden: Este método es más conveniente que los anteriores, ya que
mejora la precisión y se pueden considerar intervalos mayores de tiempo, lo que significa resolver una menor
cantidad de flujos de potencia en el sistema. Aunque la cantidad de cálculos es mayor en cada punto, la
disminución en el número de veces que hay que resolver las ecuaciones algebraicas, lo hace de todas formas
muy conveniente Las expresiones quedan:
i1  i  1 (D1  2D2  2D3  D4 )  t
6
(5.49)
'i1  'i  61 (W1  2W2  2W3  W4 )  t

con:
D1  'i


W  s P  P ( )  D ' 
1 m g i i
Ta
W t
D  '  1
2 i
2
s  t   t
W2  P mP g( i D 1 )   D' iW 1 

T  2  2
a  
t (5.50)
D  ' W
3 i 2
2
s  t   t
W3  P m Pg ( i  D 2 )   D 'i W 2 
T  2  2
a  
D4  'i W3  t


W  s P  P (  D  t)  D' W  t  
4 m g i 3 i 3
Ta
En las expresiones anteriores se ha incluido la constante asociada a los roces mecánicos y efecto de
los enrollados amortiguadores D y el tiempo de arranque del generador o constante de aceleración T a,
definido como el tiempo necesario para alcanzar la velocidad síncrona, al partir con torque nominal
constante. El tiempo de arranque Ta equivale a 2 veces la constante de inercia H.
5.2.8. Criterio de Areas Iguales

La curva de oscilación permite determinar si un sistema formado por un generador conectado a una
barra infinita (o una situación equivalente) es estable o inestable después de una perturbación brusca e
importante. Si la curva muestra que el ángulo tiende a crecer sin límite, el sistema es inestable, por otra parte,
si después de todas las perturbaciones (incluyendo desconexión y reconexión de líneas, por ejemplo), el
ángulo alcanza un máximo y luego disminuye, puede decirse que el sistema es estable.

En general, la estabilidad o inestabilidad del sistema depende en forma decisiva del hecho de si la
falla es sostenida (permanente) o si se elimina en un tiempo determinado. Existe un tiempo máximo
permisible para eliminar la falla tal que el sistema sea estable. Este tiempo medido respecto a la iniciación de
la falla, se denomina "tiempo crítico de aclaramiento".

El emplear la curva de oscilación en forma exclusiva para determinar el tiempo crítico de despeje, por
ejemplo, es un método de aproximaciones sucesivas, por lo que no es eficiente. El método consiste en
calcular la curva de oscilación para falla sostenida y luego para distintos tiempos de despeje supuestos. Del
análisis de todas las diferentes curvas se determina el tiempo crítico de despeje.

El criterio de áreas iguales, permite determinar fácilmente el ángulo crítico de despeje y a partir de la
curva de oscilación calculada para falla sostenida (permanente) o empleando un procedimiento indirecto,
calcular el tiempo correspondiente a este ángulo crítico.

Consideremos nuevamente la ecuación de oscilación de un generador conectado a una barra infinita.

d 2
 Pm  Pg (5.51)
M
2
dt
Multiplicando toda la expresión por 2(d/dt)/M y arreglando adecuadamente se puede escribir:

d 2 
M 0
(PmgP )d (5.52)
dt 

donde 0 es el ángulo inicial de la perturbación y d/dt = ' (exceso de la velocidad angular respecto
a la síncrona). Cuando d/dt = 0 significa que ' = 0; es decir, la velocidad del generador es igual a la
velocidad síncrona y se ha alcanzado el valor máximo de la curva de oscilación. Por supuesto la máquina no
permanecerá en reposo respecto a la barra infinita , la primera vez que d/dt se haga cero, pero el hecho que
 deje de variar momentáneamente puede interpretarse como una indicación de estabilidad.

Para d/dt = 0, la expresión (5.52) queda:


 (P  P ) d = 0 (5.53)
0 m g

Ecuación que puede interpretarse gráficamente mediante la curva Pg = f () del generador conectado
a la barra infinita y considerando una perturbación o falla específica.

Supongamos que la línea que une al generador con la barra infinita sufre un cortocircuito trifásico y
que después de un cierto intervalo de tiempo se la reconecta cuando la falla ha sido eliminada. La situación
puede ilustrarse como se indica en la Figura 5.20.

En la Figura 5.20, el área A1 es acelerante (A1 0) y A2 es área desacelerante (A2  0). La ecuación
(5.53) aplicada a este caso resulta:

 (Pm  Pg ) d = A1  (A 2 )  0  A1  A 2 (5.54)
0
 Pg

Pg=Pmáxsin
A2

Pm= Pg0

A1

0 r /21 lím 

Figura 5.20.- Aplicación del criterio de áreas iguales a un sistema generador-línea-barra infinita, falla trifásica

Esto significa que el sistema será estable si después de reconectar la línea en un tiempo
correspondiente a r, el generador oscila hasta un ángulo máximo 1, tal que A1 = A2. Lo anterior constituye
el criterio de áreas iguales para determinar la estabilidad transiente.

Por otra parte si r es tal que 1  lím, el sistema será inestable y A 1  A2, de modo que el valor
máximo de r o sea el ángulo crítico de aclaramiento (de reconexión, en este caso) rc se podrá determinar
aplicando el criterio de áreas iguales según se muestra en la Figura 5.21, de la siguiente forma:

Pg

Pg=Pmáx sin

A2
Pm= Pg0

A1

0 rc lím 

Figura 5.21.- Determinación del ángulo crítico de aclaramiento

 lím  rc  lím
 (P  P ) d =  d   (P  sin) d  0 (5.55)
P P
rc m máx
0 m g 0 m

de donde se obtiene que:

1 P m (lím  0 )  Pmáx cos lím 

 rc  cos  (5.56)
Pmáx 



Ejemplo 5.3. Un generador conectado a una barra infinita tiene las siguientes curvas de P g=Pg(): PgAF=2,5
sin , PgF=0,8 sin , PgFD=1,5 sin . El ángulo de torque antes de falla es de 30º y todos los datos están en pu
base común.
a. Determinar si el sistema es estable, al despejar la falla en un tiempo tal que d=60º, sin reconexión.
b. Si el sistema se reconecta 20º después que la falla ha sido despejada, calcular el ángulo crítico de despeje.
a) Sin reconexión: Las curvas Pg
de Pg() se muestran en la
Figura 5.22, donde: 2,5 P gAF

0  30º  0,5236 rad 2,0

P gFD
Pm  Pgo  2,5 * sin 30º  1,25 1,5
Pm=P g0
1,25
1,0
1,25 
  180ºsin 1
lím  1,5  0,5 P gF
 
0,0
lím  123,56º  2,1565 rad 30 60 90  lím 180 
0 d
d  60º  1,0472 rad Figura 5.22
Aplicando el criterio de Areas iguales se puede escribir:
d  lím
 (1,25  0,8 * sin ) d   (1,25  1,5 * sin ) d   el sistema es inestable
0,169  0
0 d

a) Con reconexión: Las Pg

curvas de Pg() se muestran en
la Figura 5.23, donde: 2,5 P gAF

 0  30º  0,5236 rad

2,0
'
  180º30º  150º P gFD
lím
1,5
lím  2,618 rad 1,25 Pm=P g0

r  cd  20º 1,0

Para determinar el ángulo 0,5 P gF

crítico de despeje se debe 20º
cumplir que: 0,0
30  cd 90 r ' 180 
 'límP  P () d 0

0 lím
m g Figura 5.23
0
cd r lím
Es decir:  (1,25  0,8 * sin ) d   (1,25  1,5 * sin ) d 
cd r
(1,25  2,5 * sin ) d  0
0

0,2399  0,7 * cos cd  cos (cd  20º )  0

F()  0,2399  1,6397 * cos cd  0,342 sin cd  0  cd=86,45º; con F()=8,574 *10-5

5.2.9. Reducción a un sistema simple

Ya se ha dicho que las distintas características eléctricas y mecánicas de las máquinas de un sistema
hacen necesario tratarlas en forma individual. Sin embargo, en algunos casos excepcionales (por ejemplo,
máquinas de una misma central o eléctricamente cercanas entre si y ubicadas lejos de la perturbación) es
posible reemplazar dos o más máquinas por una sola que les sea equivalente. En tal caso, la energía cinética
de la máquina equivalente debe corresponder a:
n n
ECeq   ECi   (5.57)
HiGi
i1 i1

Como las máquinas están en paralelo, sus potencias se suman. Luego, la potencia de la máquina
equivalente Geq es:
n
G eq   (5.58)
Gi
i1

Dividiendo ambos miembros de (5.57) por Geq se obtiene la constante de inercia equivalente, Heq:
n

ECeq  Hi Gi
Heq  n  i1n (5.59)
Gi Gi
i1 i1

para que esta máquina oscile junto con las máquinas que reemplaza. Como las máquinas están en
paralelo en alguna barra del sistema, la reactancia de la máquina equivalente será igual a la combinación en
paralelo de las reactancias individuales hasta esa barra y la fem E' será igual a la tensión transitoria que se
tendría en vacío en la barra de unión.

5.2.10. Ecuación de oscilación de dos máquinas finitas

Dos máquinas síncronas conectadas por medio de una reactancia pueden ser reducidas a una máquina
equivalente alimentando una barra infinita a través de la reactancia. Este modelo se utiliza en algunos casos
con propósitos demostrativos. Al respecto, consideremos la Figura 5.24, en que además supondremos que
existen consumos en las barras, para hacer más general el análisis.

E1' 1 E'2 2 M2
M1

Figura 5.24.- Dos máquinas finitas

Las ecuaciones de oscilación son:

d 1 Pm1  Pg1 Pa1

2
 
dt2 M1 M1 (5.60)
d 2 Pm2  Pg2 Pa 2
2
 
dt2 M2 M2

Es decir, varían los dos ángulos. Es posible reducir estas dos ecuaciones a una sola considerando que
sólo interesa la diferencia entre ellos. Si ésta aumenta indefinidamente, el sistema es inestable. Restando las
ecuaciones (5.60), se tiene:
 d 21 d 22 Pa1 Pa 2
(5.61)
d t 2  d t 2  M1 M 2

haciendo  = 1 - 2 en la expresión (5.61) y amplificándola por M1 M2 / (M1+M2), se puede escribir:

d
2
 Pa1M 2  Pa 2 M1
M  M1  M 2 (5.62)
dt2

donde M = M1 M2 / (M1+M2) es el momento angular equivalente a las dos máquinas

Desarrollando la expresión (5.62) considerando (5.60) se obtiene:
d 2
M  Pm  Pg (5.63)
dt
2

con:
Pm1M 2  Pm2 M1
Pm 
M1  M 2
Pg1M 2  Pg2 M1 (5.64)
Pg 
MM
1 2

La ecuación (5.63) es equivalente a la de un generador conectado a una barra infinita y por lo tanto se
puede aplicar todo lo estudiado en 5.2.6.

5.2.11. Factores que condicionan la estabilidad transitoria

El análisis hecho hasta el momento, con todas las simplificaciones que contiene, permite determinar
los factores que más influyen en la estabilidad transitoria de un sistema, los que se pueden visualizar mejor a
partir del criterio de áreas iguales y considerando el sistema generador-línea-barra infinita de la Figura 5.1.

– Potencia eléctrica inicial: Del criterio de áreas iguales se deduce que mientras mayor sea la potencia
eléctrica que está entregando la máquina antes de ocurrir la falla, mayor será el área acelerante. Ello
incrementará la velocidad relativa que adquiere el rotor, aumentado de esta forma la posibilidad de
que el sistema sea inestable. Un criterio empírico suele ser el de limitar la potencia eléctrica inicial al
80% del máximo teórico E'V/X.

– Tiempo de operación de las protecciones (t p): Mientras más rápidos son los esquemas de protección
y los interruptores utilizados, menor es el ángulo p alcanzado, y por lo tanto, menor el área
acelerante y mayor la potencia que se podría haber estado transmitiendo. Intimamente ligada con este
aspecto está la conveniencia de realizar una apertura simultánea de ambos extremos de la línea
fallada, esto es, de usar equipos de onda portadora.

– Tipo de falla: El valor máximo de la curva de potencia eléctrica transferible durante la falla depende
del tipo y ubicación de ésta. La falla más rigurosa es el cortocircuito trifásico, siguiendo el bifásico a
tierra, bifásico, monofásico y las fases abiertas. La ubicación más desfavorable depende del sistema,
pero corresponde generalmente a puntos cercanos al generador. La poca probabilidad de ocurrencia
de fallas trifásicas hace que en Chile se use como criterio de estabilidad el cortocircuito bifásico a
tierra.

– Interruptores de acción monopolar: La mayoría de las fallas que ocurren en un sistema eléctrico es
de tipo monofásico y se puede eliminar abriendo solamente la fase fallada. Ello exige disponer de
 interruptores de accionamiento monopolar, lo que implica un ligero sobrecosto en esos equipos y en
las protecciones, que se compensa por el hecho de mantener una mayor capacidad de transferencia
durante la etapa que sigue a la eliminación de la falla.

– Reconexión automática: El hecho de que la mayoría de las fallas sea de carácter fugaz y la ventaja
que significa para la estabilidad el volver a la curva de potencia eléctrica prefalla, hace muy atractiva
la reconexión automática. El tiempo de reconexión está normalmente condicionado por la necesidad
de esperar la desionización del espacio en torno al punto de falla, donde se estableció el arco. Este
tiempo de espera es variable con la tensión, condiciones atmosféricas, etc.

– Reactancia del sistema eléctrico: La potencia eléctrica que se puede transmitir es inversamente
proporcional a la reactancia total que une las máquinas. En lo posible habrá entonces que mantener
valores bajos de esas reactancias (líneas fasciculadas, condensadores serie, máquinas con alta razón
de cortocircuito, etc.), condiciones que por otra parte llevan a subir el nivel de cortocircuito. En casos
particulares puede incluso resultar conveniente construir subestaciones intermedias de
seccionalización, que reduzcan la longitud del tramo afectado por la falla.

– Inercia de los generadores: La aceleración de la máquina será menor mientras mayor su inercia
mecánica. Sin embargo, el alto costo ligado normalmente al aumento de la inercia de una máquina y
el poco efecto relativo en la estabilidad, hacen poco atractiva su modificación más allá de los valores
naturales. Durante la operación del sistema es posible elevar la inercia total conectada, poniendo en
servicio más generadores que los estrictamente necesarios. Sin embargo, ello va en desmedro de la
operación más económica, por lo que representa una medida que se emplea con reticencia.

– Tensión interna de los generadores: Aumentar las fuerzas electromotrices implica también
aumentar la estabilidad. En condiciones normales se operará entonces con una fem alta, la que se hará
más alta mientras más fuerte es la transmisión, para así mantener relativamente constante la tensión
en el consumo. Sin embargo, al ocurrir una falla y crecer la corriente del estator, crece también el
flujo desmagnetizante, produciéndose una paulatina reducción de la fem en el entrehierro. La
constante de tiempo es comparativamente alta, de manera que en gran medida es válida la hipótesis de
fem constante.

– Frenado de los generadores: Es posible pensar en mejorar la estabilidad, sometiendo los

generadores a un frenado de acción rápida cuando comienzan a acelerarse. Ello puede hacerse en las
turbinas Pelton, intercalando deflectores de chorro. Para las restantes turbinas se puede recurrir al
procedimiento indirecto de acrecentar la carga, conectando resistencia a los bornes del generador.

5.3. Estabilidad estática o permanente

5.3.1. Introducción

Un SEP está formado por un grupo de elementos electromecánicos interconectados, cuyo movimiento
puede ser representado por ecuaciones diferenciales apropiadas. Cuando las perturbaciones son grandes, el
sistema de ecuaciones es no lineal pero si ocurren cambios pequeños, las ecuaciones pueden ser linealizadas.

En los estudios de estabilidad estática se pretende aclarar la duda de si acaso la operación

cuasiestacionaria de un sistema con n máquinas es estable frente a los cambios en la magnitud de los
consumos. Ya se ha dicho que existe una potencia máxima con la cual se puede operar en forma estable, que
constituye el límite de estabilidad permanente. Este límite depende no sólo de las características constructivas
del sistema, sino también en gran medida de los reguladores de tensión, ya que si el regulador incrementa la
excitación en función de la carga, el ángulo límite se extiende más allá de los valores puramente
constructivos.

Para el análisis se recurre a las mismas ecuaciones diferenciales usadas en los estudios de estabilidad
transitoria, pero admitiendo algunas simplificaciones derivadas de la lentitud y pequeñez de las oscilaciones.
Por ejemplo, la máquina puede ser representada en su forma normal, con la fem permanente E y la reactancia
permanente Xd, ya que para frecuencias de oscilación del orden de 1 a 3 Hz no vale la pena considerar X' d. En
todo caso, se estará en el lado seguro al usar el mayor valor posible de la reactancia.

5.3.2. Generador conectado a una barra infinita

a. Despreciando el efecto de amortiguación: La ecuación de oscilación es

d 2
M  Pm  Pg (5.65)
2
dt
Como las perturbaciones son pequeñas, es posible linealizar la ecuación de la potencia eléctrica P g
alrededor del punto de operación 0, considerando pequeñas variaciones 1 del ángulo, es decir: (t)=0+1(t)
y por lo tanto: Pg ()=Pg (0+1). Desarrollando Pg en serie de Taylor se puede escribir:

 P 
Pg  Pg ( 0 )  1   2 
g 1  Pg
2  ········ (5.66)
 
1  2 
2! 
  0  0

Despreciando los términos de orden igual o superior a 2, dado que 1 es pequeño; se tiene:
d 21
 Pg 
1  0 (5.67)
M  
dt2
   0
Cuya solución es del tipo:
t t
1 (t)  C1e p1  C 2 e p 2 (5.68)

en que C1, C2 son constantes que dependen de las condiciones iniciales y p 1, p2 son las raíces de la
ecuación característica (5.69), siguiente:

2  Pg 
Mp   (5.69)
 0

  0

p1 , p2     P
g
 1
(5.70)
   M
0

Sea
K ) el coeficiente de potencia sincronizante, evaluado en =0. Para la ecuación
(Pg
0
(5.70) tenemos 3 situaciones:

- Si K  0: ambas raíces son imaginarias y el movimiento es oscilatorio no amortiguado. El sistema es

estable considerando el hecho de que se ha despreciado PD.

- Si K  0: ambas raíces son reales, una positiva y una negativa, y el sistema es inestable.

- Si 0 = 90º  K = 0 y el sistema está en el límite de estabilidad ya que las raíces son cero.

De lo anterior se deduce que basta con calcular el coeficiente de potencia sincronizante K para
determinar si el sistema es estable o inestable.
b. Con efecto de Amortiguación: En este caso, la ecuación diferencial, incluyendo D, queda:
d 21
 Pg  d1
 (5.71)
M D   1  0
2
  
dt
dt2  0

Cuya ecuación característica es:

2 D 1  Pg 
 0 (5.72)
p  p    
M M 
0


Las raíces de (5.72) son:

D  D 2 1  Pg 
p,p    
1 2  2M M    (5.73)
2M 0

- Si K  0, el sistema es inestable

- Si K  0 y además (K/M)  (D/2M)2 el sistema es estable.

5.3.3. Estabilidad de un sistema de "n" máquinas

a. Despreciando el efecto de los reguladores y del roce: La ecuación de oscilación de la máquina i

es: d 2 i
Mi  Pmi  Pgi (5.74)
dt2
en que, la potencia eléctrica de la máquina i se puede escribir como:
n
Pgi   Ei E jYij cos(ij   j  i ) (5.75)
j1

Yij incluye la reactancia síncrona de eje directo; es decir, la matriz de admitancia de barras se calcula
referida a las tensiones internas de las máquinas. La ecuación (5.75) se puede escribir:
P  E 2G  n E E Y cos(     )
gi i ii i j ij ij j i (5.76)
j1
ji

Si j = j0+j, en que j0 es el ángulo de régimen permanente y j es un pequeño incremento del
ángulo desde j0; (5.76) se puede desarrollar en serie de Taylor en torno de j0, con j=1,2,.......,n.
Despreciando los términos de orden superior a 1 y reemplazando en (5.74) se puede escribir:
d 2 i
Mi n
d t 2   Kij j  0 (5.77)
j1

donde:

 Pgi 
K ij  
 j  (5.78)
 0
n  P   Pgi 
Kii    gi    (5.79)
j1  j  i 
ji
  0
0
 Para las n máquinas y en forma matricial se puede escribir:

M1 0 · 0
  p 21 K11 K12 · K1n  1 0
 
   
K 
0 M2 0 p 2 K · K 2n 2  
· 21
22 0
 
 2   (5.80)

 · · · ·  ·   · · · ·   ·    · 
0 0 · M p 2  K K · K K  
0
 n   n   n1 n nn   n   
2
o bien:

Mp 2
 K    0 (5.81)

Pre-multiplicando (5.81) por M 1 se obtiene:

 I  p 2
 M
1
K    0 (5.82)

que corresponde al problema de valores propios:

detI A 0 (5.83)

donde A  M K es la matriz de estado del

1
I es la matriz unitaria y  = p2 corresponde
sistema, a los valores propios.

Si  es un real negativo  p es imaginario puro, por lo tanto, la respuesta es oscilatoria y el sistema

es estable.

Si  es un real positivo o un número complejo, una raíz p tendrá parte real positiva y el sistema es
inestable

En las expresiones anteriores se está trabajando con ángulos absolutos en cuyo caso, un valor propio
debe ser cero ya que la matriz A es singular. Por esta razón conviene trabajar con ángulos relativos,
considerando una de las máquinas como referencia.

b. Considerando el efecto del amortiguamiento: Si se agrega el efecto del roce y de los enrollados
amortiguadores, la ecuación (5.81) queda de la forma:

M p 2
 D  p   K    (5.84)
0
donde:

D1 0 · 0 
D2 · 0 
D  0 
(5.85)
 · · · · 
 
0 0 · Dn 

y en que D1, D2, . . ., Dn son los coeficientes de amortiguación de cada una de las máquinas.

Problemas propuestos

5.1. En un sistema eléctrico de potencia, se han establecido las siguientes ecuaciones de P g=Pg(): PgAF=1,8
sin  (pu); PgF=0,5 sin  (pu) y PgFD=1,2 sin  (pu). Determinar el tiempo crítico de despeje de la falla, si P m=1
(pu) y M=2,5x10-4 (pu).
5.2. Dadas las ecuaciones Pg=Pg() y las curvas de oscilación para falla sostenida, de los sistemas 1 y 2
siguientes, determinar:
a. Las áreas acelerantes y desacelerantes para despeje de la falla en 0,4 y 0,25 seg, indicando, en ambos
casos, si el sistema es estable o inestable
b. El tiempo máximo de despeje para que el sistema sea estable

Sistema 1: PgAF=2 sin  (pu); PgF=0,75 sin  (pu) y PgFD=1,5 sin  (pu); Pm=1,0 (pu)

t (seg) 0 0,05 0,10 0,15 0,20 0,25 0,30 0,35 0,40

 (º eléctricos) 15 17 20 25 32 40 50 65 77

Sistema 2: PgAF=2 sin (+15º) (pu); PgF=0,5 sin  (pu) y PgFD=1,5 sin  (pu); Pm=1,0 (pu)

t (seg) 0 0,05 0,10 0,15 0,20 0,25 0,30 0,35 0,40

 (º eléctricos) 30 33,13 42,18 56,19 73,97 94,54 113,9 127,6 134,7

5.3. Un generador que tiene una constante de inercia de 2,5 seg, está suministrando una potencia de 0,8 (pu) a
una barra infinita, a través de una red completamente reactiva cuando se produce una falla que reduce la
máxima potencia posible de salida a 0,5 (pu). La potencia máxima que se puede transmitir antes de la falla es
de 1,5 (pu) y de 1,3 (pu) después del despeje de la falla. Si el despeje se produce a los 0,2 seg: Dibujar
a. La curva de oscilación hasta 1 seg.
b. Las curvas de Pg =Pg() indicando claramente las áreas acelerantes y desacelerantes.

5.4. El sistema de la Figura 5.25, está operando en V 0º

condiciones normales, esto es, el generador está G Línea 1
entregando 125 MVA, factor de potencia 0,8
inductivo a la barra infinita; E'=1,8 (pu); X’ d=0,9
E'  Línea 2
(pu); V=10º (pu), H=5 seg, reactancia de cada línea 1 2
0,4 (pu); cuando ocurre un cortocircuito trifásico en Barra inf.
el medio de la línea inferior. Los interruptores 1 y 2
despejan simultáneamente la falla en un tiempo t 1 tal Figura 5.25
que (t1)=90º. Los interruptores permanecen abiertos
hasta t2 seg, donde (t2)=120º. En t2, los interruptores reconectan la línea y la falla ha desaparecido.
Considerando que todos los datos están en una base común (100 MVA), determinar si en estas condiciones, el
sistema es estable o inestable.

5.5. El sistema de la Figura 5.25, está operando en condiciones normales, esto es, el generador entregando
120 MW a la barra infinita, E'=1,65 (pu), X'd=X1=0,6 (pu); X2=0,3 (pu); X0=0,1 (pu), V=1 0 o (pu), H=5
seg, f=50 Hertz cuando ocurre un cortocircuito monofásico en la línea 2, al lado de la barra del Generador.
Los interruptores 1 y 2 (tripolares) despejan simultáneamente la falla. Los datos de la línea son: X 1=X2=0,3
(pu); X0=0,5 (pu). Todos los parámetros están en tanto por unidad, base común 100 MVA. Determinar el
tiempo crítico de despeje de la falla

5.6.- Un generador de 60 Hertz está suministrando el 60% de la potencia eléctrica máxima a través de una red
completamente reactiva, cuando ocurre una falla que incrementa la reactancia de la red entre el voltaje interno
del generador y la barra infinita en un 400%. Cuando la falla es despejada, la potencia máxima que se puede
suministrar es el 80% del valor máximo original. Determine el ángulo crítico de despeje de la falla.

5.7.- Si el generador del problema anterior tiene una constante de inercia H de 6 [MJ/MVA] y la potencia
mecánica es de 1 (pu), determinar el tiempo crítico de despeje de la falla. ¿A qué velocidad gira la máquina
en ese instante?
5.8.- Un generador síncrono está entregando el 25 % de la potencia máxima que es capaz de entregar a una
barra infinita. Si la potencia eléctrica entregada por el generador a la barra infinita se triplica bruscamente,
determinar el valor máximo alcanzado por el ángulo  durante las oscilaciones alrededor del nuevo punto de
equilibrio. ¿Cuál es el ángulo  de régimen permanente en el nuevo punto de equilibrio?

5.9.- Un generador que está operando a 50 Hz, entrega una potencia de 1 (pu) a una barra infinita a través de
una red completamente reactiva, cuando ocurre una falla que reduce la máxima potencia transferible a 0,4
(pu), mientras que la máxima potencia transferible antes de falla era de 1,8 (pu) y de 1,25 (pu) después que la
falla es despejada. La constante de inercia es de 4 seg. La falla se despeja a los 0,2 seg y el sistema se
reconecta exitosamente 0,15 seg después que la falla ha sido despejada. Utilizando exclusivamente la curva
de oscilación, determinar si el sistema es estable en las condiciones planteadas.

5.10.- En el sistema de la Figura 5.26, V=1,1 Barra inf.

todos los parámetros en pu están en base . j0,4
100 MVA. La potencia compleja en el G S
punto indicado es S&  (80  j40) j0,4
13,8 KV 1 2
MVA ,
cuando ocurre un cortocircuito trifásico F
Xd' =0,25 j0,1
en el punto F (al medio de la línea). La
falla es despejada simultáneamente por 40 MW, F. de P.=1
ambos interruptores. Determine las
ecuaciones de Pg() antes de falla, en Figura 5.26
falla y en falla despejada.

5.11.- Repetir el Problema 5.10, considerando que la falla ocurre al comienzo de la línea; es decir, al lado del
interruptor Nº 1).

5.12.- Repetir el Problema 5.11, considerando que la carga es un reactor de 10 MVAr.

5.13.- Las ecuaciones de Pg=Pg() de un sistema Generador-Barra infinita son: P gAF=2,5 sin ; PgF=PgFD=0. La
potencia que el generador entregaba a la barra infinita en el momento de producirse la falla era de 1,0 (pu). La
falla es despejada a los 0,1 seg y 0,1 seg más tarde se hace una reconexión exitosa. Todos los datos están en
una base común. H=5 seg., f=50 Hz. Determine el valor máximo alcanzado por el ángulo  durante las
oscilaciones.

5.14.- Un generador de 50 Hz con H=5 (MJ/MVA) está conectado a través de un transformador elevador a
una línea de transmisión. En el otro extremo de la línea hay un transformador reductor que une la línea a una
barra infinita. Las reactancias en pu del generador son: X 1=0,3; X2=0,15; X0=0,05, las de los transformadores
X1=X2=X0=0,1 y las de las línea de transmisión son: X 1=X2=0,25 y X0=0,70. Los transformadores están
conectados en delta en el lado de baja tensión y en estrella con el neutro a tierra en el lado de alta tensión. Se
produce una falla monofásica en el lado de alta tensión del transformador conectado al generador, en el
momento que éste está suministrando una potencia de 1,0 (pu). La tensión interna del generador es 1,24 (pu)
y la tensión de la barra infinita es de 1,05 (pu). La falla es aislada simultáneamente por interruptores
monopolares a ambos lados. Estos la aislan en 0,15 seg. y se reconectan 25 ciclos después de la apertura.
Determinar si el sistema permanece estable en estas condiciones.

5.15.- El sistema de la Figura 5.27 está trabajando en las condiciones indicadas. Todos los parámetros en
porcentaje están en base 100 MVA. Determinar el tiempo máximo de despeje de la falla si ocurre un
cortocircuito monofásico a tierra en el punto F, al comienzo de la línea, y ésta tiene interruptores
monopolares.

5.16.- Resuelva el problema 5.15, considerando que los interruptores son tripolares y que se produce
reconexión automática exitosa 0,1 seg después que la falla ha sido despejada.
 X1=X2=40%  V4 =100%
G 1 2 3
 X0=100% 100 MW 4
100 MVA 1
15 KV
75 MVA 2 75 MVA Barra
X'd =X1 =X2
15/115 KV F 115/15 KV Infinita
=30%
X0=15% X1
X1 =X2=X0=5%
=X2=X0=5%
E'=140%
H=1,875 seg

Figura 5.27

5.17.- En el sistema de la Figura 5.28 ocurre un cortocircuito bifásico a tierra en el punto F, al comienzo de
la línea, cuando la potencia transmitida en el punto indicado es S&  (175  j20)MVA . Todos los datos están
en base común 100 MVA.
a. Si los interruptores son tripolares y despejan la falla en forma simultánea, calcular el tiempo de
aclaramiento de la falla (sin reconexión)
b. Si los interruptores reconectan la línea 10º después del despeje, calcular el tiempo de aclaramiento de
la falla en esta nueva condición.

 V=100%
.
G S
X1=X2=35%
200 MVA X0=45% c/u
200 MVA 1 2
X'd =X1 =X2 =25% 0 Barra
X =X =X =20%
F
X0 1 2 infinita
=25%
H=3 seg
Figura 5.28

5.18.- El sistema de la Figura 5.29 está entregando la potencia indicada a la barra infinita, cuando se abre la
fase “a” de la línea en el punto P, próximo a la barra 2. Determinar el máximo valor alcanzado por el ángulo 
al producirse la oscilación del rotor de generador. Todos los datos en % están en base 50 MVA.

 2
G 1
X1=X2=20%

3 (50+j20) MVA 4
X0=50%

X'd =X1 P
=30% V3 =100% Barra
X2=15%
X1 =X2=X0=5% X1 Infinita
X0=10%
=X2=X0=5%

Figura 5.29

5.19.- Determinar si los sistema de las Figuras 5.30 y 5.31 siguientes, son estables o inestables para los
condiciones permanentes indicadas y suponiendo perturbaciones pequeñas. Los datos en porcentaje están en
base común.

G1 =30% E1=12 1=30º

H=3 seg 0%
X
 X=60%
G2
Figura 5.30

X=60% X=40%
X=50% X=50% H=2 seg
E2=140%
2=45º

G1 X=60%
G2

X=30% X=60% X=60%

X=40%
H=3 seg X=50% X=50% H=2 seg
E1=120% E2=140%
1=30º 2=45º
X=50%
X=30%
G3 H=3 seg
E3=100%
3=0º

Figura 5.31

5.20.- En el sistema de la Figura 5.32, ocurre un cortocircuito trifásico en el punto F (mitad de la línea
inferior), cuando el sistema estaba trabajando en las condiciones indicadas. Todos los datos están en base 100
MVA. La falla es despejada por ambos interruptores en forma simultánea y luego el sistema se reconecta.
Determinar:
a. Las ecuaciones de PgAF, PgF y PgFD.
b. La potencia eléctrica que estaba entregando el generador y su ángulo de torque en el momento en que
ocurrió la falla

V=100%
G (40+j10) MVA
X=40% c/u
100 MVA 1 2
15 kV 100 MVA F Barra
100 MVA
115/15 kV Infinita
X'd=30% 15/115 kV
X=5%
X=10%
(15+j5) MVA (25+j10) MVA

Figura 5.32
 BIBLIOGRAFIA

[1] Grainger John J., Stevenson Williams D., Jr. “Elements of Power System Analysis”, Editorial
McGraw-Hill, 1997.

[2] Sanhueza H., Hernán “Análisis de Sistemas Eléctricos de Potencia”, Departamento de Ingeniería
Eléctrica, Universidad de Santiago de Chile, 1990.

[3] Brockering Ch., Walter. “Sistemas Eléctricos de Potencia”, Pontificia Universidad Católica de
Chile, 1982.

[4] Bergen, Arthur R. “Power System Analysis”, Editorial Prentice-Hall, Inc., 1986.

[5] Gonen, Turan “Modern Power System Analysis” Editorial John Wiley & Sons, 1987.

[6] Elgerd, Olle I. “Electric Energy Systems Theory: An Introduction”, Editorial McGraw-Hill, Inc;
1971.

[7] Wadhwa, C. L. “Electric Power Systems”, India, 1985.

[8] Stagg, Glenn W.; El-Abiad, Ahmed H. “Computer Methods in Power System Analysis”, Editorial
McGraw-Hill Book Company, 1968
c