Universidad Técnica de Cotopaxi

• Asignatura: Sistemas Eléctricos de Potencia 1

• Tema: Análisis de Flujos de Potencia

 SISTEMAS ELÉCTRICOS DE POTENCIA I

Análisis nodal. Matriz de admitancias de nudos

Dada la red pasiva, determinar la matriz de admitancias cuando a=1

2
 SISTEMAS ELÉCTRICOS DE POTENCIA I

Análisis nodal. Matriz de admitancias de nudos

•Opción A: paso a paso

1 2 3

 1 1 y pl1− 3 1 1 
 + + − − 
 z sl1− 2 z sl1−3 2 z sl1− 2 z sl1− 3 
 
 1 1 1 1 
YB = − + − 
z sl1− 2 z sl1− 2 z cc z cc
 
 
 1 1 1 y pl1− 3 1 
− − + + + jB 
 z sl1− 3 z cc z sl1− 3 2 z cc 

3
 SISTEMAS ELÉCTRICOS DE POTENCIA I

Análisis nodal. Matriz de admitancias de nudos

•Opción B: Adición matriz de admitancia de elementos

Línea 1-2: línea corta Línea 1-3: línea media Trafo 3-2
1 2 1 3 3 2
 1 1   1 y pl1− 3   1 1 
−  + −
1
  2 −
 z sl1− 2  a z cc 
=  z sl1− 2 =  z sl1− 3 2 z sl1− 3  Y t 3− 2 =  a z cc 
1 
Y l1− 2 Y l1−3
− 1  1 1 y pl1− 3  − 1 1 
 z sl1− 2 z sl1− 2   − +   a z cc z cc 
 z sl1− 3 z sl1− 3 2 
a =1

1 2 3

 1 1 y pl1− 3 1 1 
 + + − − 
 z sl1− 2 z sl1−3 2 z sl1− 2 z sl1− 3 
 
 1 1 1 1 
YB = − + − 
z sl1− 2 z sl1− 2 z cc z cc
 
 
 1 1 1 y pl1− 3 1 
 − − + + + jB  Batería de
z sl1− 3 z cc z sl1− 3 2 z cc  condensadores

4
 SISTEMAS ELÉCTRICOS DE POTENCIA I

Análisis nodal. Matriz de admitancias de nudos

1 2 3
 1 1 yj0,04
pl1− 3 1 1 
 + + − − 
z sl1− 2 0,04+j0,2
0,05+j0,2 z sl1−3 2 z sl1− 2
0,05+j0,2 z sl1− 3
0,04+j0,2 
 
 1 1 1 1 
YB = − + − 
z sl1− 2 z sl1− 2 z cc ....z cc
 0,05+j0,2 0,05+j0,2 0,03+j0,1 0,03+j0,1

 
 1 1 1 yj0,04
pl1− 3 1 
 − − + + + jB
0,5

z sl1− 3
0,04+j0,2 ....z cc
0,03+j0,1 z sl1− 3
0,04+j0,2 2 .... z cc
0,03+j0,1

 2,1380 − j 9,4936 − 1,1765 + j 4,7059 − 0,9615 + j 4,8077

YB =  − 1,1765 + j 4,7059 3,9288 − j13,8802 − 2,7523 + j9,1743
− 0,9615 + j 4,8077 − 2,7523 + j9,1743 3,7138 − j13,9620 

9,7313 − 77,30º 4,8507104,04º 4,9029101,31º 

YB =  4,8507104,04º 14,4255 − 74,19º 9,5783106,70º 

 4,9029101,31º 9,5783106,70º 14,4475 − 75,10º 

5
 SISTEMAS ELÉCTRICOS DE POTENCIA I

Análisis nodal. Matriz de admitancias de nudos

Realizar los mismos cálculos si a = 1,2

1 2 3
 1 1 y pl1− 3 1 1 
 + +
j0,04
− − 
 z0,05+j0,2
sl1 − 2 z sl1 − 3
0,04+j0,2 2 z sl1− 2
0,05+j0,2 z sl1− 3
0,04+j0,2 
 
 1 1 1 
 − + −
1 
YB = z sl1− 2 z sl1− 2 z cc 
0,05+j0,2 0,05+j0,2 0,03+j0,1 .... zcc
1,2 (0,03+j0,1)
 
 1 1 yj0,04 
 1 pl1− 3 1
− − + + 0,5 
+ jB
 z sl1− 3
0,04+j0,2 ....(0,03+j0,1)
1,2 zcc z sl1− 3
0,04+j0,2 2 .... 2 z
1,2 (0,03+j0,1)
cc 
 

 2,1380 − j9,4936 − 1,1765 + j 4,7059 − 0,9615 + j 4,8077

