INSTITUTO TECNOLÓGICO DE NACIONAL DE MÉXICO

CAMPUS OAXACA

DEPARTAMENTO DE INGENIERÍA ELÉCTRICA Y

ELECTRÓNICA

MODELADO DE SISTEMAS ELÉCTRICOS DE POTENCIAS

TAREA 1 DE LA UNIDAD 3

FLUJOS DEMPOTENCIA

CARRERA DE INGENIERÍA ELÉCTRICA

ALUMNO: CAMPECHANO ANTONIO ZULEYMA / 21160261

OCTAVO SEMESTRE

DOCENTE: ING. SOMOHANO GONZALEZ BERNARDO.

FECHA DE ENTREGA : 01/08/2025.

 INTRODUCCION

El análisis de flujos de potencia es fundamental en la operación y planificación

de sistemas eléctricos de potencia. Este estudio permite conocer el
comportamiento de las tensiones, potencias activas y reactivas en los nodos,
así como los flujos de potencia en las líneas de transmisión. En esta tarea se
aplica el método iterativo de Gauss-Seide durante tres iteraciones, utilizando
la matriz de admitancias nodal (Ybus) como base de cálculo Los resultados
incluyen los flujos de potencia en las líneas y la evaluación de las pérdidas
totales en la red.
 TAREA DE FLUJOS DE POTENCIA

Para el sistema mostrado en la figura con los datos indicados, resolver para
tres iteraciónes usando el método de Gauss-Seidel con la matriz de
admitancias nodal (Ybus), indicando flujos de potencia en las líneas y las
pérdidas totales en la red.

Datos de línea
Resistencia Reactancia
Línea Entre nodos
(p.u) (p.u)
1 1–2 0.08 j0.4
2 1–3 0.12 j0.5
3 2–3 0.10 j0.4
4 3–4 0.0 j0.3

Datos de nodos
Bus
tipo voltaje Activa Reactiva
No.
1 Compensador 1.05∠0° - -
2 Carga 0.98∠0° -0.55 -j0.13
3 Carga 1.0∠0° -0.30 -j0.18
4 Generación 1.1∠0° 0.5 -

Pb = 100 MVA
KVb = 230 kV
 COMENZAMOS HACIENDO LA MATRIZ DE YBUS, CON LOS DATOS DE LA
LINEA
Cálculo paso a paso de la matriz Ybus

PRIMERO CALCULAMOS AFUERA DE LA DIAGONAL PRINCIPAL Y HACEMOS

EL CAMBIO DE SIGNO EN LA MATRIZ.

1
𝑌12 = 𝑌21 = = 0.4808 − j2.4038
0.08 + 𝑗0.4
1
𝑌13 = 𝑌31 = = 0.4539 − j1.8911
0.12 + j0.5
1
𝑌23 = 𝑌32 = = 0.5882 − j2.3529
0.10 + j0.4
1
𝑌43 = 𝑌34 = = 0.0000 − j3.3333
0.0 + j0.3

HACEMOS LOS CALCULOS PARA LA DIAGONAL PRINCIPAL

1 1
𝑌11 = 𝑌11 = + = 0.9346 − j4.2949
0.08 + 𝑗0.4 0.12 + j0.5
1 1
𝑌22 = 𝑌22 = + = 1.0690 − j4.7568
0.08 + 𝑗0.4 0.10 + j0.4
1 1 1
𝑌33 = 𝑌33 = + + = 1.0421 − j7.5773
0.12 + j0.5 0.10 + j0.4 0.0 + j0.3
1
𝑌44 = 𝑌44 = = 0.0000 − j3.3333
0.0 + j0.3

Resultado final: Matriz Ybus (compleja)

0.9346 - j4.2949 -0.4808 + j2.4038 -0.4539 + j1.8911 0.0000 + j0.0000

-0.4808 + j2.4038 1.0690 - j4.7568 -0.5882 + j2.3529 0.0000 + j0.0000
-0.4539 + j1.8911 -0.5882 + j2.3529 1.0421 - j7.5773 0.0000 + j3.3333
0.0000 + j0.0000 0.0000 + j0.0000 0.0000 + j3.3333 0.0000 - j3.3333
 Hallar las tensiones en cada nodo

