Hidráulica Básica [2015961]

Tema # 2: Análisis de sistemas de tuberı́as

Luis Alejandro Morales (Ph.D)
Profesor Asistente

Universidad Nacional de Colombia-Bogotá

Facultad de Ingenierı́a
Departamento de Ingenieria Civil y Agrı́cola

Contents
1 Sistemas de tuberı́as simples 2
1.1 Tipos de problema en sistemas de tuberı́as . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Ecuaciones para la solución de problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3 Solucioón de la ecuación de Colebrook-White . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3.1 Metodo de punto fijo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3.2 Método de Newton-Raphson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.4 Comprobación de diseño . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.5 Cálculo de la potencia requerida . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.6 Disenõ de la tuberı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

2 Sistemas de tuberı́as en serie 7

2.1 Ecuaciones para la solución de problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.2 Comprobación de diseño . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.3 Cálculo de la potencia requerida . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.4 Disenõ de la tuberı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

3 Tuberı́as con perdida uniforme de caudal por unidad de longitud (Tubos porosos) 11

4 Sistemas de tuberı́as en paralelo 14

4.1 Ecuaciones para la solución de problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
4.2 Comprobación de diseño . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
4.3 Cálculo de la potencia requerida . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
4.4 Disenõ de tuberı́as en paralelo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

5 Redes abiertas de tuberı́as 17

5.1 Ecuaciones para la solución de problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
5.2 Comprobación de diseño . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
5.3 Disenõ de redes abiertas de tuberı́as . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20

6 Redes cerradas 21
6.1 Principios fundamentales del análisis de redes cerradas . . . . . . . . . . . . . . . . . . . . . 21
6.2 Método de Hardy-Cross con corrección de caudales . . . . . . . . . . . . . . . . . . . . . . . 24
6.3 Método de Hardy-Cross con corrección de cabezas . . . . . . . . . . . . . . . . . . . . . . . 25
6.4 Comparación de los métodos de Hardy-Cross . . . . . . . . . . . . . . . . . . . . . . . . . . 27
6.5 Método de análisis lineal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

1
1 Sistemas de tuberı́as simples

Figure 1: Acueductos por gravedad.

Una tuberı́a simple es aquella que tiene un diámetro y esta hecha de un solo material a lo largo de
su longitud (ver figura 2). La energı́a que mueve el flujo dentro de la tuberı́a es gracias a la acción de la
gravedad (tanque a la entrada) o a un máquina (sistema de bombeo a la entrada). Dichas tuberı́as pueden
tener cualquier tipo de accesorio a lo largo de su longitud lo que implica unas pérdidas menores. Las
ecuaciones de Prandl, Von-Karman y Darcy-Weisbach vistas en la Unidad 1, son utilizadas para el diseño
de tuberı́as simples. Note que existe cierta dificultad para el disenño teniendo en cuenta que la ecuación
de Colebrook-White para calcular el coeficiente de rugosidad f es implicita y requiere un proceso iterativo
para su solución. Los algoritmos que aquı́ se discutirán, constituyen las bases para el análisis y diseño de
tuberı́as más complejos.

Figure 2: Tuberı́a simple alimentada por un tanque de nivel constante y con descarga a la atmosfera
(tomado de [1]).

1.1 Tipos de problema en sistemas de tuberı́as

Los problemas en sistemas de tuberı́as se clasifican de acuerdo con las variables desconocidas. Las variables
involucradas en estos problemas se pueden clasificar como:
ˆ Caracterı́sticas la tuberı́a: Diámetro (D), longitud (L), rugosidad absoluta (ε).

ˆ Propiedades del fluido: Densidad (ρ) y viscosidad dinámica (µ) o cinemática (ν).

ˆ Variables relacionadas con el esquema del sistema: Coeficientes de pérdidas menores (K) de todos
los accesorios en el sistema.

ˆ Variables relacionas con la energı́a impulsora del sistema: Cabeza de energı́a (H = E1-E2), entre la
energı́a en el embalse de entrada (E1) y la energı́a salida del sistema (E2), o potencia de la bomba
(P ).

2
 ˆ Propiedades del flujo: Caudal (Q) y velocidad (V ) del flujo.

ˆ Otras variables: Aceleración de la gravedad (g).

De acuerdo con las variables involucradas en sistemas de tuberı́as, existen tres tipos de problemas:

1. Comprobación de diseño: En este tipo de problemas la tuberı́a existe y se conoce su longitud,

su diámetro, su rugosidad absoluta (material), al igual que todos los accesorios y sus coeficiente de
pérdidas menores. También se conoce la energı́a impulsora, ya sea una cabeza de energı́a (grav-
itacional por diferencia de niveles) o una energı́a mecánica (suministrada por una bomba). Las
propiedades del fluido como la densidad y la viscosidad absoluta son también conocidas. La incógnita
es entonce el caudal o la velocidad del flujo en el sistema.

Variables conocidas Incógnita

P
D, ε, H (o P), K, ρ, µ, g, L Q (o V)

2. Cálculo de la potencia requerida: En este tipo de problemas, el sistema existe por lo que se
conocen su longitud, su diámetro, su rugosidad absoluta (material), al igual que todos los accesorios
y sus coeficiente de pérdidas menores. Las propiedades del fluido como la densidad y la viscosidad
dinámica ası́ como el caudal (o velocidad) que fluye por el sistema son también conocidas. La
finalidad es determinar la potencia, ya sea mecánica o gravitacional, requerida para mover cierto
caudal a través de la tuberı́a dada.

Variables conocidas Incógnita

P
D, ε, Q (o V), K, ρ, µ, g, L P (o H)

3. Diseño de la tuberı́a: En este tipo de problemas se conoce el caudal o la velocidad de flujo y la

potencia disponible (mecánica o gravitacional), algunas caracterı́sticas de la tuberı́a como la longitud,
los accesorios y sus coeficientes de pérdida y las propiedades del fluido como la densidad y la viscosidad
dinámica. Se desconoce el diámetro necesario para permitir el paso de el caudal demandado. En
cuanto a la rugosidad absoluta, se debe cambiar el tipo de tuberı́a (rugosidad absoluta) con el fin de
obtener la mejor opción.

Variables conocidas Incógnita

P
ε, Q (o V), P (o H), K, ρ, µ, g, L D

1.2 Ecuaciones para la solución de problemas

A continuación se presentan las ecuaciones necesarias para resolver los tres problemas ya mencionadas.
Estas ecuaciones fueron discutidas en el capitulo anterior.
Si se tiene una tuberı́a simple cuya entrada es en la sección 1 y cuya salida es en la sección, aplicando
la ecuación de Bernoulli entre 1 y 2, de manera general, se tiene:

V1 2 p1 V2 2 p2 X
+ z1 + = + z2 + + hf + he − hb + ht (1)
2g γ 2g γ
donde hb es la cabeza de energı́a suministrada por la bomba y ht es la cabeza de energı́a sustraı́da por
2
la turbina. La energı́a total en una sección (e.g. 1 o 2) de flujo se puede expresar como E = V2g + z + γp ,
por lo tanto la ecuación 1 se puede expresar como:

E1 − E2 = hf + he − hb + ht (2)

3
 en donde hf son las pérdidas por fricción estimadas con la ecuación de Darcy-Weisbach como:

L V2
hf = f (3)
D 2g

y he son las perdidas por accesorios, las cuales se pueden calcular como:

X V2
he = K (4)
2g

El factor de fricción f en la ecuacioón 3, se calcula usando el diagrama de Moody o numéricamente

usando la ecuación Colebrook-White como:
 
1 ε 2.52
√ = −2 log + √ (5)
f 3.7D Re f

donde el número de Reynolds (Re) se calcula como

VD
Re = (6)
ν

Si se reemplaza las ecuaciones 3 y 4 en la ecuación 2, se tiene:

L V2 X V2
E1 − E2 = f + K − hb + ht (7)
D 2g 2g
despejando V en la ecuación 7, se tiene:
s
E1 − E2 + hb − ht
V = 2g L P (8)
fD + K

En los tres tipos de problemas, el objetivo es usar las ecuaciones 3, 4, 5, 6 y 8 para su solución. Note
que la ecuación 5 es una ecuación implı́cita que requiere del uso de algún método iterativo o numérico para
su solución. Note que la ecuación 8 se usa en particular para la solución de problemas de comprobación de
diseño y de diseño de tuberı́as.

1.3 Solucioón de la ecuación de Colebrook-White

1.3.1 Metodo de punto fijo
Consiste en el siguiente procedimiento:

1. Leer la información de entrada: ε, ρ, µ o ν, V o Q, D

2. Calcular el Re usando la ecuación 6