Y B =  − 1,1765 + j 4,7059 3,9288 − j13,8802 − 2,2936 + j 7,6453
− 0,9615 + j 4,8077 − 2,2936 + j 7,6453 2,8728 − j10,6587 

9,7313 − 77,30º 4,8507104,04º 4,9029101,31º 

YB =  4,8507104,04º 14,4255 − 74,19º 7,9819106,70º 

 4,9029101,31º 7,9819106,7º 11,0390 − 74,92º 

6
 SISTEMAS ELÉCTRICOS DE POTENCIA I

Flujos de potencia

Realizar el flujo de cargas en el sistema de la figura para determinar

tensiones en todos los nudos y determinar flujos de potencia por líneas
y transformadores
2 ?

3 incógnitas
2 U3 ? 3 ecuaciones
3 ?
3 P2
U3 P3
Q3

Nudo Tipo Tensión Pg Qg Pd Qd Pe = Pg-Pd Qe= Qg-Qd

1 Balance 10º ? ? 0,5 0,25 ? ?

2 PU 12 ? 0,75 ? 0 0 0,75 ?

3 PQ U3 ?3 ? 0 0 0,5 0,5 -0,5 -0,5

7
 SISTEMAS ELÉCTRICOS DE POTENCIA I

Flujos de potencia

Planteamiento ecuaciones flujo de cargas

P2 = 0,75 = U2U1Y21 Cos(2 - 1 - 21 ) + U2U2Y22 Cos( - 22 ) +U2U3Y23 Cos(2 - 3 - 23 )
P3 = -0,5 = U3U1Y31 Cos(3 - 1 - 31 ) + U3U2Y32 Cos(3 - 2 - 32 ) + U3U3Y33 Cos( - 33 )
Q3 = -0,5 = U3U1Y31 Sen(3 - 1 - 31 ) + U3U2Y32 Sen(3 - 2 - 32 ) + U3U3Y33 Sen( - 33 )
Resolución de ecuaciones mediante Newton-Raphson
θ2 θ3 U3

P2  P2 P2
 e c 
P2   Δθ 2   P2 − P2 = ΔP2 
Ángulos en
 U3    radianes
Matriz  θ 2 θ 3 U3     e 
P3  P3 P3 P  Δθ   P − P3 = ΔP3 
c
jacobiana J
 θ U3 3    3  =  3 
θ 3 U3  
 2   
Q 3  Q 3 Q 3 Q 3   ΔU3  Q e − Q c = ΔQ 
U3  3 3 3
 θ 2 θ 3 U3   U3   

P2 P2
= - U2U1Y21 Sen(2 - 1 - 21 ) - U2U3Y23 Sen(2 - 3 - 23 ) = U2U3Y23 Sen(2 - 3 - 23 )
 2 3
P3 P3
= U3U2Y32 Sen(3 - 2 - 32 ) = -U3U1Y31 Sen(3 - 1 - 31 ) - U3U2Y32 Sen(3 - 2 - 32 )
 2 3
Q 3 Q 3
= -U3U2Y32 Cos(3 - 2 - 32 ) = U3U1Y31 Cos(3 - 1 - 31 ) + U3U2Y32 Cos(3 - 2 - 32 )
 2 3
P2
U3 = U2U3Y23 Cos(2 - 3 - 23 )
U3
P3
U3 = U3U1Y31 Cos(3 - 1 - 31 ) + U3U2Y32 Cos(3 - 2 - 32 ) + 2  U3U3Y33 Cos( - 33 )
U 3
Q 3
U3 = U3U1Y31 Sen(3 - 1 - 31 ) + U3U2Y32 Sen(3 - 2 - 32 ) + 2  U3U3Y33 Sen( - 33 )
U 3
8
 Datos
1 = 0 U1 = 1 U2 = 1
9,7313 − 77,30º 4,8507104,04º 4,9029101,31º 
YB =  4,8507104,04º 14,4255 − 74,19º 7,9819106,70º 
 4,9029101,31º 7,9819106,7º 11,0390 − 74,92º 

Iteración 1
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0 0 1 0,75-0,4587= -0,5+0,3823= -0,5+1,7942= 0,0417 -0,0101 0,1315
0,2913 -0,1177 1,2942

P2c = U2U1Y21 Cos(2 - 1 - 21 ) + U2U2Y22 >

Cos( - 22 ) +U2U
Tolerancia 3Y23 Cos(2 - 3 - 23 ) =
(1e-4)
= 114,8507 Cos(-104,04º ) + 1114,4255 Cos( 74,19º) +117,9819 Cos(-106,7º) = 0,4587
P3c = -0,3823
Sigue proceso iterativo
Q3c = -1,7942  
 12,35114 − 7,64526 − 2,29358   ΔP2   Δθ2 
     
