Problema 27

Código de programación Matlab

% PROBLEMA 27

clc
clear
format compact; format short g;

% Se va a vaciar un carro que contiene 10 000 gal de bencenos a 80°F en 3
% hr, el sistema de tuberías de la planta es como se indica adelante. Se
% tiene disponible para el bombeo una bomba centrifuga y los siguientes
% datos. Calcule
% a)Cuanto tiempo tardara en vaciarse el carro tanque
% b)La cantidad de trabajo y tiempo real
% c)El NPSH disponible

Di=input(‘Valor del diámetro interno ‘);
%Di=.2557
S=input(‘Valor del área de la tubería ‘);
%S=0.0513
Q=input(‘Valor del gasto o flujo ‘);
%Q=0.1238
Z1=input(‘Valor de la altura 1 ‘);
%Z1=4.5
Z2=input(‘Valor de la altura 2 ‘);
%Z2=50
P1=input(‘Valor de la presión 1 ‘);
%P1=11.3
P2=input(‘Valor de la presion 2 ‘);
%P2=11.3
gc=32.2;
g1=32.2;
d=input(‘Valor de la densidad ‘),
%d=54.1
v=input(‘Valor de la viscosidad ‘);
%v=4.03e-4
U=input(‘Valor de la U ‘);
%U=0
A=input(‘Valor de la constante A ‘);
%A=15.908
B=input(‘Valor de la constante B ‘);
%B=2788.51
C=input(‘Valor de la constante C ‘);
%C=52.36
T=input(‘Valor de la temperatura ‘);
%T=300
t=input(‘Valor del tiempo en segundos ‘);
%t=10800
n=input(‘Cuantos datos de Q hay en el sistema ‘)
g=zeros(n,4)
for j=1:n

g(j,1)=input(‘Valor de Q[gal/min] ‘);
g(j,2)=input(‘Valor de H [ft] ‘);
g(j,3)=input(‘Valor del porcentaje ‘);
g(j,4)=g(j,1)*2.228e-3
fprintf(‘\n Q[gal/min] H[ft] n% Q[ft^3/s] \n’)

g
end

nn=input(‘Cuantos accesorios tiene 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 Lei ‘);
h(j,5)=h(j,2)*h(j,4)
fprintf(‘\n Accesorio Cantidad L/D Lei Let\n’)
h
end

ktot=h(:,5);
Lt1=sum(ktot)
Lt=145
L=Lt+Lt1

m=input(‘Cuantos datos de Q hay en el sistema ‘)
p=zeros(m,9)
a=(Z2-Z1)*(g1/gc);
b=(P2-P1)/d;
for j=1:n

p(j,1)=input(‘Valor de Q’)
p(j,2)= p(j,1)/S
p(j,3)=(p(j,2)^2)/(2*gc)
p(j,4)=(d*Di*p(j,2))/v
p(j,5)=input(‘Valor del factor de fricción’)
p(j,6)=(p(j,5)*L*p(j,2)^2)/(2*gc*Di)
p(j,7)=(Z2-Z1)
p(j,8)=(P2-P1)
p(j,9)=[a]+[b]+[(p(j,2)-U)/(2*gc)]+p(j,6)

fprintf(‘\n Q[ft^3/s] U U^2/2*gc Nre f
hfs DZ DP -wf \n’)

p
end

x=g(:,4);
y1=g(:,3)
plot(x,y1,’r’);
hold on;
y2=g(:,2);
plot(x,y2,’g’);
grid on;
y3=p(:,9)
plot(x,y3,’K’)

title(‘Curvas del sistema’);
xlabel(‘Gasto de operación del sistema’);
ylabel(‘H y Eficiencia de la bomba’);
axis([.14 .18 45 55])

disp(‘Los valores que se piden a continuación los tomaras de la grafica’)
wfgraf=input(‘Valor de wf ‘);
%wfgraf=50
Qopt=input(‘Valor de Qopt ‘);
%Qopt=0.163
ngraf=input(‘Valor de ngraf ‘);
%ngraf=.47
f=input(‘Valor del factor de fricción ‘);
%f=0.0207
disp(‘El valor de w esta dado en [lb/seg]’)
w=Qopt*d
disp(‘El valor de la potencia esta dado en [HP]’)
Ph=(wfgraf*w)/550
Pr=Ph/ngraf

Po=(exp((A-(B/(T-C)))))*(2116.2/760)
v=t*Q
disp(‘Las unidades del tiempo estan dadas en [s]’)
tiempo=v/Qopt

Patm=P1*144
f1=0.0207
Lt=168.64
u=Qopt/S;
hfs=(Lt*u^2*f)/(2*gc*Di);
disp(‘Las unidades del NPSH estan dadas en [lbfft/lb]’)
NPSH=Z1+((Patm-Po)/d)-hfs