《《輸油管道設(shè)計(jì)與管理》課程設(shè)計(jì)水力計(jì)算程序4頁》由會(huì)員分享,可在線閱讀,更多相關(guān)《《輸油管道設(shè)計(jì)與管理》課程設(shè)計(jì)水力計(jì)算程序4頁(4頁珍藏版)》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
1、#include#include#define Q 1002.9 /流量,m3/h;#define D 0.559 /管外徑,m; #define d 0.5384 /管內(nèi)徑,m; #define v 9.028e-6 /運(yùn)動(dòng)粘度,m2/s; #define L 1953 /管道全長,Km;#define Zz 1530 /終點(diǎn)海拔高度Zz=1530m;#define ZQ 800 /起點(diǎn)海拔高度ZQ=0m;#define rowt 853.5 /油品在計(jì)算溫度下的密度,Kg/m3;#define Hsq 0.3e6/(9.81*rowt) /首站進(jìn)站壓力,m水柱;#define Hsz 0.
2、2e6/(9.81*rowt) /末站剩余壓頭,m水液柱;#define hm 10 /站內(nèi)損失,m;main() double Re, Re1, Re2, E, Qs; double lamda, hl; double A, N; double ai; double H, Hc,V; Qs=Q/3600; V=4*Qs/(3.14*d*d); printf(n 流速V=%f(m/s)nn,V); printf( 流量Qs=%f(m3/s)nn,Qs); Re=4*1002.9/(3.14*d*3600*v); printf( 雷諾數(shù)Re=%0.2Enn,Re); E=2*0.1/(d*100
3、0); printf( 相對(duì)當(dāng)量粗糙度E=%0.2Enn,E); Re1=59.5/pow(E,8.0/7); printf( 第一邊界雷諾數(shù)Re1=%0.2Enn,Re1); Re2=(665-765*log10(E)/E; printf( 第二邊界雷諾數(shù)Re2=%0.2Enn,Re2); /*判斷流態(tài),選用對(duì)應(yīng)摩阻計(jì)算公式;*/*如果流態(tài)在層流區(qū)*/if(Re2000) printf( /*/n); printf( Re2000 層流區(qū) n); printf( 層流區(qū)水力摩阻系數(shù)lamda=64/Re;nn); lamda=64/Re; printf( 水力摩阻系數(shù) lamda=%0.4f
4、nn,lamda); printf( 選用達(dá)西公式計(jì)算沿程阻力損失hl(m)n); hl=(lamda*(L*1000)*(V*V)/(d*2*9.81); printf( 沿程阻力損失hl hl=%0.2Enn,hl); printf( 選用列賓宗公式計(jì)算沿程阻力損失hl(m)n); hl=4.15*Qs*v*(L*1000)/pow(d,4); printf( 沿程阻力損失hl hl=%0.2Enn,hl); ai=4.15*Qs*v/pow(d,4); printf( 水力坡降ai=%0.5fnn,ai); /*如果流態(tài)在水力光滑區(qū)*/ if(3000Re&ReRe1) printf(
5、/*/n); printf( 3000ReRe1 水力光滑區(qū) n); if(Re10e5) printf( 3000Re10e5n); printf( 水力摩阻選用lamda=0.3164*pow(Re,-0.25);nn);lamda=0.3164*pow(Re,-0.25); printf( 水力摩阻系數(shù) lamda=%0.4fnn,lamda);printf( 選用達(dá)西公式計(jì)算沿程阻力損失hl(m)n); hl=(lamda*(L*1000)*(V*V)/(d*2*9.81); printf( 沿程阻力損失hl hl=%0.2Enn,hl); printf( 選用列賓宗公式計(jì)算沿程阻力損
6、失hl(m)n); hl=(0.0246*pow(Qs,1.75)*pow(v,0.25)*(L*1000)/pow(d,4.75); printf( 沿程阻力損失hl hl=%0.2Enn,hl);ai=0.0246*pow(Qs,1.75)*pow(v,0.25)/pow(d,4.75); printf( 水力坡降ai=%fnn,ai); else lamda=(1/(1.81*log10(Re)-1.53)*(1/(1.81*log10(Re)-1.53); printf( 水力摩阻系數(shù) lamda=%0.4fnn,lamda); printf( 選用達(dá)西公式計(jì)算沿程阻力損失hl(m)n
7、); hl=(lamda*(L*1000)*(V*V)/(d*2*9.81); printf( 沿程阻力損失hl hl=%0.2Enn,hl); printf( 選用列賓宗公式計(jì)算沿程阻力損失hl(m)n); hl=(0.0246*pow(Qs,1.75)*pow(v,0.25)*(L*1000)/pow(d,4.75); printf( 沿程阻力損失hl hl=%0.2Enn,hl); ai=0.0246*pow(Qs,1.75)*pow(v,0.25)/pow(d,4.75); printf( 水力坡降ai=%fnn,ai); /if(3000Re&ReRe1); /*如果流態(tài)在混合摩擦區(qū)
8、*/if(Re1Re&ReRe2) printf( /*/n); printf( Re1ReRe2 混合摩擦區(qū) n); printf( 水力摩阻選用lamda=0.11*pow(e/(d*1000)+68/Re),0.25);nn); lamda=0.11*pow(E/(d*1000)+98/Re),0.25); printf( 水力摩阻系數(shù) lamda=%0.4fnn,lamda); printf( 選用達(dá)西公式計(jì)算沿程阻力損失hl(m)n); hl=(lamda*(L*1000)*(V*V)/(d*2*9.81); printf( 沿程阻力損失hl hl=%0.2Enn,hl); prin
9、tf( 選用列賓宗公式計(jì)算沿程阻力損失hl(m)n); A=pow(10,(0.127*log10(E/(d*1000)-0.627); hl=(0.0802*A*pow(Qs,1.877)*pow(v,0.123)*(L*1000)/pow(d,4.877); printf( 沿程阻力損失hl hl=%0.2Enn,hl); ai=0.0802*A*pow(Qs,1.877)*pow(v,0.123)/pow(d,4.877); printf( 水力坡降ai=%fnn,ai); /if(Re1Re&ReRe2) printf( /*/n); printf( ReRe2 粗糙區(qū) n); pri
10、ntf( 水力摩阻選用lamda=1/pow(1.74-2*log10(E),2);nn); lamda=1/pow(1.74-2*log10(E),2); printf( 水力摩阻系數(shù) lamda=%0.4fnn,lamda); printf( 選用達(dá)西公式計(jì)算沿程阻力損失hl(m)n); hl=(lamda*(L*1000)*(V*V)/(d*2*9.81); printf( 沿程阻力損失hl hl=%0.2Enn,hl); printf( 選用列賓宗公式計(jì)算沿程阻力損失hl(m)n); N=0.11*pow(E/(d*1000),0.25); hl=(0.0826*N*(Qs*Qs)*(
11、L*1000)/pow(d,5); printf( 沿程阻力損失hl hl=%0.2Enn,hl); ai=0.0826*N*Q*Q/pow(d,5); printf( 水力坡降ai=%fnn,ai); /*管線全線所消耗的壓力能頭*/ printf( /*/n); printf( 管線全線所消耗的壓力能頭n);printf( /*/nn);H=hl+(Zz-ZQ)+Hsz; printf( 末站剩余壓頭Hsz=%0.2f(m)nn,Hsz);printf( 管線全線所消耗的壓力能頭H=%0.2f(m)nn,H); /* *泵站數(shù)量計(jì)算 */ printf( /*/n); printf( 所需
12、泵站數(shù)量n); printf( /*/nn); Hc=2*(386.43-0.000176*pow(Q,1.75); printf( 泵站的揚(yáng)程Hc(m) Hc=%0.2fnn,Hc); H=ai*L*1000+(Zz-ZQ)+Hsz; printf( 全線所需壓頭H(m)H=%0.2fnn,H); N=H/(Hc-hm); printf( 全線布站個(gè)數(shù)N N=%0.1fnn,N); N=ceil(H/(Hc-hm); printf( 向上取整站數(shù)N N=%0.0fnn,ceil(H/(Hc-hm); /*工藝方案經(jīng)濟(jì)計(jì)算,確定最優(yōu)方案*/ double Yguan, Yzhan, Yyun,
13、 YY, P; printf( /*/n); printf( 工藝方案經(jīng)濟(jì)計(jì)算,確定最優(yōu)方案n); printf( /*/nn); Yguan=2.4*(7.85*3.14*(D*D-d*d)*L*1000/4)*5603; printf( 管線投資費(fèi)用 Yguan=%0.1f(元)nn,Yguan); printf( 管線投資費(fèi)用 Yguan=%0.3f(億元)nn,Yguan/1e8); Yzhan=1.8*(N*3000e4); printf( 站場(chǎng)投資費(fèi)用 Yzhan=%0.1f(元)nn,Yzhan); printf( 站場(chǎng)投資費(fèi)用 Yzhan=%0.3f(億元)nn,Yzhan/1e8); /*計(jì)算輸油泵功率*/ P=1.1*rowt*Qs*Hc/(102*0.807); printf( 輸油泵站功率 P=%0.2E(Kw)nn,P); Yyun=1.0*0.55*N*P*15*350*24; printf( 運(yùn)營費(fèi)用 Yyun=%0.1f(元)nn,Yyun); printf( 運(yùn)營費(fèi)用 Yyun=%0.3f(億元)nn,Yyun/1e8); YY=Yguan+Yzhan+Yyun; printf( 方案總的經(jīng)濟(jì)費(fèi)用YY=%0.1f(元)nn,YY); printf( 方案總的經(jīng)濟(jì)費(fèi)用YY=%0.3f(億元)nn,YY/1e8);/main()