油氣儲運(yùn)網(wǎng)

標(biāo)題: 求解下面的編程求助紅色區(qū)域是 [打印本頁]

作者: sun7854 時間: 2017-3-12 17:46
標(biāo)題: 求解下面的編程求助紅色區(qū)域是

d1=0.001,d2=1.5;

while(1)

{

d3=(d1+d2)/2.0;

Re=1.534*Q*del*1000/(d3*u);

R1=0.001;

R2=1.0;

while(1)

{

R3=(R1+R2)/2.0;if(fr(k,d3,Re,R3)>1e-6)R1=R3; if(fr(k,d3,Re,R3)<-1e-6)R2=R3;if(fabs(fr(k,d3,Re,R3))<1e-6)break;

}

R=R3; if(fd(R,z,d3)>1e-6)d1=d3;

if(fd(R,z,d3)<-1e-6) d2=d3;if(fabs(d3-d1)<1e-6&&fabs(d3-d2)<1e-6)break;

}

d=d3; printf("摩阻系數(shù)為R=%e,管徑為d=%em\n\n\n",R,d); printf("雷諾數(shù)為Re=%e\n\n\n",Re);

}

[qq]785430264[/qq]

作者: sun7854 時間: 2017-3-12 17:47
截取了自己不懂得地方，發(fā)布出來全部的編程點(diǎn)我QQ 謝謝前輩大神們

作者: hqbtc 時間: 2017-3-13 17:16
兩個while循環(huán)是死循環(huán)
因?yàn)闂l件為1永遠(yuǎn)為真
if語句的break條件滿足時，跳出循環(huán)
具體結(jié)果自己分析

作者: sun7854 時間: 2017-3-13 23:28

hqbtc 發(fā)表于 2017-3-13 17:16
兩個while循環(huán)是死循環(huán)
因?yàn)闂l件為1永遠(yuǎn)為真
if語句的break條件滿足時，跳出循環(huán)

程序沒問題，調(diào)用子函數(shù)算的fr和fd。就是上面的管徑d1和d2的賦值0.001和1.5，這是輸氣管道管徑的編程。發(fā)不了全部的只發(fā)了我不懂的這個

歡迎光臨油氣儲運(yùn)網(wǎng) (http://black-hills-tours.com/)