Código de programación Matlab
% PROBLEMA 13
clc
clear
format compact;
% ¿Qué diámetro de tubería nueva de fierro fundido se necesita para
% transportar agua a 25°C a lo largo de 1 km con una pérdida de carga de 2m?
d=input(‘Cual es la densidad del fluido ‘);
%d=997
Q=input(‘Cual es la carga [m^3/s] ‘);
%Q=0.02298157
Lt=input(‘Valor del tramo de la tubería ‘);
%Lt=1000
v=input(‘Valor de la viscosidad de la sustancia ‘);
%v=8.8e-4
disp(‘Este problema es iterativo por lo cual usaremos sólo 5 iteracciones para llegar al resultado’)
n=input(‘Valor de iteracciones’)
h=input(‘Cual es la pérdida de la carga’);
%h=2
a=9.81
q2=Q^2
p=pi^2
m=(4*d*Q)/(pi*v)
disp(‘El valor del diámetro esta dado en [m]’)
g=zeros (n,5)
g(1,1)=input(‘Valor propuesto para el diámetro ‘);
g(1,2)=m/g(1,1)
g(1,3)=input(‘Valor de la E/D ‘);
g(1,4)=input(‘Valor del factor de fricción ‘);
g(1,5)=((8*q2*Lt*g(1,4))/(p*h*a))^(1/5)
g(2,1)=g(1,5)
g(2,2)=m/g(2,1)
g(2,3)=input(‘Valor de la E/D ‘);
g(2,4)=input(‘Valor del factor de fricción ‘);
g(2,5)=((8*q2*Lt*g(2,4))/(p*h*a))^(1/5)
g(3,1)=g(2,5)
g(3,2)=m/g(3,1)
g(3,3)=input(‘Valor de la E/D ‘);
g(3,4)=input(‘Valor del factor de fricción ‘);
g(3,5)=((8*q2*Lt*g(3,4))/(p*h*a))^(1/5)
Flujo de Fluidos 