Redes de tuberías

Leyes básicas
Las leyes físicas básicas que rigen el comportamiento del fluido en una red de tuberías son las
siguientes:

1. Ley de los nodos. La suma de los flujos entrantes en cualquier nodo debe ser 0. Esta ley
es consecuencia de aplicar la ecuación de continuidad en los nodos.
2. Ley de las mallas. La pérdida de carga total alrededor de cualquier malla cerrada debe
ser cero. Una vez definido el sentido de la malla se considerarán pérdidas de carga
positivas en las que coincida el sentido del flujo con el sentido de la malla y negativas en
caso contrario.
3. Ley de la pérdida de carga. Todas las pérdidas deben satisfacer el diagrama de Moody
(o sus ecuaciones equivalentes) y las correlaciones experimentales de las pérdidas
singulares.

Ley de nodos
Al aplicar la ecuación de continuidad en un nodo, para un fluido incompresible obtenemos que
el sumatorio de caudales entrantes en el nodo es 0. Establecemos como criterio que los caudales
entrantes en el nodo son positivos y los salientes son negativos.

∑ 𝑄𝑖 = 0

Q2 Q3

Q1 Q4

nodo

Ley de mallas
Sea ℎ𝑟 = ℎ𝑓 + ℎ𝑚 las pérdidas de carga en una tubería. Una vez definido un sentido para la
malla, se consideran pérdidas de carga positivas en aquellas tuberías donde el sentido del flujo
coincide con el sentido de la malla y negativas en el caso contrario. Entonces la suma de las
pérdidas de carga en cada una de las tuberías de la malla será 0.

∑ ℎ𝑟𝑖 = 0

t2
t1
malla

t3
t4

Con la ley de pérdidas de carga, podemos definir las pérdidas de carga en la tubería:
 8𝑓𝐿𝑄 2 8𝑄 2 8 𝑓𝐿
ℎ𝑟 = ℎ𝑓 + ℎ𝑚 = 2
+ ∑ 𝐾 · 2 4
= 2 4
( + ∑ 𝐾) 𝑄 2 = 𝑟 · 𝑄 2 , 𝑑𝑜𝑛𝑑𝑒 𝑟
𝑔𝜋 𝐷5 𝑔𝜋 𝐷 𝑔𝜋 𝐷 𝐷
8 𝑓𝐿
= ( + ∑ 𝐾)
𝑔𝜋 2 𝐷4 𝐷
Para poder darle sentido a la pérdida de carga, en lugar de utilizar el cuadrado del caudal,
utilizaremos el producto del módulo del caudal por el caudal.

ℎ𝑟 = 𝑟|𝑄|𝑄
Así la ley de mallas la podríamos enunciar como:

∑ 𝑟𝑖 |𝑄𝑖 |𝑄𝑖 = 0

Procedimiento para la solución de redes de tuberías

Para resolver los problemas de redas de tuberías seguiremos los siguientes pasos:

1. Trazado de la red.
2. Numerar mallas y nodos.
3. Establecer las ecuaciones de la red mediante la aplicación de la ley de nodos y la ley de
mallas.
4. Solución de sistema de ecuaciones no lineales mediante métodos iterativos, por
ejemplo: Hardy-Cross, Newton-Raphson…

Ejemplo
Calcular los caudales que circulan por la red de la figura, con los datos en ella reflejados.
Calcúlese después la presión en los nudos.

a) Resolver con un factor de fricción constante f=0,02.

b) Resolver con un factor de fricción dependiente del caudal. Sea la rugosidad de las
tuberías k = 0,046 mm.

L=800 m
85 l/s 56 m 48 m 100 l/s
D=250 mm

L=420 m
D=400 mm L=410 m
93 m D=300 mm

L=500 m L=810 m
55 m 57 m
D=600 mm D=400 mm
435 l/s
110 l/s
60 l/s
L=450 m L=920 m
D=400 mm 52 m D=300 mm

80 l/s
Una vez tenemos trazada la red, definimos los nodos y el sentido de los caudales.

85 l/s 2 100 l/s

2 3

1 1 3

6 435 l/s 1 4 4
7
60 l/s 110 l/s
2
6
5
5
80 l/s

Aplicamos la ley de nodos a cada uno de los nodos. En el nodo uno, tenemos un flujo entrante
de 435 l/s y 3 flujos salientes, Q1, Q4, Q6. Así la ecuación para el nodo 1 es:

−𝑄1 − 𝑄4 − 𝑄6 + 0,435 − 0,06 = −𝑄1 − 𝑄4 − 𝑄6 + 0,375 = 0 [𝐸𝑐. 𝑛𝑜𝑑𝑜 1]

Repetimos el proceso para el resto de los nodos:

𝑄1 − 𝑄2 − 0,085 = 0 [𝐸𝐶. 𝑛𝑜𝑑𝑜 2]

𝑄2 + 𝑄3 − 0,1 = 0 [𝐸𝑐. 𝑛𝑜𝑑𝑜 3]
−𝑄3 + 𝑄4 + 𝑄5 − 0,11 = 0 [𝐸𝑐. 𝑛𝑜𝑑𝑜 4]
−𝑄5 + 𝑄6 − 0,08 = 0 [𝐸𝑐. 𝑛𝑜𝑑𝑜 5]
De las 5 ecuaciones anteriores solo 4 de ellas son linealmente independientes, porque
necesitamos las 2 ecuaciones de las leyes de mallas. Para la malla 1 tenemos que el flujo en las
tuberías 1 y 2 tienen el mismo sentido que le hemos asignado a la malla y en las tuberías 3 y 4,
el sentido contrario por lo que se cumple:

