全国大学生数学建模论文A题79686

1、2010高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置

2、报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 2010 年 9 月 13 日赛区评阅编号（由赛区组委会评阅前进行编号）：2010高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：地下储油罐的变位分析与罐容表标定摘要加油站地下储油罐在使用一段时间后，由于地基变形等原因会发生纵向倾斜及横向偏转，导致与之配套的“油位计量管

3、理系统”受到影响，必须重新标定罐容表。本文即针对储油罐的变位时罐容表标定的问题建立了相应的数学模型。首先从简单的小椭圆型储油罐入手，研究变位对罐容表的影响。在无变位、纵向变位的情况下分别建立空间直角坐标系，在忽略罐壁厚度等细微影响下，运用积分的方法求出储油量和测量油位高度的关系。将计算结果与实际测量数据在同一个坐标系中作图，经计算得误差均保持在3.5%以内。纵向变位中，要分三种情况来进行求解，然后将三段的结果综合在一起与变位前作比较，可以得到变位对罐容表的影响。通过计算，具体列表给出了罐体变位后油位高度间隔为1cm的罐容表标定值。进一步考虑实际储油罐，两端为球冠体顶。把储油罐分成中间的圆柱体和

4、两边的球冠体分别求解。中间的圆柱体求解类似于第一问，要分为三种情况。在计算球冠内储油量时为简化计算，将其内油面看做垂直于圆柱底面。根据几何关系，可以得到如下几个变量之间的关系：测量的油位高度 实际的油位高度 计算体积所需的高度于是得到罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。再利用附表2中的数据列方程组寻找a与b最准确的取值。一、问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐

5、内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。题目给出了一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。并给出了罐体纵向倾斜变位的示意图和罐体横向偏转变位的截面示意图。请用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用给出的小椭圆型储油罐（两端平头的椭圆柱体）示意图，分别对罐体无变位和倾斜角为a=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模

6、型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。（2）对于实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。二、问题分析本题是一个在罐体变位后重新标定罐容表的问题，就是需要得出变位后油位高度与油料体积的关系，然后在油料高度间隔为1cm或10cm的情况下，算

7、出所有高度所对应的体积值，即可得到新的罐容表标定值。第一问中共做了两次实验，分别为罐体无变位与纵向变位。对于无变位的情况，可以选择合适的体积微元，在油位高度方向积分即可算出油体积与油位高度的关系；对于倾斜角为=4.1º的纵向变位，我们采用二重积分的方法，分三种情况进行计算。先在油位高度方向积分得到任意处油截面的面积，再积分得到体积公式。最后利用附件1中的实际数据对公式的准确度进行检验，并对比变位前后储油量与油位高度关系的差别。第二问中，将储油罐分成三部分进行计算：中间的圆柱体和两端的球冠体。对于与的处理问题，对、已经确定的静态储油罐建立空间直角坐标系，根据几何关系得出测得的油位高度与

8、实际油位高度的关系（含有参数），实际油位高度与计算体积所需的高度、的关系（含有参数），并计算得到储油量关于、的表达式，于是便得到了储油量与测量油位高度及变位参数、的关系式，代入若干组附表2中的实际数据，即可确定与，之后用实际检测数据检验所建模型的正确性与方法的可行性。三、模型假设（1）忽略油罐厚度对油罐容积的影响，认为由图中数据得到的容积即为油罐的标准容积；（2）忽略油罐内各种管道如进出油管道，油位探针所占的体积；（3）不计油浮子的厚度、大小等，认为实验中测得的高度即为油罐底部沿探针到油面的距离；（4）假设油浮子到达最高处时便不再加油。四、符号说明：储油罐任一位置平行于罐底方向实际油位高度；：

9、问题一中建立空间直角坐标系后轴方向上油料宽度的一半；：建立空间直角坐标系后轴方向上的油料长度；：建立空间直角坐标系后轴方向上的变量；：问题一纵向变位第种情况下相应某一高度时的油的体积；：问题一中变位后测得的油料高度；：问题一变位时油料平行于罐底方向的最大高度；：问题一变位情况下用任意平行于罐底平面截得的油料面积；：实际储油罐球冠内储油量；：实际储油罐中间圆柱部储油量；：附表2中编号为的流水号所对应的出油量。五、模型的建立与求解5.1小椭圆型储油罐无变位情况首先以一侧罐底中心为原点，建立如图所示的空间直角坐标系，其中下部阴影部分为油料：xYZ图1 无变位情况下建立空间直角坐标系yXh图2 截面椭