J = − 7,64526 12,45295 2,49058  J−1   ΔP3  =  Δθ3 
 2,29358 − 3,25512 8,86453  ΔQ3   ΔU3 
 U 
 3 

9
 Datos
1 = 0 U1 = 1 U2 = 1
9,7313 − 77,30º 4,8507104,04º 4,9029101,31º 
YB =  4,8507104,04º 14,4255 − 74,19º 7,9819106,70º 
 4,9029101,31º 7,9819106,7º 11,0390 − 74,92º 

Iteración 1
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0 0 1 0,75-0,4587= -0,5+0,3823= -0,5+1,7942= 0,0417 -0,0101 0,1315
0,2913 -0,1177 1,2942

Iteración 2
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0,0417 -0.0101 1.1315 0,75-0,8064= -0,5+0,5051= -0,5+0,2870= -6,7e-3 -5e-4 -0,0165
-0,0564 0,0051 -0,2130

P2c = U2U1Y21 Cos(2 - 1 - 21 ) + U2U2Y22 Cos( - 22 ) +U2U3Y23 Cos(2 - 3 - 23 ) =
= 114,8507 Cos(2,39º-104,04º ) + 1114,4255
> ToleranciaCos( 74,19º) +11,13157,9819 Cos(2,39º+0,58º-106,7º) = 0,8064
(1e-4)
P3c = -0,5051 Sigue proceso iterativo
Q3c = -0,2870  
 13,5242 − 8,7734 − 2,1432  ΔP2   Δθ2 
     
J = − 8,5043 13,9328 − 3,1728 J−1   ΔP3  =  Δθ3 
 3,0401 − 4,1831 13,3589  ΔQ3   ΔU3 
 U 
 3 

10
 Datos
1 = 0 U1 = 1 U2 = 1
9,7313 − 77,30º 4,8507104,04º 4,9029101,31º 
YB =  4,8507104,04º 14,4255 − 74,19º 7,9819106,70º 
 4,9029101,31º 7,9819106,7º 11,0390 − 74,92º 

Iteración 1
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0 0 1 0,75-0,4587= -0,5+0,3823= -0,5+1,7942= 0,0417 -0,0101 0,1315
0,2913 -0,1177 1,2942

Iteración 2
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0,0417 -0.0101 1.1315 0,75-0,8064= -0,5+0,5051= -0,5+0,2870= -6,7e-3 -5e-4 -0,0165
-0,0564 0,0051 -0,2130

Iteración 3
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0,035 -0.0106 1.1150 0,75-0,7509= -0,5+0,4999= -0,5+0,4967= -2e-4 -2e-5 -2e-4
-0,0009 -0,0001 -0,00330

P2c = U2U1Y21 Cos(2 - 1 - 21 ) + U2U2Y22 Cos( - 22 ) +U2U3Y23 Cos(2 - 3 - 23 ) =
= 114,8507 Cos(2º-104,04º ) + 1114,4255 Cos( 74,19º)(1e-4)
> Tolerancia +11,11507,9819 Cos(2º+0,61º-106,7º) = 0,7509
P3c = -0,4999 Sigue proceso iterativo
Q3c = -0,4967  
 13,3763 − 8,6322 − 2,1666   ΔP2   Δθ2 
     
J = − 8,3994 13,7484 − 3,0718 J−1   ΔP3  =  Δθ3 
 2,9428 − 4,0716 12,7550  ΔQ3   ΔU3 
 U 
 3 
11
 Datos
1 = 0 U1 = 1 U2 = 1
9,7313 − 77,30º 4,8507104,04º 4,9029101,31º 
YB =  4,8507104,04º 14,4255 − 74,19º 7,9819106,70º 
 4,9029101,31º 7,9819106,7º 11,0390 − 74,92º 

Iteración 1
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0 0 1 0,75-0,4587= -0,5+0,3823= -0,5+1,7942= 0,0417 -0,0101 0,1315
0,2913 -0,1177 1,2942

Iteración 2
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0,0417 -0.0101 1.1315 0,75-0,8064= -0,5+0,5051= -0,5+0,2870= -6,7e-3 -5e-4 -0,0165
-0,0564 0,0051 -0,2130

Iteración 3
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0,035 -0.0106 1.1150 0,75-0,7509= -0,5+0,4999= -0,5+0,4967= -2e-4 -2e-5 -2e-4
-0,0009 -0,0001 -0,00330

Iteración 4
2 rad 3 rad U3 P2 P3 Q3 2 3 U3
0,0348 -0.0106 1.1148 0,75-0,7500= -0,5+0,4999= -0,5+0,4999=
-2,1e-7 -6,5e-8 -8,7e-7
P2c = 0,7500
P3c = -0,4999 < Tolerancia (1e-4)
Q3c = -0,4999 FIN proceso iterativo 12
 SISTEMAS ELÉCTRICOS DE POTENCIA I