ℎ𝑟 1 + ℎ𝑟2 − ℎ𝑟3 − ℎ𝑟4 = 0

Tal y como vimos anteriormente, matemáticamente podemos expresar las pérdidas de carga de
la siguiente manera:

𝑟1 |𝑄11 |𝑄11 + 𝑟2 |𝑄12 |𝑄12 + 𝑟3 |𝑄13 |𝑄13 + 𝑟4 |𝑄14 |𝑄14 = 0 [𝐸𝑐. 𝑚𝑎𝑙𝑙𝑎 1]
Donde:

𝑄11 = 𝑄1 , 𝑄12 = 𝑄2 , 𝑄13 = −𝑄3 , 𝑄14 = −𝑄4

Repitiendo el proceso para la malla 2 tenemos:
 𝑟4 |𝑄24 |𝑄24 + 𝑟5 |𝑄25 |𝑄25 + 𝑟6 |𝑄26 |𝑄26 = 0 [𝐸𝑐. 𝑚𝑎𝑙𝑙𝑎 2]
Donde:

𝑄24 = 𝑄4 , 𝑄25 = −𝑄5 , 𝑄26 = −𝑄6

Para resolver el problema y obtener el valor de los caudales que fluyen por cada una de las
líneas, debemos resolver el sistema de ecuaciones.

Método de Hardy-Cross
En el método de Hardy-Cross se parte de unos valores iniciales de caudales compatibles con las
ecuaciones de los nodos. Para establecer estos valores de caudal iniciales, con ayuda del trazado
de la red, distribuimos los caudales de una manera razonable.

85 l/s 2 100 l/s

2 3

35 l/s
120 l/s

65 l/s
1 1 3

6 435 l/s 1 4 4
7
60 l/s 120 l/s 110 l/s
2
135 l/s

6
5
5
80 l/s

Al nodo 1 llegan 435 l/s y tenemos un consumo en el nodo de 60 l/s, para cumplir con la ecuación
del nodo 1 tendremos que repartir 375 l/s entre las líneas 1, 4 y 6, se puede realizar de manera
más o menos uniforme. En este caso, al definir los caudales 1, 4 y 6, las ecuaciones de los nodos
fijan el valor del resto de caudales. Por ejemplo, si estudiamos el nodo 2, al que llega un caudal
de 120 l/s y existe un consumo de 85 l/s, la ecuación del nodo 2 fija el caudal por la tubería 2 en
35 l/s. Evaluando el nodo 3, tendríamos que el caudal por la tubería 3 debe ser 65 l/s. Evaluando
el nodo 4, el caudal en la tubería 5 debe ser de 55 l/s. Y podemos comprobar que con estos
valores de caudal se cumple la ecuación en el nodo 5. Es decir, se cumplen todas las ecuaciones
de los nodos, esta es la condición necesaria para resolver el problema utilizando el método de
Hardy-Cross.

Estos caudales no cumplen con las ecuaciones de las mallas:

L [m] d [m] Q [m3/s] r

1 420 0.4000 0.1200 67.7798
2 800 0.2500 0.0350 1.3538e+03
3 410 0.3000 0.0650 278.8231
 4 810 0.4000 0.1200 130.7183
5 920 0.3000 0.0550 625.6519
6 450 0.4000 0.1350 72.6213

𝑟1 |𝑄11 |𝑄11 + 𝑟2 |𝑄12 |𝑄12 + 𝑟3 |𝑄13 |𝑄13 + 𝑟4 |𝑄14 |𝑄14 = −0,4260 𝑚

𝑚3 𝑚3 𝑚3 𝑚3
𝑐𝑜𝑛: 𝑄11 = 0,120 , 𝑄12 = 0,035 , 𝑄13 = −0,065 , 𝑄14 = −0,120
𝑠 𝑠 𝑠 𝑠
𝑟4 |𝑄24 |𝑄24 + 𝑟5 |𝑄25 |𝑄25 + 𝑟6 |𝑄26 |𝑄26 = −13,38 𝑚
𝑚3 𝑚3 𝑚3
𝑐𝑜𝑛: 𝑄24 = 0,120 , 𝑄25 = −0,055 , 𝑄26 = −0,135
𝑠 𝑠 𝑠
En ambos casos las pérdidas de carga en los tramos en sentido contrario de la malla son
superiores a las pérdidas de carga en los tramos en el sentido de la malla. Para conseguir nuestro
objetivo debemos aumentar los caudales en los tramos donde el caudal circula en el sentido de
la malla y disminuir los otros.

¿Cuánto debemos variar los caudales para llegar a la solución?

Las ecuaciones de las mallas son unas funciones multivariente de la forma:

𝑓(𝑄1 , 𝑄2 , … , 𝑄𝑛 ) = ∑ 𝑟𝑖 |𝑄𝑖 |𝑄𝑖

Buscamos un valor de Δ𝑄 tal que:

𝑓(𝑄1 + Δ𝑄, 𝑄2 + Δ𝑄, … , 𝑄𝑛 + Δ𝑄) = 0

Aplicando el desarrollo de Taylor para una función multivariante tenemos que:
𝜕𝑓
𝑓(𝑄1 + Δ𝑄, 𝑄2 + Δ𝑄, … , 𝑄𝑛 + Δ𝑄) = 𝑓(𝑄1 , 𝑄2 , … , 𝑄𝑛 ) + ∑ Δ𝑄 + 𝑅2 = 0
𝜕𝑄𝑖
Si Δ𝑄 es lo suficientemente pequeño 𝑅2 ≈ 0. Entonces:
𝜕𝑓 𝜕𝑓
𝑓(𝑄1 , 𝑄2 , … , 𝑄𝑛 ) + ∑ Δ𝑄 = 0 ⇒ 𝑓(𝑄1 , 𝑄2 , … , 𝑄𝑛 ) + Δ𝑄 ∑ =0⇒
𝜕𝑄𝑖 𝜕𝑄𝑖
−𝑓(𝑄1 , 𝑄2 , … , 𝑄𝑛 )
Δ𝑄 =
𝜕𝑓
∑
𝜕𝑄𝑖
𝜕𝑓 𝜕(∑ 𝑟𝑖 |𝑄𝑖 |𝑄𝑖 )
= = 2𝑟𝑖 |𝑄𝑖 |
𝜕𝑄𝑖 𝜕𝑄𝑖
∑ 𝑟𝑖 |𝑄𝑖 |𝑄𝑖
Δ𝑄 = −
∑ 2𝑟𝑖 |𝑄𝑖 |

Volviendo a nuestro ejemplo:

𝑟1 |𝑄11 |𝑄11 + 𝑟2 |𝑄12 |𝑄12 + 𝑟3 |𝑄13 |𝑄13 + 𝑟4 |𝑄14 |𝑄14 = 0

∑ 𝑟𝑖 |𝑄𝑖 |𝑄𝑖 𝑟1 |𝑄11 |𝑄11 + 𝑟2 |𝑄12 |𝑄12 + 𝑟3 |𝑄13 |𝑄13 + 𝑟4 |𝑄14 |𝑄14
Δ𝑄1 = − =−
∑ 2𝑟𝑖 |𝑄𝑖 | 2𝑟1 |𝑄11 | + 2𝑟2 |𝑄12 | + 2𝑟3 |𝑄13 | + 2𝑟4 |𝑄14 |
 𝑚3 𝑚3 𝑚3 𝑚3
𝑐𝑜𝑛: 𝑄11 = 0,120 , 𝑄12 = 0,035 , 𝑄13 = −0,065 , 𝑄14 = −0,120
𝑠 𝑠 𝑠 𝑠
𝑟4 |𝑄24 |𝑄24 + 𝑟5 |𝑄25 |𝑄25 + 𝑟6 |𝑄26 |𝑄26 = 0
∑ 𝑟𝑖 |𝑄𝑖 |𝑄𝑖 𝑟4 |𝑄24 |𝑄24 + 𝑟5 |𝑄25 |𝑄25 + 𝑟6 |𝑄26 |𝑄26
Δ𝑄2 = − =−
∑ 2𝑟𝑖 |𝑄𝑖 | 2𝑟4 |𝑄24 | + 2𝑟5 |𝑄25 | + 2𝑟6 |𝑄26 |
𝑚3 𝑚3 𝑚3
𝑐𝑜𝑛: 𝑄24 = 0,120 , 𝑄25 = −0,055 , 𝑄26 = −0,135
𝑠 𝑠 𝑠

Con los valores del caudal obtenemos:

Δ𝑄1 = 0,0024, Δ𝑄2 = 0,0111

Con estas variaciones de caudales calcularemos los proporcionados por la primera iteración:

𝑚3
𝑄1 𝑘 = 𝑄11𝑘 = 𝑄11𝑘−1 + Δ𝑄1 = 0,120 + 0,0024 = 0,1224
𝑠
𝑚3
𝑄2𝑘 = 𝑄12𝑘 = 𝑄12𝑘−1 + Δ𝑄1 = 0,035 + 0,0024 = 0,0374
𝑠
𝑚3
−𝑄3𝑘 = 𝑄13𝑘 = 𝑄13𝑘−1 + Δ𝑄1 = −0,065 + 0,0024 = −0,0626
𝑠
La línea 4 se ve afectada por Δ𝑄1 y Δ𝑄2 . Aplicaremos primero Δ𝑄1 y sobre el caudal obtenido
aplicaremos y Δ𝑄2 .

−𝑄4′ 𝑘 = 𝑄14𝑘 = 𝑄14𝑘−1 + Δ𝑄1 = −𝑄4𝑘−1 + Δ𝑄1 ⇒ 𝑄4′ 𝑘 = 𝑄4𝑘−1 − Δ𝑄1