Primera iteración
Como dato importante, tomamos el nodo 1 como el Slack, ya que es el nodo
del compensador.
𝑉1 ´ = 1.05 < 0° = 1.05 + 𝑗0
FORMULA
𝟏 𝑷𝟐 − 𝒋𝑸𝟐
𝑽𝟐 ´ = [ − (𝒀𝟐𝟏 𝑽𝟏 ° + 𝒀𝟐𝟑 𝑽𝟑 ° + 𝒀𝟐𝟒 𝑽𝟒 °)]
𝒀𝟐𝟐 (𝑽𝟐 °)∗

1 −0.55 + 𝑗0.13
𝑉2 ´ = [
1.06016 − 𝑗4.75674 0.98 − 𝑗0
− ((−0.48077 + j2.40385)(1.05 + 𝑗0)
+ (−0.58839 + j2.35289)(1 + j0)
+ (0.0000 + j0.0000)(1.1 + j0))] = 0.97334 − 𝑗0.10694

𝟏 𝑷𝟑 − 𝒋𝑸𝟑
𝑽𝟑 ´ = [ − (𝒀𝟑𝟏 𝑽𝟏 ° + 𝒀𝟑𝟐 𝑽𝟐 ° + 𝒀𝟑𝟒 𝑽𝟒 °)]
𝒀𝟑𝟑 (𝑽𝟑 °)∗

1 −0.30 + 𝑗0.18
𝑉3 ´ = [
1.04239 − j7.57728 1
− ((−0.45400 + j1.89106)(1.05 + 𝑗0)
+ (−0.58839 + j2.35289)(0.97334 − 𝑗0.10694)
+ (0.0000 + j3.3333)(1.1 + j0))] = 1.02242 − 𝑗0.07495

Tenemos que encontrar el valor de jQ4

Usando la formula
(𝑃4 − 𝑗𝑄4 ´) = (𝑌44 𝑉°4 + 𝑌41 𝑉´1 + 𝑌42 𝑉´2 + 𝑌43 𝑉´3 )(𝑉°4 ) ∗

(𝑃4 − 𝑗𝑄4 ´) = [(−𝑗3.3333)(1.1) + (0)(1.05) + (0)(0.97334 − 𝑗0.10694)

+ (𝑗3.3333)(1.02242 − 𝑗0.07495)](1.1) = 0.27481 − 𝒋𝟏𝟎. 𝟐𝟖𝟒𝟒𝟓

𝟏 𝑷𝟒 − 𝒋𝑸𝟒
𝑽𝟒 ´ = [ − (𝒀𝟒𝟏 𝑽𝟏 ° + 𝒀𝟒𝟐 𝑽𝟐 ° + 𝒀𝟒𝟑 𝑽𝟑 °)]
𝒀𝟒𝟒 (𝑽𝟒 °)∗
 1 0.5 − 𝑗10.28445
𝑉´4 = [
−j3.3333 1.1
− ((0)(1.05 + 𝑗0) + (0)(0.97334 − 𝑗0.1069)
+ (0.0000 + j3.3333)(1.02242 − j0.0795))]
= 1.09973 + 𝑗0.06141

Realizamos la segunda iteración

FORMULA
𝟏 𝑷𝟐 − 𝒋𝑸𝟐
𝑽𝟐 ´´ = [ − (𝒀𝟐𝟏 𝑽𝟏 ´ + 𝒀𝟐𝟑 𝑽𝟑 ´ + 𝒀𝟐𝟒 𝑽𝟒 ´)]
𝒀𝟐𝟐 (𝑽𝟐 ´)∗

1 −0.55 + 𝑗0.13
𝑉2 ´´ = [
1.06016 − 𝑗4.75674 0.97334 + 𝑗0.10694
− ((−0.48077 + j2.40385)(1.05 + 𝑗0)
+ (−0.58839 + j2.35289)(1.02242 − 𝑗0.07495)
+ (0.0000 + j0.0000)(1.09973 + 𝑗0.06141))]
= 0.97403 − 𝑗0.13773

𝟏 𝑷𝟑 − 𝒋𝑸𝟑
𝑽𝟑 ´´ = [ − (𝒀𝟑𝟏 𝑽𝟏 ´ + 𝒀𝟑𝟐 𝑽𝟐 ´´ + 𝒀𝟑𝟒 𝑽𝟒 ´)]
𝒀𝟑𝟑 (𝑽𝟑 ´)∗

