在上一篇文章中,我们得到了一维弹簧运动的方程
m2Xt2=k2Xξ2


并且得到了通解
X=F(u)+H(v)=F(ξ+βt)+H(ξβt)

或者
X(ξ,t)=12[X0(ξ+βt)+X0(ξβt)]+12βξ+βtξβtX1(s)ds

在文章的末尾,提到过这个解是有些问题的。现在让我们来详细分析它。

弹指神通
让我们用一个例子来演示一下上面的解,不妨设初始状态为
ξ={0,x<0sinx,x[0,π2]1,x>0


反解得
x=arcsinξ

代入通解表达式
X(ξ,t)=12[arcsin(ξ+βt)+arcsin(ξβt)]

一切看起来都没有问题,但是上述解确实有不当之处的。注意,ξ的取值是0到1,比如我们考虑ξ=1,t>0,那么ξ+βt>1,此时arcsin(ξ+βt)就不是实数了!可是好端端一个通解,怎么会出现这样的情况呢?

原因在于,当我们用对ξ的积分代替了求和时候,实际上使用了一个隐藏的假设:连续性!

也就是说,我们假设初始状态xi位于xi+1的左边,经过一段时间运动之后,xi还是位于xi+1的左边。这不一定成立,比如我们只考虑两个物体的弹簧,初始状态如下

弹簧演示

弹簧演示

其中x1一端绷紧,而x2一端被拉开,只要x2x1的距离足够长,也就是远远超过了原长,那么释放之后,x1向右运动,x2向左运动,而我们的假设之中它们都是质点,没有碰撞这一概念,而且在我们的方程之中,两个质点只是被一个与距离平方成正比的势能联系着,并不是由一根具体的弹簧拉扯着,所以可以使得x2在左边,x1在右边,并且两者在一段时间内越来越远。这样子就破坏了连续性假设。在无限多个质点的情况也是类似的。

于是,总的来说,我们得到的“通解”,实际上只适合描写0<ξ+βt<1,0<ξβt<1这个区域内的运动。超出了这个区域,也就需要修改了。至于怎么修改,我也还在思考当中。

转载到请包括本文地址:https://spaces.ac.cn/archives/2434

更详细的转载事宜请参考:《科学空间FAQ》

如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。

如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!

如果您需要引用本文,请参考:

苏剑林. (Mar. 13, 2014). 《一维弹簧的运动(下) 》[Blog post]. Retrieved from https://spaces.ac.cn/archives/2434

@online{kexuefm-2434,
        title={一维弹簧的运动(下)},
        author={苏剑林},
        year={2014},
        month={Mar},
        url={\url{https://spaces.ac.cn/archives/2434}},
}