′
𝑚3
𝑄4𝑘 = 𝑄24𝑘 = 𝑄24 + Δ𝑄2 = 𝑄4𝑘−1 − Δ𝑄1 + Δ𝑄2 = 0,120 − 0,0024 + 0,0111 = 0,1287
𝑘
𝑠
𝑚3
−𝑄5𝑘 = 𝑄25𝑘 = 𝑄25𝑘−1 + Δ𝑄2 = −0,0439
𝑠
𝑚3
−𝑄6𝑘 = 𝑄26𝑘 = 𝑄26𝑘−1 + Δ𝑄2 = −0,1239
𝑠
Los caudales en m3/s resultantes de es esta primera iteración son:

Q1 Q2 Q3 Q4 Q5 Q6
0.1224 0.0374 0.0626 0.1287 0.0439 0.1239

𝑟1 |𝑄1 |𝑄1 + 𝑟2 |𝑄2 |𝑄2 + 𝑟3 |𝑄3 |𝑄3 + 𝑟4 |𝑄4 |𝑄4 = −0,3528 𝑚

𝑚3 𝑚3 𝑚3 𝑚3
𝑐𝑜𝑛: 𝑄11 = 0,1224 , 𝑄12 = 0,0374 , 𝑄13 = −0,0626 , 𝑄14 = −0,1287
𝑠 𝑠 𝑠 𝑠
𝑟4 |𝑄4 |𝑄4 + 𝑟5 |𝑄5 |𝑄5 + 𝑟6 |𝑄6 |𝑄6 = −0,1514 𝑚
𝑚3 𝑚3 𝑚3
𝑐𝑜𝑛: 𝑄24 = 0,1287 , 𝑄25 = −0,0439 , 𝑄26 = −0,1239
𝑠 𝑠 𝑠
Estos caudales todavía no cumplen con las ecuaciones de mallas, repetiremos el proceso hasta
que los valores de Δ𝑄 sean lo suficientemente pequeños. En la siguiente tabla podemos
observar la evolución del cálculo iterativo.

Q1 Q2 Q3 Q4 Q5 Q6 ΔQ1 ΔQ2
It. hr1 [m] hr2 [m]
[m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s]

0 0.1200 0.0350 0.0650 0.1200 0.0550 0.1350 -0.4260 -13.338 0.0024 0.0111

1 0.1224 0.0374 0.0626 0.1287 0.0439 0.1239 -0.3528 -0.1514 0.0019 0.0014

2 0.1243 0.0393 0.0607 0.1283 0.0424 0.1224 -0.0438 -0.0651 0.0002 0.0006

3 0.1245 0.0395 0.0605 0.1287 0.0418 0.1218 -0.0209 -0.0079 0.0001 0.0001

4 0.1246 0.0396 0.0604 0.1286 0.0417 0.1217 -0.0026 -0.0037 0.0000 0.0000

85 l/s 2 100 l/s

2 3

39,6 l/s
124,6 l/s

60,4 l/s
1 1 3

6 435 l/s 1 4 4
7
60 l/s 128,6 l/s 110 l/s
2
121,7 l/s

6
5
5
80 l/s

Una vez tenemos los caudales en cada una de las líneas, podremos establecer la presión en los
nodos, conocida la presión en uno de ellos.

Aplicamos la ecuación de Bernoulli generalizado entre el depósito y el nodo 1.

𝑝0 𝑉02 𝑝1 𝑉12 𝑝𝑚 1 𝑉12

+ + 𝑧0 = + + 𝑧1 + ℎ𝑟01 ⇒ 𝑧0 = + + 𝑧1 + ℎ𝑟01 ≈ 𝐻1 + ℎ𝑟01
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔 𝜌𝑔 2𝑔
Donde 𝐻1 es la altura piezométrica en el nodo. Si despreciamos las pérdidas de carga
secundarias, tenemos:
2
𝑚3
8𝑓7 𝐿7 𝑄72 8 · 0,02 · 500 𝑚 · (0,435 𝑠 )
ℎ𝑟01 = = 𝑚 = 2,01 𝑚 ⇒
𝑔𝜋 2 𝐷75 9,81 2 · 𝜋 2 · (0,6 𝑚)5
𝑠
𝑝𝑚 1
𝐻1 = 𝑧0 − ℎ𝑟01 = 93 𝑚 − 2,01 𝑚 = 90,99 𝑚 = + 𝑧1 ⇒
𝜌𝑔
 𝑝𝑚 1
= 𝐻1 − 𝑧1 = 90,99 𝑚 − 55 𝑚 = 35,99 𝑚 ⇒ 𝑝𝑚1 = 353,05 𝑘𝑃𝑎
𝜌𝑔
Desde el nodo 1, aplicando la ecuación de Bernoulli generalizado podemos ir obteniendo la
altura piezométrica de cada uno de los nodos y su presión manométrica. Por ejemplo, para el
nodo 2:
𝑉12 𝑉22
≈ ≈0
𝑝1 𝑉12 𝑝2 𝑉22 𝑝1 2𝑔 2𝑔𝑝2 𝑝𝑚 1
+ + 𝑧1 = + + 𝑧2 + ℎ𝑟12 ⇒ + 𝑧1 = + 𝑧2 + ℎ𝑟12 ⇒ + 𝑧1
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔 𝜌𝑔 𝜌𝑔 𝜌𝑔
𝑝𝑚 2
= + 𝑧2 + ℎ𝑟12 ⇒ 𝐻1 = 𝐻2 + ℎ𝑟12
𝜌𝑔
8𝑓1 𝐿1 𝑄12
ℎ𝑟12 = = 1,05 𝑚
𝑔𝜋 2 𝐷15
𝑝𝑚 2
𝐻2 = 𝐻1 − ℎ𝑟12 = 90,99 𝑚 − 1,05 𝑚 = 89,94 𝑚 ⇒ = 𝐻2 − 𝑧2 = 89,94 𝑚 − 56 𝑚
𝜌𝑔
= 33,94 𝑚 ⇒ 𝑝𝑚2 = 332,92 𝑘𝑃𝑎

