Flujo en Tuberías PDF

UNIVERSIDAD NACIONAL DE HUANCAVELICA
FACULTAD DE CIENCIAS DE LA INGENIERÍA

ESCUELA PROFESIONAL DE INGENIERÍA CIVIL-HUANCAVELICA
"Año de la lucha contra la corrupción y la impunidad"
RESOLUCIÓN DE EXÁMENES DE MECÁNICA

DE FLUIDOS I Y LABORATORIO
CÁTEDRA: MECÁNICA DE FLUIDOS I
CATEDRÁTICO: Ing. R. Kennedy Gomez Tunque

ESTUDIANTE: Angel Sullcaray Ichpas
CICLO: V
SECCIÓN: B
HUANCAVELICA 2019
1 PROBLEMA 1.
Determinar el gasto que fluye en cada uno de los ramales del sistema
mediante la ecuación de Hazen-William. CH =100 para todas las tu-
berı́as. Obviar la carga cinética en la ecuación de Bernoulli.
SOLUCIÓN:
· utilizando la ecuación de Hazen-Williams
Q = 0,000426CH D2,63 S 0,54
como:
hf
S=
L
0,000426CH D2,63 h0,54
f
Q=
L0,54
0,000426CH D2,63
K= (coeficiente de resistencia)
L0,54
Q = Kh0,54
f
MECANICA DE FLUIDOS I 2 ING.CIVIL

hallando los coeficientes de resistencia
0,000426(100)(10)2,63
K1 = = 12,4984
20,54
0,000426(100)(8)2,63
K2 = = 8,11796
1,50,54
0,000426(100)(12)2,63
K3 = = 20,1883
20,54
0,000426(100)(8)2,63
K4 = = 7,8399
1,50,54
•Asumiendo nueva altura piezométrica Zp = 3540m
Hallando las pérdidas de carga
hf 1 = Z1 − Zp = 20m
hf 2 = Z1 − Zp = 20m
hf 3 = Zp − Z3 = 40m
hf 4 = Zp − Z2 = 30m
hallando los caudales
Q1 = 12,4984(10)0,54 = 63,0102 lt/s

Q2 = 8,11796(10)0,54 = 40,9264 lt/s
Q3 = 20,1883(50)0,54 = 147,9831 lt/s
Q4 = 7,8399(40)0,54 = 49,19897 lt/s
verificando la ecuación de continuidad

Q1 + Q2 = Q3 + Q4
Error=(Q3 + Q4 ) − (Q1 + Q2 )
E0 = 93,2455
NOTA: Debe aumentar Q1 y Q2 , por lo tanto hf en las tuberı́as 1

y 2 aumenta.

•Asumiendo nueva altura piezométricaZp = 3520m
hf 1 = 40m −→ Q1 = 91,615 lt/s

hf 2 = 40m −→ Q2 = 59,5058 lt/s
hf 3 = 20m −→ Q3 = 101,7785 lt/s
hf 4 = 10m −→ Q4 = 27,1838 lt/s

E1 = −22,1585
hf 1 = 30m −→ Q1 = 78,4332 lt/s

hf 2 = 30m −→ Q2 = 50,9439 lt/s
hf 3 = 30m −→ Q3 = 126,6908 lt/s
hf 4 = 20m −→ Q4 = 39,5245 lt/s

E2 = 36,8382
hf 1 = 35m −→ Q1 = 85,2415 lt/s

hf 2 = 35m −→ Q2 = 55,3661 lt/s
hf 3 = 25m −→ Q3 = 114,8120 lt/s
hf 4 = 15m −→ Q4 = 33,8376 lt/s

E3 = 8,042

Utilizando la interpolación de lagrange se halla la altura piezométrica

Zp con los datos obtenidos
(Zp − E1 )(Zp − E2 )(Zp − E3 )

L0 (Zp ) = Z0
(E0 − E1 )(E0 − E2 )(E0 − E3 )
(Zp − E0 )(Zp − E2 )(Zp − E3 )
L1 (Zp ) = Z1
(E1 − E0 )(E1 − E2 )(E1 − E3 )
(Zp − E0 )(Zp − E1 )(Zp − E3 )
L2 (Zp ) = Z2
(E2 − E0 )(E2 − E1 )(E2 − E3 )
(Zp − E0 )(Zp − E1 )(Zp − E2 )
L3 (Zp ) = Z3
(E3 − E0 )(E3 − E1 )(E3 − E2 )
n
X
P (Zp ) = Li (Zp)fi (Zp)
i=0
Cuando E = 0, Zp = 3523,64m

Hallando el gasto que fluye por cada una de las tuberı́as con la nueva
altura piezométrica
hf 1 = 36,36m −→ Q1 = 87,014 lt/s
hf 2 = 36,36m −→ Q2 = 56,518 lt/s
hf 3 = 23,64m −→ Q3 = 111,396 lt/s
hf 4 = 13,64m −→ Q4 = 32,145 lt/s
2 PROBLEMA 2.
Se quiere calcular el caudal del canal de irrigación artesanal que es
abastecida desde un tanque elevado, si la tuberı́a es de fierro fundido
nuevo de 4” de diámetro. Considerar las pérdidas menores.
Viscosidad cinemática υ = 10−6 m2 /s, ks = 0,00025m PVC
Usar la ecuación de Colebrook-White (pérdidas por fricción).
SOLUCIÓN:
Aplicando Bernoulli entre los puntos A y B.
PA VA2 P V2
+ + ZA = B + B + ZB + hf + hk VB = V
γ 2g γ 2g
V 2
9,15 + 3,05 = + hf + hk .........(1)
2g
Por Darcy-Weisbach.
L V2
hf = f
D 2g

