油库上机编程实验报告

1、-1-油库设计与管理实验报告小组成员李伟王文康陈万贵陈实曹增辉赖家鑫学号100102031610010203151001020319100102031410010203171001020318实验名称油库设计与管理指导教师宋晓琴单位西南石油大学-2-输入参数0HilidiLiD0L0Dn0Vv0Z计算常数4/2d4/2D4/20Dgdl2/gDL2/002/gDL1题目来源某油库的自流装卸系统，如图2所示，装油管设置在集油管中部，且两侧鹤管同时装油，两侧鹤管对称布置。油库油料为车用汽油，操作条件下车用汽油的粘度为0.6×10-6m2/s，鹤管采用Φ108×4的钢管，鹤管计算长度为25m，鹤管间距为12.5m，集油管为Φ273×6，输油管为Φ159×5，输油管计算长度为1500m，油库库容为30000m3，单个罐的容量为3000m3，油罐内径为16m，罐壁高度为15.85，罐出油口与装油鹤管出口之间的高差为30m。（1）计算鹤分别管数为20和30，液位高度为3m和10m时所对应的各鹤管流量？（2）操作鹤管数为30，发油量为1500m3，开始发油时液位高度为12m，计算发油时间？。。

2、2程序框图：YESin/2开机)1(1Q=)1(1Q+1tV(1)(1)hf1(1)u(1))1(hf2(1)i=2ton/21Q(i)=1Q(i-1)+2t)(ihf1(i)1Q(i)giV2)(2+)(igiVdl2)(2giV2)1(2+)1(igiVdl2)(2+)1(iDLgiu2)1(2u(i-1)(i)hf2(i)-3-T=2V0/（Q1+Q0）YES3程序代码3.1定义数据和数据输入3.2系统常数3.3定义数组3.4计算过程3.5计算第一根鹤管OptionBase1PrivateSubCommand1_Click()'定义数据和数据输入Dimdd,D,ll,L,D0,ν,εH,εJ,aa,A,C1,C2,C3,ε,Δ,t,S,E,Fdd=Text5D=Text4D0=Text3N=Text1ν=Text2ll=Text80Q=2()1(1Q+………+)2/(1nQ)；0u=2004DQ00hf输出Q1(1)……Q1（n/2）；0Q输入发油量V0算Z1=Z0-V0/F；调用过程再算Q1gV2)1(2+hf1(1)+hf2(1)+………+hf2(n。

3、/2)+0hf0H+Z结束-4-L=Text7L0=Text6S=Text10'系统常数εH=2*0.00015/ddεJ=2*0.00015/Daa=3.14*dd^2/4A=3.14*D^2/4C1=ll/(19.6*dd)C2=L/(19.6*D)C3=L0/(19.6*D0)'定义数组DimQ0(50)DimV0(50)Dimλ0(50)DimRe(50)DimH0(50)AsString*5Dimvv(50,50)DimV(50,50)Dimλ(50,50)DimΛ(50,50)Dimqq(50,50)DimQ(50,50)DimReH(50,50)DimReJ(50,50)DimB(50,50)'计算过程qq(1,1)=0.01Δ=0.000001-5-Fori=1To30'计算第一根鹤管Ifi1Thenqq(i,1)=qq(i-1,1)+(i-1)*ΔEndIfQ(i,1)=qq(i,1)vv(i,1)=qq(i,1)/aaReH(i,1)=vv(i,1)*dd/νIfReH(i,1)2000Thenλ(i,1)=64/ReH(i,1)ElseIfReH(i,1)3000。