10、圆示意图xXZh从侧面观察得到如下示意图：根据题目中的已知数据，得到椭圆截面的方程式为：于是有 取从上到下叠加的矩形薄片为体积微元，得到体积微元公式：体积微元在轴方向进行积分，得到体积公式： 将该结果与实际测量数据在同一以高度为横坐标，体积为纵坐标的坐标系中作图，得到如下曲线：图3 计算曲线与实际数据对比图从图像上可以看出，计算得到的数据与实际测量数据吻合较好，相对误差始终很小，实际数据稍小可能是由于探针，进出油罐管道等占一定体积及罐壁厚度造成的，为简化模型，本文忽略这部分影响。º纵向变位以椭圆罐底中心为原点，轴，轴平行于罐底，轴平行于油罐侧壁方向建立空间直角坐标系：图4 纵向变位情

11、况下建立空间直角坐标系ZXYH油位探针h0由图4可知：接下来分三种情况进行讨论，通过二重积分即可求得油料体积。第一种情况：当（单位：mm ）时，只有一端罐底接触油面，如图5：图5 第一种情况ZXYH油位探针h0先在轴方向上定积分，得到任意位置油料截面面积：再将视为变量，在轴方向上定积分：其中，代入后解得：第二种情况：当 （单位：mm）时，两端罐底都接触油面，如图6：图6 第二种情况ZXYH油位探针h0代入得：第三种情况：当，其中 （单位：mm）时，一端罐底已经完全被油浸没，如图7：图7 第三种情况ZXYH0油位探针h0且，其中，代入上式解得：将上述三种情况得到的方程式分区间画在同一坐标系中，并

12、与实际测量的数据做对比，得到如下关系图（图8）：图8 变位后储油量与油位高度关系图从图8可以看出，计算得到的公式基本符合实际检测数据。通过代入数据，误差保持在3%以内。因此，在标定罐容表时，我们以得到的公式为基础，代入数据计算即得。将变位前后储油量与油位高度关系图画在同一坐标系中，得到图9：变位前变位后图9 变位前后储油量与油位高度关系曲线对比 结合公式以及图9可以看出罐体变位对罐容表产生如下影响：变位后在油位液面到达探针之前，测量高度始终为0，刚好接触油浮子时，将数据代入公式可计算得此时储油量约为1.75L；在变位后的第一阶段内，曲线斜率小于变位前，这个阶段内储油量变化较慢；第二阶段内，曲线

13、增长趋势与变位前基本一致，即上升相同的高度，储油量增加值基本相等，但由于第一阶段储油量较少，这是储油量比变位前小220L左右；第三阶段曲线变化率逐渐降低，当油浮子的高度为1200mm时，油罐还没有装满，此时的储油量比变位前少约100L。根据假设，为使油位高度与储油量是一一对应的关系，此时不再加油，认为该值即为储油最大值。从0到1200mm每间隔10mm取一数值代入公式得到如下罐容表的标定值：表1 纵向变位后的罐容表标定值油位高度（mm）罐容量(升)油位高度（mm）罐容量(升)油位高度（mm）罐容量(升)油位高度（mm）罐容量(升)01.754310628.8476201883.55930318

14、8.83103.6066320664.2626301926.949403227.36206.33542330700.1886401970.359503265.493010.0434340736.6036502013.89603303.214014.8218350773.4856602057.259703340.515020.7535360810.8146702100.719803377.356027.9142370848.5716802144.159903413.737036.3738380886.7356902187.5610003449.628046.1975390925.28970022

15、30.9410103484.999057.4464400964.2157102274.2710203519.831007057202317.5410303554.1111084.44544201043.117302360.7310403587.81120100.3014301083.057402403.8410503620.88130117.7924401123.37502446.8410603653.32140136.9664501163.837602489.7410703685.08150157.274601204.647702532.5110803716.131

