ケプラーの法則
第1法則: 楕円軌道 $r = \dfrac{a(1-e^2)}{1+e\cos\theta}$
第2法則: 面積速度一定 $\dfrac{dA}{dt}= \text{const}$
第3法則: $T^2 = a^3$(AU/yr単位)
$GM_\odot = 4\pi^2\ \text{AU}^3/\text{yr}^2$
黄円 = 太陽(焦点) 水色三角形 = 等面積(第2法則) 矢印 = 速度ベクトル
離心率・長半径を操作して楕円軌道をリアルタイムにアニメーション。ケプラーの等面積法則(第2法則)と T²=a³(第3法則)を視覚的に確認しよう。
第1法則: 楕円軌道 $r = \dfrac{a(1-e^2)}{1+e\cos\theta}$
第2法則: 面積速度一定 $\dfrac{dA}{dt}= \text{const}$
第3法則: $T^2 = a^3$(AU/yr単位)
$GM_\odot = 4\pi^2\ \text{AU}^3/\text{yr}^2$
黄円 = 太陽(焦点) 水色三角形 = 等面積(第2法則) 矢印 = 速度ベクトル
惑星の位置は、太陽を1つの焦点とする楕円軌道の極座標形式で表されます。これがケプラーの第1法則です。
$$r(\theta) = \frac{a(1-e^2)}{1 + e \cos \theta}$$ここで、$r$は太陽からの距離、$\theta$は近日点からの角度(真近点角)、$a$は軌道長半径、$e$は離心率です。$e=0$で円、$0 < e < 1$で楕円を表します。
惑星の運動は、面積速度が一定(ケプラーの第2法則)という条件と、万有引力の法則から導かれます。シミュレーターでは、この運動を数値的に解いています。
$$\frac{dA}{dt}= \frac{1}{2}r^2\dot{\theta}= \text{const.}$$$\frac{dA}{dt}$は面積速度、$\dot{\theta}$は角速度です。この式から、太陽に近い($r$が小さい)ときは角速度$\dot{\theta}$が大きく、遠いときは小さくなることがわかります。
宇宙探査機の軌道設計:火星や小惑星への探査機は、地球の重力を利用したスイングバイを行います。この軌道計算の基礎がケプラー軌道力学です。近日点での高速移動を利用して加速し、燃料を節約します。
人工衛星の配置:通信衛星や気象衛星は、多くの場合、地球を焦点とする楕円軌道または円軌道に投入されます。特に地球観測衛星は、特定の地域を定期的に観測できるよう、軌道長半径と離心率が精密に設計されます。
彗星・小惑星の運動予測:ハレー彗星のような長周期彗星や、地球に接近する可能性のある小惑星(NEO)の将来の位置を予測するために、その軌道要素(長半径、離心率など)を測定し、ケプラーの法則に基づいて計算します。
CAEにおける軌道解析:自動車の衝突試験や構造物の動的解析で用いる数値積分法(Verlet法、Newmark-β法等)は、この惑星運動のシミュレーションで使われる手法と同系統です。宇宙機の構造CAEでは、軌道力学と連成させた解析が行われます。
まず、「シミュレーターの軌道は永遠に安定」という思い込みに注意だ。現実の太陽系では、他の惑星の重力(摂動)が働くため、軌道は完全な楕円からわずかにずれていく。このツールは「二体問題」だけを扱っているから、地球と太陽だけの理想的な世界だと思って使おう。例えば、木星の巨大な重力は小惑星帯の軌道に大きな影響を与えているんだ。
次に、パラメータ設定の落とし穴。離心率eを1以上に設定すると、軌道は楕円ではなく放物線や双曲線になってしまい、惑星は太陽に戻ってこない。あくまで惑星の「閉じた軌道」を学ぶツールだと理解しておこう。また、長半径aを大きくしすぎると、アニメーションで惑星が画面から一瞬で消えちゃうことがある。初期値の1(天文単位)を基準に、まずは火星の1.5倍くらいまでで遊ぶのがコツだよ。
最後に、「等面積表示」の本当の意味。扇形の面積が同じなのは、時間間隔が一定だからだ。でも、この「面積速度一定」は、太陽からの引力(中心力)だけが働いている時にだけ成り立つ法則だ。実務で衛星を設計する時、大気抵抗や太陽光圧が加わると、この法則は成り立たなくなる。シミュレーターの理想環境と現実の違いを意識することが、次のステップへの第一歩だ。
このシミュレーターの根幹をなす計算手法は、「数値積分」と「常微分方程式の初期値問題」そのものだ。惑星の位置と速度を少しずつ更新して軌道を描く方法は、自動車の衝突シミュレーション(クラッシュCAE)で車体の変形を追う「動的解析」と数学的に同じ土台にある。例えば、衝突時のダミー人形の挙動も、ニュートンの運動方程式を短い時間ステップで解き続けているんだ。
また、制御工学とも深く関わる。静止衛星を所定の位置に保つ「軌道制御」や、探査機が小惑星にランデブーするための軌道修正は、まさにこのケプラー軌道を出発点として設計される。軌道力学の予測と、スラスタ噴射によるフィードバック制御が組み合わさることで初めて実現する技術だ。
さらに意外なところでは、半導体設計の分野でも類似の計算が登場する。真空中で帯電した粒子(イオン)を電磁場で操り、ウェハ上に正確に打ち込む「イオン注入」プロセスでは、粒子の軌道が放物線や双曲線になることがあり、その計算基礎は万有引力の場と非常に似通っている。異なる分野でも、同じ物理法則が顔を出す好例と言えるね。
このツールに慣れたら、次は「なぜ」楕円軌道になるのかを数学的に追ってみよう。キーワードは「エネルギー保存則」と「角運動量保存則」だ。惑星の運動エネルギーと位置エネルギー(重力ポテンシャル)の和が一定であること、そして太陽の周りの角運動量が保存されることから、先ほどの極座標の式 $r(\theta)$ を導出できる。大学初年度の物理の教科書の「中心力場での運動」の章が、次の最適な教材だ。
数値計算にもっと踏み込みたいなら、このシミュレーターが背後で使っている「オイラー法」や「ルンゲ=クッタ法」を自分で実装してみることを強く勧める。例えば、ExcelやPythonで簡単なプログラムを書いて、万有引力の式 $F = GMm/r^2$ から直接、位置と速度を更新してみるんだ。時間ステップdtを粗くすると軌道が螺旋を描いて太陽に落ちてしまう(エネルギーが保存されない)など、数値解法の難しさと面白さが一気に体感できる。
最後のステップは、「三体問題」への挑戦だ。月・地球・太陽のように、3つの天体が互いに引力を及ぼし合う系は、単純な楕円軌道では記述できず、カオス的な振る舞いを示す。ここまで来れば、最新の宇宙探査機の軌道設計や、系外惑星系の安定性の研究の最前線に触れることになる。まずは、このシミュレーターで「二体問題」の完全な理解を固めることが、すべての基礎になるんだ。