1 −0.30 + 𝑗0.18
𝑉3 ´´ = [
1.04239 − j7.57728 1.02242 + 𝑗0.07495
− ((−0.45400 + j1.89106)(1.05 + 𝑗0)
+ (−0.58839 + j2.35289)(0.97403 − 𝑗0.13773)
+ (0.0000 + j3.3333)(1.09973 + 𝑗0.06141))]
= 1.02557 − 𝑗0.05511

Tenemos que encontrar el valor de jQ4

Usando la formula
(𝑃4 − 𝑗𝑄4 ´´) = (𝑌44 𝑉´4 + 𝑌41 𝑉´1 + 𝑌42 𝑉´´2 + 𝑌43 𝑉´´3 )(𝑉´4 ) ∗
 (𝑃4 − 𝑗𝑄4 ´´) = [(−𝑗3.3333)(1.09973 + 𝑗0.06141) + (0)(1.05)
+ (0)(0.97403 − 𝑗0.13773)
+ (𝑗3.3333)(1.02557 − 𝑗0.05511)](1.09973 − 𝑗0.06141)
= 𝒋𝟎. 𝟐𝟗𝟔𝟕𝟏

𝟏 𝑷𝟒 − 𝒋𝑸𝟒
𝑽𝟒 ´´ = [ − (𝒀𝟒𝟏 𝑽𝟏 ´ + 𝒀𝟒𝟐 𝑽𝟐 ´´ + 𝒀𝟒𝟑 𝑽𝟑 ´´)]
𝒀𝟒𝟒 (𝑽𝟒 ´)∗

1 0.5 − 𝑗0.29671
𝑉4 ´´ = [
−j3.3333 1.09973 − 𝑗0.06141
− ((0)(1.05 + 𝑗0) + (0)(0.97403 − 𝑗0.13773)
+ (0.0000 + j3.3333)(1.02557 − 𝑗0.05511))]
= 1.09865 + 𝑗0.08533

Realizamos la tercera iteración

FORMULA
𝟏 𝑷𝟐 − 𝒋𝑸𝟐
𝑽𝟐 ´´´ = [ − (𝒀𝟐𝟏 𝑽𝟏 ´ + 𝒀𝟐𝟑 𝑽𝟑 ´´ + 𝒀𝟐𝟒 𝑽𝟒 ´´)]
𝒀𝟐𝟐 (𝑽𝟐 ´´)∗

1 −0.55 + 𝑗0.13
𝑉2 ´´´ = [
1.06016 − 𝑗4.75674 0.97403 + 𝑗0.13773
− ((−0.48077 + j2.40385)(1.05 + 𝑗0)
+ (−0.58839 + j2.35289)(1.02557 − 𝑗0.05511)
+ (0.0000 + j0.0000)(1.09865 + 𝑗0.08533))]
= 0.97258 − 𝑗0.12537

𝟏 𝑷𝟑 − 𝒋𝑸𝟑
𝑽𝟑 ´´´ = [ − (𝒀𝟑𝟏 𝑽𝟏 ´ + 𝒀𝟑𝟐 𝑽𝟐 ´´´ + 𝒀𝟑𝟒 𝑽𝟒 ´´)]
𝒀𝟑𝟑 (𝑽𝟑 ´´)∗

1 −0.30 + 𝑗0.18
𝑉3 ´´´ = [
1.04239 − j7.57728 1.02557 + 𝑗0.05511
− ((−0.45400 + j1.89106)(1.05 + 𝑗0)
+ (−0.58839 + j2.35289)(0.97258 − 𝑗0.12537)
+ (0.0000 + j3.3333)(1.09865 + 𝑗0.08533))]
= 1.02623 − 𝑗0.04138
Tenemos que encontrar el valor de jQ4
Usando la formula
(𝑃4 − 𝑗𝑄4 ´´´) = (𝑌44 𝑉´´4 + 𝑌41 𝑉´1 + 𝑌42 𝑉´´´2 + 𝑌43 𝑉´´´3 )(𝑉´´4 ) ∗

(𝑃4 − 𝑗𝑄4 ´´´) = [(−𝑗3.3333)(1.09865 + 𝑗0.08533) + (0)(1.05)

+ (0)(0.97258 − 𝑗0.12537)
+ (𝑗3.3333)(1.02623 − 𝑗0.04138)](1.09865 − 𝑗0.08533)
= 𝒋𝟎. 𝟑𝟎𝟏𝟐𝟓

