椭圆形封头卧式容器不同液面高度的容积计算

PAGEPAGE1椭圆形封头卧式容器不同液面高度的容积计算新疆工学院孟永彪在设计卧式容器时，常常要计算不同液面高度所对应的容积，有时还需列出容积—液位高度对照表或图。例如，在盛装有毒有害介质的卧式储罐设计中，要根据体积充装系数确定最高液面高度并加以标识。在一般资料中仅能查到容器的全容积计算公式，而要计算不同液面高度下的容积则需设计者自行推导公式计算。本文以标准椭圆形封头卧式容器为例介绍不同液面高度下的容积计算方法，并以液化石油气储罐为例编制了QUICKBASIC程序，此法仅供大家参考。卧式容器的组成卧式容器是由筒体和两封头组焊而成（如图1），常用的封头为标准椭圆封头。卧式容器计算简图及说明计算简图如图2。L———筒体长度（两封头切线间的距离，含直边段长度）Di———封头及筒体内直径hi———封头曲面深度不同液面高度下封头的容积计算如图2，可假想将卧式容器两端的曲面部分合并，则形成一个完整的椭球面。其中，a=b=Ri因此，椭球面的方程为：推导出：当容器内的液面高度为h时（如图3所示）。封头的容积公式推导：对其积分得从上式可看出，h变化，V1也随之变化。不同液面高度筒体的容积计算在计算筒体的容积时，忽略尺寸公差及制造误差等因素，可将其断面方程为x2+y2=a2的一圆柱体进行计算，那么如图3所示液面高度的筒体容积为：令：y=acosθdy=-asinθdθ当y=-a时，θ=π；当y=h时，代入公式积分得：卧式容器在不同液面高度下的容积通过以上V1,V2的计算公式，可计算出卧式容器在不同液面高度下的容积之和V：3利用QUICKBASIC语言进行卧式容器的容积计算要计算不同液位高度下的容积以表格、曲线的形式列出是很麻烦的，因此本文利用简便易行的QB编制程序，当然编程语言可以有多种，本文愿起到抛砖引玉的作用。例如：一台内径与长度任意的液化石油气（以丙烷为主）储罐，再不同特征温度下，根据体积充装系数，得出的最大液位高度；程序中从h=0.1m开始，每隔1mm计算一次容积，循环比较，得出最佳值。5DIMp(7),H(7),T(7)8INPUT"请输入储罐公称容积(立方米):";GC10INPUT"请输入罐体内半径(m):";R20INPUT"请输入罐体(两封头切线间)长度(m):";L30INPUT"请输入最高工作温度下允许充装系数:";c40INPUT"请输入罐体实际全容积(立方米):";v60INPUT"请输入-20℃时丙烷饱和液体的密度(千克/立方米):";p(1)70INPUT"请输入-10℃时丙烷饱和液体的密度(千克/立方米):";p(2)80INPUT"请输入0℃时丙烷饱和液体的密度(千克/立方米):";p(3)90INPUT"请输入10℃时丙烷饱和液体的密度(千克/立方米):";p(4)100INPUT"请输入20℃时丙烷饱和液体的密度(千克/立方米):";p(5)110INPUT"请输入30℃时丙烷饱和液体的密度(千克/立方米):";p(6)120INPUT"请输入45℃时丙烷饱和液体的密度(千克/立方米):";p(7)130FORI=1TO7140k=c*p(7)/p(I)145T(I)=k150VK=v*k160FORH=.1TORSTEP.001170v1=3.1416*(R^2*H-H^3/3+2*R^3/3)/2180VA=(R^2-H^2)^.5182VB=VA/H183VC=ATN(VB)184VD=H*SQR(1-H^2/R^2)/R185v2=L*R^2*(VD-VC+3.1416)210VH=v1+v2220DH=ABS(VK-VH)230IFDH<=.1GOTO300240NEXTH300H(I)=H+R310NEXTI312LPRINTGC;"立方米液化石油气储罐各工况条件下最高充装液位高度和充装系数"313LPRINT314LPRINT""315LPRINT316LPRINT"计算条件"317LPRINT328LPRINT"储罐公称容积(立方米):";GC330LPRINT"罐体内半径(m):";R340LPRINT"罐体(两封头切线间)长度(m):";L350LPRINT"最高工作温度下允许充装系数:";c360LPRINT"罐体实际全容积(立方米):";v370LPRINT"-20℃时丙烷饱和液体的密度(千克/立方米):";p(1)380LPRINT"-10℃时丙烷饱和液体的密度(千克/立方米):";p(2)390LPRINT"0℃时丙烷饱和液体的密度(千克/立方米):";p(3)400LPRINT"10℃时丙烷饱和液体的密度(千克/立方米):";p(4)410LPRINT"20℃时丙烷饱和液体的密度(千克/立方米):";p(5)420LPRINT"30℃时丙烷饱和液体的密度(千克/立方米):";p(6)430LPRINT"45℃时丙烷饱和液体的密度(千克/立方米):";p(7)440LPRINT450LPRINT""460LPRINT"计算结果"470LPRINT520LPRINT"-20℃时最高允许充装液位高度(m)";H(1);"最高允许充装系数:";T(1)530LPRINT"-10℃时最高允许充装液位高度(m)";H(2);"最高允许充装系数:";T(2)540LPRINT"0℃时最高允许充装液位高度(m)";H(3);"最高允许充装系数:";T(3)550LPRINT"10℃时最高允许充装液位高度(m)";H(4);"最高允许充装系数:";T(4)560LPRINT"20℃时最高允许充装液位高度(m)";H(5);"