資料紹介
連続体の解析力学
ひもに解析力学が使えるか。
前回と同じモデル
前回はニュートンの運動方程式からひもの運動を論じた。 では、ラグランジアンを使った形式でひもの運動を論じる事が出来るか、というのが今回のテーマである。
まずは前回と同様、ひもは質点の集まりだというモデルから始めよう。 ラグランジアン L は L = T - V で表せたのだから、
と置けば良さそうだ。 第 1 項は「全質点の運動エネルギー」であり、そこから第 2 項の「全ての質点の間のポテンシャルエネルギー」が引いてある。 一つのラグランジアンのみで複数の質点の運動を全て表し切れてしまうというのが、この理論形式の強みの一つであるのだろう。 まぁ、問題の本質をすべてラグランジアンとして暗号化して押し込めてしまっただけだという見方もできなくもない。
ところで、この L の第 2 項目は、なぜこのような形式で書き表すことが出来ているか分かるだろうか。 あまり力学の計算に慣れていないと、こんなことでもつまづきそうだ。 前回の議論では、隣どうしの質点の、それぞれの変位の差を Y としたとき、 Y に比例した復元力 -kY が働くことを確認した。 念のため言っておくと、ここで書いた係数 k の正体は T/a なのだった。 T は張力で、a は質点間の距離である。 この力 -kY を Y で積分したものはエネルギーを意味するが、そうやって今回の式の形を得ているのである。 こまごまとした話を前回で済ませておいたお陰で、説明が楽に進んで助かる。
この L を次のような N 個のラグランジュ方程式に当てはめてやれば、 N 個の質点の運動方程式が得られてくるはずだ。
簡単な事なので念のため具体的に計算して確かめておこうか。 いや本当の事を言うと、今から出す結果を今回の話の後の方で使いたいから今のうちにやっておくのである。 まず、第 1 項の (d/dt) の後にある偏微分であるが、これは一般化運動量 pi と呼ばれている部分だった。 この計算は簡単だ。
そして (1) 式の第 2 項にある偏微分は一般化力 Fi と呼ばれているのだった。 yi で偏微分すればいいだけだが、 L の中に和の記号があって yn と yn+1 があるので、 n = i と n = i - 1 の場合の 2 つの項が関係する事に気を付けないといけない。
これらを組み合わせれば、前回作った運動方程式と同じものが出来上がる。
しかしここまでの話はまだ何ら新しい領域には突入していないことに気を付けて欲しい。 バネでつながれた質点の運動について論じているだけであり、ひもの運動ではない。
ラグランジアン密度
今回興味があるのは、これをどう改良したら「ひも」の運動を表すような方程式が得られるかということだ。 しかしいきなり魔法のような方法を期待してはいけない。 ひとまず前回と同じことを試してみようではないか。
つまり、l = na を保ったまま、n→∞ に、a→0 に、そして、それだけでは質量が無限大になって困るので、一個の質点の質量を M/n と置いてやることにしよう。 またここまでは、多数あるそれぞれの質点の上下の変位を yi として表してきたが、これからは質点が密に繋がっていると考えるので、 y ( x, t ) という関数として扱うことになる。 この考えで、先ほどの L を書き換えてやろう。
この式は試しに書き換えてみただけのものであり、まるで意味を成していない。 考え方の過程を示したいと考え
タグ
All rights reserved.
資料の原本内容
連続体の解析力学
ひもに解析力学が使えるか。
前回と同じモデル
前回はニュートンの運動方程式からひもの運動を論じた。 では、ラグランジアンを使った形式でひもの運動を論じる事が出来るか、というのが今回のテーマである。
まずは前回と同様、ひもは質点の集まりだというモデルから始めよう。 ラグランジアン L は L = T - V で表せたのだから、
と置けば良さそうだ。 第 1 項は「全質点の運動エネルギー」であり、そこから第 2 項の「全ての質点の間のポテンシャルエネルギー」が引いてある。 一つのラグランジアンのみで複数の質点の運動を全て表し切れてしまうというのが、この理論形式の強みの一つであるのだろう。 まぁ、問題の本質をすべてラグランジアンとして暗号化して押し込めてしまっただけだという見方もできなくもない。
ところで、この L の第 2 項目は、なぜこのような形式で書き表すことが出来ているか分かるだろうか。 あまり力学の計算に慣れていないと、こんなことでもつまづきそうだ。 前回の議論では、隣どうしの質点の、それぞれの変位の差を Y としたとき、 Y に比例した復元力 -kY が働くことを確認した。 念のため言っておくと、ここで書いた係数 k の正体は T/a なのだった。 T は張力で、a は質点間の距離である。 この力 -kY を Y で積分したものはエネルギーを意味するが、そうやって今回の式の形を得ているのである。 こまごまとした話を前回で済ませておいたお陰で、説明が楽に進んで助かる。
この L を次のような N 個のラグランジュ方程式に当てはめてやれば、 N 個の質点の運動方程式が得られてくるはずだ。
簡単な事なので念のため具体的に計算して確かめておこうか。 いや本当の事を言うと、今から出す結果を今回の話の後の方で使いたいから今のうちにやっておくのである。 まず、第 1 項の (d/dt) の後にある偏微分であるが、これは一般化運動量 pi と呼ばれている部分だった。 この計算は簡単だ。
そして (1) 式の第 2 項にある偏微分は一般化力 Fi と呼ばれているのだった。 yi で偏微分すればいいだけだが、 L の中に和の記号があって yn と yn+1 があるので、 n = i と n = i - 1 の場合の 2 つの項が関係する事に気を付けないといけない。
これらを組み合わせれば、前回作った運動方程式と同じものが出来上がる。
しかしここまでの話はまだ何ら新しい領域には突入していないことに気を付けて欲しい。 バネでつながれた質点の運動について論じているだけであり、ひもの運動ではない。
ラグランジアン密度
今回興味があるのは、これをどう改良したら「ひも」の運動を表すような方程式が得られるかということだ。 しかしいきなり魔法のような方法を期待してはいけない。 ひとまず前回と同じことを試してみようではないか。
つまり、l = na を保ったまま、n→∞ に、a→0 に、そして、それだけでは質量が無限大になって困るので、一個の質点の質量を M/n と置いてやることにしよう。 またここまでは、多数あるそれぞれの質点の上下の変位を yi として表してきたが、これからは質点が密に繋がっていると考えるので、 y ( x, t ) という関数として扱うことになる。 この考えで、先ほどの L を書き換えてやろう。
この式は試しに書き換えてみただけのものであり、まるで意味を成していない。 考え方の過程を示したいと考えている。 まだ和の記号 Σ を使っていたりして、質点の一つ一つを数えながら足している考えが残っているのが見て取れるだろう。 しかし今からは y ( x, t ) は連続関数であると考えようとしているのであり、 x 軸に沿って連続的に和を取ってやらないといけない。 つまり、積分記号に置き換えるべきなのである。
ところで、ただ Σ を ∫ に置き換えるだけでいいだろうか。 積分には dx という記号が付き物であるはずだが、これは単なる飾りではなくて、勝手に付け加えれば良いというものではない。 dx には「これから微小範囲 dx ごとに分けて和を取ろうと思うのだけれど、その範囲 dx はこれから無限に 0 に近付けるつもりです」という意味があるのだった。 ああ、なーんだ！ それはまるっきり a のことではないか。 ということで、a を dx と全く同じものだと見なす事で対応する。
前回と同じく、 lim 記号を付けないで変形するという無作法なことをしてしまったが、思考過程が良く分かる解説のためにそうしたのであり、そこに込められた精神だけを受け止めて欲しい。 このようにして、ラグランジアンの意味を何も変えることなく、連続的な場合へと拡張してやることが出来た。
ところで、この最後の式の積分の中身は「ラグランジアン密度」と呼ばれている。 これは今後も使うつもりがあるので、次のような記号で表しておくことにしよう。
ラグランジアン L の物理的意味がはっきりしないのだから、ラグランジアン密度についても、具体的な意味を考えて悩むだけ無駄というものである。 まさに、単位長さ辺りのラグランジアンという意味で納得しておくしかない。
汎関数微分の紹介
さて、ここからどうやったら、前回と同じ波動方程式を導いてやることができるだろうか。 同じ結果が得られてこそ、安心できると言うものだ。
まぁ、ラグランジュ方程式に代入すれば良い、と考えるのがお決まりのパターンであろう。 しかし、何の変更も無しにすんなり行くのかというのが大変心配になるところだ。 何しろ、ラグランジアン L は今や積分で表されているのである。 もはや多数の yi や i の関数ではなく、関数 y ( x , t ) の関数になっている。 このように関数全体の形によって値が決まるような関数を「汎関数」と呼ぶ。
一方、ラグランジュ方程式というのは、先ほども書いた (1) 式のことである。 これは L を yi や i で偏微分している形式になっているが、これをただ y ( x , t ) などで置き換えただけのものには意味があるだろうか？ L は x で積分してある形なので x の関数ではなくなっているはずだが、それを x や t を含むような関数で偏微分するなんて、一体それはどう計算すればいいというのだろう！？
実は、汎関数を関数で微分する方法が存在するのである。 それは「汎関数微分」と呼ばれている。 いや何のことは無い、それは第 3 部の「 ベルヌーイの問題提起 」のところで説明したのと全く同じ内容なのである。 「変分法」という言葉が「汎関数微分」の別名として使われることは普通に良く見かけるのだが、それは論理的に全く同じだからだ。 汎関数微分の詳しい説明は次回に予定している。
とりあえずはラグランジュ方程式に話を戻して、その意味を再確認しておこうか。 そしてその本質を壊さぬように、無限自由度の場合に拡張したものを作ってやればいいだろう。
まず (1) 式の第 1 項では L を i で偏微分しているが、これは i 番目の質点以外のことには目もくれないという意味が込められている。 そしてその結果を時間微分したものは、第 2 項の「 L を yi で偏微分したもの」に等しいというのである。 ところで、yi を yi + dyi に微小変化させれば、同時に i も i + d i に微小変化する。 ラグランジュ方程式というのは yi と i の 2 つの偏微分が含まれていてややこしく見えるかも知れないが、yi が微小変化すれば i はその変化に合わせて付いて来るのであり、要するに yi の微小変化さえ考えればいいのである。 つまり、どの yi を微小変化させても、この等式は成り立っているべきだという宣言になっているのだ。
では yi(t) を x 軸に沿った連続的なもの、すなわち y ( x , t ) へと拡張してやったらどうなるだろうか。 y ( x , t ) のある一部が微小変化することが、ある質点の位置が微小変化することに対応している。 その微小変化によってラグランジュ方程式の第 1 項と第 2 項で計算される量にはそれぞれ変化が生じるのだろうが、それらの変化量は互いに打ち消し合っていないといけない。 そのことが、y ( x , t ) のあらゆる地点をどんな微小変化させた場合においても成り立っていなければならないのである。
最速降下線問題とちょっと似た話になってきただろう。 つまり、関数 y をわずかにずらして y + δy に変更した時、 L にもわずかな変化 δL が生じるだろうが、δy がどんな形になっていようとも、第 1 項と第 2 項で計算される変化量は打ち消し合っていなくてはならないのである。 そのことを次のような式で表現することにする。
全くラグランジュ方程式と同じ形のままだ！ 第 1 項の時間微分が偏微分で表されていたり、ドットで表されていたものが偏微分になっていたりはするが、これは y が多変数関数になったからには従わなくてはならない表記法の理由であり、意味はあまり変わってはいない。 その他は、∂ が δ に変わっただけの違いである。
さあ、こうしてめでたく、ラグランジュ方程式に相当するものを式で表現する事だけは出来た。 しかしこの式を使って具体的に計算しろと言われても、どう手を付けたらいいのか、...
