Problema 26

Código de programación Matlab

% PROBLEMA 26

clc
clear
format compact

% Calcule el NPSH disponible para el siguiente sistema si el NPSHreal=5ft
Di=input(‘Valor del diámetro interno de la tubería ‘);
%Di=0.1722
S=input(‘Valor del área de la tubería ‘);
%S=.0233
d=input(‘Valor de la densidad de la sustancia ‘);
%d=54.1
A=input(‘Valor de la constante de Antoine para calcular la presión de vapor ‘);
%A=15.9008
B=input(‘Valor de la constante de Antoine para calcular la presión de vapor ‘);
%B=2788.51
C=input(‘Valor de la constante de Antoine para calcular la presión de vapor ‘);
%C=52.35
T=input(‘Valor de la temperatura ‘);
%T=300
Q=input(‘Valor del flujo volumétrico del sistema ‘);
%Q=.3343
Z1=input(‘Valor de la altura 1 ‘);
%Z1=15
Pabs=input(‘Valor de la presión absoluta ‘);
%Pabs=2635.2
g1=32.2;
gc=32.2;
n=input(‘Cuantos accesorios tiene el sistema ‘)
g=zeros(n,4)
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 ‘);
fprintf(‘\n Accesorio Cantidad L/D Kt\n’)
g
end
ktot=g(:,4);
kt1=sum(ktot)
e=2.71828182818281828182;
disp(‘Las unidades de la presión de vapor son [lbf/ft^2]’)
Po=e^(A-(B/(T-C)))*(2116.2/760)
disp(‘Las unidades de la velocidad son [ft/seg]’)
U=Q/S
disp(‘Las unidades del Hfs son [lbf ft/lb]’)
Hfs=(kt1*U*U)/(2*gc)
disp(‘Las unidades del NPSH son [ft]’)
NPSH=(Z1*(g1/gc))+((Pabs-Po)/d)-Hfs