Pérdidas locales
•Codo de 90◦
V2
hk = K K = 0,9
2g
•Válvula globo medio abierto
Le
= 160 → K = 160f
D
Reemplzando en (1).
V2 927,2 V 2 V2 V2 V2
12,2 = +f + 160f + 0,9 +
2g (0,102) 2g 2g 2g 2g
12,2 = V 2 (0,147854 + 471,612f )
s
12,2
V =
0,147854 + 471,612f
utilizando la ecuación de Colebrook-White
1 k 2,51
√ = −2log ( s + √ )
f 3,7D Re f
Ks 2,51 1
F (f ) = 2log ( + √ )+ √
3,7D Re f f
−2,51
√
0 Re f f 0,5
F (f ) = − √
ks 2,51 f f
LN (10)( + √ )
3,7D Re f
utilizando el método de Newton Raphson
F ( f0 )
f = f0 −
F 0 ( f0 )
• Iteración 1.
asumiendo f0 = 0,02
s
12,2
V = = 1,1285m/s
0,147854 + 471,612(0,02)
πV D2 π (1,1285)(0,102)2
Q= = = 0,0092m3 /s
4 4

VD 1,1285(0,102)
Re = = = 115107
υ 10−6
• Iteración 2.
0,00025 2,51 1
F (f ) = 2log ( + √ )+ √
3,7(0,102) 115107 0,02 0,02
−2,51
√
0 115107(0,02) 0,02 0,5
F (f ) = − √
0,00025 2,51
LN (10)( + √ ) 0,02 0,02
3,7(0,102) 115107 0,02
f = 0,0249487
s
12,2
V = = 1,0119m/s
0,147854 + 471,612(0,0249487)
π (1,0119)(0,102)2
Q= = 0,0083m3 /s
4
Re = 103213,8
• Iteración 3.
f = 0,02613834
s
12,2
V = = 0,9889m/s
0,147854 + 471,612(0,02613834)
π (0,9889)(0,102)2
Q= = 0,0081m3 /s
4
Re = 100867,8
• Iteración 4.
f = 0,02621284
s
12,2
V = = 0,9875m/s
0,147854 + 471,612(0,02621284)
π (0,9889)(0,102)2
Q= = 0,0081m3 /s
4
Error = (f − f0 ) = 0,02621284 − 0,026113834 = 0,000099

3 PROBLEMA 3.
Se desea calcular el caudal de una acantarilla (D=120cm) a presión
que pasa por debajo de una carretera, sabiendo que la pérdida de car-
ga a la entrada puede ser a 0.8 veces la carga de velocidad,considerar
la alcantarilla de tuberı́a PVC. Sabiendo: h1 = 3m h2 = 1,6m
L = 50m. ks = 0,0000015m −6
υ = 10 m /s2
Usar la ecuación de D-W para las pérdidas principales.
SOLUCIÓN:
Aplicando Bernoulli entre los puntos A y B.

PA VA2 P V2
+ + ZA = B + B + ZB + hf + hk VB = V
γ 2g γ 2g
ZA − ZB = hf + hk .........(1)
Por Darcy-Weisbach.
L V2
hf = f
D 2g
Pérdidas locales
• En la entrada
hk (entrada) = 0,8(carga cinética)
V2
hk (entrada) = 0,8
2g

• En la salida
Ensanchamiento brusco
K=1
V2
hk (salida) =
2g
reemplazando en la ecuación (1)
50 V 2 V2 V2
1,4 = f + 0,8 +
(1,2) 2g 2g 2g
1,4 = V 2 (2,124f + 9,177 ∗ 10−2 )
s
1,14
V =
2,124f + 9,177 ∗ 10−2
1 k 2,51
√ = −2log ( s + √ )
f 3,7D Re f
• Iteración 1.
asumiendo f0 = 0,02
s
1,14
V = = 3,229m/s
2,124(0,02) + 9,177 ∗ 10−2
πV D2 π (3,229)(1,2)2
Q= = = 3,652m3 /s
4 4
VD 3,229(1,2)
Re = = = 3875146,42
υ 10−6
• Iteración 2.
f = 0,00200315
s
1,14
V = = 3,818m/s
2,124(0,00200315) + 9,177 ∗ 10−2
π (3,818)(1,2)2
Q= = 4,319m3 /s
4

3,818(1,2)
Re = = 4582146,139
10−6
• Iteración 3.
f = 0,0041701
s
1,14
V = = 3,730m/s
2,124(0,0041701) + 9,177 ∗ 10−2
π (3,730)(1,2)2
Q= = 4,219m3 /s
4
3,730(1,2)
Re = = 4475972,14
10−6
• Iteración 4.
f = 0,0069275
s
1,14
V = = 3,626m/s
2,124(0,0069275) + 9,177 ∗ 10−2
π (3,626)(1,2)2
Q= = 4,101m3 /s
4
3,626(1,2)
Re = = 4351141,38
10−6
• Iteración 5.
f = 0,008801
s
1,14
V = = 3,560m/s
2,124(0,008801) + 9,177 ∗ 10−2
π (3,560)(1,2)2
Q= = 4,026m3 /s
4
3,560(1,2)
Re = = 4272048,85
10−6

