• 検索結果がありません。

x = r sin θ cos ϕ , y = r sin θ sin ϕ , z = r cos θ (1)

N/A
N/A
Protected

Academic year: 2021

シェア "x = r sin θ cos ϕ , y = r sin θ sin ϕ , z = r cos θ (1)"

Copied!
3
0
0

読み込み中.... (全文を見る)

全文

(1)

1

問題(「解析力学と相対論」から)

3次元極座標のニュートン方程式

ポテンシャルが

V

で与えられているときのニュートン方程式を極座標で表せ.ただし,極座 標への座標変換は

x = r sin θ cos ϕ , y = r sin θ sin ϕ , z = r cos θ (1)

で与えられる.

まず,極座標のパラメータへの変化方向の単位ベクトルを次のように定義する.

e

r

r r =

 sin θ cos ϕ sin θ sin ϕ

cos θ

, e

θ

= ∂e

r

∂θ =

 cos θ cos ϕ cos θ sin ϕ

sin θ

, e

ϕ

= ∂e

r

sin θ∂ϕ =

sin ϕ cos ϕ

0

 (2)

これらは正規直行基底をなす.基底のパラメータによる微分は

∂e

θ

∂θ = e

r

, ∂e

θ

∂ϕ = cos θe

ϕ

, ∂e

ϕ

∂θ = 0 , ∂e

ϕ

∂ϕ = sin θe

r

cos θe

θ

(3)

さらに,パラメータが時間によるとして基底の時間微分を計算すると

e ˙

r

= ˙ θ ∂e

r

∂θ + ˙ ϕ ∂e

r

∂ϕ = ˙ θe

θ

+ ˙ ϕ sin θe

ϕ

(4)

˙

e

θ

= ˙ θ ∂e

θ

∂θ + ˙ ϕ ∂e

ϕ

∂ϕ = θe ˙

r

+ ˙ ϕ cos θe

ϕ

(5) e ˙

ϕ

= ˙ ϕ ∂e

ϕ

∂ϕ = ϕ(sin ˙ θe

r

+ cos θe

θ

) (6)

これらを使うと,極座標での位置ベクトルは

x = re

r

(7)

速度ベクトルは

˙ x = d

dt (re

r

) = ˙ re

r

+ r e ˙

r

= ˙ re

r

+ r θe ˙

θ

+ r ϕ ˙ sin θe

ϕ

(8)

加速度ベクトルは

x ¨ = d

dt ( ˙ x) = d

dt ( ˙ re

r

+ r θe ˙

θ

+ r ϕ ˙ sin θe

ϕ

)

= re ¨

r

+ ˙ r( ˙ θe

θ

+ ˙ ϕ sin θe

ϕ

)

+ ˙ r θe ˙

θ

+ r θe ¨

θ

r θ ˙

2

e

r

+ r θ ˙ ϕ ˙ cos θe

ϕ

+ ˙ r ϕ ˙ sin θe

ϕ

+ r ϕ ¨ sin θe

ϕ

+ r ϕ ˙ θ ˙ cos θe

ϕ

r ϕ ˙

2

sin θ(sin θe

r

+ cos θe

θ

)

= (¨ r r θ ˙

2

r ϕ ˙

2

sin

2

θ)e

r

+ (r θ ¨ + 2 ˙ r θ ˙ r ϕ ˙

2

sin θ cos θ)e

θ

+ (r ϕ ¨ sin θ + 2 ˙ r ϕ ˙ sin θ + 2r θ ˙ ϕ ˙ cos θ)e

ϕ

(9)

一方ポテンシャル項を求めるためには,極座標による勾配

が必要になる.座標変換の式から

∂r = ∂x

∂r

∂x + ∂y

∂r

∂y + ∂z

∂r

∂z = ∂x

∂r · ∇ = e

r

· ∇ (10)

(2)

2

同様に

∂θ = ∂x

∂θ · ∇ = r ∂e

r

∂θ · ∇ = re

θ

· ∇ (11)

∂ϕ = ∂x

∂ϕ · ∇ = r ∂e

r

∂ϕ · ∇ = r sin θe

ϕ

· ∇ (12)

よって,

の極座標基底の方向への成分が

e

r

· ∇ =

∂r , e

θ

· ∇ =

r∂θ , e

ϕ

· ∇ =

r sin θ∂ϕ (13)

とわかる.つまり,

= e

r

∂r + e

θ

r∂θ + e

ϕ

