数学建模第12次课作业-微分模型

问题给出

美国原子能委员会以往处理浓缩的放射性废料的方法， 一直是把它们装入密封的圆桶里，然后扔到水深为90多米的海底。生态学家和科学家们表示担心，怕圆桶下沉到海底时与海底碰撞而发生破裂，从而造成核污染。原子能委员会分辩说这是不可能的。为此工程师们进行了碰撞实验。发现当圆桶下沉速度超过12.2 m/s与海底相撞时，圆桶就可能发生碰裂。这样为避免圆桶碰裂，需要计算一下圆桶沉到海底时速度是多少?

已知圆桶质量为 239.46 kg，体积为 0.2058，海水密度为 1035.71kg/ ，如果圆桶速度小于12.2 m/s就说明这种方法是安全可靠的，否则就要禁止使用这种方法来处理放射性废料。假设水的阻力与速度大小成正比例，其正比例常数 k = 0.6 。现要求建立合理的数学模型，解决如下实际问题：

（1）判断这种处理废料的方法是否合理?

（2）一般情况下， v 大， k 也大； v 小， k 也小。当 v 很大时，常用 kv 来代替 k ， 那么这时速度与时间关系如何?并求出当速度不超过 12.2 m/s，圆桶的运动时间和位移应不超过多少? ( k 的值仍设为 0.6)

模型分析

木桶受力：

得出微分方程：

第（2）问中，可得微分方程：

模型求解

问题（1）

clc,clear
syms m V rho g k %符号化
s=dsolve('m*D2s-m*g+rho*g*V+k*Ds','s(0)=0,Ds(0)=0');%微分方程（1）
s=subs(s,{m,V,rho,g,k},{239.46,0.2058,1035.71,9.8,0.6});%代入数值
s=vpa(s,10) %求位移函数
v=dsolve('m*Dv-m*g+rho*g*V+k*v','v(0)=0');%微分方程（2）
v=subs(v,{m,V,rho,g,k},{239.46,0.2058,1035.71,9.8,0.6});
v=vpa(v,7) %求速度函数
y=s-90;
tn=solve(y) %求到达海底 90 米处的时间
vn=subs(v,tn) %求到底海底 90 米处的速度

求解后可知速度为13.772m/s，超过了安全速度，因此不合理。

问题（2）

clc,clear
syms m V rho g k t real
v=dsolve('m*Dv-m*g+rho*g*V+k*v^2','v(0)=0');%微分方程（3）
v=subs(v,{m,V,rho,g,k},{239.46,0.2058,1035.71,9.8,0.6});%代入数值
v=vpa(v,7) %求速度函数
v1=real(v),v2=imag(v) %求速度函数的实部和虚部
tnn=solve(v1-12.2); %求时间的临界值
s=int(v1,0,tnn) %求位移的临界值
tnn=double(tnn)

时间不超过13.00s,
位移不超过84.84m.