自流系统装车工艺计算

油库设计与管理题目自流系统装车工艺计算专业年级指导教师姓名学号姓名学号姓名学号姓名学号姓名学号姓名学号第1页共14页目录1、题目.......................................................................22、原理分析...................................................................22.1自流装车系统的数学模型分析.............................................22.2自流装车时间的计算.....................................................33、程序框图...................................................................54、程序变量声明...............................................................64.1鹤管变量...............................................................64.2集油管变量.............................................................64.3输油管变量............................................................64.4储油罐变量............................................................64.5其他变量..............................................................65、程序代码...................................................................65.1窗体设置程序..........................................................65.2第一根鹤管流速v1(1)和λ1(1)的计算....................................75.3第一根集油管流速v2(1)和λ2(1)的计算..................................85.4第i根鹤管流速v1(i)和λ1(i)的计算....................................85.5第i根集油管流速v2(i)、λ2(i)和每根鹤管流量qq(j,i)的计算............95.6输油管流量Q0(j)、流速v0和λ0的计算.................................105.7所需油罐液位高度z（j）的计算.........................................105.8插值法找到初始液位高度z1与发油后的液位高度z2下所对应的第一根鹤管的流量.........................................................................115.9平均流量与装车时间的计算.............................................115.10程序的输出..........................................................116、运行界面..................................................................126.1鹤管数20、液位高度为3m时所对应的各鹤管流量..........................126.2鹤管数30、液位高度10m时所对应的各鹤管流量...........................136.3鹤管数为30，发油量为1500m3，开始发油时液位高度为12m所对应的发油时间.137、上机实验总结..............................................................14第2页共14页1、题目装卸系统，如图2所示，装油管设置在集油管中部，且两侧鹤管同时装油，两侧鹤管对称布置。油库油料为车用汽油，操作条件下车用汽油的粘度为0.6×10-6m2/s，鹤管采用Φ108×4的钢管，鹤管计算长度为25m，鹤管间距为12.5m，集油管为Φ219×6，输油管为Φ159×5，输油管计算长度为1500m，油库库容为30000m3，单个罐的容量为3000m3，油罐内径为16m，罐壁高度为15.85，罐出油口与装油鹤管出口之间的高差为30m。（1）计算鹤分别管数为20和30，液位高度为3m和10m时所对应的各鹤管流量？（2）操作鹤管数为30，发油量为1500m3，开始发油时液位高度为12m，计算发油时间？2、原理分析2.1自流装车系统的数学模型分析为求得在罐内液位高度一定时各鹤管流量，需建立自流装车系统水力计算的数学模型。如图1所示的自流装车系统，设装车系统的位能为Z+H0，Z为罐内液位高度，H0为罐出口与装油鹤管出口之间的高差。qi、Qi；di、Di；li、Li；λi、Λi、和ui、Ui分别为第i根鹤管和i段集油管的流量、管径、长度、摩阻系数和流速。Q0、D0、L0、λ0和U0分别为输油管的流量、管径、管长、摩阻系数和流速。图2自流装车系统第3页共14页如图1所示，从流体力学的基本方法可以得出，鹤管2和集油管的交汇点2的能量，必然是供给流经鹤管2和鹤管1及集油管1的流量q2和q1的能量，在此点上两个系统的摩阻和动能之和必须相等，即gudlgugUDLgudlgu222222222222211112111121（2）以此类推gudlgugUDLgudlguiiiiiiiiiiiiii2222222211112111121（3）因各鹤管流速分别为u1，…，ui，…，un，其流量分别为q1，…，qi，…，qn，则可求得各段集油管的流速为21114DqU，22212)(4DqqU，221)(4iiiDqqqU输油管的流速为20210)(8DqqqqUni设计自流装车系统时，一般可先假设鹤管和集油管的直径和长度。这样在方程（2）中，若知道u1就能求出u2，由方程（3）类推下去，便可求出ui直到un，同时也求出了U1，…，Ui，…，Un和U0。因此，鹤管、集油管和输油管的直径和长度已确定后，对应油罐内的液位高度Z，方程式（1）成为关于u1的一元非线性方程，可运用弦截法求解方程（1）可得出某时刻一定油罐液位高度Z对应的第一根鹤管流速u1，再由方程（3）求得其余鹤管的流速，从而计算出各段集油管及输油管的流速和流量。2.2自流装车时间的计算在设计自流装车系统时，必须满足装车时间的要求。即在一定的位差（H0+Z）下，所选装油系统的管径必须满足预定的流量。由于油罐的液位高度的变化，油品装车时间是油品液位高度的函数。设油罐的横截面积为F，输油管流量为Q0，在dτ时间内，流向油罐车的油品为Q0dτ，油罐液位下降为dZ，流出油品为FdZ，它们必然相等，即dZd0FQ（4）在自流装车系统的各管段直径和长度、鹤管数确定后，输油管流量Q0与油罐内的液位高度Z的函数关系可用下式表达。)(0ZfQ（5）将式（5）代入式（4）积分得第4页共14页ZZfFZQFZZd)(d2121ZZ0（6）式中τ—自流装车时间，s；Z1—装车前罐内的液位高度，m；Z2—装车后罐内的液位高度，m。如果需要装车的油品体积为Vr，Z2可按下式计算：FVZZr12（7）由于式（5）无法用解析式表达，式（6）可采用数组插植计算自流装车时间τ。第5页共14页3、程序框图YesYesNo程序启动计算v1(1),q1(1),λ1(1),v2(1),q2(1),λ2(1)i=2计算v1(i),q1(i),λ1(i),v2(i),q2(i),λ2(i)i=i+1j=n计算q0,v0,λ0,q=q+Δqq=mNo计算装车后液位Z2=Z1-V/A,用插值法求出Z1,Z2对应的末端鹤管的流量q1,q2利用无穷分割累加求其平均流量Qe，进而计算出装车时间T程序结束输入基本数据l1,d1,l2,d2,l0,d0….输出结果第6页共14页4、程序变量声明4.1鹤管变量鹤管内径d1，长度l1,总根数n，流速v1()、v1(n),假设的摩阻系数λ1()、λ1(n),雷渃数Re1()、Re1(n),相对粗糙度ξ1()、ξ1(n),计算的摩阻系数λ10()、λ10(n),划界雷渃数Re11()、Re11(n),划界雷渃数Re12()、Re12(n)，第一根鹤管的假设流量q，每根鹤管流量qq（j，i）4.2集油管变量集油管内径d2，长度l2,流速v2()、v2(n),摩阻系数λ2()、λ2(n),雷渃数Re2()、Re2(n),相对粗糙度ξ2()、ξ2(n),划界雷渃数Re21()、Re21(n),划界雷渃数Re22()、Re22(n)，总的摩阻损失hl4.3输油管变量输油管内径d0，长度l0,流速v0，假设的摩阻系数λ0，雷渃数Re0，相对粗糙度ξ0，划界雷渃数Re01,划界雷渃数Re02，平均流量Qe4.4储油罐变量储油罐直径D，液位初始高度z1,发油后高度z2，第一根鹤管流量下所对应的液位高度z（j）4.5其他变量发油量V，油品黏度ν，管道绝对粗糙度Δ，油罐出口与鹤管出口间的高差H05、程序代码5.1窗体设置程序PrivateSubForm_Load()Label14.Visible=FalseText13.Visible=FalseText15.Visible=FalseIfText12.Text=12Then'注Text12.Text的输入需在属性窗口输入Label14.Visible=TrueText13.Visible=TrueText15.Visible=TrueText14.Visible=FalseEndIf第7页共14页EndSub5.2第一根鹤管流速v1(1)和λ1(1)的计算PrivateSubCommand1_Click()Dimv1(),λ1(),Re1(),ξ1(),λ10(),Re11(),Re12(),v2(),λ2(),Re2(),ξ2(),Re21(),Re22()DimQ0(20000),z(20000),qq(20000,30)AsSingleν=Val(Text1.Text):ν=ν*10^-6'对输入数据赋值，并化为国际单位Δ=Val(Text2.Text):Δ=Δ*10^-3n=Val(Text3.Text)d1=Val(Text4.Text):d1=d1*10^-3l1=Val(Text5.Text)l2=Val(Text6.Text)d2=Val(Text7.Text):d2=d2*10^-3d0=Val(Text8.Text):d0=d0*10^-3l0=Val(Text9.Text)D=Val(Text10.Text)H0=Val(Text11.Text)z1=Val(Text12.Text)V=Val(Text13.Text)ReDimv1(n),λ1(n),Re1(n),ξ1(n),λ10(n),Re11(n),Re12(n),v2(n),λ2(n),Re2(n),ξ2(n),Re21(n),Re22(n)π=3.14Forq=0.000001To0.02Step0.000001'假设第一根鹤管流量为0.000001m3/s,该鹤管每次增量也为0.