《数学建模讲义3.ppt

5. 放射性核废料处理问题 以前美国原子能委员会将放射性核废料装在封闭的圆桶里，扔到水深为91米的海底。生态学家和科学家担心这种做法不安全，而原子能委员会向他们保证，圆桶绝不会破漏。经过周密的试验，证明了圆桶的封闭性确实很好，但工程师们提出疑问，圆桶是否会因与海底相撞而发生破裂，原子能委员会却说决不会。 数学与信息科学学院 于是，工程师们通过大量试验发现：当圆桶在海水中的速度超过 12.2 米时，圆桶就会因碰撞破裂。现在的问题是，圆桶到达91米的海底时的速度是不是超过了 12.2 米的问题。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 首先来建立直角坐标系。设海平面为 x 轴，y 轴的方向向下为正。 数学与信息科学学院 根据牛顿第二定律：F=ma，其中圆桶的质量为W=239.456 kg，a 为圆桶下沉的加速度，假设 t 时刻圆桶下沉到 y 处，则 。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 圆桶在向海底沉的过程中，受到三个力的作用： 一是重力，W=239.456 kg； 二是圆桶受海水的浮力 B，已知海水的浮力为1025.94 kg/m3，圆桶的体积为V=0.208 m3，那么，B=1025.94×V=213.396 kg ； 数学与信息科学学院 三是圆桶下沉过程所受的阻力 D，工程师们做了大量试验得出结论：阻力与圆桶的方位无关，而与下沉的速度成正比，比例系数 k=0.12，于是， Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 数学与信息科学学院 根据牛顿第二定律得微分方程： 其中 为重力加速度。方程（1）是 t 为自变量（但不显含 t），y 为未知函数的二阶微分方程，我们要求的是当圆桶下沉到 91 米的海底时圆桶下降的速度，因此，把速度 v 看成 y 的函数 v=v (y)。因为 ，所以 ，于是方程（1）化为一阶微分方程 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 数学与信息科学学院 或 求解微分方程（3）得其通解 将初始条件 v=v(0) 代入（4）式得 从而，方程（3）的特解为 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 数学与信息科学学院 算到这里，我们似乎感到失望，因为不能从（5）式中解出显函数 v=v (y)。但是，这个困难是可以克服的。借助于计算机很容易计算出 v(91)，我们只需要为计算机提供一个 v(91)的很好的近似表达式，近似表达式可按下述方法获得。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 在微分方程（3）中，暂时设 k=0 （无阻力），从而得新的微分方程 为了区别于方程（3）的未知函数，方程（6）用u表示未知函数。直接对（6）式积分得 数学与信息科学学院 整理上式得 特别有 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 现在来说明 u(91) 就是 v(9