0
0

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?

二重振り子でカオスを体感!— 初期値敏感性をシミュレーターで実証する

0
Posted at

▶ 今回のシミュレーター: 二重振り子(カオス)シミュレーター (ブラウザで動作・登録不要)

突然ですが、「予測不可能」って考えたことありますか?

天気予報が外れたり、将棋のAIが人間を驚かせる一手を指したり。世の中、「ちょっとした違い」で結果が大きく変わること、ありますよね。

実は、たった2本の棒と2つのおもりをヒモでつなぐだけの単純な装置が、この「予測不可能さ」の本質を体現しているんです。その名も二重振り子

「え、振り子って規則正しく揺れるだけじゃないの?」と思ったあなた。その通り、普通の振り子は規則的です。でも、それが2つ連なると、話は一変。まるで生き物のようにクネクネと複雑な動きを始め、やがて完全にカオス(混沌)に陥ります。

ざっくり言うと、こういう話

二重振り子を、上司と部下の関係に例えてみましょう。

上司(上の振り子)がちょっと方針を変える(揺れる)と、それが部下(下の振り子)に伝わり、部下は大きく動揺します(大きく揺れる)。その部下の動揺が、今度は上司への報告(フィードバック)となって返ってきて、上司もさらに動揺する…。

こうして互いが互いに影響を与え合い、増幅し合う状態が生まれます。これが「非線形の連成振動」というやつで、カオス運動の根源です。

一言で表すと「二重振り子の動きは、2つの振り子が互いに“力の綱引き”をしながら、エネルギーをやり取りするバランスで決まる」

数式を読み解く(怖くない)

この複雑な動きをコンピューターに計算させるには、やはり数式が必要です。でもご安心を。ここで活躍するのは、物理エンジニアの必須ツールラグランジュ方程式です。

難しい力のやり取り(どのヒモがどれだけ引っ張っているか)を直接考えず、運動エネルギーT位置エネルギーVという「エネルギーの差」からスマートに運動を導き出せる、魔法のような式です。

\frac{d}{dt}\left( \frac{\partial L}{\partial \dot{\theta}_i}\right) - \frac{\partial L}{\partial \theta_i}= 0 \quad (i=1,2), \quad L = T - V

この式が言っていること:

  • $L = T - V$ (ラグランジアン): 要するに「動きの勢いから、位置の高さによるブレーキを引いたもの」。
  • 左辺全体: 「その“勢いの差”の時間変化が、今の姿勢に対してどう働くか」を計算しています。
  • 式が=0: エネルギー保存則(外部から力が加わらない限り、このバランスは保たれる)を表しています。

つまり、運動エネルギーと位置エネルギーの綱引きの結果として、次の瞬間の動きが決まる、というシンプルな原理なんです。この式を2つの振り子($\theta_1$, $\theta_2$)について立てると、連立方程式が出来上がります。

シミュレーターで遊んでみよう

数式より百聞は一見に如かず。さっそくツールを開いて、実験してみましょう。

🔬 実験1: 「初期値敏感性」を体感せよ

  1. 適当な設定でシミュレーションをスタート。
  2. 「2軌道比較」ボタンを押す。
  3. 下の振り子の初期角度θ₂を、たった0.1度だけ変えてみる(例: 90.0° → 90.1°)。

どうですか?最初は同じ軌道をたどっていた2つの振り子が、数秒後にはまったく別の動きを始めます。これがカオス理論の核心「バタフライ効果」の体感です。ほんのわずかな違いが、時間とともに指数関数的に増幅されるんです。

🔬 実験2: 重さのバランスで動きがガラリと変わる

  1. 下のおもりの質量m₂小さく(例: 1.0 → 0.2)してみる。
  2. 次に、大きく(例: 1.0 → 5.0)してみる。

m₂が小さいと、下の振り子は上の振り子に振り回されるだけの“お飾り”状態。動きは比較的単調です。しかし、m₂を大きくすると、下の重りが主導権を握り、上の振り子を激しく引きずり回します。質量バランスがシステムの“主従関係”を決めるのです。

🔬 実験3: エネルギーは本当に保存してる?

  1. 「減衰」パラメータを0に設定。
  2. シミュレーションを走らせ、グラフエリアの「エネルギー」タブを確認。

理論上、減衰(空気抵抗など)がなければ、運動エネルギーと位置エネルギーの和(全エネルギー)は一定です。グラフを見ると、確かにほぼ一定の値をキープしていますね。しかし、「数値誤差」 によってごくわずかに増減しているのが見えるはず。これがシミュレーションの精度の限界であり、現実の計算で常に付きまとう課題です。

現場でハマるポイント

この計算、実務(例えばロボットアームの設計)で使う時、一番気をつけることは?

  • 落とし穴1: 数値積分法の選択
    このツールで使われているRK4(4次ルンゲ・クッタ法) は精度が高いですが、計算コストも高い。リアルタイム制御が必要な現場では、精度と速度のトレードオフを考え、もっと軽量な積分法を選ぶことがあります。
  • 落とし穴2: 特異点への対応
    振り子が真上(θ=180°)に来るような姿勢では、式の中にゼロで割る項が現れ、計算が爆発(数値不安定)することがあります。実用的なシミュレーターには、こうした特異点を回避・処理するロジックが必須です。
  • 落とし穴3: エネルギーのドリフト
    実験3で見たように、長時間シミュレーションすると数値誤差でエネルギーが少しずつ増えたり減ったりします(ドリフト)。特にゲーム物理などでは、このドリフトが蓄積してオブジェクトが暴走しないよう、“エネルギー補正” を行う処理が入ることもあります。

もっと深く知りたい人へ

二重振り子は、カオス理論・非線形力学への最高の入り口です。次に調べるなら「リアプノフ指数」(カオスの“予測不可能さ”を数値化する指標)や「ポアンカレ写像」(複雑な運動を簡潔に分析する手法)がオススメ。この小さなモデルが、実は宇宙スケールの軌道計算から経済モデリングまで幅広く応用される基礎になっていることに驚くはずです。

まとめ

今回のポイント:

  • 二重振り子のカオスは、2つの振り子が互いに力をフィードバックし合う「非線形性」 が生み出す。
  • ラグランジュ方程式は、面倒な力を直接計算せず、エネルギーの差からスマートに運動を導く便利な道具。
  • シミュレーターで実際にパラメータをいじる体験は、教科書を読むだけの10倍の理解をもたらしてくれる。

「カオス」や「非線形」と聞くと難しそうですが、その本質はこのように単純なモデルに宿っています。まずは遊びながら、その面白さに触れてみてください。

二重振り子(カオス)シミュレーター — ブラウザで即座に動作、登録不要

NovaSolverでは700以上の工学シミュレーターを無料公開中 👉 一覧はこちら

0
0
0

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
0
0

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?