En Hidráulica se entiende por tubería cualquier conducto

cerrado que transporta agua a presión. Por lo general son de

sección circular. Los conductos cerrados en que el agua circula
sin llenarlos completamente se clasifican en su totalidad como
canales abiertos, como las alcantarillas y tuberías de
avenamiento. Como los factores de que dependen las perdidas
de carga en los conductos son independientes de la presión, se
aplica la misma ley a la corriente en tuberías y en canales
abiertos, y las formulas para cada caso tienen la misma forma.
Se aplicara el principio de la energía a la solución de
problemas prácticos de flujos en tuberías, que frecuentemente
se presentan en las diversas ramas de la ingeniería. El flujo de
un fluido real es mucho más complejo que el de un fluido
ideal. Debido a la viscosidad de los fluidos reales, en su
movimiento aparecen fuerzas cortantes entre las partículas
fluidas y las paredes del contorno y entre las diferentes capas
de fluido. Las ecuaciones diferenciales en derivadas parciales,
que resolverían de forma general el problema del flujo
(Ecuaciones de Euler CAPIII), no admiten, por lo común, una
solución. Como consecuencia, los problemas de flujos reales se
resuelven aprovechando datos experimentales y usando
métodos semi-empiricos. Existen dos tipos de flujos

226
 permanentes en el caso de flujos reales, que es necesario
considerar y entender. Estos son el laminar y turbulento,
ambos están gobernados por leyes distintas, como se hace
referencia en el capitulo anterior en forma general.

V.1 OBJETIVOS DEL CAPITULO

Dentro del desarrollo del capitulo presente, se darán a conocer

las distintas situaciones de flujo permanente y las diferentes
consideraciones; El alumno será capaz de:
i. Conocer la resistencia al flujo en las
tuberías (perdidas por: rozamiento,
ensanchamiento, contracciones, válvulas,
etc.).
ii. Corrientes de flujo en tuberías laminar y
turbulenta.
iii. Efectos del fluido dentro de las tuberías
“golpe de Ariete”.
iv. Análisis de tuberías simples y múltiples.
v. Análisis de tuberías en un sistema de redes
y métodos de resolución,”Método de
Cross, Método de Newton-Raphson”.

227
 V.2 PRINCIPIOS FUNDAMENTALES DE FLUJO EN
TUBERIAS

Supóngase que la figura 5-1 representa la sección longitudinal

de una parte de tubería, si imaginamos una partícula cualquiera
P de liquido, puede suceder dos cosas respecto al movimiento
de esa partícula, puede que tenga un movimiento dirigido en
igual dirección que el eje del tubo es decir paralelo(Fig5-1a) o
que tenga un movimiento en dirección indefinida (Fig5-1b).
Cuando todas las partículas se van moviendo
longitudinalmente sin acercarse ni alejarse del eje del tubo, van
formando filetes líquidos con una determinada velocidad y
entonces el liquido circula en lo que se llama régimen
tranquilo.

228
 Puede suceder que las partículas tengan moviemientos
indefinidos, en este caso el liquido ya no circula formando
filetes paralelos, sino trayectorias mas o menos irregulares y el
régimen es llamado turbulento.
Si hacemos un corte transversal A-A (Fig5-2), veremos que los
filetes líquidos que tienen igual velocidad están distribuidos en
cilindros concéntricos de tal modo que los mas rápidos serán
los del centro y los mas lentos los cercanos a las paredes del
tubo, ello es debido al efecto del rozamiento. Como el gasto
depende de la velocidad de las partículas, en una sección dada,
para gastos grandes el régimen es turbulento y para gatos
pequeños se tendrá régimen laminar.

Imaginémonos un tubo por el que va a circular agua con

velocidades que podemos hacer variar. Al principio vamos a
hacer que el líquido contenido en el tubo, y que esta en reposo,
circule pasando por una sección dada, con una velocidad

229
 progresivamente mayor. cuando el líquido comienza a
moverse, como la velocidad es pequeña, el régimen es laminar;
al ir aumentando paulatinamente la velocidad llega a cierto
valor a partir del cual el régimen se vuelve turbulento. Este
valor que tiene que alcanzar la velocidad para pasar del
régimen laminar al turbulento, se llama Velocidad Critica
Alta.
Supongamos ahora que en la tubería esta circulando el agua a
gran velocidad, de tal manera que el régimen es turbulento y
que le hace disminuir paulatinamente; el régimen permanece
turbulento hasta que llega un momento que si se disminuye
mas la velocidad pasa a laminar. A esta velocidad para la cual
pasa de turbulento a laminar, es llamada Velocidad Crítica
Baja.

V.2.1 CORRIENTE LAMINAR EN TUBERIAS-ECUACION

DE DARCY-WEISBACH

Las leyes de la corriente laminar en las tuberías rectas se

determinaron experimentalmente, en forma independiente por
Hagen y Poiseville. La ley de Hagen-Poiseville puede
deducirse a partir de los principios fundamentales,
considérense las fuerzas que obran sobre un cilindro de

230
 Líquido de longitud “l” y radio “y”, como se ilustra en la
figura 5-3. Si el movimiento es uniforme, la fuerza originada
por la diferencia de presión sobre los extremos del cilindro,
tiene que equilibrarse exactamente con la resultante de los
esfuerzos cortantes en su superficie lateral, como lo expresa la
ecuación (5.2.1):

( p1  p 2 )y 2   2yl 5.2.1

Simplificando e introduciendo el valor de τ:
du
( p1  p 2 ) y  2 l
dy 5.2.2
Donde:
u; es la dirección según figura 5-3
μ; viscosidad (Kg-s/m2)
El segundo miembro es negativo a causa de que du y dy son de
signos contrarios (Fig5-3),. despejando la ecuación 5.2.2
tendremos:
( p1  p 2 ) y 2
u   du 
4ul
C [Kg-s/m5.2.3
2
]

231
El valor de C se calcula por la condición límite de que y =d/2
cuando u=0, tendremos:

( p1  p 2 )d 2
C
16ul
5.2.4

( p1  p 2 )  d2 
u
4ul 
  y 2 
4


5.2.5

La ecuación 5.2.5 indica que la distribución de la velocidad

para la corriente laminar en tuberías circulares es un
paraboloide de revolución cuyo eje es el del tubo. El valor
máximo de “u” se determina haciendo y= 0,
( p1  p 2 )d 2
u max 
16ul
5.2.6
El gasto se obtiene por integración:
Q  dQ   udA   u 2ydy 5.2.7

Sustituyendo en la ecuación 5.2.5:

 ( p1  p 2 ) d / 2  d 2 ydy 
Q 0    
2
y dy
2l 4 

 ( p1  p2 ) d 4
Q
128l [m3/s] 5.2.8

232
 La velocidad media, “v”, se obtiene de la ecuación 5.2.8
dividiendo por el área del tubo:

 ( p1  p 2 )d 4
v
32 l [m/s] 5.2.9

Puede verse por las ecuaciones 5.2.7 y 5.2.8 que en corriente