3. Si Re < 2000 (Flujo laminar), calcular f como:

f = Re/64 (9)
y luego ir a 6. Si Re > 2000 continuar.

4. Asumir un valor inicial de fi (e.g. f = 0.01).

4
 5. Usando la siguiente forma de la ecuación 5, calcular un valor fi+1 :
  −2
ε 2.52
fi+1 = −2 log + √ (10)
3.7D Re fi

6. Si |fi − fi+1 | ≤ η, donde η es un error (e.g. η = 1x10−6 ), ir a 6. Si |fi − fi+1 | > η, hacer fi = fi+1 e
ir a 6, para calcular un nuevo valor fi+1 .

7. Imprimir f

1.3.2 Método de Newton-Raphson

El método de Newton-Rapson es un método numérico para la solución de ecuaciones implı́citas. Converge
más rápido que el método de punto fijo. Algunas condiciones para poder aplicar el método son para un
intervalo l del dominio:
ˆ f (x) debe estar definida en l.

ˆ La función de iteración de f (x) deber ser continua en l.

ˆ La función implı́cita f (x) debe ser diferenciable (f ′ (x)) en l.

En general la ecuación de Colebrook-White cumple estas condiciones. El método consiste en el siguiente
procedimiento:

1. Leer la información de entrada: ε, ρ, µ o ν, V o Q, D

2. Calcular el Re usando la ecuación 6

3. Si Re < 2000 (Flujo laminar), calcular f como:

f = Re/64 (11)
y luego ir a 6. Si Re > 2000 continuar.

4. Asumir un valor inicial de fi (e.g. f = 0.01).

5. Calcular

1
xi = √ (12)
f
6. Calcular
 
ε 2.52xi
f (xi ) = −2 log + (13)
3.7D Re
7. Calcular

 " 2.52
#
′ −2 Re
f (xi ) = ε 2.52xi (14)
ln 10 3.7D + Re

8. Calcular el nuevo valor de xi+1 como:

f (xi ) − xi
xi+1 = xi − (15)
f ′ (xi ) − 1

5
 9. Si |xi − xi+1 | ≤ η, donde η es un error (e.g. η = 1x10−6 ), ir a 6. Si |xi − xi+1 | > η, hacer xi = xi+1
e ir a 6, para calcular un nuevo valor xi+1 .

10. Imprimir f , donde:

1
f= (16)
x2i+1

1.4 Comprobación de diseño

A continuación se describe el proceso para solucionar este tipo de problemas que consiste en determinar el
valor de V o Q.
P
1. Leer la información de entrada: ε, ρ, µ o ν, L, D, K, E1, E2, hb y ht .

2. Aplicar el procedimiento descrito en la sección 1.3.1 o en la sección 1.3.2. Note que lo único que
cambia dentro de estos procedimientos es que la velocidad V es calculada usando la ecuación 8; ya
no es un dato de entrada.

3. Imprimir V y f .

1.5 Cálculo de la potencia requerida

A continuación se describe el proceso para solucionar este tipo de problemas que consiste en determinar la
cabeza de energı́a de la bomba (hb ) y su potencia (P ).
P
1. Leer la información de entrada: ε, ρ, µ o ν, L, D, K, Q, E1, E2, ht y η.

2. Calcular la velocidad como V = Q/A.

3. Aplicar el procedimiento descrito en la sección 1.3.1 o en la sección 1.3.2 para calcular f .

4. Calcular hf usando la ecuación 3.

5. Calcular he usando la ecuación 4.

6. Calcular la cabeza de energı́a de la bomba (hb ) despejándola de la ecuación 2.

7. Calcular la potencia P como:

ρQghb
P = (17)
η
donde η es la eficiencia de la bomba.

8. Imprimir hb y P .

1.6 Disenõ de la tuberı́a

A continuación se describe el proceso para solucionar este tipo de problemas que consiste en determinar el
diámetro óptimo comercial de la tuberı́a.
P
1. Leer la información de entrada: ε, ρ, µ o ν, L, D, K, Q, E1, E2, ht y η.

2. Asumir un diámetro comercial inicial Di para la tuberı́a. El Di inicial debe ser pequeño (e.g. 1 pulg.)
pero no tan pequeño ya que el procedimiento no converge.

6
 3. Aplicar el procedimiento descrito en la sección 1.3.1 o en la sección 1.3.2. Note que lo único que
cambia dentro de estos procedimientos es que la velocidad V es calculada usando la ecuación 8; ya
no es un dato de entrada. De aquı́ sale un valor de f y V .
πDi2
4. Calcular el nuevo caudal Qn = V A usando Di , donde A = 4 .

5. Si Qn >= Q ir a 6. Si Qn < Q tomar el siguiente diámetro comercial superior Di+1 e ir a 3.

6. Imprimir D.

El diseño óptimo busca, en la mayorı́a de los casos, la solución más económica. Por lo tanto es necesario
muchas veces diseñar con otros materiales (diferente valor de ε) para encontrar la mejor solución.

2 Sistemas de tuberı́as en serie

Las tuberı́as en serie son dos o mas tuberı́as conectadas una tras de otra con diferente diámetro o rugosidad
(material) o ambos (ver figura 3). Estas tuberı́as son muy comunes en sistemas de riego localizado de
alta frequencia o en lı́neas de coducción para acueductos veredales. Al igual que en la sección 1, aquı́
se describirán las ecuaciones generales para resolver problemas relacionados con tuberı́as en serie y se
explicarán los procedimientos para resolver los tres problemas en tuberı́as: 1) Comprobación de diseño, 2)
cálculo de la potencia y 3) diseño de las tuberı́as.

Figure 3: Tres tuberı́as en serie conectando dos tanques en donde QLi representa un caudal lateral de
salida al final de la tuberı́a i (tomado de [1]).

2.1 Ecuaciones para la solución de problemas

Teniendo en cuenta las tuberias en serie de la figura 3, se plantean las siguientes ecuaciones:

ˆ Conservación de la energı́a

∆E = E1 − E2 = HT = z1 − z2 = he + hf1 + hm1 + hf2 + hm2 + hf3 + hm3 + hs (18)

donde HT es la diferencia de niveles entre los dos tanque (la energı́a total disponible en el sistema),
he pérdidas menores de entrada, hfi pérdidas por fricción en el tubo i, hmi pérdidas menores (por
válvulas, uniones, etc) en la tuberı́a i y hs pérdidas por salida.

7
 Para un numero n de tuberı́as en serie la ecuación 18 se puede escribir como:
n
X n X
X m
HT = he + hfi + hmi,j + hs (19)
i=1 i=1 j=1

donde m es numero de accesorios en la tuberı́a i; m puede ser variable. La ecuación 19 establece que
la energı́a disponible en el sistema se disipa en pérdidas a lo largo de la tuberı́a. Puede haber el caso
en el cual tengamos, por ejemplo, una tuberı́a horizontal por lo cual una bomba serı́a necesaria para
impulsar el flujo por lo que en la parte izquierda de la ecuación 18 tendrı́amos HT = hb . También se
puede dar el caso de una turbina en alguna sección de la tuberı́a por lo que habrı́a que sumarle a la
derecha de la ecuación ht .
Desarrollando la ecuación 19, se tiene:

n  n  2 Xm
V2 L V2
    2
X X V V
HT = Ke + f + Kj + Ks (20)
2g 1 D 2g i 2g i 2g n
i=1 i=1 j=1

ˆ Conservación de la masa o ecuación de continuidad:

QT = Q1 = Q2 + QL1 = Q3 + QL1 + QL2 (21)

Donde Qi es el caudal que viaja por la tuberı́a i y QLi es el caudal derivado de la tuberı́a i. La
ecuación 21 se puede expresar de forma mas compacta como:

i−1
X
QT = Qi + QLj (22)
j=1

En caso de que no existiera derivaciones (QLi = 0) la ecuación 22 queda como:

QT = Q1 = ... = Qi = ... = Qn (23)

2.2 Comprobación de diseño

Se desea calcular el valor de QT = Q1 . Antes de describir el proceso de cálculo, deduciremos una ecuación
que se utilizará para resolver este tipo de problemas.
De la ecuación de Darcy-Weisbach, despejando el factor de fricción se tiene que:

hf D2g
f= (24)
LV 2
Sacando raı́z cuadrada e invirtiendo los términos a ambos lados se tiene:
√
1 V L
√ =p (25)
f hf D2g
Igualando a la ecuación de Colebrook-White, se tiene:
√  
V L ε 2.52
p = −2 log + √ (26)
hf D2g 3.7D Re f

8
 VD
Reemplazando Re = ν y la ecuación 25 en la ecuación 26 y despejando la velocidad, se tiene:
p √ !
−2 2gDhf ε 2.52ν L
V = √ log + p (27)
L 3.7D D 2gDhf

