Capitu lo 14 / Est abilidad transitoria 267

PROBLEMAS RESUELTOS
14.1 . Una unidad de generación inyecta una potencia activa de I p.lI. a través de un tran sformador elevador
y dos líneas en paralelo a una barra de potencia infinita de tensión 1 p.lI. Además, se sabe que la
tensión en el lado de alta del transformador se regula a un valor de 1.05 p.lI. Todo el sistema trabaja a
50 Hz.
Se produce en esta situación un fall o trifásico de las líneas que parten del transformador, junto al
embmT<:\do de alta del mismo, que es despejado por las protecciones mediante la apertura de la línea
en cuestión. Adoptando el modelo clásico para la máqu ina, encontrar el tiempo crítico de apertura del
defecto.
Nota: Los parámetros del sistema se indican en la Tabla 14.1. todos bajo una misma base.

Elemento Parám etros asociados

Ge nerado r X~ =0. 1 5 , H=3s
Transformador X cc = 0.1
Líneas XL = 0.6 cada una
Tabla 14.1 Datos de la red del Problema 14.1.

Resolución
El esque ma en p.u. del sistema previo al fallo es el mostrado en la Figura 14.3.

Figura 14.3 Esquema p revio al fallo para el Problema 14.1.

De los datos de partida y planteando las ecuaciones de flujo que llegan desde las líneas al nudo 1, se
tiene:
V, V, j . 1.05
P = ~ Sen(e2 - e,) = 0. 3 sen(e2 - O) = 1

1.05 · 1
= e2 = 16.602 0
}

cos(e, - O) - - Q = 0.021
0.3 , 0.3

Con estas potencias es posible determinar la intensidad asociada a dicho flujo en el nudo 1,

J _ P - Q} _ 1 - 0021}
- (11 ,) ' - - - :I-l!l.-O-'- I l=!.:..'.2.

Se deduce así la te nsión interna del generador:

E" i.82. = 111 + ( Xcc + Xci, + '2

XL) j.J = 1.151 1 28.537
268 Sistemas eléctricos de potencia

La curva de potencia eléctrica suministrada por el generador ideal frente a la evolución del ángulo del
rotor para el período anterior a la falta es:

E'V
P = 1 sen(o _ el) = 2.093 seno
, XL)
( X,c +Xd + T
que se mu estra en la Figura 14.4, correspondiendo el punto de funcionamiento inicial a las coordenadas
(0 0 = 28 .537. P = 1) .

Pm = 1

00 01 Oc

Figura 14.4 Curvas de potencia·ángulo para el Problema 14.1.

Cuando se produce el fallo la potencia eléctrica se hace cero (debido a que la tensión se anula en
el punto de falta y a que se están despreciando las resi~tencias del sistema), manteniéndose la potencia
mecánica en el eje del generador al valor de la potencia eléctrica demandada por el sistema previo al fallo.
es decir, Pm = l . Al despejar el cortocircuito mediante la apertura de una de las líneas, Figura 14.5 , la
nueva curva de potencia eléclrica p!! vale:

+
E' ~ O.6} U I = I l!!.

Figura 14.5 Esquema unifilar tras el despeje de la fa lta para el Problema 14.1.

Aplicando el criterio de las áreas e igualando la positiva con la negativa (ver Figu ra 14.4) se determina
el ángulo crítico de despeje del cortocircuito que garantiza el sincronismo del generador.
Las áreas aceleradoras y desaceleradoras vienen dadas matemáticamente por las siguientes expresio-
nes:

A+ = (o, - oo)P",

con 0111 = 1T ~ o) = 1T ~ P",-)