Siguiendo el procedimiento para el resto de los nodos obtenemos los siguientes resultados:

𝒑𝒎𝒋 /𝝆𝒈
Nodo Hj [m] Zj [m] 𝒑𝒎𝒋 [kPa]
[m]
1 90,99 55 35,99 353,06
2 89,94 56 33,94 332,92
3 87,81 48 39,81 390,54
4 88,83 57 31,83 312,21
5 89,91 52 37,91 371,93

Si en lugar de considerar el factor de fricción constante, se calculara para las condiciones de flujo
en cada tubería, el procedimiento sería igual al desarrollado con el factor de fricción constante,
con la salvedad de que para cada iteración deberemos calcular el factor de fricción
correspondiente al caudal que fluye por cada tubería y las características de ésta.

En nuestro ejemplo, con los caudales definidos para el inicio del cálculo iterativo tenemos los
siguientes factores de fricción y valores de las contantes “r”.

L [m] d [m] Q [m3/s] Re f r

1 420 0,4 0,12 3,82E+05 0,0151 51,16
2 800 0,25 0,035 1,78E+05 0,0173 1169,91
3 410 0,3 0,065 2,76E+05 0,0161 224,18
4 810 0,4 0,12 3,82E+05 0,0151 98,67
5 920 0,3 0,055 2,33E+05 0,0164 514,13
6 450 0,4 0,135 4,30E+05 0,0149 54,05

𝑟1 |𝑄11 |𝑄11 + 𝑟2 |𝑄12 |𝑄12 + 𝑟3 |𝑄13 |𝑄13 + 𝑟4 |𝑄14 |𝑄14 = −0,1982 𝑚

𝑚3 𝑚3 𝑚3 𝑚3
𝑐𝑜𝑛: 𝑄11 = 0,120 , 𝑄12 = 0,035 , 𝑄13 = −0,065 , 𝑄14 = −0,120
𝑠 𝑠 𝑠 𝑠
 𝑟4 |𝑄24 |𝑄24 + 𝑟5 |𝑄25 |𝑄25 + 𝑟6 |𝑄26 |𝑄26 = −11,95 𝑚

𝑚3 𝑚3 𝑚3
𝑐𝑜𝑛: 𝑄24 = 0,120 , 𝑄25 = −0,055 , 𝑄26 = −0,135
𝑠 𝑠 𝑠
∑ 𝑟𝑖 |𝑄𝑖 |𝑄𝑖 𝑟1 |𝑄11 |𝑄11 + 𝑟2 |𝑄12 |𝑄12 + 𝑟3 |𝑄13 |𝑄13 + 𝑟4 |𝑄14 |𝑄14
Δ𝑄1 = − =− = 0,001348
∑ 2𝑟𝑖 |𝑄𝑖 | 2𝑟1 |𝑄11 | + 2𝑟2 |𝑄12 | + 2𝑟3 |𝑄13 | + 2𝑟4 |𝑄14 |
∑ 𝑟𝑖 |𝑄𝑖 |𝑄𝑖 𝑟4 |𝑄24 |𝑄24 + 𝑟5 |𝑄25 |𝑄25 + 𝑟6 |𝑄26 |𝑄26
Δ𝑄2 = − =− = 0,001875
∑ 2𝑟𝑖 |𝑄𝑖 | 2𝑟4 |𝑄24 | + 2𝑟5 |𝑄25 | + 2𝑟6 |𝑄26 |
Con los Δ𝑄 calculados, se corrigen los caudales y se vuelve a calcular los factores de fricción y el
valor de las constantes “r”.

L [m] d [m] Q [m3/s] Re f r

1 420 0,4 0,1213 3,86E+05 0,0151 51,09
2 800 0,25 0,0363 1,85E+05 0,0172 1163,74
3 410 0,3 0,0637 2,70E+05 0,0161 224,78
4 810 0,4 0,1305 4,15E+05 0,0149 97,69
5 920 0,3 0,0432 1,83E+05 0,0170 531,88
6 450 0,4 0,1232 3,92E+05 0,0150 54,64

Realizando el proceso de cálculo iterativo, obtendremos la siguiente evolución:

Q1 Q2 Q3 Q4 Q5 Q6 ΔQ1 ΔQ2
It. hr1 [m] hr2 [m]
[m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s]

0 0,1200 0,0350 0,0650 0,1200 0,0550 0,1350 -0,1982 -1,1195 0,001348 0,011805

1 0,1213 0,0363 0,0637 0,1305 0,0432 0,1232 -0,2833 -0,1592 0,001875 0,001875

2 0,1232 0,0382 0,0618 0,1305 0,0413 0,1213 -0,0608 -0,0573 0,000394 0,000690

3 0,1236 0,0386 0,0614 0,1308 0,0406 0,1206 -0,0204 -0,0135 0,000132 0,000164