Note que la ecuación 27 es explicita para la velocidad V .

A continuación se describe el proceso para solucionar este tipo de problemas que consiste en determinar
QT = Q1 .

1. Leer: ρ, µ o ν, Ke (coef. de pérdidas a la entrada),

P Ks (coef. de pérdidas a la salida), E1 y E2.
También, para cada tuberı́a i, leer: εi ,Li , Di , Ki , QLi , hbi , hti y ηi (eficiencia de la bomba o
turbina). Esta lectura se debe hacer para las n tuberı́as del sistema.

2. Para llevar a cabo este proceso es necesario, para la primera iteración, asumir un valor inicial de hf1 .
Se ha encontrado que hf ∝ DL5 , de acuerdo con esto, [1] ha establecido que este valor inicial de hf1
se puede expresar como:

L1 /D15
hf1 = HT Pn 5 (28)
i=1 Li /Di

3. Para la tuberı́a i = 1: calcular el valor de V1 usando la ecuación 27, calcular el caudal QT = Q1 =

V1 A1 , calcular el valor de f1 usando uno de los dos procedimientos descritos en la secciones 1.3.1 y
1.3.2, calcular hf1 utilizando la ecuación 3 y hm1 utilizando la ecuación 4.

4. Para el resto de tuberı́as i = 2...n: calcular el caudal Qi usando la ecuación 22, calcular la velocidad
Vi = Qi /Ai , fi usando uno de los dos procedimientos descritos en la secciones 1.3.1 y 1.3.2, calcular
hfi utilizando la ecuación 3 y hmi utilizando la ecuación 4.

5. Calcular la pérdida de energı́a total estimada HˆT usando la ecuación 19. Si |HˆT − HT | ≤ η donde η
es un error (e.g. η = 1x10−6 ), ir a 6. Si |HˆT − HT | > η, actualizar el valor de hf1 , como:

htf1 = ht−1
f1 + ∆hf1 (29)

donde ∆hf1 , se calcula como

L1 /D15
∆hfi = (HT − HˆT ) Pn 5 (30)
i=1 Li /Di

Una véz calculado el nuevo hf1 ir a 3.

6. Imprimir Qi donde i = 1...n.

2.3 Cálculo de la potencia requerida

A continuación se describe el proceso para solucionar este tipo de problemas que consiste en determinar la
cabeza de energı́a de la bomba (hb ) y su potencia (P ).

1. Leer: ρ, µ o ν, Ke (coef. de pérdidas a la entrada), Ks P (coef. de pérdidas a la salida), E1, E2 y

QT = Q1 . También, para cada tuberı́a i, leer: εi , Li , Di , Ki , QLi , hti y ηi (eficiencia de la bomba
o turbina). Esta lectura se debe hacer para las n tuberı́as del sistema.

2. Utilizando la ecuación 22, calcular el valor de Qi para i = 2...n.

9
 3. Para cada una de las tuberı́as i = 1...n, calcular: la velocidad como Vi = Qi /Ai , fi usando uno de
los dos procedimientos descritos en la secciones 1.3.1 y 1.3.2, calcular hfi utilizando la ecuación 3 y
hmi utilizando la ecuación 4.

4. Calcular la cabeza de energı́a de la bomba (hb ) despejandola de la ecuación 2. Note que en esta
ecuación, hf es la sumatoria de todas hfi y hm es la sumatoria de todas las hmi .

5. Calcular la potencia P usando la ecuación 17.

6. Imprimir hb y P .

2.4 Disenõ de la tuberı́a

Como se habı́a mencionado anteriormente, el diseño, y en este caso, de tuberı́as en serie es un problema
complejo por que eso implica buscar, en la mayorı́a de los casos, la solución óptima desde el punto de vista
económico. Por esto, no solo basta con determinar el diámetro si no también el material de la tuberı́a.
A esto hay que sumarle que los coeficientes de perdida dependen del material de la tuberı́a, esta es otra
variable que hay que tener en cuenta. Esto quiere decir que pueden existir multiples soluciones para un
mismo problema.
A continuación se describe el proceso para solucionar este tipo de problemas que consiste en determinar
el diámetro óptimo comercial Di de la tuberı́a para i = 1...n. En este proceso se asume que el material es
un dato dado.

1. Leer: ρ, µ o ν, Ke (coef. de pérdidas a la entrada),P Ks (coef. de pérdidas a la salida), E1, E2 y

QT = Q1 . También, para cada tuberı́a i, leer: εi , Li , Ki , QLi , hbi , hti y ηi (eficiencia de la bomba
o turbina). Esta lectura se debe hacer para las n tuberı́as del sistema.

2. Utilizando la ecuación 22, calcular el valor de Qi para i = 2...n.

3. Inicializar hfi para cada tuberı́a i = 1...n. I-pai Wu en 1975 encontró la siguiente ecuación para
inicializar hfi antes de iniciar las iteraciones:

Li cos θi
hfi = HT Pn (31)
i=1 Li cos θi

donde θi es el angulo formado entre la horizontal y la tuberı́a. Usualmente este angulo es el mismo
z1 −z2 Pn
para todas las tuberı́as ya que todas tienen la misma pendiente, por lo tanto θ = / i=1 Li ,
donde z1 es la altura a la entrada del sistema y z2 es la altura a la salida del sistema.

4. Para cada tuberı́a i = 1...n

(a) Asumir un diámetro comercial inicial Di . El Di inicial debe ser pequeño (e.g. 1 pulg.) pero no
tan pequeño ya que el procedimiento no converge.
(b) Aplicar el procedimiento descrito en la sección 1.3.1 o en la sección 1.3.2. Note que lo único que
cambia dentro de estos procedimientos es que la velocidad V es calculada usando la ecuación 8;
ya no es un dato de entrada. De aqui sale un valor de f y V .
πDi2
(c) Calcular el nuevo caudal Qni = Vi Ai usando Di , donde Ai = 4 .
(d) Si Qni >= Qi ir a 5. Si Qni < Qi tomar el siguiente diámetro comercial superior Di+1 e ir a
4a.

5. Para cada tuberı́a i = 1...n, calcular la velocidad real como VRi = Qi /Ai

6. Para cada tuberı́a i = 1...n, y con el Di y VRi , calcular fi usando el procedimiento descrito en la
sección 1.3.1 o en la sección 1.3.2.

10
 7. Para cada tuberı́a i = 1...n, calcular hfRi utilizando la ecuación 3 y hmRi utilizando la ecuación 4.
Usar en estas ecuaciones Di y VRi .

8. La energı́a del sistema a la salida de la tuberı́a n suele ser superior a la energı́a requerida del sistema
allı́ (E2), por lo que usualmente se utiliza una valvula para hacer caer la energı́a a E2. Dicha caida
hmv , se calcula como:
Xn X n
hmv = HT − hfRi − hmRi (32)
i=1 i=1

9. Si hmv > 0 y hmv ≈ E2, ir a 10. Si hmv >> E2, para cada tuberı́a i = 1...n, actualizar el valor de
hfi para la siguiente iteración como:
hfi = hfRi + hmv (33)
Luego ir a 4.

10. Imprimir Di para i = 1...n.

Es importante revisar que hmv sea positivo, de lo contrario esto supondrı́a presiones negativas que
podrı́an generar cavitación en el sistema. El procedimiento anterior puede ser repetido para diferentes
materiales de tuberı́a (diferentes valores de ε) con el fin de escoger el diseño más económico.

3 Tuberı́as con perdida uniforme de caudal por unidad de longitud

(Tubos porosos)
En sistemas de riego moderno existen sistemas de tuberı́as que estan diseñados para entregar un caudal
por unidad de longitud con el objetivo de lograr uniformidad en el riego (ver figura 5). Estas tuberı́as se
conocen como tubos porosos o mangueras exudantes.

Figure 4: Sistema de riego con tuberı́a porosa.

Si se analiza la tuberı́a de la figura 5, se sabe que la presión disminuye a lo largo de la tuberı́a debido
a las pérdidas asi como la velocidad ya que el caudal disminuye por unidad de longitud. De la figura 5,
tenemos:

ˆ Conservación de la masa
Qx = Q1 − qx (34)
donde Qx es el caudal en la tuberı́a a una distancia x de la sección de entrada 1, q es el caudal
descargado por la tuberı́a por unidad de longitud, Q1 es el caudal a la entrada (1) de la tuberı́a.

11
Figure 5: Tuberı́a porosa con entreaga de caudal por unidad de longitud (q) y su lı́nea de gradiente
hidráulico (tomado de [1]).