at-csen = 2.311
(-1.354
rad. Igualando las anteriores expresio nes se obtiene
un áng ulo crítico de valor Oc = 0.845 rad.
 Capitulo 14 1 Estabilidad transitoria 269

La evol ución del ángulo del rotor del generador vie ne gobernada, para cada intervalo de tie mpo. por
su ec uación dinámica. Para el intervalo corre spond iente a la falta resulta:
H .. .. JIT 50·][ _
- .- 0 = Pm - Pe ==> 8 = H ( P", - P,) = ~( I - O) = ~2.36
JIT
y con las co ndic iones ini c iales 80 = 0.-+98 rad y [, = O. se obtiene:

l ,
8 = 0.498 + 2"52.36'-

El tie mpo crítico asociado al {lIlgulo crítico. Oc = 0.845 rad. se obtiene así de esta últ ima ec uación
result ando ser de le = 0.115 s.

14.2. En la red de la Figura 14.6 y par~íll1etros los mostrados en la Tabla 14.2. la balTa 3, nudo de potencia
infinita y ten sión I p.lI., demanda una potencia compleja de I + 0.5) cuando se produce una maniobra
que desconecta la línea que pal1e de dicho nudo. Al cabo de 0 .2 s se desconecta también la carga del
embarrado 2.

Figura 14.6 Esquema eléctrico inicial para el Problema 14.2.

Determinar el ángulo crítico y el tiempo crítico de puesta en servicio de la línea pa ra no perder la

estabilidad dinámica del sistema (en ausencia de la carga en 2).

Elemento Parámetros asociados

Generador x;, =0.16.H =4 s
Transformador Xcc = 0.1
Línea XL =0.1 4
Carga en 2 So = 0.5 + 0.2j
Tabla 14.2 Datos de la red del Problema 14.2.

Re solución
El esquema e n p.U, del sistema previo al fallo es el mostrado en la Figura 14.7,
En el problema es posible d istinguir c uatro períodos:
J. Previo al fallo.
2. Desconexión de la línea eléc trica ,
3. Desconexión de la carga en el nudo 2,
4. Reconexión de la línea.

O.16 j O.lj 2 O.14 j 3

'\liJ\f'--l---,(lOOi'--f--'~C3

Figura 14.7 Es quema previo al fa llo para el Problema 14.2.

270 Sistemas eléctricos de potencia

Para el sistema inicial, de los datos de partida, y tomando origen de fases en ti 3 , es posible determinar
la tensión del embarrado 2,

8 el ) ' 1 - 0.5}
11, =U3 + X ¿jJ3 = 113 + X ¿j ( _ . = 1 lo + 0.14} = 1.079 17.454"
113 - I l() --

y la tensión interna del generador, previa determinación de las intensidades por el mismo:

Se2)* 0 .5 + 0.2j
JI2 = J e2 + )3 = ( -U 2 + J3 = + (1 - O.5} ) = 1.6091 - 22.802"
1.079 1-7.454

La potencia mecánica suministrada por el generador en estas circunstancias, puesto que no existen
pérdidas en ninguno de los elementos del sistema, coincide con la potencia eléctrica que demandan las
cargas:

Al producirse la desconexión de la línea, el esquema unifilar que se obtiene es el que se representa en

la Figura 14.8.

0.16) 0.1 ) 2

~e2

Figura 14.8 Esquema durante la desconexión de la línea del Problema 14.2.

La carga se modela ahora como carga de admitancia constante, obteniendo dicha admitancia de las
condiciones de tensión existentes e n la carga previas al fallo:

s;
~ C2 = ~ = 0.4624 1- 21 .796 =} Z C 2 = 2.008 + 0.803}
V-
2

Durante este período la potencia mecánica se mantiene mientras que la potencia eléctrica demandada
cOlTesponde a la consumida por la carga conectada en 2:

= ~¡(ZC2) ¡2 = :H(Zu) E' )2= 2.008· (0.59)2 = 0.698

P,>
( I(X".1, . .
+ X + ZC2) 1
cc )

Puesto que la potencia eléctrica demandada es menor a la potencia mecánica del generador se pro-
duce una aceleración del mismo, ob[Cniéndose la evolución del ángulo del rotor de la ecuación dinámica
cOlTespondiente, esto es:

H.. .. fJr
- 8 = PI/! - Pe ==> 8 = - ( Pm - Pe) = 31.491
frr H

Inicialmente, el rotor se encuentra en estado de régimen permanente, es decir, la velocidad de cambio

del ángulo del rotor es cero, 8(0) = O, lo cual implica, tras integrar la ecuación dinál1lÍca anterior, una
evolución dada por:
1
8 = 80 + 2:31.49 112
 Capítulo 14 / Estabilidad transitoria 271
El final de esta fase tiene lugar 0.2 segu ndos después de iniciarse. En ese instante el ángulo de l rotor
llega a alcanzar el valor:
1 ?
.I, ="0 + 231.49 1(0.2)- = 1.033 rad

Procede ahora analizar el tercer intervalo de tiempo. en el cual se produce la desconex ión de la carga
del nudo 2. Esto impl ica que el generador queda aislado del sistema, cayendo la potencia eléctri ca a cero.
La ec uac ión de evolución del rotor es ahora

.. frr 50rr
" = - (P", - Pe) = - ( 1.5 - O) = 58.905
H 4
con co ndiciones iniciales obtenid as del período inmediatamente anterior. es dec ir, ángu lo y velocidad
inicial de valores:
"1 = 1.033 8, = 8(02) = 3 1.491·0.2 = 6.298
La dinámica del áng ul o del rolor viene dada pues por la ecuación:

1 ?
1.033 + 6.298/ + :258.905/-
Cuando filialmente se co nccta de nuevo la línea, permaneciendo la carga del embarrado 2 fuera de
se rvic io y suponiendo que el nudo 3 sigue siendo de potenc ia infinita, la palenci a e léctri ca de mandada al
generador vale:
E'V3
Pe = 7(-:----·----c) sen" = 3.349 sen"
X;,+X cc+ XL
Para determinar el tiempo máximo que se puede esperar para la fecone xión de la línea, tiempo crítico,
en aras a evitar la pérdida de es tab ilidad del siste ma, se apli ca e l criterio de las á reas al igual que se
ha venido haciendo has ta ahora. Se trata de igualar e l afea aceleradora, A +. con la desacele radora, A - ,
indicadas en la Figura 14.9 .

