ÍNDICE

1. DESCRIPCIÓN DEL SISTEMA........................................................................................3

1.2. PLANTEAMIENTO DEL PROBLEMA....................................................................4
1.3. OBJETIVOS................................................................................................................4
2. MATRIZ DE ADMITANCIA ANTES DE LA FALLA....................................................4
3. MATRIZ REDUCIDA DURANTE LA FALLA...............................................................7
4. MATRIZ REDUCIDA DESPUÉS DE LA FALLA.........................................................10
4.1. FEM´s INTERNAS DE LOS GENERADORES.......................................................13
5. ECUACIONES DE OSCILACIÓN DURANTE LA FALLA.........................................14
5.1 GENERADOR 1........................................................................................................14
5.2. GENERADOR 2........................................................................................................14
6. ECUACIONES DE OSCILACIÓN DESPUÉS DE LA FALLA.....................................15
6.1. GENERADOR 1........................................................................................................15
6.2. GENERADOR 2........................................................................................................16
7. TABLA Δ-T MÉTODO PUNTO POR PUNTO DE LOS GENERADORES.................16
7.1. GENERADOR 1........................................................................................................16
7.2. GENERADOR 2........................................................................................................18
8. TABLA Δ-T MÉTODO EULER MEJORADO DE LOS GENERADORES..................19
8.1. GENERADOR 1........................................................................................................19
8.2. GENERADOR 2........................................................................................................22
9. ANÁLISIS DE RESULTADOS.......................................................................................24

1. DESCRIPCIÓN DEL SISTEMA

1
El sistema consta de 5 barras: Tres de generación y dos de carga. En condiciones de estado
estable los generadores G1 Y G3, están en las barras 1 y 2, respectivamente, ambos
generadores entregan potencia a la barra 3, esta última puede considerarse que representa a
un sistema infinito.
Los datos de los generadores son
- Generador 1: 400(MVA); 20(KV); X ' d=0.067( p .u .); H=11.2 (MJ/MVA)
- Generador 2: 250(MVA); 18(KV); X ' d=0.1( p . u .); H=8.0 (MJ/MVA)
Se desprecian las resistencias de los generadores, pero se consideran las resistencias de las
líneas y admitancias shunt de las líneas.
Los datos de las líneas y transformadores están dados en la tabla1.
Las barras 4 y 5 son de carga, las potencias de éstas se dan en la tabla 2 y, consisten de
cargas constantes y motores de inducción. Para el presente estudio, los motores se
representan con cargas estáticas.

Tabla 1.
Barra a Barra Z – Serie Y – Paralelo
R X B
Transformador 1 – 4 - 0.022 -
Transformador 2 – 5 - 0.04 -
Línea 3 – 4 0.007 0.04 0.082
Línea 3 – 5 (1) 0.008 0.047 0.098
Línea 3 – 5 (2) 0.008 0.047 0.098
Línea 4 – 5 0.018 0.110 0.2226

Tabla [Link] de líneas y transformadores.

Barra Voltaje Generación Carga

2
 P Q P Q
1 1.03 ∠ 8.88 ° 3.500 0.71 - -
2
2 1.02 ∠ 6.88° 1.850 0.29 - -
8
3 1 ∠0° - - - -
4 1.018 ∠ 4.68 ° - - 1.00 0.44
5 1.011∠2.27 ° - - 0.50 0.16
Tabla 2. Datos de barra y valores de flujos de potencia anteriores a la falla.
Todos los valores están en por unidad sobre la base de 230 KV y 100 MVA.

Fig. 1.

2. PLANTEAMIENTO DEL PROBLEMA

Se produce un cortocircuito trifásico transitorio a tierra en medio de la línea 4-3. Se conoce
que la línea 4-5 no tiene capacidad suficiente para transportar toda la potencia generada por
el generador G1, por lo que es necesario evitar la salida prolongada de la línea 4-3.

3. OBJETIVOS
3.1. Objetivo General

Estudiar el efecto de una falla trifásica transitoria a tierra en la línea 4-3.

3.2. Objetivo Específico

Determinar el tiempo crítico de desconexión, el correspondiente ángulo crítico y la

curva δ -t.

3
4. ANÁLISIS DE LAS OSCILACIONES

MATRIZ DE ADMITANCIA ANTES DE LA FALLA

-Admitancias de carga

P 4− j Q 4 1,00− j 0,44
Y L4 = 2
= 2
=0,9649− j 0,4246 [ p .u . ]
|V 4| 1,018

P5 − j Q 5 0,5− j 0,16
Y L5= 2
= 2
=0,4892− j 0,1565 [ p .u . ]
|V 5| 1,011

-Matriz de admitacia

Fig. 2.

Generador 1; x ´ d = 0.067 (p.u.)

Generador 2; x ´ d= 0.10 (p.u)

4
 Datos de línea y transformadores
Z serie Y paralelo
Barra a barra R X B
Transformador 1-4 - 0,022
Transformador 2-5 - 0,040
Linea 3-4 0,007 0,040 0,082
Linea 3-5(1) 0,008 0,047 0,098
Linea 3-5(2) 0,008 0,047 0,098
Linea 4-5 0,018 0,110 0,226
Tabla 1

1
Y 11= =− j11,2360
j 0,067+ j0,022

Y 12=0

Y 13 =0

1
Y 14 = =− j 11,2360
j 0,067+ j 0,022

Y 15 =0

Y 21 =0

1
Y 22= =− j7,1429
j 0,1+ j 0,04

Y 23 =0

Y 24=0

1
Y 25 = =− j 7,1429
j 0,1+ j 0,04

Y 31=0

Y 32=0

1 1 1 j0,082 j 0,098 j 0,098