r sin θ∂ϕ (14)

である.よって,ニュートン方程式

m x ¨ = −∇ V (15)

を,極座標の基底の成分に分解すると極座標表示のニュートン方程式が求まる.

m(¨ r r θ ˙

2

r ϕ ˙

2

sin

2

θ) = ∂V

∂r m(r θ ¨ + 2 ˙ r θ ˙ r ϕ ˙

2

sin θ cos θ) = ∂V r∂θ m(r ϕ ¨ sin θ + 2 ˙ r ϕ ˙ sin θ + 2r θ ˙ ϕ ˙ cos θ) = ∂V

r sin θ∂ϕ (16)

babababababababababababababababababababab

このように,ニュートン方程式を極座標で書くと非常に複雑になる.ラグランジアンを 使った方法ではこの方程式が直接得られるわけなので,ラグランジュに感謝するべきだろ う.2次元の場合と同様に

ϕ

方向には角運動量が現れている.実際方程式は

d

dt (m ϕr ˙

2

sin

2

θ) = ∂V

∂ϕ (17)

と書くことができる.これは,ポテンシャル

V

ϕ

によらないならば(

ϕ

が循環座標な らば)

L

z

= m ϕr ˙

2

sin

2

θ (18)

が保存することを意味する.

L

zは,

z

軸周りの角運動量になっている.

3次元極座標のラグランジアン

ポテンシャルが

V (r)

で与えられているときの粒子のラグランジアンを極座標で表せ.ただし,

極座標への座標変換は

x = r sin θ cos ϕ , y = r sin θ sin ϕ , z = r cos θ (19)

で与えられる.

速度ベクトルは

˙

x = ˙ re

r

+ r θe ˙

θ

+ r ϕ ˙ sin θe

ϕ

(20)

(3)

3

で与えられ,基底ベクトルの直交性を使うと

˙

x

2

= ˙ r

2

+ r

2

θ ˙

2

+ r

2

ϕ ˙

2

sin

2

θ (21)

よってラグランジアンは

L = 1

2 m( ˙ r

2

+ r

2

θ ˙

2

+ r

2

ϕ ˙

2

sin

2

θ) V (r) (22)

運動方程式は

d dt

∂L

r ˙ = r = ∂V

∂r (23)

d dt

∂L

θ ˙ = d

dt (mr

2

θ) ˙

= ∂L

∂θ = mr

2

ϕ ˙

2

sin θ cos θ (24) d

dt

∂L

ϕ ˙ = d

dt (mr

2

ϕ ˙ sin

2

θ)

= ∂L

∂θ = 0 (25)

L

2

= (mr

2

θ) ˙

2

+ (mr

2

ϕ ˙ sin θ)

2

= m

2

r

4

( ˙ θ

2

+ ˙ ϕ

2

sin

2

θ) (26)

L = 1

2 m r ˙

2

+ 1 2

L

2

r

2

V (r) (27)

L = 1 2 m

( u ˙ u

2

)

2

+ 1

2 L

2

u

2

αu (28)

mr

2

ϕ ˙ = L

z

(29)

r

2

= kdt , 1 k r ˙ = 1

r

2

r

= u

, (30)

˙

r = r

ϕ ˙ = kr

/r

2

= ku

(31) L = ( 1

2 mk

2

u

2

+ 1

2 L

2

u

2

αu)

u

2

(32)

参照

関連したドキュメント

The maximum likelihood estimates are much better than the moment estimates in terms of the bias when the relative difference between the two parameters is large and the sample size

There is a bijection between left cosets of S n in the affine group and certain types of partitions (see Bjorner and Brenti (1996) and Eriksson and Eriksson (1998)).. In B-B,

Abstract: In this paper, we proved a rigidity theorem of the Hodge metric for concave horizontal slices and a local rigidity theorem for the monodromy representation.. I

(The Elliott-Halberstam conjecture does allow one to take B = 2 in (1.39), and therefore leads to small improve- ments in Huxley’s results, which for r ≥ 2 are weaker than the result

“Breuil-M´ezard conjecture and modularity lifting for potentially semistable deformations after

S., Oxford Advanced Learner's Dictionary of Current English, Oxford University Press, Oxford

Keywords Catalyst, reactant, measure-valued branching, interactive branching, state-dependent branch- ing, two-dimensional process, absolute continuity, self-similarity,

At the end of the section, we will be in the position to present the main result of this work: a representation of the inverse of T under certain conditions on the H¨older