Pm = 1.5
0.698

Figura 14.9 Curvas de potencia-ángu lo para las distintas etapas del Problema 14.2

Se tiene. por lanto.

A + = (.I, - .I0)( PIIl - 0.698) + (.Ir - .I,) P",

A- =
1" 8,
(3.349 sen " - P", )d.l

con ,h- = ir - arcse-l1

(-
P",-) = 2.677 rad.
3.349
Igualando ambas áreas se obtiene un ángulo críti co de valor Oc = 1.578 fado
272 Sistemas eléctricos de potencia

Sustituyendo este valor en la ecuación de evolución del ángulo del rotor correspondiente al tercer
período de tiempo analizado, se obtiene el instante de tiempo crítico para el reenganche de la línea:

1
O, = 1.578 = 1.033 + 6.2981 e + 2:58.9051} ==} le = 0 .066 S

A este valor hay que ai'íadirle los 0.2 s que le preceden desde el estado de partida, lo cual hace un
tiempo total de 0. 266 s.

14.3. Sea la red de la Figura 14.10. de 50 Hz, cuyos datos (bajo una misma base y correspondientes al
modelo transitorio para el generador) se indican en la Tabla 14.3. Con la línea L2 fuera de servicio se
sabe que la tensión en el nudo I es U 1 = 1.07 + O.28j. Y que el embarrado 3 es un nudo de potencia
infinita de tensión I p.ll.

0 1
0)
Figura 14.10 Esquema unifilar de la red del Problema 14.3.

Elemento Parámetros asociados

Generador x~ = O.05,H = 5 s
Transformador X cc = 0.2
Línea Ll X = 0 . 15

Tabla 14.3 Datos de la red de la Figura 14.10.

l. Demostrar que el sistema es inestable ante UBa maniobra que desconecta transitoriamente el
generador durante un tiempo de 0 .37 segundos.
2. Con la línea L2 en servicio, y suponiendo en el generador las mismas condiciones previas a
la maniobra del apartado anterior. esto es, pOlencia mecánica y tensión interna. determinar la
reactancia que debe poseer la línea L2 para que el sistema sea estable frente a la misma maniobra
de desconexión del generador durante 0.37 segundos.

Resolución
El esquema con !11i.1gnitudes en p.ll. de la red de partida es el que se muestra en la Figura 14.11 .

(0 0.05) 0) 0.2) 0) 0.15) 8

Figura 14.11 Esquema anterior al fallo para el primer apartado de la red del Problemu 14.3.

Primera cuestión. Es pos ible distinguir tres intervalos de tiempo:

I.a) Antes del fallo: conexión del generador al embarrado de potencia infinita (lÍnea L2 abierta).
 Capitulo 14 / Estabilidad trans itoria 273

J.b) Durante el faJlo: desconexión del grupo generador.

l.c) Posteri or al fallo: reco nexión de l generador a la red,
Previamente a la desconexión del generador se ha de determinar la tensión interna de l mismo, E' ~ .
Puesto que se conocen las tensiones en los nudos 1 (U 1 = 1.07 + O.28j = 1.106 11HM ) Y 3 (U3 = I l".),
es posible determinar las potencias acti va y reacti va inyectadas en 3. esto es,

De esta fo rma, la tens ión interna del generado r, nudo O. vale:

83 ) ' 0.8 - 0.2j

E'I," = V3+XQ3j* - = 1 + (0.05 +0.2 +0.15)j-----,-"-
.- ( V3 l
= 1.08 + 0.32j = 1. 126 116.5().1
y la pote ncia mecánica suministrada por e l mismo es la inyectada en 3 puesto que las resistenc ias se han
desp reciado.
Al desconectarse e l generador, resulta una potencia eléctri ca nula, p;lf = O, por 10 que la ecuación de
osci lación del generador se reduce a:
5 ..
-8=0.8
507r
La integración de esta ec uac ión din ámica, co n la co nd ició n inicial 80 = 16.504° = 0.288 rad, propor-
ciona la siguiente evolución del ángulo del rotor:

l .. l ?
+ 'ir
?
8 = 80 = 0.288 + i5.133t-
En t = 0.37 s se alcanza una ángulo 81 2.008 rad. siendo este instante el de reconexión de l
generador a la red. En el último período (posfallo), la curva de potencia eléctrica coincide con la anterior
al fallo puesto que se retoman las condicio nes iniciales, y viene dada por:

E' V3
?:! = - - s en 8 = 2.8 16sen 8
0.4

En la Figura 14. 12 se muestra la c urva p!! con el punto inicial de eq uilibrio 80 a part ir del cual se
considera la desconex ión del grupo generador, así como el instante, equivalente al ángulo 81, en el que se
vuelve a reCQnectar a la red.

Pm = 0. 8

Figura 14.12 Va riac ión de la potencia co n el áng ul o para el Apartado 1 de la red del Problema 14.3.
274 Sistemas eléctricos de potencia

Si se lÍene en cuenta el criterio de las áreas, el área ace leradora A + no debe superar al área desacclc-
radora A -. Atendiendo a la Figura 14.12 estas áreas vienen dadas por:

A+ =1" 00
( Pm - P<~!)d8 = r" O.8d8 = 1.376
loo
= " 1"-'"
A-
1 o]
( p f! - Pm )d8 =
01
(2 .8 16 sen 8 - 0.8 )do = 0.83 1

