SISTEMAS DE POTENCIA

GAUSS SEIDEL Y NEWTON RAPHSON

SEGUNDO CORTE

PRESENTADO POR:
JOSE ARMANDO JAUREGUI ZIPACON

DOCENTE:
ING. JESUS ABELARDO VELAZCO OCHOA

UNIVERSIDAD DE PAMPLONA
FACULTAD DE INGENIERÍA Y ARQUITECTURA
INGENIERÍA ELÉCTRICA
CAMPUS VILLA DEL ROSARIO
2020-2
Dado el siguiente sistema de potencia MVA base=100

Se conoce la siguiente información

DATOS DE BARRA
BARRA |V| p.u α (°) Pg ( MW ) Q g ¿) P LOAD ( MW ) Q LOAD (MVAR )
1 1.0 0.0 ------- ------ 0.00 0.00
2 0.98517 ------ 0.0 ------ 80.00 60.0
3 1.0 ------- 150.0 ------- 200.00 150.0

DATOS DE LINEA
LINEA RESISTENCIA (p.u) REACTANCIA (p.u) SUSCEPTANCIA (p.u)
1-2 0.07448 0.03720 0.02000
1-3 0.01272 0.06360 0.03000
2-3 0.01008 0.05048 0.02500

Determinar:
1. Perfil de voltajes por Gauss Seidel y por Newton Raphson
2. Potencia reactiva que entrega el Capacitor
3. Perdidas de linea
4. Balance de potencias

SOLUCION
Identificación de barras
BARRA 1: barra de referencia V 1SP=1.0 ∠ 0 °
BARRA 2: barra PV V 2SP=0.98517
−80
P2SP =P G 2−PD 2=0.0−80= =−0.8
100
BARRA 3: barra PV V 3SP =1.0
−50
P3SP =P G 2−P D 3=150−200= =−0.5
100
Matriz de admitancias

Y 1−2=¿
1
=10.74575− j5.36711
0.07448+ j 0.03720

Y 1−3=¿
1
=3.02371− j15.11853
0.01272+ j 0.06360
Y 2−3=¿
1
=3.80401− j19.05023
0.01008+ j 0.05048

13.76946− j 20.46064 −10.74575+ j5.36711 −3.02371+ j15.11853

[
Y BARRA = −10.74575+ j5.36711 14.54976− j 24.39484 −3.80401+ j 19.05023
−3.02371+ j 15.11853 −3.80401+ j 19.05023 6.82772− j34.14126 ]
METODO GAUSS SEIDEL
sp r+1
1 Pi −Qi
V ir +1=
Y ii V i¿[ r+1
k≠ i ]
−∑ Y ik V k Qir +1=−I mag V i¿ ∑ Y ik V k
[ k ≠i ]
r +1 Vi
V icorregido =V i
V ir+1
Errores
Tensión:
∆ V i =¿ V ik+1−V ik ∨¿
Potencia:
n
sp
∆ Pi =P i −V i ∑ (Gik cos θik + Bik sin θ ik )V k
k=1
n
∆ Qi =Qisp −V i ∑ (Gik sin θ ik −B ik cos θik )V k
k=1
Potencias
Potencias inyectadas:
n
2
Pi=V i ∗Gii +V i ∑ V k (Gik cos ϴik ¿ + Bik sin ϴik ) ¿
k=1
j ≠i
n
Qi=−V i2∗B ii +V i ∑ V k (Gik cos ϴik ¿ + Bik sin ϴik ) ¿
k=1
j ≠i
Balance de potencias
Si=S Gi−S Di
Iteración No 1
Valores iniciales
V 10=1.0 ∠0 ° V 20 =0.98517 ∠0 ° V 30 =¿1.0∠ 0 °
BARRA 2
Q21=−I mag [ V 2¿ [ Y 21 V 1+Y 22 V 2+ Y 23 V 3 ] ]