𝟏 𝑷𝟒 − 𝒋𝑸𝟒
𝑽𝟒 ´´´ = [ − (𝒀𝟒𝟏 𝑽𝟏 ´ + 𝒀𝟒𝟐 𝑽𝟐 ´´´ + 𝒀𝟒𝟑 𝑽𝟑 ´´´)]
𝒀𝟒𝟒 (𝑽𝟒 ´´)∗

1 0.5 − 𝑗0.30125
𝑉4 ´´´ = [
−j3.3333 1.09865 − 𝑗0.08533
− ((0)(1.05 + 𝑗0) + (0)(0.97258 − 𝑗0.12537)
+ (0.0000 + j3.3333)(1.02623 − 𝑗0.04138))]
= 1.09745 + 𝑗0.10068

Calculamos los flujos

Usamos la formula:
𝑽´𝒑𝒒
𝑷𝑷𝒒 − 𝒋𝑸𝒑𝒒 = 𝑽𝒑 ∗ 𝒚𝒑𝒒 (𝑽𝒑 − 𝑽𝒒 ) + 𝑽𝒑 ∗ 𝑽𝒑
𝟐

𝑽´𝒑𝒒
En esta ocasión vamos a omitir 𝑽𝒑 ∗ 𝑽𝒑 ya que no tenemos la
𝟐
admitancia

• 𝑷𝟏𝟐 − 𝒋𝑸𝟏𝟐 = 𝑽𝟏 ∗ 𝒚𝟏𝟐 (𝑽𝟏 − 𝑽𝟐 )

= (1.05)( 0.4808 − j2.4038)((1.05) − (0.97258 − 𝑗0.12537))

= 0.35552 − 𝑗0.13212
• ∗
𝑷𝟐𝟏 − 𝒋𝑸𝟐𝟏 = 𝑽𝟐 𝒚𝟐𝟏 (𝑽𝟐 − 𝑽𝟏 )

= (0.97258 + 𝑗0.12537)(0.4808 − j2.4038)((0.97258 − 𝑗0.12537)

− (1.05)) = −0.34508 + 𝑗0.07993
• ∗
𝑷𝟏𝟑 − 𝒋𝑸𝟏𝟑 = 𝑽𝟏 𝒚𝟏𝟐 (𝑽𝟏 − 𝑽𝟑 )
 = (1.05)(0.4539 − j1.8911)((1.05) − (1.02623 − 𝑗0.04138))
= 0.09349 − 𝑗0.02747
• ∗
𝑷𝟑𝟏 − 𝒋𝑸𝟑𝟏 = 𝑽𝟑 𝒚𝟑𝟏 (𝑽𝟑 − 𝑽𝟏 )

= (1.02623 + 𝑗0.04138)(0.4539 − j1.8911)((1.02623 − 𝑗0.04138)

− (1.05)) = −0.09246 + 𝑗0.02316
• ∗
𝑷𝟐𝟑 − 𝒋𝑸𝟐𝟑 = 𝑽𝟐 𝒚𝟐𝟑 (𝑽𝟐 − 𝑽𝟑 )

= (0.97258 + 𝑗0.12537)(0.5882 − j2.3529)((0.97258 − 𝑗0.12537)

− (1.02623 − 𝑗0.04138)) = −0.23253 + 𝑗0.04597
• ∗
𝑷𝟑𝟐 − 𝒋𝑸𝟑𝟐 = 𝑽𝟑 𝒚𝟑𝟐 (𝑽𝟑 − 𝑽𝟐 )

= (1.02623 + 𝑗0.04138)(0.5882 − j2.3529)((1.02623 − 𝑗0.04138)

− (0.97258 − 𝑗0.12537)) = 0.23837 − 𝑗0.06934
• ∗
𝑷𝟑𝟒 − 𝒋𝑸𝟑𝟒 = 𝑽𝟑 𝒚𝟑𝟒 (𝑽𝟑 − 𝑽𝟒 )

