極坐標系 下的 牛頓第二定律


在 一個 極坐標系 里,  可以這樣 來 描述 牛頓第二定律, 如圖,  圖中 是一個 極坐標系,  質點 從 A 點 運動到 C 點,   極角 的 變化量 就是 ⊿θ  ,   極徑 的 變化量  ⊿ρ = BC 。    其中  OA = OB,  AH 垂直於 OB,  AB 間 的 弧 記為  弧AB,   AB 間 的 線段 仍然 稱為 AB    。

 

 

可以這樣 來 表達 牛頓第二定律 :

d²ρ / dt²  =  Fρ / m                  (1) 式

ρ  d²θ / dt²  =  Fθ / m              (2) 式

 

m  是 質點 的 質量  ,   Fρ  是 極徑方向 的 力,  Fθ 是 切線方向 的 力  。

 

我們說一下   (2) 式   的 推導過程 :

(1) 式   是  質點 在 極徑 ρ  方向 的 加速度,  簡稱 徑向加速度,       (2) 式 是  切線方向 的 加速度,  簡稱 切向加速度,

切線方向 就是 與  極徑 ρ  正交 的 方向   。

切向運動 總 和 徑向運動 正交,    切向位移 總是 和 徑向位移 正交,    質點 從 A 運動 到 B,  ρ 不變,   θ 變, 那么 切線方向 的 路程  s切 = 弧AB,  因為 弧AB 上每個 點 都 和  極徑 ρ  正交  。

如果 質點 從 A 運動 到 C,  ρ 變,   θ 變,  則 弧 AB 只有在 A 點 才和 ρ 正交,  弧 AB  不是 切線方向 的 路程  s切    。

此時,  切線方向 的 路程  s切  需要  通過 積分 求得 :

 

s切  =   ∫  ρ ( t )  dθ                (3) 式

 

ρ = ρ ( t )   的 時刻 t 的 ρ  ,    dθ  是  時刻 t 的 dθ       。

 

對 (3) 式 兩邊 微分,

ds切  =  ρ ( t )  dθ

兩邊除以   dt ,

ds切 / dt  =  ρ ( t )  dθ / dt

ds切 / dt   就是 切向速度 v切,  即

v切  =  ρ ( t )  dθ / dt

兩邊 對 dt 求導,

dv切 / dt = d ( ρ ( t )  dθ / dt ) / dt

dv切 / dt   就是 切向加速度 a切,  即

a切  =  d ( ρ ( t )  dθ / dt ) / dt

因為 現在是 在 計算 切向變量,  所以  徑向變量 可以看作 常量 (這是 神馬 邏輯 ?) , 所以  ρ ( t )  可以 提到 微分符號 外面 來,  於是,

a切  =  ρ ( t )  *  d ( dθ / dt ) / dt

a切  =  ρ ( t )  *  d²θ / dt²

因為 ρ = ρ ( t )  ,  ρ ( t )  可以 寫成  ρ ,  於是,

a切  =  ρ  *  d²θ / dt²

a切  =  ρ  d²θ / dt²

根據 牛頓第二定律    a切  =  Fθ / m    ,  於是,

Fθ / m  =  ρ  d²θ / dt²

ρ  d²θ / dt²  =  Fθ / m            此即  (2) 式

 

以上 推導 稱為 推導 1   。     

 

其實 還可以用 微元法 或者說 微分 的 方法 來 推導,   稱為 推導 2 :

 

設 質點 從 A 運動 到 C,  在 時刻 t 時 位於 A  。  雖然 質點 從 A 運動 到 C,  ρ 在 變化,   但是 當 ⊿θ -> 0 時,   ⊿ρ -> 0,    弧AB -> 0 ,    此時 的 弧AB 可以 認為 和  OA  正交,  OA 即 時刻 t 時 的   ρ  ,

所以 可以 用    弧AB -> 0   來 表示 時刻 t 時 的 切向路程 s切,  於是,  t 時 的 切向速度 v切 = s切 / ⊿t = 弧AB / ⊿t  ,    弧AB -> 0  ,   ⊿t -> 0     。