• Iteración 6.
f = 0,0092587
s
1,14
V = = 3,544m/s
2,124(0,0092587) + 9,177 ∗ 10−2
π (3,544)(1,2)2
Q= = 4,009m3 /s
4
3,544(1,2)
Re = = 4253373,51
10−6
• Iteración 7.
f = 0,0092827
s
1,14
V = = 3,5437m/s
2,124(0,0092827) + 9,177 ∗ 10−2
π (3,544)(1,2)2
Q= = 4,0078m3 /s
4
3,544(1,2)
Re = = 4252400,996
10−6
Error = (f − f0 ) = 0,0092827 − 0,0092587 = 0,000024

4 PROBLEMA 4.
Se tiene un sistema en serie y paralelo del cual se desea determinar
los caudales y pérdidas de carga en cada uno de las tuberı́as mediante
Hazen-William.
SOLUCIÓN:
· utilizando la ecuación de Hazen-Williams
Q = 0,000426CD2,63 S 0,54
hf
como: S=
L
0,000426CD2,63 h0,54
f
Q=
L0,54
QL0,54
h0,54
f =
0,000426CD2,65
hf = KQ1,851852
L
donde: K=
5,73 ∗ 10−7 C 1,851852 D4,90740741
Por el principio de conservación de masa
Q1 + Q2 + Q3 − 20 = 0............(1)
Q1 + Q4 + Q5 − 20 = 0............(2)

Las pérdidas de carga
hf 2 − hf 3 = 0.............(3)
hf 1 − hf 2 − hf 4 = 0.........(4)
hf 4 − hf 5 = 0..........(5)
Hallando los coeficientes de resistencia

1,2
K1 = = 0,015328
5,73 ∗ 10−7 (100)1,851852 (8)4,90740741
0,8
K2 = −7 1,851852
= 0,000243
5,73 ∗ 10 (120) (16)4,90740741
0,9
K3 = −7 1,851852
= 0,000526
5,73 ∗ 10 (120) (14)4,90740741
0,5
K4 = −7 1,851852
= 0,002137
5,73 ∗ 10 (100) (10)4,90740741
0,82
K5 = = 0,001432
5,73 ∗ 10−7 (100)1,851852 (12)4,90740741
reemplazando en las ecuaciones (3),(4) y (5) tenemos un sistema de
ecuaciones
f1 (Q) = Q1 + Q2 + Q3 − 20
f2 (Q) = Q1 + Q4 + Q5 − 20
f3 (Q) = 0,015328Q1,851852
1 − 0,000243Q1,851852
2 − 0,002137Q1,851852
4
f4 (Q) = 0,000243Q1,851852
2 − 0,000526Q1,851852
3
f5 (Q) = 0,002137Q1,851852
4 − 0,001432Q1,851852
5

Empleando el método de Newton-Raphson para la solución de ecuacio-

nes no lineales
[Xi+1 ] = [Xi ] − [Z ]−1 [fi ]

 ∂f (Q ) ∂f1 (Q2i ) ∂f1 (Q3i ) ∂f1 (Q4i ) ∂f1 (Q5i )
−1
1 1i
 ∂Q1 ∂Q2 ∂Q3 ∂Q4 ∂Q5 
     ∂f2 (Q1i ) ∂f2 (Q2i ) ∂f2 (Q3i ) ∂f2 (Q4i ) ∂f2 (Q5i )   
Q1i+1 Q 1i   f1 i
 ∂Q1 ∂Q2 ∂Q3 ∂Q4 ∂Q5
Q2i+1 Q2i f2i

   ∂f3 (Q1i ) ∂f3 (Q2i ) ∂f3 (Q3i ) ∂f3 (Q4i ) ∂f3 (Q5i )  
Q3i+1  = Q3i  −
 ∂Q1
 f3i 

Q 4i + 1
 
Q 4i

 ∂f (Q ) ∂Q2 ∂Q3 ∂Q4 ∂Q5  
f4 i

∂f4 (Q2i ) ∂f4 (Q3i ) ∂f4 (Q4i ) ∂f4 (Q5i )

Q 5i + 1 Q 5i  4 1i  f5 i
 ∂Q1 ∂Q2 ∂Q3 ∂Q4 ∂Q5 
∂f5 (Q1i ) ∂f5 (Q2i ) ∂f5 (Q3i ) ∂f5 (Q4i ) ∂f5 (Q5i )
 
∂Q1 ∂Q2 ∂Q3 ∂Q4 ∂Q5
1 1 1 0 0
 
1 0 0 1 1
0,851852
−0,00045Q0,851852 −0,00396Q0,851852
 
[Z ] = 0,02839Q1 2 0 4 0 
0,00045Q0,851852 −0,00097Q0,851852
 
0 2 3 0 0
0 0 0 0,00396Q0,851852
4 −0,00265Q0,851852
5
p
Error = (Q1i+1 − Q1i )2 + (Q2i+1 − Q2i )2 + (Q3i+1 − Q3i )2 + (Q4i+1 − Q4i )2 + (Q5i+1 − Q5i )2
• Iteración 1.
asumiendo Q1i = Q2i = Q3i = Q4i = Q5i = 10lt/s
     −1  
Q1i+1 10 1 1 1 0 0 10
Q2i+1  10   1 0 0 1 1   10 
Q3i+1  = 10  − 0,20181 −0,00320 0 −0,02814 0   0,92057 
−0,00693
Q  10   0 0,00320 0 0
 −0,02012 
4i+1
Q 5i + 1 10 0 0 0 0,02814 −0,01885 5,01235
   
Q1 4,98813
 i+1
Q2i+1 8,57535
  
 
   
Q3
 i+1

 =
