1.

INTRODUCCION

El método de Hardy Cross es un método iterativo el cual permite estimar caudales

en arreglos de tuberías en los cuales es complicado definir si se encuentran en serie
o paralelo.

Para el desarrollo de éste método se siguen los siguientes pasos:

1. Dibujar el diagrama de la red de tuberías y definir los nodos.

Figura 1. Elementos del diagrama de tuberías.

Nota: Al definir los nodos de la red de tuberías, se deben establecer los

anillos y las secciones correspondientes a cada uno de ellos.

2. Calcular las resistencias hidráulicas en cada sección.

3. Suponer los caudales de acuerdo a las resistencias hidráulicas calculados.

Nota: Tener en cuenta que en la sección donde se tiene mayor resistencia

hidráulica pasará un menor caudal.

4. Organizar los datos en una tabla que contenga las secciones de cada anillo.

5. Reemplazar en la tabla los valores de resistencia hidráulica calculados, asumiendo

que la resistencia hidráulica será constante debido a que el flujo es turbulento y
totalmente desarrollado.
6. Teniendo en cuenta los caudales supuestos y la convención definida, alimentar la
tabla con los valores de los caudales
7. Calcular 2|𝑅𝐻∀|.

8. Calcular el producto de 𝑅𝐻∀2, considerando que el signo de la expresión será el mismo

que tenga el caudal ∀.

9. Calcular la corrección de caudal (∆∀) de acuerdo a la convención. La expresión

para esto es:
∑𝑅𝐻∀2
∆∀= −
∑ 2|𝑅𝐻 ∀|

10. Para iniciar con las iteraciones, se calculan los nuevos caudales de cada sección de
acuerdo a la expresión:

∀𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜= ±(±∀𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟 ± ∆∀𝑖 ± ∆∀𝑗)

En donde:

∀𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟 : Caudal asumido en el paso3.

∆∀𝑖 : Corrección del anillo 𝑖.
∆∀𝑗 : Corrección del anillo 𝑗.

Nota: los signos ± son definidos de acuerdo a la convención.

El procedimiento se realiza hasta que el ∆∀ se encuentre en el rango de error

deseado.

Es importante aclarar que al finalizar el procedimiento se debe verificar que el flujo

permanezca turbulento, es decir, que al calcular el número de Reynolds y la relación
𝑒 ⁄𝑑 se encuentre en la zona de flujo turbulento para garantizar que las
suposiciones de los factores de fricción son las correctas. De lo contrario, se
procede a realizar
de nuevo todo el procedimiento desde el paso 1, para una segunda gran iteración,
asumiendo unos nuevos valores para el factor de fricción en la zona que turbulento.
 Figura 2. Convención de signos.

Teniendo en cuenta el diagrama anterior si la dirección del caudal va en el mismo

sentido que el establecido (verde) el signo será positivo. Si la dirección del caudal
va en sentido contrario será negativo.

Para la primera corrección de caudal el signo dependerá del sentido del anillo en
estudio, por ejemplo, la sección 8-3 en relación con el anillo I el caudal va en el
mismo sentido (amarillo), pero con relación al anillo II el caudal se desplaza en
sentido contrario (negro).

Para obtener una mayor precisión en los cálculos empleados al realizar las
iteraciones, se pueden emplear las formulas establecidas por distintos autores. Al
suponer un factor de fricción para un flujo totalmente turbulento y desarrollado se
puede emplear la segunda ecuación de Karman-Prandtl:

1 𝐷
= 2 log10 ( + 1,74)

√𝑓 2𝑒

Para un flujo en la zona de transición se pude emplear la ecuación de Colebrook-

White:

1 2,51
𝑒⁄
𝐷
= −2 log10 ( + )
√𝜆 3,7
𝑅𝑒√𝜆
 2. DESARROLLO DE LOS PROBLEMAS

2.1 PROBLEMA 1:

Por la red de tuberías circula agua 𝑄 = 1500 𝐿/𝑚𝑖𝑛. Las tuberías de 600 y 400 𝑚𝑚 son de
cemento alisado (𝑒 = 0.5 𝑚𝑚) y las tuberías de 300 y 450 𝑚𝑚 de fundición (𝑒 = 1.2 𝑚𝑚). La
presión en 𝐴 es 5𝑏𝑎𝑟. Se tiene que 𝐷1 = 300 𝑚𝑚,𝐷2=400𝑚𝑚,
𝐷3 = 450 𝑚𝑚 y 𝐷4 = 600 𝑚𝑚.

