二重振り子

二重振り子の運動 が知りたい。

多体系の拘束力 を求めたい

二重振り子（右図）は、前章で扱った振り子に、もう1つおもりを追加したものである。おもりが2つあるので、第6章の場合と同様に、運動方程式をまとめて大文字で表す：（ は重力、 は拘束力） 二重振り子の拘束条件は、「ひもの長さ （同図）が一定」である： これは、「2つのおもりと2つの拘束条件」があるような問題となっている。

運動方程式()を解くために知る必要があるのは、これまで同様、拘束力 と、初期値 に対する拘束条件である。 に対する拘束条件は、今回も となる。問題は、拘束力 である。第6章でも述べたが、複数の物体がある系（＝多体系）では、各物体に対してダランベールの原理（＝拘束力は拘束面に垂直）を適用しただけでは、 の向きが決まらない。第5章では、作用・反作用の法則を用いてダランベールの原理を多体系に拡張したわけだが、今回も、それと同じような考察が必要となるはずである。

そこでこの章では、おもりの運動 を計算するために、以下の3つの節に分けて議論を行う： 拘束力の向きが決まらない多体系のダランベールの原理の導出物体の運動の計算 なお、8.2節では、時間に依存する拘束条件 の場合も考慮する。

8.1拘束力 の向きが決まらない

第6章でボール同士の衝突を考えた際、各々のボールにダランベールの原理を適用すると、未定乗数の数が多くなりすぎるという問題があった。この節では、二重振り子の場合にも同じ問題が生じることを確認する。

に対する拘束条件：式()

まず、初期値 に対する拘束条件は、これまでと同様に である。これらは、前章と同様に となる。

後は、拘束力 が決まれば、運動方程式から運動が計算できる。 は、 の2階微分： を満たしさえすれば、拘束条件 は常に成立する。式()の成分は2つしかないので（式()）、これだけでは を完全に決めることはできない（ は6成分なので、6つの条件が必要）。よって、拘束条件以外の物理的な条件を追加する必要がある。ここまでは、前章と同じである。

各々の物体に対するダランベールの原理より式()→未定乗数が多すぎる

拘束力 に課すべき条件()に追加するものとして、各々の物体にダランベールの原理を適用することを考える。

まず、拘束力 は、 番目の拘束条件 に由来する拘束力 を足し合わせたものになる（力の加法性による）： ここで、物体ごとにダランベールの原理を適用することを考える。 番目の物体が 番目の拘束条件により受ける拘束力 の向きは、ダランベールの原理により拘束面に垂直、即ち、 と平行である： は未定乗数であり、 は での偏微分 である。

これらをまとめて書くと、式()は以下のようになる： ただし、 は、 を で微分したものである： 最後の式は、括弧をつけて としたほうが誤解がないが、後で見づらくなるので括弧は省略する。

式()の未定乗数 を全て決めれば拘束力 が決まる。しかし、式()の の数は「物体の数」×「拘束条件の数」なのに対し、式()は「拘束条件の数」しか存在しないため、 を一意的に決めるには足りない。

ダランベールの原理は式()のように拡張できるのではないか

この状況は、ボール同士の衝突（第6章）の場合と同じである。ボール同士の衝突では、作用・反作用の法則により、両方のボールの未定乗数が等しくなり（ ）、拘束力は という形に書けたのだった。もし今回も、同様の関係式 が成り立つとすれば、式()は、以下の形で書けることになる： そうすれば、未知数は のみとなり、その成分数は、拘束条件()の数と一致するので、 が決まることになる。

式()は成立するだろうか。確かに、拘束運動と衝突は似ており、実際、衝突が連続的に生じているのが拘束運動であると見なすこともできるので、同様の式が成立している可能性はある。しかし、第6章で扱ったのは、球対称なボールの衝突という特殊な状況に過ぎない。式()は、考え得る中で最も単純な形をしているが、果たしてそう上手くいくのだろうか。これを次節で考える。式()の微分を実際に実行して、振り子の各ひもに対して作用・反作用の法則を使うと、確かに式()が成り立つことが言える。しかしここでは、二重振り子に限らない一般論として示したい。

