Problema 15

Código de programación Matlab

% PROBLEMA 15

clc
clear
format compact

% Calcúlese el gasto descargado por un tubo de la figura. El fluido es agua
% a 150°F

d=input(‘Cual es la densidad del fluido ‘);
%d=61.188
Ltr=input(‘Valor del tramo de la tubería ‘);
%Ltr=240
v=input(‘Valor de la viscosidad de la sustancia ‘);
%v=2.8896e-4
Di=input(‘Valor del diámetro interno de la tubería ‘);
%Di=.1722
s=input(‘Valor del área de la tubería ‘);
%s=.0233
z1=input(‘Valor de la altura 1’ );
%z1=0
z2=input(‘Valor de la altura 2 ‘);
%z2=260
ed=input(‘Valor de la E/D ‘);
%ed=0.009
disp(‘Para hacer el cálculo de kt introdúcelos siguientes valores’)
n=input(‘Valor de interacciones ‘)
%n=1
h=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.’)

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)=h(j,3)*Di
h(j,5)=h(j,2)*h(j,4)
fprintf(‘\n Accesorio Cantidad L/D LE LET\n’)
h
Lt=h(:,5)
LT=sum(Lt)+Ltr
end

gc=32.2;
a=8/(pi*pi*gc);
b=LT/Di^5;
c=1/Di^4;
y=z2-z1;
x=4*d/(pi*v*Di);
g=ones (3,4)

g(1,1)=input(‘Valor propuesto para Q ‘)
g(1,2)=x*g(1,1)
g(1,3)=input(‘Valor del factor de fricción ‘)
g(1,4)=sqrt(y/(a*g(1,3)*b+c))

g(2,1)=g(1,4)
g(2,2)=x*g(2,1)
g(2,3)=input(‘Valor del factor de fricción ‘)
g(2,4)=sqrt(y/(a*g(2,3)*b+c))

g(3,1)=g(2,4)
g(3,2)=x*g(3,1)
g(3,3)=input(‘Valor del factor de fricción ‘)
g(3,4)=sqrt(y/(a*g(3,3)*b+c))
fprintf(‘\n Qpropuesta Nre f Qcalculada\n’)
g(4,1)=g(3,4)
g(4,2)=x*g(4,1)
g(4,3)=input(‘Valor del factor de fricción ‘)
g(4,4)=sqrt(y/(a*g(4,4)*b+c))
fprintf(‘\n Qpropuesta Nre f Qcalculada\n’)
g