Figura 3. Diagrama del problema 1.

a) Calcular los caudales que circulan por las diferentes ramas y la presión en B. Cualquier
suposición debe ser explicada y verificada al final de proceso de iteración.

b) Implementar una simulación en EPANET u otra herramienta para el estudio de redes

de tuberías. Los resultados obtenidos compárelos con los del punto A y explique las
diferencias presentadas.
Solución:

Figura 4. Distribución de caudales del problema 1

En primer lugar, se calculan los valores correspondientes para 𝑒/𝐷, el

área 𝐴 y la resistencia hidráulica 𝑅. También se estima el valor del
factor de fricción, para lo cual se asume que el flujo es turbulento y
completamente desarrollado; se emplea la ecuación de Karman-
Prandtl. Estos resultados se encuentran en la Tabla 1.
 SECCION Dint (mm) e (mm) e/D f (TURB) A (m2) L (m) R (s2/m5)
1-2 428.70 1.200 0.0028 0.0257 0.144343373 350.0 51.2
2-3 428.70 1.200 0.0028 0.0257 0.144343373 175.0 25.6
3-4 428.70 1.200 0.0028 0.0257 0.144343373 175.0 25.6
4-5 381.00 0.500 0.0013 0.0210 0.114009183 250.0 54.0
5-6 428.70 1.200 0.0028 0.0257 0.144343373 250.0 36.6
6-7 574.60 0.500 0.0009 0.0190 0.25931111 350.0 8.8
7-8 574.60 0.500 0.0009 0.0190 0.25931111 350.0 8.8
7-9 428.70 1.200 0.0028 0.0257 0.144343373 250.0 36.6
9-10 303.30 1.200 0.0040 0.0283 0.072249476 175.0 159.5
5-10 303.30 1.200 0.0040 0.0283 0.072249476 175.0 159.5
2-9 381.00 0.500 0.0013 0.0210 0.114009183 250.0 54.0
3-10 381.00 0.500 0.0013 0.0210 0.114009183 250.0 54.0
1-8 574.60 0.500 0.0009 0.0190 0.25931111 500.0 12.5

1 ITERACIÓN TIENE EL SIGNO DE Q

ANILLO SECCIÓN R Q ABS(2R*Q) R*Q2 ∆Q
1-2 51.23 3 307.4 461.0570329
2-9 54.01 3 324.1 486.1248781
7-9 36.59 1 73.2 36.59182801
I
7-8 8.76 -4 70.1 -140.2037789
1-8 12.52 7 175.3 -613.3915327
SUMA 950.0 230.1784274 -0.2423
2-3 25.61 2 102.5 -102.4571184
3-10 54.01 1 108.0 54.01387535
II 9-10 159.49 -2 638.0 637.9676437
2-9 54.01 -1 108.0 54.01387535
SUMA 956.5 643.5382759 -0.6728
3-4 25.61 1 51.2 25.6142796
4-5 54.01 -1 108.0 54.01387535
III 5-10 159.49 2 638.0 637.9676437
3-10 54.01 1 108.0 -54.01387535
SUMA 905.3 663.5819233 -0.7330
9-10 159.49 2 638.0 637.9676437
5-10 159.49 -2 638.0 637.9676437
5-6 36.59 -3 219.6 329.3264521
IV
6-7 8.76 -3 52.6 -78.86462564
7-9 36.59 4 292.7 -585.4692481
SUMA 1840.8 940.9278657 -0.5112
 23 ITERACIÓN TIENE EL SIGNO DE Q