Flujos de potencia

SOLUCIÓN
U1 = 10º p.u. U2 = 11,9938º p.u. U3 = 1,1148-0,6073º p.u.

Reparto de cargas

Nudo 1 balance
P1= U1U1Y11 Cos( - 11 ) + U1U2Y12Cos(1- 2 - 12 ) +U1U3Y13 Cos(1- 3 - 13 ) = -0,2166 p.u.
Q1= U1U1Y11 Sen( - 11 ) + U1U2Y12Sen(1- 2 - 12 ) +U1U3Y13 Sen(1- 3 - 13 ) = -0,5389 p.u.

Pg1 = P1 + Pd1 = -0,2166 + 0,5 = 0,2834 p.u.

Qg1 = Q1 + Qd1 = -0,5389 + 0,25 = -0,2889 p.u.

Nudo 2 PU
Q2= U2U1Y21 Sen(2 - 1 - 21 ) + U2U2Y22 Sen( - 22 ) +U2U3Y23 Sen(2 - 3 - 23 ) = 0,5063 p.u.

Qg2 = Q2 + Qd2 = 0,5063+0 = 0,5063 p.u.

Pérdidas
Ppérdidas = P1+ P2 + P3 = -0,2166 + 0,75 – 0,5 = 0,0334 p.u.

Qpérdidas = Q1+ Q2 + Q3 + BU32= -0,5389 + 0,5063 – 0,5 + 0,6213 = 0,0887 p.u.

13
 SISTEMAS ELÉCTRICOS DE POTENCIA I

Flujos de potencia

SOLUCIÓN U1 = 10º U2 = 11,9938º U3 = 1,1148-0,6073º

Flujos por líneas y transformador

Línea 1-2
S12= U1I12* = -0,1632 + j0,0438 p.u.
S21= U2I21* = 0,1646 - j0,0381 p.u.
 1 1 
−
I 12   z sl1−2 z sl1−2   U 1  =  1,1765 − j4,7059 − 1,1765 + j4,7059   1 
 =    
1  U 2  − 1,1765 + j4,7059 1,1765 − j4,7059  11,9938º
 
 I 21  − 1   
 z sl1−2 z sl1−2  Y l1−2
 
Y l1−2
Línea 1-3
S13= U1I13* =-0,0534 - j0,5828 p.u.
S31= U3I31* = 0,0662 + j0,6018 p.u.
 1 y pl1− 3 1 
 + − 
I13   zsl1− 3 2 z −   U1  =  0,9615− j4,7877 − 0,9615+ j4,8077   1 
 =
sl1 3
y       
I 31 1 1 pl1− 3 U3   − 0,9615+ j4,8077 0,9615− j4,7877  1,1148 − 0,6073º
 − +            
 zsl1− 3 zsl1− 3 2  Y l1−3

Y l1−3
Trafo 3-2
S32= U3I32* =-0,5662 - j0,4805 p.u.
S23= U2I23* = 0,5854 + j0,5444 p.u.
 1 1 
 −
I 32  2
a zcc azcc  U 3   1,9113 − j6,3710 − 2,2936 + j7,6453 1,1148 − 0,6073º
 = 1  = 
1  U 2  − 2,2936 + j7,6453 2,7523 − j9,1743   
I 23  − 11,9938º 
 azcc            
zcc  Y t 3 −2
  
Y t 3 −2 14
 SISTEMAS ELÉCTRICOS DE POTENCIA I

Flujos de potencia

Línea/trafo Pij Qij Pji Qji Ppérdidas Qpérdidas

Línea 1-2 -0,1632 0,0438 0,1646 -0,0381 0,0014 0,0057

Línea 1-3 -0,0534 -0,5828 0,0662 0,6018 0,0128 0,0190

Trafo 3-2 -0,5662 -0,4805 0,5854 0,5444 0,0192 0,0639

Total 0,0334 0,0886

U1 = 1 0º p.u. U2 = 1 1,9938º p.u.
P1+ P2 + P3 Q1+ Q2 + Q3 + BU32
-0,2889

0,5063
0,2834

0,75
G G
-0,1632 0,1646

Nudo 1 0,0438 -0,0381 Nudo 2

Pd3 = 0.5 p.u.

Qd3 = 0.25 p.u.
54
-0
,0

58
5

0,
34

44
-0

54
,5
8

0,
28

Nudo 3 1,2:1
62
0,
06

6
,5
62

-0

05
0,
60

U3 = 1,1148 -0,6073º p.u.

,4
18

-0

Pd3 = 0.5 p.u.

0,6213

Qd3 = 0.5 p.u.

B = j0.5 p.u.

15