Problema 17

Código de programación Matlab

% PROBLEMA 17

clc
clear
format compact

% Una tubería de 10in de diámetro se divide paralelamente en 4 como se
% muestra, uno de los ramales tiene un diámetro de 6 in y 2000 ft de
% longitud. La tubería se une de nuevo en B en una tubería de 10 in de
% diámetro. En el punto A, la altura es de 200 ft y presión es de 30 psig,
% mientras que en el punto B la altura es de 50 ft y presión es de 10 psig.
% Determinar la descarga en la tubería de 10 in; el fluido es agua a 60°F.

disp(‘Ahora ingresaras los datos correspondientes dependiendo cuantos diámetros diferentes existan en el sistema’)

n=input(‘Cuantos diámetros diferentes hay en el sistema ? ‘)
g=zeros(7,n)
for j=1:n
g(1,j)=input(‘Valor de la tubería ‘)
g(2,j)=input(‘Valor del tramo recto por tubería ‘)
g(3,j)=input(‘Valor del Diámetro interno por tubería ‘)
g(4,j)=input(‘Valor del área por tubería ‘)
g(5,j)=g(2,j)/g(3,j)
g(6,j)=input(‘Valor del factor de fricción ‘)
g(7,j)=g(6,j)*g(5,j)
end

v=input(‘Valor de la viscosidad [lb/ft.s]= ‘);
%v=8.04e-4
d=input(‘Valor de la densidad [lb/ft^3]= ‘);
%d=62.371
Pa=input(‘Valor de la Presión en el punto a [lb/ft^2]= ‘);
%Pa=4320
Pb=input(‘Valor de la Presión en el punto b [lb/ft^2]= ‘);
%Pb=1440
Za=input(‘Valor de la altura a [ft]= ‘);
%Za=200
Zb=input(‘Valor de la altura b [ft]= ‘);
%Zb=50
e=input(‘Valor de e= ‘);
%e=0.00015
L1=input(‘Valor del tramo recto de la tubería 1 [ft]= ‘);
%L1=3000
L2=input(‘Valor del tramo recto de la tubería 2 [ft]= ‘);
%L2=2000
L3=input(‘Valor del tramo recto de la tubería 3 [ft]= ‘);
%L3=200
L4=input(‘Valor del tramo recto de la tubería 4 [ft]= ‘);
%L4=50
gc=32.2;
g1=32.2;

disp(‘***Procedimiento para calcula Q en la primera tubería***’)

hfs1=(Za-Zb)*(g1/gc)+(Pa-Pb)/d
Nka=(g(3,1)*d/v)*sqrt((2*gc*g(3,1)*hfs1)/L1)
ed=e/g(3,1)
x=1/sqrt(g(6,1));
y=sqrt((2*gc*g(3,1)*hfs1)/L1);
disp(‘Las unidades de la velocidad son [ft/s]’)
u=x*y
disp(‘Las unidades del gasto son [ft^3/s]’)
Q1=u*g(4,1)

disp(‘***Procedimiento para calcula Q en la segunda tubería***’)

Nka=(g(3,2)*d/v)*sqrt((2*gc*g(3,2)*hfs1)/L2)
ed=e/g(3,2)
x=1/sqrt(g(6,2))
y=sqrt((2*gc*g(3,2)*hfs1)/L1)
disp(‘Las unidades de la velocidad son [ft/s]’)
u=x*y
disp(‘Las unidades del gasto son [ft^3/s]’)
Q2=u*g(4,2)
Qa=Q1+Q2