卧式椭圆形封头贮罐容积的计算.docx

卧式椭圆形封头贮罐容积的计算

摘要：本文采用几何推导的方式，分析了卧式椭圆形封头的贮罐不同液位高度下的容积的计算方法，并采用matlab编程和多项式拟合，推导出其相应公式。
关键词：卧式椭圆形封头matlab容积

CalculationofLiquidLevelforHorizontalTankwithEllipticalHead

Abstract:Geometricaldeductionisusedtocalculatethevolumeofthehorizontaltankwithellipticalheadfordifferentlevel,andMATLABandfittingofapolynomialmethodareintroduced,theformulaisexpressed.
Keywords:HorizontalTank,EllipticalHead,MATLAB,Volume

改良西门子法制备多晶硅工艺过程中[1]，需要大量的容器贮存三氯氢硅及其它液体副产品。容器的数量较多，体积较大，一般采用卧式椭圆形封头的贮罐。及时和准确地掌握贮罐的存量及其变化，对于多晶硅生产过程的物料平衡和动态分析具有重要的意义。
本文通过对卧式椭圆形封头贮罐几何结构的分析，采用分别计算正圆筒部分和椭圆形封头部分任意液位高度下的容积的方法，推导出任意液位高度的容积计算公式，为多晶硅及其他化工企业中所采用椭圆形封头贮罐的容积计量提供准确而快捷的方法。
1正圆筒部分容积
设椭圆形封头贮罐正圆筒部分的内半径QN长为R，母线UN长为L（L包括两边封头直边长度），装入液体高度NL高为H，如图1。

图1正圆筒部分结构示意图
正圆筒端头截面为半径为R的圆，沿垂直方向y位置取dy高度的面微元，可知此微元的面积为ds，

则相对于此面微元的液体体积微元为dV，

在区间进行定积分，得到高度为H时正圆筒部分液体的体积为V1，

2椭圆形封头部分容积
设封头长半轴为a，短半轴为b。由于两端封头是旋转体并且结构形式完全相同，呈中心对称，故可以将两端封头合并成一个椭球体有助分析，椭球体结构如图2所示。

图2椭球体结构示意图
垂直方向坐标为y时水平面与椭球面相切得到FGCEF面，取面FGCEF的包络线上任意点D，令线段OA长度为m，线段AB长度为z，线段BD长度为x，则OB长度为。由于D点同时在长半轴为a，短半轴为b的椭圆上，故有

从上式可见点D在面FGCEF的包络线上的坐标满足椭圆方程，故封头部分液体液面构成的包络线FGCEF为椭圆。AE长度为e，AC长度为c，易知，

其面积为

在y位置取体积微元，即图2中阴影部分所示，其体积微元之体积dV，

由于是正圆形椭形封头，故a=R。
在区间进行定积分，得到高度为H时正圆筒部分液体的体积为，

3matlab编程计算贮罐容积
经过上面的分析，卧式椭圆形封头贮罐的体积。以下以圆筒半径R＝1.500m，直筒长L=7.080m，封头矢高b＝0.764m的卧式椭圆形封头贮罐为例，采用matlab软件进行计算不同液位高度时贮罐内液体体积，level2volume.m为本计算的m文件，内容如下：
level2volume.m
clear;
clc;
R=1.5;
L=7.08000;
b=0.76400;
H=0:.05:2*R;
x=H-R;
y=-R;
V1=2*L*(x/2.*sqrt(R^2-x.^2)+R^2/2*asin(x./R)-(y/2.*sqrt(R^2-y.^2)+R^2/2*asin(y./R)));
V2=pi*b*H.^2/R.*(R-H/3);
V=V1+V2;
pp=polyfit(H,V,5);
hi=[0,R,2*R]
Vi=polyval(pp,hi)
ph=poly2str(pp,'h')
plot(H,V)
通过对有限数据的多项式拟合，得到卧式椭圆形封头贮罐的体积V与液位高H关系式如下：

当R=3时，即贮罐最大容积为57.3m3左右。
应当说明的，当在实际应用的时候，由于液位计的零位安装位置比罐底部要高出一定数值，故采用液位计计数计算液体体积时，要加上零位以下的液体体积。
图3为不同液位时，实例采用的贮罐所对应液体体积曲线。

图3体积VS液位
综上所述，卧式椭圆形封头贮罐在任意液位下的容积是可以通过相应液位计算出来的。
参考文献：
1、梁骏吾，电子级多晶硅的生产工艺，中国工程科学.2000,2(12).-34-39