6,43652


Q4i+1 6,93191
   
 
Q5i+1 8,07995
Error = 7,27644

• Iteración 2.
   
Q1 3,23795
 i+1
Q2i+1 10,14425
  
 
   
Q3
 i+1

 =
 6,61780


Q4i+1  7,45619
   
 
Q5i+1 9,30586
Error = 2,708351
• Iteración 3.
   
Q1 2,91530
 i+1
Q2i+1 10,29739
  
 
   
Q3
 i+1

 =
 6,78730


Q4i+1  7,62296
   
 
Q5i+1 9,46173
Error = 0,45650
• Iteración 4.
   
Q1 2,90282
 i+1
Q2i+1 10,30563
  
 
   
Q3
 i+1

 =
 6,79156


Q4i+1  7,62816
   
 
Q5i+1 9,46902
Error = 0,0179438
• Iteración 5.
     
Q1 2,90280 Q1 = 2,9 lt/s
 i+1
Q2i+1 10,30564 Q2 = 10,3 lt/s
    
  
     
=  −→ los caudales son :  Q3 = 6,8
 6,79157
Q3  lt/s 
 i+1    
Q4i+1  7,62817  Q4 = 7,6 lt/s
     
  
Q5i+1 9,46903 Q5 = 9,5 lt/s
Error = 2,828 ∗ 10−5

5 PROBLEMA 5.
Calcular los caudales que llegan o parten de cada uno de los depósitos
mediante Hazen-William.
SOLUCIÓN:
· utilizando la ecuación de Hazen-William
Q = 0,000426CH D2,63 S 0,54
como:
hf
S=
L
0,000426CH D2,63 h0,54
f
Q=
L0,54
0,000426CH D2,63
K= (coeficiente de resistencia)
L0,54
Q = Kh0,54
f

hallando los coeficientes de resistencia
0,000426(100)(15,75)2,63
K1 = = 54,3876
1,20,54
0,000426(100)(7,87)2,63
K2 = = 12,7533
0,60,54
0,000426(130)(9,84)2,63
K3 = = 18,1903
1,50,54
0,000426(90)(11,81)2,63
K4 = = 9,6265
60,54
Como la elevación de la lı́nea de alturas piezométricas en E no pue-
de determinarse, por ser desconocidos todos los caudales, el problema
se resolverá por tanteos. En el primero es conveniente elegir como al-
tura piezométrica en E, 84 m. Con esto, el caudal que sale o entra en
el reservorio B será nulo, lo que reduce el número de cálculos.
• Para una altura piezométrica ZE = 84,0 m
hf 1 = ZA − ZE = 6m
hf 2 = ZE − ZB = 0m
hf 3 = ZE − ZC = 15m
hf 4 = ZE − ZD = 33,5m
Q1 = 54,3876(6)0,54 = 143,073 lt/s

Q2 = 12,7533(0)0,54 = 0 lt/s
Q3 = 18,1903(15)0,54 = 78,574 lt/s
Q4 = 9,6265(33,5)0,54 = 64,134 lt/s
Evaluando la ecuación de continuidad

Q1 = Q2 + Q3 + Q4
Error=Q1 − Q3 − Q4 − Q2
E0 = 0,365

De los valores de estos caudales se infiere que la altura piezométrica en

E debe ser mayor, de forma que se reduzca el caudal en 1, aumente el
que va a B y circule cierto caudal hacia C Y D.
hf 1 = 5,9m
hf 2 = 0,1m
hf 3 = 15,1m
hf 4 = 33,6m
Q1 = 54,3876(5,9)0,54 = 141,780 lt/s

Q2 = 12,7533(0,1)0,54 = 3,683 lt/s
Q3 = 18,1903(15,1)0,54 = 78,860 lt/s
Q4 = 9,6265(33,6)0,54 = 64,240 lt/s

E0 = −5,003
hf 1 = 5,8m
hf 2 = 0,2m
hf 3 = 15,2m
hf 4 = 33,7m
Q1 = 54,3876(5,8)0,54 = 140,477 lt/s

Q2 = 12,7533(0,2)0,54 = 5,355 lt/s
Q3 = 18,1903(15,2)0,54 = 79,138 lt/s
Q4 = 9,6265(33,7)0,54 = 64,341 lt/s


E0 = −8,357
hf 1 = 5,95m
hf 2 = 0,05m
hf 3 = 15,05m
hf 4 = 33,55m
Q1 = 54,3876(5,95)0,54 = 142,427 lt/s

Q2 = 12,7533(0,05)0,54 = 2,534 lt/s
Q3 = 18,1903(15,05)0,54 = 78,715 lt/s
Q4 = 9,6265(33,55)0,54 = 64,186 lt/s

E0 = −3,007
Utilizando la interpolación de lagrange se halla la altura piezométrica

ZE con los datos obtenidos

Cuando E = 0, ZE = 84,003m
Q1 = 143,034 lt/s
Q2 = 0,554 lt/s
Q3 = 78,582 lt/s
Q4 = 64,137 lt/s

6 PROBLEMA 6.
En la figura mostrada se representa el abastecimiento mediante dos
depósitos escalonados a dos redes de riego.La toma A abastece una red
de riego por aspersión para lo que se requiere una altura piezométrica
de 35m.La toma G abastece una red una red de riego por goteo para lo
que se requiere una altura piezométrica de 20m.El depósito S mantiene
su nivel constante alimentado por el depósito T, que a su vez es de
grandes dimensiones. Se pide:
- Caudales circulantes por el sistema de tuberı́as en litros/segundo.