ˆ Pendiente de la lı́nea de gradiente hidráulico Aplicando la ecuación de Darcy-Weisbach (para pérdidas

por fricción), e ignorando pérdidas menores, para una longitud infinitesimal de tuberı́a dx, en donde
no existe salida de caudal, se tiene:
dx Vx2
dhf = fx (35)
D 2g
donde dhf es la pérdida de energı́a a lo largo de dx, fx es el factor de fricción en dx y Vx es la
velocidad en dx. Expresando la ecuación 35 en términos del caudal Qx , se tiene:

dhf fx Q2x
= (36)
dx 2g DA2

dhf 8fx Q2
= 2 x5 (37)
dx π gD
reemplazando la ecuación 34 en la ecuación 37, se tiene:

dhf 8
= 2 5 fx (Q1 − qx)2 (38)
dx π gD

Integrando la ecuación 38 para calcular las pérdidas a lo largo de la longitud de la tuberı́a (L), se
tiene:

Z L
8
hf = 2 5 fx (Q1 − qx)2 dx (39)
π gD 0

donde fx esta dado por la ecuación de Colebrook-White:

 
1 ε 2.52ν
√ = −2 log + √ (40)
fx 3.7D Vx D fx

12
 Note que la ecuación 39 no se puede resolver analiticamente ya que el valor de fx , coeficiente de
fricción a una distancia x, se calcula a través de la ecuación implicita 40. Sin embargo la ecuación 39
se puede resolver aproximadamente utilizando los métodos siguientes:

1. Factor de fricción constante

Si se supone que fx es constante y se calcula como el promedio del factor de fricción a la entrada
(f1 ) y a la salida (fL ):

f1 + fL
fx = (41)
2
integrando la ecuación 39 para x, se tiene:

x=L
q 2 x3
 
8fx 2 2
hf = 2 5 Q1 x − Q1 qx + (42)
π gD 3
x=0

en donde se obtiene:

q 2 L3
 
8fx 2 2
hf = 2 5 Q1 L − Q1 qL + (43)
π gD 3

2. Factor de fricción constante tramo a tramo

Para utilizar este método se supone que la tuberı́a de longitud L está dividida en n tramos de
igual longitud l = L/n. Por lo tanto se calculan las pérdidas para cada tramo de la tuberı́a
suponiendo que fx es constante en el tramo, y luego se suman. Las pérdidas de energı́a hfi para
un tramo i donde i = 1 · · · n son:

q 2 l3
 
8fx 2 2
hfi = 2 i 5 Q1i l − Q1i ql + (44)
π gD 3

donde Q1i es el caudal transportado por el tramo i e igual a Q1i = Q1 − ql(i − 1). Note que fxi
es el promedio entre tramos:

fi + fi+1
fxi = (45)
2
Las pérdidas totales a lo largo de la tuberı́a se calculan como:
n
X
hf = h fi (46)
i=1

El procedimiento para el calculo de la pérdida de energı́a en una tuberı́a porosa de describe a contin-
uación.

1. Calcular el caudal de entrada a la tuberı́a (Q1 ) de longitud L como:

Q1 = qL (47)

2. Cálculo de la zona laminar. Debido a que la manguera o tuberı́a va perdiendo caudal su velocidad
disminuye y el flujo, en algún punto, se convierte en laminar. Para calcular ese punto, se asume un
número de Reynolds crı́tico (donde cambia a laminar) igual a Rec = 2000. Dicho caudal, para el cual
ocurre Rec , se calcular como:

13
 Rec πDν
Qc = (48)
4
la distancia desde el final de la tuberı́a en donde el flujo cambia a laminar, se calcula como:

Lc Qc
Lc = (49)
Q1

Teniendo en cuenta que en los últimos tramos de la tuberı́a (x < xc ) el número de Reynolds tiende a
cero, el f tiende a infinito. Por esto, es necesario dejar cierta longitud de la tuberı́a en donde no se
cálculan pérdidas. En estos ultimos tramos, las pérdidas son muy pequeñas por lo tanto no influyen
en el resultado final. Dicho valor puede ser el 0.3% de L.

3. Implementar cualquiera de los dos métodos descritos: 1) factor de fricción constante para la tuberı́a
y 2) factor de fricción constante tramo a tramo.

4 Sistemas de tuberı́as en paralelo

Las tuberı́as en paralelo son un conjunto de dos tuberı́as o mas las cuales se derivan de un nodo común
(nodo de entrada) y convergen hacia otro nodo (nodo de salida). Esto implica que, de acuerdo a la ley
de conservación de la masa, la suma de los caudales que viajan por cada una de las tuberı́as debe ser
igual aguas arriba del nodo de entrada y aguas abajo del nodo de salida; esto implica que el caudal de
entrada y de salida de las tuberı́as es le mismo. En casos prácticos y cotidianos, los sistemas de tuberı́as
en paralelo están limitados a máximo 3 o 4 tuberı́as, y son mas comunes los sistemas de dos tuberı́as. En
estos sistemas, cada tuberia puede tener longitudes, diámetros, accesorios y materiales diferentes.

Figure 6: Aumento de la capacidad de un acueducto veredal mediante la construcción de una tercera

tuberı́a.

4.1 Ecuaciones para la solución de problemas

Teniendo en cuenta las tuberı́as en paralelo de la figura 7, se plantean las siguientes ecuaciones:

14
Figure 7: Esquema tridimensional de dos tuberı́as en paralelo y sus lineas de energı́a, con un nodo de
entrada 1 y un nodo de salida 2. (tomado de [1]).

Conservación de la energı́a
∆E = E1 − E2 = HT = h1f1 + h1m1 + h1f2 + h1m2 + h1f3 + h1m3 + h1f4 + h1m4 + h1f5 (50)
donde HT es la diferencia de energı́a entre el node de entrada (E1) y el nodo de salida (E2); la energı́a
total disponible en el sistema, hfi ij son las pérdidas por fricción en el tramo i = 1...5 de la tuberı́a j = 1
y hjmi son las pérdidas menores (por válvulas, uniones, etc) en el tramo i = 1...4 de la tuberı́a j = 1.
La ecuación 50 puede expresarse de forma más compacta como:
nj mj
hjfi
X X
HT = + hjmi (51)
i=1 i=1

donde nj es el número de tramos de la tuberı́a j = 1 y mj es el número de accesorios en la tuberı́a j = 1.

Note que la ecuación 51 aplica también para la tuberı́a j = 2.
De acuerdo con la ley de conservación de la energı́a, esto implica que las pérdidas de energı́a en la
direccón del flujo a través de la tuberı́a 1 y 2 es la misma. De acuerdo con esto, la ecuación 51, se escribe
como:
Xn1 Xm1 n2
X m2
X
1 1 2
HT = hfi + hmi = hfi + h2mi (52)
i=1 i=1 i=1 i=1

Conservación de la masa
De acuerdo con la ley de conservación de la masa, el caudal que entra al nodo de entrada (QT ), debe ser
igual a la suma de los caudales que se distribuyen en las tuberı́as e igual al caudal en el nodo de salida.
Para el esquema de la figura 7, se tiene que:

QT = Q1 + Q2 (53)
donde Q1 es el caudal que transporta la tuberı́a 1 y Q2 es el caudal que transporta la tuberı́a 2. En
general para un numero n de tuberı́as en paralelo, se tiene:
n
X
QT = Qi (54)
i=1

15
4.2 Comprobación de diseño
Para este caso se conocen las caracterı́sticas fı́sicas de n tuberı́as, los accesorios en cada una de ellas, la
energı́a disponible en el sistema y las propiedades del fluido circulante. Se desea entonces conocer los
caudales individuales en cada tuberı́a (Qi ). A continuación se presenta el procedimiento de cálculo:
1. Leer:
P ρ, µ o ν, hb (Cabeza de la bomba), E1 y E2. También, para cada tuberı́a i, leer: εi , Li , Di ,
Ki , hti y ηi (eficiencia de la turbina). Esta lectura se debe hacer para las n tuberı́as del sistema.

2. Para cada una de las n tuberı́as en paralelo, hacer la comprobación de diseño para tuberı́as simples
descrito en la sección 1.4. Se obtendrán n caudales Q y n velocidades V para cada tuberı́a.

3. El caudal total QT se obtendra usando la ecuación 54.

4.3 Cálculo de la potencia requerida

