3

力学変数の変換とEuler-Lagrange方程式

386
1

私は最近、「基幹講座 解析力学」という本で解析力学を学んでいるのですが、その途中で出てきた次の問いに苦戦していました。一昨日解決したので、思考の整理がてら記事を生やそうと思い筆を執っています。

さて、次のような問題です。

LagrangianL(q,q˙,t)で記述されるN自由度系において、力学変数qの代わりにN個のN+1変数関数fiを用いて、
qi(t)=fi(Q1(t),Q2(t),,QN,t)(i=1,2,,N)
で関係した新力学変数Qをとる。qのLagrangianL(q,q˙,t)を新力学変数Qで表した、
LQ(Q,Q˙,t)=L(f(Q,t),ddtf(Q,t),t)に対応したEuler-Lagrange方程式

LQQiddtLQQi˙=0
が、元のEuler-Lagrange方程式と等価であることを示せ。

問題文の主張は?

問の主張は、ざっくり説明すると、

例えば2次元のある系の運動に関しての記述で、力学変数として、直交座標(x,y)を用いて、Euler-Lagrange方程式をたてようと、極座標(r,θ)を用いて、Euler-Lagrange方程式をたてようと同じ結果が得られる。
これは次元が上がろうと、良くわからない頭のおかしい力学変数の取り方をしようと、同じ系の運動に関しての記述であれば結果は同じになる。

というものです。

一般化座標を用いての議論なのだから、座標の取り方によらず同一の結果が得られることはなんとなくわかりますが、定量的に議論していきましょう。

さて、いきなり多変数(高次元)で話をしても理解しにくいので、まずは2次元で議論しましょう。

例えば、2次元で(x,y)(r,θ)の変換を考えてみます。

x=rcosθなので、rx=1cosθであり、また、
dx=cosθdrrsinθdθ なので、両辺を微小時間で割って、
x˙=r˙cosθθ˙rsinθr˙x˙=1cosθ

よって、
LxddtLx˙=0LrrxddtLr˙r˙x˙=0Lr1cosθddtLr˙1cosθ=01cosθ(LrddtLr˙)=0

できましたね。
θについても全く同様に可能です。

ここで重要だったのが、rx=r˙x˙です。
これが成り立っていたから最後に1cosθでくくれた訳ですから。

だんだん答えが見えてきましたね。一般化していきましょう。

(q1,q2,,qN)(Q1,Q2,,QN)の変換。

これは、点変換(point transformation)と呼ばれる、Lagrangianについての変換です。

qi(t)=fi(Q1,Q2,,QN,t)である。
Qi˙qj˙=(dQidt)(dqjdt)=Qiqj=Qifjを踏まえて考えていく。

LqiddtLqi˙=0LQQiQiqjddtLQQi˙Qi˙qj˙=0Qiqj(LQQiddtLQQi˙)=0

よって、LQQiddtLQQi˙=0である。(証明終了)

いやあ疲れましたね。偏微分の記号が非常に打つのに時間がかかる。
解析力学を学び始めて5日目にしてこの問題にぶち当たった訳ですからそれはもう心が折れそうでした。
問題の解釈、それを当然の主張に見ること、そして多変数関数や偏微分の考え方など問題が山積みの状態だったので、理解までに10日もかけてしまいました。
ですがこの問題を解決したことで、解析力学の基本的なもののとらえ方の片鱗が理解できたので、時間をかけた価値は十分にあったと思います。
また、具体例で考えることの大切さも再認識しました。
拙い文章でしたが、最後まで読んでいただきありがとうございました。

2020/11/20

投稿日:20201119
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

はてな
はてな
23
1694
今は 解析力学、複素解析、Γ関数周り について学んでいます

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中