Sist 04
Sist 04
Fernando Ballesteros,2
ESCUELA TÉCNICA SUPERIOR DE 37700 BÉJAR
INGENIERÍA INDUSTRIAL ÁREA DE INGENIERÍA ELÉCTRICA
1 CONSIDERACIONES GENERALES
El estudio de sistemas eléctricos de potencia en situación estable analiza el comportamiento
de los sistemas en régimen permanente, es decir, en ausencia de fenómenos transitorios. Di-
chos fenómenos transitorios surgen de maniobras, descargas atmosféricas o alteraciones de
la red de algún otro tipo, como cortocircuitos, fallos a tierra o desconexiones bruscas de algún
elemento por actuación de sus protecciones.
El análisis de situaciones estables se ocupa del comportamiento de los sistemas en estado
estacionario, que es aquél caracterizado por el mantenimiento a largo plazo de las condicio-
nes iniciales. Esto no significa que no se evalúen los cambios en esas condiciones, sino que se
omiten las fluctuaciones transitorias para considerar solo las nuevas condiciones a que se da
lugar. Así se estudian las repercusiones para un sistema de cambios estacionarios como la
conexión o desconexión de líneas, de generadores o de cargas, pero no los fenómenos pasa-
jeros que ocurren justo en el momento de la variación.
El problema comienza planteándose como qué potencia deben inyectar los generadores a la
red en determinados nudos, para satisfacer la demanda absorbida por los consumidores en
esos u otros nudos. Se conocen la estructura de la red, los parámetros característicos de los
generadores, consumidores, líneas y transformadores, y los límites de funcionamiento de to-
dos ellos. Así, las fluctuaciones de tensión suelen restringirse a un cierto porcentaje alrededor
de la tensión nominal (por ejemplo, ±7%), la capacidad de transporte de las líneas está limi-
tada por un determinado valor característico de los materiales empleados, la potencia de ge-
neración tiene un máximo o un mínimo según el tipo de central de que se trate, la carga puede
oscilar entre ciertos límites, etc. Los cambios son suficientemente lentos como para que cada
situación constituya una solución estable, y la potencia se reparte de forma estacionaria a lo
largo de toda la red.
Se trata de una página web interactiva, gracias a las posibilidades que ofrecen HTML5 y Ja-
vascript, que llama a un programa escrito en R y que se ejecuta en el servidor de la página. La
interface entre el navegador y el intérprete de R es openCpu, de JEROEN OOMS
([Link]
A continuación siguen algunos ejemplos de cómo utilizarla y cómo interpretar sus resultados.
El modelo eléctrico de la línea es el de la fig. 2, que se ajusta a las condiciones de los modelos
de carga impuestas en el tema 2, y que eran:
- Una red de admitancias en paralelo con fuentes de intensidad independientes, salvo
la del nudo 1 que es una fuente de tensión cuya misión es fijar los potenciales de la
red y entregarle toda la potencia que esta necesita.
- ( )
Las fuentes de nudo entregan a la red una intensidad I 1 = S1 U 1 * , I 2 = S2 U 2 * , ( )
proporcional a la potencia que entregan a la red.
- Todas las fuentes tienen un terminal común.
Los parámetros iniciales se introducen en los campos siguientes:
- Admitancia en paralelo (fig. 3): cada uno de los dos fasores de las impedancias de 1
km que en el modelo de línea larga hemos denotado como y = G + BC i , en Ω/km. Su
parte real es la conductancia, de valor G = 0 Ω/km, y la imaginaria es la susceptancia
capacitiva, de valor BC = C ω , con C = 8,49275x10-9 F/km y ω = 2π 50 = 314 (para una
frecuencia de red de 50 Hz): BC = + 2,668076x10-6 i Ω/km.
Por defecto aparecen, al cargar la página, dos nudos y una línea que los une:
Fig. 4: Datos del problema por defecto, que aparecen con cada recarga de la página.
- Nudo (fig. 4): el número de cada nudo de la red, que debe ser único.
- Carga L (fig. 4): potencia aparente absorbida de la red por la fuente colocada en cada
nudo, en MVA. La carga del ejemplo en el nudo 2 significa una demanda de 20 MW,
con una componente reactiva de 15 MVA reactivos. En el modelo el nudo 1 no admite
carga y por eso su valor está atenuado. Tampoco admite generación, pese a represen-
tar una central, porque es el programa el que calcula su potencia.
- Gener G (fig. 4): potencia aparente entregada a la red por la fuente colocada en cada
nudo, en MVA. El nudo 1 no admite este dato por lo dicho antes, ya que por defecto es
un nudo que siempre entrega potencia y no la puede absorber. Y en este caso es nula
para el nudo 2, lo que lo convierte en demandante puro de energía.
- Tensión de nudo U (fig. 4): potencial de la fuente colocada en ese nudo, respecto del
terminal común a todas las fuentes de la red, en V. Se puede fijar arbitrariamente para
el nudo 1, y al resolver el problema aparecen las tensiones que impone la red a los
demás nudos. En este caso, el nudo 2 estará a 132 0º kV si el nudo 1 se sitúa en
142.251 2.893º kV = (142 070 + 7 180 i) kV:
Fig. 5. Solución del problema con los valores calculados de las tensiones,
intensidades y potencias aparentes de las fuentes de nudo.
- Intensidad de las fuentes del nudo I (fig. 5): intensidad, en A, que entregan las fuentes
al nudo de la red al que están conectadas. Su valor aparece al resolver el problema: la
fuente del nudo 1 entrega 93.1 −21.4º A y la del nudo 2 entrega 109.3 143.1º A.
- Potencia de las fuentes del nudo I (fig. 5): potencia aparente, en kVA, que entregan las
fuentes al nudo de la red en la que están conectadas. Es otro valor que aparece des-
pués de resolver el problema. La potencia que entrega el nudo 2 aparece confirmada:
-(20 000+15 000 i) kVA son los 20+15 i MVA de carga L introducidos en el campo de
la izquierda. La que entrega la fuente del nudo 1, que hasta ahora no se conocía, es la
necesaria para satisfacer la demanda más las pérdidas de la red. La potencia aparente
entregada a la red por la fuente del nudo 1 es 20 900 + 9 430i kVA, es decir, 20 900
kW de potencia activa, con componente reactiva de 9 340 kVA reactivos. Como la
única demanda activa son los 20 000 kW del nudo 2, las pérdidas de la red serán, re-
dondeando decimales, 20 900 + (-20 000) = 900 kW.
- Rendimiento de la red (%) (fig.
6): relación entre la potencia
activa retirada de la red (signo
negativo) y la entregada a la
red (signo positivo), en tantos
por ciento y con hasta tres de-
cimales de precisión: 95.696%.
Las pérdidas aparecen justo
encima con hasta cuatro deci-
males. Son resultados que se
obtienen al resolver el pro-
Fig. 6: Pérdidas, rendimiento y tabla de resultados para las blema.
ramas de la red
- Potencias entregadas a las líneas (fig. 6): tabla de la potencia aparente, en kVA, que
entregan las fuentes del nudo-fila al nudo-columna a través de la línea que los une, si
existe. En la diagonal no aparecen valores porque no tiene sentido que un nudo se
entregue potencia a sí mismo a través de la red.
- Intensidades de las líneas (fig. 6): intensidad que sale del nudo-fila hacia el nudo-co-
lumna, en A. No hay valor si la línea no existe, ni tampoco en la diagonal, por lo mismo
que en la tabla de potencias anterior.
SIGNOS ALGEBRAICOS:
La potencia que la fuente de la de-
recha entrega a la red (fig. 8) es
positiva cuando lo son su parte
real y su parte imaginaria, esto es: Fig. 8: Sentidos positi-
vos de la tensión, la in-
Sentr = vAB iBA * = Pentr + j Qentr tensidad y la potencia
entregada a la red de
Sentr > 0 ⇒ Pentr > 0 y Qentr > 0 una fuente del modelo.
Decir que la parte real y la imaginaria de la potencia aparente son positivas es lo mismo que
decir que su argumento está comprendido entre 0º y 90º. Puede verse que, en efecto, la po-
tencia que entrega el nudo 1 es positiva, ya que
Puede resultar llamativo el valor eficaz de las intensidades que entregan ambas fuentes a la
línea: un amperímetro situado junto al nudo 1, que entrega la potencia que absorben la línea
y la carga del nudo 2, marcará 93.1 A, mientras que otro colocado junto al nudo 2, marcaría
I j = I A + I jk ⎫⎪
( ) (
⎬ I = IA + IB − Ik ; Ik = IA + IB − I j
I jk = I B − I k ⎪ j
)
⎭
Teniendo en cuenta que las intensidades que se derivan a través de las capacidades de la línea
( I A e I B ) tienen un argumento de 90º en adelanto respecto a sus respectivas tensiones, la
ecuación anterior corresponde al siguiente diagrama fasorial, donde se aprecia claramente
que el módulo de I k es bastante mayor que el de I j .
Fig. 12. Si las potencias introducidas como dato en los nudos no coinciden
con las potencias calculadas para ellos, el problema no está resuelto.
Fig. 1. Ejemplo de hoja de cálculo para resolver sistemas de hasta 6 nudos (hasta 15 líneas).
De los dos factores anteriores, es el empleo de números complejos el que más dificulta el
cálculo, pues para resolver sistemas indeterminados se pueden utilizar simples hojas de
cálculo. Por este motivo, y para facilitar al alumno la posibilidad de resolver ejercicios por sí
mismo sin necesidad de utilizar aplicaciones informáticas no tan habituales, en adelante se
trabajará solo con números reales, de acuerdo con los planteamientos simplificados recogi-
dos en la sección 0, p. 94 ss. para el modelo simplificado, y en la sección 7, p. 104 ss. para el
modelo mejorado que trabaja con redes de dos puertas en "pi". Esto significa que el modelo
⎛ 6⎞ 6! 720
NL = ⎜ ⎟ = = = 15 líneas.
(
⎝ 2⎠ 2! 6 − 2 ! 48 )
Con este modelo se pueden hacer ejercicios sobre un pequeño ordenador en el que funcione
una hoja de cálculo dotada de la función de cálculos iterativos para resolver las referencias
circulares que surgen al programar la solución del sistema de ecuaciones del método de los
nudos basada en el método de GAUSS-SEIDEL.
La generación conectada en cada nudo se indica en las celdas C10:H10. La hoja permite la
asignación de cualquier potencia en las celdas amarillas, pero en la C10 se obtiene el resultado
del proceso iterativo correspondiente a la potencia de generación del nudo 1, para el que se
puede fijar la tensión en la celda C40. Fijada la tensión de la fuente, la variable que se calcula
como resultado es la intensidad, que a su vez determina la potencia de esa fuente. Por ello no
es dato (celda amarilla) sino resultado.
Las potencias siempre se introducen con signo positivo. La herramienta corrige interna-
mente, asignando signo positivo a la potencia inyectada por los generadores, y negativo a la
absorbida por los centros de demanda.
Cuando figure una "H" en la celda B11 la hoja calculará el estado de la red tomando los valores
de la máxima carga (D08:H08), y cuando figure "L" tomará el valor reducido (D09:H09). El
valor seleccionado en cada caso aparece en el rango de celdas D11:H11. De esta forma es po-
sible alternar entre los dos estados de la red rápidamente, para evaluar los límites de funcio-
namiento de ésta en una situación dada.
La configuración topológica de la red y los parámetros de las líneas se indican en el rango de
celdas de la fig. 3.
Los tramos están predefinidos en las columnas marcadas como INI y FIN, así como en la co-
lumna LÍNEA. Estos códigos no pueden modificarse.
Si la línea 12 existe, es decir, si los nudos 1 y 2 de la red están conectados, en la celda corres-
pondiente (en este caso la H14) se indicará el valor de la admitancia serie en Ω-1 de esa rama
bajo la leyenda y2. En las celdas anterior y siguiente (G14 e I14 respectivamente) aparecen
las impedancias transversales y1 e y3 en Ω-1. Estos valores son fruto de un cálculo, por lo que
el contenido de la celda no puede ser modificado.
Si la línea no existe, porque los nudos no están conectados, la celda amarilla correspondiente
quedará en blanco. La herramienta interpretará que la admitancia es cero (impedancia infi-
nita).
Anotar, por el contrario, un cero en el valor de la impedancia de una línea significaría que la
línea existe y que es un cortocircuito de impedancia nula, lo que causaría la inestabilidad de
la herramienta por admitancia infinita y el fracaso de la solución por el método iterativo. Por
eso la hoja está protegida frente a este descuido, y un valor cero no afecta a la impedancia de
la línea, que se considerará a pesar de ello infinita, y su admitancia nula.
Las líneas que parten de los nudos 2 a 6 pueden ser convertidas en transformadores regula-
dores de tensión indicando la relación de transformación en la columna B de la hoja de cálculo
y una 'longitud' equivalente a la impedancia serie del correspondiente modelo en 'pi'.
Si la relación de transformación se deja en blanco o se pone a 1, la fila vuelve a corresponder
a una línea normal.
son los opuestos de los valores que figuran en la matriz (en este caso el valor tomado auto-
máticamente de la celda H14 y cambiado de signo - ver fig. 5 -).
El resultado del cálculo se muestra en el bloque de celdas de la fig. 6:
Fig. 7. Cuando no se parte de un conjunto de valores que fueran solución de un análisis anterior, la iteración no
tiene éxito y aparecen errores irreversibles como el que se muestra.
La iteración está bien terminada cuando al forzar nuevos cálculos sobre los mismos paráme-
tros no se percibe que los valores cambien. Si, por el contrario, el programa aún muestra in-
tentos de búsqueda de una solución, es que la iteración no ha concluido. Revísense las condi-
ciones de partida, en especial la proximidad de la situación inicial a la solución buscada. Si la
distancia fuera mucha, ajústense los parámetros gradualmente para ir aproximando la
La suma total de las potencias de fuente son las pérdidas totales en la red, que se indican en
la celda I51. El rendimiento de la red aparece en la celda I52.
La celda C53 es un parámetro que evita el error de división por cero en el cálculo automático.
Este error se produciría al calcular la red cada vez que se llegase a dos nudos no unidos por
una línea, donde la admitancia sería cero. La hoja se ha programado para que dicha admitan-
cia nunca sea exactamente cero, sino al menos el valor indicado en C53, que al ser muy pe-
queño, no altera apreciablemente el resultado pero evita el error.
Por último, conocidas las potencias entregadas a las líneas y las tensiones de las fuentes en
los nudos, se calculan las intensidades de cada línea en una matriz idéntica a la de las poten-
cias entregadas a las líneas, según la expresión propia de redes trifásicas
Pjk
I jk =
3Uj
Los valores obtenidos para cada rama son comparados automáticamente con el parámetro
de intensidad máxima introducido en la celda B16, y se resaltan aquellos valores que vulne-
ran el límite prefijado. En la fig. 9 la línea 1-2 soporta más del límite y por ello su intensidad
aparece en número rojo sobre fondo naranja.
13 Agradecimiento al profesor D. Roberto Carlos Redondo Melchor por la ayuda prestada para poner a punto y depurar.
1 ⎛ Sj ⎞
* k< j k =n
U (j m+1) = ⎜ * ( m) − ∑ Y jkU k − ∑ Y jkU k( m) ⎟
( m+1)
Y jj ⎝ U j k =1 k> j ⎠
( )
n
I j = S j U j * = ∑U k Y jk
k =1
⎛ Uk − U j ⎞
S kj = U k I kj * = U k ⎜ + ytkjU k ⎟ *
⎝ Z jk ⎠
Fig. 9. Línea en proceso de duplicación. Fig. 10. Protecciones en un tramo de línea durante el
El apoyo admitirá tres cables más. tendido de los nuevos cables sobre una carretera.
- Si la nueva línea no es eficiente, mejorará a pesar de ello la tensión del nudo de llegada
y también el rendimiento de la red, pero en tan baja proporción que la rentabilidad
de su inversión quedará en entredicho.
14 Definición de MPE de acuerdo a la definición del Reglamento (UE) 2016/631 que establece un código de red sobre requisitos
de conexión de generadores a la red: “una unidad o un conjunto de unidades que genera electricidad, que está conectado de
forma no síncrona a la red o que está conectado mediante electrónica de potencia, y que además dispone de un solo punto de
conexión a una red de transporte, una red de distribución, incluidas las redes de distribución cerradas, o un sistema HVDC”.
los elementos no afectados por el incidente y considerar si las protecciones automáticas de-
ben permitirlas o no. Con ello se puede configurar un sistema muy robusto frente a gran tipo
de sucesos, en beneficio de la estabilidad del sistema y de los usuarios suministrados.
5 PROBLEMAS RESUELTOS
En las páginas que siguen se plantean y resuelven ejercicios de análisis simplificado de redes
malladas de hasta 6 nudos con líneas puramente resistivas, empleando la herramienta de
cálculo que se ha explicado en la sección 2, p. 116.
1. Sea la red de 45 kV de la figura en la que se han representado 6 nudos y las líneas que
actualmente los vinculan. Los valores en los recuadros son las longitudes de cada una (km) y
abajo se indica la potencia que se extrae de algunos de ellos. Suponiendo una impedancia
serie de 5 mΩ/km, una admitancia transversal de 5.50x10-8 S/km, y haciendo las suposicio-
nes habituales se pide:
a) Resolver la red, calculando las tensiones de cada nudo y la potencia transportada por
cada línea.
b) Hallar la potencia que se debe inyectar en el nudo 1 y el rendimiento de la red.
c) Calcular la tensión que ha de mantenerse en el nudo 1 para que la tensión de todos
los nudos quede dentro del margen del ±1,50% de la tensión nominal de la red. Recalcular los
parámetros anteriores si fuera necesario.
Solución:
La matriz de admitancias de nudo es:
1 2 3 4 5 6
1 0,004667 0,0 -0,004 0,0 0,0 -0,000667
2 0,0 0,002 0,0 -0,002 0,0 0,0
3 -0,004 0,0 0,004 0,0 0,0 0,0
4 0,0 -0,002 0,0 0,006 0,0 -0,004
5 0,0 0,0 0,0 0,0 0,002 -0,002
6 -0,000667 0,0 0,0 -0,004 -0,002 0,006667
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,004667 0,002 0,004 0,006 0,002 0,006667
TENSIÓN NUDO (V) 45 000 44 167 44 933 44 201 44 184 44 235
INT. FTES. DE NUDO (A) 448,9 -39,2 -154,2 -39,2 -58,8 -156,6
POT. FTES. DE NUDO (kW) 34 985 - 3 000 - 12 000 - 3 000 - 4 500 - 12 000
Luego las potencias entregadas por cada nudo a cada rama valen, en kilovatios:
1 2 3 4 5 6
1 0 12 023 0 0 22 962
2 0 0 - 3 000 0 0
3 - 12 000 0 0 0 0
4 0 3 013 0 0 - 6 013
5 0 0 0 0 - 4 500
6 - 22 539 0 0 6 023 4 516
Estos datos arrojan unas pérdidas totales de 485 kW, que comparadas con la potencia inyec-
tada en el nudo 1 de 34 985 kW, significan un rendimiento de la red del 98,613%.
La tensión de todos los nudos es menor que la del nudo 1, y la más pequeña se da en el nudo
2, donde se calculan 44 167 V. La diferencia respecto de la tensión nominal, a la que se ha
supuesto el nudo 1, es del -1,85%, fuera del margen, por tanto, del ±1,50% deseado.
Para resolver esta situación, la tensión del nudo 1 debe subir 155 V, de manera que la nueva
solución de la red pasaría a ser
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,004667 0,002 0,004 0,006 0,002 0,006667
TENSIÓN NUDO (V) 45 155 44 325 45 088 44 359 44 342 44 393
INT. FTES. DE NUDO (A) 447,3 -39,1 -153,7 -39,0 -58,6 -156,1
POT. FTES. DE NUDO (kW) 34 983 - 3 000 - 12 000 - 3 000 - 4 500 - 12 000
Ahora la caída de tensión en el nudo 2 es justo del -1,50% que es ya un valor admisible.
El nuevo rendimiento se extrae de la matriz de potencias entregadas a las líneas desde cada
nudo en kW, que es la siguiente:
1 2 3 4 5 6
1 0 12 023 0 0 22 959
2 0 0 - 3 000 0 0
3 - 12 000 0 0 0 0
4 0 3 013 0 0 - 6 013
5 0 0 0 0 - 4 500
6 - 22 539 0 0 6 023 4 516
Las pérdidas son ahora 483 kW, por lo que el rendimiento ha subido ligeramente, hasta el
98,620%.
2. Con el paso del tiempo los aisladores de las líneas se ven cubiertos con una fina capa de
suciedad procedente de la fuerte contaminación ambiental, lo que unido a la humedad por las
lluvias hace que haya que considerar una admitancia transversal 10 veces mayor, es decir,
que pasa de 5,50x10-8 S/km a 5,50x10-7 S/km. Calcular ahora las nuevas pérdidas y el rendi-
miento de la red en las condiciones del ejercicio 1.
Solución:
La matriz de admitancias de nudo es ahora distinta (no se representa), y aunque las diferen-
cias afectan a los últimos decimales porque los números son muy pequeños, determinan una
solución de la red diferente:
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,004667 0,002 0,004 0,006 0,002 0,006667
TENSIÓN NUDO (V) 45 169 44 325 45 103 44 359 44 343 44 394
INT. FTES. DE NUDO (A) 454,9 -39,1 -153,6 -39,0 -58,6 -156,1
POT. FTES. DE NUDO (kW) 35 588 - 3 000 - 12 000 - 3 000 - 4 500 - 12 000
No solo hay que elevar la tensión del nudo 1 algo más para mantener la del nudo 2 en el mí-
nimo aceptable, sino que las potencias entregadas a las líneas han cambiado también:
1 2 3 4 5 6
1 0 12 074 0 0 23 515
2 0 0 - 3 000 0 0
3 - 12 000 0 0 0 0
4 0 3 111 0 0 - 6 111
5 0 0 0 0 - 4 500
6 - 22 783 0 0 6 169 4 614
Estos datos arrojan unas pérdidas totales de 1088 kW, que comparadas con la potencia in-
yectada en el nudo 1 de 35 588 kW, significan un rendimiento de la red del 96,942%.
Solución:
Utilizaremos el valor de la admitancia transversal de 5,50x10-7 S/km del ejercicio 2.
Se observa que el nudo de menor tensión deja de ser el 2 y pasa a ser ahora el 5, aunque la
caída de tensión es de solo el -0,71% de la nominal de la red. Nótese que estos valores se han
obtenido recuperando, para el nudo 1, dicha tensión nominal de 45 kV, pues ya no es necesa-
rio elevarla para mantener los niveles dentro del ±1,50%.
El rendimiento se obtiene de la matriz de potencias entregadas a las líneas en esta situación,
en kW:
1 2 3 4 5 6
1 15 121 12 074 0 0 8 168
2 - 15 037 0 12 037 0 0
3 - 12 000 0 0 0 0
4 0 - 11 891 0 0 8 891
5 0 0 0 0 - 4 500
6 - 7 789 0 0 - 8 826 4 615
Para la opción b), es decir, la línea entre 3 y 5 de 50 km de longitud también, se obtienen los
siguientes resultados:
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,004667 0,002 0,008 0,006 0,006 0,006667
TENSIÓN NUDO (V) 45 000 44 613 44 857 44 648 44 782 44 682
INT. FTES. DE NUDO (A) 454,4 -38,8 -154,4 -38,8 -58,0 -155,1
POT. FTES. DE NUDO (kW) 35 416 - 3 000 - 12 000 - 3 000 - 4 500 - 12 000
1 2 3 4 5 6
1 0 25 709 0 0 9 706
2 0 0 - 3 000 0 0
3 - 25 572 0 0 13 572 0
4 0 3 112 0 0 - 6 112
5 0 0 - 13 494 0 8 994
6 - 9 307 0 0 6 171 - 8 864
Ahora las pérdidas son de 916 kW, por lo que el rendimiento es peor que antes, del 97,415%,
y la opción preferible es la a).
4. Con el paso de los años, la red de los problemas anteriores se ve completada con la línea
1-2 y la línea 3-5, pero debe hacer frente a la nueva potencia suministrada que se indica más
abajo. Resolver la red y determinar el margen de regulación aplicable al nudo 1 para que las
tensiones de la red se mantengan dentro del límite del ±1,50%.
Nudo 2: 10.5 MW Nudo 3: 12.0 MW Nudo 4: 25.0 MW Nudo 5: 4.5 MW Nudo 6: 24.0 MW
Solución:
La mínima tensión admisible es un 1,50% inferior a la nominal de 45 kV, es decir, 675 V me-
nos ó 44 325 V. Esta condición se da en el nudo 4 cuando la tensión del nudo 1 es de 44 769 V:
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,008667 0,006 0,008 0,006 0,006 0,006667
TENSIÓN NUDO (V) 44 769 44 582 44 598 44 325 44 494 44 338
INT. FTES. DE NUDO (A) 997,8 -136,0 -155,3 -325,6 -58,4 -312,5
POT. FTES. DE NUDO (kW) 77 369 - 10 500 - 12 000 - 25 000 - 4 500 - 24 000
1 2 3 4 5 6
1 33 630 30 708 0 0 13 032
2 - 33 434 0 22 934 0 0
3 - 30 535 0 0 18 535 0
4 0 - 22 694 0 0 - 2 306
5 0 0 - 18 438 0 13 938
6 - 12 580 0 0 2 361 - 13 781
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,008667 0,006 0,008 0,006 0,006 0,006667
TENSIÓN NUDO (V) 45 675 45 491 45 507 45 240 45 405 45 252
INT. FTES. DE NUDO (A) 978,1 -133,3 -152,2 -319,1 -57,2 -306,2
POT. FTES. DE NUDO (kW) 77 377 - 10 500 - 12 000 - 25 000 - 4 500 - 24 000
1 2 3 4 5 6
1 33 630 30 709 0 0 13 038
2 - 33 437 0 22 937 0 0
3 - 30 539 0 0 18 539 0
4 0 - 22 698 0 0 - 2 302
5 0 0 - 18 441 0 13 941
6 - 12 578 0 0 2 359 - 13 782
De manera que el margen de regulación para mantener la tensión de la red en los límites
deseados, actuando exclusivamente desde el nudo 1, es de
45 675 - 44 769 = 906 V (2.02% de la tensión de red).
Se comprueba que, en efecto, la línea 1-2 supera su capacidad, marcada por el límite térmico
de 426 A. Se verifica, además, que la línea está más sobrecargada en el origen, junto al nudo
1, que en el final, donde la intensidad es ya ligeramente menor (433,6994 - 432,9901 = 0,709
A de diferencia). Este fenómeno es debido a la influencia de las corrientes de fuga a lo largo
de la línea, introducidas en el modelo mediante las admitancias transversales del esquema en
"pi". Pero recuérdese que en una red real (véase el ejemplo de cálculo fasorial de una línea en
la sección 2.1, p. 117) ocurre justo al contrario, y por la misma razón.
A medida que sube la tensión de red la línea soporta algo menos de intensidad, de manera
que, cuando la tensión del nudo 1 está en 45 578 V, del nudo 1 salen hacia el nudo 2 exacta-
mente 426 A:
1 2 3 4 5 6
1 426,0000 389,0 0,000 0,0 165,1
2 -425,2778 0,0 291,7 0,0 0,0
3 -388,3 0,0 0,0 235,7 0,0
4 0,000 -290,3 0,0 0,0 -29,5
5 0,0 0,0 -235,0 0,0 177,6
6 -160,8 0,0 0,0 30,2 -176,2
Efectivamente, la línea 1-2 está sobrecargada cuando la red soporta la carga máxima y la ten-
sión del nudo 1 cae por debajo de esos 45 578 V.
Esa sobrecarga también se produce cuando la red suministra el 98,2% de la carga máxima
estadística y la tensión es la mínima posible (44 761 V en el nudo 1 para que el nudo 4 no baje
de 44 325 V).
b) Para simular en la herramienta de cálculo que la línea 1-2 se duplica puede multipli-
carse por dos la admitancia en serie (y2) de esa línea y multiplicarse también por dos sus
admitancias en paralelo (y1 e y3). O bien puede dividirse por dos su longitud, lo que es equi-
valente. Esto segundo es lo que vamos a hacer.
Con la nueva configuración la tensión más baja a la que puede trabajar la red cuando
Nótese que al variar las características de la línea 1-2 varía la matriz de admitancias de la red
y las admitancias propias de los nudos 1 y 2 que se ven en la tabla de arriba.
La nueva matriz de intensidades hacia las líneas, para esa tensión más baja, es ahora:
1 2 3 4 5 6
1 470,0340 374,1 0,000 0,0 153,4
2 -469,6794 0,0 333,8 0,0 0,0
3 -373,4 0,0 0,0 217,9 0,0
4 0,000 -332,4 0,0 0,0 6,8
5 0,0 0,0 -217,2 0,0 158,8
6 -149,1 0,0 0,0 -6,1 -157,4
Solución:
La matriz de admitancias de nudo de esta red es:
1 2 3 4 5 6
1 0,002486 -0,000816 -0,0008 -0,00087 0,0 0,0
2 -0,000816 0,002066 0,0 -0,00125 0,0 0,0
3 -0,0008 0,0 0,004988 0,0 -0,003077 -0,001111
4 -0,00087 -0,00125 0,0 0,004786 0,0 -0,002667
5 0,0 0,0 -0,003077 0,0 0,004327 -0,00125
6 0,0 0,0 -0,001111 -0,002667 -0,00125 0,005028
En la condición de carga máxima la tensión de la red no debería ser, en ningún punto, inferior
en más de un 1,0% a la nominal de la red, es decir, no bajar de 65 340 V. Pero en condiciones
de carga mínima, ningún punto debería superar la tensión en un 1,0%, esto es, pasar de 66
660 V.
Para la carga máxima, la tensión del nudo 1 puede oscilar entre el máximo de 66 660 V y bajar
hasta 66 208 V. Por debajo de ese valor la tensión en el nudo 6 vulneraría el límite del -1.0%:
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,002486 0,002066 0,004988 0,004786 0,004327 0,005028
TENSIÓN NUDO (V) 66 208 65 371 65 447 65 405 65 361 65 340
INT. FTES. DE NUDO (A) 1149,4 -418,2 -130,7 -278,1 -135,9 -183,6
POT. FTES. DE NUDO (W) 131 806 - 47 350 - 14 820 - 31 500 - 15 390 - 20 780
Para la carga mínima (se observa que es un 12% de la máxima), las tensiones a las que trabaja
la red tienden a ser las más altas, luego hay que verificar cuándo se sobrepasa el límite del
+1.0% en algún nudo. Evidentemente, al ser una red alimentada solo desde el nudo 1, basta
De forma que es sí es posible gestionar correctamente esta red simplemente manteniendo fija
la tensión en el nudo de generación en cualquier valor entre 66 660 V y 66 208 V, pues nunca
se vulnerarán los límites bajo cualquiera de los dos extremos de carga estudiados y, conse-
cuentemente, tampoco bajo cualquier caso intermedio.
7. Obtener las máximas intensidades que soportan las ramas de la red para un nuevo flujo
de cargas en el que se debe tener en cuenta un incremento de la carga máxima en el nudo 4
de 10 MW, quedando la mínima en 4,98 MW. Determinar si alguna supera el límite térmico
del conductor, establecido en 426 A, y en qué condiciones lo hace.
Solución:
Las intensidades serán máximas cuando las tensiones sean lo más pequeñas posible. Pero el
valor mínimo lo fija el umbral de -1.0% de variación respecto de la tensión nominal, que debe
respetarse en todos los nudos. De manera que primero hay que encontrar cuál es la tensión
mínima admisible de los nudos de la red en las condiciones de carga pedidas:
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,002486 0,002066 0,004988 0,004786 0,004327 0,005028
TENSIÓN NUDO (V) 66 276 65 389 65 466 65 390 65 375 65 340
INT. FTES. DE NUDO (A) 1237,6 -418,1 -130,7 -366,4 -135,9 -183,6
POT. FTES. DE NUDO (W) 142 064 - 47 350 - 14 820 - 41 500 - 15 390 - 20 780
Con la tensión del nudo 1 a 66 276 V la caída de tensión en el nudo 6, el más desfavorable, es
del -1.0%. Ahora las potencias en kW entregadas a las líneas desde cada nudo:
1 2 3 4 5 6
1 48 007 42 967 51 090 0 0
2 - 47 306 0 - 44 0 0
3 - 42 383 0 0 18 344 9 219
4 - 50 353 81 0 0 8 771
5 0 0 - 18 303 0 2 913
6 0 0 - 9 159 - 8 747 - 2 874
Se comprueba que con la red trabajando a plena carga y las tensiones al nivel más bajo posi-
ble, la intensidad de la rama 1-4 supera el límite térmico de 426 A.
Si las tensiones fueran, por el contrario, las más altas posibles (lo que de nuevo se produce
cuando el nudo 1 trabaja a 66 660 V) la línea 1-4 soportaría también una intensidad superior
a su intensidad nominal:
1 2 3 4 5 6
1 415,7369 372,1 442,438 0,0 0,0
2 -415,2217 0,0 -0,4 0,0 0,0
3 -371,6 0,0 0,0 160,8 80,8
4 -441,954 0,7 0,0 0,0 77,0
5 0,0 0,0 -160,7 0,0 25,6
6 0,0 0,0 -80,4 -76,8 -25,2
Una primera solución con la tensión del nudo 1 al nivel más bajo (66 276 V según el ejercicio
anterior) indica que la central en el nudo 5 debería tener una potencia de 6 836 kW:
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,002486 0,002066 0,004988 0,004786 0,004327 0,005028
TENSIÓN NUDO (V) 66 276 65 413 65 534 65 428 65 464 65 397
INT. FTES. DE NUDO (A) 1176,4 -417,9 -130,6 -366,2 -75,4 -183,5
POT. FTES. DE NUDO (W) 135 047 - 47 350 - 14 820 - 41 500 - 8 554 - 20 780
Y aunque se comprueba que la intensidad de la línea 1-4 ha sido corregida hasta 426 A, ocurre
que la tensión de la red podría ser algo más baja, pues en el nudo 6 ya no estamos a 65 340 V
sino a 65 397 V. Esto es debido a la presencia de la central del nudo 5.
Es preciso volver a ajustar la solución, en un segundo barrido, que primero reconsidere la
tensión en 1, y luego establezca la nueva potencia de la central en 5. En este caso, y con los
6 836 kW de generación ya en el nudo 5, la tensión en 1 puede bajar hasta 66 221 V.
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,002486 0,002066 0,004988 0,004786 0,004327 0,005028
TENSIÓN NUDO (V) 66 220 65 355 65 477 65 371 65 407 65 340
INT. FTES. DE NUDO (A) 1177,5 -418,3 -130,7 -366,5 -75,5 -183,6
POT. FTES. DE NUDO (W) 135 050 - 47 350 - 14 820 - 41 500 - 8 554 - 20 780
Para esta nueva tensión, la línea 1-4 presenta de nuevo una ligera sobrecarga, pues su inten-
sidad es de 426.4 A. Si la central en 5 proporcionara 6 969 kW esa sobrecarga se corregiría.
Pero ahora la solución de la red para esa nueva inyección de potencia vuelve a ser distinta,
aunque la diferencia es mucho menor que en el barrido anterior: se trata de que la central del
nudo 1 podría fijar su tensión en 66 220 V, apenas 1 V menos que antes:
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,002486 0,002066 0,004988 0,004786 0,004327 0,005028
TENSIÓN NUDO (V) 66 219 65 355 65 477 65 371 65 408 65 340
INT. FTES. DE NUDO (A) 1176,3 -418,3 -130,7 -366,5 -74,3 -183,6
POT. FTES. DE NUDO (W) 134 911 - 47 350 - 14 820 - 41 500 - 8 419 - 20 780
Y la central en el nudo 5 puede subir su potencia un par de kW más, quedando en 6 971 kW,
que era el resultado anunciado.
Por tanto, el problema queda resuelto en tres barridos de búsqueda de la solución: la central
a autorizar en el nudo 5 debería tener al menos la capacidad de incorporar 6 971 kW a la red
en su nudo, lo que aliviaría el problema de la línea 1-4 aunque se dieran las condiciones más
adversas posibles.
Pero aún hay que establecer la potencia máxima a autorizar para esa central. En este caso la
fija la condición de carga opuesta: la red trabaja al mínimo estadístico de carga y suministra
baja potencia a los nudos demandantes. Las tensiones tienden a ser las más altas y una nueva
central no hace sino elevarlas todavía más, incrementando el riesgo de sobretensión por en-
cima del +1.0% admisible.
La red se estudia ahora para su menor carga estadística (un 12% solo de la máxima carga), y
para las máximas tensiones posibles, que las fija el nudo 1 al alcanzar los 66 660 V permitidos.
En esas condiciones, la central en el nudo 5 solo puede inyectar 8 519 kW porque con ellos ya
9. Comprobar que construir una nueva línea de 200 km de longitud, del mismo tipo que el
resto de la red, para unir el nudo 3 con el 4, no soluciona el problema de los dos ejercicios
anteriores.
Comencemos por obtener la nueva matriz de admitancias de la red:
1 2 3 4 5 6
1 0,002486 -0,000816 -0,0008 -0,00087 0,0 0,0
2 -0,000816 0,002066 0,0 -0,00125 0,0 0,0
3 -0,0008 0,0 0,005988 -0,001 -0,003077 -0,001111
4 -0,00087 -0,00125 -0,001 0,005786 0,0 -0,002667
5 0,0 0,0 -0,003077 0,0 0,004327 -0,00125
6 0,0 0,0 -0,001111 -0,002667 -0,00125 0,005028
Y las soluciones:
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,002486 0,002066 0,005988 0,005786 0,004327 0,005028
TENSIÓN NUDO (V) 66 278 65 398 65 450 65 403 65 364 65 340
INT. FTES. DE NUDO (A) 1237,9 -418,0 -130,7 -366,3 -135,9 -183,6
POT. FTES. DE NUDO (W) 142 109 - 47 350 - 14 820 - 41 500 - 15 390 - 20 780
En efecto, la línea nueva 3-4 apenas soporta carga, aportando tan solo 3 072 kW al nudo 4
mediante una corriente de poco más de 27 A de intensidad, luego no alivia apenas el problema
de la 1-4, que sigue sobrecargada.
10. Calcular, para el caso de la red del ejercicio 6, gestionada de la manera simple que se
describió en él, la máxima potencia de generación que pudiera llegar a inyectarse, por sepa-
rado (caso A), en cada uno de los nudos de la red excluido el 1. Calcular, después, el valor
máximo de generación que puede admitir la red en dichos nudos inyectada en todos a la vez
(caso B).
La gestión simple de la red consiste en fijar, para cualquier caso de carga, el valor de la tensión
en el nudo 1 al máximo admisible, es decir, 66 660 V (un 1.0% superior al valor nominal de
66 kV).
En el procedimiento de cálculo que seguimos la potencia de generación inyectada en un nudo
se resta aritméticamente de la potencia suministrada a la carga en ese nudo, para así obtener
la potencia neta demandada por el nudo de la red. Hallando este valor en el caso más desfa-
vorable, quedan estudiadas todas las demás situaciones, porque en ellas la demanda no será
tan notable. Dicho caso más desfavorable es aquel en que la red suministra el mínimo de po-
tencia a las cargas, pero debe absorber toda la potencia nominal de generación, pues la ten-
sión en los nudos puede subir por encima del límite admisible. Como la red ya trabaja a la
tensión más alta posible, por efecto de la gestión simple adoptada, este límite se alcanza
pronto.
La solución de la red en la situación de partida sería esta (red cargada al 12%, tensiones lo
más altas posible, generación solo en el nudo 1):
NUDO 1 2 3 4 5 6
ADMITANCIA PROPIA 0,002486 0,002066 0,004988 0,004786 0,004327 0,005028
TENSIÓN NUDO (V) 66 660 66 560 66 569 66 564 66 559 66 556
INT. FTES. DE NUDO (A) 138,0 -49,3 -15,4 -32,8 -16,0 -21,6
POT. FTES. DE NUDO (W) 15 937 - 5 682 - 1 778 - 3 780 - 1 847 - 2 494
Estudiemos ahora el caso A). Se trata de ver, nudo a nudo, cuánta generación podemos admi-
tir en cada uno, comenzando por el nudo 2.
Inyectando 9 713 kW en el nudo 2 (rendimiento de la red 97,851%):
NUDO 1 2 3 4 5 6
TENSIÓN NUDO (V) 66 660 66 660 66 597 66 612 66 590 66 596
INT. FTES. DE NUDO (A) 53,8 34,9 -15,4 -32,8 -16,0 -21,6
POT. FTES. DE NUDO (W) 6 210 4 031 - 1 778 - 3 780 - 1 847 - 2 494
En caso de admitirse generación por todos los nudos (caso B) sin orden o criterio alguno, es
decir, de forma aleatoria, caben infinitas soluciones. Las restricciones serían que en ningún
caso se violase el límite de tensiones en todos los nudos de la red, y que el nudo 1 (una central
pura) se mantenga generando o como mucho se pare, pero no puede convertirse en un nudo
que absorba potencia de la red.
Un resultado de los infinitos posibles puede ser este:
NUDO 1 2 3 4 5 6
TENSIÓN NUDO (V) 66 660 66 623 66 660 66 639 66 659 66 648
INT. FTES. DE NUDO (A) 28,9 -28,1 9,0 -14,0 7,2 -0,2
POT. FTES. DE NUDO (W) 3 335 - 3 238 1 044 - 1 615 833 - 23
El límite de generación lo ha impuesto el nudo 3, que llega al máximo valor admisible para su
tensión, aunque no es el que más generación ha admitido. En otras soluciones puede ser otro
nudo distinto, y la cifra final deberá variar. En esta, por ejemplo, el límite lo ha impuesto el
nudo 4:
OTRA SOLUCIÓN CUALQUIERA (kW)
MÁXIMA GENERACIÓN ADMISIBLE POR TODOS LOS NUDOS A LA VEZ (CASO B)
NUDO 1 2 3 4 5 6
TENSIÓN NUDO (V) 66 660 66 621 66 656 66 660 66 658 66 660
INT. FTES. DE NUDO (A) 21,1 -45,8 -8,8 28,8 3,8 3,8
POT. FTES. DE NUDO (W) 2 440 - 5 282 - 1 020 3 327 436 436
El rendimiento en todas las soluciones a las que se pueda llegar será siempre algo superior al
caso anterior, donde sólo se permitía a un nudo, además del nudo 1, inyectar potencia. Y mu-
cho mejor al obtenido cuando toda la carga se suministraba desde el nudo 1.
Se comprueba otra vez, por tanto, la conclusión que se anunció en el apartado -, p. 132: la
generación distribuida mejora el rendimiento de la red. Sin embargo, creamos una duda:
¿cuál, de las infinitas soluciones posibles, será el reparto de generación óptimo para esta red?
¿Cuál es la solución para la que la generación total admisible es la más grande posible?
11. Sea la red de la figura cuya tensión nominal son 132 kV. Se indican las potencias su-
ministradas por cada nudo y las longitudes de las líneas en km. Considere que el conductor
es del tipo LA-280 y tome una impedancia en serie de 2x10-3 Ω/km y una admitancia paralelo
de 1x10-7 S/km.
Como se trata de un modelo en "pi", ha sido muy fácil adaptar la herramienta .XLSX para con-
siderar que una rama entre dos nudos es un transformador. Una vez adaptada, basta incluir
la relación de transformación t en la celda correspondiente de la columna B y simular la im-
pedancia serie de la máquina indicando una 'longitud' para esa rama.
En este caso, el transformador se intercalará entre los nudos 5 y 6, por lo que hay que prepa-
rar un poco el problema: primero estudiaremos qué ocurre cuando dichos nudos están sim-
plemente unidos por una línea de una longitud simbólica (un kilómetro, es decir, los nudos 5
y 6 son prácticamente el mismo nudo), para luego sustituirla por un transformador entre 5 y
6.
Así, para empezar hay que resolver la red en la que los nudos 5 y 6 están unidos por una línea
de 1 km, como muestra la imagen de la página anterior.
[Nota: como la red está escasamente mallada, el procedimiento de iteración utilizado re-
quiere forzar el cálculo repetidas veces hasta que la solución converja en los valores correc-
tos. Verifique que las potencias de las fuentes de nudo coinciden con los datos de partida del
problema para comprobar que el cálculo ha terminado, y de lo contrario, insista en las itera-
ciones hasta llegar a él].
Situación de partida
Toma: ! 6 7 8 9 10 11 12 13 14 15 16 !
99 99, 25 99, 50 99,75 100 100, 25 100, 50 100,75 101 101, 25 101, 50
t: ! !
100 100 100 100 100 100 100 100 100 100 100
Si la tensión del nudo 1 se mantiene en el máximo posible (132 660 V), la posición más ade-
cuada del cambiador de tomas del transformador es la [13] (obsérvese la relación t = 1,0075
en la fila de los nudos 5-6), pues se consigue que la tensión en el nudo 6 sea prácticamente la
nominal (132 014 V), como muestra el cálculo siguiente:
Si la tensión del nudo 1 fuera solo la nominal (132 000 V) la relación de transformación más
adecuada para el transformador entre 5 y 6 sería la [15] (t = 1,0125), que permite tener
132 006 V en el nudo 6 (véase la imagen en la página siguiente).
El problema es que con esta tensión aparece sobrecarga en la línea 1-2 (564,4 A > 563 A), por
lo que quizá resulte más conveniente la segunda opción que indica el enunciado.
6 PROBLEMAS PROPUESTOS
4.A. Complete la tabla adjunta teniendo en cuenta que la impedancia de las líneas es
únicamente resistiva y de valor z = 5.0 Ω/km, y que la admitancia es el valor inverso Y = 1/Z.
Represente gráficamente la red.
4.B. Suponga que en la red anterior todos los nudos son demandantes de potencia ex-
cepto el primero, que es de generación pura y que se encarga de mantener la tensión de la red
en el valor nominal de 45 kV. De las dos soluciones al problema de flujo de cargas que se
presentan más abajo, ¿cuál cree que es la que corresponde a este caso? Razone brevísima-
mente la respuesta.
Potencias suministradas:
Nudo 1 2 3 4 5 6
kW 29 400 35 280 14 700 21 560 43 120
NUDO 1 2 3 4 5 6
TENSIÓN DE NUDO (V) 45 000
INTENSIDAD DE FUENTES DE NUDO (A)
POTENCIA DE FUENTES DE NUDO (kW)
4.D. Tras el análisis de flujos de carga de una red de 6 nudos, se obtiene la siguiente
matriz de potencias entregadas desde cada nudo a las líneas. Represente gráficamente la red
y calcule las potencias de las fuentes de nudo cumplimentando la tabla adjunta.
Potencias entregadas a las líneas (kW):
1 2 3 4 5 6
1 28 239 46 245 48 086 23 076 0
2 - 28 042 0 0 - 1 358 0
3 - 45 928 0 10 648 0 0
4 - 47 630 0 - 10 620 0 43 550
5 - 22 918 1 358 0 0 0
6 0 0 0 - 43 120 0
4.E. A partir de los valores anteriores, calcule las pérdidas de potencia en la red y el ren-
dimiento de ésta en tantos por ciento.
4.F. Suponga que las intensidades que circulan desde los nudos a las líneas en la red que
estamos estudiando en los ejercicios anteriores resultan ser las de la matriz adjunta, en la que
se observa que la línea 1-4 está sobrecargada porque su intensidad límite son 600 A y está
soportando 616.9 A. A partir del gráfico que obtuvo en el ejercicio D, proponga dos amplia-
ciones distintas de la red que contribuyan a disminuir la carga soportada por dicha línea 1-4.
Represente ambas soluciones y razone brevísimamente por qué las propone.
Intensidades hacia las líneas desde cada nudo (A):
1 2 3 4 5 6
1 362,3 593,3 616,9 296,1 0,0
2 -362,3 0,0 0,0 -17,5 0,0
3 -593,3 0,0 137,6 0,0 0,0
4 -616,9 0,0 -137,6 0,0 564,1
5 -296,1 17,5 0,0 0,0 0,0
6 0,0 0,0 0,0 -564,1 0,0
4.G. En una red de 6 nudos como la que venimos estudiando en los ejercicios anteriores,
las potencias entregadas a las líneas bajo cierta hipótesis de carga resultan ser las de la matriz
a). Se considera la ampliación de la red y se barajan la opción 1 y la opción 2, cuyas matrices
de potencias se adjuntan en b) y en c). Represente gráficamente las tres redes, calcule sus
rendimientos y determine cuál es la opción de ampliación más beneficiosa.
1 2 3 4 5 6
1 28 239 37 662 56 847 23 076 0
2 - 28 042 0 0 - 1 358 0
3 - 37 452 0 22 752 0 0
4 - 56 208 0 - 22 623 0 43 551
5 - 22 918 1 358 0 0 0
6 0 0 0 - 43 120 0 0
1 2 3 4 5 6
1 38 888 29 861 41 265 35 474 0
2 - 38 514 0 0 9 114 0
3 - 29 729 0 15 029 0 0
4 - 40 929 0 - 14 972 0 20 621
5 - 35 101 - 9 106 0 0 22 647
6 0 0 0 - 20 525 - 22 595
1 2 3 4 5 6
1 28 239 47 190 46 949 23 076 0
2 - 28 042 0 0 - 1 358 0
3 - 46 860 0 9 180 0 22 980
4 - 46 514 0 - 9 159 0 20 393
5 - 22 918 1 358 0 0 0
6 0 0 - 22 821 - 20 299 0