储油罐的变位识别与罐容表标定数模竞赛.doc
《储油罐的变位识别与罐容表标定数模竞赛.doc》由会员分享,可在线阅读,更多相关《储油罐的变位识别与罐容表标定数模竞赛.doc(48页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、【精品文档】如有侵权,请联系网站删除,仅供学习与交流储油罐的变位识别与罐容表标定数模竞赛.精品文档.2010高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃
2、处理。我们参赛选择的题号是(从A/B/C/D中选择一项填写): 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名): 参赛队员 (打印并签名) :1 2. 3. 指导教师或指导教师组负责人 (打印并签名): 日期: 2010 年 9 月 13 日赛区评阅编号(由赛区组委会评阅前进行编号):2008高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号储油罐的变位识别与罐容表标定摘要本文是与储油罐的变
3、位识别与罐容表标定相关的数学问题。问题一中,我们运用高等数学中的微元的思想,分别建立了储油罐在变前位、变位后的数学模型,并对其求解,得出变位前与变位后,在相同油高下,间隔为1cm的罐容表标定值(详见39页)。油体高度变位后的体积变位前的体积121.6016.1718222.2024.46211093658.103786.66201103687.703809.9516通过对变位前、变位后标定值的对比,在时得出倾斜对罐容表的影响为:倾斜之后的油罐,油体体积的标定值在任意高度下都要比倾斜之前标定体积要小。当储油罐内装油过多或过少时,倾斜对油体体积的标定值影响较小;当储油罐内油面深度接近罐体中线时影响
4、最大。问题二中,储油罐发生纵向变位的同时还发生横向变位。通过几何分析找出倾角、偏角与油体高度间的函数关系。按照问题一中的分析方式,利用微元法得到体积与深度、纵向倾角和横向偏角的关系式,但方程比较复杂,不能够得到解析解。因此我们采用近似算法,找出相邻油层间的体积,利用最优化理论得到、的最优解,再将代入上述关系式中求解,即得到罐体变位后油位高度间隔为的罐容表标定值。(详见41页)高度标定值高度标定值403.68183626058.4239505.41510627060.2533607.35832928061.851我们用和高度算出体积,计算出理论出油体积,与实际出油体积作对比,算出相对误差均值在0
5、.1%左右,证明我们模型的准确性较高。最后,我们用蒙特卡洛法对模型中的高度进行扰动时,发现对高度添加一个服从正态分布随机扰动值时产生的相对误差与的值产生的相对误差在一个数量级。从而证明了我们所建立的模型具有较好稳定性。一、问题重述(略)二、问题分析加油站的地下储油罐通过预先标定的罐容表(即罐内油位高度与储油量的对应关系)进行实时计算,以得到罐内油位高和储油量的变化情况。罐体发生变位后会发生纵向倾角和横向偏角,影响油浮子的位置,致使读出的油位高度与实际高度不符,从而导致罐容表发生改变。通过对附表1数据的分析,可以看出进油曲线和出油曲线完全重合,即没有漏油存在,且流量计读数完全正确。我们不考虑温度
6、、压力及不计浮盘上下浮动对罐内油量体积的影响,但要考虑到油罐内有注油管道和出油管道(未达油罐底部)以及油位探测管道(焊接至油罐底部),会占一部分体积。考虑到在实际情况中,压力及温度对油罐变形的影响,为确保安全储油罐不可能完全装满,也不可能在油量很少的情况下供油,且储油罐的变位不可能很大,纵向倾角和横向偏角都在小范围内。简化的储油罐如下图,我们只考虑在油量在阴影部位的情况。 简化的储油罐2.1对问题一的分析:我们首先找出油罐倾斜前、倾斜后体积与倾角,高度间的函数关系,然后通过互相对比得出倾角对罐容表的影响程度。未倾斜的油罐:如图:对于横截面为椭圆形的储油罐,封闭的罐体内油面为水平面,在统一的坐标
7、系下,油罐内处处油体深度相等,且处处横截面面积相等。只需计算出任一截面面积,并对其长度方向积分即可求出体积。通过得到的体积与高度间的关系式,即可得到任意高度时储油罐中的理论储油量。但实际中储油罐内含有输油、进油管道等其他影响油体体积的因素。其标定值可能不等于理论值。我们通过修正得到它们之间的差值,从而找到正确的标定值。倾斜的油罐:罐体纵向倾斜时,油深不再是处处相等,各截面面积也此时倾角、油体高度对体积都有影响。我们打算通过微元法建立他们的函数关系,对其进行积分求解。同样找出修正值,从而得出油罐倾斜时的正确标定值。2.2对问题二的分析:附件2中给出的数据是在倾角和偏角同时影响下的实际检测数据。同
8、时考虑倾角和偏角对罐容表标定的影响,为了简化数学建模的求解,我们分为两步求解,为此我们首先考虑横向偏转的情况,找出偏转之后油体的实际高度与读数高度、偏转角之间的函数关系。将读数高度用偏转角和油体实际高度的函数表达,即建立:这样既可将偏转后实际油体高度转为读数高度,再考虑纵向变位时,就可将问题转化为问题一的情况。从而得到的关系式。但问题二中的储油罐的主体是圆柱体,两端为球冠体。这就使得建立积分方程比较复杂,无法得出解析解。我们采用两种方法来解决这个问题。1. 对积分方程进行数值解;2. 采用近似法求得两油层间体积的近似解,步骤如下:1)首先找出液面所在的截面面积与高度、变位角度之间的函数关系;2
9、)取上下两层距离为的两油层,相邻高度的体积可以用一个近似柱体来代替,其计算公式为:表示上层油面面积,表示下层油面面积,表示两油层间的间距 3)由题中所给数据,得到相邻高度间的采集出油量,用最小二乘法思想建立优化模型,可得到的估计值。把高度h、代入积分函数,便可以算出任意高度体积。从而给出油位高度间隔为的罐容表标定值。相比之下,第二种方法比第一种方法更易编程实现,因此我们采用第二种方式进行建模。三、模型假设1、流量计所测数据准确;2、储油罐为刚体,即油罐不会变形3、忽略因温度变化而对油体容积的影响;4、不计油品的静压力作用而引起得油体的膨胀和收缩的容积;5、油浮子与标尺为点接触,不会因为油浮子影
10、响油体体积;6、罐体发生变位时纵向倾角的取值范围为,横向偏角的取值范围为四、符号说明-出油罐内油体的读数高度-油体的实际高度-椭圆的长半轴-椭圆的短半轴-储油罐的长度-按照读数高度得出的油罐内油体体积-纵向倾角-横向偏角-椭圆横截面面积-两相邻为的油层间的体积-上层油面面积-下层油面面积-两油层间的间距-真实体积部分变量在使用时说明 五、模型的建立与求解5.1问题(一)微分方程模型5.1.1油罐未变位时的模型根据对前面对问题一的分析,首先建立了未变位时,储油罐内油量与油位高度的数学程模型1。建立如图5.1的坐标系,设油罐内液面高度为,为恒定值;椭圆的长半轴为,短半轴为,xyh简化小椭圆储油罐A
11、AAA图5.1椭圆轨迹方程为: (5.1)对半弓横截面(阴影部分)进行积分得出面积: (5.2)由(1)式得出为: (5.3)将(3)式代入(2)式得: (5.4)积分后可得: (5.5)面积再乘以长度即为体积: (5.6)其中,半弓椭圆横截面面积,表示油罐长度,表示油罐内油体体积。代入已知值以及变值,用求解得出结果部分见表1,详见附录1。累加进油量/L实际累加体积/L油位高度/mm油高/m计算体积/m3实际体积/m3差值50312159.020.159020.322882590.3120.01088259100362176.140.176140.3746330340.3620.0126330
12、34365739191152.361.152364.0556137683.918910.136703768370739691193.491.193494.1073621283.968910.138452128累加出油量/L实际累加体积/L油位高度/mm油高/m计算体积/ m3实际体积/ m3差值52.7239161150.721.150724.0527979523.916190.136607952102.738661123.991.123994.0010492213.866190.1348592213653316.2160.480.160480.327208490.316190.0110184
13、93703266.2142.620.142620.2754733340.266190.009283334表1注:实际累加体积是在罐内油量初值的基础上累次加上每次进油量得出,计算体积是按照上述函数式(5)代值求解的。不考虑其他任何影响因素的情况下,油深应该与高度的关系为(6)式所得积分值。则计算体积应该与实际体积相等,没有差值产生。但由结果可见计算体积与实际体积之间有差值。对题中所给附件1中数据分析可知,储油罐在13:20以前一直不间断进油,储存油量到达最高值,之后一直出油,油量减到最低值。差值按照先增后减的趋势排列。我们考虑差值与油深间存在着对应关系,利用拟合函数找出它们间的关系式: (5.7
14、)并求出相关系数,相关系数越接近1,说明线性关系很强。将计算值减去差值得出正确值,该值与真实值十分贴近。这也验证了我们前面的分析,即不可忽视油罐内的注油管道和出油管道(未达油罐底部)以及油位探测管道(焊接至油罐底部)占去的体积。图5.2 出油量和进油量误差排序后的拟合曲线算出的油体体积与真实值不符,它们之间有一个差值即真实体积=计算体积-差值,联系(6)(7)式:根据此式可计算出未变位时,油位高度间隔为的罐容表标定值(见表2详见附录4):h(m)V/m3v(L升)0.010.01617216.171750.020.02446224.462110.030.03554535.544591.073.
15、7370553737.0551.083.7623433762.3431.093.7866623786.6621.13.8099523809.952表2 未倾斜时罐容表标定值5.1.2 油罐发生纵向变位时的模型: 图5.3建立如图5.3的三维坐标。设通过油位计读出的油深为;某一微小截面处油深为,值不恒定,随着油深改变;已知椭圆的长半轴为;短半轴为;油罐逆时针变位纵向倾角。如图5.3在区间内任意去某一处微小段,同时用与底面平行的平面与液体面相截,得到图中阴影部分微元面积元的体积: (5.8)由上述(5)中推出的弓形面积共识式可知: (5.9),之间的几何关系见下图, 为已知数据:fdZhh (5.
16、10)联立(7)(8)(9)式得: (5.11)对上式进行积分: (5.12)再代入各项已知数据,利用MATLAB求出结果如下表。由于篇幅有限,在此仅列出部分结果(见表3),详见附录5。累加进油量(升)油位高度实际体积计算体积实际与计算体积的差值747.86411.29962.861.01000.0472797.86423.451012.861.05830.04553297.731034.243512.733.5694574790.05673299.741035.363514.743.5732283220.0585累加出油量(升)油位高度实际体积计算体积实际与计算体积的差值501020.653
17、464.743.5231136560.05837371001007.733414.743.4780872070.06334722500425.831014.741.0678373580.05309742550411.73964.741.0117861410.0470461表3由(11)式就可得出倾角为时,小椭圆储油罐在不同高度下的计算体积。同样,油罐倾斜时的计算体积与实际体积间仍有差值,差值=计算体积-实际体积。实际体积是在罐内油量初值的基础上累次加上每次进油量得出,计算体积将值代入函数式(11)中求出的。同样按照5.1.1的方法将差值按照先增后减的趋势排列(见附录6),油罐倾斜时的差值与读数
18、高度同样存在对应关系。利用MATLAB拟合工具箱CFTOOL对差值和高度进行分析(程序代码见附录7),得出差值与油面深度呈非线性关系为(见图5.4): (5.13) ,拟合度好。图5.4利用函数式(11)计算出的油体体积与真实值不符,它们之间有一个差值即:真实体积=计算体积-差值 (5.14)按照题目要求即可得到油位高度间隔1cm的罐容表标定值(见表3,详见附录8):h(cm)h(m)v(m3)v(L升)180.180.2136213.6190.190.2378237.8200.200.2630263.01071.073.59653596.51081.083.62773627.71091.09
19、3.65813658.11101.103.68773687.7表3 倾斜之后的罐容表标定值说明:标定高度范围在之间。原因:由前面的假设,纵向倾角和横向偏角都,讨论极限情况,即油体液面刚好位于线处,假设这时的倾角为最大值。几何关系见图5.5f图5.5为已知数据,则则油面深度的标定区间为。5.1.3倾斜前后罐容表标定值比较(详细见附录8):油体高度未倾斜时罐容表标定值倾斜之后的罐容表标定值同等高度下倾斜与未倾斜油量差值)0.0116.1717521.6-5.40.0224.4621122.22.30.0335.5445923.811.71.083762.3433627.7134.61.093786
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 储油罐 变位 识别 罐容表 标定 数模 竞赛
限制150内