二阶微分方程及其模型幻灯片.ppt

过去一段时间，美国原子能委员会为了处理浓缩的放射性废物，他们把废物密封在圆桶后，扔到水深100米以上的海中。为此，许多科学家表示担心，特别是圆桶在运动过程中与海底碰撞，是否会因速度过快而破裂，从而导致放射性废物对大海的污染。这种担心是否会发生呢？首先需要描述圆桶在深海中运动的过程。 请建立描述该运动过程的数学模型。 问题的提出 模型一 无阻力运动 假设：圆桶运动过程中，所受阻力忽略不计，圆桶在水中进行直线运动，不产生旋转运动，即圆桶只受重力G和浮力e影响。 在上述假设条件下，圆桶所 受合力F=G-e。 由牛顿第二定律：F=ma 如图，设圆桶入海时记t=0，位移s=0。若s(t)为t时刻对应的位置函数，则a(t)=s’’(t)。 由此可得微分方程（1）： m s’’(t)= mg-pgV 且满足 s(0)=0, v(0)=v0 . 模型二 阻力运动 假设：圆桶运动过程中，若圆桶除受重力G和浮力e以外，还受阻力 f影响(和速度成正比) ，其它情况同模型一。 在此假设条件下，圆桶所受合力F=G-e-f。 由牛顿第二定律：F=ma 如图，设圆桶入海时计t=0，位移s=0。若s(t)为t时刻对应的位置函数，则 v(t)=s’(t), a(t)=s’’(t)。 由此可得微分方程（2）： m s’’(t)= mg-pgV-ks’(t) 且满足s(0)=0, v(0)=v0 . 模型求解并分析 假设我们已测得下面数据，试利用所建立的模型，求解微分方程并分析结果，从而判断圆桶是否会因碰撞海底而发生破裂。 测得圆桶质量m≈240kg，体积V ≈0.2立方米，海水比重p ≈ 1025kg/立方米。 通过大量科学试验表明：当圆桶的速度超过12m/s，与海底碰撞就容易破裂；另外，当圆桶在水中运动时，所受阻力和速度成正比，比例常数k约为0.08g~0.12 g （牛 米/秒），即f= kv。 模型一的求解 s’’(t)= g-pgV/m （1） 且满足s(0)=0, v(0)=v0 利用积分，可求解方程(1) s’(t)=Kt+c，K=g-pgV/m. 由v(0)=s’(0)=v0, v(t)=s’(t)=Kt+v0, 再次积分， s(t)=Kt2/2 +v0 t+ c’ 由s(0)=0, c’=0, 即解为 s(t)= Kt2/2 +v0 t，K≈1.43 s(t)=0.72t2 +v0 t. 设v0=0， s(t)=0.72 t2 . 当s=100米时，t =11.8秒， 而v(t)=s’(t)=1.43 t v=1.43*11.8=16.8 m/s. 即当圆桶到达100处的海底时，其速度将超过12m/s，因而容易破裂。 由速度函数v(t)=s’(t)=Kt+v0 ，海底愈深，其速度愈大。 若v0≠0，可得圆桶到达海底的速度更大，也就愈容易破裂。 模型二的求解 s’’(t)= mg-pgV-k s’(t) (2) 且满足s(0)=0, v(0)=0. 利用降阶法，设 v=s’(t) v’=mg-pg V-k v 再利用变量分离法或常数变易法，可求得： v(t)=s’(t)= [1-exp(-kt/m)] *(m-pV)g/k 再求积分，s(t)=…….。 同模型一，先求出到达海底的时刻t0，然后求得v(t0)与12m/s比较， 从而得到问题的答案。 不难发现，在求解t0时刻时，将会遇到求解超越方程的困难（可利用计算机求解，略）。 由模型可知，在有阻力的情况下，取k=0.12g, v(t)=290 e-0.005t，当t→∞时， v→ 290，其变化趋势如图所示： 模型二的改进 事实上，对原问题的判断，我们仅须知道速度和位置间的关系，即速度关于位置的函数v=v(s)，为此我们建立该函数的微分方程。 由于dv/dt = dv/ds * ds/dt=v dv/ds，代入方程（2）得： v dv/ds = g-pVg/m-k/m v，且满足v(0)=0. 利用变量分离法或常数变易法，可得 s= -m v/k - mg (m-p V)/k2 ln[(m-p V- k v /g)/(m - p V)] 当v=12 m/s，代入上式s =52m. 由于v=v(s)是增函数，知s=1