储油罐的游标设计.docVIP

实验十三 标尺的刻度设计 13.1 实验目的 本实验涉及微积分和线性代数的若干基本知识。通过实验复习函数的反函数、定积分的应用和计算、解微分方程和解代数方程组等内容；并且介绍了连续型数学问题近似化处理的若干基本技巧，例如插值法、二分法和两点边值问题的差分法等。 13.2 实际问题 在石油的生产地和加工厂，为了储存原油，经常使用大量的储油罐。油罐的外形为一个圆柱体和两个圆锥体的组合，上端有一注油孔（见图13.1）。由于经常注油和取油，有时很难知道油罐中剩油的数量。这给现有储油量的统计带来很大的麻烦。显然，将剩油取出计量是不现实的。因此，希望能设计一个精细的标尺：工人只需将该尺垂直插入使尺端至油罐的最底部，就可以根据标尺上的油痕位置的刻度获知剩油量的多少。这是一个来自油田的问题。 图13.1 13.3 数学模型 设圆柱的底面半径为，长度为；而圆锥的底面半径也是，高为。若标尺被油浸湿位置的高度为，而此时罐内的油量为。那么我们的问题归结为求出函数 （13.1） 其中 （13.2） 换言之，要求油量函数 （13.3） 的反函数。由于对称性，我们只要研究的情形就行了。 13.4 问题的解法 Ⅰ．求油量函数的解析方法 记 （13.4） 其中和分别为相应圆柱和圆锥（一侧部分）中的储油量。 先求。设圆柱体截面中储油部分对应的弓形区域面积为，弓形对应的圆心角的一半为（见图13.2）， 图13.2 那么易得 利用的三角函数表达式，就有 （13.5） 于是 （13.6） 再求。设与圆锥体底面平行且距底面处的截面上表示储油部分的弓形区域面积为，那么和前面不同的是：不仅与有关，而且与有关。设该弓形的半径为高为，则由几何关系（见图13.3）不难看出： 从而得到 （13.7） 图13.3 得 （13.8） 于是类似于的求法，有 通过作变换 即 （13.9） 将有 （13.10） 这个积分虽然比较复杂，但还是可以积出来，建议读者自己进行计算，也可以借助数学软件。具体结果为 从而可知 （13.11） 综合（13.4）、（13.6）和（13.11）三式，就有 （13.12） Ⅱ．刻度位置函数的求法 显然，要由油量函数的解析表达式（13.12）来解析地求出反函数的显式表示是不可能的。这里，我们介绍两种近似方法。 1．插值法 将区间[0，]作如上等分： ，那么由式（13.12）可以求得在各分点的值： （13.13） 现在的问题是：如果给出油量，如何确定（求出）相应的标尺刻度位置？而这实际上是求方程 的根。 （1）线性插值 设，在此区间上用线性函数近似，则有 （13.14） （2）二次插值 先找出最接近的的，然后在用二次函近似，则有 （13.15） 插值方法还有多种，这里介绍的是最简单的。希望读者由此了解插值的思想。 2．二分法 对给定的，先估计的范围，例如可以从函数在等分点的值（13.13）来判断。假设，当然。取，由式（13.12）求出相应的。若，则偏大，取；若，则偏小，取，；再取，由式（13.12）求出相应的。依照上述原则，可求得一系列的；直至；（是给定的误差界）。此时可取 综上所述，可以将标尺刻度设计的过程归纳为： 步骤一 把区间[0，]依精度需要分成等份，得剖分点，，；其中。利用公式（13.12）算出刻度位置处相应的油量。 步骤二 根据所需的各油量读数所在上一步骤中的某一油量区间，通过插值或二分逼近法获得相应的刻度位置值。 13.5 数值实例 现在我们以（单位：）为例进行实际计算。 （1）求油量函数的值 取，由油量表达式（13.12）可得计算结果如表13.1所示。 表13.1 0 0.0 0.000 0 0.000 0 0.000 0 1 0.1 0.293 6 0.002 3 0.298 2 2 0.2 0.817 5 0.012 8 0.843 1 3 0.3 1.477 5 0.034 3 1.546 1 4 0.4 2.236 5 0.068 2 2.372 9 5 0.5 3.070