Al ser A + > A - se tendría un sistema inestable donde el ángulo del rotor aumentaría sin control
perdiendo el generador su si ncronismo (lógicame nte actuarían las protecciones evitando tal extremo).
Segunda cuestión. También ahora es posible di stinguir tres intervalos de tiempo:
2.a) Antes del fallo: de acuerdo co n el enunciado se mantienen en este caso la potencia mecánica del
generador, Pm = 0.8, Y su tensió n interna, E ' = 1.126. El nudo 3 continúa siendo un nudo de
potencia infinita de tensión V3 = 1, variando la tensión en el nudo l.
Lb) Durante el fallo: se produce de nu evo la desconexión del grupo generador con lo que p(~/f = 0,
exactamenle igual a lo que ocurría e n el primer apartado.
2.e) Posterior al fallo: la reconexión del generador a la red recupera la red inicial , donde la curva de
potencia queda en función de la reactancia de la línea L2 desconocida.
El área aceleradora que se obtiene en esta ocasión vale:

De la integración de la ecuación dinámica, dado que el instante de reconexión del generador es el

mismo qu e en el caso anterior (0.37 segundos después del fallo) se deduce la siguiente relación :

1.. ? d! rrf ?
O, - 00 = - o¡- ( Pm - P, ) = - ¡-. 0. 8 = 1.72 = K
2 2H
qu e llevada a la ecuación inmediatame nte anteli or a esta última, nos define un área ace leradora
A+ = 1.376. Obsérvese que, a pesar de que 00 tiene un valor distinto, aún por encontrar, el área ace-
leradora coincide con la del primer apanado.
Po!' otra parte. la Curva posterior al fall o sigue la expresión :

'1 E' h 1.1 26

Pe1 = - sen o = --seno
0.25 + c-0_.lc-5_._X_ ) XT
( 0. 15 + X
siendo X la reactancia de la línea L2 que se desea determinar.
Si se mantiene el área aceleradora del primer caso y se desea que el.sistema sea estable para al menos
un tiempo crítico de 0.37 s, esto implica que se ha de aume ntar el área desaceleradora, y. por tanto, te ner
una reactancia total X T me nor que la lOtal del caso anterior (OA), para ob tene r así una c urva de potencia
con un máximo más elevado. Matemáticamente esto se traduce en la igualdad de áreas aceleradoras y
desaceleradoras. El área desaceleradora queda como:

A- =
1', o]
( pI! - Pm )d8 = 1"-8"(1.126)
o]
- - sen o - 0.8 do =
XT

=
[
-1.126
cos8 - 0. 88
]"-'0
XT K+oo
donde son incógnitas X T yoo .
De la igualdad de las áreas, A + = A - , y poniendo el ángulo 80 en func ión de X T a partir de las
condiciones iniciales previas al fallo, esto es,
E' V3
- - senoa = sen 00 = 0.8
XT
 Capítulo 14 / Estabilidad transitoria 275

se obtie ne la siguie nte ecuación no lineal:

1.126
1- - 0.8
- X T ) ' ( l +cosK)-0.8senK - 0.8iT + 1.6 arcsen (0. 8
- -XT ) =0
Xr ( I.I~ I . I~

Resolviendo iterativamente se obtie ne:

XT = 0.318 = X = 0.1 23
es dec ir. el valor máximo para la reacta ncia X de la línea L2, e n aras a garanti zar la estabilidad para un
ti empo crítico de recollcxión de 0.37 s. es de 0.123 p.lI.

14.4. En la red de la Figura 14.1 3 y con los parámetros mostrados en la Tabla 14.4. se conocen las tensiones
de los embarrados 3 (nudo de potencia infinita) y 2:

U3 = 1.0 + O.Oj
U2 = 1.01 + 0.2j

sie ndo la potencia demandada por la carga local co nectada en el nudo 2 de I + OAj.

Figura 14.13 Esquema eléctrico inicial para el Problema 14.4.

Se picle:
l. Demostrar que e l sistema es inestable si se desconecta instantáneamente locla la carga de la barra
2, quedando el trafo y la línea conectad os.
2. Si en la situación de partida se desco necta permane ntemente la línea. ¿e s e stable e l sistema?

Elemento Parámetros asociados

Generador X;, = 0.35
Transformador XC(' = 0.15
Línea XL = 0.4

Tabla 14.4 Datos de la red del Problema 14.4.

Resolución
El esq uema unifi lar previo a la perturbació n corresponde al de la Figura 14. 14 .

O. 35j O.15j 2 3
113 = l l~

Figura 14.14 Esqu ema previo al fallo para el Problema 14.4.

276 Sistemas eléctricos de potencia

Con las tensiones de los nudos 2 y 3 es posible determinar la intensidad por la línea J L

"J L -_ (U 2 -'U 3
. )_- O'
.) _ OO?"_O
. -)i - .)' OI 1 - .,. 86,,
XLi ----
y la inte nsi dad qu e demanda la carga e n 2:

