| 5 | 1/1 | 返回列表 |
| 查看: 2059 | 回復(fù): 4 | ||
yuan020545金蟲 (初入文壇)
|
[求助]
關(guān)于c++中計(jì)算結(jié)果出現(xiàn)-1.#IND的問題
|
|
流體力學(xué)中用c++編程計(jì)算串聯(lián)管路的流量q,完整的程序如下,但是計(jì)算結(jié)果出現(xiàn)-1.#IND,不知道什么情況,請大神們幫看看吧,應(yīng)該怎么寫才對? #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 ] |
木蟲 (小有名氣)

金蟲 (著名寫手)
新蟲 (初入文壇)
金蟲 (著名寫手)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)工程321分求調(diào)劑 +18 | 大米飯! 2026-03-15 | 22/1100 |
|
|
[考研] 326求調(diào)劑 +5 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 8/400 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗(yàn)) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|