油库上机编程实验报告(共9页).doc
精选优质文档-倾情为你奉上 油库设计与管理 实验报告小组成员李伟王文康陈万贵陈实曹增辉赖家鑫 学号 实验名称油库设计与管理指导教师宋晓琴单位西南石油大学1题目来源某油库的自流装卸系统，如图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程序框图：开 机输入参数 n 计算常数 =+V(1) (1) hf1(1) u(1) hf2(1)i=2 to n/2(i)= (i-1)+ hf1(i) (i)+<+ YESu(i-1) (i) hf2(i) i<n/2=2(+)； = +hf1(1)+hf2(1)+hf2(n/2)+ <+Z YES输出Q1(1)Q1（n/2）； 输入发油量V0 算Z1=Z0-V0/F；调用过程再算Q1T=2V0/（Q1+Q0）结 束3程序代码3.1定义数据和数据输入3.2系统常数3.3定义数组3.4计算过程3.5计算第一根鹤管Option Base 1Private Sub Command1_Click()'定义数据和数据输入Dim dd, D, ll, L, D0, , H, J, aa, A, C1, C2, C3, , , t, S, E, Fdd = Text5D = Text4D0 = Text3N = Text1 = Text2ll = Text8L = Text7L0 = Text6S = Text10'系统常数H = 2 * 0.00015 / ddJ = 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)'定义数组Dim Q0(50)Dim V0(50)Dim 0(50)Dim Re(50)Dim H0(50) As String * 5Dim vv(50, 50)Dim V(50, 50)Dim (50, 50)Dim (50, 50)Dim qq(50, 50)Dim Q(50, 50)Dim ReH(50, 50)Dim ReJ(50, 50)Dim B(50, 50)'计算过程qq(1, 1) = 0.01 = 0.For i = 1 To 30'计算第一根鹤管 If i <> 1 Then qq(i, 1) = qq(i - 1, 1) + (i - 1) * End If Q(i, 1) = qq(i, 1) vv(i, 1) = qq(i, 1) / aa ReH(i, 1) = vv(i, 1) * dd / If ReH(i, 1) < 2000 Then (i, 1) = 64 / ReH(i, 1) ElseIf ReH(i, 1) < 3000 And ReH(i, 1) > 2000 Then (i, 1) = 0.16 / ReH(i, 1) 0.25 ElseIf ReH(i, 1) < 59.7 / H (8 / 7) And ReH(i, 1) > 3000 Then (i, 1) = 0.3164 / ReH(i, 1) 0.25 Else (i, 1) = (1 / (-1.8 * Log(6.8 / ReH(i, 1) + (H / 7.4) 1.11) / Log(10) 2 End If'计算第一段集油管 V(i, 1) = qq(i, 1) / A ReJ(i, 1) = V(i, 1) * D / If ReJ(i, 1) < 2000 Then (i, 1) = 64 / ReJ(i, 1) ElseIf ReJ(i, 1) < 3000 And ReJ(i, 1) > 2000 Then (i, 1) = 0.16 / ReJ(i, 1) 0.25 ElseIf ReJ(i, 1) < 59.7 / J (8 / 7) And ReJ(i, 1) > 3000 Then (i, 1) = 0.3164 / ReJ(i, 1) 0.25 Else (i, 1) = (1 / (-1.8 * Log(6.8 / ReJ(i, 1) + (J / 7.4) 1.11) / Log(10) 2 End If '计算后面鹤管 For j = 2 To N qq(i, j) = qq(i, j - 1) Do qq(i, j) = qq(i, j) + 0.00001 vv(i, j) = qq(i, j) / aa ReH(i, j) = vv(i, j) * dd / If ReH(i, j) < 2000 Then (i, j) = 64 / ReH(i, j) ElseIf ReH(i, j) < 3000 And ReH(i, j) > 2000 Then (i, j) = 0.16 / ReH(i, j) 0.25 ElseIf ReH(i, j) < 59.7 / H (8 / 7) And ReH(i, j) > 3000 Then (i, j) = 0.3164 / ReH(i, j) 0.25 Else (i, j) = (1 / (-1.8 * Log(6.8 / ReH(i, j) + (H / 7.4) 1.11) / Log(10) 2 End If E = vv(i, j - 1) 2 / 19.6 + (i, j - 1) * C1 * vv(i, j - 1) 2 + (i, j - 1) * C2 * V(i, j - 1) 2 F = vv(i, j) 2 / 19.6 + vv(i, j) 2 * (i, j) * C1 If Abs(E - F) < 0.05 Then Exit Do Loop'计算后面集油管 Q(i, j) = Q(i, j - 1) + qq(i, j) V(i, j) = Q(i, j) / A ReJ(i, j) = V(i, j) * D / If ReJ(i, j) < 2000 Then (i, j) = 64 / ReJ(i, j) ElseIf ReJ(i, j) < 3000 And ReJ(i, j) > 2000 Then (i, j) = 0.16 / ReJ(i, j) 0.25 ElseIf ReJ(i, j) < 59.7 / J (8 / 7) And ReJ(i, j) > 3000 Then (i, j) = 0.3164 / ReJ(i, j) 0.25 Else (i, j) = (1 / (-1.8 * Log(6.8 / ReJ(i, j) + (J / 7.4) 1.11) / Log(10) 2 End If Next j'计算输油管 Q0(i) = 2 * Q(i, N) V0(i) = 4 * Q0(i) / (3.14 * D0 2) = 2 * 0.00015 / D0 Re(i) = V0(i) * D0 / If Re(i) < 2000 Then 0(i) = 64 / Re(i) ElseIf Re(i) < 3000 And Re(i) > 2000 Then 0(i) = 0.16 / Re 0.25 ElseIf Re(i) < 59.7 / (8 / 7) And Re(i) > 3000 Then 0(i) = 0.3164 / Re(i) 0.25 Else 0(i) = (1 / (-1.8 * Log(6.8 / Re(i) + ( / 7.4) 1.11) / Log(10) 2 End If H0(i) = 0(i) * C3 * V0(i) 2 For j = 1 To N H0(i) = H0(i) + (i, j) * C1 * vv(i, j) + (i, j) * C2 * V(i, j) 2) * 2 Next j'显示结果 List2.AddItem Q0(i) List1.AddItem H0(i)Next i'计算输油时间 Text11 = H0(25) Text12 = H0(5) t = 2 * S * (H0(25) - H0(5) / (Q0(5) + Q0(25) * 3600) Text9 = tEnd Sub4程序说明4.1显示封面窗体,启动4.2单击“计算”命令按钮,封面窗体消失,显示基本参数设置界面;此时如果单击“退出”命令,则程序结束.4.3修改各个参数后,单击“计算发油时间”命令按钮,进入计算时间窗体;单击“计算各鹤管流量”命令按钮则进入鹤管界面.单击“返回”按钮,回到启动窗体.4.4进入计算时间窗体后,再输入各个变量.单击确定后显示发油时间.单击“返回”按钮,回到上一级窗体.4.5进入鹤管窗体后, 再输入各个变量.单击“输出图形”命令按钮后显示平均流量与液面高度图.单击“返回”按钮,回到该窗体.4.6在鹤管窗体里,单击“显示各鹤管流量”命令按钮,则显示鹤管数一半的各管流量.单击“返回”按钮,回到该窗体.5运行结果及结果分析专心-专注-专业
