Problema 16

Código de programación Matlab

% PROBLEMA 16

clc
clear
format compact

% Los tubos que se muestran en el sistema siguiente son de acero comercial.
% Si el gasto es de 20 ft^3/s determine la perdida de carga y la
% distribución de flujos de los tubos entre A y B.

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(‘ Diámetro de la tubería en pulgadas ‘) ;
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

%K1=g(5,1)*g(6*1)
% K2=g(5,2)*g(6*2)
% K3=g(5,3)*g(6*3)
%K4=g(5,1)*g(6*4)

Qp=input(‘Valor propuesto de la carga’)
gc=32.2;
Q=input(‘Valor del gasto inicial’)
%Q=20
Hfs1=(g(7,1)*Qp^2)/(2*gc*g(4,1)^2)
Qp2=sqrt((2*gc*g(4,2)^2*Hfs1)/g(7,2))
Qp3=sqrt((2*gc*g(4,3)^2*Hfs1)/g(7,3))
Qpt=Qp+Qp2+Qp3
Q1=(Qp/Qpt)*Q
Q2=(Qp2/Qpt)*Q
Q3=(Qp3/Qpt)*Q
Hfs1=(g(7,1)*Q1^2)/(2*gc*g(4,1)^2)
Hfs2=(g(7,2)*Q2^2)/(2*gc*g(4,2)^2)
Hfs3=(g(7,3)*Q3^2)/(2*gc*g(4,3)^2)
disp(‘Ahora calcularemos el Hfs para las distancias A-N, M-B y A-B’)
L=input(‘valor de la longitud en el punto A-N) ‘)
%L=2000
L2=input(‘valor de la longitud en el punto M-B) ‘)
%L2=4000
Q=input(‘Valor del gasto en la tubería de 24in’ )
%Q=20
k1=g(6,4)*(L/g(3,4))
k2=g(6,4)*(L2/g(3,4))
hfsAN=(k1*Q*Q)/(2*gc*g(4,4)*g(4,4))
hfsMB=(k2*Q*Q)/(2*gc*g(4,4)*g(4,4))
hfsAB=hfsAN+hfsMB+Hfs3