Y 33 = + + + + + =11,2842− j65,4732
0,007+ j0,040 0,008+ j 0,047 0,008+ j 0,047 2 2 2

1
Y 34= =4,2450− j 24,2571
0,007+ j 0,040

1 1
Y 35= + =7,0392− j 41,3550
0,008+ j0,047 0,008+ j 0,047

5
 Y 41=Y 14 =− j 11,2360

Y 42=Y 24 =0

Y 43=Y 34=4,2450− j 24,2571

1 1 1 j 0,082 j 0,226
Y 44 = + + + + + 0,9649− j 0,4246=6,6587− j 44,61
j 0,067 + j 0,022 0,007+ j0,040 0,018+ j 0,110 2 2

1
Y 45= =1,4488− j 8,8538
0,018+ j0,11

Y 51=Y 15=0

Y 52 =Y 25=− j7,1429

Y 53 =Y 35=7,0392− j 41,3550

Y 54=Y 45=1,4488− j 8,8538

1 1 1 1 j 0,226 j 0,098 j 0,098

Y 55 = + + + + + + + 0,4892− j0,1
j 0,1+ j 0,04 0,018+ j 0,11 0,008+ j 0,047 0,008+ j0,047 2 2 2

Y BUS =¿

[ ]
− j11,2360 0 0 j 11,2360 0
0 − j 7,1429 0 0 j 7,1429
¿ 0 0 11,2842− j 65,4732 −4,2450+ j 24,2571 −7,0392+ j 41,3550
j11,2360 0 −4,2450+ j 24,2571 6,6587− j 44,6175 −1,4488+ j8,8538
0 j 7,1429 −7,0392+ j 41,3550 −1,4488+ j8,8538 8,9772− j57,2971

MATRIZ REDUCIDA DURANTE LA FALLA

Para hallar la matriz correspondiente se añade un nodo auxiliar 6 donde ocurre la falla
trifásica.

6
 Fig. 3

1
Y 11= =− j11,2360
j 0,067+ j0,022

Y 12=0

Y 13 =0

1
Y 14 = =− j 11,2360
j 0,067+ j 0,022

Y 15 =0

Y 16 =0

Y 21 =0

1
Y 22= =− j7,1429
j 0,1+ j 0,04

Y 23 =0

Y 24=0

1
Y 25 = =− j 7,1429
j 0,1+ j 0,04

Y 26=0

Y 31=0

7
 Y 32=0

1 1 1 j 0,082 j 0,098 j0,098

Y 33 = + + + + + =15,5291− j 89,7508
0,007+ j 0,040 0,008+ j0,047 0,008+ j 0,047 4 2 2
2

Y 34 =0

1 1
Y 35= + =7,0392− j 41,3550
0,008+ j0,047 0,008+ j 0,047

1
Y 36 = =8,4900− j 48,5143
0,007+ j 0,040
2

Y 41=Y 14 =− j 11,2360

Y 42=Y 24 =0

Y 43=Y 34=0

1 1 1 j 0,082 j 0,226
Y 44 = + + + + +0,9649− j 0,4246=10,9037− j 68,
j 0,067 + j 0,022 0,007+ j 0,040 0,018+ j0,110 4 2
2

1
Y 45= =1,4488− j 8,8538
0,018+ j0,11

1
Y 46 = =8,4900− j 48,5143
0,007+ j0,040
2

Y 51=Y 15=0

Y 52 =Y 25=− j7,1429

Y 53 =Y 35=7,0392− j 41,3550

Y 54=Y 45=1,4488− j 8,8538

1 1 1 1 j 0,226 j 0,098 j 0,098

Y 55 = + + + + + + + 0,4892− j0,1
j 0,1+ j 0,04 0,018+ j 0,11 0,008+ j 0,047 0,008+ j0,047 2 2 2

Y 56 =0

Y 61=Y 16=0

Y 62=Y 26 =0

8
 Y 63 =Y 36=8,4900− j 48,5143

Y 64=Y 46=8,4900− j 48,5143

Y 65 =Y 56=0

1 1
Y 66 = + =16,9800− j 97,0285
0,007+ j 0,040 0,007+ j 0,040
2 2

Se tiene la siguiente matriz de admitancia Y BUS =¿

