2003 年高一力學課程大要 陳義裕
Lecture 6
本次演講是個高等課題,它會讓你們驚嘆力學理論的美妙!
回顧:
哈米頓原理:(Hamilton's principle)
給定一段固定的時間T,同時也把起始以及最後的狀態都固定住。則一個日常 力學系統中的各個粒子會自動選擇自己的軌跡,使得
成為一個極值!
換言之,你如果把真實軌跡做一些微小的改變而得到一些「測試軌跡」,則沿 著「測試軌跡」所算出來的積分值和沿著真實軌跡所算出來的積分值是一樣的。
{ }
0T ( dt
∫
總動能)-(總位能)以下是哈米頓原理的一點應用。
例:
在一度空間中有兩個粒子在交互作用,它們之間的交互作用位能V 只和它們間的 相對位置有關:V =V x( 1−x )2 。請利用這個系統所具有在空間中的平移對稱性來 證明系統的總動量是守恆的。
解:
把時間 分成很多 小段 :
T
∆t
T = ∆N t。在右圖中 我們以黑色曲線畫 出粒子真實的軌跡 ABCDE 及 LFGHK 隨時間 t 改變的情 形。
現在我們考慮兩條粒子的「測試軌跡」AB’C’D’E 及 LF’G’H’K。(真實的粒子當 然
不是
走這些軌跡!) 這兩條軌跡是把 BCD 段在空間中整個平移 成為 B’C’D’,同時也把 FGH 段在空間中整個平移 成為F’G’H’而得。δx δx
在哈米頓原理中我們必須考慮以下的積分:
2 2
1 1 2 2
1 2
0
2 2
( ) ( )
( ) ( )
1 1 2 2
1 2
( )
2 2
( )
2 2
T
j j
j j
j
m dx m dx
V x x dt
dt dt
m dx m dx
V x x t
dt dt
⎧ ⎫
⎪ ⎛ ⎞ + ⎛ ⎞ − − ⎪
⎨ ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ ⎬
⎪ ⎪
⎩ ⎭
⎧ ⎛ ⎞ ⎛ ⎞ ⎫
⎪ ⎪
= ⎨ ⎜ ⎟ + ⎜ ⎟ − − ⎬
⎝ ⎠ ⎝ ⎠
⎪ ⎪
⎩ ⎭
∫
∑
∆根據該原理,沿著「測試軌跡」所算出來的積分值和沿著真實軌跡所算出來的積 分值必須是一樣的。 然而很明顯的是:在真實軌跡( ,x x1 2)上的位能值
和「測試軌跡」對應位置
1 2
( ) V x −x
1 2
(x ',x ')上的位能值V x( 1'−x2')都是一樣的——這是因 為這個系統具有在空間上的平移對稱性之故——所以我們便發現:只需考慮以下 的積分便可:
2 2
1 1 2 2
0
2 2
( ) ( )
1 1 2 2
2 2
2 2
T
j j
j
m dx m dx
dt dt dt
m dx m dx
dt dt t
⎧ ⎫
⎪ ⎛ ⎞ + ⎛ ⎞ ⎪
⎨ ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ ⎬
⎪ ⎪
⎩ ⎭
⎧ ⎛ ⎞ ⎛ ⎞ ⎫
⎪ ⎪
= ⎨ ⎜ ⎟ + ⎜ ⎟ ⎬
⎝ ⎠ ⎝ ⎠
⎪ ⎪
⎩ ⎭
∫
∑
∆但是我們製作出來的測試軌跡還有另外一個特性就是:沿著 B’C’D’段上之測試 粒子速度與在BCD段上之粒子原始速度根本是一樣的;同理,沿著 F’G’H’段上 之測試粒子速度與在FGH段上之粒子原始速度也是一樣的。(這是因為大家都平 移了δ 距離之故。) 因此,哈米頓原理就告訴了我們: x
2 2 2 2
1 1 2 2 1 1 2 2
2 2 2 2
1 1 2 2 1 1 2 2
' ' '
2 2 2 2
2 2 2 2
AB LF DE HK
AB LF D E H K
m dx m dx m dx m dx
dt dt dt dt t
m dx m dx m dx m dx
dt dt dt dt t
⎧ ⎫
⎪ ⎛ ⎞ + ⎛ ⎞ + ⎛ ⎞ + ⎛ ⎞ ⎪
⎨ ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ ⎬
⎪ ⎪
⎩ ⎭
⎧ ⎫
⎪ ⎛ ⎞ ⎛ ⎞ ⎛ ⎞ ⎛ ⎞ ⎪
=⎨ ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ ⎬
⎝ ⎠ ⎝ ⎠ ⎝ ⎠ ⎝ ⎠
⎪ ⎪
⎩ ' ⎭
∆
∆
接著注意到:
1 1 '
1 1
'
2 2
'
2 2
'
AB AB
D E DE
LF LF
H K HK
dx dx x
dt dt t
dx dx x
dt dt t
dx dx x
dt dt t
dx dx x
dt dt t
δ δ δ δ
= +
∆
= −
∆
= +
∆
= −
∆
把這些全部代入上式,同時只保留到δ 的一次方項(因為x δ 很小之故),於是我x 們就得到
1 2 1 2
1 2 1 2
AB LF DE HK
dx dx dx dx
m m m m
dt dt dt dt
⎧ ⎛ ⎞ + ⎛ ⎞ ⎫ ⎧= ⎛ ⎞ + ⎛ ⎞
⎨ ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠ ⎬ ⎨ ⎜⎝ ⎟⎠ ⎜⎝ ⎟⎠
⎩ ⎭ ⎩
⎫⎬
⎭ 這道式子不偏不倚正是動量守恆定律!
例:
一個質量是一個單位的粒子在一度空間中受到彈簧的回復力在做運動。已知彈簧 的彈性位能是
2
( ) 2 V x = x
若粒子一開始的時候是位於原點,且希望它在時間t=2的時候恰好跑到 之 位置。那麼請利用哈米頓原理近似地算出粒子於任一時刻 之位置。(0 2)
2 x=
t ≤ ≤t
解:
如果我們根據牛頓運動定律來處理這個問題,則所列出之方程式就應該是
1 dV
dv F dx x
dt m v dx
dt
= = = − = −
= = 加速度 速度
(請注意第一道式子告訴我們粒子受到的力量遵守虎克定律。) 不過我們到目前 為止所學到的微積分知識仍不足以應付這兩道聯立方程式。唉! 所以我們才要 退而求其次,利用哈米頓原理來近似地算出粒子於任一時刻 t 之位置!
我們的猜測是:
( ) (2 ) x t ≈ +t at − t
(注意到這個猜測確實滿足我們要求的條件:粒子一開始的時候是位於原點,且 它在時間t=2的時候恰好跑到x=2之位置。) 其中 a 是個未定的係數。現在把
( )
x t 代入哈米頓原理的積分式中:
( )
{ }
2 2
0
2 2 2
0
2 2 4 3 2 2
0 2
1 ( ) 2
1 1
1 2 2 ( (2 ))
2 2
1 2 (1 2 ) (4 1) 4 (1 2 ) (1 2 ) 2
4 4 1 5 3 3
dx V x dt dt
a at t at t dt
a t a a t a t a a t a dt
a a
⎧ ⎫
⎪ ⎛ ⎞ − ⎪
⎨ ⎜⎝ ⎟⎠ ⎬
⎪ ⎪
⎩ ⎭
⎧ ⎫
= ⎨ + − − + − ⎬
⎩ ⎭
= − + + − + − + + +
= − −
∫
∫
∫
哈米頓原理說我們一定要讓這個積分值取極(小)值才會越接近真實的軌跡,所以 我們可以從
4 2 4 1 5 3 3 0
d a a
da
⎛ ⎞
− − =
⎜ ⎟
⎝ ⎠
求出來 5
a= 。所以 6
5 8 2
( ) (2 )
6 3
5 x t ≈ +t t − =t t−6t 就是一個近似解。事實上,這題的答案是
2 sin
( ) 2.2 sin sin 2
x t = t ≈ t
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4
x(t)
0.2 0.4 0.6 0.8 時間 t1 1.2 1.4 1.6 1.8 2
右方這個圖把這兩個答案畫 在一起供你比較。(黑色線是 精確解,藍色線是近似解) 你 大概會很驚訝說我們這麼粗 糙的近似竟然和真實的答案 如此接近吧!