2010年全国数学建模A题答案

1 储油罐的变位识别与罐容表标定 摘要摘要 加油站都有若干个储存燃油的地下储油罐 需要采用流量计和油位计来测 量进 出油量与罐内油位高度等数据 通过预先标定的罐容表 即罐内油位高度 与储油量的对应关系 进行实时计算 以得到罐内油位高度和储油量的变化情 况 但是许多储油罐在使用一段时间后 由于地基变形等原因 使罐体的位置 会发生纵向倾斜和横向偏转等变化 从而导致罐容表发生改变 需要定期对罐 容表进行重新标定 在求解过程中 我们对于罐体无变位 罐体产生纵向变位 罐体在水平和纵向都产生变位三种情况 利用解析几何的方式计算出体积与变 位参数之间的关系 同时应用契比雪夫多项式对体积值进行近似多项式展开用 以对标高和出油量的关系进行拟合表示 得到较为满意的效果 第一问 1 针对无变位情况 我们计算得到椭圆油罐容积表达式为 利用契比雪夫多项式方法在提高拟合精 abl v bh v bh v bh V arcsin 1 2 2 椭 度的前提下用 5 阶多项式拟合处标高和容量之间的函数关系 2 对于纵向变位 的情况 当椭圆型罐体发生变位纵向变位角度时 我们利用体积等效思 1 4 想 讲上述罐内不规则油量容积的计算转为 1 中规则油容进行计算 利用附件 1 中数据利用最小二乘拟合方法算出油位高度的真实值 继而利用拟合多项式 408 5976 H 395 774852 5322H 13 2498H 1 1361H 0 0320H 2345 椭变 V 进行间隔为 1cm 的此罐容表进行标定 得出的表标定值如下 1cm2cm3cm4cm 118cm119cm120cm 0L0L0L0L 4017 26L4050 08L4082 80L 第二问 1 利用第一问中等体积的思想 我们同样可以对纵向倾斜角度 和横向倾斜角度时进行数学模型的建立 2 在模型的建立过程中得到一个 关于浮游子高度和偏转角 以及等效高度之间的一个表达式 从而利H h 用最小二乘拟合确定变位参数 3 利用已给数据求得表达式 继而再次利用拟合拟合多项式得出间隔为 10cm tan2tan 10 hRz 值 10cm20cm30cm 280cm290cm300cm 1352 172223 853082 13 60303 8860690 9861807 31 利用附表 2 中的数据进而进行模型正确性与可靠性的检验 2 关键词 关键词 储油罐 罐容表 变位 契比雪夫公式 等效高度 最小二乘拟合 一 问题的重述一 问题的重述 加油站一般都有若干个储存燃油的地下储油罐 并且也都有与之配套的 油位计量管理系统 采用流量计和油位计来测量进 出油量与罐内油位高度 等数据 通过预先标定的罐容表 即罐内油位高度与储油量的对应关系 进行 实时计算 以得到罐内油位高度和储油量的变化情况 许多储油罐在使用一段时间后 由于地基变形等原因 使罐体的位置会发 生纵向倾斜和横向偏转等变化 以下称为变位 从而导致罐容表发生改变 按照有关规定 需要定期对罐容表进行重新标定 储油罐的主体为圆柱体 两 端为球冠体 用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题 1 为了掌握罐体变位后对罐容表的影响 利用的小椭圆型储油罐 两端 平头的椭圆柱体 分别对罐体无变位和倾斜角为的纵向变位两种情况 1 4 做了实验 实验数据已给出 请建立数学模型研究罐体变位后对罐容表的影响 并给出罐体变位后油位高度间隔为 1cm 的罐容表标定值 2 对于实际储油罐 试建立罐体变位后标定罐容表的数学模型 即罐内 储油量与油位高度及变位参数 纵向倾斜角度 和横向偏转角度 之间的一 般关系 请利用罐体变位后在进 出油过程中的实际检测数据 已给出 根据 所建立的数学模型确定变位参数 并给出罐体变位后油位高度间隔为 10cm 的罐 容表标定值 进一步利用已给出的实际检测数据来分析检验模型的正确性与方 法的可靠性 二 问题的分析二 问题的分析 综述 该题要解决的问题是储油罐的变位识别和罐容表的标定 显然这是 几何求解的应用问题 可以通过所学过知识求解 对于问题一 题目中给出一个小椭圆型的平头油罐的模型和纵向的倾斜角的 值 并且也给出了对这个罐体无变位和有纵向变位两种情况的实验数据 所以我 们可以把这个问题分成两部分来解决 首先解决当罐体无变位时罐中储油量与油 液面高度的关系 其次解决当罐体发生倾斜时罐体中油的储油量随油液面高度的 变化 然后对这两种情况求出的结果进行比较 判断出油罐变位后对罐容表的影 响 最后重新确定变位后罐容表的标定值 分析发现当油罐无变位时 问题比较 好解决 我们可以利用积分的方法求出罐中的储油量和油液面高度的关系 当油 罐变位时 首先要选好坐标系 确定倾斜后液面高度是怎样变化的 也是通过积 分的方法进行求解 但是我们也可以寻找一点 可以把它等效成无变位的情况去 解决 对于问题二 给的是一个实际储油罐 其主体为圆柱体 两端为球冠体 该问题主要是要解决罐内储油量与油位高度及变位参数 纵向倾斜角度 和横 向偏转角度 之间的一般关系 对这个问题可以看成两部分求解 首先也是 要选好坐标系 确定好油面高度和变为参数的关系 分别求解主体圆柱体和两 端球冠体中的储油量和油位高度的关系 建立出一定的数学模型 然后通过题 目给出的数据 确定数据的可用性求出模型中的变位参数 最后通过误差分析 检验模型的正确性和方法的可靠性 3 三 问题的假设三 问题的假设 1 忽略油罐中温度 压强对储油量的影响 2 假设油罐的形状是规则的 3 假设储容表的读数是无读数上的技术误差 4 假设油罐变位的角度不会太大 5 假设罐中的油分布均匀 6 假设油罐不会因为出油的变化而发生位置变化 四 模型的建立与求解四 模型的建立与求解 问题一 问题一 对于本问题 首先可以把它分解成两个问题来考虑 第一个问题是确定小 椭圆型储油罐没有变位时储油量与油面高度的关系 利用解析几何的方法确定 标高和油量之间的关系 第二个问题确定油体容量和标高及变位角之间的关系 利用等效体积替换的思想把它转化为上一问中未变位的问题进行求解 同时我 们分析发现变位数据的油高不是真实值 所以考虑利用已有读数寻求真值 依 据上述思想分别对这两个小问题建立合理的模型 1 1无变位时 无变位时 小椭圆型储油罐的横截面很明显是一个标准的椭圆 首先建立合适的坐标系 并且椭圆的图形如图 1 图 1 1 符号声明 a 椭圆罐横截椭圆面的半长轴 b 椭圆罐横截椭圆面的半长轴 h 无变位椭圆油罐中的油位高度 面积无变位时椭圆的横切面 椭 S 小椭圆油罐的母线长l xx 0 h 2b a 4 方差和 2 R 根据模型得出的数据 i V 实测数据 i V 2 2 模型建立 模型建立 椭圆的方程为 1 2 2 2 2 b by a x 1 式中 a 椭圆罐横截椭圆面的半长轴且为 0 89 米 b 椭圆罐横截椭圆面的半长轴且为 0 6 米 油面高度为 h 时 储油罐横截面的面积可求得为 dybyb b a S h 0 2 2 2 椭 2 式中 h 无变位椭圆油罐中的油位高度 可推导得 ab b bh b bh b bh s arcsin1 2 2 椭 3 我们可以设 由可以得 b bh t bh20 11 t 可求得小椭圆油罐未变位时储油量随液面高度的变化为 abltttV arcsin1 2 2 椭 4 式中 小椭圆油罐的母线长且为 2 45 米 l 在求解时我们可以利用最佳平方逼近多项式 即契比雪夫多项式 求出储 油量的近似值 根据公式 1 我们可以准确的计算出储油量 但是公式里有反 三角函数没有多项式计算方便 所以我们用多项式逼近方法解决此问题 具体 理论证明见附录二 求得 5 abltttV arcsin1 2 2 椭 5 abltTatTa n nn 1 00 2 1 2 当 k 1 时 求得近似是一个 2 阶的多项式 当 k 2 时 求得的近似 椭 V 椭 V 是一个五阶的多项式 由于拟合时考虑的误差 选择五阶多项式拟合效果较好 所以我们用当 k 2 时的五阶契比雪夫多项式逼近 椭 V 椭 V abl xxx 53 4880615 1575 16 2 6 由于储油罐的形状以及温度压强的关系这里的系数不太准确 我们可以根 据题目 附件 1 中变位进油量数据 给出的累加进油量和油位高度的数据进行 拟合 用 MATLAB 软件拟合后得到的结果是 3 9321 99 4693h76 1928h 9 2692h 0 6284h 0 0209h hV 23 45 椭 对该模型进行误差分析 误差分析的数据源 附件 1 位变位出油量数据 将数据代入到上述 5 阶拟合多项式中与已知油量数据对比验证 验证数据方差和为 2 R 1 0517 1 1 1 2 2 n i ii VV n S 7 1 1 1 2 1 2 2 n i ii n i ii VV VV R 8 由以上可知拟合的效果比较好 有变位时 椭圆油罐有纵向倾斜且倾斜角 选定的坐标系及有关参数如图 2 所 1 4 示 分析发现变位后油高数据由于变位后产生影响 而油量数据我们认为是真 实的 因此首先利用几何关系求出油量容积的解析表达式 利用最小二乘逼近 6 的方法求出真实反映油量的油高读数 然后利用拟合的方法建立油高和油量的 模型 从而解决标高问题 图 2 1 符号说明 标准差 2 S 纵向倾角 变位后液面在 Y 坐标轴上的分量H 图中 c 0 4m d 2 05m 设浮标的距离罐底的距离为 可以知道浮标的坐标为 c 直线 AB f h f h 的斜率为 所以可求得 AB 的直线方程为 tan L2 tantanchxy f 9 显然油液面的高度随横坐标 x 的变化为 tantan chxH f 10 y L2 x 0水平线 Q c d 7 此时储油罐横截面的面积为 dybyb b S H 0 22 a2 椭 11 ab b bH b bH b bH arcsin1 2 2 所以变位时小椭油罐中油容量与油位高度的关系是 l dxSV 0 椭椭 12 根据公式 2 我们可以准确的计算出不同油位高的储油量 但是 公式中 还有反三角函数和开根号 我们计算起来很不方便 所以我们可以寻找一点做 体积上的等效转换仿效上一问的方法 我们设那一点为 Q 那么它的坐标可设 为 见图 2 H 并且该点 Q 也在直线 l1 上 根据 l1 的直线方程可以求得为 H tantan chH f 13 对于一个固定外壳的容器 其中装定量体积的液体 无论如何对容器发生偏 转 其液体的体积始终是一个固定值 根据等效的原理可以肯定点 Q 可以把楔 状油形体积转化为标准的椭圆柱的体积 并且此时的油位高度为 所以这样 H 就可以转化成没变位的情况去解决问题 得到转化后的椭圆面积为 dybyb b a S H 0 22 1 2 椭 14 ab b bH b bH b bH 2 arcsin1 2 可以求得体积为 abl b bH b bH b bH V 2 1 arcsin1 2 椭 15 8 我们可以设 得到 b bH t abltttV 2 1 arcsin1 2 椭 16 得到的体积公式的形式和未变位时得到的公式形式相似 所以我们同样可 以 利用契比雪夫多项式来近似地求解储容量与油位高度的关系 求得 abltttV 53 48 80 615 1575 16 2 椭 17 题目中附录的高度是变位后的读数 为了求 我们可以用最小二乘法 H 2 1 min n i i VV椭 算出的值 讲求出来的带入中可以得到真实等价的油位高度 我们用这个真实 H H 等价的油位高度和题目给出的累加进油量的数据进行 5 阶多项式拟合 求得 H 的近似的为 椭变 V 408 5976 H 395 774852 5322H 13 2498H 1 1361H 0 0320H 2345 椭变 V 对这个模型进行误差分析根据上面的公式求得 标准差为 1602 2 2 S R 检验 1 2 R 可见拟合的效果很好 符合要求 然后根据这个模型就可以算出罐体变位后油位高度间隔为 1cm 的罐容表标 定值如下 详细结果见附录 1cm2cm3cm4cm 118cm119cm120cm 0L0L0L0L 4017 26L4050 08L4082 80L 利用附件 1 中的有变位的出油数据分别带入上述两个模型中 做出的图如 下 从图中可以看出 当油位高度一样时变位后的油罐中储油量的值比没变位 的储油量低 说明罐容表实际的读数比理论的要偏大 9 问题二问题二 该问题是问题一的变形 是研究实际储油罐变位后储容量与油位高度之间 的关系 储油罐的变位是横向和纵向都有一定的偏角 并且储油罐两端不是平 头而是球冠体 首先建立合适的坐标系 并且要设定一定的参数 如图 4 R fh h z x y 水平线 图 3 图 4 10 1 1 符号声明 符号声明 实际测量的油位高度 f h 储油罐横向偏转角度 R 圆柱截面的半径 r 球缺的半径 b 1m 球缺的顶高 柱体的体积 柱 V 储油罐的中油的体积为 冠柱 V VV 18 一 一 对储油罐圆柱体部分的计算 对储油罐圆柱体部分的计算 如图在 XOZ 坐标面内 Z 轴方向的水深为 19 RRhh f cos 1 式中 实际测量的油位高度 f h 储油罐横向偏转角度 z 1 h x 图 5 11 R 圆柱截面的半径 浮游子在 ZOY 面投影的点的坐标为 2 R f h ZOY 坐标面与油液面的交线的斜率为 将浮游子在 ZOY 坐标面投2l tan 影的点的坐标带入的方程得2l 20 2tan 1 yRhz 即求得的方程得 2l tan2tan 1 Rhyz 对于这个计算我们也可以寻找一个点 P 使得 R 变为没有变位时 0 z 0 z 油面的高度 即 tan2tan 10 Rhz 在 YOZ 平面内截圆柱的方程为 222 Ryz 它的截面积为 21 dzzRS z R 0 22 2 2 000 arctan1 2 2R R z R z R z SlV 柱 lR R z R z R z 2 000 arctan1 2 2 22 由第一问可知它可以用 5 次多项式近似代换 2 2 对储油罐两段球缺部分的计算 对储油罐两段球缺部分的计算 由球的关系有 23 2 2 2 2 R brr 式中 r 球缺的半径 b 1m 球缺的顶高 解得 r 1 625m 在 XOZ 截面内用平面 z 截到一个平面 0 z 12 面积为 24 dyyzrS zr br z 22 222 2 另 22 yrp 有dyypS p br z 22 brpbr p brpp 2 22 2 1 arcsin 242 dzSV z z 0 5 1 球 25 0 5 1 2 22 2 1 arcsin 242 z dzbrpbr p brpp 这样储油罐的体积就可以表示为 26 球柱 V VV V 也可以用多项式拟合 不过不能直接用附录里的数据 哪里的体积是不准 确的 它是没有考虑变位时对体积的影响 我们就用第一个体积 把它作为储 油罐的初始体积 分别用当前行的体积减去下一行的出油量 这样我们认为此 时的体积即为实际体积 同第一问中问题二的解决方法一样通过最小二乘法计 算出 和 求得 7 2 6 3 然后反向应用拟合出 5 次多项式 即 27 54 2 3 3 2 4 1 5 0 ahahahahahaV 实 800 8409 479 2440h 178 2758h 6 3084h 0 1206h 0 0015h 2345 对这个模型进行误差分析求得 标准差为 8848 0 2 S R 检验 1 2 R 符合要求 然后根据这个模型就可以算出罐体变位后油位高度间隔为 1cm 的罐容表标定值如下 13 显示油高 cm 显示油量容积 L 显示油高 cm 显示油量容积 L 101352 1716031568 85 202223 8517033378 42 303082 1318036168 99 403195 81