経緯
物理以外の参考書等(例えば制御工学)で拘束条件付きの問題に対して最初から一般化座標系を導入して拘束条件なしのLagrangeの運動方程式を立てている様子を見て、それが直角座標系に於けるNewtonの運動方程式と等価であることを納得したくて暫く考えていた。その結果をここに書き残しておく。
はじめに
まず直角座標系に於いて拘束条件付きのLagrangeの運動方程式がNewtonの運動方程式と同値であることを理解する。そして、可逆な座標変換、所謂「一般化座標系」への変換を行ってもLagrangeの運動方程式が表面的には形を変えない(Lagrange関数の数式は変わるが、それを相変わらずと書く限り形を変えない)ことを示す。そして拘束条件を表す制約式に現れる一般化座標変数の個数が最小となるような変換を選ぶと、一般化座標変数の一部が拘束条件から容易に求まることと、残りの一般化座標変数について拘束条件が消える可能性について述べる。最後に、そのことを確認できる例題を一つ紹介し、数値計算結果を示す。
問題設定
を質点の個数とし、を番目の質点(以下「質点」と呼ぶ)の質量とする。各質点に働く力は保存力、および質点の速度と直角に働く(つまり質点に対して仕事をしない)拘束力のみであるとする。ここでの拘束力とは、質点がある曲面や曲線に沿ってしか運動できず、かつそれら曲面,曲線を構成する物体(拘束物)と質点の間に摩擦力が生じない状況で、質点が拘束物から受ける力を想定する(例えば滑らかなワイヤに通されたビーズの運動や、半径が異なる、中心を共有する2つの滑らかな球面の間に挟まった鉄球の運動などを想像するとよい)。
直角座標系では質点1個あたり3つの座標変数があるので、系全体では個の座標変数がある。これらをとする。質点の座標はとなる。個々の質点の座標を分けて考えるのではなくにまとめて表現することで、質点系を次元空間を運動する1個の高次元の質点として扱える。このようにしても数学的には元の運動方程式と変わらない。
保存力
直角座標系に於いて、質点系に保存力をもたらすポテンシャルをとする(CはCartesianのC)(空間に固く固定された電荷などポテンシャルの原因が複数あるときはそれらの和をとすればよい)。質点系が受ける保存力はとなる。
Lagrange関数
直角座標系に於けるLagrange関数をとすると
ここには「質量行列」と呼ばれる対角行列であり、である(直角座標系では対角であるが、他の座標系ではそうとも限らない)。
拘束条件
質点系の拘束条件は複数あり得る。をその個数とする(CはConstraintのC)。これらの制約が、直角座標系において級の関数を用いてと表される場合を考える(の添え字のCはCartesianのC)。但し各制約式は、それが満足される領域でとなるように選ぶものとする(後で判るが、勾配が消えると計算上都合が悪い)。例えば質点1に対する制約条件として直線への拘束を考えると、制約条件を満たす領域で左辺の勾配がになる。しかし同じ拘束条件を、平面を表す2つの制約式の連立として表現すれば勾配はそれぞれとなる。同様に曲線も曲面の交線として表現する。今回は不等式制約(物理的には壁に進入しない等に相当する)は扱わない。冒頭で仮定したように、質点が曲面から受ける拘束力は曲面に垂直、すなわちに平行である。
直角座標系に於けるLagrangeの運動方程式
以上に述べたことから、直角座標系に於ける拘束条件付きのNewtonの運動方程式は次のようになり、これは拘束条件付きのLagrangeの運動方程式(Lagrange’s equations of the first kind (第一種Lagrange方程式)とも呼ばれる)と同値である。
ここには質点系の位置と加速度, 質点系に働く保存力に依存する、究極的にはの関数である。質点が拘束条件のレールからはみ出ないようにに平行な力が働いて、質点系に働く慣性抗力(加速度と質量の積の-1倍)と質点系に働く保存力をキャンセルしていることは確かだが、いつどんな大きさで働くのかは自明ではない。そのためはと同時に解かねばならない。
可逆な座標変換をしてもLagrangeの運動方程式の形が変わらないこと
産業応用上解く必要に迫られる力学の問題は解析的に解けないものが殆どであり、数値解析に頼ることになるのだが、その場合でも前出の式,をできるだけ解きやすくなるように座標変換することは必要である。拘束条件を利用して解くべき座標変数の個数を減らせれば問題の次元が下がって扱いやすくなる。
との間の可逆な級の座標変換を考える。つまりはの級関数として表せて、かつも級関数として表せる状況を考える。
まずはの関数になる。なぜなら となり、はの関数なのでがの関数となるからである。全く同様にしてはの関数となる。後々の計算の為に座標系上でのLagrange関数を次式で定義しておく。 とは関数値は等しいが、写像としての定義が異なる。物理学ではそういうものを区別せずに、異なる座標系で同じ記号を使いまわすのが慣例だが、ここでは式変形から手品を排除して全てをはっきりと見る為に敢えて執拗に数学の形式に拘ることにする。(「いやしかし、は既にの関数としてと表しているのに、の関数としてとして表すのはいい加減ではないか?」という指摘もあろう。その通りで、有言実行するには新たに記号を定義して区別するしかないが、それでは却って式の可読性が下がる。結局、厳密さによる混乱の回避の恩恵と式の可読性のトレードオフで戦略を決めている。)をで表したのと同じ理屈で次式も成り立つ。
座標変換によってがどうなるか考える。
ここで次正方行列をで定義すると、これは関数のJacobianであり、関数は可逆だからは正則である。同様に次正方行列(正則とは限らない)をで定義すると次式を得る。
次に座標変換によってがどうなるか考える。 ここで より であるから次式を得る。 これをさらにで微分すると
次に座標変換によってがどうなるか考える。を定義したときと同じように、とするとである。 よって次式が成り立つ。
以上の式,,をに適用すると次式を得る。 となることを示す。まず 次に である。左から2つ目の等号は、がの関数である(はの関数なのでそれをで偏微分してもの関数である)ことから従う。右から2つ目の等号はがの級関数であと仮定しておいたことから従う。以上2つの式からであるので結局次式を得る。 先述の通りは正則なので両辺に左からを掛けて次式を得る これと制約条件を合わせて、座標変換後の拘束条件付きのLagrangeの運動方程式は次式である。
Lagrange関数と拘束条件式の記号が変わった以外は直角座標系のときと同じ形の式である。
座標変換を工夫すると制約なしのLagrangeの運動方程式が得られる可能性について
式は行の連立微分方程式であるが、の取り方を工夫することでの幾つかの行を0にできる。のうちに現れるものの個数が最小になるようなを見つけると都合が良い。もしがに現れなければの行目は0である。そしてに現れるの成分が簡単に求まる場合、残った成分について拘束なしのLagrangeの運動方程式(Euler–Lagrange方程式とかLagrange’s equations of the second kind (第二種Lagrange方程式)と呼ばれる)を解けば良いので問題が格段に扱いやすくなる。
物理学を道具として用いる分野の技術書等(例えば制御工学)では普通、最初から拘束条件が消えるような一般化座標系への変換を考えてLagrange関数を構築して拘束なしのLagrangeの運動方程式を鮮やかに立てているが、これは「拘束条件が消えるような一般化座標系への変換を考え」る過程で自動的に「に現れるの成分が簡単に求ま」って、そして「がに現れなければの行目は0」となって、最初から拘束なしのLagrangeの運動方程式に辿り着ているのである。
問題例
とし、質点の質量をとする。拘束条件として潰れた楕円形の閉曲面を考える。軸の負の向きに重力加速度が働くとする。
拘束条件の可視化
Lagrange関数は次式である。
直角座標系に於ける拘束条件付きのLagrangeの運動方程式を具体的に書き下すと次のようになる。
人間に解けるような問題ではないから計算機で数値計算するとして、拘束条件が厄介である。そこで次のような座標変換を考える。
そうすると拘束条件は
となりが決まる。そしては次式になる。
そしてLagrange関数をの関数で表してで偏微分してLagrangeの運動方程式を立てるのだが、実はここから先の計算は恐ろしく大変である。筆者の手作業では1日掛かっても終わりそうにないのでMathematicaにやらせたが、Lagrangeの運動方程式に至ってはここに書けないほど長い数式になった。簡単に表せる部分のみ具体的に表示し、他は記号的に書くと次のようになる。
式の3行目はにを代入する前に計算しなければならないことに注意する(1,2行目はを先に代入してから計算しても構わない)。そもそも式はが変数であるとして導出されたのだから、の計算前にを特定の値で固定してはいけない。例えばという写像があって、何らかの条件からであることが解っているからといってとするのはおかしい。なのでであってもである。「関数の定義」と「関数値」は別物であって、「関数の定義」は「関数値」よりも多くの情報を保持している。偏微分は「関数の値」ではなく「関数の定義」に対して定められている。
式の上2行はの要素が0になっており、の影響を受けない。既には解っていて未知の座標変数はの2つだけであるから式の上2行からこれらが唯一に決定されるように思えるし、そうであって欲しい。そして実際にそうでなくてはならない。仮に式の上2行だけでが唯一に決定されないとする。これ以外に利用できる情報は3行目の式しか残されていなが、それを考慮に入れたところで、が同時に解くべき未知変数に追加されるため、結局解けないのである。古典力学的に正しく定義された問題に対しては解が唯一に決まることは疑いようのない事実であり、「が唯一に決まらない」状況は物理的に許されず、式の上2行だけでが唯一に決定されるより他にあり得ないのである。そして3行目からが決定される。
先にも述べた、解析力学の理解を前提とした技術書では座標変換で一般化座標変数の一部を制約条件から事前に求めてしまい(制約条件と同じ個数の座標変数が代数的に決まることを想定する)、それらを最初から定数とした上で、Lagrange関数を残りの一般化座標の関数と見做して拘束なしの運動方程式を解いている。こうしてを表に出さずに解軌道に辿り着く。
この一般化座標系に於けるLagrangeの運動方程式をMathematicaで解いた結果を元の座標系に変換して可視化したものが次の図である。
数値計算で求めた軌道
初期条件はすなわちである。この位置に固定されていた質点が時刻0に解放されて楕円曲面に拘束されながら重力によって落下し、底に辿り着くと勢いで壁を登っていき、速度が落ちて再び落下していくのを繰り返す様子が見える。
以下は、このシミュレーションを行ったMathematicaノートブックである。