Problema 19

Código de programación Matlab

% PROBLEMA 19

clc
clear
format compact

% Agua fluye a 75°F por gravedad de un deposito A hacia un depósito B a
% través del sistema de tuberías mostrado. Las longitudes de tramo recto de
% la tubería, cuyo diámetro es de 4 in, es 200ft, mientras que la de la
% tubería, cuyo diámetro es de 6 in, es 100ft, ambas tuberías son de hierro
% fundido. Calcular el flujo volumétrico a través del sistema.

dz=input(‘valor de la diferencia de alturas ‘);
%dz=50
Di1=input(‘Valor del diámetro interno de la tubería 1 ‘);
%Di1=.3355
S1=input(‘Valor del área de la tubería correspondiente a la cédula 2 ‘);
%S1=.0884
Di2=input(‘Valor del diámetro interno de la tubería 1 ‘);
%Di2=.5054
S2=input(‘Valor del área de la tubería correspondiente a la cédula 2 ‘);
%S2=.2006
L1=input(‘Valor del tramo recto de la tubería 1 ‘);
%L1=200
L2=input(‘Valor del tramo recto de la tubería 2 ‘);
%L2=100
f1=input(‘Valor del factor de fricción 1’)
%f1=0.017
f2=input(‘Valor del factor de fricción 2’)
%f2=0.015

disp(‘Este problema solo puede hacer el cálculo si existen 2 tuberías con diferente diámetro’)

disp(‘Ahora introduce los accesorios para la primera tubería ‘)

n=input(‘¿Cuántos accesorios tienes en el sistema? ‘);
g1=32.2;
gc=32.2;
g=zeros(n,5)
for j=1:n

disp(‘Antes de introducir los datos recuerda asignar un número determinado para cada accesorio por ejemplo 1=entrada, 2=codo de 90°, 3= VCTA, etc.’)

g(j,1)=input(‘Nombre del accesorio ‘);
g(j,2)=input(‘Cantidad del accesorio anterior ‘);
g(j,3)=input(‘L/D del accesorio anterior ‘);
g(j,4)=input(‘Valor de Ki Tabla A-50 Crane ‘);
g(j,5)=g(j,2)*g(j,4)
fprintf(‘\n Accesorio Cantidad L/D Ki Kt\n’)
g
end

disp(‘Ahora introduce los accesorios para la segunda tubería ‘)

nn=input(‘Cuantos accesorios tienes en el sistema ‘);
h=zeros(nn,5)
for j=1:nn
disp(‘Antes de introducir los datos recuerda asignar un número determinado
para cada accesorio por ejemplo 1=entrada, 2=codo de 90°, 3= VCTA, etc.’)
h(j,1)=input(‘Nombre del accesorio ‘);
h(j,2)=input(‘Cantidad del accesorio anterior ‘);
h(j,3)=input(‘L/D del accesorio anterior ‘);
h(j,4)=input(‘Valor de Ki Tabla A-50 Crane ‘);
h(j,5)=h(j,2)*h(j,4)
fprintf(‘\n Accesorio Cantidad L/D Ki Kt\n’)
h
end
ktot=g(:,5);
ktot3=sum(ktot)
ktot1=h(:,5);
ktot4=sum(ktot1)
Ktr1=f1*(L1/Di1)
Kta=Ktr1+ktot3
Ktr2=f2*(L2/Di2)
Ktb=Ktr2+ktot4
a=dz*(g1/gc);
b=1/(2*gc);
c=(Kta*S2*S2)/(S1*S1*2*gc);
x=Ktb/(2*gc);
disp(‘Las unidades para la velocidad son [ft/s]’)
Ub=sqrt(a/(b+c+x))
disp(‘Las unidades para el gasto son [ft^3/s]’)
Q=Ub*S2