ANILLO SECCIÓN R Q ABS(2R*Q) R*Q2 ∆Q
1-2 51.23 2.744789825 281.2 385.949366
2-9 54.01 2.744789825 296.5 406.933579
7-9 36.59 -2.273573438 166.4 189.1481419
I
7-8 8.76 -6.844975366 120.0 -410.566505
1-8 12.52 6.744789825 168.9 -569.480082
1033.0 1.984499883 -0.0019
2-3 25.61 5.018363263 257.1 -645.069245
3-10 54.01 0.561682279 60.7 17.04067458
II 9-10 159.49 -1.571401928 501.3 -393.8340169
2-9 54.01 -4.273573438 461.7 986.4786275
1280.7 -35.38395972 0.0276
3-4 25.61 4.456680983 228.3 508.7509595
4-5 54.01 2.456680983 265.4 -325.9889401
III 5-10 159.49 1.133084208 361.4 -204.7684462
3-10 54.01 0.561682279 60.7 -17.04067458
915.8 -39.04710142 0.0426
9-10 159.49 2.428598072 774.7 940.6974205
5-10 159.49 -2.428598072 774.7 -940.6974205
5-6 36.59 -0.410234809 30.0 -6.158133815
IV
6-7 8.76 -0.410234809 7.2 -1.474703641
7-9 36.59 1.155024634 84.5 -48.81649561
1671.1 -56.44933306 0.0338

Qreal Qreal Qreal Qreal Qreal

4.1172 0.07068583 58.246 0.3 2.18E+07
4.1172 0.03141593 131.054 0.2 3.28E+07
3.4104 0.03141593 108.555 0.2 2.71E+07
10.2675 0.03141593 326.823 0.2 8.17E+07
10.1172 0.04908739 206.106 0.25 6.44E+07

7.5275 0.03141593 239.609 0.2 5.99E+07

0.8425 0.01767146 47.677 0.15 8.94E+06
2.3571 0.04908739 48.019 0.25 1.50E+07
6.4104 0.01767146 362.752 0.15 6.80E+07

6.6850 0.03141593 212.791 0.2 5.32E+07

3.6850 0.01767146 208.530 0.15 3.91E+07
1.6996 0.04908739 34.625 0.25 1.08E+07
0.8425 0.04908739 17.164 0.25 5.36E+06
e (mm) 0.15
f (TOTAL
Sección Dint (mm) e/D TURB) A (m2) L (m) R (s2/m5)
A-B 381 0.0004 0.0158 0.114009 457.21 74.50
B-C 381 0.0004 0.0158 0.114009 457.21 74.50
C-F 381 0.0004 0.0158 0.114009 609.61 99.33
F-I 202.7 0.0007 0.0183 0.032270 1219.21 5376.90
I-H 202.7 0.0007 0.0183 0.032270 457.21 2016.34
H-G 303.3 0.0005 0.0167 0.072249 457.21 245.09
G-D 333.3 0.0005 0.0163 0.087249 1219.21 399.34
D-E 381 0.0004 0.0158 0.114009 457.21 74.50
E-F 303.3 0.0005 0.0167 0.072249 457.21 245.09
A-D 428.7 0.0003 0.0154 0.144343 609.61 53.68
B-E 303.3 0.0005 0.0167 0.072249 609.61 326.78
E-H 303.3 0.0005 0.0167 0.072249 1219.21 653.57

1 ITERACIÓN TIENE EL SIGNO DE Q

ANILLO SECCIÓN 𝑅 𝑚 Q ABS(2R*Q) 𝑅 𝑄 ∆Q
A-B 74.50 0.21237435 31.6 3.359999
B-E 326.78 0.07079145 46.3 1.63765869
I D-E 74.50 -0.1557412 23.2 -1.8069328
A-D 53.68 -0.2265326 24.3 -2.754653
SUMA 125.4 0.43607193 -0.0035
B-C 74.50 0.1415829 21.1 1.49333289
C-F 99.33 0.09910803 19.7 0.97564416
II E-F 245.09 -0.1132663 55.5 -3.1443047
B-E 326.78 -0.0707914 46.3 -1.6376587
SUMA 142.6 -2.3129863 0.0162
H-G 245.09 0.09910803 48.6 2.40735827
G-D 399.34 0.01415829 11.3 0.08005059
III D-E 74.50 -0.0849497 12.7 -0.5375998
E-H 653.57 0.11326632 148.1 8.38481249
SUMA 220.6 10.3346215 -0.0468
F-I 5376.90 0.01415829 152.3 1.07783845
I-H 2016.34 -0.0707914 285.5 -10.104736
IV E-H 653.57 0.15574119 203.6 15.8525361
E-F 245.09 0.08494974 41.6 1.76867138
SUMA 683.0 8.59431044 -0.0126
 2 ITERACIÓN TIENE EL SIGNO DE Q