12 = (S,)*
U: = 1.046 1- 10.601

Así. la inte nsidad por el tran sformador corresponde a la sum a rl e las anteriores:

Finalmente, se determina la ten sión interna del ge nerador

y la potencia eléctrica sumini strada :

Pe = P", = ~/(E' L& · 1~) = 1.5

Primera cuestión. En esta süuación, al desco nectar repentiname nte la carga conectada al nud o 2. la curva
de án gulo· potenc ia eléctrica asociada val e:

dI E ' V)
P, = , sen o = 1.641 sen o
(X cc + X d + X L!
Siguiendo COIll O en ejemplos anteriores el criterio de las áreas, para que el sistema fuera estable se
habría de cumplir que e l área acelerad ora no superase al área desaceleradora máxima (ver Figura 14.15).

Pm = 1.5

Figura 14.15 Curva de potencia-á ngulo tras la pérdida de carga en el Prob lema 14.4

Estas áreas valen:

A += r
Joo
O
' (P", -I.64l sen 8) d o

A- = r
Jo,
O
'( 1.64/ Sen o- P",)do

con 8 1 = arc sen ( P", / 1.641) = 1.1 53 rad, y 82 = 7r - 81 = 1.988 rad.

Tras sustitu ir va lores en es tas ec uaciones se deduce que las protecciones del generado r, para evi tar
la ace leración incontrolada del mismo, dispararían. ya que e l área acelerado ra, A + = 0.085, supera a la
desace le radora máx ima, A + = 0.078.
 Capitulo 14 1 Estabilidad transito ria 277

Segunda cuestión. Se parte de nuevo de la si tuac ión inicial. y. como es habitual pa ra este tipo de estudios,
se supone, para el período de falta . un modelo de admitancia constante para la carga. Así. la carga en 2
viene definida por una impedancia que se obtiene con las condiciones sobre la mi sma prev ias al fallo:

(1.01 2 + 0.2 2 )
0.913879 + 0.36555j
( 1 - OAj )

El c ircuito equivalente desde que se produce la pérdida permanente de la línea se indi ca en la Figura
14.16. La intens idad que circula por e l mi smo va le :

,
E' l! =
o
1=1.1 704
J=
(Xdj + X,d + Z,)

0.35) 0.15) 2

Figura 14.16 Modelo del sistema para el Problema 14.4. tras la perdida de la linea

La única potencia activa demandada en el sistema corresponde a la consumida por la parte real de la
impedancia q ue modela la carga e n 2. y corre sponde pues a la que suminis tra el generador. S u c uan tía es
de:

De este valor se deduce fácil me nte la inestabilidad del sistema. ya que la potencia mecá nica del
ge ne rador. PIII = 1.5. supera la eléctrica demandada. P;'j = 1.258 (ésta permanece constante debido a l
modelo clásico adoptado para el generador y a l model o de admitancia cons tante para la carga). con lo q ue
se produce un área aceleradora que crece indefinidamen te. En la Figura 14.17 queda de manifiesto dicha
si tu ación.

Pm = 1.5
I

Figura 14.17 Diag ramas de potencia -ángulo tras la pérdida de la línea en el Problema 14.4.

14.5. En la red de la Figura 14.18 y con los panímetros mostrados en la Tabla 14.5, el embarrado 5, de
potencia de cortocircuito infinita y tensió n I en p.lI .. cons ume una potencia S5 = 1 + O.2j p.u.
Determinar:
l . Tensión interna del generador.
2. El ángulo Clítico de eliminació n de un fallo trifásico rígido en el embarrado 3, fallo que es
eliminado mediante la actuación simultánea de los interruptores B 1 Y B2.
278 Sistemas eléctricos de potencia

Figura 14.18 Esquema eléctrico inicial para el Problem a 14.5.

Elemento Parámetros asociados

Generador de rotor liso x;, = X;, = O.IS .
Transformador 1-2 X cc = 0.05
Transformador 4-5 X cc = O.OS
Línea 2-4 XLI = 0.3
Línea 2-3 XL2 = O.1
Línea 3-4 XL 3 = 0.2