A continuación se describe el proceso para solucionar este tipo de problemas que consiste en determinar la
presión en el nodo de entrada (1) o en el nodo de salida (2) de las tuberı́as en paralelo. Con esto es posible
calcular la potencia consumida por el flujo a través de las tuberı́as en paralelo.
Uno de los problemas es que, a pesar de que el caudal total QT es un dato dado en el problema, no se
sabe como éste se distribuye a través de las tuberı́as. Para esto se ha encontrado la siguiente ecuación a
partir de la ecuación de Darcy-Weisbach:
5/2 √
d1 / l1
Qk1 = QT P 5/2 √
(55)
n
i=1 di / li
en donde l1 es la londitud de la tuberı́a 1 y d1 es el diámetro de la tuberı́a 1. La ecuación 55 sirve
para inicializar el caudal Qk1 en la tuberı́a 1. Con este caudal Qk1 es posible calcular la energı́a (presión)
en el node 1 o nodo 2 siguiendo el algoritmo de la sección 1.5. Una vez se tiene la energı́a disponible en
el sistema (E1 − E2 ), es posible calcular el caudal en cada una de las otras n − 1 tuberı́as siguiendo el
algoritmo de la sección 1.4. Sin embargo, es posible que la sumatoria de los n caudales no se sea igual a
QT , por lo que es necesario corregir la suposicion inicial de Q1 como:

QT
Qk+1
1 = Qk1 (56)
Q∗T
donde
n
X
Q∗T = Qki (57)
i=1
es el caudal total estimado para la iteración k. El proceso de iteracion termina cuando QT = Q∗T .
A continuación se describe el algoritmo para encontrar la presión (energı́a) en el nodo de entrada o de
salida al sistema y la potencia disipada.

1. Leer: ρ,P µ o ν, hb (Cabeza de la bomba), QT y E1 o E2. También, para cada tuberı́a i, leer: εi ,
Li , Di , Ki , hti y ηi (eficiencia de la turbina). Esta lectura se debe hacer para las n tuberı́as del
sistema.

2. Inicializar el caudal transportado por la tuberı́a 1 (Q1 ) usando la ecuación 55.

3. Con base en el caudal Q1 determinar la energı́a en el nodo 1 (E1 ) o en el nodo 2 (E2 ) siguiendo el
algoritmo de la sección 1.5.

4. Una vez calculada la energı́a en el nodo 1 (E1 ) o en el node 2 (E2 ) para cada una de las n − 1 tuberı́as
restantes, calcular el caudal Q siguiendo el procedimiento descrito en la sección 1.4.

16
 5. Una vez calculado el caudal Q en cada tuberı́a, calcular el caudal total estimado Q∗T para la iteración
k. Si |QT − Q∗T | ≤ η, donde η es un error (e.g. η = 1x10−6 ), ir a 6. Si |QT − Q∗T | > η, corregir Q1
usando la ecuación 56 e ir 3.

6. Imprimir E1 o E2 y Qi donde i = 1 · · · n. Con base en las pérdidas a lo largo de la tuberı́a i, calcular

la potencia disipada en el sistema.

4.4 Disenõ de tuberı́as en paralelo

En la practica en ingenierı́a no es común diseñar sistemas de tuberı́as en paralelo debido a su ineficiencia
a nivel hidráulico ya que el perı́metro mojado de 2 o mas tuberı́as es mayor que el perı́metro mojado de
una sola tuberı́a. Este diseño solo se justifica si el diámetro de diseño es mayor al diámetro comercial
disponible. De ser ası́, cada tuberı́a se diseña como una tuberı́a simple para un caudal QT /n, donde n
es el número de tuberı́as en paralelo. Sin embargo, lo que si es común es la ampliación de un sistema
conformado por una tuberı́a simple ya existente debido a un incremento en la demanda, lo cual es común
en acueductos y sistemas de riego. Para esta ampliación se diseña y se construye una tuberı́a paralela a
la ya existen para la cual se conocen sus caracterı́sticas. Para este diseño se conoce el nuevo caudal total
QT , la energı́a en el nodo de entrada al sistema y las caracterı́sticas de las tuberı́as. Las incógnitas son el
diámetro de la nueva tuberı́a y la energı́a en el nodo de salida del sistema. Note que la energı́a en el nodo
final será diferente a la energı́a inicial para este nodo (antes de la ampliación) ya que la hidráulica de la
tuberı́a existente se modifica también.
A continuación se presenta un algoritmo para el diseño de la nueva tuberı́a (número 2), en donde la
tuberı́a 1 es existente.

1. Leer: ρ, µ o ν, hb (Cabeza de la bomba), Q1 , QT , E1 (energı́a del nodo a la entrada) y E2 (energı́a

del nodo Pa la salida. Esta energı́a cambia durante el proceso). También, para cada tuberı́a i, leer:
εi , Li , Ki , hti y ηi (eficiencia de la turbina) donde i = 1 · · · 2. Leer el diametro D de la tuberı́a
existente.

2. Calcular el caudal por la tuberı́a 2: Q2 = QT − Q1 .

3. Diseñar la tuberı́a 2 como una tuberı́a simple siguiendo el procedimiento descrito en la sección 1.6.
Teniendo en cuenta que, debido a que el diámetro comercial siempre deja pasar mas caudal, Q2 es
mayor al caudal deseado (caudal inicial), la energı́a en el nodo 2 aumenta y por lo tanto el caudal
Q1 cambia también.

4. De aquı́ en adelante se vuelve un problema de cálculo de la potencia requerida para dos tuberı́as en
paralelo partiendo del numeral 2 en el procedimiento descrito en la sección 4.3.

5 Redes abiertas de tuberı́as

En la practica, los sistemas de tuberı́as suelen ser mas complejos que tuberı́as simples, en serie o en
paralelo, por lo que las tuberı́as suelen estar conectadas unas con otras formando diferentes configuraciones.
Las tuberı́as ramificadas o redes abiertas de tuberiı́as suelen ser las tuberı́as principales (lı́neas madres o
expresas) en sistemas de acueductos, sistemas de irrigación o sistemas de tratamiento de agua, que están
compuestas ademas por tanques de almacenamiento, desarenadores, plantas de tratamiento y obras de
captación (ver 9). En muchos casos estas redes abiertas estń conectadas a redes cerradas.
Los problemas que se estudiaran para redes cerradas son el problema de comprobación de diseño, el cual
busca determinar los caudales en cada una de las tuberı́as, y el de diseño de la red en donde se determina
el diámetro de cada tuberı́a. Note que el problema de cálculo de la potencia no es pertinente ya que la
energı́a del sistema siempre es un dato dado que es determinado por la localización topográfica de la red
o por la energı́a suministrada por un bomba en una tuberı́a.

17
 Figure 8: Red de embalses de para el sumunistro de agua en Bogotá D.C.

5.1 Ecuaciones para la solución de problemas

Teniendo en cuenta la figura 10, se plantean las siguientes leyes de conservación para la solución de prob-
lemas de redes abiertas.

Conservación de la masa
En una red abierta es necesario conservar la masa a lo largo de toda la red. Se plantea entonces que en cada
nodo se conserva la masa, teniendo en cuenta las posibles demandas desde los tanques de almacenamiento
y desde los nodos. Analizando la red de la figura 10 se tiene que el tanque de distribución A conduce
el caudal total demandado a través de la tuberı́a madrea AU hacia el nodo U de donde se distribuye el
caudal hacia los tanques B,C y D, de donde se derivan unas demandas. De acuerdo con esto, la ecuación
de conservación de la masa en el nodo U es:
n
X
QiU − QLU = 0 (58)
i=1

donde QiU es el caudal en la tuberı́a iU , QLU es el caudal demandado en el nodo U que puede ser cero,
y n es el número de tuberı́as conectadas con el nodo U .

Conservación de la energı́a
De acuerdo con la figura 10, la ecuacion de conservación de la energı́a entre un tanque i = A, B, C, D y un
nodo U , se expresa como:
  2
LiU X ViU
|Zi − ZU | = HtiU = fiU + KiU (59)
DiU 2g

18
 Figure 9: Red abierta compuesta por tuberı́as, tanques de almacenamiento y nodos (tomado de [1]).

Figure 10: Red abierta conectando un tanque de suministro A, al cual llega un caudal total demandado,
con los tanques de almacenamiento B,C y D, de los cuales se demandan caudales (tomado de [1]).

donde HtiU es la energı́a disponible entre iU , ZU es la cabeza piezométrica en la unión (siempre es

desconocida), Zi es la cabeza piezométrica en el tanque i, ViU es la velocidad de flujo en la tuberı́a iU , diU
es el diámetro en laPtuberı́a iU , LiU es la longitud de la tuberı́a iU , fiU es el factor de fricción de Darcy
en la tuberı́a iU y KiU es el coeficiente de global de pérdidas menores en la tuberı́a iU .
La ecuación 59 se puede expresar en términos del caudal como:

KiU Q2iU
 P 
LiU
|Zi − ZU | = HtiU = fiU + (60)
2gDiU 2g A2iU

