Código de programación Matlab
% PROBLEMA 18
clc
clear
format compact
% Encuentre el caudal volumétrico de cada uno de los ramales. El fluido es
% benceno a 68°F. La caída de presión es igual a cero.
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=55
dp=input(‘Valor de la caída de presión’)
%dp=0
dz=input(‘Valor de la diferencia de alturas’)
%dz=23
e=input(‘Valor de e= ‘);
%e=0.00015
L1=input(‘Valor del tramo recto de la tubería 1 [ft]= ‘);
%L1=30
L2=input(‘Valor del tramo recto de la tubería 2 [ft]= ‘);
%L2=30
gc=32.2;
g1=32.2;
disp(‘Procedimiento para calcula Q en la primera tubería’)
hfs1=(dz)*(g1/gc)+(dp)/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
ir al siguiente problema