Tabla 14.5 Datos de la red del Problema 14.5.

Resolución
Atendiendo a los dmos asociados a cada ele me nto, el esquema unifilar previo al fall o con'esponde al
indicado en la Figura 14, 19,

0.3)

@ 0. 15) (0 0.05)
.-----~~~--~
0.05 ) 0
,
cp E'~

0) S = 1 + 0 .2)

Figura 14.19 Esquema previo al fallo para el Problema 14.5.

La reactancia equi vale nte ent re el nudo O y e l e mbalTado 5 es:

, X LI (X L2 + X u ) )
X = X,
,
+ X,.,. + ( XLI + XL2 + XL3
+ X,.,. = 04 .
.' J

De la pote nc ia compleja consumida e n 5 y de la ten sión de dicho nudo se ded uce la tensión interna
del ge nerador así como el resto de las condiciones previas al fallo:

p = E' v, se n (80 - es) =

X
.E..
DA
sen (80 - O) = I 1= 80 = 20.323°
E' . V V2
E' I { E' = 1.152
Q= ' cos(80 - 05) - í = - cos(80 - O) - - = 0.2
X X OA 0.4
Durante e l cortocircuito se tie ne el primer esque ma de la Fig ura 14.20. donde se muestra n ade más
los distintos un ifila rcs simpli ficados que se van obteniendo cua ndo se reali zan conversiones estrella-
triá ngulo. La primera equivalencia consiste en pasar la estrell a constit ui da por los extremos O. 4 y 3 a
 Capitulo 14 / Estabilidad transitoria 279

un triángulo (desaparece el nud o 2). con lo que se obtiene el segundo esqu ema de la Figura 14.20, y
la segunda tran sformación es pasar de nuevo a un triángu lo la estre lla form ada por los nudos 0, 3 Y 5
(desaparece el nudo 4). obte ni éndose el últ imo circu ito monofásico de la Figura 14 .20. En es te último
esquema se ti e ne una curva de potencia-ángulo dada por:

d! E/V)
P, = ----;:-::sen(8 - e5) = 0.755 sen.
1.)2)

0.3)

0.15) 0) 0.( 5) 005) 0

0.1)

E'~

1.1 ) 0.05) 0
E'~ 0.55)
0.2)

- -
0) - -

íf)\
\V
1.525)
'00()'
0
o' 1 l'!.

p,----.--- ~ w';¡
r -·" ,'- "'"'; ~ I
o
Figura 14.20 Esquem a del fa llo para el Proble m a 14.5 y reducción del circuito.

La curva de potencia-ángulo Iras el despeje de la fa lla se determina del ci rcui to monofásico corres-
pondiente a este período y que se indica en la Figura 14.21. Así. resulta:

pf E ' Vs
P,. = sen(8 - e5) = 2.095 sen 8
X;J+ X{,("+ XL I +X cc
Todas la curvas de potencia. las ante riores al fallo. durr.nte la falta y la poste rior al despeje de la misma
se muestran e n la Figura J 4.22, donde se ap lica de lluevo el criterio de líJS áreas para determinar el ángulo
crítico de apertura de los interruptores .
Los á ngu los 01 y 01 valen:

P", )
8, = aresen ( 2.095 = 0.498 rad

82 =" - 8, = 2.644 rad

280 Sistemas electricos de pot encia

0. 15) (í) 0.05) 0) O.3j 0 0 .05 ) 0

-E'~

Figura 14.2 1 Esquem a tras el despeje de la f alta en el Problema 14.5.

El área ace leradora y desacelerad ora vienen dadas por:

A+ =
ro, ( Pm -
J.. 0.755sen8) d8

A- =
1"
0,
(2.095 5en8- Pm )d8

De la igualdad entre ambas se obtiene un ángulo crítico de valor 1.766 rad. equ ivalentes a 10 1. 19° . La
obtenc ión del tie mpo crítico requeriría en este caso la integración numérica de la ec uación de oscilac ión.
al no ser la ace leración constante.

PI/! = I

Figura 14.22 Cu rvas de potencia-ángulo para las d istinta s et apas del Problema 14. 5,

14.6. En el sistema de la Figura 14.13 y con los parámetros most rados en hI, Tabla 14.6. se produce un fallo
en la línea L3, cerca del nudo 3. que se desp'eja en 0 .1 segundos abriendo las protecciones de dicha
línea.

0 0)
L2

~I
GI G2

o L4
L3
0
Fig u ra 14.23 Esquema unifi lar del Problema 14.6.

Tomando una potencia base de 100 MVA y conociendo el estado de la red previo al fallo, Ta-
bla 14.7, determinar numéricamente la evolución de los ángu los de los rotores de los generadores.
Nota: modelar las cargas como impedancias constantes.
 Capitulo 14 / Estabilidad transitoria 281

Elemento Parámetros asociados

