基于重积分与MATLAB.doc

基于重积分与MATLAB的油罐注油体积计算 摘要 本文运用重积分及立体几何的有关理论，根据油罐的位置变化对实际储油罐的罐容表的标定值进行相应修改。这样，可使加油站有效利用燃油及提高加油站对燃油的有效管理。平头小椭圆型储油罐问题中，油罐的形状为侧面椭圆，两端平头。油罐的倾斜只在纵向上而没有横向倾斜，同时根据已知的倾斜角度要求重新给出罐容表标定值。而圆头圆柱型储油罐问题中，油罐形状变为中间罐身为圆柱体，两端为球灌体。同时油罐既有纵向倾斜和横向倾斜，要求根据所给数据确定油罐位置变化的横向及纵向角度。根据这些特点及要求，我们对平头小椭圆型储油罐问题运用重积分对倾斜后的燃油体积进行积分，同时对不同液面分情况进行讨论；对圆头圆柱型储油罐问题，我们先运用油罐的空间几何特点对油罐进行合适的分割，再运用三重积分的换元的方法分别对燃油体积进行积分，最后相加解决问题。最后利用最小二乘拟合法求出油罐倾斜的横向及纵向角度参数。 在平头小椭圆型储油罐问题中，我们建立适当坐标系，利用油罐倾斜角度给出油面倾斜斜率及几何分析，对油的体积进行积分。在积分过程中再讨论不同液面高度（主要分为液面左右两端分别位于椭圆上或下两部分）下积分后表达式的不同。最后利用MATLAB软件，对积分表达式编写程序，求出具体不同液面高度下实际燃油的体积，从而得到修改后的罐容表标定值。 在圆头圆柱型储油罐问题中，由于油浮子不是在油罐的中心位置，因而需要考虑油罐的向左向右倾斜。而左右倾斜的区别只在于油浮子高度的影响，只需稍加变化即可。我们将油罐依圆柱及两个球灌体分割成三个部分，从而将燃油依此分割分为三个部分。这里用到了“分而治之”的思想。由于油罐的横向倾斜对三个部分中燃油体积没有影响，因此只需将油浮子读数做相应修改得到实际燃油液面高度即可。然后考虑纵向倾斜的情况，此时，对罐身圆柱体中燃油体积的计算依照平头小椭圆型储油罐问题的模型求解；对两个球灌体中燃油体积，同样分别建立适当坐标系，给出体积的三重积分公式，这里运用积分的极坐标变换，以及几何知识，分析得出积分限，从而得到积分结果。最后将得到的三个体积公式相加即为燃油的实际体积公式。最后利用MATLAB软件运用最小二乘拟合法编写程序，拟合出油罐的横向及纵向倾斜角度，得出以下结果： 右倾（右端在下、左端在上）：； 左倾（左端在下、右端在上）：。 关键词：重积分 MATLAB 体积 倾斜 偏转 分而治之 问题重述 对于储存燃油的地下储油罐，运用的是“油位计量管理系统”，该系统采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以此得到罐内油位高度和储油量的变化情况。 在初始时刻，储油罐处于水平状态，而使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。因此需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。 在本文中我们要研究解决储油罐的变位识别与罐容表标定的以下问题。 （1）利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为(=4.10的纵向变位两种情况做了实验，数据见附件1，并建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 （2）对于图1所示的实际储油罐，建立罐内储油量与油位高度及变位参数（纵向倾斜角度(和横向偏转角度( ）之间的一般关系，数据见附件2，并根据所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验模型的正确性与方法的可靠性。 符号说明 符号 说明 备注 倾斜的角度的大小 在问题1中为 偏转的角度的大小 a 椭圆的短轴长的一半 b 椭圆的长轴长的一半（问题1） 图1中油罐圆柱部分的半径（问题2） 在问题2中b为1.5米 L 液面离椭圆中心的距离 h 油浮子的高度 油浮子高度为h时，油的体积 表1：符号说明表 模型建立 在本模型中，我们的中心任务是根据油罐中油浮子的高度来得到油罐中注油的体积。因此下文中会用到一系列积分来得到体积。 图5：情况1侧面透视图 问题1： 情况1 如图5所示，将油罐竖直放置，则AC为油的液面，此时液面的最高处不超过短轴长的一半，最低处高度为0（相对稍有倾斜的横放椭圆底部而言），为的大小，AC离底面椭圆圆心的距离为L（L根据x轴 的方向取正值或负值），ABC面垂直于底面，并设底部椭圆的方程为： 其中a为椭圆短轴长的一半，b为椭圆长轴长的一半。 在被积区域内取被积△A