20
Jan
《方程与宇宙》：三体问题和它的初积分(六)
By 苏剑林 | 2011-01-20 | 61892位读者 | 引用The Three Body Problem and its Classical Integration
很多天文爱好者都已经接触到了“二体问题”（我们在高中学习到的“开普勒三定律”就是内容之一），由于在太阳系中行星质量相对较小而且距离相对较远，应用“二体问题”的解对天体进行计算、预报等能够满足一定的近似需求。不过，如果需要更高精度的计算，就不能把其他行星的引力给忽略掉了，于是就产生了所谓N体问题（N-Body Problem），即N个质点尽在它们各自引力的相互作用下的运动规律问题。最简单的二体已经被彻底解决，而三体或更多体的问题则与二体大相径庭，因为庞加莱证明了，三体问题不能严格求解，而且这是一个混沌系统，任何微小的扰动都会造成不可预期的效果。
根据牛顿力学，选择惯性参考系，设三个质点分别为$M_1,M_2,M_3$，向径分别为$\vec{r_1},\vec{r_2},\vec{r_3}$，可以列出运动方程（以下的导数都默认是对时间t求导）
26
Jan
唠叨下，关于三体问题周期轨道
By 苏剑林 | 2011-01-26 | 32547位读者 | 引用
4
Feb
[春礼]《方程与宇宙》:圆形限制性三体问题(七)
By 苏剑林 | 2011-02-04 | 23725位读者 | 引用平面圆形限制性三体问题运动方程及能量积分
plane circular restricted three-body problem
02.04有重要修正！！
寒假一个很大的目标就是能够在三体问题的周期轨道上有点突破，于是就出动了“向量”、“复分析”、“微分方程”等理论“核武”，遗憾的是，“有心栽花花不开”，到今天还是没有多少进展。不过俗语也说“无心插柳柳成荫”，也不错。今天回看《天体力学引论》中的“圆形限制性三体问题”，经过一番思考，利用这些天的思考方法重新推导出了其运动方程和能量积分，也算是“意外收获”在此作为春节礼物与大家分享。
平面圆形限制性三体问题
所谓“圆形限制性三体问题”，就是指两个大质量天体（质点）在它们相互引力作用下做圆周运动，假设第三天体（质量趋于0）只受到这两个天体的引力作用而不影响两个天体运行的一种运动情况。由于普通三体问题无法积分，而这个“限制性模型”能够把问题化简不少（不过还是不能积分出来的），因此也得到了一定应用。它的应用条件是：第三体质量小（如当前航天器与地球、太阳）、短程。注意短程也是相当重要的条件之一，注意短程也是相当重要的条件之一，质量越小应用范围越大。要是质量大的话，就不能计算太长的路程。
19
Feb
《方程与宇宙》:一种有趣的三体问题坐标
By 苏剑林 | 2011-02-19 | 22875位读者 | 引用通常来说，选取惯性系为参考系，列出的三体问题方程为
$$\ddot{\vec{r}}_k=\sum_{i=1,i != k}^{n} Gm_i\frac{\vec{r}_i-\vec{r}_k}{|\vec{r}_i-\vec{r}_k|^3}$$
历史上出现过很多不同形式的变换，使得三体问题的运动方程有了各样的形式，如Lagrange形式、Jacobi形式、Hamilton形式等。这些变换形式都各有特点，都能够在一定程度上化简三体问题。BoJone在研究摆弄等质量型三体问题的运动方程时，也发现了一种很有趣的变换，在此贴出与大家分享。
设$\vec{R}_1=\vec{r}_1-\vec{r}_2,\vec{R}_2=\vec{r}_2-\vec{r}_3,\vec{R}_3=\vec{r}_3-\vec{r}_1$，则三体问题的运动方程变为
8
Mar
沐浴问题——调控水温
By 苏剑林 | 2011-03-08 | 25151位读者 | 引用载入正题之前，不妨闲扯一下BoJone的家...
BoJone在一些文章中已经提到过，我是一个来自农村的孩子，目前我的家也在农村。虽然生活并不能说“贫困”，家中也添置了不少电器，不过一直没有购置的就是洗衣机和热水器。洗衣机嘛，我觉得衣服自己动手洗是很好的，至少不让自己偷懒。至于热水器，因为家在农村，所以能够比较方便地弄到一些柴草，而且稻谷收割完后的桔梗也可以当燃料用，平时烧菜一般都用烧柴草，因此热水器实在没有多大必要。（很遗憾，沼气池没有能够在这里普及起来，大家可不要责怪我排放温室气体哦...^_^）
烧柴草的炉灶
既然没有热水器，那只能人工烧水了。往往是烧好一大锅水，洗澡时盛一盆子，然后加水降温，接着就可以洗白白了。本文的问题正是来源于调水温。当水很热时，为了加快降温，我们往往“双管齐下”：一边向盆子注入冷水，一般从盆子放出热水。于是就有了一个问题：水的温度与时间成什么关系？
5
Apr
重提“旋转弹簧伸长”问题（变分解法）
By 苏剑林 | 2011-04-05 | 19993位读者 | 引用感谢Awank-Newton读者的来信，本文于2013.01.30作了修正，主要是弹性势能的正负号问题。之前连续犯了两个错误，导致得出了正确答案。现在已经修正。参考《平衡态公理的修正与思考》
在下面的两篇文章中，BoJone已经介绍了这个“旋转弹簧伸长”的问题，并从两个角度提供了两种解答方法。前者列出了一道积分方程，然后再转变为微分方程来解；后者直接从弹性力学的角度来列出一道二阶微分方程，两者殊途同归。
http://kexue.fm/archives/782/
今天，再经过一段时间的变分法涉猎后，BoJone尝试从变分的角度（总能量最小）来给出一种新的解法。同样设r为旋转达到平衡后弹簧上一点到旋转中心的距离，该点的线密度为$\lambda =\lambda (r)$，该点到中心的弹簧质量为$m=m(r)$，旋转前的长度为$l_0$，旋转平衡后的长度为$l_1$。由于弹簧旋转后已经达到了平衡状态，由平衡态公理（参看《自然极值》系列），平衡意味着总能量“动能-势能”取极值。
我们在研究地球附近的小天体运动时，如果把天体和地球看作一个二体系统，那最多只能算上一个零级近似，如果使用“地球+月球+小天体”组成的圆形限制性三体问题模型，那可以算上一个二级近似了。那么，一级近似又是什么了。BoJone认为，它就是本文将要讲的“双固定引力中心问题”了，也叫“双不动中心问题”，英文名是two fixed-center problem。这是一种特殊的限制性三体问题。在这个三体系统中，两个主天体（或称有限质量天体）固定不动，第三个小天体在两个固定的主天体吸引下运动。欧拉、拉格朗日、勒让德、雅可比等人很早就研究过这个问题。其中，欧拉最先成功地求出了这个系统的积分。[引用]
另外，双固定引力中心问题还有另外一个应用，在研究人造卫星的运动时，可以只考虑地球引力，但是由于地球不是完美的球体，把其看成一个质点其实不十分精确，要是把它拆分为两个引力源，就可以很大程度上提高精确度。毕竟双固定引力中心问题是完全可以积分的，可以作为一个比较好的中间轨道（介乎圆锥曲线和精确轨道之间的）。
看完了“双不动中心”问题，我们不妨再来看一个貌似简单一点的力学问题，在一个固定质点的引力吸引的基础上，增加一个恒力作用，研究这样的力场中小天体的运动。
咋看上去这个问题比“双不动中心”简单多了，至少运动方程也显得更简单：
$$\ddot{vec{r}}=-GM\frac{\vec{r}}{|\vec{r}|^3}+\vec{F}$$
其中$\vec{F}$是一个常向量。不过让人比较意外的是，这个问题本质上和“双不动中心”是一样的，它可以看作是双不动中心问题的一个极限情况。而且它们的解法也是惊人地相似，下面我们就来分析这一个过程。
首先很容易写出这个方程的能量守恒积分：
$$1/2 \dot{vec{r}}^2-GM\frac{1}{|\vec{r}|}-\vec{F}\cdot \vec{r}=h$$
最近评论