Generador I X~ = 0.28 p.u .. S = 50 MVA. H = 5 p.u.
Generador 2 X~ = 0.25 p.u., S = 120 MVA, H = 4 p.ll.
Transformador 1 X cc = 0.08 p.U .. S = 50 MVA
Transformador 2 X cc = 0.07 p.u., S = 120 MVA
Lín ea LI XLI = 0.08 (Ss = 100 MVA)
Línea L2 X L2 = 0.06 (S8 = 100 M VA)
Línea L3 X L3 = 0.06 (S8 = 100 MVA)
Línea L4 XL 4 = 0 . 13 (S8 = 100 MVA)

Tabla 14.6 Datos de la red del Problema 14.6.

Nudo Tensión Angula Potencia Sc Potencia Sc

(p.u.) (Grados) consumida ( MVA) generada ( MVA)
1.03 O O+ Oj 30 + 23.l j
2 1.02 - 0.5 80 + 40j 100 + 37. 8}
3 1.018 - 1.0 50 + 20j O+Oj
Tabla 14.7 Estado de la red del Problema 14.6 previo al fallo.

Resolución
Pasando el sistema a p.U. en la base de pote ncias indi cada. se obtie ne el esque ma unifi lar mostrado en la
Figura 14.24. En é l se han mode lado las cargas como admitan cias constantes de valor:

lJC2 = ( :¡
S ) ' = 0.768 - 0. 384j

¡icJ ,)*
= ( S~; = 0.482 - 0.193j

O.266 j

O.1 3j

Figura 14.24 Esquem a unifilar de l Problema 14.6 previo al fallo.

De las condic iones iniciales se determin a el estado inic ial de los generadores, fuerza electromotri z
tras la reactancia transitoria y án gulo inicial:

E: lb. = UI +(XCI +XTl )j(~II)' = 1.03 ~ +(0. 72j )(0.291 - 0.233j) = 1. 209 19.95

E; liL = U2 + (XC2 + Xnlj(~n * = 1.021 -05 + (0.266j )(1.048 1- 21.2) = 1.15211345

282 Sistemas eléctricos de potencia

Las potencias mecánicas en cada generador coinciden en régimen permanente con las eléctricas su-
ministradas por los mismos:
P ml = 0.3 Pm 2 = I

Transformaodo las fuentes de tensión en fuentes de intensidad, el esquema de red previo al rallo es el
mostrado en la Figura 14.25. La malriz de admitancias que lo representa resulta:

-21.58j 7.69j 12.5 j ]

IYa! J = [ 7 .69j 0 .769 - 45.15j 33 .32j
12.5 j 33.32j 0.482 - 46.026j

(í) 0.06)
J1

E; 181 0.72} E; 101

0.266)
0.72) 0.266}

0.13) ~ C2

Figura 14.25 Esquema unifilar del Problema 14.6 previo al fallo usando fuentes de intensidad.

Cuando se produce la falta, la tensión del nudo 3 se hace cero con lo que la matriz de admitancias del
sistema resulta:
-21.58j 7.69j 12.5 j ]
[Yd¡ l = 7 .69j 0.769 - 45 . 15j ..,., .,? .
).) . .)-}
[ 12.5j 33.32j Yoe
donde Y00 representa una impedancia nula a tierra.
El término Y00 que apararece en la matriz no se tiene e n cuenta. puesto que la tensión del nudo 3 es
nul a
7.69j 12.5j
0.769 - 45.15j 33.32j
33.32) Yoo
con lo que realmente se utilizaría el sistema de dos por dos:

7 .69 j ' ] [ 111 ]

0.769 - 45. 15¡ 112

Finalmente, tras despejar la falta mediante la actuac ión de las protecciones, que en este caso se trad uce
en la apertura de la línea L3 , resulta la matriz de adtuitancias :

- 2 1.58} 7.69j 12.S j ]