16、60179.4944701245.717802575.1410903746.44170203.1564801287.037902617.6211003775.96180228.0064901328.588002659.9411103804.66190253.9375001370.348102702.0811203832.49200280.8675101412.318202744.0311303859.39210308.7325201454.468302785.7811403885.3220337.4755301496.798402827.3111503910.19230367.04654015

17、39.298502868.611603934.01240397.4025501581.938602909.6511703957.74250428.5015601624.718702950.4411803958.07260460.3085701667.628802990.9511903978.33270492.7885801710.638903031.1712003996.88280525.9115901753.759003071.09290559.6466001796.959103110.68300593.9676101840.229203149.935.2实际储油罐变位分析我们将储油罐分成三

18、段来考虑，两端为球缺，中间为圆柱体。中间部分采用类似第一题的积分方法求解。对于两端的球冠体，若直接积分，结果将十分复杂，为方便计算，同时使误差尽量小，本文把球冠内油液面看做与轴平行。对于纵向与横向都已经变化好的静态储油罐来说，我们以中间圆柱体一侧底面圆心为原点，平行于罐体的轴为轴，平行于油面的轴为轴建立空间直角坐标系。 水平线油位探针油油浮子图10 储油罐纵向变位示意图h0H190º90ºZYH2O根据图10可以得到以下关系式：用垂直于轴的平面去截油罐得到图11所示的储油罐的横向变位截面示意图，图中两个油液面是指将横向变位前后的截面图画在一个图中，并使油位探针方向相同，以方

19、便计算，此时前后液面形成夹角：横向变位后油液面横向变位前油液面h0h图11 储油罐横向变位示意图油位探针为测量值，实际油位高度，根据图像可得如下关系式：综合上面几个式子，可得、与的关系式：球冠体内储油量的计算XY1.625m1mxO根据已知数据容易解得球冠所在球的半径为1.625m，球过球心的截面图如下，以圆心为原点，平行于空间坐标系轴的轴为轴，建立新的平面直角坐标系，阴影部分为储油部分：图12 球冠还原为球后截面图该圆的方程为：表示圆上一点到轴距离，所以：以平行于空间坐标系Y轴的平面去截球冠，得到如下所示截面图：x0.625m图13 球冠体截面图可以得知：所以球冠内油料截面面积为：当球冠内油

20、位高度为时，球冠内储油量为：在计算两端球冠内储油量时，分别用、代替即可求出结果。5.2.2中间圆柱体内储油量的计算计算方法与第一问中类似，用垂直于轴的平面去截得到如下截面示意图：XZh图14 圆柱部截面示意图截面圆的方程为：于是得到：又有：即：于是该截面面积：由于有转折点，又要分三种情况讨论，分别求解。当（单位：m）时当（单位：m）时当（单位：m）时其中为圆柱体的总体积用Matlab积分得到的结果过于冗长，不便于写在正文中，具体结果见附录。5.2.3参数，的确定由于第二种情况的可能性最大，数据最多，所以在求解参数与时，利用附表2中显示油高值在中间部分的值进行计算。由于显示的油量容积是利用没有变

21、位情况下的公式计算得到的，不是真实值，故不能加以利用。附表2给出了出油量与显示油高的对应数据，我们用差值计算，即利用累计出油量与油高的变化值的对应关系求解a、b。取流水号分别为323、337、351的三组数据，令：于是得到如下方程组：用Matlab7.0求解该方程组，得到一组解a=1.6º，b=0º于是便得到了变位后储油量与油位高度的关系式，间隔10cm取值代入得到如下罐容表标定值：表2 变位后实际储油罐罐容表标定值油位读数(m)罐容量(L)油位读数(m)罐容量(L)0.0 25.67711.6 33580.30.1 331.11141.7 36401.20.2 1123.

22、381.8 39199.90.3 2352.81.9 41960.70.4 3883.32.0 44667.70.5 5656.842.1 47304.20.6 7634.742.2 49853.20.7 9787.132.3 52296.50.8 12089.12.4 54614.50.9 14518.72.5 56785.71.0 17056.32.6 58785.41.1 19683.52.7 60584.71.2 223832.8 62146.21.3 25138.32.9 63414.91.4 27933.43.0 63856.41.5 30752.5将得到的关系曲线六、模型的评价与推广6.1模型的评价本题主要运用微积分的方法与立体几何的相关知识建立数学模型，进而求出罐内油料体积与测量