ANILLO SECCIÓN 𝑅 𝑚 Q ABS(2R*Q) 𝑅 𝑄 ∆Q
A-B 74.50 0.20889783 31.1 3.25089483
B-E 326.78 0.05109154 33.4 0.85302111
I D-E 74.50 -0.1123699 16.7 -0.9406658
A-D 53.68 -0.2300091 24.7 -2.839851
106.0 0.32339914 -0.0031
B-C 74.50 0.1578063 23.5 1.85516991
C-F 99.33 0.11533143 22.9 1.32120164
II E-F 245.09 -0.0845 41.4 -1.7482886
B-E 326.78 -0.0904914 59.1 -2.6759377
147.0 -1.2478548 0.0085
F-I 245.09 0.05573679 27.3 0.76138952
I-H 399.34 -0.0201054 16.1 -0.1614236
III E-H 74.50 -0.1317975 19.6 -1.2940448
E-F 653.57 0.06641856 86.8 2.88317343
149.8 2.18909451 -0.0146
H-G 5376.90 -0.0146492 157.5 -1.1538785
G-D 2016.34 -0.0833755 336.2 -14.016536
IV D-E 653.57 0.14315709 187.1 13.3942231
E-H 245.09 0.05068608 24.8 0.62965183
705.7 -1.1465397 0.0016
 11 ITERACIÓN TIENE EL SIGNO DE Q
ANILLO SECCIÓN 𝑅 𝑚 Q ABS(2R*Q) 𝑅 𝑄 ∆Q
A-B 74.50 0.2138826 31.9 3.4078929
B-E 326.78 -4.417E-05 0.0 6.3747E-07
I D-E 𝑅 𝑚74.50 -0.0973636 14.5 -0.7062014
A-D 53.68 -0.2250244 24.2 -2.7180942
70.6 -0.0164021 0.0002
B-C 74.50 0.21392676 31.9 3.40930052
C-F 99.33 0.17145189 34.1 2.91983531
II E-F 245.09 -0.0505 24.7 -0.6239216
B-E 326.78 -0.1416271 92.6 -6.5547224
183.2 -0.8495081 0.0046
F-I 245.09 0.04073048 20.0 0.40659501
I-H 399.34 -0.0522435 41.7 -1.0899512
III E-H 74.50 -0.141819 21.1 -1.4983184
E-F 653.57 0.05639702 73.7 2.07875864
156.5 -0.102916 0.0007
H-G 5376.90 -0.0486531 523.2 -12.727801
G-D 2016.34 -0.061259 247.0 -7.5666378
IV D-E 653.57 0.16527365 216.0 17.8524995
E-H 245.09 0.01854797 9.1 0.08431717
995.4 -2.3576217 0.0024
Q 0.4389
Qreal A (m2) V (m/s) Dint (m) Re
0.0939 0.1140 0.8234 0.3810 V*D/(0,8*10^-6)
0.0000 0.0722 0.0003 0.3033 1.02E+02
0.0427 0.1140 0.3748 0.3810 1.79E+05
0.0988 0.1443 0.6842 0.4287 3.67E+05

0.0939 0.1140 0.8236 0.3810 3.92E+05

0.0753 0.1140 0.6600 0.3810 3.14E+05
0.0221 0.0722 0.3065 0.3033 1.16E+05
0.0622 0.0722 0.8604 0.3033 3.26E+05

0.0179 0.0323 0.5540 0.2027 1.40E+05

0.0229 0.0323 0.7106 0.2027 1.80E+05
0.0622 0.0722 0.8615 0.3033 3.27E+05
0.0248 0.0722 0.3426 0.3033 1.30E+05

0.0214 0.0722 0.2956 0.3033 1.12.E+05

0.0269 0.0872 0.3082 0.3333 1.28.E+05
0.0725 0.1140 0.6363 0.3810 3.03.E+05
0.0081 0.0722 0.1127 0.3033 4.27E+04

Hr A-H
Caudal (pie3/s) 15.5 CIRCUITO A, D, G, H 10.10036 m
PA (psig) 700 CIRCUITO A, B, E, H 2.15337 m
CIRCUITO A, D, E, H 2.16057 m
PH 4.8.E+06 Pa
PH 693.1652 psig