流体力学matlab例题,流体力学简单计算MATLAB程式
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
流体⼒学matlab例题,流体⼒学简单计算MATLAB程式
《流体⼒学简单计算MATLAB程式》由会员分享,可在线阅读,更多相关《流体⼒学简单计算MATLAB程式(2页珍藏版)》请在⼈⼈⽂库⽹上搜索。
1、⽤matlab进⾏编程计算第⼀问:
z=30;p1=50*9.8*104;p2=2*9.8*104;jdc=0.00015;gama=9800;d=0.257;L=50000;mu=6*10(-6);hf=z+(p1-
p2)/(0.86*gama)xdc=2*jdc/d;beta=4.15;m=1;Q=(hf*d(5-m)/(beta*mum*L)(1/(2-
m);v=4*Q/(pi*d2);Re=v*d/mu;Re1=59.7/xdc(8*xdc/7);Re2=(665-765*log(xdc)/xdc;i=hf/L;if Re3000Q=Q;elseif
3000ReRe1m=0.25;be。
2、ta=0.0246;Q=(hf*d(5-m)/(beta*mum*L)(1/(2-m);v=4*Q/(pi*d2);Re=v*d/mu;elseif
Re1ReRe2m=0.123;A=10(0.127*log(jdc/d)-0.627);beta=0.0802*A;Q=(hf*d(5-m)/(beta*mum*L)(1/(2-
m);v=4*Q/(pi*d2);Re=v*d/mu;else m=0;langda=1/(2*log(3.7*d/jdc)2;beta=0.0816*langda;Q=(hf*d(5-m)/(beta*mum*L)(1/(2-m);v=4*Q/(pi*d2);Re=v。
3、*d/mu;endiQRevhf =588.1395i =0.0118Q =0.0915Re =7.5526e+004v =1.7632利⽤IF 语句对四种流态⼀⼀进⾏试算,最终的结果⽔⼒损失为588.1395m , ⽔⼒坡降0.0118第⼆问:将其中10km换成直径305mm的管⼦
z=30;p1=50*9.8*104;p2=2*9.8*104;jdc=0.00015;gama=9800;mu=6*10(-6);d=0.257 0.305;L=40000
10000;x=zeros(4,1);Q1=0.0915;hf=z+(p1-p2)/(0.86*gama);xdc=2*jdc./d;bet。
4、a=4.15;m=1;Q=(hf/(beta*mum*(L(1)/d(1)(5-m)+L(2)/d(2)(5-m)(1/(2-
m);v(1)=4*Q/(pi*d(1)2);v(2)=4*Q/(pi*d(2)2);Re(1)=v(1)*d(1)/mu;Re(2)=v(2)*d(2)/mu;Re1(1)=59.7/xdc(1)
(8/7);Re1(2)=59.7/xdc(2)(8/7);Re2(1)=(665-765*log(xdc(1)/xdc(1);Re2(2)=(665-765*log(xdc(2)/xdc(2);if Re(1)3000 & Re(2)3000 Q=Q;elseif 3000。
5、Re(1)Re1 & 3000Re(2)Re1 m=0.25;beta=0.0246;Q=(hf/(beta*mum*(L(1)/d(1)(5-m)+L(2)/d(2)(5-m)(1/(2-
m);v(1)=4*Q/(pi*d(1)2);v(2)=4*Q/(pi*d(2)2);Re(1)=v(1)*d(1)/mu;Re(2)=v(2)*d(2)/mu;elseif Re1Re(1)Re2 &
Re1Re(2)Re2m=0.123;A=10(0.127*log(jdc/d)-0.627);beta=0.0802*A;Q=(hf/(beta*mum*(L(1)/d(1)(5-m)+L(2)/d(2)。
6、(5-m)(1/(2-m);v(1)=4*Q/(pi*d(1)2);v(2)=4*Q/(pi*d(2)2);Re(1)=v(1)*d(1)/mu;Re(2)=v(2)*d(2)/mu;else
m=0;langda=1/(2*log(3.7*d/jdc)2;beta=0.0816*langda;Q=(hf/(beta*mum*(L(1)/d(1)(5-m)+L(2)/d(2)(5-m)(1/(2-
m);v(1)=4*Q/(pi*d(1)2);v(2)=4*Q/(pi*d(2)2);Re(1)=v(1)*d(1)/mu;Re(2)=v(2)*d(2)/mu;endQ Revq=Q-Q1baifenbi=q/Q1Q
=0.0978Re =1.0e+004 *8.0795 6.8080v =1.8863 1.3393q =0.0063baifenbi =0.0694先假设流态均为层流,再进⾏试算,最终可得输量可提⾼6.94个百分点。