《储油罐的变位识别与罐容表标定模型数学建模》由会员分享，可在线阅读，更多相关《储油罐的变位识别与罐容表标定模型数学建模（26页珍藏版）》请在金锄头文库上搜索。
1、2010高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置
2、报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 2010 年 9 月 13 日赛区评阅编号（由赛区组委会评阅前进行编号）：2010高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：储油罐的变位识别与罐容表标定摘 要本文旨在对储油罐的变位情况进行分析，并建立体积积分模型对罐容表的标定值进行求解计算。根据题意要求对发生变位
3、，即纵倾与横滚情况下的储油罐罐容表进行正确的标定定位，需要标定的油罐分为两种，通过分析我们考虑到，对于椭圆柱体油罐和实际油罐其初步的分析情况是一致的，都需求解出两侧的面积，并在此基础上对圆柱的变化情况进行积分，并根据其特有的结构进行一定的近似积分求解。本题所涉及到主要工作有体积积分模型的建立，数据的拟合，误差的分析以及罐容表的标定。对于第一问，要求对椭圆型的储油罐进行罐容表标定，由于椭圆油罐的倾斜角是已知的，故我们可以根据正切关系求解出测量高度下的左侧面的液面高度，为了积分求解的方便，我们建立椭圆曲线方程式，并移动坐标轴至椭圆低端，进一步在转化为极坐标的情况下，对相应液面高度的椭圆面积进行积分
4、，从而得到左侧液面高度的面积值。通过分析我们考虑到，由于液面是按一定的规律从左至右倾斜而下，以及考虑到斜液面的积分分为三段高度，故我们在求解出左侧面积的基础上对斜率变化下的椭圆面积进行分段积分，由此得到椭圆柱体储油罐相应高度的体积情况，通过对原始数据的比较，发现有一定的误差，并按油面的高度成二次曲线的关系，最后我们根据体积积分模型以及结合误差关系曲线得出了120个间隔为1cm高度的椭圆柱油罐体体积。对于第二问，要求对实际的球冠柱体的储油罐进行罐容表的标定，由于题中所给的纵倾角和横滚角是未知的，故我们需要在建立模型的基础上用最小二乘法拟合得到角度的参数。通过分析我们考虑到附件2中的油量累加体积数
5、据可能存有很大的误差，而两个时间点下少量的出油体积数据是准确无误的，故拟合参数时，必须要把前一个高度下计算得出的体积量减去下一个的体积量。在建立模型时，我们使用了坐标变换的方法并结合使用两个切面的二重积分方法得出了中间圆柱体的体积，进一步的我们对两个球冠头依次按照横切和纵切的一重积分法得出其相应体积，由于切割面积是建立在相对于油罐的水平情况下，而液面的高度是水平于地平线的，故还存有一个小切块体积，我们同样使用三角形近似切割积分法对两个球冠头下剩余的切块体积进行求解。综合以上三个积分体积之和以及在考虑到液面的分段情况下，我们使用逐步搜索法对角度参数进行精度搜索并结合最小二乘法得出了的值为。最后我
6、们根据得出的角度值计算出了30个间隔为10cm高度的储油罐体积值。本文最后对模型的优缺点进行了分析，并对在相互变化的情况下进行了灵敏度分析。文章中我们使用了大量的图表使抽象的理论情况更为形象化。关键词：最小二乘拟合 体积积分 纵倾横滚 坐标变换 逐步搜索法1 问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏
7、转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。请利用罐体变位后在进/出
8、油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。2 问题分析本文主要是让我们根据小椭圆型模型求得实际罐容与高度的关系式，进而标定罐容表。对于小椭圆模型，想到应该先求出无变位时的模型，这部分是规则题，应该能求出具体表达式。在得到无变位模型后，我们分析变位和无变位之间的区别仅仅在于变位后，罐内油位高度不是定值，而是沿着一个平面线性变化。这时计算体积，想到只要把每个地方的椭圆面面积都求出来，进而积分就可以得到。而椭圆面的截面积在无变位模型时已求解出，
9、我们只需将始终的高度用现在的高度表达式替换即可。在得到第一问的模型后，分析第二问中实际储油罐与小椭圆储油罐的区别在于中间部分为圆柱体，两边为球罐体。中间部分与第一问中的椭圆柱体本质上是一致的，只是加上横向偏转后，实际各点高度与油位计的高度关系更为复杂。两边的球罐体可以采用分割近似求解来求。3 模型假设1 假设在理论计算体积时，不考虑罐壁厚度；2 假设理论值与实际值的误差主要来自于罐体内部结构； 3 假设给出的实验数据与真实值很接近，可以作为检验数据；4 假设罐体结构不会随温度、储油量或其它环境因素的影响而发生变形。4 符号说明无变位时油罐体罐容； 变位时油罐体罐容；油位计测量罐内油位高度； 截
10、面面积。5 小椭圆储油罐罐容标定模型（问题一）为了掌握罐体变位后对罐容表的影响，本节利用小椭圆型储油罐（两端平头的椭圆柱体）来模拟实际储油罐，分别针对罐体无变位和倾斜角为a=4.100的纵向变位两种情况建立数学模型，与实验数据进行比对，分析模型的正确性，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。本节内容安排如下：l 罐体无变位模型l 罐体变位模型l 结果与误差分析l 罐容表标定5.1 罐体无变位模型问题要求我们建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。考虑到罐体倾斜后，涉及到倾角，在用几分求体积表达式时比较复杂，为此我们先研究简单情
11、况罐体无变位时罐容（体积）与油位计测量罐内油位高度的关系式，即为罐体无变位时的数学模型。下面分析求解过程及计算方法。5.1.1 罐容与油位高度的关系式无变位模型为使问题数学化，我们采用解析几何的方法，首先建立如图所示三维直角坐标系，如图5.1所示。图5.1 三维直角坐标系我们要求罐容，也就是图中虚线以下部分体积，与油浮子高度的关系式，显然需要用到多重积分。已知三重积分的积分函数时，其在空间有界闭区域上得三重积分即为有界闭区域的体积，即 （5-1）根据上述知识，我们只需要确定有界闭区域的空间表达式，即可在上对函数积分得到的体积。（1）积分空间的表达式图5.1给出的是椭圆形储油罐的正面图，我们画出
12、其侧面图（椭圆），将其放在与图5.1一致的直角坐标系内，如图5.2所示。图5.2 储油罐侧面图我们可以清楚地看出侧面是一个椭圆，已知其半长轴和半短轴，可以很快写出其表达式如下： （5-2）在确定侧面椭圆表达式后，我们来确定上面表达式，上面应是平行与平面的平面，其表达式为： （5-3）在确定这两个面以后，还需知道的范围，从图中可以看出 （5-4）由此，我们可以得到空间有界闭区域应为曲面，平面，围成的闭区域，即： （2）闭区域的体积无变位时罐容表达式三重积分求解过程如下：将投影到面上，得投影区域如图5.3所示。图5.3 投影区域则应有 （5-5）只需求解双重积分即可得出体积。显然双重积分应为图中阴影部分的面积，记为。由于满足椭圆关系式，这里采用极坐标计算。先进行坐标变换：双重积分 （5-6）令 ，由积分限 ，得 ，只取到椭圆的右半部分，故 （5-7）最终得到无变位时罐容表达式为： （5-8）5.1.2 分析理论与实验结果我们将附件1中无变位进油和无变位出油时测得的油位高度带入式（5-8），求出对应的罐容理论计算值，如附表1.1所示，MATLAB程序见附录2.1。这里给出部分理论值表，如表5.1所示。表5.1 无变位时部分罐容值油位高度/mm实际罐容/L理论罐容/L误差/L100.00163.594159.02312322.8810.88176.14362374.63