laminar, la velocidad media es la mitad de la velocidad
máxima. La expresión para la variación de la carga
piezomaetrica en una longitud “l” se la obtendrá de la siguiente
manera:

( p1  p 2 ) 32lv
h 2 5.2.10
w d g

Introduciendo el número de Reynolds (Capitulo IV), la

expresión se transforma en:

64 l v 2
hf 
Re d 2 g 5.2.11

Si en la ecuación 5.2.11 se usa el coeficiente de rozamiento en

lugar de 64/R, tendremos:

233
 l v2
hf  f
d 2g 5.2.12
Ec. Darcy-
Weisbach
Donde:
f; es el factor de rozamiento
v; es la velocidad
d; diámetro de la tubería
Re; número de Reynolds
hf; perdida por fricción

Conocida como la formula de Darcy-Weisbach para la

corriente en tuberías, siendo h la perdida de energía en l metros
de tubería en kilogramos por kilogramo. Lo anterior demuestra
que en corriente laminar, “f” es completamente independiente
de la rugosidad, pues solo varía con el número de Reynolds, es
decir, con la resistencia relativa de las fuerzas viscosa y de
inercia.

V.2.2 CORRIENTE TURBULENTA EN TUBERIAS

Cuando la corriente ocurre a números de Reynolds mayores de

4000, los valores de “f” en la ecuación 5.2.12 de Darcy-
Weisbach, varían con la rugosidad, la viscosidad y la densidad.
La corriente turbulenta se divide en tres grupos: en tubos lisos,
relativamente rugosos a velocidades grandes, y en la zona de
transición comprendida entre las dos primeras mencionadas.
234
 En la corriente de tubos muy lisos, los valores de “f” varían
con R (ANEXO). Se puede observar que nunca llega a
convertirse en una recta horizontal, lo cual demuestra que las
propiedades de los líquidos influyen sobre la corriente en todo
el intervalo de números de Reynolds. Los tubos de vidrio y los
de metal estirado con superficie muy lisa se encuentran en esta
categoría.
En la figura 5-4 se ilustra la corriente en tuberías rugosas a
valores grandes de R por medio de la zona situada arriba y a la
derecha de la línea de trazos, llamada zona de turbulencia
completamente establecida o plena. En esta zona, las curvas
“f” se vuelven horizontales, demostrando así que la corriente
es completamente independiente de las propiedades de los
líquidos. La tercera categoría de corriente turbulenta se
presenta cuando los valores de “f” se encuentran en la zona
comprendida entre la curva para tuberías lisas y la línea de
trazos (Fig5-4). La corriente en los tubos comerciales se
produce generalmente dentro de esta categoría. En esta zona,
las curvas de “f” para los diversos valores de la rugosidad
relativa se separan en puntos sucesivos de aquellas para
tuberías lisas y se vuelven horizontales cuando entran en la
zona de turbulencia plena.
Todas estas curvas presentadas en la figura 5-4 fueron
presentadas por Moody. Por lo que la figura 5-4 es conocida

235
como el Diagrama de Moody, con el cual conocido el número
de Reynolds y la rugosidad equivalente podemos encontrar el
correspondiente factor de rozamiento. Sin embargo si no se
conoce el número de Reynolds, teniendo de incógnitas el gasto
o el diámetro, el factor de rozamiento se obtiene por
aproximaciones sucesivas.

Fuente: Mecánica de Fluidos Aplicada-

Robert L. Mott
Coeficientes de rugosidad para algunos materiales

236
 Fuente: Mecánica de Fluidos Aplicada- DIAGRAMA DE
Robert L. Mott
MOODY

V.3 USO DEL DIAGRAMA DE MOODY

Como ya se menciono en el anterior subtitulo, el diagrama de

Moody se usa utiliza como una ayuda para determinar el valor
del factor de fricción, “f”, para flujo turbulento. Deben
conocerse los valores del número de Reynolds y de la
rugosidad relativa. Por consiguiente, los datos básicos
requeridos son el diámetro interior del conducto, el material

237
 con que el conducto esta hecho, la velocidad de flujo y el tipo
de fluido y su temperatura, con los cuales se puede encontrar la
viscosidad.

EJEMPLO DE APLICACIÓN Determine la fricción

o rozamiento “f” si el agua a 160ºF esta fluyendo a 30
pies/seg en un conducto de hierro forjado no
recubierto cuyo diámetro interior es de 1 pulg.
Solución. Primero, se debe calcular el número de
Reynolds para determinar si el flujo es laminar o
turbulento:
vD
NR 


Pero D=1 pulg. = 0.0833 pies, y υ = 4.38x10-6 pies2/s,

entonces tendremos:

(30)(0.0833)
NR   5.70  10 5
4.38  10 6

Por lo tanto el flujo es turbulento. Ahora se debe

calcular la rugosidad relativa. De la tabla 5.1 tenemos
Є= 8x10-4 pies. Entonces, la rugosidad relativa será:

D 0.0833
  1.04  10 2  104
 8  10  4

238
 Obsérvese que, con el fin de que D/Є sea un cociente
sin dimensiones, tanto D como Є deben estar en las
mismas unidades.
Por ultimo se debe localizar el número de Reynolds
sobre la abscisa del Diagrama de Moody (Fig5-4)
N R  5.70  10 5

proyectamos verticalmente hasta que se alcance la

curva correspondiente a D/Є=104, podemos usar la
curva numero 100 ya que es cercana. Proyectamos
horizontalmente hacia la izquierda, se lee el valor que
en este caso es aproximadamente f=0.038.

V.4 ECUACIONES DEL FACTOR DE FRICCION

El diagrama de Moody es un medio conveniente y lo

suficientemente preciso para determinar el factor de fricción
cuando se resuelven problemas mediante cálculos manuales.
Sin embargo, si los cálculos deben ser algo automáticos para
poder obtener la solución en una computadora o con una
calculadora programable, es necesario tener ecuaciones para el
factor de fricción.
La ecuación que se utiliza en el trabajo hecho por Moody
cubre tres diferentes zonas del diagrama de Moody. En la zona

239
 de flujo laminar, para valores del número de Reynolds por
debajo de 2000, “f” puede encontrarse con la ecuación:

f  64 / N R 5.4.1
Esta relación, esta graficada en el diagrama de Moody como
una línea recta en el lado izquierdo del diagrama. Desde luego,
para los números de Reynolds desde 2000 hasta 4000, el flujo
se encuentra en la región crítica y es imposible predecir el
valor de “f”.
Por encima del número de Reynolds de 4000, por lo general el
flujo se conoce como turbulento. Sin embargo, en esencia
existen dos zonas de interés en este punto. Hacia el lado
derecho del diagrama, el flujo esta en la zona de completa
turbulencia. Se puede observar que el valor de “f” no depende
del número de Reynolds, sino solo de rugosidad relativa. En
este intervalo, se aplica la ecuación:
1
 2 log10  3.7 D / 
f 5.4.2
La frontera de esta zona es la línea punteada que corre, por lo
general, de la parte superior izquierda a la parte inferior
derecha del diagrama, la ecuación de esta línea es:

1 NR