8.2多体系のダランベールの原理→ の導出

この節では、式()が確かに成立することを見た後、拘束力 を与える公式()を示す。

拘束力 は、運動エネルギー()を変化させない

拘束された運動においても、ボール同士の衝突（第6章）における弾性衝突と同じように、「運動の勢い」が保存していると考えられる。実際、2重振り子の例をイメージしても分かるように、全体の運動がどんどん激しくなっていったり、減衰して止まったりすることは考えづらい。とりあえず簡単のため、外力 （＝拘束力以外の力、例えば重力や摩擦力）がゼロで、拘束条件も時間に依存しない場合を考えている。拘束条件()を満たすような勝手な拘束力 を取ると「運動の勢い」は保存しなくなるので、「運動の勢い」の保存則は、何らかの新しい条件を与えることになる。

では、「運動の勢い」とはどのように定量化されるのだろうか。それを探るため、ダランベールの原理()が成り立つと仮定してみて、候補となる保存量が存在するかを調べることにしよう。まず、式()より、運動方程式()は以下の形になる：（ としている） これが正しいのであれば、この式から、「運動の勢い」に対応するような保存量が導けるはずである。「運動の勢い」は、拘束条件 とは無関係に定義出来るはずなので、上式から を消すことを考える。これは簡単で、 を左から掛ければよい。実際、拘束条件 の微分（ は時間によらないとする） を用いれば、式()の右辺が消える：（公式 を使う） 右辺 よって、左辺のみが残る： この式の左辺を変形すると 左辺 となる。これが となるのだから、赤字部分を ： とおくと、 は保存する： を運動エネルギーという。

それらしき量が出てきたが、この は本当に「運動の勢い」を表しているのだろうか。少なくとも、一方のおもりが静止している場合には、 が保存すれば、もう一方のおもりの速度の大きさが保存することになるのでもっともらしい。第6章でボール同士の衝突を扱った際、弾性衝突は「運動の勢い」が保存するような衝突であった。従って、ボール同士の衝突の際にも が保存することが期待されるが、これは実際に保存していることを示すことができる（以下の【8.2-注1】）。そこで、 は「運動の勢い」を表していることを認めることにしよう。つまり、拘束力 は、運動エネルギー を保存するように働く。今は物体の個数が2つだけだが、式()の項を増やすことで、より多くの物体がある場合にも自然に一般化できる。

ところで、例えば、重力下での振り子の場合、おもりの速度の大きさは時間とともに変化するので、運動エネルギーは保存していない。しかし、全体としては周期的な運動を繰り返すのだから、運動が減衰しているわけではないように見える。この場合でも、何らかの保存則が成り立っているではないだろうか。これは実際に正しく、重力下の運動では、運動エネルギーに「ポテンシャルエネルギー」と呼ばれる項を加えたものが保存する（以下の【8.2-注2】）。

拘束力 は拘束条件のミクロな破れによって生じる

上述のように、拘束力 は、運動エネルギー を保存するように働く。しかし、 という1つの条件だけでは、式()からダランベールの原理()：（再掲） は導けない。 の性質について、もう少し考察する必要がある。

まず、拘束条件 は、ミクロに見れば、運動に応じてほんのわずかだけ実際には破れているはずである。拘束力は本来、そのミクロな破れを押し戻そうとする力である（右図の ）。これは例えば、硬いバネを変形させた時、変形に応じて元に戻ろうとする力が働くのと同様である（バネを非常に硬くしていったものが拘束だと考えてもよいだろう）。よって、拘束面上の同一点において 拘束条件のミクロな破れ方が同じであれば、は等しい と言えそうである。これを認めることにしよう。これは、 が、拘束のミクロな破れの量にのみ依存し、物体の速度や外力には直接依存しないということである。もちろん、速度や外力が大きいと も大きくなる傾向があるが、それは、その方が拘束がより大きく破れるからであって、速度や外力に直接依存しているからではない。

