耦合振动与04年全国物理决赛题.docVIP

谈谈耦合振动 作为引子，也作为例子我们先来看看2004年全国中学生物理决赛最后一题。原题是这样叙述的： 如图，a b c d 是位于光滑水平桌面上的四个小物块，它们的质量均为m .a 、b 间有一自然长度为l,劲度系数为的弹簧联结。c、d之间有一自然长度为l,劲度系数为的弹簧联结四个物块的中心在同一直线上。如果b、c发生碰撞，碰撞是完全弹性的，且碰撞时间极短。开始时，两个弹簧都处在自然长度状态，物块c、d静止，物块a、b 以相同的速度向右运动，试定量论述 若 = ,四个物块相对于桌面怎样运动？ 若 =4 ,四个物块相对于桌面怎样运动？ 图如下： （摘自《2005全国中学生物理竞赛专辑》 北京教育出版社） 下面是原书给出的解答（有删动）： 在a b 两物块向右运动过程中，b 与c 第一次相遇并发生碰撞，在极短的时间内，两根弹簧的长度都来不及发生变化，处在原长状态，故a,d的运动未发生改变 。容易写出四个物块的运动方程： (1) ( k1=k2 时不再讨论) k1=4k2 时 （2） 第一次碰撞后 物块相对于桌面坐标系的位置由 （1） 确定 各物块速度为： （3） 如果在以后的过程中 b c 会相遇利用（1）则有： （4） 因得 sinω2t(cosω2t+1)=0 2 （5） 其解为 ω2t=nπ n为整数 （6） 从而使t异于0的最小正整数解满足 ω2t=π (7) 此时 b c 相对于质心系的速度为： vb’=-v0/2; vc’=-v0/2; (8) 这表明，b c 相遇时不相撞。 于是b c 之间无相互作用，各自独立运动… … (以后部分是按照b c 无相互作用解的，暂略) 问题是,此时b c 真的仅仅是相遇而不相撞且将保持各自独立运动吗？从速度上看似乎是这样，但是如果我们分析一下振动图像，就会得出不同的结论。注意到Tc=2π/ω2 , Tb=2π/2ω2, Tc=2Tb 在t0时刻,b c 有相同的位置和速度，但在t=t1时刻，XcXb 也就是说 c 在b 的左边！要做到这点，c 必须穿过b !!这当然是不可能的。因此，上面得出的b c 各自独立运动的结论是错的。虽然，b, c 在t=t0时刻速度相等，不相撞，但在t=t0+ε(ε0,且任意小)马上会产生相互作用。 值得一提的是，从图上看，在t=t2时刻，好像b c 发生碰撞，但实际上，此时vb0 相左运动, vc0 向右运动; 也就是说，b c 发生分离而不是相撞。当然，在t=t0时刻以后按原模型的任何分析已经没有意义。 下面，尝试给出此题目的正确解法。 如图， 我们在质心系中考虑此问题，并选取t0时刻为时间的起点 即t=0. 因此可以将b c 在有相互作用的这段时间内看成是一个整体。 这就得到了一个三自由度的耦合振动（如图）： 开始时，速度为（质心系中） Va=Vc=V0/2;Vb,c=-V0/2, (9) 各个弹簧都处在自然伸长状态，即: Xa=Xbc=Xd=0; (10) 参照下图，列出动力学方程： (11) (12) (13) 其中， Xi 表示i(i=1,2,3) 相对于平衡位置的偏移量。 记 方程组可化为： （14） （15） (16) 由于我们对运动方程感兴趣，故需解出Xi来； 记 方程可化为： （17） 其中， 欲对X 实行变换T 使在变换T 下A 为对角阵，即： ，而 为对角阵； 这样（17）可化为： （18） 这个方程是容易解的。 首先求特征值： |A-ΛI|=0;由此得方程： (19)解得三个特征值: 于是解得 利用初值,最终可解得: (19) 式中 而b c之间的作用 当时 b c 分离，故可求得当ωt=1.31时 b c 将分离。 为避免犯原解的错误，还可以检验T(ωt)在ωt=1.31处单调减，这说明 b c 确实在此刻分离了。 进一步，运用质心系中运动方程，可以求得物块 a b c d 最终运动状态。这里略去冗杂的推导直接给出我们最关心的b c的末态运动（在系统质心系中，且以b c 分离时刻为时间零点）： （20） （21） 并且可以验证，此后b c 不