5.2 Comprobación de diseño

En este caso se conocen las caracterı́sticas del fluido transportado, las caracterı́sticas de las tuberı́as
(diámetros, rugosidades y longitudes), las cotas de los tanques que integran el sistema y las posibles
demandas de caudal en los nodos. Se desea conocer el caudal total demandado y las demandas de caudal
en cada uno de los tanques QLi . El proceso iterativo se basa en el proceso de comprobación de diseño para
tuberı́as simples aplicado a cada tuberı́a del sistema. Sin embargo como las cabeza de energı́a (ZUj ) en
los nodos no son conocidas estas se deben ir corrigiendo en el proceso de iteración hasta cumplir (aproxi-
madamente) con la ecuación 58 en cada nodo. A continuación se presenta un algoritmo para resolver este
tipo de problemas:
P
1. Leer: ρ y µ o ν. Para cada tuberı́a i = 1 · · · n, leer: εi , Li , Di , Ki , hpi y ηi (eficiencia de la
bomba). Para cada embalse o tanque de almacenamiento i leer la cota piezométrica Zi . Para cada
nodo j = 1 · · · m, leer la demanda de caudal en el nodo QLUj .

19
 2. Estimar los valores iniciales de ZUj . Esta estimaci´’on se debe hacer cuidadosamente ya que el algo-
ritmo puede converger hacia una solución errada (quedar atrapado en un minimo local de la función
objetivo). En el caso en donde no existan bombas en las tuberı́as los valores de ZUj pueden asumirse
por debajo del tanque principal y por encima de los tanques a los que se distribuye el flujo. En el
caso de un tramo con un sistema de bombeo ZUj puede asumirse por debajo del tanque al cual se
esta bombeando.

3. Para cada tuberı́a i = 1 · · · n, calcular el caudal siguiendo el algoritmo de comprobación de diseño

para tuberı́as simples descrito en la sección 1.4.

4. Para cada nodo j = 1 · · · m comprobar que | ni=1 Qij − QLUj | ≤ η, donde η es un error (e.g. η =
P

1x10−6 ). Note que Qij es positivo si entra al nodo y es negativo si sale del nodo. Si esta condición
no se cumple se debe corregir ZUj de la siguiente manera:

ZUt j = ZUt+1
j
+ ∆ZUj (61)

donde ∆ZUj es la correción que se debe calcular siguiendo la metodologı́an de R.J. Cornish (1940)
como:
P 
n
2 i=1 Qij − QLUj
∆ZUj = Pn Qij
(62)
i=1 Zi −ZUj

Note que en la medida que la condicioón | ni=1 Qij − QLUj | se va acercando a cero, el valor de ∆ZUj
P

también se acerca a cero. El proceso terminar cuando m

P
j=1 |∆ZUj | ≤ η.

5. Imprimir Qij , donde i = 1 · · · n y j = 1 · · · m. También es conveniente imprimir ZUj .

5.3 Disenõ de redes abiertas de tuberı́as

En este caso se conocen las caracterı́sticas del fluido transportado, las caracterı́sticas de las tuberı́as (ru-
gosidades y longitudes), las cotas de los tanques que integran el sistema y las posibles demandas de caudal
en los nodos. Se desea conocer el diámetro de cada tuberı́a Di . El proceso iterativo se base en el proceso
de comprobación de diseño para tuberı́as simples aplicado a cada tuberı́a del sistema, asumiendo previa-
mente los diámetros de cada tuberı́a. Sin embargo como las cabeza de energı́a (ZUj ) en los nodos no son
conocidas estas se deben ir corrigiendo en el proceso de iteración hasta cumplir (aproximadamente) con la
ecuación 58 en cada nodo. A continuación se presenta un algoritmo para resolver este tipo de problemas:
P
1. Leer: ρ y µ o ν. Para cada tuberı́a i = 1 · · · n, leer: εi , Li , Ki , hpi , ηi (eficiencia de la bomba) y
Qi . Para cada embalse o tanque de almacenamiento i leer la cota piezométrica Zi . Para cada nodo
j = 1 · · · m, leer la demanda de caudal en el nodo QLUj .

2. Estimar los valores iniciales de ZUj . Esta estimaci´’on se debe hacer cuidadosamente ya que el algo-
ritmo puede converger hacia una solución errada (quedar atrapado en un mı́nimo local de la función
objetivo). En el caso en donde no existan bombas en las tuberı́as los valores de ZUj pueden asumirse
por debajo del tanque principal y por encima de los tanques a los que se distribuye el flujo. En el
caso de un tramo con un sistema de bombeo ZUj puede asumirse por debajo del tanque al cual se
esta bombeando.

3. Disenar cada tuberı́a i = 1 · · · n de la red abierta siguiento el algoritmo de diseño para tuberias
simples descrito en la sección 1.6.

20
 4. Para cada nodo j = 1 · · · m y con base en los caudales (Qij ) calculados (generalmetne diferentes a
los leidos en el paso 1) para cada tuberı́a en el paso anterior, determinar la corrección de cabeza de
energı́a para cada nodo con base en la ecuación 62.
5. Verificar que la condición m −6
P
j=1 |∆ZUj | ≤ η donde donde η es un error (e.g. η = 1x10 ). Note que
|∆ZUj = ZUt j − ZUt+1
j
|. Si la condición anterior se cumple ir a 6. Si la condición no se cumple ir a 3.
6. Imprimir Di , donde i = 1 · · · n.

6 Redes cerradas

Figure 11: Red cerrada del laboratorio de hidráulica.

Las redes cerradas son también conocidas como redes de distribución (en sistemas de abastecimiento de
agua). Estas redes se caracterizan por tener al menos un circuito cerrado garatizando que cualquier zona
de la red sea alcanzada por una tuberı́a, lo cual incrementa la confiabilidad del abastecimiento. Al igual
que en los casos de sistemas de tuberı́as ya vistos, las ecuaciones para el análisis de redes cerradas se basan
en las ecuaciones de conservación de la energı́a y de la masa, ası́ como de las ecuaciones para el calculo
de las pérdidas. Se presentan aquı́ los diferentes métodos de análisis para la comprobación de diseño, en
donde se calculan los caudales en las tuberı́as y las presiones en los nodos de las redes. Los métodos se
presentan en orden cronológico, en donde los métodos más antiguos son menos eficientes desde el punto de
vista computacional.

6.1 Principios fundamentales del análisis de redes cerradas

Ecuaciones de conservación en los nodos
Conservación de la masa
Si se considera la red cerrada de la figura 13 en donde QD1 , QD2 , QD3 , · · · , QDN u son los caudales
demandados en los nodos y Qe1 , Qe2 , · · · , Qem son los caudales que alimentan la red, la ecuacion de
continuidad se puede escribir como:
m
X Nu
X
Qe = QD (63)
i=1 i=1

21
 Figure 12: Red cerrada de tuberı́as de acuerdo con la norma Colombiana (Insfopal) (tomado de [1]).

donde m es el número de nodos con entrada de caudales a la red y N u es el número de nodos con
demandas. La ecuación de continuidad para cada nodo, se puede expresar como:
N Ti
X
Qij − QDi = 0 (64)
j=1

donde N Ti es el número de tubos que llegan al nodo i y Qij es el caudal que pasa por la tuberı́a ij
hacia el nodo i desde el nodo j. Qij es positivo si entra a i y negativo si sale de este.

Conservación de la energı́a

Para cada tramo ij, se puede expresar la ecuación de conservación de la energı́a como:

X Vij2
Hj − Hi = Kmij + hfij (65)
2g
En términos de caudal, la ecuación 65 se expresa como:

Q2ij
X 
Lij
Hj − Hi = Kmij + fij (66)
2gA2ij Dij
despejando Qij de la ecuación 66, se tiene:
 1
2
Hj − Hi p
Qij =  P L
 2gAij (67)
Kmij + fij Dijij

Con el fin de tener en cuenta el signo de los caudales, la ecuación anterior se expresa como

22
Figure 13: Red cerrada de tuberı́as en donde se indican los caudales demandados y los caudales que la
alimentan (tomado de [1]).

 
Hj − Hi 1
 |Hj − Hi |− 2
p
Qij =  P L
2gAij (68)
Kmij + fij Dijij
Reemplazando la ecuación 67 en la ecuación 64, se tiene:

 1
N Ti 2

 P Hj − Hi
X p
L
 2gAij − QDi (69)
j=1 Kmij + fij Dijij

Las ecuaciones 69 se conocen como las ecuaciones de cabeza. Como generalmente es conoce una cabeza
de energı́a en uno de los nodos, se tienen N u − 1 ecuaciones de cabeza en la red, donde N u es el numero
de nodos de la red.

