关于阻尼振动

  任何振动都不可避免地会受到阻力的作用，这时的振动叫做阻尼振动。一般而言，如果物体所受阻力$F_f$与其速度$v$成正比，也即有：
$$F_f=-\gamma v=-\gamma \frac{\mathrm{d}x}{\mathrm{d}t}\text{\hspace{1em}(1)}$$

  $\gamma$是阻力系数。那么，在弹力、阻力的共同作用下，弹簧振子的运动方程就是：
$$m\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{t}^2}=-kx-\gamma \frac{\mathrm{d}x}{\mathrm{d}t}\text{\hspace{1em}(2)}$$

  记${\omega }_0=\sqrt{\frac{k}{m}},\beta =\frac{\gamma }{2m}$.则式$(2)$可化为：
$$\frac{{\mathrm{d}}^{2}x}{\mathrm{d}{t}^2}+2\beta \frac{\mathrm{d}x}{\mathrm{d}t}+{\omega }_0^{2}x=0\text{\hspace{1em}(3)}$$

  这个微分方程的特征方程是
$${\lambda }^2+2\beta \lambda +{\omega }_0^2=0\text{\hspace{1em}(4)}$$

  是一个一元二次方程。根据判别式$\Delta =4({\beta }^2-{\omega }_0^2)$的不同，微分方程$(3)$有三种形式的解。

弱阻尼$(\beta <{\omega }_0)$

  这种情况下，$\Delta <0$，记$\omega =\sqrt{{\omega }_0^2-{\beta }^2}$，式$(4)$有共轭复根
$${\lambda }_{1,2}=-\beta \pm \omega \mathrm{i}\text{\hspace{1em}(5)}$$

  从而得到$(3)$的解为$x={\mathrm{e}}^{-\beta t}(C_1\cos \omega t+C_2\sin \omega t)$.这一形式也可以写成式$(6)$所示的形式：
$$x=A_0{\mathrm{e}}^{-\beta t}\cos (\sqrt{{\omega }_0^2-{\beta }^2}t+{\phi }_0)\text{\hspace{1em}(6)}$$

  取$x_0=1\mathrm{m},v_0=-1m/s,{\omega }_0=10rad/s,\beta =0.5{\mathrm{s}}^{-1}.$画出的$x-t$图象如下所示：

  使用的画图方式为 GeoGebra .后面的图也是用它画的

  蓝色的包络线为$x=\pm A_0{\mathrm{e}}^{-\beta t}$，红色的是弹簧振子$x-t$图象。

过阻尼$(\beta >{\omega }_0)$

  在这种情况下，$\Delta >0$，记${\omega }^{\prime }=\sqrt{{\beta }^2-{\omega }_0^2}$，式$(4)$有两根
$${\lambda }_{1,2}=-\beta \pm {\omega }^{\prime }\text{\hspace{1em}(7)}$$

  此时运动方程$(3)$的解的形式为
$$x={\mathrm{e}}^{-\beta t}(A_1{\mathrm{e}}^{{\omega }^{\prime }t}+A_2{\mathrm{e}}^{-{\omega }^{\prime }t})\text{\hspace{1em}(8)}$$

  取$x_0=1\mathrm{m},v_0=-1m/s,{\omega }_0=10rad/s,\beta =12.5{\mathrm{s}}^{-1}.$画出的$x-t$图象如下所示：

  $x$和$v$都迅速衰减到$0$.

临界阻尼$(\beta ={\omega }_0)$

  在这种情况下，$\Delta =0$，此时$(4)$有重根
$${\lambda }_{1,2}=-\beta \text{\hspace{1em}(9)}$$

  也可以由此得到运动方程$(3)$的解为
$$x={\mathrm{e}}^{-\beta t}(A_0+A_{1}t)\text{\hspace{1em}(10)}$$

  取$x_0=1\mathrm{m},v_0=-1m/s,\beta ={\omega }_0=10rad/s.$画出的$x-t$图象如下所示：

  $x$和$v$也是很快地衰减到$0$.

总结

  把上面三张图画在一起可以得到下图：

  其中，蓝色的是包络线$x=\pm A_0{\mathrm{e}}^{-\beta t}$，红色的是欠阻尼情况下的$x-t$图象，橙色的是过阻尼情况下的$x-t$图象，紫色的是临界阻尼情况下的$x-t$图象。在系统只改变$\beta$的情况下，发现临界阻尼情况下弹簧振子最快回到平衡状态。
  实际上，可以统一$(3)$的解的形式为$x={\mathrm{e}}^{-\beta t}f(t)$.临界阻尼情况下的$f(t)$是多项式，过阻尼情况下的$f(t)\sim {\mathrm{e}}^{|{\omega }^{\prime }|t}$ 是指数阶，所以临界阻尼衰减得比过阻尼快。