ケプラー軌道シミュレーター 戻る
宇宙力学・天体物理

ケプラー軌道シミュレーター — 惑星の運動と万有引力

離心率・長半径を操作して楕円軌道をリアルタイムにアニメーション。ケプラーの等面積法則(第2法則)と T²=a³(第3法則)を視覚的に確認しよう。

軌道パラメータ
長半径 a 1.00 AU
離心率 e 0.017
シミュレーション速度
表示設定
プリセット
計算結果
公転周期 T
T²/a³(=1.0で正確)
現在の距離 r
現在の速度 v
近日点距離 r_peri
遠日点距離 r_aph

ケプラーの法則

第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$

数値積分との関係: Verlet法は構造動解析の時間積分(Newmark-β法等)と同系統の手法です。宇宙機CAEでも軌道解析が基礎になります。
公転周期 T (年)
速度 (km/s)
距離 (AU)
T²/a³

黄円 = 太陽(焦点) 水色三角形 = 等面積(第2法則) 矢印 = 速度ベクトル

ケプラー軌道とは

🧑‍🎓
「惑星の軌道は楕円」って教科書で見ますけど、どうやってシミュレーターで確かめられるんですか?
🎓
ざっくり言うと、太陽を中心にした「長半径」と「離心率」で軌道の形が決まるんだ。このシミュレーターの上のスライダーで「離心率」を0から0.9くらいに動かしてみて。e=0なら真円、eが大きくなるほど細長い楕円になるのが一目瞭然だよ。
🧑‍🎓
え、そうなんですか!離心率を大きくすると、惑星の動く速さもガラッと変わるんですね。太陽に近いところではめっちゃ速くて、遠いところではすごくゆっくり。
🎓
その通り!これがケプラーの第2法則「面積速度一定」の直感的な表現だ。シミュレーターの「等面積表示」をオンにすると、等時間ごとに描かれる扇形の面積が常に同じなのがわかる。実務では、探査機が近日点でスイングバイする際の速度計算にこの原理が使われてるんだ。
🧑‍🎓
プリセットで「ハレー彗星」ってありますね。これ、離心率が0.967って…めちゃくちゃ細長い軌道なんですか?
🎓
そうだね。ハレー彗星を選んでアニメーションを再生すると、太陽の近くを猛スピードで通り過ぎて、遠日点ではほとんど動いてないように見えるはず。このシミュレーターでは「長半径」も変えられるから、例えば地球(a=1, e≈0)と火星(a≈1.5)を比べると、第3法則「周期の2乗は長半径の3乗に比例」が成り立って、火星の方がゆっくり公転するのが確認できるよ。

物理モデルと主要な数式

惑星の位置は、太陽を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つの天体が互いに引力を及ぼし合う系は、単純な楕円軌道では記述できず、カオス的な振る舞いを示す。ここまで来れば、最新の宇宙探査機の軌道設計や、系外惑星系の安定性の研究の最前線に触れることになる。まずは、このシミュレーターで「二体問題」の完全な理解を固めることが、すべての基礎になるんだ。