IY"f 1= 7.69j 0769 - 28.5j 16.66j
[ 12.5j 16. 66j 0.482 - 29.35j

Considerando invaria ntes las potencias mecánicas y las tensiones internas que representan a los gene-
radores, el sistema de ecuaciones y el proceso iterativo para resolverlas sería:
l . Obtener las intensidades inyec t ad<~s en los nudos utilizando la matriz de admüanc ias adecm!da que
modela el sistema según el instante de tiempo de integración.

[Jl = rlJl 1111 (14.10)

 Capitu lo 14 / Estabilidad t ran sito ri a 283

2. Determ inar las potencias e léctricas su ministradas por cada generador.

Pd = ~¡(é:'IJI)* = ~¡( 1. 209 ~ . J[)*

Pe' = m([~J2)* = ~¡(1.I521 8, . J2)*

3. Resolve r las ecuaciones dinámicas de cada generador.

4. Incrementar el interva lo de integración, {JI = f v- I + 61 , e ir al paso l.

Comenzamos aplicando el algoritmo con las condic iones iniciales 8? = 0 .174, o~ = 0.221 rad y
8? = 0, ¿g = o. Se co nsidera qu e el instante 1=0 es aque l en el que se produce el fallo y se toma un paso
de integ racion 6.1 = 0.0 I s.

Primera Iteración:
1. Util izando la matri z de ad mitanc ias que caracteri za el sistema durante el fallo, ?JdJ, resultan las
intens idades y te nsiones:
I
J[ = 0.261-1.492}
.
Jl = 0.859 - 3.779}
U: = 0.117 +0 .022} ul = 0.113 + 0.023}
2. Las potencias eléctri cas sum inistradas por cada generador, aplicando las Ecuaciones ( 14. I 1) son:

3. Se procede ahora a in teg rar las ecuaciones dinámicas (14. 11 ), para lo cua l se ap lica en este ca-
so el método de Euler y las condiciones in iciales indicadas an te riOlTllente, obteniéndose para el
generador 1:

01.. I
.
= --;¡-::-(0.3
50· 1f
-.)
.o
'1
+ 1.5 79 · 10- ') = 18.95

o: = 18.95 · I'>t + D[ = 0.189

I
o: =0. 189I'>t+o?=0. 175

y para el generador 2
.. I = --¡:s-(I
02 SO· 1f - 0.0 12) = 32.33
. .o
oi = 32.33 . I'>t + 02 = 0.323
oj = 0.3231'>t + og = 0.222

4 . Finalmente se pasa a la segunda iteración, 1 2 = 1

1
+ 6.1 = 0.02, volv iendo de nuevo al paso lo

Segunda Iteración:
1. Puesto que aún no se ha alcanzado el instante de despej e de la falta (t=O. I s) se continú a utili zando
la misma matriz de admitanci as durante el fallo, lddf. resultando. para esta segu nda iteración las
intensidades y tensiones:

J; = 0.263 - 1.492} J~ = 0.865 - 3.777}

U; = 0.117 + 0.022} 1l~ =0.1 13 +0.023}
284 Sistemas el éct rico s de potencia

2. Con los nuevos valores de los ángulos e intensidades, las potencias e léclricas suministradas por
cada generador valen ahora:

P}I = -1.619· 1O- 3

3. De la integració n de las ecuaciones dinámicas resulta:

.. ') 50· J[ 3
81 - = - -(0.3 + 1.619·10- ·) =
2.5
8·¡ = 18.95 ·1'-.1 + 8·¡ = 0.379
= 0.379/'-.1 + 8¡ = 0. 178

I
8¡

y para el generador 2

.. •- = -;¡:s(l
8, 50·rr - 0.0 12) = 32.33

s"i = 32.33 . 1'-.1 + 81 = 0.647

, 1
8, = 0.6471'-.1 + 8, = 0.227
4. Se incremema el tiempo de integrac ión. f3 = [2 + ~l = 0.03 Y comenzaríamos la tercera iteración.
Sigui endo este proceso. se llega al instante t::::O.1 s. correspondiendo a la décima iteración. con los
valores:
8:°
= 0.269 = 1.902 8:°
8iO = 0.382 ~ IO
02
. . , T'"
= .J. _ .J

En este insta nte se produce el de speje de la falta, con lo que el primer paso de la sigui ente iteración ha de
uti lizar la llueva matriz de admitancias }J flf.
Tras co mpl etar la iteración 11. de forma análoga a como se ha hecho anteri orm ente, res ulta:

8¡ 1 = 0.288 8¡ =
1 1.948
8i 1 = 0.415 8i l = 3.207

En las Figuras 14.26 y 14.27 se muestran la evolució n temporal de los ángulos y ve locidades de los
generadores del siste ma. Debido a que no se ha te nido en cue nta ning ún factor de amortiguamjento el
resultado es una oscilación permanente de los generadores; rea lmente ésta acabaría anulándose debido a
los pares amortiguadores que en realidad ex isten.
 Capitu lo 14 / Estabil idad transitoria 285

4 ,-----------------~

Ángulo del generador 1

~ An9UIO del generador 2

3 ~--~-~~_f------------~~---------

2 f-+--?'--------------------""'~:$"=-._-

~ I Velocidad del generador 1

~ Velocidad del generador 2
~H~--------~~~-----------------

o
O 0.2 04 0.6 0.8 1.2

Figura 14.26 Evolución de los ángulos y veloc idades de los generadores del Problema 14.6.

315

314 i-------------------------------------------------------~

313 ~~ -------------------------------------------------------
Velocidad relativa del generado r 2
Velocidad relativa del generador 1

312 t- - '\------~l>-----------~/-------
311 +----'\b-"/'-----',~
-----"-~_'----""
_ _ _ _:7"/'_-~
-------
310 ~-------,--------,_------_.--------._--------------_.

O 0.2 0.4 0.6 0.8 1.2

Figura 14.27 Ev o lución de las velocidades, relativas a la de sincronismo, de los g en erado res del
Problema 14.6.