f 200( D / ) 5.4.3

240
La tercera zona del diagrama de Moody, que se conoce como
zona de transición, se encuentra entre la zona de completa
turbulencia y la línea que se identifica como conductos lisos.
La línea de “conductos lisos” tiene la ecuación:

1 N f 
f
 2 log10  R

 2.51



5.4.4

Siendo lisos, estos conductos no presentan irregularidades

superficiales al flujo, de modo que el factor de fricción solo es
función del número de reynolds. Los conductos hechos de
vidrio o de cobre tiene un valor de rugosidad relativa que los
acerca a la línea de conductos lisos.
En la zona de transición, el factor de fricción es función tanto
del número de Reynolds, como de la rugosidad relativa. C.F.
Colebrook desarrollo la relación para el factor de fricción en
esta zona:
1  1 2.51 
f
 2 log10 

 3.7 ( D / )

NR f



5.4.5

La ecuación siguiente, permite el cálculo directo del valor del

factor de fricción, esta fue desarrollada por P.K. Swamee y
A.K. Jain,

241
 0.25
f  2
  1 5.74  5.4.6
log  

  3.7( D / ) N R
0.9


La ecuación 5.4.6 produce valores para “f” que se encuentran

entre  1.0 % del valor de los correspondientes a la ecuación
5.4.5, dentro del intervalo de rugosidad relativa, D/Є,
comprendido entre 1000 y 1x106, y para números de Reynolds
que van de 5x103 hasta 1x108. Esta es virtualmente la zona de
turbulencia completadle diagrama de Moody.

EJEMPLO DE APLICACIÓN Calcule el valor del

factor de fricción si el numero de Reynolds para el
flujo es de 1x105, y la rugosidad relativa es de 2000.

Solución. Debido a que se encuentra en la zona

turbulenta, usamos la ecuación 5.4.6

0.25
f  2
  1 5.74 
log  
  3.7(2000) (1  10 )
5 0. 9


f  0.0204

V.5 PÉRDIDAS DE FRICCION EN SECCIONES

242
 TRANSVERSALES NO CIRCULARES

La ecuación de Darcy para la perdida por fricción puede ser

utilizada para secciones transversales no circulares si la
geometría esta representada por el radio hidráulico en lugar del
diámetro del conducto, como se utiliza en secciones circulares.
Después de calcular el radio hidráulico, podemos calcular el
número de Reynolds sustituyendo 4R en lugar del diámetro, D.
La ecuación, entonces, queda como:

 (4 R )   ( 4 R )
NR 


v 5.5.1

Una sustitución parecida se puede hacer para calcular la

perdida por fricción. En la ecuación de Darcy:
l v2
hf  f
4R 2 g 5.5.2

La rugosidad relativa, se convierte en 4R/Є. el factor de

fricción se puede encontrar en el diagrama de Moody.

EJEMPLO DE APLICACIÓN Determine la caída de

presión para una longitud de 50 m de un conducto
cuya sección transversal muestra la figura 5-5. Esta
fluyendo etilenglicol a 25ºC y con una rapidez de

243
 0.16 m3/s. La dimensión interna del cuadrado es de
250 mm y el diámetro exterior es de 1500 mm.
Utilice Є=3x10-5 m, un poco más liso que el conducto
de acero comercial.

Solución. El área, la velocidad, el radio hidráulico y

el número de Reynolds fueron calculados y se obtuvo
los siguientes resultados:
A  0.0448m 2

v  3.7 m / s

RH  0.0305m

N R  2.96 x10 4

El flujo es turbulento y se puede usar la ecuación de

Darcy para calcular la perdida de energía entre dos
puntos separados 50m. Con el fin de determinar el
244
 facto de fricción, primero debemos encontrar la
rugosidad relativa:
4 R /  ( 4)(0.0305) /(3 x10 5 )  4067

De acuerdo al diagrama de Moody, f  0.0245 entonces

tendremos:
l v2 50 (3.57) 2
hf  f  0.0245  
4R 2 g ( 4)(0.0305) ( 2)(9.81)

h f  6.52m

Si el conducto esta colocado horizontalmente,

hL  p /  p  hL 

En la que p es la caída de presión ocasionada por la

perdida de energía. entonces tendremos:
KN
p  10.79  6.52m p  70.4 Kpa
m3

V.6 RESISTENCIA AL FLUJO EN TUBERIAS

(ECUACIONES)

V.6.1 FORMULA DE HAZEN-WILLIAMS.

Esta formula, establecida tanto para tuberías como para canales

abiertos, ha sido muy usada para tuberías. La selección de los
exponentes se hizo con la idea de conseguir una variación
mínima del coeficiente para todos los conductos del mismo
grado de rugosidad.

245
La formula de Hazen-Williams puede ser aplicada a conductos
libres o conductos forzados, ha sido usada para tuberías de
agua y alcantarillados.
En el sistema técnico de unidades se expresa como:

v  1.538C HW Rh 0.63 S 0.54 5.6.1

Donde:
Rh; radio hidráulico (m)
S; pendiente de la línea de energía
CHW; coeficiente de Hazen-Williams, adimensional

Para tuberías funcionando a presión tenemos:

D hf
Rh 
4 ; S
l

Por lo que la ecuación de Hazen-William, se escribirá:

0.63 0.54
D  hf 
v  1.538C HW  
 4

 l


5.6.2
Despejando,
v1.85l
h f  2.271 1.85 5.6.3
C HW D1.16

Si multiplicamos y dividimos el segundo miembro por l/D,

tendremos:
v1.85 l
h f  2.271 1.85 5.6.4
C HW D1.16 D

246
de la misma manera, multiplicamos y dividimos el segundo
v2
miembro por 2g , tenemos:
 44.56  l v2  44.56  l v2
hf   1.85 1.16 0.15 
 0.15  5.6.5
 C HW D v  D 2 g  C HW D  Dv   D 2 g
1.85 0.01

multiplicando y dividiendo por v 0.15 ,

 44.56  l v2
hf   1.85 0.01 0.15 0.15  5.6.6
 C HW D Re v  D 2 g

Igualando esta ecuación a la ecuación de Darcy-Weisbach,

tendremos:

 44.56 
f  1.85 0.01 0.15 0.15  5.6.7
 C HW D Re v 

Para el CHW, existen distintos valores dependiendo de los

materiales, en la tabla 5.2 se muestran valores para algunos
tipos de materiales.

247
 Fuente: Introducción a la Hidráulica-
Texto auxiliar de la materia

V.6.2 FORMULA DE MANNING

Esta es una formula empírica, obtenida por experimentación, y

es otra alternativa de dar solución a problemas de flujos en
conductos cerrados y canales abiertos, esta es considerada para
tuberías de 1 metro de diámetro, siendo muy fiable para la
gama de diámetros comprendidos entre 0.40 y 1.30 m.
1
v
n
Rh 2 / 3 S 1 / 2 5.6.8
Donde:
Rh; radio hidráulico
S; pendiente de la línea de energía
n; coeficiente de Manning, adimensional

Para tuberías a presión, es decir a sección llena tendremos:

248
 2/3 1/ 2
1D  hf 
v  
n 4 

 l


5.6.9
Despejando:
n 2 lv 2
h f  6.34
D4/3
5.6.10
Multiplicando y dividiendo el segundo miembro por D y por
2g, tendremos:
 n2  l v2
h f  124.4 1 / 3 
 D  D 2g 5.6.11
igualando esta ecuación a la ecuación de Darcy-Weisbach,
tendremos:

 n2 
f  124.4 1 / 3 
 D  5.6.12
El coeficiente “n” depende del material de que esta construido
el tubo, el la tabla 5.3 se muestran algunos materiales con sus
respectivos valores de “n”.

Fuente: Abastecimiento de Agua

Potable-Ing. Juan Carlos Rocha C.

V.7 CONCEPTO DE LONGITUD VIRTUAL

249
 La longitud virtual es un concepto que se utiliza para
determinar el gasto que circulará por cada tramo de tubo, a este
se le denomina gasto propio.

De la figura se deduce que el tramo B-A abastece por un solo

lado, mientras que el tramo D-C que abastece por los dos lados
deberá conducir más gasto. Si se tratara de una zona con la
población uniformemente distribuida el tramo D-C conducirá
el doble de gasto que el tramo B-A. Relacionado gastos con
longitudes, parecería que el tramo D-C mide el doble de la
longitud que el tramo B-A, pero en realidad esto no es cierto,
por lo tanto se considera que su longitud virtual es de 200 y su
longitud real de 100. El tramo B-A tiene una longitud virtual
igual a su longitud real de 100.
Para tuberías de alimentación LVirtual =0

250
 Para predios que abastecen a un solo lado LVirtual =
LReal
Para predios que abastecen a los dos lados LVirtual
= 2*LReal
Sumando estas longitudes virtuales (de todo el sistema de red)
se podrá encontrar, se obtiene el coeficiente de caudal por
metro de tubería:
Qmax h
q [l/s*m]
 LVirtual
Donde:
q ; Coeficiente de caudal por metro de tubería
Qmax-h ; Caudal máximo horario (l/s)
L ; Sumatoria de longitudes virtuales (m)
Virtual

V.8 METODO DE LA LONGITUD EQUIVALENTE

Este método consiste en sumar a la extensión real de tubo, una

longitud equivalente a las pérdidas de carga que se ocasionan
en las piezas especiales existentes en la tubería. A cada pieza
especial corresponde una longitud ficticia y adicional (1). La
tabla siguiente muestra la longitud equivalente para diferentes
piezas.

LONGITUDES EQUIVALENTES A PERDIDAS LOCALES (EXPRESADAS

EN METROS DE TUBERIA RECTILINEA) (1)

251
 Diámetro Codo Codo 90° Codo 90° Codo 45°
90° radio me- radio corto
radio largo dio
mm pulg

13 01-feb 0.3 0.4 0.5 0.2

19 03-abr 0.4 0.6 0.7 0.3
25 1 0.5 0.7 0.8 0.4
32 1 1/4 0.7 0.9 1.1 0.5
38 1 1/2 0.9 1.1 1.3 0.6
50 2 1.1 1.4 1.7 0.8
63 2 1/2 1.3 1.7 2.0 0.9
75 3 1.6 2.1 2.5 1.2
100 4 2.1 2.8 3.4 1.5
125 5 2.7 3.7 4.2 1.9
150 6 3.4 4.3 4.9 2.3
200 8 4.3 5.5 6.4 3.0

EJEMPLO DE APLICACIÓN Se tiene una tubería de PVC de

250 mm (10“) de diámetro, clase 10. La tubería tiene
una longitud total de 1,150 m y conduce agua de
un tanque elevado 80 metros sobre el nivel de la descarga (ver
figura). Determinar el gasto que circula por la tubería,
si se tienen en el recorrido 3 codos de 22.5º,
6 de 45º y 2 de 90º, además de tener dos válvulas de compuerta
totalmente abiertas.

252
 Solución.

(a) Usando Darcy -Weisbach:

Aplicado del punto A al B, el término P/γ se elimina en ambos
miembros, ya que están expuestos a la
presión atmosférica, la velocidad en los tanques es muy
pequeña siendo despreciable por lo que la ecuación queda de
la siguiente manera al despejar para
pérdidas de carga. (Se asume una Tª de 20º C)

Z A  Z B  ht A B siendo ht A B  hf A B  hx A B
Z A  Z B  hf A B  hx A B

L V2 V2 m
f   Kx  Z A  Z B
D 2 g 2 g x 1
factorizando
V2 L m
(f   Kx)  Z A  Z B
2g D x 1

V2 8Q 2 8Q 2 L m

como  2 4  ( f   Kx)  Z A  Z B
2g  D g  2D4 g D x 1
al despejar el gasto

g Z A  ZB
Q 
Datos: 8 L 1 m
f  ( ) Kx
ZA- ZB= 80 m, D 5 D 4 x 1

L = 1150 m,
D N = 250 mm,

253
 D I = 231.7 mm
g= 9.81 m/s2
π= 3.1416
Kcodo 22.5º=0.2 (VER ANEXOS-PERDIDAS LOCALES),
Kcodo 45º =0.4 ,
Kcodo 90º = 0.9,
KVálvula de comp.= 0.2
ΣKx = 3(0.2)+6(0.4)+2(0.9)+2(0.2)=5.2, ΣKx/D4 =
5.2/(0.2317)4 =1,804.26
L/D5 = 1150/(0.2317)5 = 1.7221 x106 ,
La fórmula para calcular Q queda:

De Reynolds
Re = (1.2606 x 106 ) (Q/0.2317 m) = 5.44 x 106 Q , con Q en
m3 /s
ε/D = 0.0015 mm / 231.7 mm = 6.474 x 106
Como se puede observar el Número de Reynolds queda en
función del gasto, por lo tanto la f
también queda en función del gasto, para resolver este tipo de
problemas se hacen tanteos (a prueba y error). Regularmente
se tienen entre cinco y seis iteraciones.

254
 No. Q asumido Número de f de Swamee
(m³/s) Tabla de
Reynolds Jain Q
iteración
calculado

(m³/s)
1 0.100 5.4E+05 0.013006 0.200
2 0.200 1.1E+06 0.011614 0.211
3 0.211 1.1E+06 0.011520 0.212
4 0.212 1.2E+06 0.011514 0.212
5 0.212 1.2E+06 0.011513 0.212
6 0.212 1.2E+06 0.011513 0.212

En la 6ª iteración la f se mantiene constante por lo que la solución al

problema es:
Solución (a)

Q = 212 lps.

V= 5.02 m/s
V.9 PÉRDIDAS MENORES

En la mayor parte de los sistemas de flujo, la perdida de

energía primaria se debe a la fricción de conducto, como se vio
en el subtitulo anterior. Los demás tipos de perdidas
generalmente son pequeñas en comparación, y por
consiguiente se hace referencia a ellas como perdidas menores.
Las perdidas menores ocurren cuando hay un cambio en la
sección cruzada de la trayectoria de flujo o en la dirección de
flujo, o cuando la trayectoria de flujo se encuentra obstruida,