4 0,1237 0,0387 0,0613 0,1308 0,0405 0,1205 -0,0053 -0,0041 0,000034 0,000050

5 0,1238 0,0388 0,0612 0,1308 0,0404 0,1204 -0,0013 -0,0009 0,000000 0,000000

Para los caudales correspondientes a la solución los valores del número de Reynolds y factores
de fricción son:

L [m] d [m] Q [m3/s] Re f r

1 420 0,4 0,1238 3,94E+05 0,0150 50,97
2 800 0,25 0,0388 1,98E+05 0,0170 1153,43
3 410 0,3 0,0612 2,60E+05 0,0162 225,91
4 810 0,4 0,1308 4,16E+05 0,0149 97,66
5 920 0,3 0,0404 1,72E+05 0,0172 537,14
6 450 0,4 0,1204 3,83E+05 0,0151 54,80
Método de Newton-Raphson
Para resolver el problema anterior debemos resolver el siguiente sistema de ecuaciones no
lineales:

𝑄1 − 𝑄2 − 0,085 = 0 [𝐸𝐶. 𝑛𝑜𝑑𝑜 2]

𝑄2 + 𝑄3 − 0,1 = 0 [𝐸𝑐. 𝑛𝑜𝑑𝑜 3]
−𝑄3 + 𝑄4 + 𝑄5 − 0,11 = 0 [𝐸𝑐. 𝑛𝑜𝑑𝑜 4]
−𝑄5 + 𝑄6 − 0,08 = 0 [𝐸𝑐. 𝑛𝑜𝑑𝑜 5]
𝑟1 |𝑄1 |𝑄1 + 𝑟2 |𝑄2 |𝑄2 − 𝑟3 |𝑄3 |𝑄3 − 𝑟4 |𝑄4 |𝑄4 = 0 [𝐸𝑐. 𝑚𝑎𝑙𝑙𝑎 1]
𝑟4 |𝑄4 |𝑄4 − 𝑟5 |𝑄5 |𝑄5 − 𝑟6 |𝑄6 |𝑄6 = 0 [𝐸𝑐. 𝑚𝑎𝑙𝑙𝑎 2]

85 l/s 2 100 l/s

2 3

1 1 3

6 435 l/s 1 4 4
7
60 l/s 110 l/s
2
6
5
5
80 l/s

Sea 𝑄⃗ 𝑖 = (𝑄1𝑖 , 𝑄2𝑖 , … , 𝑄𝑛𝑖 ) un conjunto de valores que no satisfacen el sistema de ecuaciones.
El objetivo del método es encontrar un vector 𝑄 ⃗ 𝑖+1 que satisfaga el sistema de ecuaciones.

⃗ 𝑖+1 = 𝑄
Sea 𝑄 ⃗ 𝑖 + Δ𝑄
⃗ 𝑖 , 𝑑𝑜𝑛𝑑𝑒 Δ𝑄
⃗ 𝑖 = (Δ𝑄1𝑖 , Δ𝑄2𝑖 , … , Δ𝑄𝑛𝑖 ), entonces 𝑄
⃗ 𝑖+1 debe cumplir que:
𝑛
⃗ 𝑖)
𝜕𝑓𝑗 (𝑄
⃗ 𝑖+1 ) = 0 ⇒ 𝑓𝑗 (𝑄
𝑓𝑗 (𝑄 ⃗ 𝑖+1 ) = 𝑓𝑗 (𝑄
⃗ 𝑖 + Δ𝑄
⃗ 𝑖 ) ≈ 𝑓𝑗 (𝑄
⃗ 𝑖 ) + ∇𝑓𝑗 · Δ𝑄
⃗ 𝑖 = 𝑓𝑗 (𝑄
⃗ 𝑖) + ∑ · Δ𝑄𝑘𝑖
𝜕𝑄𝑘
𝑘=1

Si aplicamos la expresión para la función j sobre las n funciones obtendríamos la siguiente

ecuación matricial:
 ⃗ 𝑖)
