1、1储油罐的变位识别与罐容表标定的数学模型摘要本文解决了储油罐罐容表变位后标定的问题。通过把实际的储油罐抽象成直角坐标系下的几何柱体，然后从区分不同的油面高度入手建立了几何柱体体积的积分模型。再通过合理运用所给数据进行数据拟合，得出了油量体积与油面高度之间的函数关系，进而进行理论与实际体积之间的误差分析和模型可行性分析。针对问题一，首先对于无变位的小椭圆柱体建立了直角坐标系下的容积积分模型（见第 4 页） 。通过 Minitab15 软件对实验数据进行曲线拟合，得出一个油量作为高度的函数关系。利用这个函数关系计算出相应罐容表高度的实际油量容积，对比理论积分模型的容积值，计算出误差值（见表 3 和
2、表 5） 。观察知误差属于正常范围内，则得出通过理论模型来标定的标准罐容表（见第 7 页表 6） 。然后当只有纵向倾斜的变位时，根据柱体内的倾斜油面将柱体容积分为三个部分，分段计算出相对应部分中的容积积分，建立了变位后的分段容积积分模型，通过 Matlab7.0 编程得出容积积分函数（见第 9 页） 。而这个模型是与纵向倾斜角度和油高两个因素有关的。当倾斜角一定时，代入条件数据进行拟合对比，得出模型是合理有效的，从而得出变位后的罐容表（见第12 页表 7） 。最后将每变化 0.01m 的油量变化量与标准罐容表作比，得出比例系数。针对问题二，将储油罐分割成两个球冠和一个圆柱三部分，并将其截面放入
3、平面直角坐标系下建立容积积分模型，分别求出各个部分的油量容积，再相加求总容积（见第 15 页） 。而当纵向倾斜和横向偏转都存在时，考虑将空间直角坐标系作一个相应变换，即把轴乘以相应的三角函数得到新的坐标系，此时积分模型得出的是关于两个倾斜角度和高度的函数。然后根据所给数据作拟合计算出实际油量，且分别选取两个倾斜角度的合理范围，固定高度后代入容积积分函数，将得到的油量与拟合出的实际油量作比较，利用最小二乘的方法从两边逐步逼近，最终得出最优的倾斜角度（见第 17 页）和倾斜后的罐容表（见第 17 页表 8） 。最后，利用附件 2 中的实际检测数据对所建模型的正确性与所运用方法的可靠性进行了分析。关
4、键词 椭圆柱；储油罐；积分模型2一、问题的重述对于加油站的地下储油罐都有与之配套的“油位计量管理系统” ，工作人员采用流量计和油位计来测量进/出油量与罐内油位高度等数据。再通过预先标定的罐容表进行实时计算，以得到罐内油位高度和储油量的变化情况。但由于地形变化而导致储油罐的位置会发生纵向倾斜和横向偏转等变位，使得罐容表的标度不再准确。故要研究变位与罐容表之间的对应关系。其中的条件包括：条件一：实际圆柱型（两端为球冠体）储油罐的尺寸及变位后的各项数据见附录一图 1 至图 3，图 4 是小椭圆型储油罐模型（两端平头的椭圆柱体） 。条件二：利用小椭圆型储油罐模型，分别对罐体无变位和倾斜角为 =4.10
5、的纵向变位两种情况做了实验，实验数据见附件 1 所示。条件三：问题只考虑纵向倾斜和横向偏转的变位，且纵向倾斜角度为 和横向偏转角度为 。要讨论的问题有：问题一：针对小椭圆型储油罐，建立数学模型来研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为 1cm 的罐容表标定值。问题二：针对实际储油罐，建立罐体变位后标定罐容表的数学模型，并用附录二的数据确定 和 的值，给出罐体变位后油位高度间隔为 10cm 的罐容表标定值。进一步利用附件 2 中的实际检测数据来分析检验模型的正确性与方法的可靠性。二、问题的分析对于问题一，首先计算无变位时的罐容表。已知罐容表就是罐内油位高度与储油量的对应关系，
6、故要根据已给小椭圆型的数据，积分得出其体积与高度之间的关系，即建立积分模型。具体做法是取从左往右的左视图，将所截半椭圆面函数作为被积函数。然后乘以罐长即得体积与高度的理论函数。再利用无变位时进油出油的实验数据，拟合出油量作为高度的函数关系。然后令高度为标准罐容表的刻度，通过编程计算出高度间隔为 1cm 的油量体积数，即为标准罐容表。而用拟合出的关系函数，则可以计算出相应高度的计算油量。将理论油量与计算油量做比较，可得出罐容表计量的误差。再考虑变位后的罐容表。由于小椭圆型油罐的两端是椭圆形，故在研究积分时先分段，以到一半长轴为界限，分别计算半椭圆的面积，从而得出被积函数的表达式。然后乘以模型的长
7、，就可以表达出变位时的体积公式。用编程积分得体积与高度的函数，再与拟合的函数比较，得出变位后的罐容表。将标准罐容表与变位后的罐容表数值关于纵向倾斜角度 作比值，就可以得出变位后对刻度的影响。对于问题二，类似于对问题一的分析，首先计算实际标准储油罐的体积与高度的函数关系，同样要建立积分模型。由于标准储油罐的两端是球冠体，因此可以将储油罐分成两部分分别积分，即分为圆柱体和两个球冠体。为了方便积分，可以取从上往下的视角，即俯视图，来对所截平面积分。故应设油高为积分变量，用高度表示不同油量时球冠弧对应的半径，从而计算出球冠处的面积，最后编程积分得理论函数。当有变位时，加入两个角度值为自变量，用实验数据
8、二来代入计算，得出角度与罐容表的函数。3三、模型的假设（1）不考虑储油罐内外温度的变化；（2）油位探针的压力传感器在最下端；（3）出油管等管的体积忽略不计；（4）储油罐在进油后没有其他外在的抽油，即进油后的总量与刚开始出油时等量；四、符号说明误差： ； 高端液高为 时的体积： （单位：m 3） ； HHV椭圆长半轴： ； 椭圆短半轴： ；a b油高： （单位：m） ； 油高坐标值： ；Hh储油罐长度： （单位：m） ； 理论油量： （单位：L） ；L 1V实际进油油量： （单位：L） ； 实际出油油量： （单位：L） ；2V 3储油罐的一半高度： （单位：m） ； 储油罐球冠半径： （单位：m
9、） ；Rr储油罐球冠高： （单位：m） ； s五、对问题一求解的小椭圆柱体容积积分模型5.1 无变位积分模型5.1.1 模型的建立对于小椭圆型罐模型，把它转化成坐标系中的椭圆柱，如图 1；当向里面进油时，油面是水平的，即从左往右看的左视图如图 2 y图 1 椭圆柱 图 2 左视图根据数学分析中的定积分应用知识 1，将储油罐储油量计算问题转化成求椭圆柱容积计算问题来处理。当油面水平时，罐长是一定的，且左面与右面的半椭圆面大小相同。因此只需要计算侧平面图形椭圆的面积。故在图 2 所建的直角坐标系中，椭圆的方程为hx4，221ybxa为了将面积表示成高度的函数，其中横坐标可表示为。 (1)222ya
10、xby取一个高度为 的积分元，即用 作为积分变量。值得注意的是，上式不是标准hd的椭圆方程，其目的是使得油高与纵坐标相一致。由定积分知识得，以椭圆方程为被积函数，分别以油高坐标值和短半轴为上、下限，积分得图 2 中油所覆盖半椭圆面积为S, (2)220Habyd而要求的理论柱体体积为, (3)2210VL通过换元计算积分式(3)得理论体积为2210Habyd。 (4)22 21arcsinHbLb 上式就是表示油高与体积之间的理论函数。根据这个理论函数，令 且按0,1.2H1cm 的增量递增，代入式(4)就可以计算出相应高度的油量体积，即刻画出标准罐容表。5.1.2 模型的求解对于小椭圆型的储
11、油罐模型，已有数据为， ， ，1.780.92a1.206b2.45L将其代入(4)式，得理论容积作为油高的函数关系为 2 21. 1450.6arcsin0.66 .HVH 。 (5)23.0.10.3arcsin1.5.为了先标刻出标准罐容表的刻度，利用问题的约束条件，令 且按 1cm 的0,1.2H增量递增，代入式(6)计算得出理论油量体积，即罐容表。部分数据见如下的表 1：表 1 理论罐容表高度 （m）h理论油量（L）高度 （m）h理论油量（L）高度 （m）h理论油量（L）50.31 839.5066865 0.36 1034.824992 0.41 1238.0842150.32 8
12、77.8381279 0.37 1074.904214 0.42 1279.5213150.33 916.5526829 0.38 1115.284782 0.43 1321.1867510.34 955.6331116 0.39 1155.951504 0.44 1363.0666950.35 995.0626542 0.40 1196.889508 0.45 1405.1475285.1.3 模型的检验为了检验理论油量函数的可实用性，需要将理论值与试验值做比较。有了这个目的，就可以利用对问题所给附件 1 的实验数据中无变位进油的数值，通过 Minitab 软件统计中的数据拟合进行函数拟合，
13、得出拟合油量体积作为油位高度的函数和图 3 如下：。 (6)23249.67043607Vhh图 3 进油时油量体积与油高的函数图令其中的油位高度为标度的 0.01cm，代入式(6)得出对应的实际油量，部分数据如表 2表 2 由拟合函数得的进油时油量高度 （m）h实际油量（L）高度 （m）h实际油量（L）高度 （m）h实际油量（L）0.31 801.5761212 0.36 1034.824992 0.41 1238.0842150.32 839.5066865 0.37 1074.904214 0.42 1279.5213150.33 877.8381279 0.38 1115.284782
14、 0.43 1321.1867510.34 916.5526829 0.39 1155.951504 0.44 1363.0666950.35 955.6331116 0.40 1196.889508 0.45 1405.147528观察表 1 和表 2 的对应高度的油量容积，由此来计算相应于表 1、表 2 中理论进油6量与实际进油量之间的误差。计算误差公式为， (7)210%V对于上述部分数据，在 Excel 中计算得误差值，见表 3表 3 理论油量和实际进油油量的误差高度 （m）h误差（%） 高度 （m）h误差（%） 高度 （m）h误差（%）0.31 3.336009583 0.36 3.
15、462651596 0.41 3.4745182240.32 3.373321679 0.37 3.472456814 0.42 3.4676039140.33 3.404129678 0.38 3.478150707 0.43 3.4583831760.34 3.428956075 0.39 3.480140554 0.44 3.4471521250.35 3.448303383 0.40 3.478809692 0.45 3.434187644以上讨论的是进油时的理论与实际油量，还要研究出油时的误差，从而确定误差的来源与理论的可靠性。由假设（4）可知，进油后的一个小时中并没有其它的出油，即
16、出油初始时的燃油量就是进油后的总量。故首先要将附件一中无变位出油中的累计出油量转化为出油后剩余量，即，30V其中 是已给的累加出油量， 是进油最后的总量，本问题中 （L） 。类似0V03968.1V于表 1 对出油后剩余量和油高通过 Minitab 曲线拟合，得到出油时的油量与油高的函数关系和图 4。 (8)2336.2175436Vhh7图 4 出油时油量体积与油高的函数图将罐容表的刻度值为 0.01m 代入式(8)，对应上面的部分数据的实际出油量见表 4表 4 由拟合函数得的出油时油量高度 （m）h实际油量（L）高度 （m）h实际油量（L）高度 （m）h实际油量（L）0.31 3156.0
17、78446 0.36 2968.163136 0.41 2772.0202260.32 3119.236608 0.37 2929.537418 0.42 2731.9565280.33 3082.010322 0.38 2890.596432 0.43 2691.6467420.34 3044.413424 0.39 2851.354014 0.44 2651.1047040.35 3006.45975 0.40 2811.824 0.45 2610.34425由此来计算相应于表 3、表 4 中理论出油量与实际出油量之间的误差，其表达式为，310%V代入数据计算得表 5表 5 理论油量和实际
18、出油油量的误差高度 （m）h误差（%） 高度 （m）h误差（%） 高度 （m）h误差（%）0.31 4.412699048 0.36 4.508145615 0.41 4.6256322310.32 4.430623746 0.37 4.529521478 0.42 4.6526631130.33 4.449013578 0.38 4.551881496 0.43 4.6810151720.34 4.467990574 0.39 4.575305523 0.44 4.7107390450.35 4.487666891 0.40 4.599866768 0.45 4.7418805348比较表
19、3 和表 5 在相同高度处的误差，发现出油时的理论计算油量与实际油量的差距更大。这种现象可以认为是由于出油时油是通过出油管抽出的，且随着燃油的减少，储油罐中的油面高度也在下降。这样就会有部分燃油因油的粘稠性而附着在储油罐内壁，且出油越多，油面高度下降越大，附着的油量就越多，从而产生的误差就越大。据此分析依据这个角度看这个理论油量的积分模型所得出的油量与油高的函数关系是可行的，故可以使用积分来标定标准罐容表，具体数据见表 6。表 6 标准罐容表高度（m） 油量（L）高度（m） 油量（L）高度（m） 油量（L）高度（m） 油量（L）0.01 4.24803983 0.31 839.506687 0
20、.61 2095.23398 0.91 3339.287240.02 13.8803461 0.32 877.838128 0.62 2138.77986 0.92 3376.360430.03 26.3015684 0.33 916.552683 0.63 2182.30153 0.93 3412.975830.04 40.9482806 0.34 955.633112 0.64 2225.78685 0.94 3449.113050.05 57.4904368 0.35 995.062654 0.65 2269.22365 0.95 3484.750890.06 75.7023752 0.
21、36 1034.82499 0.66 2312.59974 0.96 3519.867250.07 95.4162230 0.37 1074.90421 0.67 2355.90286 0.97 3554.439010.08 116.500171 0.38 1115.28478 0.68 2399.12070 0.98 3588.441970.09 138.846744 0.39 1155.95150 0.69 2442.24085 0.99 3621.850660.1 162.365829 0.40 1196.88951 0.7 2485.25085 1 3654.638240.11 186
22、.980216 0.41 1238.08422 0.71 2528.13812 1.01 3686.776330.12 212.622609 0.42 1279.52132 0.72 2570.88999 1.02 3718.234760.13 239.233522 0.43 1321.18675 0.73 2613.49365 1.03 3748.981370.14 266.759755 0.44 1363.06669 0.74 2655.93617 1.04 3778.981730.15 295.153260 0.45 1405.14753 0.75 2698.20447 1.05 380
23、8.198740.16 324.370271 0.46 1447.41583 0.76 2740.28531 1.06 3836.592250.17 354.370628 0.47 1489.85835 0.77 2782.16525 1.07 3864.118480.18 385.117245 0.48 1532.46201 0.78 2823.83069 1.08 3890.729390.19 416.575673 0.49 1575.21388 0.79 2865.26779 1.09 3916.371780.2 448.713758 0.5 1618.10115 0.8 2906.46
24、249 1.1 3940.986170.21 481.501343 0.51 1661.11115 0.81 2947.40050 1.11 3964.505260.22 514.910033 0.52 1704.23130 0.82 2988.06722 1.12 3986.851830.23 548.912987 0.53 1747.44914 0.83 3028.44779 1.13 4007.935780.24 583.484751 0.54 1790.75226 0.84 3068.52701 1.14 4027.649630.25 618.601107 0.55 1834.1283
25、5 0.85 3108.28935 1.15 4045.861560.26 654.238949 0.56 1877.56515 0.86 3147.71889 1.16 4062.403720.27 690.376170 0.57 1921.05047 0.87 3186.79932 1.17 4077.050430.28 726.991569 0.58 1964.57214 0.88 3225.513872 1.18 4089.4716540.29 764.064764 0.59 2008.11802 0.89 3263.845313 1.19 4099.103960.30 801.576
26、121 0.6 2051.676 0.9 3301.775879 1.2 4104.392645.2 变位后的容积积分模型5.2.1 模型的建立9当小椭圆柱有纵向变位 时，可将其分成、三个部分 2，其中、两部分为楔形体，部分是水平液面为梯形的中部区，具体见图 5。hxLdlH2ahA面 A面图 5 下部楔形体的部分容积计算 图 6 左视图（1）部分下部楔形体积分模型为计算此部分的容积，就是要先计算图 6 的左视图中半椭圆的面积，再对垂直于罐底的各个面积分，即得容积。故在图 5 中，当 ，即 时，为下部dL0tanHL楔形地区。取在距离液面高端（液面高端就是指从左边观察油面的高度） 处，用一垂l
27、直于椭圆柱体的截面，截得 A 面。该区内通过此面的任一高度的油，在 A 面上为一弓形面积 ，其值为Sl2 22 20 0.6.6h haabydbyd220rcsinh 22 .60.60.6ari2ahbhb。22 .rcsin由上式可以看出，这个截面面积跟油高 有关。而要求容积的积分式中应该是罐H长 为积分变量，即高端液高为 时的部分容积为l， (9)0dHVSl观察图 5，寻求 与 之间的关系，有h，得 ，cotlcotldh所以积分变量可换为，tdl当 时， ；当 时， 。0lanhHl0将换好的积分变量代入(9)式，即得出部分的容积积分公式0 0cottHHVSldhSldhy102
28、20 0.6cot.60.6arcsin2Hahhbhb 。26.arcsinbd利用 Matlab7.0 软件编程（见附录二）就可以将上式的积分式解出，即322 20.13cot590459HaVHbHbb2 2375rsin0. arcsin.2 33 2 25900.675ri59bbbHb。 (10)2230.65arcsin7（2）部分中部区积分模型如图 7 所示，当 ，即 时，为中部区，该区高端起始高 可dLta2Hb0H按 计算。0tanHLH2aMhHb1b2图 7 中部区的部分容积计算 图 8 高端 图 9 低端当液面处于此区域时，其液面的形状为两端分别为弦长 b1和 b2截
29、去部分的椭圆面。令此时高端液高为 ，低端液高为 。在下楔形积分模型中， 。在此M0tanHL模型中通过部分总截面是平行四边形这一条件，可以将高端液面处的容积计算转化成部分的最上限，从而由式(10)得液面处于中部区的部分容积计算公式为。 (11)HMV（3）部分上部楔形体积分模型最后一部分是上部楔形体，如图 10 所示，当高端液高 ，即高端液面进入第2Hb部分时，其部分容积按下式计算， (12)2HbKV11也就是总容积与类似部分容积之差，其中 和 都用(10)式表示。2bVHKLGHH00HkM2R图 10 上部楔形体的部分容积计算 图 11 油浮子处高与高端液高关系（4）油浮子标度与高端液面
30、的关系上面分情况讨论的三种情况都是以高端液面 为自变量的积分模型，而问题要求的是油浮子处的刻度高度 与标准罐容表之间的关系。如图 11，利用三角形中的三角H函数关系可得，tanGL其中 是油浮子到左端点的固定距离。这样把 用 表示，再代入上面的式(10)等H三个分段函数，就可以得出油量体积作为油浮子处高度的函数表达式。5.2.2 模型的求解根据式(10)的公式，代入已知数据，即 ， ，得部分下1.780.92a1.206b部楔形体的部分容积函数为32cot0.396.105HVH63263.24.78.10H。550.3arcsin10.4arcsin.80.53对于本题中附件 1 中的变位实
31、验数据知 ，同无变位时的自变量 的取值相.1同，令 且按 1cm 的增量递增，代入上式用 Excel 计算得出变位后相应刻度处0,.2H的理论油量。同理将附件 1 中变位后的进油数据，通过 Minitab 软件进行数据拟合，得出变位后实际的油浮子量出的油面高度，见图 12，其油量与油高的函数关系为，2324.9.53497Vhh12再令自变量 ，得出变位后罐容表中的油量数。0,1.2H图 12 变位后油量与油高的拟合图但由于本问题中在进油前有一部分初始油，则要先算出此部分初始油量的高端液面。即直接令理论油量积分式中的容积 ，计算得出对应的高端高度125V（m） 。而在分段函数中部分的上限值为
32、，当 时，上限值就是0.24oH tanL4.1（m） ，tan0.425tan4.08L比较初始油量高端液面与这个上限值，即 ，这说明由于变位的倾斜角.240.18H度过小，导致初始油量已经处于部分。所以在标油浮子理论刻度时，要根据上限值来分段求解。利用 Excel 计算得出变位后理论与实际油量之间的误差，其值大约都为0.34%左右。这说明得出的理论函数较好，由此求变位后的罐容表如下表 7表 7 变位后罐容表高度（m） 油量（L）高度（m） 油量（L）高度（m） 油量（L）高度（m） 油量（L）0.01 -0.5353702 0.31 535.858544 0.61 1749.58882 0
33、.91 3060.997410.02 -0.1660991 0.32 570.400631 0.62 1793.76760 0.92 3102.143940.03 0.86877809 0.33 605.544492 0.63 1838.03415 0.93 3142.968680.04 2.6928779 0.34 641.263373 0.64 1882.37562 0.94 3183.454450.05 5.4164598 0.35 677.531919 0.65 1926.77920 0.95 3223.583530.06 9.13760072 0.36 714.326017 0.66
34、 1971.23217 0.96 3263.337590.07 13.9446918 0.37 751.622661 0.67 2015.72182 0.97 3302.697640.08 19.9183457 0.38 789.399833 0.68 2060.23549 0.98 3341.64392130.09 27.1327811 0.39 827.636410 0.69 2104.76052 0.99 3380.155810.1 35.6568416 0.40 866.312072 0.7 2149.28427 1 3418.211720.11 45.5547663 0.41 905
35、.407229 0.71 2193.79409 1.01 3455.788950.12 56.8867848 0.42 944.902950 0.72 2238.27730 1.02 3492.863490.13 69.7095893 0.43 984.780911 0.73 2282.72121 1.03 3529.409840.14 84.0767139 0.44 1025.02334 0.74 2327.11308 1.04 3565.400720.15 100.038847 0.45 1065.61295 0.75 2371.44012 1.05 3600.806700.16 117.
36、644087 0.46 1106.53295 0.76 2415.68946 1.06 3635.595770.17 136.938164 0.47 1147.76694 0.77 2459.84818 1.07 3669.732680.18 157.964622 0.48 1189.29891 0.78 2503.90325 1.08 3703.178050.19 181.300351 0.49 1231.11324 0.79 2547.84152 1.09 3735.887190.2 205.544967 0.5 1273.19459 0.8 2591.64977 1.1 3767.808
37、150.21 230.975368 0.51 1315.52795 0.81 2635.31458 1.11 3798.878940.22 257.504134 0.52 1358.09858 0.82 2678.82244 1.12 3829.022810.23 285.055635 0.53 1400.89199 0.83 2722.15963 1.13 3858.140220.24 313.564935 0.54 1443.89392 0.84 2765.31226 1.14 3886.093990.25 342.975370 0.55 1487.09033 0.85 2808.2662
38、4 1.15 3912.679050.26 373.23671 0.56 1530.46738 0.86 2851.00724 1.16 3937.553040.27 404.303829 0.57 1574.01138 0.87 2893.52069 1.17 3960.047140.28 436.135744 0.58 1617.70883 0.88 2935.79175 1.18 3978.491370.29 468.694892 0.59 1661.54636 0.89 2977.80527 1.19 3986.167520.30 501.946577 0.6 1705.51073 0
39、.9 3019.54577 1.2 3992.392645.2.3 模型的检验在计算误差时发现其值基本上是保持在一个数值上，因此可以说明这个理论函数是可用性较高的。六、对问题二求解的实际储油罐容积积分模型6.1 模型的建立对于实际储油罐的特殊形态，为了求其中油面对应的容积，要把储油罐分解成两个球冠与一个圆柱，分别求出各自的容积后求和。其中从正视图分解的各部分如图 13图 13 储油罐截面分解图 图 14 左视图6.1.1 中间圆柱的容积积分模型xyhR-ho14当无变位时，以球冠的中心点为原点建立平面直角坐标系，而实际容积积分模型要根据储油罐的油高分为两类：（1）油面高度 的情况hR当 时，储
40、油罐中间圆柱体的左视图见图 14。为了避免油高出现负值给积分和解释带来不便，我们将 区域，即把图 14 中阴影区域翻折到 时的阴影区域，0y 0y则可求出阴影面积为，222arcsinRh RhSdRh 因此中间这个圆柱容积就可以求出为。 (13)22rsiVL L 柱（2）油面高度 的情况hR当 时，中间圆柱体的左视图中阴影面积为22hSyd，2234arcsinRhh0 xyr图 15 求半径图 图 16 球冠图因此对应容积为。 (14)22234arcsinRhVSLRhRhL 柱6.1.2 两端球冠的容积积分模型如图 15 所示，先求出球冠的半径 ，即r， (15) 22rsRr之后以
41、垂直于地平面为截面画平面图，其中以右弧与中间长方形接线处中点为圆心建立平面直角坐标系。则截球冠的圆弧所在的方程为，22xrsyr在立体图中球冠任意一截面所截得的圆半径为xy015， (16)22yrxs这个就是油面高度。而作为两端的球冠容积也要分成两种情况讨论：（1）油面高度 的情况R若 ，则按式(13)有y2 2222 2SxRrxsRrxsrxs，222arcsinrxsR则球冠的容积为 。 (17)0sVSxdz冠（2）油面高度 的情况y若 ，则按(14)式有yR2 2 222 2 234SxrxsRrxsRrxs,22arcsinRrxs则球冠的容积为 。 (18)0sVSxdz冠综上
42、所述，得储油罐的容积积分模型为或者 。2冠柱 2V冠柱6.2 模型的求解先计算中间圆柱体容积，即已知 ，代入(13)式和(14)式，得8L2 21.581.5.1.5.arcsinhVShh 柱2 2 21.53.4.arcsinh柱16对于球冠，根据附录一图 1 中的数据，可得 ， ，代入式(15)得出半径1s3R。再代入式(16),得出1.625r22.0.65yx110 2 222 21arctnhoVSdHRhRHrhRrh d冠 2231111211323123211 21312()63arctnrt2arctn2hRhRhrHhRHhRrhhHrR17220 2 222 22232
43、222323221arctn()63arctnhoVSdHRHhRrhRrh dhRrhRHRhRhrHHhRrH 冠 3222 2322rtarctnhRhrhRHrh而 和 中的 、 和 的关系如下1Vh21h2， ，tanltanl代入亦可得出 、 。所以，纵向变为 角度时，油浮子测量12V高度所对应的油的理论容积值为：。 (19)(tan()tan( ) 321 hVlhlhV6.3 带倾斜角度的容积积分模型已知原来的空间直角坐标系是以罐长方向为 轴，垂直地平面为 轴，指向纸外为zy轴的坐标系。在上面无变位的容积积分模型基础上，将空间直角坐标系作变换如下x， (20)cosxyz即由
44、建立了新的沿罐的空间直角坐标系。之后将无变位的容积公式中的坐标用,x新坐标系的坐标表示，即得出容积关于两个倾斜角度和高度的函数关系 。由,Vh(20)式的变换得出 ，代入(19)式即cosh18。 (21)11223,(costan)(costan)(cos)VhlVhlVh又由附录二的数据，拟合曲线发现倾斜角度的值不是很大。因此我们规定一个准则，即 。令 是一个固定值，就可以得出 。分别将2,5,0,的上下限代入 中，得到的容积量分别为 。而附录二数据中对应高度,1,的容积是 。若 ，则取 。依次这样计算，即是最小二hV11hhV2hV乘法的思想。通过从两边逐步逼近的过程，可以近似求出 的值
45、，使得 的,Vh数值最接近 。通过 Matlab7.0 软件编程，求得最优的倾斜角度为h。2.3,.1最后，将最优的倾斜角度都代入(21)式，令 且每隔 0.10m 递增。这0.1,3h样就得出了倾斜时的罐容表，如下表 8表 8 有倾斜角度的罐容表高度（m） 油量（L） 高度（m） 油量（L） 高度（m） 油量（L）0.1 0.34568 1 2.10.2 0.43505 1.1 2.12101 2.2 4.263200.3 0.51452 1.2 2.23453 2.3 4.543780.4 0.67498 1.3 2.67534 2.4 4.878670.5 0.77546 1.4 2.8
46、3412 2.5 5.534210.6 0.86512 1.5 3.05154 2.6 6.351560.7 0.96965 1.6 3.27410 2.7 7.835410.8 1.05578 1.7 3.46785 2.8 8.905560.9 1.38756 1.8 3.68615 2.9 11.534781.0 1.86435 1.9 3.94541 3 13.574386.4 模型的评价对于所给的附录二实验数据进行拟合等，我们发现油浮子的质地与油浮子种类有关。当油浮子是指针型时，影响较小。但油浮子若连接着压力传感器，就与油的深度有关。又当外界气压和温度变化较大时，对储油罐内部的气体蒸发也有影响，使得蒸发的气体积压在储油罐上方，从而导致油高的不准确。19参考文献1 华东师范大学数学系， 数学分析 ，北京
