2010年全国大学生数学建模竞赛A题全国一等奖论文.pdf

数值解法在储油罐变位标定中的应用 摘要 本文主要采用微元分析法、数值积分法和回归分析法，对储油罐油位测量高度 和 H 燃油体积 建立了模型，即问题一中椭圆型储油罐无变位模型 ,椭圆型储油罐纵 V H V 11 向变位模型H V ，问题二中实际储油罐变位模型H V 。 12 2 对于问题一中的椭圆型储油罐，探究变位对罐容表的影响，先建立无变位时的 H V 模型。以椭圆型储油罐左底面中心为原点，该面为 平面，建立空间直角坐标 xoy 系。根据简单柱体的体积计算公式，得出无变位初模型。又将由模型得到的体积理论计 算值与所给实际体积数据对比，发现两者大致呈线性关系。将两者进行线性拟合以对初 模型进行修正，最终得到椭圆型储油罐无变位模型H V 。 11 当储油罐发生纵向变位后，在上面所建坐标系下，用微元分析的思想，以积分的方 法求得假设条件下燃油体积 关于油位测量高度 的解析解。但是由于解的复杂性，很 V H 难由直接积分的方法求得给定 的燃油体积。转而利用数值积分的方法，将积分转化为 H 若干微元的和，借助MATLAB 即求得燃油体积，这样就得到了纵向变位模型。将依模型 所求的理论计算值与处理得到的实际值进行对比分析，两者关系走向基本呈线性，于是 对两者进行线性拟合，以缩小误差。最终得到精度较高的修正模型H V ，以此对储油 12 罐进行重新标定，标定结果参见附录一中的图表一。 对于问题二中的实际储油罐，将其分成左球冠、中间圆柱体和右球冠三部分，参照 问题一建立的空间直角坐标系，利用微元法，分别建立了各部分燃油体积关于油位测量 高度 和参数( ,) 的积分表达式。然后将三部分的燃油体积相加，即得到储油罐内燃 H 油体积关于变量 和参数 、 的积分形式解析解 。由于积分表达式十分复 H   V (,,H ) 2 杂，依然利用MATLAB 数值积分来处理。对于未知参数 、 ，本文用回归分析法和遍   历思想求解。即以较小间隔，从 取到 ，让 和 分别遍历所有的取值点。定义一 0 10   个函数 2 ,其中 为附件 中所给的部分 数据所组 F (, ) (V (, ,h(k )) V (k )) h 2 H 2 V ( ,