二重振り子のシミュレーション

二重振り子とは？？

物理部らしいことしてみた

二重振り子はカオス現象ですが、運動方程式をラグランジアンに変換して、Runge–Kutta法で微分方程式の近似解出したらシミュレーションできます(言ってみたかった) pic.twitter.com/fuN8Pdsiui

— 甲陽学院高等学校 物理部 (@koyobutsuri) 2018年6月25日

普通の振り子の先にもう一つ振り子をつないだだけの簡単な構造なのですが、これがまた興味深い。いわゆるカオスな運動をします。

カオス理論は、力学系の一部に見られる、数的誤差により予測できないとされている複雑な様子を示す現象を扱う理論である。ここで言う予測できないとは、決してランダムということではない。その振る舞いは決定論的法則に従うものの、初期値鋭敏性ゆえに、数値解析の過程での誤差によっても、得られる値と真の値とのずれが増幅される。そのため予測が事実上不可能という意味である。(Wikipediaより引用)

もう少し噛み砕いて説明すると、振り子の動きを運動方程式として書くことができる、つまり振り子はある法則に従って動いているのですが、その運動方程式を使って二重振り子の動きを予測することは不可能であると言うことです。

どれだけ高性能なスパコンを使っても、数十秒後には誤差が指数関数的に増加して、バラバラの動きとなってしまいます。

そのとっても無意味な運動方程式を立ててやる！！というのが今回のお題です。

どうでもいいですけど、二重振り子の動きとオタ芸の腕の動きってなんか似てません？はい、どうでもいいですね。

ひたすら計算

二つのおもりで構成される二重振り子について考えてみたいと思います。

おもり1は質量m1で紐の長さはl1で角度をθ1、座標を(x,y)
おもり2は質量m2で紐の長さはl2で角度をθ2、座標を(X,Y)とします。
そうすると

各座標の速度は時間tで微分して

おもり１の速度をv、おもり２の速度をVとすると