255
 como sucede con una válvula. La energía se pierde bajo estas
condiciones debido a fenómenos físicos bastante complejos. La
predicción teórica de la magnitud de estas perdidas también es
compleja, y por tanto, normalmente se usan los datos
experimentales.

V.9.1 COEFICIENTE DE RESISTENCIA

Las perdidas de energía son proporcionales a la cabeza de

velocidad del fluido al fluir éste alrededor de un codo, a trabes
de una dilatación o contracción de la sección de flujo, o a
través de una válvula. Los valores experimentales de perdidas
de energía generalmente se reportan en términos de un
coeficiente de resistencia “k”, de la siguiente forma:
h  K (v / 2 g )
L
2
5.9.1
donde hL es la perdia menor.

En algunos casos, puede haber mas de una velocidad de flujo,

como con las dilataciones o en las contracciones. Es de mayor
importancia que el estudiante sepa que velocidad se debe usar
con cada coeficiente de resistencia. Si la cabeza de velocidad
v / 2 g de la ecuación 5.7.1 se expresa en las mismas unidades
2

de metros, entonces la perdida de energía también estará en

metros o N*m/N de flujo de fluido. El coeficiente de

256
resistencia no tiene unidades, pues representa una constante de
proporcionalidad entre la perdida de energía y la cabeza de
velocidad. La magnitud de coeficiente de resistencia depende
de la geometría del dispositivo que ocasiona la perdida y
algunas veces depende de la velocidad de flujo.

V.9.2 DILATACION SUBITA

Al fluir un fluido de un conducto menor a uno mayor a través

de una dilatación súbita, su velocidad disminuye abruptamente
(Fig5-6), ocasionando una turbulencia, y por consiguiente, la
cantidad de perdida de energía, depende del cociente de los
tamaños de los dos conductos.

La menor perdida se calcula con la ecuación:

2
hL  K (v1 / 2 g ) 5.9.2

257
 Donde v1 es la velocidad de flujo promedio en el conducto
menor que esta delante de la dilatación. Las pruebas han
mostrado que el valor del coeficiente de perdida “k” depende
tanto de la proporción de los tamaños de los dos conductos
como de la magnitud de la velocidad de flujo (Fig5-7). Al
hacer ciertas suposiciones de simplificación respecto del
carácter de la corriente de flujo al expandirse a través de una
dilatación súbita, es posible predecir analíticamente el valor de
“k” a partir de la siguiente ecuación:
2

K  1   A1 / A2    1  ( D1 / D2 ) 2  2
5.9.3

258
 Fuente: Mecánica de Fluidos Aplicada-
Robert L. Mott

Los subíndices 1 y 2 se refieren a las secciones menores y

mayores, respectivamente (Fig5-6). Los valores para “k” de
esta ecuación concuerdan con los datos experimentales cuando
la velocidad v1 es aproximadamente 1.2m/s. A velocidades
mayores, los valores reales de “k” son menores que los valores

259
teóricos. Recomendamos que se usen los valores
experimentales si se conoce la velocidad de flujo.

EJEMPLO DE APLICACIÓN Determine la perdida

de energía que ocurrirá al fluir 10L/min de agua a
través de una dilatación súbita proveniente de un tubo
de cobre de 1 pulg. a un tubo de 3 pulg. (tipo K)

Solución. Usando el subíndice 1 para la sección

delantera de la dilatación y 2 para la sección corriente
debajo de la dilatación, obtenemos:
D1  25.3mm  0.0253m

A1  2.017 x10 4 m 2

D2  73.8mm  0.0738m

A1  4.282 x10 3 m 2
Q 100 1
v1   4
  3.32m / s
A1 5.017 x10 60000
2
v1 3.32 2
 m  0.56m
2 g 2(9.81)

Para encontrar un valor de k, se requiere la

proporción del diámetro, entonces tendremos:
D2 / D1  73.8 / 25.3  2.92mm

De la figura 5-7 obtenemos k= 0.72, por lo que

tendremos:
hL  K (v 2 / 2 g )  (0.72)(0.56)  0.40m

260
 Este resultado indica que 0.40N*m de energía se
disipa de cada newton de agua que fluye a través de
una dilatación súbita.

V.9.3 PERDIDA DE SALIDA

Durante el flujo de un fluido de un conducto hacia un gran

deposito o tanque (Fig5-8), su velocidad disminuye hasta casi
cero. En el proceso la energía cinética que el fluido poseía en
el conducto, indicada por la cabeza de velocidad, se disipa. Por
lo tanto, la perdida de energía para esta condición es:
2
hL  1.0(v1 / 2 g ) 5.9.4
Esta se denomina la pérdida de salida. El valor de k=1.0 se usa
sin importar la forma de la salida donde el conducto se conecta
con la pared del tanque.

EJEMPLO DE APLICACIÓN Determine la perdida

de energía que ocurrirá al fluir 100l/min de agua de

261
 un tubo de cobre de 1 pulg. (tipo k) en un tanque
mayor.
Solución. Usando la ecuación 5.9.4, tenemos:
2
hL  1.0(v1 / 2 g )

Sabemos que:
v1  3.32m / s y 2
v1 / 2 g  0.56m

Entonces, la perdida de energía será:

hL  1.0(0.56)  0.56m

V.9.4 DILATACION GRADUAL

Si la transición de un conducto menor a uno mayor puede

hacerse menos abrupta que la dilatación súbita de bordes
cuadrados, la perdida de energía se reduce. Esto normalmente
se hace colocando una sección entre los dos conductos (Fig5-
9). Las paredes en pendiente del cono tiende a guiar el fluido
durante la desaceleración y expansión de la corriente flujo.

La pérdida de energía para una dilatación gradual se calcula a

partir de:
2
hL  k (v1 / 2 g ) 5.9.5

262
donde v1 es la velocidad del conducto menor que está delante
de la dilatación. La magnitud de k depende tanto de la
proporción de diámetro, como del ángulo de cono “θ”, se dan
valores para diferentes casos en la figura 5-10 y la tabla 5.4
para “θ” y D2/D1.

263
 Fuente: Mecánica
de Fluidos Aplicada- Robert
L. Mott

La perdida de energía calculada de la ecuación 5.7.5 no incluye

la perdida debido a la fricción en las paredes de la transición.
Para lo ángulos de cono relativamente empinados, la longitud
de la transición es corta y por lo tanto, la perdida de fricción de
la pared es despreciable. Sin embargo, al disminuir el ángulo
del cono, la longitud de la transición se incrementa y la

264
 fricción de la pared se hace significativa. Tomando en cuenta
tanto la perdida de fricción de la pared como la pérdida debido
a la dilatación, podemos obtener la perdida de energía mínima
con un ángulo de cono de aproximadamente 7º.

Fuente: Mecánica de Fluidos Aplicada-

Robert L. Mott

EJEMPLO DE APLICACIÓN Determine la perdida