𝜕𝑓1 (𝑄 𝜕𝑓1 (𝑄⃗ 𝑖) ⃗ 𝑖)
𝜕𝑓1 (𝑄
⋯
⃗ 𝑖) 𝜕𝑄1 𝜕𝑄2 𝜕𝑄𝑛
0 𝑓1 (𝑄 Δ𝑄1𝑖
⃗ 𝑖)
𝜕𝑓2 (𝑄 ⃗
𝜕𝑓2 (𝑄𝑖 ) ⃗ 𝑖)
𝜕𝑓2 (𝑄
⃗ Δ𝑄
(0) = 𝑓2 (𝑄𝑖 ) + 𝜕𝑄1 𝜕𝑄2
…
𝜕𝑄𝑛 ( 2𝑖 ) ⇒ 0 = 𝐵 + 𝐽 · Δ𝑄
⃗ 𝑖 ⇒ −𝐵
⋮ ⋮ ⋮
0 ⃗ ⋮ ⋮ ⋱ ⋮ Δ𝑄𝑛𝑖
(𝑓𝑛 (𝑄𝑖 )) ⃗ ⃗
𝜕𝑓𝑛 (𝑄𝑖 ) 𝜕𝑓𝑛 (𝑄𝑖 ) ⃗ 𝑖)
𝜕𝑓𝑛 (𝑄
…
( 𝜕𝑄1 𝜕𝑄2 𝜕𝑄𝑛 )
⃗ −1 −1 ⃗
= 𝐽 · Δ𝑄𝑖 ⇒ −J · 𝐵 = J · 𝐽 · Δ𝑄𝑖 ⇒ Δ𝑄 ⃗ 𝑖 = −J−1 · 𝐵 ⇒
−1
⃗ 𝑖)
𝜕𝑓1 (𝑄 ⃗ 𝑖)
𝜕𝑓1 (𝑄 ⃗ 𝑖)
𝜕𝑓1 (𝑄
⋯
𝜕𝑄1 𝜕𝑄2 𝜕𝑄𝑛 ⃗ 𝑖)
Δ𝑄1𝑖 𝑓1 (𝑄
⃗ 𝑖)
𝜕𝑓2 (𝑄 ⃗ 𝑖)
𝜕𝑓2 (𝑄 ⃗ 𝑖)
𝜕𝑓2 (𝑄
Δ𝑄 𝑓2 (𝑄⃗ 𝑖)
( 2𝑖 ) = − 𝜕𝑄1 𝜕𝑄2
…
𝜕𝑄𝑛 ⇒
⋮ ⋮
Δ𝑄𝑛𝑖 ⋮ ⋮ ⋱ ⋮ ⃗
𝜕𝑓𝑛 (𝑄⃗ 𝑖 ) 𝜕𝑓𝑛 (𝑄⃗ 𝑖) ⃗ 𝑖)
𝜕𝑓𝑛 (𝑄 (𝑓𝑛 (𝑄𝑖 ))
…
( 𝜕𝑄1 𝜕𝑄2 𝜕𝑄𝑛 )
𝑄1 𝑖+1 𝑄1 𝑖 Δ𝑄1𝑖
𝑄 𝑄 Δ𝑄
( 2 𝑖+1 ) = ( 2 𝑖 ) + ( 2𝑖 )
⋮ ⋮ ⋮
𝑄𝑛 𝑖+1 𝑄𝑛 𝑖 Δ𝑄𝑛𝑖
−1
𝜕𝑓1 (𝑥𝑖 ) 𝜕𝑓1 (𝑥𝑖 ) 𝜕𝑓1 (𝑥𝑖 )
⋯
𝜕𝑥1 𝜕𝑥2 𝜕𝑥𝑛 𝑓1 (𝑄⃗ 𝑖)
𝑄1 𝑖
𝜕𝑓2 (𝑥𝑖 ) 𝜕𝑓2 (𝑥𝑖 ) 𝜕𝑓2 (𝑥𝑖 )
𝑄 … 𝑓2 (𝑄⃗ 𝑖)
= ( 2 𝑖) − 𝜕𝑥1 𝜕𝑥2 𝜕𝑥𝑛
⋮ ⋮
𝑄𝑛 𝑖 ⋮ ⋮ ⋱ ⋮ ⃗
𝜕𝑓𝑛 (𝑥𝑖 ) 𝜕𝑓𝑛 (𝑥𝑖 ) 𝜕𝑓𝑛 (𝑥𝑖 ) (𝑓𝑛 (𝑄𝑖 ))
…
( 𝜕𝑥1 𝜕𝑥2 𝜕𝑥𝑛 )
Como criterio de convergencia podemos utilizar el siguiente estimador de error:

‖Δ𝑄 ⃗ 𝑖‖ √∑𝑛𝑗=1 Δ𝑄𝑗𝑖

2

𝜀= =
⃗ 𝑖+1 ‖
‖Q
√∑𝑛𝑗=1 𝑄𝑗𝑖+1
2

Calculamos el Jacobiano del sistema de ecuaciones

𝜕𝑓1 𝜕𝑓1
⋯
𝜕𝑄1 𝜕𝑄𝑛
𝐽= ⋮ ⋱ ⋮
𝜕𝑓𝑛 𝜕𝑓𝑛
⋯
(𝜕𝑄𝑛 𝜕𝑄𝑛 )
−1 0 0 −1 0 −1
1 −1 0 0 0 0
0 1 1 0 0 0
𝐽= 1
0 0 −1 0 0
2𝑟1 |𝑄1 | 2𝑟2 |𝑄2 | −2𝑟3 |𝑄3 | −2𝑟4 |𝑄4 | 0 0
( 0 0 0 2𝑟4 |𝑄4 | −2𝑟5 |𝑄5 | −2𝑟6 |𝑄6 |)
 8 𝑓𝑖𝑘 𝐿𝑖
𝐷𝑜𝑛𝑑𝑒 𝑟𝑖𝑘 = 4 ( + ∑ 𝐾𝑖𝑗 )
𝑔𝜋 2 𝐷𝑖 𝐷𝑖
𝑗

Solución suponiendo un factor de fricción contante para todas las tuberías 𝑓𝑖𝑘 = 𝑓 = 0,02, los
valores de ri no dependen de la iteración, y tendrán los siguientes valores:

1 2 3 4 5 6
𝟐
𝒔
𝒓 [ 𝟓] 67,78 1353,76 278,82 130,72 625,65 72,62
𝒎

Fijamos unos valores iniciales para la iteración, en este caso no es necesario que cumplan con la
ecuación de los nodos, por ejemplo, definimos 𝑄 ⃗ 0 = (0.1, 0.1, 0.1, 0.1, 0.1, 0.1). Así el Jacobiano
del sistema es el siguiente:
−1 0 0 −1 0 −1
1 −1 0 0 0 0
0 1 1 0 0 0
𝐽=
0 0 −1 1 0 0
13,5560 270,7517 −55,7646 −26,1437 0 0
( 0 0 0 26,1437 −125,1304 −14,5243)

1 2 3 4 5 6
3
𝑚
⃗ [ ]
Δ𝑄 0,03597 -0,04903 -0,05097 0,00673 -0,04770 0,03230
𝑠
𝑚3
⃗1 [ ]
𝑄 0,13597 0,05097 0,04903 0,10673 0,05230 0,13230
𝑠

⃗ ‖ 0,0983
‖Δ𝑄
𝜀= = = 41,86%
‖Q1 ‖ 0,2348
Siguiendo con la iteración obtenemos los siguientes resultados:

It. 𝑸𝟏 𝑸𝟐 𝑸𝟑 𝑸𝟒 𝑸𝟓 𝑸𝟔 ε
0 0,1000 0,1000 0,1000 0,1000 0,1000 0,1000
1 0,1360 0,0510 0,0490 0,1067 0,0523 0,1323 41,86%
2 0,1250 0,0400 0,0600 0,1282 0,0417 0,1217 13,93%
3 0,1246 0,0396 0,0604 0,1287 0,0417 0,1217 0,35%
4 0,1246 0,0396 0,0604 0,1287 0,0417 0,1217 0,00%

Si en lugar de tomar un factor de fricción constante, utilizamos el factor de fricción

correspondiente al número de Reynolds en cada tubería para cada iteración, tendríamos que
recalcular el jacobiano en cada una de las iteraciones.

1 2 3 4 5 6
3
𝑚
⃗𝟎 [ ]
𝑄 0,1 0,1 0,1 0,1 0,1 0,1
𝑠
𝒇 0,0154 0,0153 0,0153 0,0154 0,0153 0,0154
 𝒔𝟐
𝒓𝟎 [ 𝟓 ] 52,36 1037,44 213,27 100,98 478,56 56,10
𝒎

Así, el Jacobiano para la primera iteración será:

−1 0 0 −1 0 −1
1 −1 0 0 0 0
0 1 1 0 0 0
𝐽=
0 0 −1 1 0 0
10,4716 207,48 −42,6540 −20,1953 0 0
( 0 0 0 20,1953 −95,7114 −11,2196)

1 2 3 4 5 6
3
𝑚
⃗ [ ]
Δ𝑄 0,03598 -0,04902 -0,05098 0,00671 -0,04768 0,03232
𝑠
𝑚3
⃗1 [ ]
𝑄 0,13598 0,05098 0,04902 0,10671 0,05232 0,13232
𝑠

⃗ ‖ 0,0983
‖Δ𝑄
𝜀= = = 41,86%
‖Q1 ‖ 0,2348

En la siguiente iteración ajustamos los factores de fricción a los nuevos caudales y recalculamos
los valores r y el jacobiano.

1 2 3 4 5 6
3
𝑚
⃗1 [ ]
𝑄 0,13598 0,05098 0,04902 0,10671 0,05232 0,13232
𝑠
𝒇 0,0149 0,0165 0,0167 0,0153 0,0165 0,0149
𝒔𝟐
𝒓𝟏 [ 𝟓 ] 50,4063 1113,5294 232,7747 100,1323 517,6346 54,1796
𝒎
−1 0 0 −1 0 −1
1 −1 0 0 0 0
0 1 1 0 0 0
𝐽=
0 0 −1 1 0 0
13,7084 113,5329 −22,8218 −21,3693 0 0
( 0 0 0 21,3693 −54,1608 −14,3376)
Realizando el proceso de cálculo iterativo, obtendremos la siguiente evolución:

⃗‖
‖Δ𝑄
It. Q1 [m3/s] Q2 [m3/s] Q3 [m3/s] Q4 [m3/s] Q5 [m3/s] Q6 [m3/s] 𝜀=
‖Q 𝒊 ‖
0 0,1000 0,1000 0,1000 0,1000 0,1000 0,1000
1 0,1360 0,0510 0,0490 0,1067 0,0523 0,1323 41,86%
2 0,1249 0,0399 0,0601 0,1287 0,0413 0,1213 14,25%
3 0,1239 0,0389 0,0611 0,1307 0,0405 0,1205 1,26%
4 0,1238 0,0388 0,0612 0,1308 0,0404 0,1204 0,08%
Anexos
Derivada de la función 𝑓 (𝑄) = 𝑟|𝑄|𝑄
La función 𝑓(𝑄) = 𝑟|𝑄|𝑄 donde r es una constante positiva, es una función monótona
creciente, por lo que su derivada es siempre positiva o 0.

𝑑𝑓 𝑑(𝑟 𝑄 𝑄)
ℎ𝑟 = 𝑓(𝑄) = 𝑟 𝑄 𝑄 = = 2𝑟 𝑄
𝑑𝑄 𝑑𝑄