一维弹簧的运动(下)
By 苏剑林 | 2014-03-13 | 28974位读者 |在上一篇文章中,我们得到了一维弹簧运动的方程
m∂2X∂t2=k∂2X∂ξ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一端被拉开,只要x2到x1的距离足够长,也就是远远超过了原长,那么释放之后,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}},
}
March 18th, 2014
超过这个区域,就有“碰撞”发生了,于是就不能视作质点了,于是就要考察弹簧本身的属性了,比如直径,杨氏模量,碰撞系数等等
我考虑过一个模型,就是它们同时还受到δ(xi−xi−1)的势能,这样子一旦碰撞,就立马反弹。想象中这样的模型应该可以成立的,估计可以给出适当的方程,但是我在数学上还不能描述它。