- Diámetro de la tuberı́a que une ambos depósitos T a S.
Despreciar sumandos cinéticos, υ = 10−6 m2 /s
SOLUCIÓN:
Por continuidad
4Q
Q = VA −→ V =
πD2
Por Darcy-Weisbach.
L V2 8f LQ2
hf = f = 2 5
D 2g π gD

v
u hf π 2 gD 5
u
−→ Q = t
8f L
VD
Re =
υ
1 k 2,51
√ = −2log ( s + √ )
f 3,7D Re f
Ks 2,51 1
F (f ) = 2log ( + √ )+ √
3,7D Re f f
−2,51
√
0 Re f f 0,5
F (f ) = − √
ks 2,51 f f
LN (10)( + √ )
3,7D Re f
utilizando el método de Newton Raphson
F ( f0 )
f = f0 −
F 0 ( f0 )
• Para una altura piezométrica ZN = 40m
−→ Las pérdidas de carga son:
hf 2 = 4m
hf 3 = 5m
hf 4 = 20m
NOTA: Utilizando el tipo II se hallan los caudales
PARA LA TUBERÍA 2
asumiendo f0 = 0,02
• Iteración 1. • Iteración 2. • Iteración 3.
f0 = 0,02 f0 = 0,0152909 f0 = 0,015571

Q2 = 0,032 m3 /s Q2 = 0,037 m3 / Q2 = 0,036 m3 /
V = 0,808 m/s V = 0,924 m/s V = 0,916 m/s
Re = 181800 Re = 207912,98 Re = 206034,62
f = 0,0152909 f = 0,015571 f = 0,0156022

• Iteración 4. • Iteración 5.
f0 = 0,0156022 f = 0,0156053
Q2 = 0,036 m3 / Q2 = 0,036 m3 /
V = 0,915 m/s V = 0,915 m/s
Re = 205828,46 Re = 205805,97
f = 0,0156053 f = 0,0156056
error ( %) = 1,97 ∗ 10−2 error ( %) = 1,92 ∗ 10−3
PARA LA TUBERÍA 3
En 5 iteraciones se optimiza el cálculo con un error porcentual de

9,39 ∗ 10−3 obteniendose los siguientes valores
Q3 = 0,64 m3 /s f = 0,013757
PARA LA TUBERÍA 4
En 5 iteraciones se optimiza el cálculo con un error porcentual de

4,365 ∗ 10−4 obteniendose los siguientes valores
Q3 = 0,005 m3 /s f = 0,0183932
Verificando la ecuación de continuidad
Error = Q2 − Q3 − Q4
E0 = −0,033
Se obtienen los siguientes caudales en cada tuberı́a
Q2 = 0,046 m3 /s f = 0,014954
Q3 = 0,048 m3 /s f = 0,0144764
Q4 = 0,005 m3 /s f = 0,0186177
E1 = −0,007

Q2 = 0,0495 m3 /s f = 0,0147176
Q3 = 0,0386 m3 /s f = 0,0150927
Q4 = 0,005 m3 /s f = 0,0187413
E2 = 0,0059
• Para una altura piezométrica ZN = 37,5m
Q2 = 0,048 m3 /s f = 0,0148305
Q3 = 0,044 m3 /s f = 0,0147484
Q4 = 0,005 m3 /s f = 0,0186784
E3 = −0,001
Utilizando la interpolación de lagrange se halla la altura piezométri-

ca ZN con los datos obtenidos

(Zp − E1 )(Zp − E2 )(Zp − E3 )

L0 (Zp ) = Z0
(E0 − E1 )(E0 − E2 )(E0 − E3 )
(Zp − E0 )(Zp − E2 )(Zp − E3 )
L1 (Zp ) = Z1
(E1 − E0 )(E1 − E2 )(E1 − E3 )
(Zp − E0 )(Zp − E1 )(Zp − E3 )
L2 (Zp ) = Z2
(E2 − E0 )(E2 − E1 )(E2 − E3 )
(Zp − E0 )(Zp − E1 )(Zp − E2 )
L3 (Zp ) = Z3
(E3 − E0 )(E3 − E1 )(E3 − E2 )
n
X
P (Zp ) = Li (Zp)fi (Zp)
i=0
Cuando E = 0, ZN = 37,42m
Q2 = 0,048 m3 /s f = 0,0148117
3
Q3 = 0,043 m /s f = 0,0147978
Q4 = 0,005 m3 /s f = 0,0186883

Por dato del problema:
El depósito S mantiene su nivel constante alimentado por el depósi-

to T, entonces se infiere que Q2 = Q1 .
Utilizando TIPO III se halla el diámetro de la tuberı́a 1.
asumiendo f = 0,02
se calcula el diámetro
v
2
5 8f LQ
u
u
D= t
π 2 ghf
luego se calcula el el número de reynold
VD
Re =
υ
Se reemplaza en la ecuación de colebrook-white y se halla el nuevo
valor del factor de frición f .
RESOLVIENDO
Datos
f0 = 0,02
hf 1 = 51,7 − 44 = 7,7m
Q1 = 0,048 m3 /s
L1 = 1200m
K = 0,0025
υ = 10−6
RESULTADOS
D = 0,26m
f = 0,0378098
Error = 9,1511 ∗ 10−7

PROGRAMA PARA HALLAR PÉRDIDAS DE CARGA, CAU-

