i am designing a 500m pipeline, but i'm facing a problem with my code, it can't seem to plot a graph. i know it's a bit long but this is the only way to post the question. i have a lot of columns due to the range of velocities used.

Question

omar mahallawy 2018-12-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/434557-i-am-designing-a-500m-pipeline-but-i-m-facing-a-problem-with-my-code-it-can-t-seem-to-plot-a-grap

评论： omar mahallawy 2018-12-8

clear,clc,clf
format long
Q=0.03;
rouh=847;
g=9.81;
Gamma=(rouh.*g);
Zs=0;
Ze=0;
L=500;
Ks=(6*10.^-5);
viscosity=(1.1*10.^-6);
V=1.34:0.5:4.5;
for I=1
D=sqrt(((4.*Q))./(V.*pi))
end
Re=(V.*D./viscosity);
Er=(1.*10.^-8);
Fo=0.01;
for I=0:1:10^6;
    Fn=(1./(-4.*log10((Ks./(3.71.*D))+((1.26)./(Re.*sqrt(Fo)))))).^2;
    E=abs((Fn-Fo)/Fn);
    if  E<=Er;
     display(Fn),display(E),break,
      end
if E>Er;
   Fo=Fn;
end
end
Hloss=((Fn.*L.*Q.^2)./(12.*D.^5))
Hpump=Hloss
%calculating Hloss every 100M
for HLC=1:100:500
    display(HLC)
    r0=HLC
HlossN=((Fn.*HLC.*Q.^2)./(12.*D.^5))
TEL=(Hpump-HlossN+Zs)
r1 = TEL(:,1)
HGL=(TEL-((V.^2)./2.*g))
r2 = HGL(:,1)
plot(r0,r1)
hold on
plot(r0,r2)
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

madhan ravi 2018-12-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/434557-i-am-designing-a-500m-pipeline-but-i-m-facing-a-problem-with-my-code-it-can-t-seem-to-plot-a-grap#answer_351178

编辑：madhan ravi 2018-12-8

在 MATLAB Online 中打开

clear,clc,clf
format long
Q=0.03;
rouh=847;
g=9.81;
Gamma=(rouh.*g);
Zs=0;
Ze=0;
L=500;
Ks=(6*10.^-5);
viscosity=(1.1*10.^-6);
V=1.34:0.5:4.5;
for I=1
D=sqrt(((4.*Q))./(V.*pi))
end
Re=(V.*D./viscosity);
Er=(1.*10.^-8);
Fo=0.01;
for I=0:1:10^6;
    Fn=(1./(-4.*log10((Ks./(3.71.*D))+((1.26)./(Re.*sqrt(Fo)))))).^2;
    E=abs((Fn-Fo)/Fn);
    if  E<=Er;
     display(Fn),display(E),break,
      end
if E>Er;
   Fo=Fn;
end
end
Hloss=((Fn.*L.*Q.^2)./(12.*D.^5))
Hpump=Hloss
%calculating Hloss every 83.33M
 HLC=linspace(0,500,7);
    r0=HLC;
HlossN=((Fn.*HLC.*Q.^2)./(12.*D.^5));
TEL=(Hpump-HlossN+Zs);
r1 = TEL;
HGL=(TEL-((V.^2)./2.*g));
r2 = HGL;
plot(r0,r1,'-ok')
hold on
plot(r0,r2,'-om')