[
− j 11,2360 0 0 j 11,2360 0 0
0 − j7,1429 0 0 j 7,1429 0
0 0 15,5291− j89,7508 0 −7,0392+ j 41,3550 −8,4900+ j 48,5
j 11,2360 0 0 10,9037− j 68,8951 −1,4488+ j8,8538 −8,4900+ j 48,5
0 j7,1429 −7,0392+ j 41,3550 −1,4488+ j8,8538 8,9772− j 57,2971 0
0 0 −8,4900+ j 48,5143 −8,4900+ j 48,5143 0 16,9800− j 97,0

[ ]
− j 11,2360 0 0 j 11,2360 0
0 − j7,1429 0 0 j7,1429
0 0 15,5291− j 89,7508 0 −7,0392+ j 41,3550
j 11,2360 0 0 10,9037− j68,8951 −1,4488+ j 8,8538
0 j7,1429 −7,0392+ j 41,3550 −1,4488+ j 8,8538 8,9772− j 57,2971

Para hallar la matriz reducida se tiene la siguiente ecuación:

−1
Y RED=[ Y nn ] −[ Y nr ]∗[ Y rr ] ∗[ Y rn ]

Realizando las operaciones matriciales se tiene:

Y RED=¿

9
 [ ]
0,2882− j9,4120 0,0271+ j 0,1796 −0,0199+ j 1,0664
0,0271+ j 0,1796 0,1387− j6,2561 −0,0708+ j5,2710
−0,0199+ j 1,0664 −0,0708+ j5,2710 9,9248− j 59,3032

MATRIZ REDUCIDA DESPUÉS DE LA FALLA

Después de ocurrida la falla se aísla la línea 3-4.

Fig. 4.

1
Y 11= =− j11,2360
j 0,067+ j0,022

Y 12=0

Y 13 =0

1
Y 14 = =− j 11,2360
j 0,067+ j 0,022

Y 15 =0

Y 21 =0

1
Y 22= =− j7,1429
j 0,1+ j 0,04

Y 23 =0

Y 24 =0

10
 1
Y 25 = =− j 7,1429
j 0,1+ j 0,04

Y 31=0

Y 32=0

1 1 j 0,098 j 0,098
Y 33 = + + + =7,0392− j 41,2570
0,008+ j0,047 0,008+ j 0,047 2 2

Y 34=0

1 1
Y 35= + =7,0392− j 41,3550
0,008+ j0,047 0,008+ j 0,047

Y 41=Y 14 =− j 11,2360

Y 42=Y 24 =0

Y 43=Y 34=0

1 1 j 0,226
Y 44 = + + + 0,9649− j0,4246=2,4137− j 20,4014
j0,067 + j 0,022 0,018+ j0,110 2

1
Y 45= =1,4488− j 8,8538
0,018+ j0,11

Y 51=Y 15 =0

Y 52 =Y 25=− j7,1429

Y 53 =Y 35=7,0392− j 41,3550

Y 54=Y 45=1,4488− j 8,8538

1 1 1 1 j 0,226 j 0,098 j 0,098

Y 55 = + + + + + + + 0,4892− j0,1
j 0,1+ j 0,04 0,018+ j 0,11 0,008+ j 0,047 0,008+ j0,047 2 2 2

[ ]
− j 11,2360 0 0 j11,2360 0
0 − j7,1429 0 0 j 7,1429
0 0 7,0392− j 41,2570 0 −7,0392+ j 41,3550
j 11,2360 0 0 2,4137− j 20,4014 −1,4488+ j 8,8538
0 j7,1429 −7,0392+ j 41,3550 −1,4488+ j 8,8538 8,9772− j 57,2971

11
 Para hallar la matriz reducida se usa la siguiente ecuación:
−1
Y RED=[ Y nn ] −[ Y nr ]∗[ Y rr ] ∗[ Y rn ]

[ ]
0,7500− j 4,6889 0,0694+ j 0,6443 −0,2333+ j 3,7988
¿ 0,0694+ j 0,6443 0,1426− j 6,2104 −0,0936+ j 5,5396
−0,2333+ j 3,7988 −0,0936+ j5,5396 1,0382− j 9,2769

FEM´s INTERNAS DE LOS GENERADORES

Para realizar el análisis de las oscilaciones es necesario encontrar las fem´s internas

Para ello usaremos la siguiente ecuación:

E ∠δ ´ = V +( Q x´d
V
+j)P x ´d
V

δ 0=δ ´ +α

Generador 1

(
E1 ∠δ ´ 1 = 1.030 ∠0 ° +
0,712∗0,067
1.030 ∠ 0 °
+j )
3,500∗0,067
1.030 ∠ 0°
=1.1001 ∠11,9436 °

δ 1=δ ´ 1+ α 1=11,9436 ° +8,88 °=20,8236°

E1 ∠ δ 1 =1.1001∠ 20,8236 °

Generador 2

(
E2 ∠δ ´ 2 = 1,020 ∠ 0 °+
0,298∗0,10
1,020 ∠0 °
+j )
1,850∗0,10
1,020 ∠ 0 °
=1,0648 ∠9,8075 °

δ 2=δ ´ 2+ α 2 =9,8075° + 6,38° =16,1875°

E2 ∠δ 2 =1,0648∠ 16,1875 °

Para la barra infinita al ser la referencia se tiene:

E3 ∠δ 3=1,0000 ∠0,00°

12
ECUACIONES DE OSCILACIÓN DURANTE LA FALLA
Para hallar la potencia inyectada en el nodo i, suministrada por el generador i usamos la
ecuación:
3
Pei =E2i Gii + ∑ Ei E j Y ij cos ( θij −δ i +δ j )
j=1
j ≠i

GENERADOR 1
Y 12 =0,0271+ j0,1796=0,1816 ∠81,4193 °

Y 13 =−0,0199+ j 1,0664=1,0666 ∠ 91,0691 °

Y 23=−0,0708+ j 5,2710=5,2715 ∠ 90,7695 °

Pe 1=E21 G11+ E1 E2 Y 12 cos ( θ 12−δ1 +δ 2 ) + E1 E3 Y 13 cos ( θ 13−δ 1 + δ 3 )

2
Pe 1=1.1001 ∗0,2882+1.1001∗1,0648∗0,1816 cos ( 81,4193° −δ 1 +16,1875 ° ) +1.1001∗1∗1,0666 cos ( 91,0691

Pe 1=0,3488+1,3849 cos ( δ 1−92,0710 ° )

Pe 1=0,3488+1,3849 sen ( δ 1−2,0710° )

Ecuación de oscilación:

2 H i d 2 δi
=Pmi−P ei
ω s dt 2

Ecuación de oscilación de G1:

2
d δ1 ωs
2
= ( P −P e1 )
2 H 1 m1
dt

d2δ 1
dt 2
=
2∗π∗60
2∗11,2
[ 3,5−0,3488−1,3849 sen ( δ 1−2,0710° ) ]
rad elec
s2 [ ]
[ ]
2
d δ1 rad elec
2
=53,0346−23,3078 sen ( δ 1−2,0710 ° ) 2
dt s

GENERADOR 2

Pe 2=E22 G 22+ E 2 E1 Y 21 cos ( θ 21−δ 2 + δ 1 ) + E2 E3 Y 23 cos ( θ 23−δ 2+ δ 3 )

2
Pe 2=1,0648 ∗0,1387 +1,0648∗1.1001∗0,1816 cos ( 81,4193 °−δ 2 +20,8236 ° ) +1,0648∗1∗5,2715 cos ( 90,7695

Pe 2=0,1573+5,8217 cos ( δ 2−91,1859° )

13
 Pe 2=0,1573+5,8217 sen ( δ 2−1,1859° )

Ecuación de oscilación de G2:

2
d δ2 ωs
2
= ( P −P e2 )
2 H 2 m2
dt

[ ]
2
d δ2 2∗π∗60 rad elec
dt
2
=
2∗8
[ 1,85−0,1573−5,8217 sen ( δ 2 −1,1859° ) ]
s
2

[ ]
2
d δ2 rad elec
2
=39,8833−137,1706 sen ( δ 2−1,1859 ° ) 2
dt s

ECUACIONES DE OSCILACIÓN DESPUÉS DE LA FALLA

Para hallar la potencia inyectada en el nodo i, suministrada por el generador i es:

3
Pei =E Gii + ∑ Ei E j Y ij cos ( θij −δ i +δ j )
2
i
j=1
j ≠i

GENERADOR 1

Y 12 =0,0694+ j 0,6443=0,6480 ∠83,8522 °

Y 13=−0,2333+ j 3,7988=3,8060 ∠93,5144 °

Y 23=−0,0936+ j 5,5396=5,5404 ∠ 90,9680°

Pe 1=E21 G 11+ E1 E2 Y 12 cos ( θ 12−δ1 +δ 2 ) + E1 E3 Y 13 cos ( θ 13−δ 1 + δ 3 )

2
Pe 1=1.1001 ∗0,7500+1.1001∗1,0648∗0,6480 cos ( 83,8522°−δ 1 +16,1875 ) +1.1001∗1∗3,8060 cos ( 93,5144−

Pe 1=0,9077+ 4,9419 cos ( δ 1−94,5146 ° )

Pe 1=0,9077+ 4,9419 sen ( δ1−4,5146 ° )

Ecuación de oscilación:
2
2 H i d δi
=Pmi −P ei
ω s dt 2

Ecuación de oscilación de G1:

14
 d2δ 1 ωs
2
= ( P −P e1 )
2 H1 m1
dt

[ ]
2
d δ1 2∗π∗60 rad elec
dt
2
=
2∗11,2
[ 3,5−0,9077−4,9419 sen ( δ 1−4,5146° ) ] s 2

[ ]
2
d δ1 rad elec
2
=43,6283−83,1720 sen ( δ 1−4,5146 ° ) 2
dt s

GENERADOR 2

Pe 2=E22 G22+ E 2 E1 Y 21 cos ( θ 21−δ 2 + δ 1 ) + E2 E3 Y 23 cos ( θ 23−δ 2+ δ 3 )

2
Pe 2=1,0648 ∗0,1426 +1,0648∗1.1001∗0,6480 cos ( 83,8522 °−δ 2+ 20,8236° ) + 1,0648∗1∗5,5404 cos ( 90,9680

Pe 2=0,1617+6,6393 cos ( δ 2−92,5205° )

Pe 2=0,1617+6,6393 sen ( δ 2−2,5205 ° )

Ecuación de oscilación de G2

Generador 2
2
d δ2 ωs
2
= ( P −P e2 )
2 H 2 m2
dt

[ ]
2
d δ2 2∗π∗60 rad elec
dt
2
=
2∗8 [ 1,85−0,1617−6,6393 sen ( δ 2−2,5205 ° ) ]
s
2

d2δ 2
dt 2
=39,7796−156,4348 sen ( δ 2−2,5205 ° )
[ rad elec
s2 ]

4.1 MÉTODO PASO POR PASO

Ecuación Después de la falla Ecuación Durante de la falla

GENERADOR 1

E1 ∠δ 1 =1.1001∠ 20,8236 °

H=11,2 [ s ]

15
Pm=3,5

Pe 1=0,3488+1,3849 sen ( δ 1−2,0710° ) Durante la falla

Pe 1=0,9077+ 4,9419 sen ( δ1−4,5146 ° ) Después de la falla

t(seg) Pe Pa(n-1) KPa(n-1) ∆δn δn(°elect)

0 20,8236
0+ 0,7940 2,7060 20,8236
0 1,3530 3,2617 3,2617 20,8236
0,05 0,8679 2,6321 6,3452 9,6069 24,0853
0,10 1,0749 2,4251 5,8462 15,4531 33,6922
0,15 1,3629 2,1371 5,1520 20,6051 49,1452
0,20 5,3951 -1,8951 -4,5686 16,0365 69,7503
0,25 5,7924 -2,2924 -5,5263 10,5102 85,7868
0,30 5,8472 -2,3472 -5,6585 4,8517 96,2970
0,35 5,8165 -2,3165 -5,5844 -0,7327 101,1488
0,40 5,8234 -2,3234 -5,6011 -6,3338 100,4161
0,45 5,8495 -2,3495 -5,6639 -11,9976 94,0823
0,50 5,7338 -2,2338 -5,3850 -17,3826 82,0846
0,55 5,1956 -1,6956 -4,0875 -21,4701 64,7020
0,60 3,9988 -0,4988 -1,2023 -22,6725 43,2319
0,65 2,2736 1,2264 2,9565 -19,7160 20,5594
0,70 0,5913 2,9087 7,0121 -12,7038 0,8435
0,75 -0,4855 3,9855 9,6080 -3,0959 -11,8604
0,80 -0,7396 4,2396 10,2204 7,1245 -14,9563
0,85 -0,1490 3,6490 8,7966 15,9212 -7,8317
0,90 1,2158 2,2842 5,5065 21,4276 8,0895
0,95 2,9964 0,5036 1,2140 22,6416 29,5171
1,00 4,5596 -1,0596 -2,5545 20,0871 52,1587

16
 Oscilación Punto por punto
120.0000

100.0000

80.0000

60.0000

40.0000

20.0000

0.0000
0 1 2

-20.0000

-40.0000

El tiempo crítico de desconexión estará en el intervalo de:

0,10<t cr <0,15 [ seg ]

El ángulo crítico de desconexión estará en el intervalo de:

49,1452< δ cr <69,7503 [ ° ]

En estas condiciones el ángulo delta, no es inestable, pero tampoco es estable si no que el ángulo
oscila, y si tomamos un intervalo mayor el ángulo delta aumenta exponencialmente, provocando
inestabilidad.

17
GENERADOR 2

E2 ∠δ 2 =1,0648∠ 16,1875 °

H=8 [ s ]

Pm=1,85

Pe 2=0,1573+5,8217 sen ( δ 2−1,1859° ) Durante la falla

Pe 2=0,1617+6,6393 sen ( δ 2−2,5205 ° ) Después de la falla

t(seg) Pe Pa(n-1) KPa(n-1) ∆δn δn(°elect)

0 16,1875
0+ 1,6642 0,1858 16,1875
0 0,0929 0,3135 0,3135 16,1875
0,05 1,7657 0,0843 0,2845 0,5980 16,5010
0,10 1,8329 0,0171 0,0579 0,6559 17,0990
0,15 1,9063 -0,0563 -0,1900 0,4659 17,7549
0,20 1,9583 -0,1083 -0,3656 0,1003 18,2208
0,25 1,9695 -0,1195 -0,4033 -0,3031 18,3210
0,30 1,9357 -0,0857 -0,2892 -0,5923 18,0180
0,35 1,8695 -0,0195 -0,0657 -0,6580 17,4257
0,40 1,7957 0,0543 0,1834 -0,4746 16,7677
0,45 1,7423 0,1077 0,3634 -0,1112 16,2931
0,50 1,7298 0,1202 0,4057 0,2945 16,1819
0,55 1,7629 0,0871 0,2938 0,5884 16,4765
0,60 1,8290 0,0210 0,0708 0,6592 17,0648
0,65 1,9028 -0,0528 -0,1783 0,4808 17,7240
0,70 1,9565 -0,1065 -0,3596 0,1212 18,2048
0,75 1,9701 -0,1201 -0,4052 -0,2840 18,3260
0,80 1,9384 -0,0884 -0,2983 -0,5823 18,0420
0,85 1,8733 -0,0233 -0,0786 -0,6608 17,4598
0,90 1,7992 0,0508 0,1715 -0,4893 16,7989
0,95 1,7442 0,1058 0,3572 -0,1321 16,3096
1,00 1,7293 0,1207 0,4074 0,2752 16,1775

18
 Oscilación Punto por Punto
18.5000

18.0000

17.5000

17.0000

16.5000

16.0000

15.5000

15.0000
0 1 2

El tiempo crítico de desconexión estará en el intervalo de:

0<t cr <0 [ seg ]

El ángulo crítico de desconexión estará en el intervalo de:

16,1875<δ cr < 16,5010 [ ° ]

En estas condiciones el ángulo delta, no es inestable, pero tampoco es estable si no que el ángulo
oscila, y si tomamos un intervalo mayor el ángulo delta aumenta exponencialmente, provocando
inestabilidad.

19
 4.2. MÉTODO EULER MEJORADO

Ecuación Después de la falla Ecuación Durante de la falla

GENERADOR 1

E1 ∠δ 1 =1.1001∠ 20,8236 °

H =11,2 [ s ]

Pm=3,5

Pe 1=0,3488+1,3849 sen ( δ 1−2,0710° ) Durante la falla

Pe 1=0,9077+ 4,9419 sen ( δ1−4,5146 ° ) Después de la falla

Pendiente prin. de
Interval
paso Valores predecidos Pendiente final de paso Valores corregidos δn°
o
dδn/dt d(∆w)n/dt δn ∆wn dδn/dt d(∆ω)n/dt δn(rad) ∆wn(rad/s)
0,00 0 0 0 0 0 0 0,3634 0 20,8236
0,01 0,0000 45,5415 0,3634 0,4554 0,4554 45,5415 0,3657 0,4554 20,9541
0,02 0,4554 45,4913 0,3703 0,9103 0,9103 45,4913 0,3725 0,9103 21,3453
0,03 0,9103 45,3409 0,3816 1,3637 1,3637 45,3409 0,3839 1,3637 21,9968
0,04 1,3637 45,0912 0,3976 1,8146 1,8146 45,0912 0,3998 1,8146 22,9073
0,05 1,8146 44,7440 0,4180 2,2621 2,2621 44,7440 0,4202 2,2621 24,0752
0,06 2,2621 44,3017 0,4428 2,7051 2,7051 44,3017 0,4450 2,7051 25,4982
0,07 2,7051 43,7678 0,4721 3,1428 3,1428 43,7678 0,4743 3,1428 27,1735
0,08 3,1428 43,1465 0,5057 3,5742 3,5742 43,1465 0,5079 3,5742 29,0978
0,09 3,5742 42,4433 0,5436 3,9987 3,9987 42,4433 0,5457 3,9987 31,2673
0,10 3,9987 41,6649 0,5857 4,4153 4,4153 41,6649 0,5878 4,4153 33,6777
0,11 4,4153 40,8193 0,6319 4,8235 4,8235 40,8193 0,6340 4,8235 36,3245
0,12 4,8235 39,9157 0,6822 5,2227 5,2227 39,9157 0,6842 5,2227 39,2025
0,13 5,2227 38,9649 0,7364 5,6123 5,6123 38,9649 0,7384 5,6123 42,3065
0,14 5,6123 37,9793 0,7945 5,9921 5,9921 37,9793 0,7964 5,9921 45,6309
0,15 5,9921 36,9729 0,8563 6,3619 6,3619 36,9729 0,8582 6,3619 49,1701
0,16 6,3619 35,9609 0,9218 6,7215 6,7215 35,9609 0,9236 6,7215 52,9182
0,17 6,7215 34,9602 0,9908 7,0711 7,0711 34,9602 0,9926 7,0711 56,8694
0,18 7,0711 33,9891 1,0633 7,4110 7,4110 33,9891 1,0650 7,4110 61,0182
0,19 7,4110 33,0669 1,1391 7,7416 7,7416 33,0669 1,1407 7,7416 65,3591
0,20 7,7416 -29,0059 1,2181 7,4516 7,4516 -29,0059 1,2167 7,4516 69,7117
0,21 7,4516 -31,8716 1,2912 7,1328 7,1328 -31,8716 1,2896 7,1328 73,8898
0,22 7,1328 -34,2129 1,3609 6,7907 6,7907 -34,2129 1,3592 6,7907 77,8786
0,23 6,7907 -36,0623 1,4271 6,4301 6,4301 -36,0623 1,4253 6,4301 81,6661
0,24 6,4301 -37,4612 1,4896 6,0555 6,0555 -37,4612 1,4878 6,0555 85,2429
0,25 6,0555 -38,4571 1,5483 5,6709 5,6709 -38,4571 1,5464 5,6709 88,6023

20
0,26 5,6709 -39,1013 1,6031 5,2799 5,2799 -39,1013 1,6012 5,2799 91,7395
0,27 5,2799 -39,4461 1,6540 4,8854 4,8854 -39,4461 1,6520 4,8854 94,6516
0,28 4,8854 -39,5434 1,7008 4,4900 4,4900 -39,5434 1,6989 4,4900 97,3375
0,29 4,4900 -39,4427 1,7438 4,0956 4,0956 -39,4427 1,7418 4,0956 99,7971
0,30 4,0956 -39,1904 1,7827 3,7037 3,7037 -39,1904 1,7808 3,7037 102,0314
0,31 3,7037 -38,8289 1,8178 3,3154 3,3154 -38,8289 1,8159 3,3154 104,0422
0,32 3,3154 -38,3964 1,8490 2,9314 2,9314 -38,3964 1,8471 2,9314 105,8318
0,33 2,9314 -37,9265 1,8764 2,5522 2,5522 -37,9265 1,8745 2,5522 107,4027
0,34 2,5522 -37,4484 1,9001 2,1777 2,1777 -37,4484 1,8982 2,1777 108,7577
0,35 2,1777 -36,9870 1,9200 1,8078 1,8078 -36,9870 1,9181 1,8078 109,8995
0,36 1,8078 -36,5633 1,9362 1,4422 1,4422 -36,5633 1,9344 1,4422 110,8305
0,37 1,4422 -36,1941 1,9488 1,0802 1,0802 -36,1941 1,9470 1,0802 111,5531
0,38 1,0802 -35,8931 1,9578 0,7213 0,7213 -35,8931 1,9560 0,7213 112,0692
0,39 0,7213 -35,6703 1,9632 0,3646 0,3646 -35,6703 1,9614 0,3646 112,3803
0,40 0,3646 -35,5330 1,9651 0,0093 0,0093 -35,5330 1,9633 0,0093 112,4874
0,41 0,0093 -35,4851 1,9634 -0,3456 -0,3456 -35,4851 1,9616 -0,3456 112,3911
0,42 -0,3456 -35,5282 1,9581 -0,7009 -0,7009 -35,5282 1,9564 -0,7009 112,0913
0,43 -0,7009 -35,6607 1,9494 -1,0575 -1,0575 -35,6607 1,9476 -1,0575 111,5875
0,44 -1,0575 -35,8784 1,9370 -1,4163 -1,4163 -35,8784 1,9352 -1,4163 110,8789
0,45 -1,4163 -36,1744 1,9210 -1,7780 -1,7780 -36,1744 1,9192 -1,7780 109,9638
0,46 -1,7780 -36,5384 1,9014 -2,1434 -2,1434 -36,5384 1,8996 -2,1434 108,8404
0,47 -2,1434 -36,9574 1,8782 -2,5130 -2,5130 -36,9574 1,8763 -2,5130 107,5064
0,48 -2,5130 -37,4147 1,8512 -2,8871 -2,8871 -37,4147 1,8493 -2,8871 105,9594
0,49 -2,8871 -37,8899 1,8205 -3,2660 -3,2660 -37,8899 1,8186 -3,2660 104,1967
0,50 -3,2660 -38,3590 1,7859 -3,6496 -3,6496 -38,3590 1,7840 -3,6496 102,2155

21
 140 Oscilación Euler

120

100

80

60

40

20

0
0.00 0.50 1.00 1.50 2.00 2.50

-20

-40

El tiempo crítico de desconexión está en el intervalo de:

0,18<t cr <0,19 [ seg ]

El ángulo crítico de desconexión está en el intervalo de:

65,3591<δ cr < 69,7117 [ ° ]

En estas condiciones el ángulo delta, no es inestable, pero tampoco es estable si no que el ángulo
oscila, y si tomamos un intervalo mayor el ángulo delta aumenta exponencialmente, provocando
inestabilidad.

GENERADOR 2

E2 ∠δ 2 =1,0648∠ 16,1875 °

H =8 [ s ]

22
 Pm=1,85

Pe 2=0,1573+5,8217 sen ( δ 2−1,1859° ) Durante la falla

Pe 2=0,1617+6,6393 sen ( δ 2−2,5205 ° ) Después de la falla

Pendiente prin. de Pendiente final de

Interval paso Valores predecidos paso Valores corregidos
δn°
o δn(rad
dδn/dt d(∆w)n/dt δn ∆wn dδn/dt d(∆ω)n/dt ∆wn(rad/s)
)
0,00 0 0 0 0 0 0 0,2825 0 16,1875
0,01 0,0000 2,8174 0,2825 0,0282 0,0282 2,8174 0,2827 0,0282 16,1956
0,02 0,0282 2,7960 0,2829 0,0561 0,0561 2,7960 0,2831 0,0561 16,2197
0,03 0,0561 2,7320 0,2836 0,0835 0,0835 2,7320 0,2838 0,0835 16,2597
0,04 0,0835 2,6259 0,2846 0,1097 0,1097 2,6259 0,2848 0,1097 16,3151
0,05 0,1097 2,4791 0,2858 0,1345 0,1345 2,4791 0,2860 0,1345 16,3850
0,06 0,1345 2,2937 0,2873 0,1574 0,1574 2,2937 0,2874 0,1574 16,4687
0,07 0,1574 2,0720 0,2890 0,1782 0,1782 2,0720 0,2891 0,1782 16,5648
0,08 0,1782 1,8173 0,2909 0,1963 0,1963 1,8173 0,2910 0,1963 16,6721
0,09 0,1963 1,5332 0,2929 0,2117 0,2117 1,5332 0,2930 0,2117 16,7890
0,10 0,2117 1,2238 0,2951 0,2239 0,2239 1,2238 0,2952 0,2239 16,9137
0,11 0,2239 0,8937 0,2974 0,2328 0,2328 0,8937 0,2975 0,2328 17,0446
0,12 0,2328 0,5478 0,2998 0,2383 0,2383 0,5478 0,2998 0,2383 17,1796
0,13 0,2383 0,1912 0,3022 0,2402 0,2402 0,1912 0,3022 0,2402 17,3167
0,14 0,2402 -0,1709 0,3046 0,2385 0,2385 -0,1709 0,3046 0,2385 17,4538
0,15 0,2385 -0,5328 0,3070 0,2332 0,2332 -0,5328 0,3070 0,2332 17,5890
0,16 0,2332 -0,8892 0,3093 0,2243 0,2243 -0,8892 0,3093 0,2243 17,7200
0,17 0,2243 -1,2346 0,3115 0,2120 0,2120 -1,2346 0,3115 0,2120 17,8450
0,18 0,2120 -1,5638 0,3136 0,1963 0,1963 -1,5638 0,3135 0,1963 17,9620
0,19 0,1963 -1,8717 0,3155 0,1776 0,1776 -1,8717 0,3154 0,1776 18,0691
0,20 0,1776 -2,1536 0,3171 0,1561 0,1561 -2,1536 0,3170 0,1561 18,1647
0,21 0,1561 -2,4049 0,3186 0,1320 0,1320 -2,4049 0,3185 0,1320 18,2472
0,22 0,1320 -2,6219 0,3198 0,1058 0,1058 -2,6219 0,3197 0,1058 18,3153
0,23 0,1058 -2,8009 0,3207 0,0778 0,0778 -2,8009 0,3206 0,0778 18,3679
0,24 0,0778 -2,9390 0,3214 0,0484 0,0484 -2,9390 0,3212 0,0484 18,4041
0,25 0,0484 -3,0340 0,3217 0,0181 0,0181 -3,0340 0,3215 0,0181 18,4231
0,26 0,0181 -3,0840 0,3217 -0,0128 -0,0128 -3,0840 0,3216 -0,0128 18,4246
0,27 -0,0128 -3,0880 0,3214 -0,0437 -0,0437 -3,0880 0,3213 -0,0437 18,4085
0,28 -0,0437 -3,0455 0,3209 -0,0741 -0,0741 -3,0455 0,3207 -0,0741 18,3747
0,29 -0,0741 -2,9569 0,3200 -0,1037 -0,1037 -2,9569 0,3198 -0,1037 18,3238
0,30 -0,1037 -2,8231 0,3188 -0,1319 -0,1319 -2,8231 0,3186 -0,1319 18,2563
0,31 -0,1319 -2,6457 0,3173 -0,1584 -0,1584 -2,6457 0,3172 -0,1584 18,1731
0,32 -0,1584 -2,4271 0,3156 -0,1826 -0,1826 -2,4271 0,3155 -0,1826 18,0754
0,33 -0,1826 -2,1702 0,3136 -0,2043 -0,2043 -2,1702 0,3135 -0,2043 17,9646
0,34 -0,2043 -1,8786 0,3115 -0,2231 -0,2231 -1,8786 0,3114 -0,2231 17,8421

23
0,35 -0,2231 -1,5562 0,3092 -0,2387 -0,2387 -1,5562 0,3091 -0,2387 17,7098
0,36 -0,2387 -1,2077 0,3067 -0,2508 -0,2508 -1,2077 0,3066 -0,2508 17,5696
0,37 -0,2508 -0,8381 0,3041 -0,2592 -0,2592 -0,8381 0,3041 -0,2592 17,4235
0,38 -0,2592 -0,4528 0,3015 -0,2637 -0,2637 -0,4528 0,3015 -0,2637 17,2737
0,39 -0,2637 -0,0575 0,2988 -0,2643 -0,2643 -0,0575 0,2988 -0,2643 17,1225
0,40 -0,2643 0,3420 0,2962 -0,2608 -0,2608 0,3420 0,2962 -0,2608 16,9720
0,41 -0,2608 0,7396 0,2936 -0,2534 -0,2534 0,7396 0,2936 -0,2534 16,8247
0,42 -0,2534 1,1292 0,2911 -0,2421 -0,2421 1,1292 0,2912 -0,2421 16,6827
0,43 -0,2421 1,5050 0,2887 -0,2271 -0,2271 1,5050 0,2888 -0,2271 16,5483
0,44 -0,2271 1,8609 0,2866 -0,2085 -0,2085 1,8609 0,2866 -0,2085 16,4235
0,45 -0,2085 2,1916 0,2846 -0,1866 -0,1866 2,1916 0,2847 -0,1866 16,3104
0,46 -0,1866 2,4916 0,2828 -0,1617 -0,1617 2,4916 0,2829 -0,1617 16,2106
0,47 -0,1617 2,7562 0,2813 -0,1341 -0,1341 2,7562 0,2814 -0,1341 16,1259
0,48 -0,1341 2,9810 0,2801 -0,1043 -0,1043 2,9810 0,2803 -0,1043 16,0576
0,49 -0,1043 3,1622 0,2792 -0,0727 -0,0727 3,1622 0,2794 -0,0727 16,0069
0,50 -0,0727 3,2968 0,2786 -0,0397 -0,0397 3,2968 0,2788 -0,0397 15,9747

Oscilación Euler
25

20

15

10

0
0.00 0.50 1.00 1.50 2.00 2.50

24
El tiempo crítico de desconexión estará en el intervalo de:

0<t cr <0 [ seg ]

El ángulo crítico de desconexión estará en el intervalo de:

16,1875<δ cr < 16,1956 [ ° ]

En estas condiciones el ángulo delta, no es inestable, pero tampoco es estable si no que el ángulo
oscila, y si tomamos un intervalo mayor el ángulo delta aumenta exponencialmente, provocando
inestabilidad.

5. COMPARACIÓN DE LOS MÉTODOS

Se analizara los resultados comparando los 2 generadores:

 Potencia eléctrica antes y después de la falla

Generador 1

Pe 1=0,3488+1,3849 sen ( δ 1−2,0710° ) Durante la falla

Pe 1=0,9077+ 4,9419 sen ( δ1−4,5146 ° ) Después de la falla

Generados 2

Pe 2=0,1573+5,8217 sen ( δ 2−1,1859° ) Durante la falla

Pe 2=0,1617+6,6393 sen ( δ 2−2,5205 ° ) Después de la falla

- Se puede ver que la potencia eléctrica de G1 Después de la Falla es mucho mayor que
Durante de la Falla, esto es correcto ya que la línea 3-4 está conectada directamente a G1, lo
que repercute y perturba directamente a G1

- En cambio la potencia eléctrica de G2 Después de la Falla es levemente mayor a la

potencia eléctrica Durante de la Falla, esto se podría interpretar que la línea 3-4, donde
ocurre la falla no está conectada directamente a G2, este es afectado levemente.

 Ángulo y tiempo crítico de desconexión

Generador G1

Método Punto por Punto

0,10<t cr <0,15 [ seg ]

49,1452< δ cr <69,7503 [ ° ]

25
En estas condiciones el ángulo delta, no es inestable, pero tampoco es estable si no que el ángulo
oscila.

Para un tiempo de desconexión 0,15<t cr < 0,20 el ángulo δ se incrementaría exponencialmente lo

que genera inestabilidad.

Método Euler Mejorado

0,18<t cr <0,19 [ seg ]

65,3591<δ cr < 69,7117 [ ° ]

En estas condiciones el ángulo delta, no es inestable, pero tampoco es estable si no que el ángulo
oscila.

Para un tiempo de desconexión 0,19<t cr < 0,2 0 el ángulo δ se incrementaría exponencialmente lo

que genera inestabilidad.

Generador G2

Método Punto por Punto

0<t cr <0 [ seg ]

16,1875<δ cr < 16,5010 [ ° ]

En estas condiciones el ángulo delta, no es inestable, pero tampoco es estable si no que el ángulo
oscila.

Método Euler Mejorado

0<t cr <0 [ seg ]

16,1875<δ cr < 16,1956 [ ° ]

En estas condiciones el ángulo delta, no es inestable, pero tampoco es estable si no que el ángulo
oscila.

- El Generador G1, por el Método Punto por Punto su intervalos es de 0,10<t cr <0,15 [ seg ] y
por el Método Euler Mejorado es 0,18<t cr <0,19 [ seg ], los intervalos de los dos métodos no
son los mismos aunque cabe decir que si se aproximan entre ello, esto de puede deber a que
los métodos que usamos son métodos aproximados y no exactos por lo que tienen errores, en
cuanto al ángulo critico estos si se encuentran en el mismo intervalo dando así una región de
donde podría estar el ángulo de desconexión.

- El Generador G2, según podemos ver su ángulo crítico es 0 y su δ cr es prácticamente δ 0 ,

esto podemos interpretar de la siguiente manera, como la falla ocurre en la línea 3-4, y esta línea está

26
alejada de G2, la perturbación en G2 es muy leve, haciendo que esta falla no provoque inestabilidad
en el generador G2.

27