Ecuaciones de conservación en los circuitos

Para los circuitos que conforman la red, se deben cumplir las leyes de conservación.
Conservación de la masa
Para cada nodo que conforma el circuito se debe cumplir la ecuación 64.
Conservación de la energı́a
Para cada circuito, las pérdidas alrededor del circuito deben ser iguales a cero, es decir:
N Ti′ N Ti′
X X
hfij + hmij = 0 (70)
j=1 j=1

donde N Ti′ es el número de tuberı́as en el circuito i. La ecuación anterior se pude expresar como:

N Ti′
X Q2ij X Lij

Kij + fij (71)
j=1
2gA2ij Dij

23
 Con el fin de tener en cuenta los signos de los caudales en las tuberı́as, la ecuación 71 se convierte en:
N Ti′
X Qij |Qij | X Lij

2 Kij + fij (72)
2gA ij Dij
j=1

Las ecuaciones 71 son conocidas como las ecuaciones de caudal. Se tienen entonces N C ecuaciones de
caudal, donde N C es el numero de circuitos en la red. Note que pueden existir circuitos que contienen
otros y en todos ellos se deben cumplir las ecuaciones de caudal. Sin embargo el análisis se hace para
circuitos independientes.
En total se tienen entonce N C + (N U − 1) ecuaciones no lineales las cuales sirven para encontrar los
caudales en las tuberı́as y las cabezas de energı́a en los N U − 1 nodos de la red dadas las caracterı́sticas
de las tuberı́as y de fluido circulante. Esto quiere decir que el problema se convierte entonce en una
comprobación de diseño más que en un diseño de la red. Es posible para diferentes combinaciones de
diámetros y mediante métodos de optimización de los costos de la red obtener los diámetros en la red.

6.2 Método de Hardy-Cross con corrección de caudales

Este método desarrolado por H. Cross en 1936, resuelve las ecuaciones 69 y 71 para una red cerrada de
tuberı́as dadas las caracterı́sticas de las tuberı́as y del fluido. Este método corrige en cada iteración los
caudales en las tuberı́as. Si se supone un error ∆Qi en el circuito i, las perdidas para un tubo se expresan
como:

(Qij + ∆Qi )2
 
Lij X
hfij + hmij = fij + Kmij (73)
Dij 2gA2ij
Si la cabeza de energı́a en la tuberı́a ij se define como hij = hfij + hmij , la ecuación 73 se expresa como:
 2
(Qij + 2Qij ∆Qi + ∆Q2i )

Lij X
hij = fij + Kmij (74)
Dij 2gA2ij
Despreciando el término ∆Q2ij en la ecuación 74 ya que es pequeño en relación con los otros dos términos,
se tiene:

(Q2ij + 2Qij ∆Qi )

 
Lij X
hij = fij + Kmij (75)
Dij 2gA2ij
Para el numero de tubos en el circuito i N Ti′ , se tiene:
N Ti′ N Ti′
X Q2ij X Lij
 X Qij X Lij

Kij + fij + 2∆Qi Kij + fij (76)
j=1
2gA2ij Dij
j=1
2gA2ij Dij

Despejando ∆Qi y simplificando la ecuación 76, tenemos:

P
(hf + hmij )
∆Qi = − P  hij +hm  (77)
fij ij
2 Qij

24
Comprobación de diseño y diseño
Pasos para la comprobación de diseño y para el diseño usando el método de Hardy-Cross con correción de
caudales.

1. Definición de la topologı́a de la red. Los circuitos están conformados por tuberı́as y la conección
de las tuberı́as en los circuitos es positiva siguiendo las manecillas del reloj y negativo en sentido
contrario. Todos los tubos deben quedar incluidos en almenos un circuito. Definir los nodos con
cabeza conocida.

2. Si existe más de un nodo con cabeza conocida es necesario conectarlos a través de una tuberı́a
imaginaria formando un circuito que los una. En estas tuberı́as imaginarias se deben suponer sus
caracterı́sticas con el fin de calcular el caudal hipotético que circuları́a por ellas. Sin embargo en las
correcciones de caudal, los tubos imaginarios no deben ser incluidos.

3. Se suponen unos caudales iniciales en cada una de las tuberı́as de la red. Recuerde que los caudales
que van en el sentido de las manecillas del reloj son positivos, de lo contrario son negativos. Recuerde
ademas que los caudales que entran a un nodo son positivos y los que salen son negativos.

4. Para cada circuito, calcular las pérdidas a lo largo de cada tuberı́a utilizando la ecuación 71.

5. Se calcula la sumatoria de las pérdidas a lo largo del circuito, teniendo en cuenta que los flujos en
sentido de las manecillas del reloj arrojan perdidas positivas y los flujos en sentido contrario dan
pérdidas negativas. Si la sumatoria es diferente de cero, se procede a calcular la corrección de caudal
para las tuberı́as del circuito con la ecuación 77. Note que al balance de perdidas se deben restar las
cabezas de bombas que existiesen en alguna tuberı́a del circuito.

6. Los pasos 4 y 5 se repiten para el número de circuitos. Note que hay tubos comunes en circuitos lo
que implica que el caudal por ellos debe estar corregido para el siguiente circuito en la iteración.

7. Los pasos 4,5 y 6 se repiten hasta que las ecuaciones de cabeza (ecuación 69) en los nodos que
conforman el circuito se cumplan, es decir, cuando su valor sea muy cercano a 0.

6.3 Método de Hardy-Cross con corrección de cabezas

Este método se utiliza para resolver las ecuaciones de cabeza (ecuación 69) para los nodos de la red. Las
ecuaciones del método son las ecuaciones 67. Por lo tanto en lugar de suponer caudales en las tuberı́as, se
supone la cabeza en cada uno de los (N u − 1) nodos de la red. Note que la cabeza en al menos uno de los
nodos debe ser conocida. Luego, se ajustan las cabezas en cada uno de los nodos y el proceso de ajuste se
repite hasta que la ecuación de continuidad en los nodos (ecuación 64) se cumpla.
Si se considera que la cabeza Hi en un nodo i esta sobrestimada o subestimada, el factor de ajuste de
cabeza en un nodo i es ∆Hi , la ecuación 67 se convierte en:
 1
2
(Hj − Hi ) − ∆Hi  p
Qij =  P L
2gAij (78)
Kmij + fij Dijij

Analizando el termino ((Hj − Hi ) − ∆Hi )1/2 y utilizando el teorema del binomio, la ecuación 78 se
expresa como:
√  
2gAij 1/2 1 −1/2
Qij =   1 (Hj − Hi ) − 2 (Hj − Hi ) ∆Hi (79)
P Lij 2
Kmij + fij Dij

25
 Teniendo en cuenta la ecuación 67 y reemplazando la ecuación 79 en la ecuación de continuidad para
un nodo i (ecuación 64), se tiene que:
N Ti NT
X ∆Hi Xi Qij
Qij − QDi − =0 (80)
2 Hj − Hi
j=1 j=1

Despejando el factor de corrección de cabezas (∆Hi ) para el nodo i, se tiene que:

P
2( Qij − QDi )
∆Hi = P Qij (81)
Hj −Hi

Note que en un nodo i, en la sumatoria de Qij se debe incluir la entrada (posible) de caudal externo
Qe . Cuando se asumen las primeras cabezas en todo los nodos se calcula Qij (utilizando la ecuación 67)
y a su vez se corrige la cabeza de tal forma que los próximos caudales se calculan con cabezas corregidas.
La cabeza para la iteración siguiente en el nodo i se calcula como:

Hit+1 = Hit + ∆Hi (82)

Si en alguna de las tuberı́as de la red se tiene una bomba rotodinamica, la cabeza de la bomba expresada,
por ejemplo, de la siguiente manera:

HB = AQ2ij + BQij + C (83)

afecta el calculo del caudal. Teniendo en cuenta esta adición de energı́a HB , la ecuación 66, queda
expresada como:

Q2ij
X 
Lij
HB = AQ2ij + BQij + CHj − Hi = Kmij + fij − (AQ2ij + BQij + C) (84)
2gA2ij Dij
Teniendo en cuenta que la ecuación 84 tiene la forma de una polinómio de grado 2, aQ2ij + bQij + c = 0,
la solución es encontrar las raı́ces reales del polinómio.

Comprobación de diseño y diseño