DALES Y DIÁMETRO MEDIANTE DARCY-WEISBACH
EXPORT FLUIDOS_I()
BEGIN
LOCAL Re,ks,vs,V,f,fo,D;
LOCAL Q,QQ,fx,dfx,hf,g;
LOCAL MENU;
LOCAL Ei,Emax;
ks:=0.003;
D:=0.3;
hf:=5;
Q:=0.125;
L:=400;
vs:=0.00000113;
fo:=0.02;
g:=9.807;//gravedad
//INICIO DEL PROGRAMA

RECT;
FOR I FROM 1 TO 240 DO
LINE_P(0,I-0.99,320,I-0.99,RGB(13,I,I))
END;
TEXTOUT_P("ESCUELA PROFESIONAL DE INGENIERÍA CIVIL",18,20+I,3,RGB(0,162,232));
TEXTOUT_P("MECÁNICA DE FLUIDOS I",80,35+I,3,RGB(0,114,168));
END;
LINE_P(40,60+I,276,60+I,RGB(138,0,138));
LINE_P(40,83+I,276,83+I,RGB(138,0,138));
END;
TEXTOUT_P("FLUJO DE FLUIDOS EN TUBERIAS",53+I,65+I,3,RGB(255,255,255));
END;
TEXTOUT_P("FLUJO DE FLUIDOS EN TUBERIAS",53,68,3,RGB(225,0,113));
FOR J FROM 1 TO 2 DO
TEXTOUT_P("Desarrollado por: ",36,90+J,2,RGB(120,24,34));
END;
TEXTOUT_P("ANGEL SULLCARAY ICHPAS",133,91,2,RGB(255,234,79));
TEXTOUT_P("Presione Enter para ingresar",146,223,2,RGB(255,0,0));

TEXTOUT_P("!!!",298,217,5,RGB(255,0,0));
TEXTOUT_P("Docente: Ing. R. Kennedy Gómez",60,200,2,RGB(130,0,130));
//LA ECUACIÓN DE DARCY
TEXTOUT_P("?",18,121,4,RGB(0,64,128));
TEXTOUT_P("hf = f",34,123,3,RGB(0,0,0));
TEXTOUT_P("L",73,117,3,RGB(0,0,0));
LINE_P(72,132,80,132,RGB(0,0,0));
TEXTOUT_P("D",73,132,3,RGB(0,0,0));
TEXTOUT_P("V",89,117,3,RGB(0,0,0));
LINE_P(86,132,99,132,RGB(0,0,0));
TEXTOUT_P("2g",86,132,3,RGB(0,0,0));
TEXTOUT_P("2",97,117,1,RGB(0,0,0));
TEXTOUT_P("Ecuación de Darcy Weisbach.",130,130,1,RGB(0,0,157));
//LA ECUACIÓN DE COLEBROOK-WHITE
TEXTOUT_P("?",18,159,4,RGB(0,64,128));
TEXTOUT_P("1",34,155,3,RGB(0,0,0));
LINE_P(32,170,47,170,RGB(0,0,0));
LINE_P(35,172,45,172,RGB(0,0,0));
LINE_P(32,185,35,172,RGB(0,0,0));
LINE_P(32,185,31,183,RGB(0,0,0));
TEXTOUT_P("f",37,172,3,RGB(0,0,0));
TEXTOUT_P("=",50,162,3,RGB(0,0,0));
TEXTOUT_P("-2log",58,162,3,RGB(0,0,0));
TEXTOUT_P("[",94,157,6,RGB(0,0,0));
//DENTRO DEL PARENTESIS
TEXTOUT_P("Ks",105,158,2,RGB(0,0,0));
LINE_P(100,170,124,170,RGB(0,0,0));
TEXTOUT_P("3.7D",101,171,2,RGB(0,0,0));
TEXTOUT_P("+",129,163,3,RGB(0,0,0));
TEXTOUT_P("2.51",145,158,2,RGB(0,0,0));
LINE_P(140,170,170,170,RGB(0,0,0));
//f
LINE_P(159,172,169,172,RGB(0,0,0));
LINE_P(156,184,159,172,RGB(0,0,0));
LINE_P(156,184,155,182,RGB(0,0,0));
TEXTOUT_P("f",161,174,2,RGB(0,0,0));
TEXTOUT_P("Re",140,172,2,RGB(0,0,0));
TEXTOUT_P("]",170,157,6,RGB(0,0,0));

TEXTOUT_P("Ecuación de Colebrook-White",182,166,1,RGB(0,0,157));
WAIT(-1);
//UTILIZANDO EL COMANDO CHOOSE
REPEAT
CHOOSE(MENU,"ELEGIR TIPO",{"TIPO I","TIPO II","TIPO III","SALIR"});
//TIPO I (CÁLCULO DE LA PÉRDIDA DE CARGA)

IF MENU==1 THEN
LOCAL MI:=[[0]];
INPUT({{fo,[0],{50,30,1}},
{D,[0],{50,30,2}},
{Q,[0],{50,30,3}},
{L,[0],{50,30,4}},
{vs,[0],{50,30,5}},
{ks,[0],{50,30,6}}},
"DATOS INICIALES",
{"fo: ","D[m]: ","Q[m3/s]: ","L[m]: ","Vs[m2/s]: ","Ks: "},
{"ingrese F.fricción inicial","ingrese Diámetro","ingrese caudal",
"ingrese longitud","ingrese viscosidad cinemática",
"ingrese rugosidad absoluta"});
//INGRESANDO AL BUCLE WHILE
V:=4*Q/(?*Dˆ2);
Re:=V*D/vs;
Ei:=0.1;
Emax:=0.0000001;
I:=0;
//LIMPIANDO LA PANTALLA
PRINT;
WHILE Ei>Emax DO
fx:=1/(foˆ0.5) +2*LOG((ks/(3.7*D))+(2.51/(Re*foˆ0.5)));
dfx:=(-2.51/(Re*fo*foˆ0.5))/((ks/(3.7*D))+(2.51/(Re*foˆ0.5))*LN(10))
-0.5/(fo*foˆ0.5);
f:=fo-(fx)/(dfx);
Ei:=ABS(f-fo)*100/f;
fo:=f;
I:=I+1;
MI(I,1):=I;
MI(I,2):=ROUND(f,7);
MI(I,3):=Ei;

