二重振り子エクスプローラー

結合運動における決定論的カオスとバタフライ効果を探る

インタラクティブシミュレーションを読み込んでいます...

二重振り子はなぜ予測不能になるのか

これはランダムではなく、決定論的カオスです。わずか0.1°だけ異なる2つのシミュレーションを始めると、しばらくすると全く無関係に見えるようになります。その差はほぼ指数関数的に拡大するからです(正のリアプノフ指数)。非線形の結合によってエネルギーが2本の腕の間を絶えず行き来し、わずかな位相差が増幅されてまったく異なる未来の状態になります。方程式自体は正確です。崩れるのは予測可能な範囲のほうです。
この振り子は左右に揺れるだけではありません。どちらの腕がダンスをリードするか、絶えず気が変わっているのです。

二本の棒と重力だけ

モーターも乱数発生器も仕掛けもなく、あるのは二本のアーム、二つの関節、そして重力だけです。あの荒々しい動きは、各アームが互いを押し合う二つ目の関節から生まれます。静止点の近くで小さく弾ませてみてください。ほぼ普通の振り子のように揺れます。アームが最高点を越えて回れるだけのエネルギーを与えて初めて、カオスが主役になります。
複雑さは組み込まれたものではなく、恥ずかしいほど単純なものから育つのです。

先端は重力より速く落ちうる

アームを水平に保って静止から放し、外側のおもりを見てください。一瞬、自由落下する球より速く、つまり g を超えて下向きに加速することがあります。剛体のアームがてこと鞭のように働き、内側のアームが外側を重力だけでは不可能なほど強く振り下ろすのです。これは、蝶番付きの棒が落ちる硬貨より先に床へ届く現象と同じです。
落下を始めるのは重力ですが、誰が先に着地するかを決めるのは幾何学です。

全プロセスの詳細解説

  1. 静止状態から放したときの2回の完全な振動の周期 7 ステップ

    m₁ = m₂ = 1、l₁ = l₂ = 1、g = 9.81、減衰 c = 0 のもと、θ₁ = 20° および θ₂ = 10° の静止状態から放つ。これほど小さな角度ではカオスは生じず、運動は2つの整った振動の重ね合わせとなる。それぞれの周期を求め、パネル上のサイクルマークがそのどちらをカウントしているかを判定せよ。

    この問題をツールで開く

    1. 上記の運動方程式を簡略化する。角度が小さいため sin θ は θ に、cos 2Δ は 1 になり、速度の2乗項は小振幅において2次の微小量となるため消去される。両方の分母は 2m₁ + m₂ − m₂ = 2 に整理され、残るのは線形でありながら連成した運動である。どちらのおもりも、他方を押し動かすことなしには運動しない。

      θ¨1=−2kθ1+kθ2,θ¨2=2kθ1−2kθ2,k=gℓ\ddot\theta_1=-2k\theta_1+k\theta_2,\qquad \ddot\theta_2=2k\theta_1-2k\theta_2,\qquad k=\frac{g}{\ell}
    2. 何をも求める前に、まず比較の基準を置く。長さ 1 m の単振り子は ω₀ = √(g/l) で振動し、パネルが周期のヒントとして表示し、各サイクルマークに刻んでいる 2.01 s という値を与える。この値を覚えておいてほしい。以下で得られる計算結果のどれも、この値とは一致しないからである。

      k=9.811=9.81 s−2,ω0=k=3.132 rad/s,T0=2πω0=2.01 sk=\frac{9.81}{1}=9.81\ \mathrm{s^{-2}},\quad \omega_0=\sqrt{k}=3.132\ \mathrm{rad/s},\quad T_0=\frac{2\pi}{\omega_0}=2.01\ \mathrm{s}
    3. 固有モードとは、両方の角度が1つの共通した振動数で振動し、一定の比を保つ運動のことである。θᵢ = aᵢ cos ωt を代入すると、2つの運動方程式のそれぞれから、互いに逆の側からたどり着く形でその比 a₂/a₁ が得られる。

      θi=aicos⁡ωt ⟹ (2k−ω2)a1=k a2,(2k−ω2)a2=2k a1\theta_i=a_i\cos\omega t\ \Longrightarrow\ (2k-\omega^{2})a_1=k\,a_2,\quad (2k-\omega^{2})a_2=2k\,a_1
    4. 1つの比を表すこれら2つの式は一致しなければならず、これらを等しいとおくことで両方の振幅が消去される。残るのは ω² に関する2次方程式であり、2次方程式には 2 つの解が存在する。第2の周期はここから導かれる。

      (2k−ω2)2=2k2⟹ω2=(2∓2) k=5.747, 33.49 s−2(2k-\omega^{2})^{2}=2k^{2}\quad\Longrightarrow\quad \omega^{2}=(2\mp\sqrt2)\,k=5.747,\ 33.49\ \mathrm{s^{-2}}
    5. したがって、この振り子には 2.62 s と 1.09 s という 2 つの周期が存在し、サイクルマークの 2.01 s はそのどちらでもない。それは両者の間に位置し、いずれのモードにも属していない。

      ω−=2.397,  ω+=5.787 rad/s⟹T−=2.62 s,  T+=1.09 s\omega_-=2.397,\ \ \omega_+=5.787\ \mathrm{rad/s}\quad\Longrightarrow\quad T_-=2.62\ \mathrm{s},\ \ T_+=1.09\ \mathrm{s}
    6. それぞれの解は独自の値と形態をもたらす。遅い解からは a₂/a₁ = +1.414 が得られ、両方のおもりが同じ方向に揺れ、下のおもりの振幅が 1.414 倍大きくなる。早い解からは −1.414 が得られ、おもりは逆方向に揺れ、その間の関節がハサミのように交差する。

      a2a1=2k−ω2k=±2=±1.414\frac{a_2}{a_1}=\frac{2k-\omega^{2}}{k}=\pm\sqrt2=\pm1.414
    7. 初期値 20° と 10° はどちらの形態とも一致しないため、両方のモードが同時に動作している。これを遅いモードの振幅 A と早いモードの振幅 B に分割すると、上のおもりは A + B = 20°、下のおもりは √2(A − B) = 10° を満たす必要があり、初期振幅のうち遅いモードが 13.54°、早いモードが 6.46° を占める。

      A+B=20∘,  2 (A−B)=10∘⟹A=13.54∘,  B=6.46∘A+B=20^\circ,\ \ \sqrt2\,(A-B)=10^\circ\quad\Longrightarrow\quad A=13.54^\circ,\ \ B=6.46^\circ

    解答

    2.62 s と 1.09 s であり、その比は厳密に 1 + √2 である。無理数の比は、これら 2 つのモードが二度と周期を揃えて同期しないことを意味し、したがって下のおもりが描く図形が完全に繰り返されることはない。増してや、これはカオスが一切含まれない、おとなしい小角度の場合である。これにより、2.01 s のサイクルマークの正体も確定する。それは 2 つの周波数を持つ運動にあてがわれた単振り子の定規であり、各刻みで記録される角度が 20.0° や 10.0° に戻らずドリフトし続けるのはそのためである。離す角度を 145° に広げると、2 つの周波数すら消滅する。なぜなら sin 145° = 0.5736 であるのに対し、角度自体は 2.531 rad であり、4.4 倍もの誤差が生じるため、この導出全体が依拠している線形化にはもはや述べるべきことが何もないからである。

    T−T+=2+22−2=1+2=2.4142…\frac{T_-}{T_+}=\sqrt{\frac{2+\sqrt2}{2-\sqrt2}}=1+\sqrt2=2.4142\ldots
参考文献 (1)