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

《《输油管道设计与管理》课程设计水力计算程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、#include#include#defineQ1002.9//流量,m^3/h;#defineD0.559//管外径,m;#defined0.5384//管内径,m;#definev9.028e-6//运动粘度,m^2/s;#defineL1953//管道全长,Km;#defineZz1530//终点海拔高度Zz=1530m;#defineZQ800//起点海拔高度ZQ=0m;#definerowt853.5//油品在计算温度下的密度,Kg/m^3;#defineHsq0.3e6/(9.81*rowt)//

2、首站进站压力,m水柱;#defineHsz0.2e6/(9.81*rowt)//末站剩余压头,m水液柱;#definehm10//站内损失,m;main(){doubleRe,Re1,Re2,E,Qs;doublelamda,hl;doubleA,N;doubleai;doubleH,Hc,V;Qs=Q/3600;V=4*Qs/(3.14*d*d);printf("流速V=%f(m/s)",V);printf("流量Qs=%f(m^3/s)",Qs);Re=4*1002.9/(3.14*d*3600*v);printf("

3、雷诺数Re=%0.2E",Re);E=2*0.1/(d*1000);printf("相对当量粗糙度E=%0.2E",E);Re1=59.5/pow(E,8.0/7);printf("第一边界雷诺数Re1=%0.2E",Re1);Re2=(665-765*log10(E))/E;printf("第二边界雷诺数Re2=%0.2E",Re2);/**********************************************判断流态,选用对应摩阻计算公式;*************************

4、*******************//************如果流态在层流区**************/if(Re<2000){printf("/***********************************/");printf("Re<2000层流区");printf("层流区水力摩阻系数lamda=64/Re;");lamda=64/Re;printf("水力摩阻系数lamda=%0.4f",lamda);printf("选用达西公式计算沿程阻力损失hl(m)");hl=(lamda*(L*10

5、00)*(V*V))/(d*2*9.81);printf("沿程阻力损失hlhl=%0.2E",hl);printf("选用列宾宗公式计算沿程阻力损失hl(m)");hl=4.15*Qs*v*(L*1000)/pow(d,4);printf("沿程阻力损失hlhl=%0.2E",hl);ai=4.15*Qs*v/pow(d,4);printf("水力坡降ai=%0.5f",ai);}/************如果流态在水力光滑区**************/if(3000

6、***********************************/");printf("3000

7、V))/(d*2*9.81);printf("沿程阻力损失hlhl=%0.2E",hl);printf("选用列宾宗公式计算沿程阻力损失hl(m)");hl=(0.0246*pow(Qs,1.75)*pow(v,0.25)*(L*1000))/pow(d,4.75);printf("沿程阻力损失hlhl=%0.2E",hl);ai=0.0246*pow(Qs,1.75)*pow(v,0.25)/pow(d,4.75);printf("水力坡降ai=%f",ai);}else{lamda=(1/(1.81*log10

8、(Re)-1.53))*(1/(1.81*log10(Re)-1.53));printf("水力摩阻系数lamda=%0.4f",lamda);printf