又因為       弧AB -> 0 = OA * ⊿θ = ρ * ⊿θ   ,      ⊿θ -> 0    ,於是,

v切 = ρ * ⊿θ / ⊿t      ,      ⊿θ -> 0  ,    ⊿t -> 0

寫成 微分 的 形式 就是 

v切 = ρ * dθ / dt          

 

注意,  推導 1  中 的 “現在是 在 計算 切向變量,  所以  徑向變量 可以看作 常量”   這個 原則 在 這里 仍然 適用,   因為 對於 dt,  ρ 是 有 dρ 的,   如果 在 計算 切向 速度 加速度 的 時候 再把  dρ  計算進來, 那就  ……  沒完沒了 了 。

 

從這里, 可以看出,   數學家 很多時候 也是 依靠 直觀 直覺 的 ,   只不過 這些 不會 寫在 論文 傳記 教科書 里  。  

 

不然的話,   請給出 為什么 此處 計算 v切 a切 不算入 dρ ,或者說 把 ρ 當作 常量 的 證明 。

我估計 對 這個 問題 較真 會 引發 第四次數學危機   ……  ~

 

書歸正傳,   對 v切 求導 就是 加速度 a切,   即

a切 = ( v切 ) ′ = d ( ρ * dθ / dt ) / dt  =  ρ *  d ( dθ / dt ) / dt  =  ρ  d²θ / dt²    ,    當然, 這里把 ρ 從 微分符號 里 提出來, 又用了一次   “現在是 在 計算 切向變量,  所以  徑向變量 可以看作 常量”   原則   。

於是,

a切 =  ρ  d²θ / dt²                此亦  (2) 式

 

大家可能會說   推導 1  其實 也是 以 推導 2 為 前提 的 ,  這樣的話,  推導 1 是 無意義 的 循環論證   ……    隨便,  你們 愛怎么想 就 怎么想 吧 ,哈哈 。

 

為什么 要 研究 極坐標系 下的 牛頓第二定律 ?    一個 直接的 需求 是 天體力學 里的 一體 二體 問題 。

 

一體 問題 就是 一個 行星 圍繞 一個 恆星 公轉,  恆星 不受 行星 的 引力 影響  。  恆星 不受 行星 的 引力 影響 是 一種 理想狀況, 所以 一體問題 又稱為 理想公轉問題   。

二體 問題 和 一體 問題 的 區別 是 兩個 天體 均受到 對方 引力 的 影響    。

見 《人造衛星軌道 和 天體軌道 原理》 https://www.cnblogs.com/KSongKing/p/11867972.html     。

在 《人造衛星軌道 和 天體軌道 原理》  中 提出了  一體 問題 在 直角坐標系 下 的  微分方程組, 由 x 坐標 和 y 坐標 的 2 個 方程 組成,

因為 兩個 方程 之間 互相關聯,  所以 方程組 難於求解   。

但是 用 極坐標系 的 話,    以 恆星 為 原點,    行星 受到 的 引力 F 永遠 在 極徑 方向,  這樣 行星 的 運動規律 表達為 牛頓第二定律 只需要 一個 方程 :

 

d²ρ / dt²  =  -  F / m         ,         m 為 行星質量

 

根據 萬有引力公式,   F = G M m / ρ²    ,    M 為 恆星質量  ,  代入 方程,

 

d²ρ / dt²  =  -  G M / ρ²                 (4) 式

 

這就是 極坐標系 下 的 一體問題 的 微分方程,   在 極坐標系 是 牛頓第二定律 用 一個 方程 就可以 表達 出來,   一個 方程 比 方程組 易於 求解   。

 

因為 行星 可能有 初速度,  初速度 可能 包含 徑向 和 切向 分量,  切向速度 會 讓 θ 發生變化,  當   θ  改變 時,   行星速度 的 徑向分量 和 切向分量  會 隨之改變,  

在 沒有 力 作用 的 情況下,   當   θ  改變 時,   行星速度 的 徑向分量 和 切向分量  會 隨之改變,  