de energía que ocurrirá al fluir 100L/min de agua de
un tubo de cobre de 1 pulg (tipo k) en un tubo de
cobre de 3 pulg (tipo k ) a través de una dilatación
gradual con un ángulo de cono de 30º.

Solución. Usando anexos y los resultados de algunos

cálculos de los ejemplos anteriores, sabemos que:
v1  3.32m / s
2
v1 / 2 g  0.56m

D2 / D1  2.92

de la figura 5-10 obtenemos k= 0.48, entonces:

265
 2
hL  0.48(v1 / 2 g )  (0.48)(0.56)
hL  0.27 m

Comparada con la dilatación súbita, la perdida de

energía disminuye en 33% cuando utiliza la dilatación
gradual de 30º.

V.9.5 CONTRACCION SUBITA

La pérdida de energía debido a una contracción súbita (Fig5-

11), se calcula a partir de:

2
hL  k (v2 / 2 g ) 5.7.6

donde v2 es la velocidad en la corriente hacia abajo del

conducto menor a partir de la contracción. El coeficiente de
resistencia k depende de la proporción de los tamaños de lo
dos conductos y de la velocidad de flujo (Fig 5-12, tabla5.5).

266
El mecanismo mediante el cual se pierde energía debido a una
contracción súbita es bastante complejo. La figura 5-13
muestra los que pasa al converger la corriente de flujo. Las
líneas de la figura 5-13 representan las trayectorias de las
diversas partes de la corriente de flujo llamadas líneas de
trayectorias. Al aproximarse las líneas de trayectoria a la
contracción, asumen una trayectoria curva y la corriente total
continúa estrechándose durante cierta distancia más allá de la
contracción. Por lo tanto, la sección de cruce mínimo del flujo
es menor que la del conducto menor. La sección donde ocurre
esta área de flujo mínimo se denomina vena contracta
(Capitulo IV). Más allá de la vena, la corriente de flujo debe
desacelerar y dilatarse nuevamente para llenar el conducto. La
turbulencia ocasionada por la contracción y la posterior
dilatación genera la perdida de energía.

267
Fuente: Mecánica de Fluidos Aplicada-
Robert L. Mott

Fuente: Mecánica de Fluidos Aplicada-

Robert L. Mott

268
EJEMPLO DE APLICACIÓN Determine la perdida
de energía que ocurrirá al fluir 100l/min de agua de
un tubo de cobre de 3 pulg (tipo k) en un tubo de
cobre de 1 pulg (tipo K) a través de una contracción
súbita.

Solución. De la ecuación 5.9.6, tenemos:

2
hL  k (v2 / 2 g )

para el tubo de cobre sabemos que D1=73.8 mm,

D2=25.3 mm y
A2=5.017x10-4m2.

Entonces podemos encontrar los siguientes valores:

D2 / D1  2.92

Q 100 1
v2     3.32m / s
A2 5.017 x10  4 60000

269
 2
v2 / 2 g  0.56m

De la figura 5-12 encontramos k=0.42, entonces:

2
hL  k (v2 / 2 g )  (0.42)(0.56)  0.24m

V.9.6 CONTRACCION GRADUAL

La pérdida de energía en una contracción puede disminuirse

sustancialmente haciendo la contracción más gradual. La
figura 5-14 muestra una contracción de este tipo, formada
mediante una sección cónica entre los dos diámetros con
cambios abruptos en las junturas. El ángulo “θ” se denomina el
ángulo del cono.

La figura 5-15 muestra los datos para el coeficiente de

resistencia contra la proporción de diámetro para varios valores

270
del ángulo del cono. La pérdida de energía se calcula a partir
de la ecuación 5.7.6, donde el coeficiente de resistencia se
basa en la cabeza de velocidad en el conducto menor después
de la contracción.

Estos datos son para números de Reynolds mayores que 1.0 x

l05, Observe que para ángulos sobre el amplio intervalo de 15 a
40º, K= 0.05 o menos, un valor muy bajo. Para ángulos de
hasta 60º, K es menor que 0.08. Al disminuir el ángulo de cono
de la contracción por debajo de l5 grados, el coeficiente de
resistencia de hecho se incrementa, como se muestra en la
figura 5-16. La razón es que los datos incluyen los efectos
tanto de la turbulencia local ocasionada por la separación del
flujo como de la fricción del conducto. Para los ángulos de
cono menores, la transición entre los dos diámetros es muy
larga, lo que incrementa las pérdidas de fricción.

271
 Fuente: Mecánica de Fluidos Aplicada-
Robert L. Mott

El redondeo del extremo de la transición cónica para juntarla

con el conducto menor puede disminuir el coeficiente de
resistencia por debajo de los valores mostrados en la figura 5-
15.

272
Fuente: Mecánica de Fluidos Aplicada-
Robert L. Mott

V.9.7 PERDIDA DE ENTRADA

Un caso especial de una contracción ocurre cuando un fluido

fluye desde un depósito o tanque relativamente grande hacia
un conducto. El fluido debe acelerar desde una velocidad
relativamente despreciable a la velocidad de flujo del
conducto. La facilidad con que se realiza la aceleración
determina la cantidad de pérdida de energía y por lo tanto, el
valor del coeficiente de resistencia de entrada depende de la

273
 geometría de la entrada. La figura 5-17 muestra cuatro
configuraciones diferentes y el valor sugerido de K para cada
una. Las líneas de corriente ilustran el flujo de fluido hacia el
conducto y muestran que la turbulencia asociada con la
formación de una vena contracta en el tubo es una causa
importante de pérdida de energía. Para una entada bien
redondeada con Dr> 0.15, no se forma ninguna vena contracta,
la pérdida de energía es bastante pequeña, y usamos K = 0.04.

274
 EJEMPLO DE APLICACIÓN Determine la perdida
de energía que ocurrirá al fluir 100L/min de agua de
un deposito hacia un tubo de cobre de 1 pulg. (tipo k)
a) hacia un tubo de proyección hacia adentro.
b) a traves de una entrada bien redondeada.

275
 Solución. Tenemos los datos de los anteriores
ejemplos de aplicación por lo tanto

Para:
a)
D2  25.3mm

A2  5.017 x10 4 m 2
Q
v2   3.32m / s
A2
2
v2 / 2 g  0.56m

Para una entrada de proyección hacia adentro, K=1.0,

entonces:
2
hL  k (v2 / 2 g )  (1.0)(0.56)  0.56m

b)
Para una entrada bien redondeada, K=0.04, entonces:
2
hL  k (v2 / 2 g )  (0.04)(0.56)  0.02m

V.9.8 PERDIDAS POR VALVULAS

Los coeficientes de perdida por válvulas varían de acuerdo con

el tipo de válvula y según las posiciones, estas deben ser
276
proporcionados por los fabricantes. A falta de estos datos, se
puede utilizar los valores medios.

V.9.8.1 VALVULAS DE COMPUERTA

Según los valores mostrados en la figura 5-18

Fuente: Introducción a la Hidráulica-Texto

auxiliar de la materia

V.9.8.2 VALVULAS ESFERICAS

Según valores de la figura 5-19

Fuente: Introducción a la Hidráulica-Texto

