四训练题2分析和总结.docx

评分：广州大学2010年数学建模暑期培训答题卡

评分：

学员姓名题目类型评语

学院及年级 组号

微分方程模型2 日期 请在表格下面书写答案

题目：

以前美国原子能委员会将放射性核废料装在封闭的圆桶里，扔到水深为91米的海底。生态学家和科学家担心这种做法不安全，而原子能委员会向他们保证，圆桶绝不会破漏。经过周密的试验，证明了圆桶的封闭性确实很好，但工程师们提出疑问，圆桶是否会因与海底相撞而发生破裂而原子能委员会却说决不会。于是，工程师们通过大量试验发现：当圆桶在海水中的速度超过12.2米时，圆桶就会因碰撞破裂。现在的问题是，圆桶到达91米的海底时的速度是不是超过了12.2米的问题。现在来计算这一问题。其中圆桶的质量（重

力）为W?239.456(N)，海水的浮力为1025.94(N/m3)，圆桶的体积为V ?0.208m3，圆桶下沉过程所受

的阻力为D，工程师们做了大量试验得出结论：阻力与圆桶的方位无关，而与下沉的速度成正比，比例系数k?0.12.

要求给出计算机仿真，拟题人：钟老师

答案： 海中物体下落问题

问题摘要问题重述

对于以前美国原子能委员会将放射性核废料装在封闭的圆桶里，扔到水深为91米的海底很多科学家对其安全问题持有怀疑，之前工程师已对圆桶经过周密的试验，证明了圆桶的封闭性确实很好，但却发现当圆桶在海水中的速度超过12.2米时，圆桶就会因碰撞破裂，而我们要解决的问题是圆桶在到达91米的海底时的速度会不会超过了12.2米的问题。

问题假设

假设圆桶在下落到海底之前不受任何因素干扰；

假设圆桶内的东西不会在下落时发生爆炸；

假设圆桶在下落过程中受到的重力加速度不会因为深度的改变而改变；

假设圆桶在下落时是垂直下落的；

假设圆桶在下落开始时速度为零，海面距离为零；

符号说明

1.W 圆桶受到的重力；

2. F 圆桶受到海水的浮力；

3.M 圆桶的质量；

s 圆桶下落过程中与海面的距离；

t 圆桶在下落过程中于海平面距离s时的时间力；

v

t

圆桶在t时刻下落时的速度；

k 阻力与速度的比例系数；

a

t

圆桶在t时刻下落过程中的加速度；

V 圆桶的体积；

D

t

圆桶在t时刻下落时受到的阻;

g 为重力加速度；

问题分析

圆桶在下沉过程中受到3个力的作用的作用：重力W垂直向下，受到海水的浮力

F=V1025.94(N/m3)垂直向上，受到的阻力D

t

且与圆桶在水中的方位无关，与下落速度

v成正比，比例系数为k；速度v

t t

积分与t时间比就是t时刻的加速度a

t

，t时刻的距离s积

分与t时间积分的比就是速度v

t

，根据牛顿第二定理可得出距离s与v

t

模型的建立

的二介的关系。

根据上述可知：D

t

=kv

t

； （1）

下落过程时垂直下落，在t时刻是下落到与海平面距离s处，可得速度与距离的函数为

v(s)，g?9.8(ms2)；根据牛顿定律可得：

M=W

g

； （2）；

v=ds

t dt

a=dv

t dt

； （3）；

； （4）；

Ma ?W?F?D

t t

； （5）；

综上（1）（2）（3）（4）（5）的到微分方程为（6）（7）：

d2s

M ?W

dt2

F?kds

dt

； （6）；

s(0)?0,v(0)?0； （7）；

模型的求解

根据（2）（3）（4）（6）（7）进行化为一介微分方程：

v dv?Mds； （8）；

W?F?kv

两边积分：

? v dv??

W?F?kv

Mds； （9）；

求的解为：

1 s?

M

v(s)k

W?F

k2

ln(W?F?kv(s))；

W?F

g v(s)

即 s?

W k

W?F

k2

ln(W?F?kv(s)) ； （10）

W?F

S=91，即v(91)通过MATLAB求的v(91)?13.64 ms；根据（10）的函数关系式大概描述如下图：

v m

s

s/m

结果分析：图像和结果所显示，圆桶在到达海底91米时，圆桶的速度已达到了和超过12.2ms，因此他们的做法是不安全的。

模型分析

对于动态变量可通过建立微分方程，建立模型，可以求出答案，但没有用到计算