= (1.02623 + 𝑗0.04138)(0.0000 − j3.3333)((1.02623 − 𝑗0.04138)

− (1.09745 + 𝑗0.10068)) = −0.49577 + 𝑗0.22402
• ∗
𝑷𝟒𝟑 − 𝒋𝑸𝟒𝟑 = 𝑽𝟒 𝒚𝟒𝟑 (𝑽𝟒 − 𝑽𝟑 )

= (1.09745 − 𝑗0.10068)(0.0000 − j3.3333)((1.09745 + 𝑗0.10068)

− (1.02623 − 𝑗0.04138)) = 0.49577 − 𝑗0.30820

Multiplicamos cada flujo de potencia por la potencia base, que son 100 MVA

➢ 𝑷𝟏𝟐 − 𝒋𝑸𝟏𝟐 = (0.35552 − 𝑗0.13212)(100) = 35.552 − 𝑗. 13.212

➢ 𝑷𝟐𝟏 − 𝒋𝑸𝟐𝟏 = (−0.34508 + 𝑗0.07993)(100) = −34.508 + 𝑗7.993
➢ 𝑷𝟏𝟑 − 𝒋𝑸𝟏𝟑 = (= 0.09349 − 𝑗0.02747)(100) = 9.349 − 𝑗2.747
➢ 𝑷𝟑𝟏 − 𝒋𝑸𝟑𝟏 = (−0.09246 + 𝑗0.02316)(100) = −9.246 − 𝑗2.316
➢ 𝑷𝟐𝟑 − 𝒋𝑸𝟐𝟑 = (−0.23253 + 𝑗0.04597)(100) = −23.253 + 𝑗4.597
➢ 𝑷𝟑𝟐 − 𝒋𝑸𝟑𝟐 = (0.23837 − 𝑗0.06934)(100) = 23.837 − 𝑗6.934
➢ 𝑷𝟑𝟒 − 𝒋𝑸𝟑𝟒 = (−0.49577 + 𝑗0.22402)(100) = −49.577 + 𝑗22.402
➢ 𝑷𝟒𝟑 − 𝒋𝑸𝟒𝟑 = (0.49577 − 𝑗0.30820)(100) = 49.577 − 𝑗30.820
 Con los datos obtenidos realizamos una tabla

GENERACION CARGA FLUJO

BUS V MW MVVAR MW MVAR BUS MW MVAR

1 1.0500 < 0.00° 44.901 15.959 0 0 2 35.552 13.212

3 9.349 2.747

2 0.98062 <-7.34518° 0 0 55.0 13.0 1 -34.508 -7.993

3 -23.253 -4.597

3 1.02706 <-2.30904° 0 0 30.0 18.0 1 -9.246 -2.316

2 23.837 6.934

4 -49.577 -22.402

4 1.10205 < 5.24163° 50.0 30.820 0 0 3 49.577 30.820

Totales del sistema 94.901 46.779 85 31 1.731 16.405

REALIZAMOS LOS CALCULOS DE PERDIDAS EN EL SISTEMA

PERDIDIAS EN LAS LINEAS

o (𝑷𝟏𝟐 ) − (𝑷𝟐𝟏 ) = 35.552 −34.508=1.044 MW
o (𝑷𝟏𝟑 ) − (𝑷𝟑𝟏 ) = 9.349 − 9.246 = 0.103 MW
o (𝑷𝟑𝟐 ) − (𝑷𝟐𝟑 ) = 23.837 − 23.253 = 0.584 MW
o (𝑷𝟑𝟒 ) − (𝑷𝟒𝟑 ) = 49.577 − 49.577 = 0 MW

PERDIDAS TOTALES
𝑷𝒕𝒐𝒕𝒂𝒍 = 1.044 MW + 0.103 MW + 0.584 MW + 0 MW = 1.731MW.
REALIZAMOS EL DIAGRAMA
CONCLUSION

Mediante el método de Gauss-Seidel se logró estimar satisfactoriamente las

tensiones nodales y los flujos de potencia en el sistema eléctrico tras solo tres
iteraciones, lo que proporciona una aproximación aceptable para propósitos
académicos. Se observaron las pérdidas de potencia activa en las líneas, reflejando
la eficiencia energética del sistema. Esta tarea reafirma la utilidad del análisis de
flujos de carga para el diseño, operación y mejora del desempeño en redes
eléctricas.

ANEXAMOS LA TAREA, QUE TAMBIEN SE HIZO MANUALMENTE.