auxiliar de la materia

277
V.9.8.3 VALVULAS MARIPOSAS
Según valores de la figura 5-20

Fuente: Introducción a la Hidráulica-Texto

auxiliar de la materia

278
 V.9.9 PÉRDIDAS POR CAMBIO DE DIRECCION

Estas perdidas son ocasionadas por distintos accesorios y estas

perdidas pueden variar según el sentido, posición y tamaño del
accesorio, así como la dirección del flujo.

V.9.9.1 BIFURCACIONES

Las bifurciones generalmente denominadas por su forma

suelen ser “tes o yes”, estás ocasionan perdidas que dependen
del ángulo que forman la tubería secundaria con la principal o
maestra, la relación entre los diámetros de ambas tuberías y de
la dirección de la corriente. esta perdida es mayor donde
sucede la unión de ambas tuberías, es decir donde esta la
bifurcación.
En la figura 5-21(a, b, c, d, e, f, g, h, i) se muestran varias
situaciones en las que se podrían presentar este tipo de
bifurcaciones.

279
280
281
282
283
284
285
 Fuente de graficas 5-21a-5-21i: Introducción
a la Hidráulica-Texto auxiliar de la materia

V.9.10 CODOS DE TUBERIA

A menudo es más conveniente curvar un conducto o tubo que

instalar un codo comercialmente hecho. La resistencia al flujo
de un codo depende de la proporción del radio “r” del codo con
el conducto dentro del diámetro D. En la figura 5-22 se
muestra que la resistencia mínima ocurre cuando la proporción
r/D es aproximadamente “3”. La resistencia se da en términos
de la proporción de longitud equivalente Le/D, y por lo tanto, la
ecuación:
k  ( Le / D) f T 5.9.7
Donde:

286
Le; longitud equivalente y es la longitud del conducto recto del
mismo diámetro nominal como el accesorio que tendría la
misma resistencia que está.
fT; factor de fricción en el conducto al cual está conectado el
accesorio.

Fuente: Mecánica de Fluidos Aplicada-

Robert L. Mott

Esta ecuación 5.7.7 debe usarse para calcular el coeficiente de

resistencia. La resistencia mostrada en figura 5-22 incluye
tanto la resistencia del codo como la resistencia debido a la
longitud del conducto en el codo. Cuando calculamos la
proporción r/D, “r” se define como el radio a la línea del centro
del conducto o tubo, denominado el radio medio (Fig5-23).

287
Esto es, si Ro es el diámetro externo del conducto o tubo y Ri
el interno:

r  Ri  Do / 2 5.9.8

r  Ro  Do / 2 5.9.9

r  ( Ro  Ri ) / 2 5.9.10

EJEMPLO DE APLICACIÓN un sistema de

distribución para propano liquido está hecho a partir
de tubería de acero de 11/4 pulg con un grosor de

288
 0.083 pulg. Se requieren varios codos de 90º para
ajustarse al otro equipo del sistema. Las
especificaciones requieren que el radio del interior de
cada codo sea de 200mm. Cuando el sistema lleva
160L/min de propano a 25ºC, calcule la pérdida de
energía en cada codo.

Solución. La ecuación de Darcy debe usarse para

calcular la perdida de energía con la proporción Le/D
para los codos (Fig5-23). Primero, determinamos r/D,
recordando que D es el diámetro interno del tubo y r
es el radio a la línea central del tubo.

De los anexos tenemos:

D  27.5mm  0.0275m

r  Ri  Do / 2

Donde Do=31.75mm, y el diámetro externo del tubo

lo obtenemos de lo anexos, tendremos:
r  200mm  (31.75mm) / 2  215.9mm

r / D  215.9mm / 27.5mm  7.85

De acuerdo con la figura 5-22, encontramos que la

proporción de longitud equivalente es
aproximadamente 23. Ahora debemos calcular la

289
 velocidad para completar la evaluación de la perdida
de energía de la ecuación de Darcy:
Q 160 L / min 1m 3 / s
v   4.48m / s
A 5.954 x10 4 m 2 60000 L / min

La aspereza relativa es:

D/Є=(0.0275m)/(4.6x10-5m) = 598
Entonces, podemos encontrar f  0.022 , del diagrama
T

de Moody visto anteriormente, en la zona de

completa turbulencia. Ahora puede calcularse la
perdida de energía:
Le 2 (4.48) 2
hL  f T (v / 2 g )  (0.022)(23)
D ( 2)(9.81)

hL  0.517m  0.517 N * m / N

V.10 GOLPE DE ARIETE

Supongamos que por un tubo con una válvula en su salida se

descarga un depósito. Si se cierra la válvula súbitamente, se
crea dentro de la tubería una presión dinámica, además de la
presión estática normal. A la presión dinámica se le llama
comúnmente “golpe de ariete” y se produce por la
transformación súbita de energía cinética a energía de presión
(Fig5-24). Si “p” es la fuerza de presión debida al golpe de
ariete, y “m” es la masa de agua contenida en el tubo, cuya
velocidad decrece en dv durante el intervalo dt,

290
 dv
pm
dt 5.10.1

Este golpe violento denominado Golpe de Ariete producido

sobre las paredes del conducto forzado, suele suceder cuando
hay una caída, o existen bombas en el sistema de las tuberías
estas paran y arrancan bruscamente.
145  v
hi 
EcuaciónEde
1 a
Alievi
d [m] 5.10.2
Et  e

Donde:
hi ; Sobre presión por golpe de ariete
v ; velocidad del agua en la tubería
Ea ; modulo de elasticidad del agua
Et MATERI E ; modulo de elasticidad del material
[Kg/cm
AL 2 de la tubería
]
d Acero 2,10*1006 ; Diámetro interior de la tubería
Hierro
e Fundido 9,3*1005
; Espesor de la tubería
Concreto
simple 1,25*1005
Asbesto-
cemento 3,28*1005
PVC 3,14*1004
291
Polietileno 5,2*1003
2,067*100
4
Agua
 Módulos de Elasticidad
para algunos materiales

Fuente: Abastecimiento de Agua

potable-Enrique Cesar Valdez

Al cerrar instantáneamente o para el equipo de bombeo, la

compresión del agua y expansión de la tubería comienzan en el
punto de cierre, transmitiéndose hacia arriba a una velocidad
determinada, conocida como velocidad de propagación de la
onda. El tiempo requerido para que la onda de presión regrese
a la válvula es:
2L
T 
a 5.10.3
Donde:
L ; Longitud de la tubería
a ; Celeridad de la onda de presión

292
 V.11 PLANTEAMIENTO PARA REDES

Es común el uso para sistemas de distribución de agua potable,

