MATLAB计算方法和技巧小球沉降的规律.doc

小球在水中沉降的规律 问题小球的质量为m，水的浮力为B，水的粘性力为R = -Kv，式中K是与水的粘性和小球的半径有关的常量，计算小球在水中竖直沉降随时间变化的速度和深度。 [] 如B6.1图所示，小球受到重力mg，方向向下；浮力B，方向向上；粘性力Kv，方向向上。取向下为正方向，根据牛顿第二定律可列小球的运动方程 mg - B - Kv = ma。 (6.1.1) 显然重力要大于浮力mg B。当速率为零时可得最大加速度aM = g - B/m。当加速度为零时可得极限速率vT = (mg - B)/K。 由于a = dv/dt，可得微分方程 ， (6.1.2) 分离变量得 ， 积分 ， 可得 ， 化得 ， 速度为 。 (6.1.3) 当t→∞时v→vT，这就是极限速度。 当t = m/K时，小球速率为v = vT(1 - e-1) = 0.6321vT。m/K称为特征时间，用τ表示；0.6321vT是特征速度。在t τ时，小球将以极限速度vT匀速下降。 在求小球沉降的深度时，则利用关系v = dx/dt，可得微分方程 ， 积分得 ， 即。 (6.1.) 当时间t足够长时，小球下沉的深度随时间直线增加 x = vT(t - τ)。 (6.1.) 这是一条直线，也是深度曲线的渐近线。 令t* = t/τ，t*就是无量纲的时间，也称为约化时间；令v* = v/vT，v* 是无量纲的速度，称为约化速度 v* = 1 – exp(-t*)。 (6.1.) 极限速度与特征时间的乘积是特征长度vTτ，取无量纲深度或约化深度x* = x/vTτ，则得 。 (6.1.) 牛顿运动方程是关于时间的二阶微分方程，当微分方程有解析解时，可深入讨论物体的运动规律，也画出运动曲线。 [程序].m部分如下。 clear %清除变量 tm=5; %最大时间 dt=0.1; %时间间隔 t=0:dt:tm; %时间向量 v=1-exp(-t); %求速度 x=t+exp(-t)-1; %深度 figure %创建图形窗口 subplot(2,1,1) %选子图 plot(t,v) %画速度曲线 grid on %加网格 fs=12; %字体大小 title(小球沉降的速度曲线,FontSize,fs)%标题 xlabel(时间\itt/\tau,FontSize,fs) %横坐标 ylabel(速度\itv/v\rm_T,FontSize,fs)%纵坐标 subplot(2,1,2) %选子图 plot(t,x) %画深度曲线 grid on %加网格 title(小球沉降的深度曲线,FontSize,fs)%标题 xlabel(时间\itt/\tau,FontSize,fs) %横坐标 ylabel(深度\itx/v\rm_T\it\tau,FontSize,fs)%纵坐标 [图示](1)如P6之上图所示，速度随时间减速增加，趋于极限速度。 (2)如P6之下图所示，深度开始时随时间加速增加，最后随时间直线增加。 方法二：用一个一阶微分方程的数值解。由(6.1.2)式得 ， 此式可化为 ， 其中，τ = m/K。利用约化时间t* = t/τ和约化速度v* = v/vT，可得 。 (6.1.) 这是一个关于约化速度的一阶微分方程，左边是导数，右边是导数的表达式。微分方程的初始条件为：当t* = 0时v* = 0。微分方程的数值解就是速度与时间的关系，P6.1图 。