有质动力:倒立单摆的稳定性
By 苏剑林 | 2013-12-29 | 49096位读者 |前几天在“宇宙的心弦”浏览网页时,发现他更新了一篇很有趣的文章,叫《倒立单摆的稳定性与Ponderomotive Force》(果然,物理系的能接触到各种各样有趣的现象),里边谈到通过施加一个运动在单摆上面,倒立的单摆也可以是稳定的。这勾起了我的兴趣,遂也计算了一番。
单摆有两个平衡点,一个是垂直向下的稳定平衡点,另一个是垂直向上的不稳定平衡点。用一个形象的说法,那就是如果你有一个垂直向下的单摆,那么你稍微晃一下它,它只是来回摆动;如果那个单摆是垂直向上的,那么如果你晃一下它,它就会迅速倒下到垂直乡下的位置了。如果给单摆施加一个运动,可以让倒立的单摆变成稳定的?咋看之下,很难判断这种运动是否存在,但是只需要想象将拉着单摆向下以$10m/s^2$的加速度向下运动,那么倒立的单摆自然就稳定了,因为这是的力已经反方向了。这其实我们,非惯性运动可能给我们带来一个额外的力场,使得倒立单摆稳定。当然,上面举的加速运动的例子没有太大用处,因为它的运动幅度太大了,以至于那已经不是我们能控制的单摆了。但是还有另外一种运动有可能使得倒立单摆稳定,那就是高频振荡。我们将详细分析它。
Physixfan的文章也对这个问题进行了一番数学分析。但是我总觉得他那样的分析过程是不能让人满意的,他往原系统中的运动方程$m\ddot{x} = -dU/dx$加了一个“高频的振荡项”$f$,变为
$$m\ddot{x} = -dU/dx+f$$
但是这里的$f$意义是不明确的,如果$x$表示位移,那么$f$可以理解为力,但是$x$也可以是角度或者其他变量,那时就不好把握了。另外,即使将$f$解释为力,那么“将一个高频振荡的力加在单摆上面”与实验而言其实也是不清楚的,更清晰的做法应该是将一个高频振荡直接加在原系统上面,换句话说,把单摆放到一个高频振荡的参考系里边去。(想象一个极端的情况,有一个以极快速度不断一上一下的电梯,单摆就放在里边。)
在这里采用最小作用量原理来分析是最为方便的,以垂直向上为x轴正方向,水平向右为y轴正方向,那么单摆的作用量为:
$$S=\int \left[ \frac{1}{2}m(\dot{x}^2+\dot{y}^2)- mgx \right]dt$$
其中还有$x^2+y^2=l^2$这一约束,由此可以导出单摆的运动方程。选取这样的坐标系对于分析倒立单摆较为方便。
现在,将一个$x=h(t)$的振荡项加到原单摆中,也就是说让单摆的固定中心以$x=h(t)$的方式运动这,这样一来,单摆的作用量还是
$$S=\int \left[ \frac{1}{2}m(\dot{x}^2+\dot{y}^2)- mgx \right]dt$$
只不过约束条件变为
$$[x-h(t)]^2+y^2=l^2$$
现在就设$x=l\cos\theta+h(t),y=l\sin\theta$,代入作用量,得到
$$S=\int \left[ \frac{1}{2}m(l^2\dot{\theta}^2+\dot{h}^2-2l\dot{h}\dot{\theta}\sin\theta)- mg\left(l\cos\theta+h\right) \right]dt$$
由欧拉-拉格朗日方程就得到运动方程
$$\frac{d}{dt}\left(ml^2\dot{\theta}-ml\dot{h}\sin\theta\right)=mgl\sin\theta-ml\dot{h}\dot{\theta}\cos\theta$$
即
$$l\ddot{\theta}-(\ddot{h}+g)\sin\theta=0$$
可见,稳定点$\theta=0$依然存在,当$h$项不存在时,方程$l\ddot{\theta}-g\sin\theta=0$在原点处是不稳定的,为了认识到这一点,只需用近似$\sin\theta\approx\theta$,得到$l\ddot{\theta}-g\theta=0$,解得
$$\theta=Ae^{t\sqrt{g/l}}$$
这是发散至无穷的指数解。
只要$(\ddot{h}+g)<0$,就可以让上述方程变为受约束的方程,从而让$\theta=0$处是稳定的。但是我们希望原来的单摆改变的“幅度”不要太大,因此$h$只能够是小幅振荡。不妨假设
$$h(t)=h_0 \cos(\omega t)$$
有
$$\ddot{h}=-h_0 \omega^2 \cos(\omega t)$$
代入上述方程得
$$l\ddot{\theta}+\left[h_0 \omega^2 \cos(\omega t)-g\right]\sin\theta=0$$
在$t=0,\theta=0$附近展开,近似有
$$l\ddot{\theta}+\left(h_0 \omega^2 -g\right)\theta=0$$
可见,为了让$\theta=0$成为稳定平衡点,必须有$h_0 \omega^2 -g \gg 0$,即
$$\omega \gg \sqrt{\frac{g}{h_0}}$$
如果$h_0=0.01m$,那么估计$\omega \gg 30 s^{-1}$,这只不过是$5Hz$左右,而我们交流电流的频率有$50Hz$,所以上述视频的实验是可能实现的。总共的来说,长期观察而言,就是高频振荡可以产生一项力来使得系统保持它的稳定性,这个力就叫做Ponderomotive Force,中文翻译是有质动力。这种现象在粒子物理中有着重要应用。
转载到请包括本文地址:https://www.spaces.ac.cn/archives/2232
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Dec. 29, 2013). 《有质动力:倒立单摆的稳定性 》[Blog post]. Retrieved from https://www.spaces.ac.cn/archives/2232
@online{kexuefm-2232,
title={有质动力:倒立单摆的稳定性},
author={苏剑林},
year={2013},
month={Dec},
url={\url{https://www.spaces.ac.cn/archives/2232}},
}
December 31st, 2013
写的很不错,我当时也觉得范翔的写的不明白,那个x和f究竟什么意思我也没弄懂。
话说你现在写作能力了得啊
谢谢^_^大概是读国外的著作比较多,喜欢模范他们的风格来写。当然,我觉得关键是把自己的想法交代清楚。
其实还有一些没有想清楚,过段时间我会接着写这个问题。
January 20th, 2014
这是信号系统,或控制理论的东西。
也有这个应用。据说搞核聚变的人有想法想用这个原理来操控里边的等离子~