4、AndReH(i,1)2000Thenλ(i,1)=0.16/ReH(i,1)^0.25ElseIfReH(i,1)59.7/εH^(8/7)AndReH(i,1)3000Thenλ(i,1)=0.3164/ReH(i,1)^0.25Elseλ(i,1)=(1/(-1.8*Log(6.8/ReH(i,1)+(εH/7.4)^1.11)/Log(10)))^2EndIf'计算第一段集油管V(i,1)=qq(i,1)/AReJ(i,1)=V(i,1)*D/νIfReJ(i,1)2000ThenΛ(i,1)=64/ReJ(i,1)ElseIfReJ(i,1)3000AndReJ(i,1)2000ThenΛ(i,1)=0.16/ReJ(i,1)^0.25ElseIfReJ(i,1)59.7/εJ^(8/7)AndReJ(i,1)3000ThenΛ(i,1)=0.3164/ReJ(i,1)^0.25ElseΛ(i,1)=(1/(-1.8*Log(6.8/ReJ(i,1)+(εJ/7.4)-6-^1.11)/Log(10)))^2EndIf'计算后面鹤管Forj=2ToNqq(i,j)=qq(i,。

5、j-1)Doqq(i,j)=qq(i,j)+0.00001vv(i,j)=qq(i,j)/aaReH(i,j)=vv(i,j)*dd/νIfReH(i,j)2000Thenλ(i,j)=64/ReH(i,j)ElseIfReH(i,j)3000AndReH(i,j)2000Thenλ(i,j)=0.16/ReH(i,j)^0.25ElseIfReH(i,j)59.7/εH^(8/7)AndReH(i,j)3000Thenλ(i,j)=0.3164/ReH(i,j)^0.25Elseλ(i,j)=(1/(-1.8*Log(6.8/ReH(i,j)+(εH/7.4)^1.11)/Log(10)))^2EndIfE=vv(i,j-1)^2/19.6+λ(i,j-1)*C1*vv(i,j-1)^2+Λ(i,j-1)*C2*V(i,j-1)^2F=vv(i,j)^2/19.6+vv(i,j)^2*λ(i,j)*C1IfAbs(E-F)0.05ThenExitDoLoop-7-'计算后面集油管Q(i,j)=Q(i,j-1)+qq(i,j)V(i,j)=Q(i,j)/AReJ(i,j)=V(i,j)。

6、*D/νIfReJ(i,j)2000ThenΛ(i,j)=64/ReJ(i,j)ElseIfReJ(i,j)3000AndReJ(i,j)2000ThenΛ(i,j)=0.16/ReJ(i,j)^0.25ElseIfReJ(i,j)59.7/εJ^(8/7)AndReJ(i,j)3000ThenΛ(i,j)=0.3164/ReJ(i,j)^0.25ElseΛ(i,j)=(1/(-1.8*Log(6.8/ReJ(i,j)+(εJ/7.4)^1.11)/Log(10)))^2EndIfNextj'计算输油管Q0(i)=2*Q(i,N)V0(i)=4*Q0(i)/(3.14*D0^2)ε=2*0.00015/D0Re(i)=V0(i)*D0/νIfRe(i)2000Thenλ0(i)=64/Re(i)ElseIfRe(i)3000AndRe(i)2000Thenλ0(i)=0.16/Re^0.25ElseIfRe(i)59.7/ε^(8/7)AndRe(i)3000Thenλ0(i)=0.3164/Re(i)^0.25Else-8-λ0(i)=(1/(-1.8*Log(6.8/Re(i)+(。

7、ε/7.4)^1.11)/Log(10)))^2EndIfH0(i)=λ0(i)*C3*V0(i)^2Forj=1ToNH0(i)=H0(i)+(λ(i,j)*C1*vv(i,j)+Λ(i,j)*C2*V(i,j)^2)*2Nextj'显示结果List2.AddItemQ0(i)List1.AddItemH0(i)Nexti'计算输油时间Text11=H0(25)Text12=H0(5)t=2*S*(H0(25)-H0(5))/((Q0(5)+Q0(25))*3600)Text9=tEndSub4程序说明4.1显示封面窗体,启动4.2单击“计算”命令按钮,封面窗体消失,显示基本参数设置界面;此时如果单击“退出”命令,则程序结束.4.3修改各个参数后,单击“计算发油时间”命令按钮,进入计算时间窗体;单击“计算各鹤管流量”命令按钮则进入鹤管界面.单击“返回”按钮,回到启动窗体.4.4进入计算时间窗体后,再输入各个变量.单击确定后显示发油时间.单击“返回”按钮,回到上一级窗体.4.5进入鹤管窗体后,再输入各个变量.单击“输出图形”命令按钮后显示-9-平均流量与液面高度图.单击“返回”按钮,回。

8、到该窗体.4.6在鹤管窗体里,单击“显示各鹤管流量”命令按钮,则显示鹤管数一半的各管流量.单击“返回”按钮,回到该窗体.5运行结果及结果分析。