利用Matlab进行弹簧振子运动状态的仿真(打印版).doc

利用Matlab进行弹簧振子运动状态的仿真 汪家培（安庆师范学院物理与电气工程学院 安徽 安庆 246001） 指导教师：张杰 摘要:　简谐振动是最基本的振动，弹簧振子是最典型的简谐振动。本文从牛顿第二定律出发推导了弹簧振子的运动学方程，从而进一步给出由两相同轻质弹簧和一个质量为m的小球组成的非线性弹簧振子系统的运动学方程，利用Matlab建立了该系统的仿真模型，研究了在不同条件下系统的的运动图像、能量图和相图，详细分析了此非线性系统的特点。 关键词：弹簧振子，能量图，相图，运动图，计算机仿真 1. 引言 普通物理力学中[1-2]讲振动时，常以水平弹簧振子和单摆为例作为问题的切入点，它们只是非线性振动的近似(胡克定律有一定的适用范围)。所谓非线性振动是指恢复力与位移不成线性比例或阻尼力与速度不成线性比例的系统的振动。一般说，线性振动只适用于小运动范围，超过此范围，就变成非线性振动。非线性系统的运动微分方程是非线性的，不能用叠加原理求解。方程中不显含时间的非线性系统称为非线性自治系统；显含时间的称为非线性非自治系统。保守非线性自治系统的自由振动仍是周期性的，但其周期依赖于振幅。对于渐硬弹簧，振幅越大，周期越短；对于渐软弹簧，振幅越大，周期越长。非保守非线性自治系统具有非线性阻尼，阻尼系数随运动而变化，因而有可能在某个中间振幅下等效阻尼为零，从而能把外界非振动性能量转变为振动激励而建立起稳定的自激振动（简称自振）。弦乐器和钟表是常见的自振系统。周期地改变系统的某个参量而激起系统的大幅振动称参变激发2. 弹簧振子系统概述 弹簧振子是一种简谐振子。它的典型结构是由一个一端固定，质量可以忽略的轻弹簧和连在它另一端（自由端）与一个跨在气垫导轨上的滑块相连，如图（1）。这样的，以及与此类似的系统称为弹簧振子。将滑块自弹簧自由伸展时的O点移动一小位移至A点，然后释放。可观察到滑块作变加速运动，由A点通过O 点到达B点，继而又经过O点回到A点。并且再度开始与上述过程完全相同的运动，如此往返。不考虑作用于滑块的空气阻力及气垫的摩擦。将滑块视作质点，弹簧自由伸展时质点的位置是平衡位置，以此为坐标原点建立坐标系Ox。X表示质点的位置坐标，又等于相对与原点的位移，也是弹簧的伸长（压缩）量。x很小是，力f与x之间成线性关系，即f=-kx， k是弹簧劲度系数，弹簧弹性力是线性回复力，弹簧振子做简谐振动。 根据牛顿第二定律有[7] （1） 用m除上式两端，并令，上式可写作 或 （2） 上式即为弹簧振子的运动方程。根据常微分方程的理论，（2）式的解可以写作 （3） A和是待定常数，需要根据初始条件来确定。（3）式就是简谐振动的运动学方程。 由于 令 代入上式后再代入（3）式，得 可见，简谐振动的运动规律也可以用正弦函数表示。正弦和余弦函数都是周期函数，因此简谐振动是围绕平衡位置的周期运动。 按简谐振动的运动学方程，物体的最大位移不能超过A，物体离开平衡位置的最大位移的绝对值叫做振幅。为了确定振幅A的大小，现对位移对时间求一阶导数，得： （4） 将初始条件t=0，，代入上式后取二式平方之和，既可求出振幅， （3）式表明，只有知道A、ω0、α，才能完全决定系统的运动状态。在（3）式中我们把时间T的线性函数φ=ω0t+α叫做简谐振动的相位。由于简谐振动的位移、速度是按余弦、正弦规律变化的，所以相位是当振幅一定时进一步决定简谐振动任何瞬时运动状态的物理量。所以我们可以用比较两个系统的相位方法来比较两个简谐振动的运动状态。 弹簧振子中线性回复力为弹簧的弹性力，是保守力。所以简谐振动系统的总机械能守恒。关于弹簧振子，应用质点动能公式 ，，将（3）式代入可得到： 因为。所以有 （5） 至于势能，，将简谐振动运动学方程代入，可得 （6） 有(5)、(6)两式可知：弹簧振子的动能和势能按余弦或正弦的平方随时间变化。显然，动能最大时势能最小，动能最小是势能最大，简谐振动的过程正是动能与势能相互转换的过程。由于弹簧振子的运动行为已经有很多文章讨论过，我们这里就不再赘述，下面我们将讨论偶合弹簧振子的运动