多体系のダランベールの原理：式()

以上で予備的な考察が終わったので、ここからは、外力 があってもよいとする。拘束面上の1点 において、ある拘束力 を与える「拘束条件の破れ」を とおく。 が等しい全ての に対し、拘束力は同じ値 を取ることになる。速度 は、拘束条件()：（再掲） を満たす全ての を取ることができる。そして、これら全ての に対して、 は運動エネルギー を変化させない。式()の導出で見たように、 が変化しないのであれば、 は と垂直となる：

以上により、拘束力 に対する条件をまとめると、「式()を満たす全ての速度 に対して、式()が成立する」、となる。そのような は、式()の形以外にはあり得ない（以下の【8.2-注3】参照）。以上により、ダランベールの原理()が導かれた。

この議論の仕組みについて、見直しておこう。そもそも、 だけからダランベールの原理()が得られなかったのは、1つの運動に対する のみを考えていたからである。性質()を使うことにより、同じ に対し、拘束条件を満たす全ての を考えることができるようになった。それら全ての に対して、 が運動エネルギーを変化させないという条件から、式()が得られたのである。

拘束条件が時間変化する場合も、ダランベールの原理()は成り立つ

ここまでは、拘束条件が時間変化しない場合を考えた。時間変化する拘束条件 の場合であっても、拘束力 は、拘束の破れにのみ依存し、拘束面の移動速度には直接依存しない。よって、拘束の時間変化を直接考える必要はない。ただしこの場合は、速度 に拘束面の移動速度が含まれるため、 と は垂直でなくなる。そのため、式()はそのままでは成り立たないが、代わりに、 を、拘束面に沿ったベクトル に置き換えたものが成立する：

また、 が拘束面上にあるという条件は、以下のように書ける： は、考えている時刻における拘束面の形状のみから決まり、移動速度とは無関係である。 が時刻に依存しない場合には、 と は、取り得る方向が一致する。

以上により、拘束力 に対する条件は、「式()を満たす全ての に対して、式()が成立する」、となる。よって結局、 を に置き換えるだけの同じ議論により、 が時間変化する場合もダランベールの原理()が成り立つことが分かる。

拘束力 ：式()

拘束力 を求めるには、ダランベールの原理()の未定乗数 を求めればよいわけだが、「 の成分数」と「拘束条件 の要素数」が等しいので、拘束条件と連立してやれば、 は一意的に決まる。よって、式()の拘束力 が確定する。

実際に を計算すると、 が時間に依存する場合も含めて、以下の【8.2-注4】の式()のようになる。これが、一般的な拘束力 の公式である。運動方程式に代入したものを書き下しても良いが、きれいになるわけでもないので省略する。前章の【7.2-注3】のように、射影行列を使って運動方程式を書くこともできる。

8.3物体の運動 の計算

以上で、必要な議論がそろった。この節では、まず、拘束条件がある場合の運動 の計算方法についてまとめる。その後、具体的な例として、ディアボロ（2物体1拘束）、2重振り子（2物体2拘束）、シーソー（2物体4拘束）について、数値計算を行う。

物体の運動 の計算方法

複数の物体 に、複数の拘束条件 が課されているとする（ は拘束条件が時間に依存する場合に必要）。物体の運動 を計算するには、まず、 において拘束条件： を満たすような初期値 を設定する。その後の時間発展は、運動方程式 によって計算できる。ただし、 は重力などの外力であり、拘束力 は式()で与えられる。

【例題1】ディアボロ：初期条件()、運動方程式()

まず、2つの物体に1つの拘束条件が課せられている例として、右図のようなディアボロの運動を考える。 は固定点であり、ひもの全長 が一定となる（ひもはたるまない）。また、重力が働いているとする。

拘束条件 およびその微分は、以下のようになる： 公式（【注】）の計算：【注】

これを用いると、初期値 が満たすべき拘束条件()は となる。運動方程式()についても、 に重力を代入し、拘束力 に式()を代入して計算していくと、以下のようになる：