储油罐的变位识别与罐容表标定.docx

2011高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 孟龙 2. XXX 3. XX 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：储油罐的变位识别与罐容表标定摘要通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”。储油罐在使用一段时间后，由于地基变形等原因，使罐体发生变位，现通过建立模型求解罐体油位高度与罐容表标定值之间的关系。本文针对涉及到的3个问题，进行了合理的解答。我们对问题一第1问使用了高等数学中求定积分的方法求解。对问题一第2问借用第1问的解可类似的使用定积分的方法求解。对第二问可使用分解求积的方法求解。对于问题一第1问我们使用高等数学中的定积分法首先建立了无变位模型。以小椭圆油罐的侧面为横截面，首先借助MATLAB软件寻找横截面积与油位高度之间的函数关系。然后对油罐的长度进行积分从而得到油罐储油量和油位高度之间的函数关系，通过和无变位的实验数据通过百分误差进行判断，可以得出最终模型为对于问题一第2问我们仍然使用定积分的方法建立纵向变位模型。已知小椭圆油罐侧面的横截面积与油位高度之间的关系，现分五种情况对油罐储油量和油位高度之间的关系进行讨论，对模型进行合理的理论证明和推导，可得到大致的关系为再通过MATLAB软件编程求得只有纵向变位的数学模型。并通过和实验数据对比给出罐体变位后油位高度间隔为1cm的罐容表标定值。对于问题二我们将整个罐体分解为左侧球冠体，右侧球冠体以及中间部分计算油罐储油量和油位高度之间的关系。在问题一的第2问建立的纵向变位模型基础上我们再次分五种情况建立双变位模型。现通过枚举法求得最佳的 和 分别为2.14和4.6，并通过求得的 和 与未变位情况进行对比并寻找差值的规律。最后给出预测罐体变位后油位高度间隔为10cm的罐容表标定值。关键词：储油罐 罐体变位 罐容表 定积分 多重积分 MATLAB目 录第一部分 问题重述.(3)第二部分 问题分析.(3)第三部分 模型的假设.(4)第四部分 定义与符号说明.(4)第五部分 模型的建立与求解(5)1问题1的模型(5) 模型I(无变位模型)(5) 模型II(纵向变位模型).(7)2问题2的模型(11)模型I(双变位模型)(11)第六部分 对模型的评价(14)第七部分 参考文献(15)一、问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。现提出下列问题：（1） 为了掌握罐体变位后对罐容表的影响，利用简化的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为a=4.1的纵向变位两种情况做了实验。根据实验结果，建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。（2） 对于实际的储油罐，建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。利用罐体变位后在进/出油过程中的实际检测数据，根据所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用给出的实际检测数据来分析检验模型的正确性与方法的可靠性。二、问题分析（一）问题一的分析（1）第1问属于使用多重积分求解体积函数的问题。解决此类问题需要首先求得横截面积与函数坐标之间的关系，即小椭圆油罐的侧面积与油位高度之间的函数关系，使用MATLAB软件积分可顺利求得。在此基础上对油罐的长度再次进行积分，从而可以得到油罐储油量和油位高度之间的关系。通过计算推导得到的关系和实验数据之间的方差判断模型的可行性。（2）第2问仍然属于使用定积分法求解体积函数的问题。但和第1问不同的是需要考虑油罐的纵向变位角 。利用第1问推导出的弓形面积公式，以油罐的长为z轴建立坐标系，分5种情况讨论油罐储油量和油位高度之间的关系。再次通过计算推导得到的关系和实际数据之间的误差判断纵向变位模型可行性，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。（二）问题二的分析问题二属于使用分解法求积分的问题。解决此类问题可将需要求积的物体分解为若干部分，此问题中我们将油罐分解为三部分，即中间的圆柱部分以及两端的球冠体部分。其中中间的圆柱部分可使用问题一第2问的结果直接得到。两端的球冠体部分使用增补法转化为简化的规则模型，在此基础上将坐标系建立在左侧球冠体与圆柱相接处，通过枚举法求得最佳的 和后再使用计算推导出油罐储油量和油位高度之间的关系。并和实际数据进行百分误差的计算验证模型的可行性。通过求得的模型给出罐体变位后油位高度间隔为10cm的罐容表标定值。利用所给数据做出油罐储油量的双变位理论值和未变位显示值，并以此判断差值的规律。三、模型假设1. 假设题目所给的数据真实可靠；2. 假设小椭圆油罐形状规则；3. 假设小椭圆油罐厚度可忽略不计；4. 假设储油罐在罐空和罐满之间均可随意进出油；5. 假设储油罐在发生变位后其形状不发生改变；6. 假设储油罐发生两种变位的角度不超过10；7. 假设计算出的储油罐容积即为储油罐的储油量；四、定义与符号说明 ：小椭圆油罐侧面长半轴； ：小椭圆油罐侧面短半轴； ：油罐侧面横坐标； ：油罐侧面纵坐标； ：油罐显示油位高度； ：油罐长度； ：小椭圆油罐侧面弓形面积； ：油罐储油量；：油罐的存储量； ：油罐液面高度到油罐顶部距离； ：油罐等效油位高度； ：未横向偏转时油位高度； ：油罐侧面横截面半径；五、模型的建立与求解第一部分：求解问题一（一） 无变位模型1.模型建立根据问题一（1）的分析，首先需要求得小椭圆油罐的侧面积与油位高度h之间的函数关系，建立如图所示坐标系现设椭圆的面积公式为 则弓形面积为 。使用MATLAB求解可以得到弓形面积，现对油罐长度进行积分，即 。将得到的解与实验数据进行作图计算百分误差并对模型进行调整得到最终结果。2.模型求解首先通过MATLAB积分工具箱计算弓形面积，解出得到弓形面积为：再对油罐长度进行积分可以得到油罐储油量和油位高度的函数为：通过MATLAB作图可以得到理论推导出的函数与实验采集数据拟合的较好（程序见附件一）。途中可以发现在增大的过程中存在一定的误差值，因此可判定为油罐体厚度等影响所致。现使用MATLAB编程计算推导数据和实验数据之间的百分误差，误差稳定在3.49%左右，因此在推导出的油罐储油量V和油位高度h的函数中除以1.0349，得到最终的模型为其中a=0.89，b=0.6，l=2.45。以下给出小椭圆油罐罐容表的理论值和实际值理论值：油位高度/mm储油量/L油位高度/mm储油量/L油位高度/mm储油量/L油位高度/mm储油量/L105.12310813.136102027.919103231.392014.43320850.226202070.039203267.253026.45330887.676302112.139303302.674040.62340925.476402154.29403337.635056.62350963.626502196.229503372.116074.243601002.086602238.189603406.087093.313701040.856702280.079703439.5280113.713801079.926802321.889803472.4290135.333901119.266902363.69903504.73100158.084001158.867002405.210003536.45110181.894101198.717102446.6910103567.54120206.694201238.87202488.0510203597.98130232.444301279.17302529.2610303627.72140259.074401319.627402570.3210403656.74150286.534501360.337502611.2110503685.01160314.84601401.227602651.9210603712.47170343.824701442.277702692.4310703739.1180373.564801483.497802732.7410803764.841904044901524.857902772.8310903789.65200435.095001566.348002812.6811003813.46210466.85101607.948102852.2811103836.21220499.125201649.668202891.6211203857.83230532.025301691.468302930.6811303878.23240565.465401733.368402969.4611403897.3250599.435501775.328503007.9211503914.92260633.915601817.348603046.0711603930.92270668.875701859.48703083.8711703945.09280704.295801901.518803121.3211803957.1290740.155901943.638903158.4111903966.42300776.446001985.779003195.112003971.54实际值：油位高度/mm储油量/L油位高度/mm储油量/L油位高度/mm储油量/L油位高度/mm储油量/L159.02312438.121312677.632312892.823168.83176.14362450.41362678.542315.83892.843168.91192.59412462.621412690.532365.83906.533218.91208.5462474.781462690.822367.06920.453268.91223.93512486.891512702.852417.06934.613318.91238.97562498.951562714.912467.06949.053368.91253.66612510.971612727.032517.06963.83418.91268.04662522.951662739.192567.06978.913468.91282.16712534.91712751.422617.06994.433518.91296.03762546.821762763.72666.981010.433568.91309.69812558.721812764.162668.831026.993618.91323.15862570.611862776.532718.831044.253668.91336.44912582.481912788.992768.831062.373718.91349.57962594.351962801.542818.831081.593768.91362.561012606.222012814.192868.831102.333818.91375.421062618.092062826.952918.831125.323868.91388.161112629.962112839.832968.831152.363918.91400.791162641.852162852.843018.831193.493968.91413.321212653.7522128663068.83425.761262665.672262879.323118.83（二） 纵向变位模型1 模型建立已知小椭圆油罐侧面的弓形面积的计算方法，现在不同情况下通过弓形面积的计算公式找到储油量和油位高度之间的函数关系。（1） 情况1：图一如图一所示，当纵向变位后，储油量很小时，油位高度不可测量，则此时的储油量无法测量，其中在0到一定的范围之间。则有其中（2） 情况2：图二当 时，储油量的表达式如下所示，示意图见图二。（3） 情况3：当 时，储油量的表达式如下所示，示意图见图三。图三（4） 情况4：当 时，储油量为整个小椭圆油罐体积减去上面空缺部分的体积，即 图四 其中a=0.89，b=0.6，l=2.45。由可推出因此等效的油位高度为图五（5）情况5：当时，因油面高度不可测量，因此储油量在一定值到 之间。2 模型求解在已建立的模型基础上，我们使用MATLAB软件编程求解所列出的公式（程序见附录二）。从而得到如下解。（1） 情况1经过计算可以得到，当油位高度较小不可测量时，储油量的范围是：（2） 情况2由题意可知，此时取值为，将其带入公式，得到油位高度的范围是现通过程序求得储油量的范围是：（3） 情况3当时，即时，求得储油量的范围是：（4） 情况4当时，即时，求得储油量的范围是：（5） 情况5当时，可以求得储油量的范围是：将模型求解所得到的值与实验测得数据使用（1）中的百分误差法可以求得模型求解所得到的数据和实验测得数据误差不超过5%，因此可认为模型有效。通过MATLAB编程将纵向变位前后储油量和油位高度的关系曲线表现在图像中，并求得两者差值最大为195.5293L。由于变位后储油量具有一定值时油位高度为0，因此变位后储油量变小。综上所述，现给出小椭圆油罐纵向变位后的罐容表的理论值理论值：油位高度/mm储油量/L油位高度/mm储油量/L油位高度/mm储油量/L油位高度/mm储油量/L001.67 310629.94 6201884.44 9303188.74 103.53 320665.36 6301927.80 9403227.22 206.26 330701.29 6401971.20 9503265.30 309.97 340737.71 6502014.62 9603302.97 4014.75 350774.60 6602058.05 9703340.22 5020.69 360811.93 6702101.48 9803377.01 6027.85 370849.69 6802144.90 9903413.34 7036.31 380887.85 6902188.29 10003449.17 8046.13 390926.41 7002231.65 10103484.49 9057.38 400965.33 7102274.95 10203519.26 10070.11 4101004.61 7202318.19 10303553.48 11084.37 4201044.23 7302361.36 10403587.11 120100.23 4301084.16 7402404.44 10503620.12 130117.71 4401124.40 7502447.41 10603652.48 140136.88 4501164.93 7602490.28 10703684.17 150157.78 4601205.74 7702533.02 10803715.14 160180.21 4701246.80 7802575.62 10903745.36 170203.94 4801288.11 7902618.07 11003774.80 180228.84 4901329.64 8002660.36 11103803.40 190254.81 5001371.40 8102702.46 11203831.11 200281.77 5101413.35 8202744.38 11303857.88 210309.66 5201455.50 8302786.09 11403883.65 220338.43 5301497.82 8402827.58 11503908.32 230368.03 5401540.30 8502868.84 11603931.80 240398.40 5501582.93 8602909.85 11703953.94 250429.52 5601625.70 8702950.60 11803976.60 260461.34 5701668.59 8802991.08 11903995.50 270493.84 5801711.59 8903031.26 12004012.704110.10 280526.97 5901754.68 9003071.13 12004012.704110.10 290560.72 6001797.87 9103110.68 300595.05 6101841.12 9203149.89 第二部分：求解问题二1. 模型建立正视图首先考虑油罐未发生横向偏转的情况（即），我们把油罐车分为主体的圆柱体和两端的球冠体三个部分，其中主体的圆柱体的计算方法可以利用上题中计算椭圆柱的方法。现在不同情况下通过计算弓形面积的计算公式找到储油量和油位高度之间的函数关系（程序见附录三）。所建立的坐标系如图所示（图中仅作出球冠体）。（1） 情况1：如图一所示，当纵向变位后，储油量很小时，油位高度不可测量，则此时的储油量无法测量，其中在0到一定的范围之间。则有（2） 情况2：当时，储油量的表达式如下所示，示意图见图二。（3） 情况3：当时，储油量的表达式如下所示，示意图见图三。（4） 情况4：左视图当时，储油量为整个小椭圆油罐体积减去上面空缺部分的体积，即 其中R=1.5，l=8。由可推出因此等效的油位高度为俯视图（5） 情况5：当达到3时，因油面高度已经达到油位的最大高度，但是此时油罐并未满，因此储油量在（即情况4中取2R时，的值）到之间。l 下面开始考虑有横向偏转的情况：（1） 当时，（2） 当时，（3） 当时2. 模型求解通过分别对和赋值对比可知，只有当，超过10时才有可能出现其他情况，由假设可排除这种可能，现通过110以步长为0.1进行枚举，发现当分别为2.1和4.6。通过已经求得的的值，我们用已经得到的模型进行求解。立体示意图当时，所能到得最大值为；此时的满足上述的情况2；故，此时的油量为55.77L.当时，此时的满足上述情况2；故，此时的油量为343.61L.当时，此时的满足上述情况2；故，此时的油量为947.22L.经过计算，当时，通过计算得到的满足上述情况3；当时，此时的满足上述情况4；故，此时的油量为63785.06L.通过模型所得到的数据与实际进出油量的数据进行对比，使用问题一第1问中的百分误差法可以求得模型求解所得到的数据和实验测得数据误差不超过5%，因此可认为模型有效。相应的油位高度间隔为10cm的油量如下：油位高度/mm储油量/L油位高度/mm储油量/L油位高度/mm储油量/L油位高度/mm储油量/L055.77 80011387.76 160032571.86 240053660.75 100343.61 90013759.66 170035378.87 250055871.85 200947.22 100016245.24 180038166.79 260057921.34 3002075.99 110018826.69 190040924.94 270059780.07 4003492.44 120021485.87 200043635.20 280061416.53 5005172.68 130024203.76 210046281.40 290062777.25 6007069.67 140026968.94 220048844.55 300063785.06 7009150.30 150029763.45 230051311.40 六、模型的评价与推广1.优点：问题一第1问中引入了误差测定的方法，并通过此方法将求解得到的模型进行了优化，使模型更加切合实际测量值。问题一第2问中将油量的可能值详细的分为5种情况，使模型所求解更加精确可靠。问题二中引入了数学中的分解求积法和增补法简化了模型，将求解过程变得方便但不失准确性。本文在计算求解时多次使用MATLAB软件，提高了计算效率。2.缺点：题目假设偏转角度不大于10，对模型的求解有了一定的限制。在考虑储油量时将储油罐容积作为储油量进行计算求解，增加了一定的误差。问题一第2问求得的模型较为复杂，因此并没有使用第1问使用的方法对所求模型进行优化，因此存在一定误差。问题二使用增补法虽然简化了模型，但在一定程度上也增加了误差。3.模型推广此模型不仅仅可以用作地下储油罐的容积测定，还可以推广到其他液体储存器械未变位及变位的储存量进行测定。同时可将实际测量出的偏转角度带入模型求出储存量与液面高度的关系，增加储存器械显示储存量的正确性。七、参考文献1 姜启源编，数学模型（第二版），北京：高等教育出版社，1993。2 卓金武，MATLAB在数学建模中的应用，北京：北京航空航天大学出版社，2011。3 韩中庚编，数学建模方法及其应用（第二版），北京：高等教育出版社，2009。4 Mark M. Meerschaert著 刘来福等译，数学建模方法与分析（第三版），北京：机械工业出版社,2009。5 许丽佳, 穆炯主编，MATLAB程序设计及应用，北京：清华大学出版社,2011。15