utilizar unas mallas de tuberías, estas consisten en tuberías
interconectadas configurando una red. Una red de tuberías es
en realidad un sistema de tuberías múltiples, por lo tanto su
solución se basa en lo siguientes principios.
La cantidad de agua que entra a la red debe ser igual a la que
sale de ella. En sistemas de agua potable, los caudales de salida
se denominan consumos, este principio corresponde al
principio de continuidad.
Q  C
a 5.11.1
Donde:
Qa; caudales que alimentan la red.
C: Consumos.
Esta ecuación (5.9.1) representa el balance global de la red.
En cada nodo debe cumplirse la ecuación de la continuidad. Al
asignarle signo a los gastos que concurren a un nodo (positivos
afluentes y negativos efluentes)
Q  C  0
ij 5.11.2
donde para un nodo determinado, el subíndice “j” indica el
nodo hacia donde va el flujo y el subíndice “i” desde donde
proviene; el termino C representa los caudales retirados.

293
en cada nodo existe un solo nivel de energía, en consecuencia,
las perdidas de un nodo a otro, independientemente del
camino de unión seleccionado, son iguales. Esto lleva a la
siguiente conclusión “ si se considera perdida de energía
positiva a aquella que ocurre en sentido contrario a las agujas
del reloj y viceversa, se podrá establecer para cada malla de
red:
5.11.3
nij
h  0  K Q
ij ij ij

Una red que tenga “n” nodos y “m” mallas independientes se

podrá, en principio, plantear n+m ecuaciones; sin embargo
como debe cumplirse la ecuación del balance global de la red,
únicamente “n-1” ecuaciones de nodo son independientes entre
si, dado lo cual solo se dispone de “m+(n-1” ecuaciones. Es
importante hacer notar que si las incógnitas son los gastos, las
ecuaciones de nodo son lineales, mientras las ecuaciones de
malla no lo son. Para la resolución de problemas de redes
existen varios métodos, plantearemos algunos de ellos.

V.11.1 SISTEMA Q

Consiste en suponer unos gastos iniciales cuales quiera en las

tuberías de la red que cumplen con las ecuaciones de nodo.
Cada gasto inicial es afectado por un incremento correctivo Q

294
que se considera constante para las tuberías de cada malla. Los
incrementos se calculan mediante las ecuaciones de malla.
El sistema de ecuaciones constará entonces de “m” ecuaciones
no lineales con un numero igual de incógnitas que son los Q
de cada malla.
K (Q  Q )   0
ij oij I
nij
I

K (Q  Q )   0
ij oij I
nij
II

……..
……..
K ij (Qoij  QI ) nij  m
0

Donde:
Qo; gasto inicial supuesto.
I, II,.. m; cantidad de mallas

V.11.2 SISTEMA H

Las incógnitas en este sistema son los niveles o alturas de

energía “H” en cada nodo. Los gastos en las tuberías se
expresan en función de estos niveles y a continuación se
establecen las ecuaciones de nodo correspondientes. Como
solo existen “n-1” ecuación es independientes de nodo, es
necesario conocer el nivel de energía en un nodo cualquiera.
Se debe tomar en cuenta que:

295
 1
 Hi  H j  nij

;
nij
hij  H i  H j  K ij Q Qij   
ij
 K 
 ij 

El sistema H es:
  
1

 H i  H j
nij
  C  0
  K ij 
 
 1

 H H 
1
nij 
 i j
  C  0
  K ij 
 
 2

……
……
  
1

 H i  H j
nij
  C  0
  K ij 
 
  n 1

Donde:
1, 2,…., n-1; representan los n-1 nodos con ecuaciones
independientes, las ecuaciones de este sistema son no
lineales.

V.11.3 SISTEMA Q

Es el planteamiento general, es decir, tanto las ecuaciones de

nodo como las de malla, considerando como incógnitas a los
gastos en las tuberías. El sistema constará de m+(n-1)
ecuaciones y se puede expresar de la siguiente manera:
Ecuaciones de nodo:
 Q  C   0 ij 1

296
 Q ij C  2
0

…..
 Q ij  C n 1
0

Ecuaciones de malla:
 K Q  ij
n ij
ij I 0

K Q  ij
nij
ij II 0

………
 K Q  ij
nij
ij m 0

En este sistema existen n-1 ecuaciones lineales y m ecuaciones

no lineales.

V.11.4 METODOS PARA SOLUCIONAR REDES DE

TUBERIAS

a) Método de Cross

Este método consiste en hacer lineales las ecuaciones del

sistema Q mediante el binomio de Newton aplicando la
siguiente forma:
n ij n 1
K ij (Qoij  Qm )  K ij (Qoij  nij Qoijij Qm )  0

con lo cual ignora todos los términos del binomio a partir del
tercero. Como Q es constante para cada malla, pude sacarse
m

fuera del paréntesis, resultando entonces que el factor de

corrección de gasto en una determinada malla es:

297
 n
K ij Qoijij
Q   n 1
nij K ij Qoijij

Para aplicar esta ecuación es oportuno recordar que

únicamente los sumandos del numerador, es decir las perdidas
de energía, tienen signo; el denominador es siempre positivo.
Los gastos corregidos se calculan mediante la expresión
siguiente:
QIij  Qoij  Qm

Donde el signo negativo corresponde a tuberías con perdidas

de carga negativas y viceversa. Si el valor de Q resulta Iij

negativo significa que el flujo tiene sentido contrario al

supuesto en la iteración previa.

b) Método de Newton- Raphson

Este método puede aplicarse a cualquier sistema, en esta

ocasión solo presentaremos su aplicación al sistema H.
El método se basa en el proceso iterativo de Newton para
encontrar el valor de la variable que anula una función de ella,
proceso que fue ampliado a funciones de varias variables por
Raspón, extendiendo de esta forma su campo de aplicación,
1
H  H  nij

uno de los cuales es la solución de redes

 C de tuberías. El sistema
i j
F  
j
 K ij

H es el formado por las ecuaciones de nodo; en consecuencia,

el problema reside en encontrar los valores de los niveles de

298
 energía o cotas piezometricas “H” que anulen la siguiente
función:

En cada nodo de la red existirá una función Fj, por lo tanto, si

lo gastos alimentadores de la tubería son conocidos, existirán
“n-1” funciones y cuando no lo sean es decir mallas
imaginarias el numero será “n”.
Se comienza a trabajar con unos valores de H cuales quiera
supuestos. Sin embargo, las suposiciones deden ser razonables,
en caso contrario, el método podría no converger. Estos
valores iniciales se corrigen sustrayéndoles un incremento
finito H que, puede calcularse mediante la siguiente ecuación,
expresada en forma matricial:

 F1 F1 F1 

 H ....
H 2 H j 
  H1   F 1 
'
 1
 F2 F2
....
F2  
H   F '2 
 H1 H 2 H j   2   
   .   . 
 . . . . *  = 
 .   . 
 . . . .   . 
   . 
 . . . .  H j   ' 
 Fj F j Fj    F j 
 
....
 H1 H 2 
H j 


Donde F1’, F2’, etc., son los valores de las funciones Fj para los
valores iniciales de H. Para los efectos de colocar signos, es

299
oportuno recordar que la función Fj representa la ecuación de
la continuidad en cada uno de los nodos, esto implica que el
gasto efluente es positivo y el gasto afluente es negativo,
además que el subíndice “i” corresponde al nodo de donde
proviene el flujo, y “j” hacia el nodo donde se dirige.

300