END;
//RESULTADOS
EDITMAT(MI,{"RESULTADOS",{},{"ITERACIÓN","f","Error(%)",""}});
hf:=f*(L/D)*(Vˆ2/(2*g));
RECT;
LINE_P(0,I-0.99,320,I-0.99,RGB(24,I,I));
END;
TEXTOUT_P("RESULTADOS",110,50+I,4,RGB(0,0,0));
TEXTOUT_P("? La pérdida de carga es:",50,80,3, RGB(0,0,135));
TEXTOUT_P("hf = ",70,100,2,RGB(0,0,0));
TEXTOUT_P(ROUND(hf,5),95,100,2,RGB(0,0,0));
TEXTOUT_P("? El factor de fricción es:",50,130,3,RGB(0,0,135));
TEXTOUT_P("f = ",70,150,2,RGB(0,0,0));
TEXTOUT_P(ROUND(f,7),90,150,2,RGB(0,0,0));
TEXTOUT_P("? El número de Reynolds es:",50,180,3,RGB(0,0,135));
TEXTOUT_P("Re = ",70,200,2,RGB(0,0,0));
TEXTOUT_P(ROUND(Re,2),95,200,2,RGB(0,0,0));
END;
WAIT(-1);
END;
//TIPO II (CÁLCULO DEL CAUDAL)