Pasos para la comprobación de diseño y para el diseño usando el método de Hardy-Cross con correción de
cabezas.
1. Definición de la topologı́a de la red. Los circuitos están conformados por tuberı́as y la conección
de las tuberı́as en los circuitos es positiva siguiendo las manecillas del reloj y negativo en sentido
contrario. Todos los tubos deben quedar incluidos en al menos un circuito. Definir los nodos con
cabeza conocida.
2. Suponer la cabeza de energı́a en cada uno de los nodos de la red. Note que al menos un nodo es de
cabeza conocida y constante en el tiempo. Entre mejor la estimación de cabezas iniciales, más rápido
convergerá el método.
3. Calcular el caudal en cada una de las tuberı́as utilizando la ecuación 67.
4. Utilizando la ecuación 64, comprobar la conservación de la masa en el nodo i. Si esta ecuación no se
cumple (> η = 1x10−6 ), se debe corregir la cabeza del nodo utilizando la ecuación 81.
5. Se repiten los pasos 3 y 4 para cada nodo.
6. El algoritmo se detendrá una vez la ecuación 64 se cumpla (con un margen de error muy pequeño)
en todos los nodos de la red.

26
6.4 Comparación de los métodos de Hardy-Cross
ˆ El método de H-C con corrección de caudale converge mas rápido ya que la corrección de caudales
en todas las tuberı́as de circuito se hace de manera simultánea en cada iteración, a diferencia del
método de correción de cabezas en donde la correción se hace nodo por nodo.

ˆ El método de H-C con correción de cabezas tiene la ventaja de que inicializar la energı́a en los nodos
es más sensillo, por lo tanto el paso previo al proceso iterativo es más rápido. Esto contrasta con la
suposición de caudales iniciales en la red, lo cual es más lento ya que es necesario que se cumpla la
conservación de la masa en cada nodo.

ˆ El diseño de redes de tuberı́as cerradas se hace mediante prueba y error. Esto quiere decir que se
asumen diferentes combinaciones de diámetros en la red hasta llegar al diseño más adecuado (óptimo)
para la cual hay que correr la red tantas veces como combinaciones de diámetros en la red tengamos.
Es por esto que el método de H-C con correción de cabezas, teniendo en cuenta la rapidez para
calcular las cabezas iniciales en los nodos, es más apropiado en general.

6.5 Método de análisis lineal

Este método fue desarrollado por D.J Wood y C.O. Charles entre 1970 y 1972. Se basa en la linearización
de las ecuaciones de conservación de la energı́a de cada una de las tuberı́as de la red. Este conjunto de
ecuaciones linealizadas junto con las ecuaciones lineales de coservación de la masa en los nodos de la red,
forman un sistema de ecuaciones lineales que se resuelve a traves del algebra lineal. El problema consiste
en formar la matriz que contiene los términos de las ecuaciones e invertirla para solucionar el sistema de
ecuaciones que son los caudales por cada tuberı́a. Se ha demostrado que este método converge más rápido
que los métodos ya explicados.

Ecuaciones usadas en el método

1. Conservación de la masa en el nodo: Para cada nodo de la red, se debe cumplir la ecuación de
continuidad:

N Ti
X
Qij − QDi = 0 (85)
j=1

Note que los caudales entrantes a la red (Qei ) en un nodo i deben tenerse encuenta en la ecuación 85.
Recuerde que los caudales que entran al nodo son positivos y que los que salen son negativos. Si N u
es el número de nodos habrá N u ecuaciones 85 y una redundande.

2. Conservación de la energı́a en el circuito: Para cada uno de los circuitos de la red se debe
cumplir la ecuación de la energı́a:

N Ti
X
hfij + hmij = 0 (86)
j=1

donde N Ti es el numero de tubos conectados en el circuito i. Si el numero de circuitos es N C, se

tendrán N C ecuaciones 86. La ecuación 86 se puede expresar como:

N Ti X 
X Qij Lij
2 Kmij + fij =0 (87)
2gA ij Dij
j=1

27
Note que la ecuación 87 es una ecuación no lineal. Para la solución de las N C ecuaciones no lineales, las
ecuaciones 87, se puede expresar como:
N Ti
X
kij Q2ij = 0 (88)
j=1
L
Kmij +fij Dij
ij
donde kij = 2gA2ij
.
Para resolver el sistema de ecuaciones 85 y 88, el método de teoria lineal propone que para un tramo
ij:

′
hfij + hmij = kij Qij (89)
donde k ′ ij = kij Q0ij . El caudal Q0ij es el caudal estimado para la primera iteración, o el caudal
corregido en la iteración previa para las demás iteraciones. Al reemplazar la ecuación 89 en la ecuación 88,
se tiene:

N Ti
X
′
kij Q2ij = 0 (90)
j=1

En el caso de existir una bomba en una de las tuberı́as del circuito, la cual sumistra una cabeza de
energı́a HB , la ecuación 90 queda:

N Ti
X
′
kij Q2ij = HB (91)
j=1

El numero total de ecuaciones a resolver es N T = (N U − 1) + N C; N U − 1 ecuaciones 85 y N C

ecuaciones 90. En N T coincide con el número de tuberı́as en la red y lo que se desea es encontrar Qij para
cada tuberı́a.
Para utilizar el método, es necesario estimar caudales iniciales en las tuberı́as (Q0ij ). Una de las
ventajas del método es que al suponer estos caudales iniciales no es necesario que se conserve la masa en
los nodos; los caudales pueden ser los mismos en todas las tuberı́as. Este método propone que el caudal
en la iteración siguiente es igual a:

Qt0ij + Qtij
Qt+1
0ij = (92)
2

Las N T ecuaciones se pueden expresar matricialmente, y el metodo consiste en invertir la matrix

N T xN T para determinar los caudales Qij . Si se tiene la red de la figura 14, se tiene las siguientes
ecuaciones que forman el sistema de ecuaciones:

ˆ Ecuaciones de continuidad en los nodos: Las siguientes ecuaciones se deducen teniendo en

cuenta que los caudales entrantes al nodo son positivos y los caudales salientes son negativos.
−Q12 + Q16 = −QE
Q12 − Q23 − Q25 = QD2
Q23 − Q34 = QD3
(93)
Q34 − Q45 = QD4
Q25 + Q45 − Q56 = QD5
Q56 − Q61 = QD6 Redundante

28
 Figure 14: Red cerrada que ilustra el uso del método de programación lineal (tomado de [1]).

ˆ Ecuaciones de conservación de la energı́a en los circuitos: Las siguientes ecuaciones se

deducen teniendo en cuenta que las pérdidas de energı́a en el circuito que siguen las manecillas del
reloj son positivas mientras que las pérdidas en sentido contrario son negativas.
′ ′ ′ ′
k12 Q12 + k25 Q25 + k56 Q56 + k16 Q16 = 0
′ ′ ′ ′ (94)
k23 Q23 + k34 Q34 + k45 Q45 − k25 Q25 = 0

Ordenando las ecuaciones 93 y 94 em forma matricial se tiene:

−1 0 −QE
    
0 0 0 0 1 Q12
 1 −1 0 0 −1 0 0  Q23   QD2 
   

 0
 1 −1 0 0 0 0  Q34   QD3 
   

 0
 0 1 −1 0 0 0  Q45  =  QD4 
    (95)
 0
 0 0 1 1 −1 0  Q25   QD5 
   

k ′ 0 0 0 k ′ k ′ ′  Q 
k61 56
 0 
12 25 56
0 k23′ ′
k34 ′
k45 −k25 ′ 0 0 Q61 0

La ecuación 95 en la forma compacta, se expresa como:

A·B=C (96)
La incognita es el vector B, el cual contiene los caudales de cada tuberı́a. La solución del sistema es:

B = A−1 · C (97)
Note que los valores de k ′ ij se calculan en función de Q0ij para la primera iteración. Para las iteraciones
siguientes Q0ij es calculado con la ecuación 92.
El algoritmo para solucionar la ecuación 96 utilizando el método de programación lineal es el siguiente:

1. Se suponen todos los caudales en cada una de las tuberı́as asi como la dirección del caudal (siguiendo
las manecillas del reloj o en sentido contrario). No es necesario cumplir con ecuación de continuidad
en los nodos.
′ .
2. Con estos caudales se calculan los valores de kij

3. Se plantean las ecuaciones lineales de continuidad (ecuación 85) y de conservación de energı́a (ecuación 90).

4. Se construye la matriz A de la ecuación 96.

29
 5. Se encuentran los caudales Qij al solucionar la ecuación 97.

6. Utilizando la ecuación 92, se corrigen los caudales Q0ij .

′ para los caudales corregidos.
7. Se calculan de nuevo los valores de kij

8. Se repiten los pasos 3 al 7 hasta que los caudales Qij cambien muy poco entre dos iteraciones sucesivas.

References
[1] Juan Guillermo Saldarriaga Valderrama. Hidráulica de tuberı́as. abastecimiento de agua, redes, riegos.
Editorial Alfaomega, Ediciones Uniandes, 2007.

30