| 5 | 1/1 | 返回列表 |
| 查看: 2058 | 回復(fù): 4 | ||
yuan020545金蟲(chóng) (初入文壇)
|
[求助]
關(guān)于c++中計(jì)算結(jié)果出現(xiàn)-1.#IND的問(wèn)題
|
|
流體力學(xué)中用c++編程計(jì)算串聯(lián)管路的流量q,完整的程序如下,但是計(jì)算結(jié)果出現(xiàn)-1.#IND,不知道什么情況,請(qǐng)大神們幫看看吧,應(yīng)該怎么寫(xiě)才對(duì)? #include #include #define pi 3.1415926 #define g 9.8 float CL(float dlt,float d,float re,int *id) { float lmd1,temp,lmd,Rec1,Rec2,Rec3; Rec1=2320; Rec2=80*d/dlt; Rec3=4160*pow(d/2.0/dlt,0.85); if(re if(re<1e5){lmd=0.3164/pow(re,0.25);*id=2;} else if(re<3e6){lmd=0.0032+0.221*pow(re,-0.237);*id=3;} } else if(re lmd=0.01;*id=4; while(fabs(lmd-lmd1)>=1e-6) { lmd1=lmd; temp=log(dlt/(3.7*d)+2.51/re/sqrt(lmd1))/log(10.0); lmd=pow(1.0/(-2.0)/temp,2); } } else if(re>=Rec3) { lmd=pow((1.74+2.0*log(d/(2.0*dlt))/log(10.0)),-2); *id=2; } return(lmd); } float CL(float,float,float,int *); void main() { int ID1=0; int ID2=0; float DLT1,LMD1,LMD11,Re1,D1,L1,V1; float DLT2,LMD2,LMD22,Re2,D2,L2,V2; float Q,NU,H,KCE; float DK,lll; printf("please input:h,l1,d1,l2,d2,nu,dlt1,dlt2,kce:\n" ;scanf("%f%f%f%f%f%f%f%f%f",&H,&L1,&D1,&L2,&D2,&NU,&DLT1,&DLT2,&KCE); LMD11=0.025;LMD22=0.015; loop: DK=D1*D1/(D2*D2); lll=KCE+LMD11*L1/D1+(1-DK)*(1-DK); lll=lll+LMD22*L2/D2*DK*DK+DK*DK; V1=sqrt(2*g*H/lll); V2=DK*V1; Re1=V1*D1/NU; Re2=V2*D2/NU; LMD1=CL(DLT1,D1,Re1,&ID1); LMD2=CL(DLT2,D2,Re2,&ID2); if((fabs(LMD1-LMD11)>=5e-6)||(fabs(LMD2-LMD22)>=5e-6)) { LMD11=LMD1; LMD22=LMD2; goto loop; } Q=pi*D1*D1*V1/4; printf("ID=%d\nv1=%7.4fm/s\nre1=%7.2f\nlmd1=%7.4f\nID2=%d\nv2=%7.4fm/s\nre2=%7.2f\nlmd2=%7.4f\nq=%7.4fm3/s\n",ID1,V1,Re1,LMD1,ID2,V2,Re2,LMD2,Q); } 運(yùn)行程序: input:H,D1,D2,DLT1,DLT2,L1,L2,NU 6.,0.6,0.9,0.0015,0.0003,300.,240.,1.E-6 正確的輸出如下: ID1=5 V1=2.8713m/s RE1=1722771.00 LMD1=0.0249 ID2=4 V2=1.2761m/s RE2=1148514.00 LMD2=0.0159 Q=0.8118m3/s [ Last edited by yuan020545 on 2011-5-15 at 22:35 ] |
木蟲(chóng) (小有名氣)

金蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (初入文壇)
金蟲(chóng) (著名寫(xiě)手)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚(yú) 2026-03-18 | 7/350 |
|
|
[考研] 354求調(diào)劑 +5 | Tyoumou 2026-03-18 | 8/400 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸! 2026-03-16 | 8/400 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變,給編輯又寫(xiě)郵件了,沒(méi)回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|