IF MENU==2 THEN
LOCAL MII:=[[0]];
INPUT({ {fo,[0],{50,30,1}},
{D,[0],{50,30,2}},
{hf,[0],{50,30,3}},
{L,[0],{50,30,4}},
{vs,[0],{50,30,5}},
{ks,[0],{50,30,6}}},
"DATOS INICIALES",
{"fo: ","D[m]: ","hf[m]: ","L[m]: ","Vs[m2/s]: ","Ks: "},
{"ingrese F.fricción inicial","ingrese Diámetro","ingrese pérdida de carga",
"ingrese longitud","ingrese viscosidad cinemática",

Ei:=0.1;
Emax:=0.0000001;
I:=0;
PRINT;
WHILE Ei>Emax DO
//CÁLCULO DE LAS VARIABLES
QQ:=((hf*?ˆ2*g*Dˆ5)/(8*fo*L))ˆ0.5;
V:=4*QQ/(?*Dˆ2);
Re:=V*D/vs;
-0.5/(fo*foˆ0.5);
f:=fo-(fx)/(dfx);
fo:=f;
I:=I+1;
MII(I,1):=I;
MII(I,2):=ROUND(f,7);
MII(I,3):=Ei;
END;
//RESULTADOS
EDITMAT(MII,{"RESULTADOS",{},{"ITERACIÓN","f","Error(%)",""}});
RECT;
LINE_P(0,I-0.99,320,I-0.99,RGB(24,I,I));
END;
TEXTOUT_P("? El caudal es:",50,80,3, RGB(0,0,135));
TEXTOUT_P("Q = ",70,100,2,RGB(0,0,0));
TEXTOUT_P(ROUND(QQ,5),95,100,2,RGB(0,0,0));
TEXTOUT_P("f = ",70,150,2,RGB(0,0,0));
TEXTOUT_P("Re = ",70,200,2,RGB(0,0,0));
END;
WAIT(-1);

END;
//TIPO III (CÁLCULO DEL DIÁMETRO DE LA TUBERÍA)

IF MENU==3 THEN
LOCAL MIII:=[[0]];
INPUT({{fo,[0],{50,30,1}},
{Q,[0],{50,30,2}},
{hf,[0],{50,30,3}},
{L,[0],{50,30,4}},
{vs,[0],{50,30,5}},
{ks,[0],{50,30,6}}},
"DATOS INICIALES",
{"fo: ","Q[m3/s]: ","hf[m]: ","L[m]: ","Vs[m2/s]: ","Ks: "},
{"ingrese F.fricción inicial","ingrese caudal","ingrese pérdida de carga"
,"ingrese longitud","ingrese viscosidad cinemática",
Ei:=0.1;
Emax:=0.0000001;
I:=0;
PRINT;
WHILE Ei>Emax DO
//CÁLCULO DE LAS VARIABLES
D:=((8*fo*L*Qˆ2)/(?ˆ2*g*hf))ˆ0.2;
V:=4*Q/(?*Dˆ2);
Re:=V*D/vs;
-0.5/(fo*foˆ0.5);
f:=fo-(fx)/(dfx);
fo:=f;
I:=I+1;
MIII(I,1):=I;
MIII(I,2):=ROUND(f,7);
MIII(I,3):=Ei;
END;
//RESULTADOS
EDITMAT(MIII,{"RESULTADOS",{},{"ITERACIÓN","f","Error(%)",""}});
RECT;


LINE_P(0,I-0.99,320,I-0.99,RGB(24,I,I));
END;
TEXTOUT_P("? La diámetro de la tuberı́a es:",50,80,3, RGB(0,0,135));
TEXTOUT_P("D = ",70,100,2,RGB(0,0,0));
TEXTOUT_P(ROUND(D,5),95,100,2,RGB(0,0,0));
TEXTOUT_P("f = ",70,150,2,RGB(0,0,0));
TEXTOUT_P("Re = ",70,200,2,RGB(0,0,0));
END;
WAIT(-1);
END;
UNTIL MENU==4;
END;
PROGRAMA PARA INTERPOLAR MEDIANTE EL POLINO-

MIO DE LAGRANGE PARA 4 DATOS
EXPORT LAGRANGE()
BEGIN
PRINT;
LOCAL Zo,Z1,Z2,Z3;
LOCAL Eo,E1,E2,E3;
LOCAL Lo,L1,L2,L3;
LOCAL fo,f1,f2,f3;
LOCAL Px;
LOCAL LA,a;
INPUT({{Zo,[0],{20,20,1}},
{Z1,[0],{20,20,2}},
{Z2,[0],{20,20,3}},
{Z3,[0],{20,20,4}},
{Eo,[0],{60,20,1}},
{E1,[0],{60,20,2}},

{E2,[0],{60,20,3}},
{E3,[0],{60,20,4}}},
"INGRESAR DATOS",{"Zo: ","Z1: ","Z2: ","Z3: ","Eo: ","E1: ","E2: ","E3: "});
//POLINOMIO DE LAGRANGE
Lo:=((X-E1)*(X-E2)*(X-E3))/((Eo-E1)*(Eo-E2)*(Eo-E3));
L1:=((X-Eo)*(X-E2)*(X-E3))/((E1-Eo)*(E1-E2)*(E1-E3));
//POLINOMIO GENERAL
Px:=(Lo*Zo)+(L1*Z1)+(L2*Z2)+(L3*Z3);
PRINT("P(x): "+Px);
END;

ANEXO

Flujo en Tuberías PDF

Cargado por

Copyright:

Formatos disponibles

Flujo en Tuberías PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Flujo en Tuberías PDF

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE HUANCAVELICA

FACULTAD DE CIENCIAS DE LA INGENIERÍA

"Año de la lucha contra la corrupción y la impunidad"

RESOLUCIÓN DE EXÁMENES DE MECÁNICA

CÁTEDRA: MECÁNICA DE FLUIDOS I

CATEDRÁTICO: Ing. R. Kennedy Gomez Tunque

· utilizando la ecuación de Hazen-Williams

Q = 0,000426CH D2,63 S 0,54

MECANICA DE FLUIDOS I 2 ING.CIVIL

hallando los coeficientes de resistencia

Hallando las pérdidas de carga

hallando los caudales

Q1 = 12,4984(10)0,54 = 63,0102 lt/s

verificando la ecuación de continuidad

NOTA: Debe aumentar Q1 y Q2 , por lo tanto hf en las tuberı́as 1

MECANICA DE FLUIDOS I 3 ING.CIVIL

•Asumiendo nueva altura piezométricaZp = 3520m

hf 1 = 40m −→ Q1 = 91,615 lt/s

verificando la ecuación de continuidad

•Asumiendo nueva altura piezométrica Zp = 3530m

hf 1 = 30m −→ Q1 = 78,4332 lt/s

verificando la ecuación de continuidad

•Asumiendo nueva altura piezométrica Zp = 3525m

hf 1 = 35m −→ Q1 = 85,2415 lt/s

verificando la ecuación de continuidad

MECANICA DE FLUIDOS I 4 ING.CIVIL

Utilizando la interpolación de lagrange se halla la altura piezométrica

(Zp − E1 )(Zp − E2 )(Zp − E3 )

MECANICA DE FLUIDOS I 5 ING.CIVIL

MECANICA DE FLUIDOS I 6 ING.CIVIL

MECANICA DE FLUIDOS I 7 ING.CIVIL

MECANICA DE FLUIDOS I 8 ING.CIVIL

Usar la ecuación de D-W para las pérdidas principales.

Aplicando Bernoulli entre los puntos A y B.

hk (entrada) = 0,8(carga cinética)

MECANICA DE FLUIDOS I 9 ING.CIVIL

MECANICA DE FLUIDOS I 10 ING.CIVIL

MECANICA DE FLUIDOS I 11 ING.CIVIL

MECANICA DE FLUIDOS I 12 ING.CIVIL

MECANICA DE FLUIDOS I 13 ING.CIVIL

Las pérdidas de carga

Hallando los coeficientes de resistencia

MECANICA DE FLUIDOS I 14 ING.CIVIL

Empleando el método de Newton-Raphson para la solución de ecuacio-

[Xi+1 ] = [Xi ] − [Z ]−1 [fi ]

MECANICA DE FLUIDOS I 15 ING.CIVIL

Error = 2,828 ∗ 10−5

MECANICA DE FLUIDOS I 16 ING.CIVIL

· utilizando la ecuación de Hazen-William

Q = 0,000426CH D2,63 S 0,54

MECANICA DE FLUIDOS I 17 ING.CIVIL

hallando los coeficientes de resistencia

• Para una altura piezométrica ZE = 84,0 m

Hallando las pérdidas de carga

hallando los caudales

Q1 = 54,3876(6)0,54 = 143,073 lt/s

Evaluando la ecuación de continuidad

MECANICA DE FLUIDOS I 18 ING.CIVIL

De los valores de estos caudales se infiere que la altura piezométrica en