《输油管道设计与管理》课程设计水力计算程序.docVIP

#includestdio.h #includemath.h #define Q 1002.9 //流量,m^3/h; #define D 0.559 //管外径,m; #define d 0.5384 //管内径,m; #define v 9.028e-6 //运动粘度,m^2/s; #define L 1953 //管道全长,Km; #define Zz 1530 //终点海拔高度Zz=1530m; #define ZQ 800 //起点海拔高度ZQ=0m; #define rowt 853.5 //油品在计算温度下的密度,Kg/m^3; #define Hsq 0.3e6/(9.81*rowt) //首站进站压力,m水柱; #define Hsz 0.2e6/(9.81*rowt) //末站剩余压头,m水液柱; #define hm 10 //站内损失,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)\n\n,V); printf( 流量Qs=%f(m^3/s)\n\n,Qs); Re=4*1002.9/(3.14*d*3600*v); printf( 雷诺数Re=%0.2E\n\n,Re); E=2*0.1/(d*1000); printf( 相对当量粗糙度E=%0.2E\n\n,E); Re1=59.5/pow(E,8.0/7); printf( 第一边界雷诺数Re1=%0.2E\n\n,Re1); Re2=(665-765*log10(E))/E; printf( 第二边界雷诺数Re2=%0.2E\n\n,Re2); /******************************************* ***判断流态,选用对应摩阻计算公式; ********************************************/ /************如果流态在层流区**************/ if(Re2000) { printf( /***********************************/\n); printf( Re2000 层流区 \n); printf( 层流区水力摩阻系数lamda=64/Re;\n\n); lamda=64/Re; printf( 水力摩阻系数 lamda=%0.4f\n\n,lamda); printf( 选用达西公式计算沿程阻力损失hl(m)\n); hl=(lamda*(L*1000)*(V*V))/(d*2*9.81); printf( 沿程阻力损失hl hl=%0.2E\n\n,hl); printf( 选用列宾宗公式计算沿程阻力损失hl(m)\n); hl=4.15*Qs*v*(L*1000)/pow(d,4); printf( 沿程阻力损失hl hl=%0.2E\n\n,hl); ai=4.15*Qs*v/pow(d,4); printf( 水力坡降ai=%0.5f\n\n,ai); } /************如果流态在水力光滑区**************/ if(3000ReReRe1) { printf( /***********************************/\n); printf( 3000ReRe1 水力光滑区 \n); if(Re10e5) {printf( 3000Re10e5\n); printf( 水力摩阻选用lamda=0.3164*pow(Re,-0.25);\n\n); lamda=0.3164*pow(Re,-0.25); printf( 水力摩阻系数 lamda=%0.4f\n\n,lamda); printf( 选用达西公式计算沿程阻力损失hl(m)\n); hl=(lamda*(L*1