Q 21=−I mag [ (0.98517 ∠0 °) [ (−10.74575+ j 5.36711 ) ( 1.0∠ 0 )+ ( 14.54976− j24.39484 )( 0.98517 ∠0 )+ (−3.8040

Q 21=−0.37858
1
V 2 1= ¿
Y 22

1 −0.8+ j0.37858
V 2 1= [
14.54976− j24.39484 0.98517− j0
−( (−10.74575+ j5.36711 ) ( 1.0 ∠ 0 ) + (−3.80401+ j 19.05023 )( 1.0

1 −0.8+ j0.37858
V 2 1= [
14.54976− j24.39484 0.98517− j 0
−(−14.54976+ j24.41734) ]
V 2=0.97442− j 0.01803=0.97458∠−1.05998 °
Se corrige V2 con la magnitud especificada y el ángulo calculado:
V 2 corr 1=0.98517 ∠−1.05998 °=0.98500− j 0.01822

Errores
∆ V 2 =¿ 0.98517−0.98500∨¿0.00017
∆ V 2 imag =¿ 0−0.01822∨¿ 0.01822
∆ ϴ2=¿1.05998−0∨¿ 1.05998
BARRA 3
Valores actualizados
V 10=1.0 ∠0 ° V 21 =0.98517 ∠−1.05998 ° V 30=¿1.0∠ 0 °
Q31=−I mag [ V 3¿ [Y 31 V 1+ Y 32 V 2+ Y 33 V 3 ] ]

Q 31=−I mag [ (1.0 ∠0 °) [ (−3.02371+ j 15.11853 ) ( 1.0∠ 0 )+ (−3.80401+ j19.05023 )( 0.98517 ∠−1.05998 ) + ( 6.827

Q 31=0.188899
1
V 3 1= ¿
Y 33

1 −0.5− j0.188899
V 3 1=
6.82772− j34.14126 [
1− j 0
−( (−3.02371+ j 15.11853 ) ( 1.0∠ 0 ) + (−3.80401+ j 19.05023 )( 0.

1 −0.5− j0.18889
V 3 1=
6.82772− j34.14126 [
1− j 0
−(−6.42348+ j33.95236) ]
V 31=0.99491− j 0.02546=0.99523∠−1.46613 °
Se corrige V3 con la magnitud especificada y el ángulo calculado:
V 3 corr 1=1.0 ∠−1.46613° =0.99967− j 0.025586

Errores
∆ V 2 =¿ 1−0.99967∨¿0.00033
∆ V 2 imag =¿ 0−0.025586∨¿ 0.025586
∆ ϴ2=¿1.46613−0∨¿ 1.46613

Iteración No 2
Valores actualizados
V 11=1.0 ∠0 ° V 21=0.98517∠−1.05998 ° V 31=¿1.0∠−1.46613°
BARRA 2
Q22=−I mag [ V 2¿ [ Y 21 V 1+Y 22 V 2+ Y 23 V 3 ] ]

Q 22=−I mag [ (0.98517 ∠1.05998° ) [ (−10.74575+ j 5.36711 ) (1.0 ∠ 0 )+ (14.54976− j 24.39484 ) ( 0.98517 ∠−1.05

Q 22=−I mag [ (0.98517 ∠1.05998° ) [−0.1741+ j 0.21432 ] ]

Q22=−¿0.20984
 1
V 2 2= ¿
Y 22

1 −0.8+ j0.20984
V 2 2= [
14.54976− j24.39484 0.98517 ∠ 1.05998
−( (−10.74575+ j5.36711 )( 1.0 ∠ 0 ) + (−3.80401+ j19.05023 ) (

1 −0.8+ j 0.207936
V 2 2= [
14.54976− j24.39484 0.98517 ∠ 1.05998
−(−14.061045+ j24.50844 ) ]
V 22=0.97392− j 0.03792=0.97465 ∠−2.22971°
Se corrige V2 con la magnitud especificada y el ángulo calculado:
V 2 corr 2=0.98517 ∠−2.22971 °=0.98442− j0.03833

Errores
∆ V 2 =¿ 0.98500−0.98442∨¿0.00058
∆ V 2 imag =¿ 0.01822−0.03833∨¿ 0.02011
∆ ϴ2=¿1.05998−2.22971∨¿ 1.16973
BARRA 3
Valores actualizados
V 11=1.0 ∠0 ° V 22=0.98517∠−2.22971 ° V 31=¿1.0∠−1.46613°

Q31=−I mag [ V 3¿ [Y 31 V 1+ Y 32 V 2+ Y 33 V 3 ] ]

Q 32=−I mag [ (1.0 ∠1.46613 ° ) [ (−3.02371+ j 15.11853 )( 1.0 ∠0 ) + (−3.80401+ j19.05023 ) ( 0.98517 ∠−2.22971 ) +

Q 32=−I mag [ ( 1.0 ∠1.46628 ° ) [ −0.10025− j 0.29186 ] ]

Q32=0.289036
1
V 3 2= ¿
Y 33

1 −0.5− j0.289036
V 3 2= [
6.82772− j34.14126 1.0 ∠ 1.46613° °
−( (−3.02371+ j 15.11853 ) ( 1.0 ∠ 0 ) + (−3.80401+ j 19.05023 ) ( 0.

1 −0.5− j0.29433
V 3 2= [
6.82772− j34.14126 1.0 ∠ 1.46628 °
−(−6.05211+ j34.01563) ]
V 32=0.99419− j 0.036743=0.994868∠−2.116558 °
Se corrige V3 con la magnitud especificada y el ángulo calculado:
V 3 corr 2=1.0 ∠−2.116558 °=0.999317− j0.036932

Errores
∆ V 3 =¿ 0.99967−0.999317∨¿0.00035
∆ V 3 imag =¿ 0.02559−0.036932∨¿ 0.0113
∆ ϴ3 =¿1.46628−2.116558∨¿ 0.6502

Iteración No 3
Valores actualizados
V 11=1.0 ∠0 ° V 22=0.98517∠−2.22971 ° V 32=1.0 ∠−2.116558 °
BARRA 2
Q23=−I mag [ V 2¿ [Y 21 V 1+ Y 22 V 2 +Y 23 V 3 ] ]

Q 23=−I mag [ (0.98517 ∠2.22971° ) [ (−10.74575+ j 5.36711 )( 1.0 ∠ 0 ) + ( 14.54976− j 24.39484 ) ( 0.98517 ∠−2.22

Q 23=−I mag [ (0.98517 ∠2.22971° ) [−0.43449− j 0.017426 ] ]

Q 23=¿0.044674
1
V 2 3= ¿
Y 22

1 −0.8− j 0.044674
V 2 3= [
14.54976− j24.39484 0.98517 ∠ 2.18782
−( (−10.74575+ j5.36711 )( 1.0 ∠0 ) + (−3.80401+ j19.05023 ) (

1 −0.8− j 0.033496
V 2 3= [
14.54976− j24.39484 0.98517 ∠ 2.18782
−(−13.84055+ j24.545325) ]
V 23=0.979524− j0.049925=0.980795 ∠−2.91776 °
Se corrige V2 con la magnitud especificada y el ángulo calculado:
V 2 corr 3=0.98517∠−2.91776 °=0.983893− j 0.050147
Errores
∆ V 2 =¿ 0.98445−0.98394∨¿ 0.00051
∆ V 2 imag =¿ 0.03761−0.04912∨¿ 0.01151
∆ ϴ2=¿2.18782−2.8581∨¿ 0.67028
BARRA 3
Valores actualizados
V 11=1.0 ∠0 ° V 23=0.98517 ∠−2.91776 ° V 32=1.0 ∠−2.116558 °

Q31=−I mag [ V 3¿ [Y 31 V 1+ Y 32 V 2+ Y 33 V 3 ] ]

Q 33=−I mag [ (1.0 ∠2.116558 ° ) [ (−3.02371+ j 15.11853 ) ( 1.0∠ 0 )+ (−3.80401+ j 19.05023 )( 0.98517 ∠−2.9177 ) +

Q 33=−I mag [ ( 1.0 ∠2.12563 ° ) [−0.27414− j0.32125 ] ]

Q33=0.326404
1
V 3 3= ¿
Y 33

1 −0.5− j0.326404
V 3 3= [
6.82772− j34.14126 1.0 ∠2.116558
−( (−3.02371+ j15.11853 )( 1.0 ∠ 0 ) + (−3.80401+ j19.05023 ) ( 0.

1 −0.5− j0.326404
V 3 3= [
6.82772− j34.14126 1.0 ∠2.116558
−(−5.83084 + j 34.04976) ]
V 33=0.997564− j0.044276=0.998546 ∠−2.54135°
Se corrige V3 con la magnitud especificada y el ángulo calculado:
V 3 corr 3=1.0 ∠−2.54135 °=0.99901− j 0.0443403

Errores
∆ V 3 =¿ 0.99931−0.99901∨¿0.00029
∆ V 3 imag =¿ 0.036932−0.0443403∨¿ 0.0074
∆ ϴ 3 =¿ 2.116558−2.54135∨¿ 0.4247

Iteración No 4
Valores actualizados
V 11=1.0 ∠0 ° V 23=0.98517 ∠−2.91776 °=V 33 =1.0 ∠−2.54135 °
BARRA 2
Q 24 =−I mag [ V 2¿ [ Y 21 V 1 +Y 22 V 2 +Y 23 V 3 ] ]

Q 24 =−I mag [(0.98517∠2.91776 ° ) [ (−10.74575+ j5.36711 )( 1.0 ∠ 0 ) + ( 14.54976− j 24.39484 ) ( 0.98517 ∠−2.91

Q23=−I mag [ (0.98517 ∠2.91776 °) [ −0.43449− j 0.017426 ] ]

Q 23=0.192291
1
V 24 = ¿
Y 22

1 −0.8− j 0.192291
V 24 = [
14.54976− j 24.39484 0.98517∠ 2.91776
−( (−10.74575+ j5.36711 ) ( 1.0 ∠0 ) + (−3.80401+ j 19.05023 )

1 −0.8− j 0.192291
V 24 = [
14.54976− j 24.39484 0.98517∠ 2.91776
−(−13.84055+ j24.545325) ]
V 24 =0.979749− j0.056353=0.98136 ∠−3.29193°
Se corrige V2 con la magnitud especificada y el ángulo calculado:
V 2 corr 4 =0.98517 ∠−3.29193° =0.983544− j 0.056572

Errores
∆ V 2 =¿ 0.98394−0.983544∨¿0.0004
∆ V 2 imag =¿ 0.04912−0.056572∨¿ 0.0074
∆ ϴ2=¿2.8581−3.29193∨¿ 0.4338
BARRA 3
Valores actualizados
V 11=1.0 ∠0 ° V 24 =0.98517 ∠−3.29193° V 33=1.0 ∠−2.54135 °

Q34 =−I mag [ V 3¿ [ Y 31 V 1 +Y 32 V 2+Y 33 V 3 ] ]

Q 34 =−I mag [ ( 1.0 ∠ 2.54135° ) [ (−3.02371+ j 15.11853 ) ( 1.0 ∠ 0 ) + (−3.80401+ j 19.05023 )( 0.98517 ∠−3.29193 ) +

Q 34 =−I mag [ ( 1.0 ∠ 2.54135° ) [ −0.27414− j 0.32125 ] ]

Q34 =0.356259 j
1
V 34 = ¿
Y 33

1 −0.5− j 0.356259
V 34 = [
6.82772− j 34.14126 1.0 ∠−2.54135
−( (−3.02371+ j15.11853 )( 1.0 ∠ 0 ) + (−3.80401+ j19.05023 ) ( 0.

1 −0.5− j 0.356259
V 34 = [
6.82772− j 34.14126 1.0 ∠2.54135
−(−5.83084+ j 34.04976) ]
V 34 =0.99808−0.048109 j=0.999238 ∠−2.75961°
Se corrige V3 con la magnitud especificada y el ángulo calculado:
V 3 corr 4 =1.0 ∠−2.75961 °=0.99884− j0.048145

Errores
∆ V 3 =¿ 0.99901−0.99884∨¿0.00017
∆ V 3 imag =¿ 0.0443403−0.048145∨¿ 0.0038
∆ ϴ 3 =¿ 2.54135−2.75961∨¿ 0.21826

Iteración No 5
Valores actualizados
1 4 4
V 1 =1.0 ∠0 ° V 2 =0.98517 ∠−3.29193° V 3 =1.0 ∠−2.75961
BARRA 2
Q 25=−I mag [ V 2¿ [Y 21 V 1+ Y 22 V 2 +Y 23 V 3 ] ]

Q 25=−I mag [ (0.98517 ∠3.29193 °) [ (−10.74575+ j 5.36711 ) (1.0 ∠ 0 )+ (14.54976− j24.39484 )( 0.98517 ∠−3.29

Q25=−I mag [ (0.98517 ∠ 3.29193° ) [−0.43449− j0.017426 ] ]

Q 25=0.273685 j
1
V 2 5= ¿
Y 22
 1 −0.8− j 0.273685
V 2 5= [
14.54976− j24.39484 0.98517 ∠ 3.29193
−( (−10.74575+ j5.36711 )( 1.0 ∠ 0 ) + (−3.80401+ j19.05023 ) (

1 −0.8− j 0.192291
V 2 5= [
14.54976− j24.39484 0.98517 ∠ 3.29193
−(−13.84055+ j 24.545325) ]
V 25=0.980998− j 0.060331=0.982851 ∠−3.519236°
Se corrige V2 con la magnitud especificada y el ángulo calculado:
V 2 corr 5=0.98517∠−3.519236 °=0.983312− j 0.060473

Errores
∆ V 2 =¿ 0.983544−0.983312∨¿0.00023
∆ V 2 imag =¿ 0.056572−0.060473∨¿ 0.0039
∆ ϴ2=¿3.29193−3.519236∨¿ 0.2273
BARRA 3
Valores actualizados
V 11=1.0 ∠0 ° V 25=0.98517 ∠−3.519236 ° V 34=1.0 ∠−2.75961

Q35=−I mag [ V 3¿ [ Y 31 V 1+ Y 32 V 2 +Y 33 V 3 ] ]

Q 35=−I mag [ (1.0 ∠2.75961 ° ) [ (−3.02371+ j15.11853 )( 1.0 ∠0 ) + (−3.80401+ j19.05023 ) ( 0.98517 ∠−3.519236

Q 35=−I mag [ ( 1.0 ∠2.75961 ° ) [ −0.27414− j 0.32125 ] ]

Q35=0.3700798
1
V 3 5= ¿
Y 33

1 −0.5− j0.3700798
V 3 5=
6.82772− j34.14126 [
1.0∠ 2.75961
−( (−3.02371+ j 15.11853 ) ( 1.0 ∠ 0 ) + (−3.80401+ j 19.05023 )( 0

1 −0.5− j0.3700798
V 3 5=
6.82772− j34.14126 [
1.0∠ 2.75961
−(−5.83084+ j34.04976) ]
V 35=0.998273−0.0504068 j=0.999544 ∠−2.890638 °
Se corrige V3 con la magnitud especificada y el ángulo calculado:
V 3 corr 5=1.0 ∠−2.890638 °=0.99873− j 0.05043

Errores
∆ V 3 =¿ 0.99884−0.99873∨¿ 0.00011
∆ V 3 imag =¿ 0.048145−0.05043∨¿ 0.00 22
∆ ϴ 3 =¿ 2.75961−2.890638∨¿ 0.131

Iteración No 6
Valores actualizados
V 11=1.0 ∠0 ° V 25=0.98517 ∠−3.519236 ° V 35=1.0 ∠−2.890638
BARRA 2
Q 26=−I mag [ V 2¿ [Y 21 V 1+ Y 22 V 2 +Y 23 V 3 ] ]

Q 26=−I mag [(0.98517 ∠3.519236 °) [ (−10.74575+ j 5.36711 ) ( 1.0∠ 0 ) + ( 14.54976− j24.39484 )( 0.98517 ∠−3.5

Q26=−I mag [ (0.98517 ∠3.519236 °) [ −0.752502− j0.28268 ] ]

Q 26=0.323469
1
V 2 6= ¿
Y 22

1 −0.8− j 0.323469
V 2 6= [
14.54976− j 24.39484 0.98517 ∠ 3.519236
−( (−10.74575+ j 5.36711 ) (1.0 ∠ 0 )+ (−3.80401+ j19.05023

1 −0.8− j 0.323469
V 2 6= [
14.54976− j 24.39484 0.98517 ∠ 3.519236
−(−13.584221+ j 24.584936) ]
V 26=0.98376∠−3.65714 °
Se corrige V2 con la magnitud especificada y el ángulo calculado:
V 2 corr 6=0.98517∠−3.65714 ° =0.98316− j 0.062839

Errores
∆ V 2 =¿ 0.983312−0.98316∨¿0.00015
∆ V 2 imag =¿ 0.060473−0.062839∨¿ 0.0023
∆ ϴ2=¿3.519236−3.65714∨¿ 0.1379
BARRA 3
Valores actualizados
V 11=1.0 ∠0 ° V 26 =0.98517 ∠−3.65714 ° V 35=1.0 ∠−2.890638

Q36=−I mag [ V 3¿ [ Y 31 V 1 +Y 32 V 2 +Y 33 V 3 ]]

Q 36=−I mag [ ( 1.0∠2.890638 ° ) [ (−3.02371+ j15.11853 )( 1.0 ∠0 )+ (−3.80401+ j19.05023 ) ( 0.98517 ∠−3.65714

Q 36=−I mag [ (1.0 ∠2.890638 ° ) [−0.469264− j0.355069 ] ]

Q36=0.378282
1
V 3 6= ¿
Y 33

1 −0.5− j 0.378282
V 3 6= [
6.82772− j34.14126 1.0∠ 2.890638
−( (−3.02371+ j 15.11853 ) ( 1.0∠ 0 ) + (−3.80401+ j 19.05023 )( 0.

1 −0.5− j 0.378282
V 3 6= [
6.82772− j34.14126 1.0∠ 2.890638
−(−5.566561+ j 34.08707) ]
V 36=0.999723∠−2.970115 °
Se corrige V3 con la magnitud especificada y el ángulo calculado:
V 3 corr 6=1.0∠−2.970115 °=0.99865− j0.051815

Errores
∆ V 3 =¿ 0.99873−0.99865∨¿0.00008
∆ V 3 imag =¿ 0.05043−0.051815∨¿ 0.0013
∆ ϴ 3 =¿ 2.890638−2.970115∨¿0.0794

Iteración No 7
Valores actualizados
V 11=1.0 ∠0 ° V 26 =0.98517 ∠−3.65714 ° V 36=1.0 ∠−2.970115
BARRA 2
Q 27=−I mag [ V 2¿ [Y 21 V 1+ Y 22 V 2 +Y 23 V 3 ] ]
Q 27=−I mag [(0.98517 ∠3.65714 ° ) [ (−10.74575+ j5.36711 )( 1.0 ∠0 ) + ( 14.54976− j 24.39484 ) ( 0.98517∠−3.65

Q27=−I mag [ (0.98517 ∠3.65714 ° ) [ −0.785732− j0.309574 ] ]

Q 27=0.353737
1
V 2 7= ¿
Y 22

1 −0.8− j 0.353737
V 2 7= [
14.54976− j24.39484 0.98517 ∠ 3.65714
−( (−10.74575+ j 5.36711 ) ( 1.0∠ 0 ) + (−3.80401+ j 19.05023 )(

1 −0.8− j 0.353737
V 2 7= [
14.54976− j24.39484 0.98517 ∠ 3.65714
−(−13.557561+ j 24.588854) ]
V 27=0.98431 ∠−3.74101°
Se corrige V2 con la magnitud especificada y el ángulo calculado:
V 2 corr 7=0.98517∠−3.74101 °=0.98307− j 0.06427

Errores
∆ V 2 =¿ 0.98316−0.98307∨¿0.00009
∆ V 2 imag =¿ 0.062839−0.06427∨¿ 0.0014
∆ ϴ2=¿3.65714−3.7410∨¿ 0.08386
BARRA 3
Valores actualizados
V 11=1.0 ∠0 ° V 27 =0.98517 ∠−3.74101 ° V 36 =1.0 ∠−2.970115

Q37=−I mag [ V 3¿ [ Y 31 V 1 +Y 32 V 2 +Y 33 V 3 ] ]

Q 37=−I mag [ ( 1.0∠2.970115 ° ) [ (−3.02371+ j 15.11853 ) ( 1.0∠ 0 ) + (−3.80401+ j 19.05023 )( 0.98517 ∠−3.74101

Q 37=−I mag [ (1.0 ∠2.970115 ° ) [ −0.48927− j 0.358404 ] ]

Q37=0.383275
1
V 3 7= ¿
Y 33
 1 −0.5− j0.383275
V 3 7= [
6.82772− j34.14126 1.0 ∠2.970115
−( (−3.02371+ j 15.11853 ) ( 1.0 ∠ 0 ) + (−3.80401+ j 19.05023 ) ( 0.

1 −0.5− j0.383275
V 3 7= [
6.82772− j34.14126 1.0 ∠ 2.970115
−(−5.538792+ j 34.090772) ]
V 37=0.99983∠−3.018455 °
Se corrige V3 con la magnitud especificada y el ángulo calculado:
V 3 corr 7=1.0 ∠−3.018455 °=0.99861− j 0.05265

Errores
∆ V 3 =¿ 0.99865−0.99861∨¿0.00004
∆ V 3 imag =¿ 0.051815−0.05265∨¿ 0.00083
∆ ϴ3 =¿ 2.970115−3.018455∨¿0.04834

Iteración No 8
Valores actualizados
V 11=1.0 ∠0 ° V 27 =0.98517 ∠−3.74101 ° V 37 =1.0∠−3.018455
BARRA 2
Q28=−I mag [ V 2¿ [Y 21 V 1+ Y 22 V 2 +Y 23 V 3 ] ]

Q 28=−I mag [(0.98517 ∠3.74101° ) [ (−10.74575+ j 5.36711 )( 1.0 ∠ 0 )+ ( 14.54976− j 24.39484 ) ( 0.98517 ∠−3.74

Q 28=−I mag [ (0.98517 ∠3.7410 °) [ −0.805973− j 0.325877 ] ]

Q28=0.372167
1
V 2 8= ¿
Y 22

1 −0.8− j 0.372167
V 2 8= [
14.54976− j24.39484 0.98517 ∠ 3.74101
−( (−10.74575+ j5.36711 )( 1.0 ∠ 0 ) + (−3.80401+ j19.05023 ) (

1 −0.8− j 0.372167
V 2 8= [
14.54976− j24.39484 0.98517 ∠3.74101
−(−13.541342+ j24.59122) ]
V 28=0.984646 ∠−3.792118 °
Se corrige V2 con la magnitud especificada y el ángulo calculado:
V 2 corr 8=0.98517∠−3.792118 ° =0.98301− j 0.06515

Errores
∆ V 2 =¿ 0.98307−0.98301∨¿0.00006
∆ V 2 imag =¿ 0.06427−0.06515∨¿ 0.00088
∆ ϴ2=¿3.74101−3.792118∨¿ 0.051

Cumple con el criterio de divergencia V 28=0.98517 ∠−3.792118 °

BARRA 3
Valores actualizados
V 11=1.0 ∠0 ° V 28 =0.98517 ∠−3.792118 ° V 37=1.0 ∠−3.018455

Q 38=−I mag [ V 3¿ [ Y 31 V 1 +Y 32 V 2 +Y 33 V 3 ] ]

Q 38=−I mag [ ( 1.0∠3.018455 ° ) [(−3.02371+ j15.11853 )( 1.0 ∠ 0 ) + (−3.80401+ j19.05023 ) ( 0.98517 ∠−3.79211

Q38=−I mag [ (1.0 ∠3.018455 ° ) [−0.48927− j0.358404 ] ]

Q 38=0.360416
1
V 3 8= ¿
Y 33

1 −0.5− j0.360416
V 3 8= [
6.82772− j34.14126 1.0∠ 3.018455
−( (−3.02371+ j15.11853 )( 1.0 ∠ 0 ) + (−3.80401+ j 19.05023 ) ( 0.

1 −0.5− j0.360416
V 3 8= [
6.82772− j34.14126 1.0∠ 3.018455
−(−5.521868+ j34.093008) ]
V 38=0.99917∠−3.03957 °
Se corrige V3 con la magnitud especificada y el ángulo calculado:
V 3 corr 8=1.0 ∠−3.03957 °=0.99859− j 0.053025

Errores
∆ V 3 =¿ 0.99861−0.99859∨¿0.00002
∆ V 3 imag =¿ 0.05265−0.053025∨¿ 0.000375
∆ ϴ3 =¿ 3.018455−3.03957∨¿ 0.021
Cumple con el criterio de convergencia
Valores finales obtenidos
V 18=1.0 ∠0 ° V 28 =0.98517 ∠−3.792118 ° V 38=1.0∠−3.03957 °

V 18=1.0 ∠0 ° V 28 =0.98517 ∠−3.86781° V 38=1.0 ∠−3.09259 °

b. Potencia reactiva que entrega el Capacitor

I 2sp=Y 21∗V 1 +Y 22∗V 2 +Y 23∗V 3

I 2sp=(−10.74575+ j 5.36711)∗(1.0 ∠0)+(14.54976− j24.39484)∗(0.98517 ∠−3.86781)+(−3.80401+ j 19.0

I 2sp=0.90671 ∠−157.26623

S2sp =V 2∗I 2¿ =( 0.98517 ∠−3.86781 )∗( 0.90671 ∠157.26623 )

S2sp=0.89326∠ 153.39842=−0.79870+ j 0.4

Q2sp=QG 2−QLOAD

Q G 2=Q 2sp+Q LOAD =0.4+ 0.6=1

Q G 2=Q C =1 POTENCIA REACTIVA DEL CAPACITOR

REACTANCIA DEL CAPACITOR

V 22 V 22 0.98517 2
QG 2=QC = X = = X C =0.970559
X C C QC 1

PERDIDAS DE LAS LINEAS

LINEA 1-2

V 1−V 2
I 1=I Z + I Ye = +V 1 Y e
Z
S1=V 1∗I 1¿

1.0 ∠0−0.98517 ∠−3.86781°

I 1= + 1.0∠0∗ j0.0 1=0.83172∠49.5019 3
0.07448+ j0.0372
S1=V 1∗I 1¿ =( 1.0 ∠0 )∗( 0.83172∠−49.50193 ) =0.83172∠−49.50193
V 2 −V 1
I 2=I Z + I Ye = +V 2 Y e
Z
0.98517∠−3.86781−1.0 ∠ 0
I 1= +0.98517 ∠−3.86781∗ j0.0 1=0.81631 ∠−131.3665
0.07448+ j0.0372
S2=V 2∗I 2¿ =( 0.98517 ∠−3.86781 )∗( 0.8 1631∠131.3665 )=0. 804204 ∠ 127.49869

S L12=S1 +S 2=( 0.83172 ∠−49.50193 ) +( 0. 804204 ∠127.49869)

S L12=0.05089 ∠ 6.27949

LINEA 2-3

V 2 −V 3
I 2=I Z + I Ye = +V 2 Y e
Z

S2=V 2∗I 2¿

0.98517∠−3.86781−1.0 ∠−3.09259 °
I 2= +0.98517 ∠−3.86781∗ j0.0 1 25
0.0 1008+ j 0.0 5048
I 2=0. 396051∠138.53761

S2=V 2∗I 2¿ =( 0.98517 ∠−3.86781 )∗( 0. 396051∠−138.53761 )

S2=0.390177 ∠−142.40542
V 3 −V 2
I 3=I Z + I Ye = +V 3 Y e
Z
1.0∠−3.09259−0.98517∠−3.86781
I 3= + 1.0∠−3.09259∗j 0.0 1
0.0 1008+ j 0.0 5048
I 3=0. 381281∠−38.52212

S3=V 3∗I 3¿=( 1.0 ∠−3.09259 )∗( 0. 381281∠38.52212 ) =0.38128 ∠35.42953

S L2 3=S2 + S3= ( 0.390177 ∠−142.40542 )+(0. 38128 ∠35.42953)

S L23=0.0 1708 ∠−84.89354

LINEA 1-3

V 38=1.0 ∠−3.0925 9

V1=1.0 ∠ 0

V 1−V 3
I 1=I Z + I Ye = +V 1 Y e
Z

S1=V 1∗I 1¿

1.0 ∠0−1.0 ∠−3.09259 °

I 1= + 1.0∠0∗ j0.0 1 5
0.0 1 272+ j 0.0 636
I 1=0. 83477 ∠ 10.77833

S1=V 1∗I 1¿ =( 1.0 ∠0 )∗( 0.83477 ∠−10.77833 )

S1=0.83477 ∠−10.77833
V 3 −V 1
I 3=I Z + I Ye = +V 3 Y e
Z
1.0∠−3.09259−1.0 ∠0
I 3= +1.0 ∠−3.09259∗ j0.0 1 5
0.0 1272+ j 0.0 636
I 3=0. 828887 ∠−171.2473

S3=V 3∗I 3¿=( 1.0 ∠−3.09259 )∗( 0. 828887 ∠ 171.2473 )=0. 828887 ∠168.15471

S L1 3=S1 + S3= ( 0.83477 ∠−10.77833 ) +(0.828887 ∠168.15471)

S L1 3=0.01 6569 ∠57.89264

POTENCIAS ENTREGADAS POR CADA GENERADOR

GENERADOR 1

I 1sp=Y 11∗V 1+ Y 12∗V 2+ Y 13∗V 3

I 2sp=(13.76946− j 20.46064)∗(1.0 ∠0)+(−10.74575+ j5.36711)∗(0.98517 ∠−3.86781)+(−3.02371+ j 15.1

I 1sp=1.57225 ∠30.10362

S1sp =V 1∗I 1¿ =( 1.0 ∠ 0 )∗(1.57225 ∠−30.10362 )

S1sp =1.57225 ∠−30.10362

SG 1=S 1+ S D 1= (1.57225 ∠−30.10362 )+ ( 0 )=1.57225 ∠−30.10362
GENERADOR 3
I 3sp=Y 31∗V 1 +Y 32∗V 2 +Y 33∗V 3

I 3sp=(−3.02371+ j15.11853 )∗(1.0 ∠0)+(−3.80401+ j 19.05023)∗(0.98517 ∠−3.86781)+(6.82772− j34.14

I 3sp=0.63528 ∠−145.0863

S3sp =V 3∗I 3¿=( 1.0 ∠−3.09259 )∗( 0.63528 ∠145.0863 )

S3sp =0.63528∠ 141.99371

SG 3 =S 3+ S D 3=( 0.63528 ∠141.99371) + ( 2+ j1.5 ) =2.41347 ∠ 51.59049

BALANCE DE POTENCIAS

∑ S G=¿ ∑ S D +∑ S L ¿
( 2.41347 ∠ 51.59049 ) + ( 1.57225∠−30.10362 ) + j 1=( 0.8+ j 0.6 ) + ( 2+ j 1.5 ) + ( 0.016569 ∠ 57.89264 ) + ( 0.01708 ∠
3.55∠36.32=3.55∠36.32

METODO NEWTON RAPHSON

ECUACIONES
Errores de potencia en barras
Barras PQ
n
ΔP i=Pisp−V i ∑ (Gik cos ϴik + Bik sinϴik )V k
k=1

n
sp
ΔQ i=Qi −V i ∑ (G ik sin ϴ ik−Bik cos ϴik )V k
k=1

Barras PV
 n
ΔP i=Pisp−V i ∑ (Gik cos ϴik + Bik sinϴik )V k
k=1

Términos y matriz jacobiana

ΔP = H N
[ ΔQ ] [ J L ][ ΔVΔϴ/V ]
n
H ii =V i ∑ (−Gik sinϴik +B ik cos ϴik ) V k =−Qi sp −V i2 Bii
j=1
k ≠i
H ik =V i (Gik sin ϴ ik −Bik cos ϴik ) V k
n
N ii =2 V i2 Gii + V i ∑ (Gik cos ϴik + Bik sin ϴik )V k =Pisp +V i2 G ii
k=1
k≠ i
N ik =V i (G ik cos ϴ ik + Bik sin ϴ ik )V k
n
J ii =V i ∑ (Gik cos ϴik + Bik sin ϴik )V k =Pisp −V i2 Gii
k=1
k ≠i
n
2 sp 2
J ik =−V i (G ik cos ϴik + Bik sin ϴ ik )V k Lii =−2 V i Bii +V i ∑ ( Gik sin ϴik−Bik cos ϴik )V k =Qi −V i B ii
k=1
k≠i
Lik =V i (Gik sin ϴik −Bik cos ϴik ) V k
SOLUCION

La barra 2 es una PQ por tanto se tendrán que determinar ΔP 2 , ΔP 3

N ecuaciones=2∗3−1−3=2

ΔP 2 H H 23 Δϴ2
[ ][
ΔP3
= 22
H 32 H 33 Δϴ3 ][ ]
1 iteración
Valores iniciales
V 1=1. 0 ∠ 0 V 2(0) =0.98517 ∠ 0 ϴ2(0)=1.0 ∠ 0

ΔP 2(1)=P2 sp−V 2 ¿

ΔP 3(1)=P 3sp−V 3 ¿

−80
P2SP =P G 2−PD 2=0.0−80= =−0. 8
100
−50
P3SP =P G 2−P D 3=150−200= =−0. 5
100
ΔP 2(1)=−0.8−0.98517 ¿

ΔP 2(1)=−0.587427 no cumple el criterio

ΔP 3(1)=P 3sp−V 3 ¿

ΔP 3(1)=−0.5−1.0 ¿

ΔP 3(1)=−0.556413no cumple con el criterio

H 22(1)=V 2 [(−G21 sin ϴ21+ B21 cos ϴ21 ) V 1+(−G23 sin ϴ23+ B23 cos ϴ23)V 3 ]
H 22(1)=0.98517 [ (−(−10.74575 ) sin ( 0−0 )+5.36711 cos ( 0−0 ) ) 1.0+ (−(−3.80401 ) sin ( 0−0 )+19.05023 cos ( 0−0

H 23=V 2 (G23 sinϴ23−B23 cos ϴ23) V 3

H 23=0.98517 ( (−3.80401 ) sin ( 0−0 )−19.05023 cos ( 0−0 ) ) 1.0=−18.767715

H 32=V 3 (G32 sin ϴ32−B32 cos ϴ32 )V 2

H 32=1.0 ( (−3.80401 ) sin ( 0−0 )−19.05023 cos ( 0−0 ) ) 0.98517=−18.767715

H 33(1)=V 3 [(−G31 sin ϴ31+ B31 cos ϴ31 ) V 1+(−G32 sinϴ32 + B32 cos ϴ32 )V 2 ]
H 33(1)=1.0 [ ( −(−3.02371 ) sin ( 0−0 ) +15.11853 cos ( 0−0 ) ) 1.0+ ( −(−3.80401 ) sin ( 0−0 ) +19.05023 cos ( 0−0 ) ) 0.98
Reemplazando valores en la matriz jacobiana e invertirla
−0.587427 = 24.05523 −18.767715 Δϴ2(1)
[ ][
−0.556413 −18.767715 33.886245 Δϴ3(1) ][ ]
−1 (1)
−0.587427 24.05523 = Δϴ2(1)
[ ][ −18.767715
−0.556413 −18.767715 33.886245 ] [ ]
Δϴ3
(1 )
−0.587427 0.0732 0.0405 = Δϴ2
[ ][
−0.556413 0.0405 0.0520 Δϴ3(1)] [ ]
Δϴ2(1) = −0. 0656
[ ][Δϴ3(1) −0.0527 ]
180 180
ϴ 2(1) =ϴ 2(0) +∆ ϴ 2(1) =0−0.0656 =−3.758603 °
π π
180 180
ϴ 3(1) =ϴ 3(0) + ∆ϴ 3(1) =0−0.0527 =−3.019487°
π π
Iteración 2
Actualizando valores
V 1=1.0 ∠ 0 V 2(1) =0.98517 ∠−3.758603 V 2(1 )=1.0 ∠−3.019487
ΔP 2(2)=P 2sp−V 2 ¿

ΔP 3(2)=P 3sp−V 3 ¿

−80
P2SP =P G 2−PD 2=0.0−80= =−0. 8
100
−50
P3SP =P G 2−P D 3=150−200= =−0. 5
100

ΔP 2(2)=−0.8−0.98517¿

ΔP 2(2)=−0. 021799 no cumple el criterio

(2) sp
ΔP 3 =P 3 −V 3 ¿

ΔP 3(2)=−0.5−1.0 ¿

ΔP 3(2)=−0.006642no cumple con el criterio

H 22(1)=V 2 [(−G21 sin ϴ21+ B21 cos ϴ21 ) V 1+(−G23 sin ϴ23+ B23 cos ϴ23)V 3 ]
H 22(1)=0.98517 [ (−(−10.74575 ) sin (−3.758603 )+ 5.36711 cos (−3.758603 ) ) 1.0+ (−(−3.80401 ) sin (−3.758603+

H 23=V 2 (G 23 sinϴ 23−B23 cos ϴ 23) V 3

H 23=0.98517 ( (−3.80401 ) sin (−3.758603+3.019487 )−19.05023cos (−3.758603+ 3.019487 ) ) 1.0=−18.99957

H 32=V 3 (G 32 sin ϴ 32−B32 cos ϴ 32 )V 2

H 32=1.0 ( (−3.80401 ) sin (−3.019487+3.758603 )−19.05023 cos (−3.019487+3.758603 )) 0.98517=−19.097715

H 33=V 3 [ (−G 31 sin ϴ 31 +B 31 cos ϴ 31 ) V 1 +(−G 32 sin ϴ 32+ B32 cos ϴ 32) V 2]
H 33=1.0 [ (−(−3.02371 ) sin (−3.019487 ) +15.11853 cos (−3.019487 ) ) 1.0+ (−(−3.80401 ) sin (−3.019487+ 3.7586
Reemplazando valores en la matriz jacobiana e invertirla
(2)
−0. 021799 = 23.299984 −18.99957 Δϴ2
[−0.006642 ][
−19.097715 33.752761 Δϴ3(2) ][ ]
−1 (2)
−0. 021799 23.299984 −18.99957 = Δϴ2(2)
[ ][
−0.006642 −19.097715 33.752761 ] [ ]
Δϴ3
(2)
−0. 021799 0.0797 0.0449 = Δϴ2
[ ][
−0.006642 0.0451 0.0550 ]
Δϴ3(2) [ ]
Δϴ2(2) = −0.002
[ ][Δϴ3(2) −0.0013 ]
 180 180
ϴ 2(2) =ϴ 2(1) +∆ ϴ 2(2) =−3.758603−0.002 =−3.87319
π π
180 180
ϴ 3(2) =ϴ 3(1) +∆ ϴ 3(2 ) =−3.019487−0.0013 =−3.093971
π π

Iteración 3
Actualizando valores
V 1=1.0 ∠ 0 V 2(2) =0.98517 ∠−3.8731 9 V 2(2 )=1.0∠−3.093971

ΔP 2(2)=P 2sp−V 2 ¿

ΔP 3(2)=P 3sp−V 3 ¿

−80
P2SP =P G 2−PD 2=0.0−80= =−0. 8
100
−50
P3SP =P G 2−P D 3=150−200= =−0. 5
100

ΔP 2(2)=−0.8−0.98517¿

ΔP 2(2)=−0. 00044 cumple con el criterio criterio de divergencia

ΔP 3(2)=P 3sp−V 3 ¿

ΔP 3(2)=−0.5−1.0 ¿

ΔP 3(2)=−0.00039 cumple con el criterio de divergencia

Por tanto, los valores son

V 1=1.0 ∠ 0 V 2(2) =0.98517 ∠−3.8731 9 V 2(2 )=1.0∠−3.093971