储油罐的变位识别与罐容表标定学生数学建模竞赛A题获奖

1、本文通过对储油罐中油位高度及变位参数之间的不同情形的储油量进行分 析并建立相应的数学模型，在该过程中先利用投影法、截面法及微元法得出储油 量与油位高度及变位参数的函数关系。再由 Matlab编程可知各高度储油量的理 论数据，最后分析误差及评价模型的合理性。对于问题一的任一种情形，我们均建立笛卡尔坐标系，当储油罐无变位时， 利用微元法得到体积关于h的公式，当储油罐发生变位时，根据储油罐中油量的 多少分成三种情形，就每一类利用微元法得到体积关于 h的公式。代人附件1实 验数据中的高度得到储油罐中的理论油量 V。根据理论油量及实际油量得出误 差，判断误差所服从的分布，再利用相对误差进行误差分析并评价

2、模型的合理性。 由上述得到储油罐发生变位时体积关于 h的公式我们给出了罐体变位后油位高 度间隔为1cm的罐容表标定值（即进/出油量与罐内油位高度的表格）。对于问题二中的储油罐，我们先将问题进行简化考虑，得出了储油罐水平卧 放时油量与浮油子高度的函数关系；再考虑储油量与油位高度及变位参数（纵向 倾斜角度a和横向偏转角度P）的一般情况，在该过程中，我们进行近似处理， 利用投影法和截面法得出了储油量关于油位高度及变位参数的函数关系；并在周定的横向偏转角度P条件下，就纵向倾斜角度a的变化进行分成三类讨论，这三 类又可以分成八种情形，得到了每一种情形下实际储油罐罐内储油量与油位高度 的函数关系。在模型的

3、改进中，我们就问题二储油量与油位高度及变位参数的一般情况进 行了仔细的考虑，将含油部分的体积分成四个部分，每一个部分将上述所提到的 积分方法相结合，得到了各个部分的储油量与油位高度及变位参数的函数关系， 从而可得总储油量与油位高度及变位参数的函数关系；并据此利用Matlab编程和实际测量的数据求得口和P值；与此同时我们可以得出在固定a、P值时各高 度下的理论储油量；根据理论油量及实际油量得出误差， 判断误差所服从的分布 再利用相对误差进行误差分析并评价模型的合理性。由上述得到储油罐发生变位 时体积关于h的公式我们给出了罐体变位后油位高度间隔为10cm的罐容表标定值。【关键词】投影法 截面法 微

4、元法 Matlab编程§ 1 问题重述通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的 “油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等 数据，通过预先标定的罐容表(即罐内油位高度与储油量的对应关系) 进行实时 计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生 纵向倾斜和横向偏转等变化(以下称为变位)，从而导致罐容表发生改变。按照 有关规定，需要定期对罐容表进行重新标定。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。(1)为了掌握罐体变位后对罐容表的影响，利用如图

5、4的小椭圆型储油罐 (两端平头的椭圆柱体)，分别对罐体无变位和倾斜角为口 =4.1。的纵向变位两种 情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容 表的影响，并给出罐体变位后油位高度间隔为 1cm的罐容表标定值。(2)对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模 型，即罐内储油量与油位高度及变位参数(纵向倾斜角度 口和横向偏转角度P) 之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据(附件2), 根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为 10cm的罐容表标定值。进一步利用附件 2中的实际检测数据来分析检验你们模

6、型的正确性与方法的可靠性。(以上涉及的图14均在附录中)§ 2模型的假设与符号的约定§ 2.1 型的假设与说明(1)在储油罐倾斜的情况下，忽略油浮子高度为 0时油所占的体积;(2)在储油罐倾斜的情况下，假设当油浮子高度达到最大后不再进油；(3)油的挥发速度很慢，忽略因油的挥发而造成储油量的减少；(4)储油罐的材料为钢体，忽略因渗出油而造成储油量的减少；(5)储油罐管理妥当，不会因特殊情况而造成储油量的变化。§ 2.2 号的约定与说明V表示储油罐中油的体积；L表示储油罐圆柱体部分的长度a表示任一椭球截面的长半轴1b表示任一椭球截面的短半轴 ；ai表示油浮子在圆柱体高

7、方向上投影至两端的较小值 ;h表示油浮子到圆柱体高方向的距离；几表示储油罐接地一端油面到地面得距离；a表示纵向倾斜角度；一：表示横向倾斜角度；ho表示球冠高；bo表示球冠底半径；§ 3问题的分析§ 3.1 题一的分析当储油罐无变位时，储油罐圆柱体的接地一端为原点，以圆柱体高方向为z轴，建立笛卡尔坐标系，利用微元法得到体积关于 h的公式，代人附件1实验数 据中的高度得到储油罐中的理论油量 V。利用附件1实验数据中得到储油罐中的 实际油量，根据理论油量及实际油量就可以得出误差，判断误差所服从的分布， 利用相对误差进行误差分析。当储油罐发生变位时，以储油罐圆柱体的接地一端为原点，

8、圆柱体高方向为 z轴，建立笛卡尔坐标系。根据储油罐中油量的多少分成三类，然后就每一类利 用微元法得到体积关于h的公式，代人附件1实验数据中的高度得到储油罐中的 理论油量V。利用附件1实验数据中得到储油罐中的实际油量， 根据理论油量及 实际油量就可以得出误差，判断误差所服从的分布，利用相对误差进行误差分析。由上述得到储油罐发生变位时体积关于 h的公式可以给出罐体变位后油位 高度间隔为1cm的罐容表标定值（即进/出油量与罐内油位高度的表格）。§ 3.2 题二的分析对于实际储油罐，我们首先将问题进行简化考虑，得出了当实际储油罐水平 卧放时实际储油罐中油量与浮油子高度的函数关系；然后我们先考

9、虑实际储油罐罐内储油量与油位高度及变位参数（纵向倾斜角 度式和横向偏转角度P）的一般情况，在该过程中，我们进行近似处理，利用投 影法和截面法得出了储油量关于油位高度及变位参数的函数关系，再在固定的横向偏转角度P条件下，就纵向倾斜角度口的变化进行分类讨论，一共有三种情形， 得到了每一种情形下实际储油罐罐内储油量与油位高度的函数关系。最后我们先利用附件2中的少量实际数据得出了附件2所处状态下的纵向倾 斜角度口和横向偏转角度P ,再利用附件2中给定各高度进行代人，得到实际储 2 油罐理论的储油量，与实际储油量进行比较，求出误差及相对误差。由上述得到储油罐发生变位时体积关于 h的公式可以给出罐体变位后

10、油位 高度间隔为10cm的罐容表标定值(即进/出油量与罐内油位高度的表格)。§ 4模型的建立与求解§ 4.1 问题一§ 4.1.1 模型一当储油罐体无变位时，储油罐圆柱体的接地一端为原点，以圆柱体高方向为 z轴，建立笛卡尔坐标系，如图4-1所示，在高度为h时，利用微元法过垂直z轴 的方向做截面S(z),对S(z)关于z进行积分，得到体积关于h的公式。油位X探针水平线一5图4-122由三七=1 = x = ： . y2 a bb2 .-y dya2.2. h - b 1 , 2_.= (h b)、, 2bh - h b arcsin bb _b 2dV =S(z)d

11、zLV = 0 S(z)dz=aL (h -b) .2bh -h2 b2 arcsinh-b 1二 b2b _b 2V =aL (h -b) , 2bh -h22 h-b 12b arcsin -二bb 2利用Matlab中的命令subs代人附件1实验数据中的各高度得到储油罐中的 理论油量V。由附件1实验数据中进油量、出油量及储油罐罐内油量初值可以得到储油罐 中的实际储油量，根据理论油量及实际油量就可以得出误差。由附录中的程序youliangl,我们得到了理论储油量，误差及相对误差。进油后理论储油量与实际储油量随高度的变化规律如图4-2所示：o Oo O5 O4 400 53o o O 口 o

12、 O 口 o o O 0 o O0 5 0 5 0 53 2 2 1 1幅匣均坦旺W24&81012进油后储油量*储油理论量 储油实际量图4- 2出油后理论储油量与实际储油量随高度的变化规律如图4-3所示:储油理忆量 储油实际量*Li， . X揖产一4 5。口4000350030002500200015001000耀内油位图度睡出油后储油量的变化情况24&91012出油后储油量无变位进油和无变位出油的储油理论量和储油实际量及误差和相对误差的 数据如表4-1所示：表4- 1无变位进油无变位出油储油量实 际值储油量 理论值.（必）相对误差储油量实际 值储油量 理论值初（误差）相对误

13、差312322.910.883.37%3916.24052.8136.613.37%362374.612.633.37%3866.24001.0134.863.37%412426.414.363.37%3816.23949.3133.133.37%462478.116.133.37%3766.23897.6131.383.37%512529.917.853.37%3716.23845.8129.653.37%562581.619.613.37%3666.23794.1127.883.37%612633.421.353.37%3616.23742.3126.163.37%662685.123.0

14、83.37%3566.23690.6124.393.37%712736.824.853.37%3516.23638.8122.653.37%762788.626.583.37%3466.23587.1120.923.37%812840.328.333.37%3416.23535.4119.173.37%862892.130.063.37%3366.23483.6117.433.37%912943.831.803.37%3316.23431.9115.703.37%962995.533.543.37%3266.23380.1113.943.37%10121047.335.303.37%3216.

15、23328.4112.183.37%10621099.137.053.37%3166.23276.6110.453.37%11121150.838.813.37%3116.23224.9108.703.37%11621202.640.553.37%3066.23173.2106.963.37%12121254.342.293.37%3016.23121.4105.223.37%12621306.044.033.37%2966.23069.7103.463.37%13121357.845.773.37%2916.23017.9101.733.37%13621409.547.493.37%2866

16、.22966.2100.003.37%14121461.249.243.37%2816.22914.498.233.37%14621513.050.983.37%2766.22862.796.503.37%15121564.752.743.37%2716.22811.094.763.37%15621616.554.493.37%2666.22759.293.023.37%16121668.256.243.37%2616.22707.591.263.37%16621720.057.983.37%2566.22655.789.513.37%17121771.759.733.37%2516.2260

17、3.987.763.37%17621823.561.463.37%2466.22552.286.033.37%18121875.263.193.37%2416.22500.584.303.37%18621927.064.953.37%2366.22448.782.533.37%19121978.766.683.37%2316.22397.080.803.37%19622030.468.433.37%2266.22345.379.073.37%20122082.270.203.37%2216.22293.577.323.37%20622134.071.953.37%2166.22241.775.

18、553.37%21122185.773.673.37%2116.22190.073.843.37%21622237.475.433.37%2066.22138.372.083.37%22122289.277.163.37%2016.22086.570.323.37%22622340.978.893.37%1966.22034.868.603.37%23122392.780.673.37%1916.21983.066.853.37%2315.82396.680.783.37%1866.21931.365.123.37%2365.82448.482.543.37%1816.21879.563.35

19、3.37%2367.12449.682.563.37%1766.21827.861.613.37%2417.12501.484.343.37%1716.21776.159.883.37%2467.12553.186.053.37%1666.21724.358.123.37%2517.12604.987.823.37%1616.21672.656.373.37%2567.12656.689.533.37%1566.21620.854.643.37%2617.12708.391.283.37%1516.21569.152.873.37%26672760.093.033.37%1466.21517.

20、351.143.37%2668.82761.993.113.37%1416.21465.649.423.37%2718.82813.794.833.37%1366.21413.847.653.37%2768.82865.496.593.37%1316.21362.145.913.37%2818.82917.298.343.37%1266.21310.444.183.37%2868.82968.9100.093.37%1216.21258.642.423.37%2918.83020.7101.843.37%1166.21206.940.683.37%2968.83072.4103.583.37%

21、1116.21155.138.943.37%3018.83124.1105.313.37%1066.21103.437.193.37%3068.83175.9107.063.37%1016.21051.635.433.37%3118.83227.6108.803.37%966.2999.933.713.37%3168.83279.4110.553.37%916.2948.231.963.37%3168.93279.5110.553.37%866.2896.430.203.37%3218.93331.2112.273.37%816.2844.628.453.37%3268.93382.9114.

22、033.37%766.2792.926.713.37%3318.93434.7115.763.37%716.2741.224.993.37%3368.93486.4117.523.37%666.2689.423.253.37%3418.93538.2119.263.37%616.2637.721.513.37%3468.93589.9121.013.37%566.2585.919.743.37%3518.93641.7122.763.37%516.2534.218.023.37%3568.93693.4124.513.37%466.2482.516.283.37%3618.93745.1126

23、.233.37%416.2430.714.503.37%3668.93796.9127.983.37%366.2379.012.773.37%3718.93848.6129.743.37%316.2327.211.023.37%3768.93900.4131.483.37%266.2275.59.283.37%3818.93952.1133.233.37%3868.94003.9134.953.37%3918.94055.6136.703.37%3968.94107.4138.453.37%#由上述的表格可以得出相对误差稳定，不会随高度发生变化。当储油罐体发生变位时，我们就油面及油浮子的位置与

24、倾斜角度的变化情况 分成将该问题三类，然后就每一类进行分析，找出了储油罐中油的体积与油浮子 位置的函数关系，其立体图如图 4-4所示：图4- 4(1)当油浮子的高度在0 MhMh -ztana时，如图4-5所示，截面图如图4-6所示图4- 5图4- 6截面图截面面积为：加 -ztanOCs(z)= 12f7b2-(y-b2)dy0 0b=ai A - ztan(: ) -b 2b(h1 - ztan(二)YA - ztan(： )2 b,a '.2. 'hi -ztan(«) -b Y 1 . 2十 b arcsin +-bb < I b J 2；所以油罐体积为

25、:hlV = f(A)=S(z)dz9 ，tan ：- a h n0 b h1-ztan(.：s) - b . 2b(hi - ztan(.：O) - (h1 - ztan(.：Q)2(1.1 )tanaa .2.一 b arcsin0 bhl -ztanC >b 1-.b2b2将卜=几-atana代入公式(1.1 ),则公式(1.1 )人 h a1tan 二 一ztan 二 一b1tbV = ab tan 二(t 1 -t2 arcsint 二)dtab2h a1tan ： - b h a1 tan : - barcsin I +1/h+a1tana-b ；ab2+tan(二)h 阚

26、tan(： ) -bH I113i'h + a1 tan(ct) b ；b3/2、(2)当油浮子的高度在Itana < h c 2b-a1tana时，如图4-7所示,图4- 7截面面积为:rhi -ztanotS(z)=2a . b2 -(y-b2)dy b(0.2)=abrLh+a1tan( a )-ztan( a )-bb2h+a1tan(二)-ztan(二)-b arcsin I 1十一冗dz2 J将h=h1 -a 1tan( a)公式(0.3),4 h+a1tan(二)-ztan(二)-b v t= -Hh -ztan(： ) -b 2b(h1 -ztan(： ) - (

27、h1 - ztan(： )2QDbh1 - ztan(： ) -b arcsin 所以油罐装油的体积为;LV = 0 S z dz#公式（2.2）化为：2, (h+a1tan(a)-Ltan(a)-b)/B-ab I 1(h+atan(a)-b)/bt、1-t 2 arcsin(t) 二 2将积分下限令为p,积分上限令为q,则p= h +a tan(口)-bh+a1tan( : )-Ltan( : )-bV(h)=ab2；,1tan。：） 3-q2- 1 -q2 -qsin( q) - 1 - p2 psin( p) 1 - p23ab2tan( -)h a1 tan(二)一 b1- . b:

28、 、y J（3）当油浮子的高度Ltana<hW 2b时，如图4-8所示,图4- 813 其中,此时截面面积为:卜 h =几 一 h =几-2b a tan ：, h2 = l tan : - h.八2 -ztanaS(z) = l2Wb2-(y-b)2dyJ 0 b=a (h -(z-a1)tan(«) -b )J2b(h -(z -a1) tan(«) - b ) -(h-(z-a1)tan(«) b )2 ) ba 12.1h-(z-a)tan(a)b 12'十一 b arcsin +nbb < I b 2 2)所以油罐装油的体积为Ltan

29、()_h!；2b 口 tan( )V(h) = ( tan(a s(z)d z令1=h 一(z -a1)tan(二)一bh -(L - a1)tan(-) b2h 2a1tan(： )-3b-Ltan(：)则上式公式变为:利用Matlab中的命令subs代人附件1实验数据中的各高度得到储油罐中的 理论油量V。由附件1实验数据中进油量、出油量及储油罐罐内油量初值可以得到储油罐中的实际油量，根据理论油量及实际油量就可以得出误差由附录中的程序youliang2,我们得到了理论油量，误差及相对误差。倾斜时进油后理论储油量与实际储油量随高度的变化规律如图4-9所示:o O o OO 5 4 300刈I

30、4 OODOO o o o o O 5 0 6 0 5 2 2 11 超帼/建苣蜜储油理论量 储油实际量一 * J> * « * . *« _* * J* ' * -，二1I1i倾斜时罐内油位高度随进油后储油量的变化情况567891011倾斜时进油后储油量图4- 9倾斜时出油后理论储油量与实际储油量随高度的变化规律如图4-10所示:储油理论量 储油实际量_ >9一 .*4:一* « .* * J /*，：. :卡 西 口 口 4*IILi倾斜时罐内油位高度随出油后储油量的变化情况400035003000250020001500100050045

31、67891011倾斜时出油后储油量图 4-10将倾斜进油和倾斜变位出油的储油理论值和储油实际值及误差和相对误差 的数据如表4-2所示：表4- 2倾斜变位进油倾斜变位出油储油量实储油量 理论值鸟（误差）相对误差储油量实际 值储油量理 论值鸟（误差）相对误 差962.91010.07.194.67%3464.73523.158.371.66%1012.91058.345.474.30%3414.73478.163.351.82%1062.91118.055.194.94%3364.73430.465.671.91%1112.91167.554.674.68%3314.73382.067.281.9

32、9%1162.91222.159.294.85%3264.73330.966.191.99%1212.91279.266.345.19%3214.73289.474.692.27%1262.91327.564.634.87%3164.73234.569.772.16%1312.81382.669.825.05%3114.73188.974.172.33%1362.81433.670.814.94%3064.73137.372.572.31%1412.71485.072.234.86%3014.73088.974.112.40%1462.71536.173.354.78%2964.73040.1

33、75.402.48%1512.71591.979.144.97%2914.72989.274.482.49%1562.71645.182.425.01%2864.72947.782.982.82%1612.71697.484.674.99%2814.72891.176.362.64%1662.71749.987.164.98%2764.72845.480.622.83%1712.71796.683.894.67%2714.72794.079.312.84%1762.71848.886.084.66%2664.72747.582.723.01%1812.71900.087.324.60%2614

34、.72695.881.073.01%1862.71952.790.014.61%2564.72644.579.783.02%1912.71999.286.484.33%2514.72598.183.333.21%1962.72052.789.974.38%2464.72547.182.383.23%2012.72103.290.504.30%2414.72500.285.473.42%2062.72148.485.723.99%2364.72451.987.113.55%2112.72202.389.544.07%2314.72402.988.143.67%2162.72252.790.014

35、.00%2264.72349.784.973.62%2212.72303.791.013.95%2214.72298.884.113.66%2262.72352.289.443.80%2164.72256.091.294.05%2312.72402.790.023.75%2114.72204.489.664.07%2362.72452.289.513.65%2064.72156.691.834.26%2412.72497.985.173.41%2014.72106.791.974.37%2462.72548.785.933.37%1964.72053.188.354.30%2512.72599

36、.686.833.34%1914.72005.690.814.53%2562.72644.781.963.10%1864.71953.588.784.54%2612.72696.383.553.10%1814.71905.190.394.74%2662.72748.886.113.13%1764.71852.888.024.75%2712.72798.986.143.08%1714.71801.586.774.82%2762.72847.284.492.97%1664.71752.787.995.02%2812.72895.983.132.87%1614.71699.384.554.98%28

37、62.72942.179.342.70%1564.71645.981.184.93%2912.72992.679.842.67%1514.71602.087.275.45%2962.73044.281.482.68%1464.71544.179.355.14%3012.73089.676.912.49%1414.71492.677.895.22%3062.73140.878.102.49%1364.71444.179.415.50%3112.73189.776.992.41%1314.71392.377.565.57%3162.73234.071.322.21%1264.71333.468.6

38、25.15%3212.73284.171.392.17%1214.71280.565.785.14%3262.73334.671.882.16%1164.71230.665.915.36%3312.73379.066.291.96%1114.71175.160.375.14%3362.73423.560.821.78%1064.71124.759.975.33%3412.73473.260.461.74%1014.71067.853.104.97%3462.73517.754.931.56%964.71011.847.054.65%3512.73569.556.731.59%3514.7357

39、3.258.491.64%由上述的程序还可以描绘出倾斜时储油量的误差随罐内油位高度的变化情况，如图4-11所示:OOJ40倾斜时进油后储油量的误差随罐内油位高度的变化情况567 S 91011罐内油位高度倾斜时出油后储油量南误差随罐内油位高度的变化情况4 o o O 8 6 4111:*>， ! . !*:* *；* !-总* i 10 *1« s iBt一麟1=*, ： -一 40 一 *1；iiii567891011罐内油位高度806G图 4- 11由上述图我们可以知道倾斜时不管是进油后储油量的误差还是出油后储油 量的误差均随罐内油位高度的变化呈正态分布，说明该模型建立得贴

40、合实际。由上述得到储油罐发生变位时体积关于 h的公式我们可以给出罐体变位后 油位高度间隔为1cm的罐容表标定值（即进/出油量与罐内油位高度的表格），如 表4-3所示：表4- 3罐体变位后油位高度间隔为1cm的罐容表标定值高度（h）cm体积（v）L高度（h）体积（v）高度（h）体积（v）高度（h）体积（v）cmLcmLcmL13.53126461.5511413.9762491.326.26352749452145677253439.974828527.1531498.4782576.6414.756329560.9541540.9792619.1520.690830595.2551583.58

41、02661.4627.854231630.1561626.3812703.6736.316332665.6571669.2822745.5846.142433701.5581712.2832787.2957.393534738591755.3842828.71070.12735774.9601798.58528701184.396836812.2611841.8862911.112100.254137850621885.1872951.813117.747538888.2631928.5882992.314136.92339926.7641971.9893032.515157.82464096

42、5.7652015.4903072.416180.493941100.5662058.891311217204.971421044.6672102.3923151.218228.9431084.5682145.7933190.119254.94411248692189.1943228.620281.9451165.3702232.5953266.721309.8461206.2712275.8963304.422338.5471247.2722319.1973341.723368.1481.2886732362.3983378.524398.5491330.1742405.4993414.92

43、5429.7501371.9752448.41003450.7§ 4.2 问题二我们一开始考虑当储油罐水平放置时油量与浮油子高度的函数关系。当储油罐水平放置无横向与纵向倾斜时，设油高为h,所建立的空间直角建立坐标系如图4-12所示：图 4-12储油罐中储油的体积可分为两部分，圆柱体中储存油的体积与两端的球罐体 中油的体积，根据对称性，两端的球罐体中油的体积相等。用垂直于z轴 的平面截储油罐，根据截面法求体积公式可求出储油的体积。第一种情况：当0Wh至凤时171、对于圆柱体部分：截面圆的方程为x2+(y b0)2 = bo2由x2 (ybo)2 =bo2= x =、.；bo2 (y b

44、o)2截面s(z) =2( Jbo2 (y bo)2dy21y =2-bo,bo2 -(y -bo)2 b2arcsiny-bo +bo /二(h-bo) , bo2 -(h-bo)2 bo2arcsinh bo1 +-nbo lb3.，2Lo体积v1 = 0 s(z)dz=Lo(h -bo)Kbo2 -(h -bo)2Lobo2arcsinh-bo+1兀bo2Lo19电2 bo22h32、对于两端的球冠体，根据对称性，两端的球冠体储油体积相等4-13所示：R2=(R -h3)2 bo2 = R =用平面z=z截球罐体，设截得小圆半径为r1 ,所图4-14所示，根据几何关系得:21_2'

45、;2=R -(d z)=d' = , R2 - bo221= R2- R -bo22z1= . R2- R2 -bo2z截面圆的油高为h截面 S(z) =2 °ri2 -(y -r1)2dy二(h -1)1 -(h -J12acsinhh333 V2=2obo S(z)dzh=2h h3 -h3 bo积分解得:二 hR2h T11二122(h-bo)jR2 TjR2_bo2+z j-(h _b°)2dzh3Kb)+ jiR2 T. JR2 - bo2 z ? |acsinoh'店-3 2o R2- 1 R2 -bo2 z dz3-h 2h,2h3R-h2 a

46、2 h3-R-R3 arctan 3h R-h3.2hhR -h2 - h32h - bo 2R2- , R2-bo2z323卜 一3k 2R arctandzh,2h3R-h2-h32)hLo,b02 -h2-2 x R2 -b02 2R2 b02 arctanV2V2=:hR219 第二种情况：当bo W h W 2bo时对于圆柱体的部分同第一种情况L0T)+*L0Vi = 0 S(z)dz-Lo(h -bo) ,bo2 _(h-bo)2Lobo2arc sin对于球罐体体积截面圆的油高为h'h'=r1 (h-b0)'S(z) =2 0 r； -(y-ri)2dy=(h -r1) r； -(h -了r；'arc sin这种情况与第