UNIVERSIDAD NACIONAL
AGRARIA LA MOLINA
ESCUELA DE POSTGRADO
MAESTRIA EN RECURSOS HIDRICOS
SEMESTRE : 2015 I
DOCENTE
: Ing. Jose Arapa
ALUMNO
: Aguedo Tahua Robert Henry.
CURSO
: Mtodos numricos en ingeniera
de recursos hdricos.
TEMA
: Hidrulica de tuberas calculo y
anlisis de redes.
CODIGO
: 20150824
La Molina - Lima
Maestra en recursos Hdricos - Universidad Nacional Agraria la Molina
ndice
MATLAB. .................................................................................................................................. 3
I.
1.1. Ejemplo 4.3. Calculo el dimetro que debe tener una tubera nueva, de cemento enlucido
(K=0.0004 m) para conducir 2 m3/s. la viscosidad del agua es de 1.2 x 10-6 m2/s. la longitud
de la tubera dela tubera es de 1000 m. la perdida de carga admisible. Es de 25 m. ............. 3
1.2. Ejemplo 4.4. Que presin se requiere impulsar 20 lps a lo largo de una tubera lisa,
horizontal, de 2 de dimetro. La longitud del tramo es 300 m. la viscosidad de agua es de 10
-6
m2/s. ............................................................................................................................................. 4
1.3. Ejemplo 4.5. Calculo el gasto sistema mostrado en la figura. La viscosidad del agua es
1.2x10-6 m2/s. la tubera es lisa. Considerar nicamente las prdidas de carga continuas. El
dimetro de ls tubera de descarga es de 2 cm. ........................................................................... 5
II.
EXCEL. .................................................................................................................................. 6
2.1. Ejemplo 4.3. Calculo el dimetro que debe tener una tubera nueva, de cemento enlucido
(K=0.0004 m) para conducir 2 m3/s. la viscosidad del agua es de 1.2 x 10^-6 m2/s. la longitud
de la tubera de la tubera es de 1000 m. la perdida de carga admisible. Es de 25 m. ............ 6
2.2. Ejemplo 4.4. Que presin se requiere impulsar 20 lps a lo largo de una tubera lisa,
horizontal, de 2 de dimetro. La longitud del tramo es 300 m. la viscosidad de agua es de 10
-6 m2/s. ........................................................................................................................................... 7
2.3. Ejemplo 4.5. Calculo el gasto sistema mostrado en la figura. La viscosidad del agua es
1.2x10-6 m2/s. la tubera es lisa. Considerar nicamente las prdidas de carga continuas. El
dimetro de la tubera de descarga es de 2 cm. .......................................................................... 8
Maestra en recursos Hdricos - Universidad Nacional Agraria la Molina
I. MATLAB.
1.1. Ejemplo 4.3. Calculo el dimetro que debe tener una tubera nueva, de cemento enlucido
(K=0.0004 m) para conducir 2 m3/s. la viscosidad del agua es de 1.2 x 10-6 m2/s. la longitud
de la tubera dela tubera es de 1000 m. la perdida de carga admisible. Es de 25 m.
clc;clear;format short, format compact
%ASIGNAR EL VALOR f hasta igualarse.
f=0.0171; %coeficiente de Darcy - weisbach
%DATOS DE ENTRADA
H=25;%perdida de carga
g=9.81; %aceleracion de la gravedad
Ks=0.0004; %Rugosidad absoluta PVC
Q=2; %Caudal
L=1000;%Longitud
v=1.2e-6;%Viscosidad cinematica
%CALCULOS PRELIMIARES
S=H/L;%pendiente
d=(0.0827*f*Q^2/S)^(1/5);
x=1/sqrt(f);
A=pi*d^2/4;%Area de tuberia
V=Q/A; %velocidad flujo
Re=V*d/v; %Nmero de reynodls
i=1;tol=0.0001;err=100;
while err>=tol
fx=-2*log10(Ks/(3.71*d)+2.51*x/Re)-x;
dfx=-(2/log(10))*((2.51/Re)/(Ks/(3.71*d)+2.51*x/Re))-1;
x1=x;
x=x-fx/dfx;
err=abs((x-x1)/x*100);
T(i)=i;X1(i)=x1;F(i)=fx;DF(i)=dfx;Err(i)=err;
i=i+1;
end
TABLA=[T' X1' F' DF' Err']
f1=(1/x)^2;
[f f1 d]
disp('
f
f
d')
TABLA =
1.0000
2.0000
ans =
0.0171
f
7.6472
7.6372
-0.0100
0.0000
0.0171
f
0.7429
d
-1.0050
-1.0050
0.1306
0.0000
Maestra en recursos Hdricos - Universidad Nacional Agraria la Molina
1.2. Ejemplo 4.4. Que presin se requiere impulsar 20 lps a lo largo de una tubera lisa,
horizontal, de 2 de dimetro. La longitud del tramo es 300 m. la viscosidad de agua es de
10 -6 m2/s.
clc;clear;format short, format compact
%ASIGNAR EL VALOR f hasta igualarse.
f=0.0136; %coeficiente de Darcy - weisbach
%DATOS DE ENTRADA
d=2*2.54/100; %Diametro tuberia
Ks=1.5e-6; %Rugosidad absoluta PVC
Q=20e-3; %Caudal
g=9.81;%Aceleracion de la gravedad
v=1e-6;%Viscosidad cinematica
L=300;%longitud
%CALCULOS PRELIMIARES
x=1/sqrt(f);
A=pi*d^2/4;%Area de tuberia
V=Q/A; %velocidad flujo
Re=V*d/v; %Nmero de reynodls
i=1;tol=0.0001;err=100;
while err>=tol
fx=-2*log10(Ks/(3.71*d)+2.51*x/Re)-x;
dfx=-(2/log(10))*((2.51/Re)/(Ks/(3.71*d)+2.51*x/Re))-1;
x1=x;
x=x-fx/dfx;
err=abs((x-x1)/x*100);
T(i)=i;X1(i)=x1;F(i)=fx;DF(i)=dfx;Err(i)=err;
i=i+1;
end
TABLA=[T' X1' F' DF' Err']
f1=(1/x)^2;
hf=f1*L*V^2/(2*g*d);
[f f1 hf]
disp('
f
f
hf')
TABLA =
1.0000
2.0000
ans =
0.0136
f
8.5749
8.5857
0.0117
0.0000
0.0136
f
397.5855
hf
-1.0855
-1.0854
0.1257
0.0000
Maestra en recursos Hdricos - Universidad Nacional Agraria la Molina
1.3. Ejemplo 4.5. Calculo el gasto sistema mostrado en la figura. La viscosidad del agua es
1.2x10-6 m2/s. la tubera es lisa. Considerar nicamente las prdidas de carga continuas. El
dimetro de ls tubera de descarga es de 2 cm.
clc;clear;format short, format compact
%ASIGNAR EL VALOR f hasta igualarse.
f=0.0186; %coeficiente de Darcy - weisbach
%DATOS DE ENTRADA
hf=5;%Perdida de carga continua
d=2/100; %Diametro tuberia
Ks=0; %Rugosidad absoluta PVC
g=9.81;%Aceleracion de la gravedad
v=1.2e-6;%Viscosidad cinematica
L=4;%longitud
%CALCULOS PRELIMIARES
x=1/sqrt(f);
A=pi*d^2/4;%Area de tuberia
V=((2*g*d*hf)/(f*L))^0.5; %velocidad flujo
Re=V*d/v; %Nmero de reynodls
i=1;tol=0.0001;err=100;
while err>=tol
fx=-2*log10(Ks/(3.71*d)+2.51*x/Re)-x;
dfx=-(2/log(10))*((2.51/Re)/(Ks/(3.71*d)+2.51*x/Re))-1;
x1=x;
x=x-fx/dfx;
err=abs((x-x1)/x*100);
T(i)=i;X1(i)=x1;F(i)=fx;DF(i)=dfx;Err(i)=err;
i=i+1;
end
TABLA=[T' X1' F' DF' Err']
f1=(1/x)^2;
Q=V*A;
[f f1 Q]
disp('
f
f
Q')
TABLA =
1.0000
2.0000
ans =
0.0186
f
7.3324
7.3347
0.0026
0.0000
0.0186
f
0.0016
Q
-1.1185
-1.1184
0.0321
0.0000
Maestra en recursos Hdricos - Universidad Nacional Agraria la Molina
II. EXCEL.
2.1. Ejemplo 4.3. Calculo el dimetro que debe tener una tubera nueva, de cemento enlucido
(K=0.0004 m) para conducir 2 m3/s. la viscosidad del agua es de 1.2 x 10^-6 m2/s. la
longitud de la tubera de la tubera es de 1000 m. la perdida de carga admisible. Es de 25
m.
Solucin
Perdida de carga.
Aceleracin de la gravedad.
Caudal.
H=
g =
Q=
25
9.81
2
m
m2/s
m3/s
Longitud.
Viscosidad cinemtica.
Rugosidad del material.
Coeficiente de Darcy - Weisbach.
L=
v=
Ks =
f=
1000
1.20E-06
0.0004
0.0171
m
m2/s
m
Pendiente
S=
0.025
D=
0.74
A=
V=
0.433
4.614
m2
m/s
Dimetro
rea.
Velocidad del flujo
Nmero de Reynolds.
Re =
2.86E+06
xi
f(xi)
F'(xi)
xi+1
7.647191129
-0.0100270
-1.01E+00
7.63721432
0.130634
7.637214317
0.0000000
-1.01E+00
7.63721432
0.000000
f = 1/x^2
X=
f=
7.637
0.0171
D=
0.74 m
er
Maestra en recursos Hdricos - Universidad Nacional Agraria la Molina
2.2. Ejemplo 4.4. Que presin se requiere impulsar 20 lps a lo largo de una tubera lisa,
horizontal, de 2 de dimetro. La longitud del tramo es 300 m. la viscosidad de agua es de
10 -6 m2/s.
Solucin
Dimetro tubera
Aceleracin de la gravedad.
d=
g =
0.0508
9.81
m
m2/s
Caudal.
Longitud.
Q=
L=
0.02
300
m3/s
m
Viscosidad cinemtica.
Rugosidad del material.
Coeficiente de Darcy Weisbach.
v=
Ks =
1.00E-06
1.50E-06
m2/s
m
rea.
Velocidad del flujo
Nmero de Reynolds.
A=
V=
Re =
f=
0.0136
0.002
m2
9.868
m/s
5.01E+05
xi
f(xi)
F'(xi)
xi+1
8.574929257
0.0117114
-1.09E+00
8.59E+00
0.125667
8.59E+00
0.0000005
-1.09E+00
8.58571914
0.000005
8.59E+00
0.0000000
-1.09E+00
8.58571914
f = 1/x^2
X=
f=
8.586
0.0136
hf = f * L * V^2 / (2 * g * d)
hf =
397.59 m
er
Maestra en recursos Hdricos - Universidad Nacional Agraria la Molina
2.3. Ejemplo 4.5. Calculo el gasto sistema mostrado en la figura. La viscosidad del agua
es 1.2x10-6 m2/s. la tubera es lisa. Considerar nicamente las prdidas de carga
continuas. El dimetro de la tubera de descarga es de 2 cm.
Solucin
Diametro tuberia
Aceleracin de la gravedad.
Perdida de carga por dif. Cota.
Longitud.
Viscosidad cinematica.
Rugosidad del material.
Coeficiente de Darcy Weisbach.
Area.
Velocidad del flujo
Nmero de Reynolds.
d=
g =
hf =
L=
v=
Ks =
f=
A=
V=
Re =
0.02
9.81
5
4
1.20E-06
0
m
m2/s
m
m
m2/s
m
0.0186
0.000314
m2
5.135
m/s
8.56E+04
xi
f(xi)
F'(xi)
xi+1
7.332355751
0.0026333
-1.12E+00
7.33E+00
0.032100
7.33E+00
0.0000000
-1.12E+00
7.3347102
0.000001
7.33E+00
0.0000000
-1.12E+00
7.3347102
f = 1/x^2
X=
f=
7.335
0.0186
Q = V*A
Q=
0.00161 m
er