也就是說,  徑向速度 的 加速度 除了 由 徑向 的 引力 F 引起 外,  還 包含 由 θ 變化 引起 的 成分,   θ 變化 的 原因 是 行星 的 速度 包含 切向分量 ,  行星 的 速度 包含 切向分量 的 原因 是 行星 的 初速度 包含 切向分量 ,  因為 如果 初速度 沒有 切向分量,   則  行星 只會 在 徑向 上 運動,  不會 產生 切向速度   。

 

所以,  還需要在 (4) 式 中 加入 θ 變化 引起 的 徑向加速度    。

 

設  行星速度 為 v,  徑向分量 為 v徑,   切向分量 為 v切, θ 變化 引起 的 徑向加速度 為  a_ρθ   ,

 

設  θ 變化 引起 的 徑向速度 變化 為 ⊿ v_θ引起  ,   微分 為  dv徑_θ引起 ,     則 :

 

dv徑_θ引起  =  v切 因 θ 變化 產生的 徑向增量 + v徑 因 θ 變化 產生的 徑向增量

 

v切 因 θ 變化 產生的 徑向增量  =  v切 sin dθ

v徑 因 θ 變化 產生的 徑向增量  =  v徑 cos dθ - v徑

 

dv徑_θ引起 = v切 sin dθ  +  v徑 cos dθ - v徑 

 

因為

v切 = ρ dθ / dt

v徑 = dρ / dt

 

所以,

dv徑_θ引起 = ρ dθ / dt sin dθ  +  dρ / dt cos dθ - dρ / dt

兩邊 除以 dt,

dv徑_θ引起 / dt = ( ρ dθ / dt sin dθ  +  dρ / dt cos dθ - dρ / dt ) / dt,

dv徑_θ引起 / dt  就是   θ 變化 引起 的 徑向加速度  a_ρθ   ,

即   a_ρθ = ( ρ dθ / dt sin dθ  +  dρ / dt cos dθ - dρ / dt ) / dt     。

 

行星 在 徑向 總的 加速度 aρ 等於  引力 引起 的 徑向加速度 a_引 加上  θ 變化 引起 的 徑向加速度  a_ρθ   。   即 :

 

aρ = d²ρ / dt² 

a_引 = - G M / ρ²

a_ρθ =  ( ρ dθ / dt sin dθ  +  dρ / dt cos dθ - dρ / dt ) / dt 

 

aρ = a_引 + a_ρθ

 

d²ρ / dt²  =  -  G M / ρ²  +  ( ρ dθ / dt sin dθ  +  dρ / dt cos dθ - dρ / dt  ) / dt       (5) 式

 

(5) 式 就是 在 (4) 式 中 加上了 a_ρθ  ,      (5) 式 就是  一體問題 在 極坐標系 下 的 微分方程 。

 

這個 微分方程  怎么解 ?    不知道   。  好處 是 把  ρ , θ   放到了 一個 方程 里,  避免了 方程組   。

這個 方程 的 解 是 ρ , θ 和 t 的 函數關系, 記為  f ( ρ , θ ) = t     ,    以 t 為 自變量 的 函數圖像 應該是一個 橢圓  。

即   一體問題 行星 公轉 的 軌道 是一個 橢圓,   橢圓形狀 和 位置 由 恆星質量 、行星質量 、行星 初始位置 初始速度 決定   。

 

解微分方程 的 過程中 會 進行 積分,   積分 會 產生 積分常數,    行星 初始位置 初始速度 就 體現 在 積分常數,   積分常數 的 具體數值 由 行星 初始位置 初始速度 決定    。

 

二體問題 可以通過 約化質量 轉化為 一體問題,  所以 也可以用 極坐標系 求解,  但 這樣得到的 解 仍然 是 不完備 的,  因為 這樣的解 只描述了 一個 質點 相對於 另一個 質點 的 運動狀況,   沒有 描述  兩個質點 在 第三方參照系 中 的 運動狀況 。

 

 

 

 

 

 

 

 

 

 

 

 

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM