A题《储油罐的变位识别与罐容表标定》论文解析.docx
《A题《储油罐的变位识别与罐容表标定》论文解析.docx》由会员分享，可在线阅读，更多相关《A题《储油罐的变位识别与罐容表标定》论文解析.docx（20页珍藏版）》请在冰豆网上搜索。
A题《储油罐的变位识别与罐容表标定》论文解析
储油罐的变位识别与罐容表标定
摘要
为了解决储油罐的变位识别与罐容表标定的问题，我们建立了积分方程模型。
从计算油面高度与储油体积的函数关系出发，考虑有无变位及纵向、水平偏转对罐容表测量的影响，建立更加合理的标定罐容表模型。
问题
（1），对于无变位且几何外形规则的小椭圆型储油罐，首先通过几何及积分运算，得到罐内储油量的函数表达式，建立模型一，然后对求得的模型进行误差分析和修正。
当=4.1时，在无变位情况下修正模型一的基础上，建立模型二，利用附件1中的实验数据，计算模型二的误差值，从而再次修正由积分得到的理论模型。
然后用实际数据做出检验，得到修正后的相对误差，说明模型的修正取得了很好的效果。
利用修正后的模型，计算并给出罐体变位后油位高度间隔为1cm的罐容表标定值（详见表1），取部分数值如下：
高度/cm
油量/L
高度/cm
油量/L
高度/cm
油量/L
高度/cm
油量/L
1
1.2
31
611.5
61
1805.2
91
3057.5
2
3.4
32
646.4
62
1847.9
92
3096.1
3
6.5
33
681.7
63
1890.7
93
3134.4
问题
（2），对于两端是球冠体，中间是圆柱体的储油罐，首先考虑无变位情况下实际储油罐两端球冠中的油量体积，继而考虑纵向倾斜角度的情况，可以通过油罐中倾斜液面的平均高度简化为水平液面的液面高度近似计算。
而中间的圆柱体，在模型二的计算基础上进行了模型的改进。
将三部分的油量体积值相加，可以得出储油罐内油量关于纵向倾斜角度和油位高度的函数表达式。
由几何关系可知储油罐内实际油位高度与显示油位高度和横向偏转角度有关，通过模型六的建立，求解出实际油位高度与显示油位高度的转化关系式，从而对罐容表重新进行标注。
再利用附件2的数据，求得出油量相对误差的平均值，以平均相对误差最小作为目标函数，用matlab编程优化出变位参数。
最后将变位参数代入模型，计算并给出罐体变位后油位高度间隔为10cm的罐容表标定值。
关键词：
罐容表标注变位识别枚举法
2问题重述
生活中常用罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。
储油罐常在使用一段时间后，发生罐体的位置发生变位，从而导致罐容表发生改变。
按照规定，需要定期对罐容表进行重新标定。
图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。
图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。
请用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。
（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。
请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。
（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度和横向偏转角度）之间的一般关系。
请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。
进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。
3问题分析
通过建立数学模型，求得罐内油量关于纵向倾斜角、横向倾斜角及油位高度的函数关系，解决罐容表的重新标定。
利用罐体变位后在进/出油过程中的实际检测数据，根据数学模型确定变位参数，解决储油罐的变位识别问题。
3.1问题一的分析
当油罐无变位时，油的体积为油的侧面积与油罐长度之积。
当油罐纵向倾斜角时，油面相对于水平面的高度为，油量的计算将随着的变化而改变。
当和将油量计算分成了三种情况。
由于探针、注油管和出油管具有一点的体积，以及随着油面的升高，油的压强的变化，会使计算量与测量量之间存在误差，综合考虑各方面的因素，可以达到误差与油面高度的关系，在可调节误差范围内，使得计算值更接近于测量值。
3.2问题二的分析：
实际的油罐与实验用的油罐主要区别在于多出两端球冠体。
将油罐的体积分为中间圆柱体和两端球冠体三个部分来计算。
圆柱体的体积计算与问题一的计算方式类似。
两端的球冠体将根据油面相对于水平面的高度的不同，计算球冠中油的体积、高度与纵向倾斜角、横向倾斜角的关系。
将球冠体积油量与圆柱体油量相加得到总的油量与高度，纵向倾斜角，横向倾斜角，根据所给的测量数据，代入不同的，得到多个的关系式，从而可以解出。
即可得到油量与的关系，取不同的值，对罐容表重新标注。
4条件假设与符号说明
4.1问题的假设
（一）假设问题
（2）中的；
（二）假设对球冠体积的近似计算产生的误差对结果不造成太大影响
（三）假设假设油面最高达到探针口
4.2符号说明
符号
含义
油面距水平高度
纵向偏移角（弧度）
横向偏移角（弧度）
椭圆或圆的弓形面积（）
椭圆长半轴长
椭圆短半轴长
储油罐的底面长度
储油体积（）
左球冠油量体积
右球冠油量体积
问题二圆柱半径，常量
邻界高度
用于说明某些问题提出的高度
用于说明某些问题提出的圆的半径，变量
5模型的建立与求解
5.1问题一模型的建立与求解
5.1.1模型一
当油罐无变位时，油罐的侧面图如图1所示，建立如图所示坐标系。
图1油罐侧面图
当油面高度为h时，侧面图上阴影部分面积表达式：
可求得储油体积表达式：
利用所建模型，并结合附表1中油位高度的数据，通过matlab编程，得到理论计算罐容量值。
将计算得到的预测值与附表中给出的测量值，代入matlab中，比较结果如图2：
图2无变位时预测值与测量值
预测值与测量值的差值，可得如图3：
图3差值（进油量/L×油位高度/mm）
模型一的误差分析：
测量值与计算值之间的误差值与高度值间的函数关系式近似为一条直线，查阅相关资料。
探针、注油管和出油管具有一点的体积。
探针的直接大约为，而注油管与出油管一般为直接，管壁厚为的空心圆柱体。
同时，随着油量高度的上升，油罐内部的温度、压强等因素将发生改变，油的实际体积也会受到影响。
综合考虑各因素，给出上述因素对计算值的影响关系为。
模型一的改进：
根据上述误差的分析，对模型进行改进，考虑温度压强，已经管道自身的体积对油量的影响，通过计算改良后的预测值与测量值的差值，可得如图4：
图4差值（进油量/L×油位高度/mm）
由相对误差，可得图5：
图5相对误差
经过改进，计算值与测量值更加接近，相对误差。
5.1.2模型二
油罐纵向变位且变位角度为。
探究罐体变位后对罐容表的影响，根据变位后油面可能处于不同的位置，考虑到当油面超过探针口时，随着油面的上升，探针已经失去了显示的意义，因此可确定三个端点、、（如图6）分三种情况进行求解。
图6油罐拋面图
（一）当油面水平高度小于时：
（如图7）。
图7油罐剖面图
将图7中油量的剖面单独取出分析，建立如图8所示直角坐标系，根据几何关系，可做出图8中的边长表达式，设在0y方向上，y点处与0z方向平行高度为h’。
图8油部结构坐标系
当Y=y时：
可求出储油体积：
表达式中的范围为：
（二）当油面处于水平高度小于时：
如图9所示：
图9油罐剖面图
与
（一）同理，建立的y0z直角坐标系（如图10），根据图9的角度、长度表达式，可做出图10中的边长表达式。
设在0y方向上，y点处与0z方向平行高度为h’。
图10油部结构坐标系
当Y=y时：
表达式中h的可行范围为：
（三）当油面水平高度大于小于时：
如图11所示：
图11油罐拋面图
取图11油罐无油的剖面，如图12建立直角坐标系，根据图11的角度、长度表达式，可做出图12中的边长表达式，设在0y方向上，y点处与0z方向平行高度为h’。
图12油部结构坐标系
可求出高度表达式：
代入图1中的表达式：
可求出油面体积：
表达式中h的可行范围为：
模型求解：
模型中所使用的高度为油面与水平面之间的高度，而实际标注的为探针所测得的油面高度，事实上，与存在着线性关系，如图13：
图13油罐剖面图
根据几何关系可以求出。
利用所建模型，并结合附表1中油位高度的数据，考虑模型一中提出的误差，通过matlab编程，得到理论计算罐容量值。
比较预测值与测量值，可得到如下折线图：
图14有变位时预测值与测量值（进油量/L×油位高度/mm）
通过图14可以发现，利用模型求出的油高和测量值非常吻合，因此可以利用模型建立新的罐容表如表一：
表1罐体变位后的罐容表标定值
高度/cm
油量/L
高度/cm
油量/L
高度/cm
油量/L
高度/cm
油量/L
1
1.2
31
611.5
61
1805.2
91
3057.5
2
3.4
32
646.4
62
1847.9
92
3096.1
3
6.5
33
681.7
63
1890.7
93
3134.4
4
10.6
34
717.5
64
1933.5
94
3172.3
5
16.0
35
753.8
65
1976.4
95
3209.8
6
22.5
36
790.6
66
2019.2
96
3246.9
7
30.4
37
827.7
67
2062.1
97
3283.5
8
39.6
38
865.3
68
2104.9
98
3319.8
9
50.3
39
903.3
69
2147.7
99
3355.5
10
62.4
40
941.6
70
2190.5
100
3390.8
11
76.1
41
980.3
71
2233.3
101
3425.5
12
91.4
42
1019.4
72
2275.9
102
3459.7
13
108.3
43
1058.7
73
2318.5
103
3493.4
14
126.8
44
1098.4
74
2361.0
104
3526.4
15
148.8
45
1138.3
75
2403.4
105
3558.9
16
170.6
46
1178.5
76
2445.7
106
3590.7
17
193.8
47
1219.0
77
2487.9
107
3621.8
18
218.1
48
1259.7
78
2529.9
108
3652.2
19
243.5
49
1300.7
79
2571.8
109
3681.8
20
269.8
50
1341.9
80
2613.5
110
3710.7
21
297.1
51
1383.2
81
2655.0
111
3738.7
22
325.3
52
1424.8
82
2696.3
112
3765.9
升级会员 