經典力學——正交曲線坐標系

本文是自己學習經典力學做的筆記的一些整理,自娛自樂向。


1、曲線坐標

位置矢量: vec{r}=vec{r}(q_1,q_2,q_3)

坐標面:確保某一個坐標為常數,其上兩個坐標可變。即 q_1=C_1 q_2=C_2 q_3=C_3 等。

坐標線:坐標面兩兩相交得坐標線,其上只有一個坐標可變。

基矢:取基矢沿坐標線切線方向,正向為相應坐標增大方向。

單位基矢:模長歸一的基矢。數學定義如下:

e_alpha(vec{q})=frac{1}{h_alpha}frac{partialvec{r}}{partial q_alpha},alpha=1,2,3

其中 h_alpha 為度規係數或拉梅係數,起模長歸一作用: h_alpha(vec{q})=|frac{partialvec{r}}{partial q_alpha}|

線元矢量的定義為:

dvec{r}equivvec{r}(vec{q}+dvec{q})-vec{r}(vec{q})\ =sum_alphafrac{partialvec{r}}{partial q_alpha}dq_alpha\ =sum_alpha h_alpha dq_alpha vec{e}_alpha

兩點間距離定義為:

ds^2=dvec{r}cdot dvec{r} =sum_{alpha,eta}e_alpha e_eta h_alpha h_eta dq_alpha dq_eta

若坐標系滿足正交歸一條件:

vec{e}_alphacdotvec{e}_eta=delta_{alphaeta}

ds^2 可進一步化簡為:

ds^2=sum_alpha h_alpha^2dq_alpha^2

以柱坐標為例:

有位置矢量 vec{r}=
ho(cos{phi}vec{e_x}+sin{phi}vec{e_y})+zvec{e_z}

各坐標偏導數有

frac{partialvec{r}}{partial
ho}=cosphivec{e_x}+sinphivec{e_y}=vec{e_
ho},h_
ho=1\ frac{partialvec{r}}{partialphi}=
ho(-sinphivec{e_x}+cosphivec{e_y})=
hovec{e_phi},h_phi=
ho\ frac{partialvec{r}}{partial z}=vec{e_z},h_z=1

位矢可重新表示為 vec{r}=
hovec{e_
ho}+zvec{e_z}

線元矢量有 dvec{r}=d
hovec{e_
ho}+
ho dphivec{e_phi}+dzvec{e_z}

兩點間距離有 ds^2=d
ho^2+
ho^2dphi^2+dz^2

2、正交曲線坐標系中的速度與加速度

位置矢量 vec{r}=r(vec{q})

對於速度:

依定義,速度矢量 vec{v}=dot{vec{r}}=sum_alphafrac{partialvec{r}}{partial q_alpha}dot{q_alpha}=sum_alpha h_alphadot{q_alpha}vec{e_alpha}

速度的平方 v^2(vec{q},dot{vec{q}})=vec{v}cdotvec{v}=sum_alpha h_alpha^2dot{q_alpha}^2

對於加速度:

依定義,加速度矢量 vec{a}=ddot{vec{r}}=sum_alpha a_alphavec{e_alpha}

加速度的各個分量滿足 a_alpha=vec{a}cdotvec{e_alpha}=frac{1}{h_alpha}left[ frac{d}{dt}frac{partial}{partialdot{q_alpha}}left( frac{v^2}{2}
ight)-frac{partial}{partial{q_alpha}}left( frac{v^2}{2}
ight) 
ight]...(1)

下面證明式 (1)

先證引理1(可以把偏微分中上下的點同時約去或添加): frac{partialdot{vec{r}}}{partial dot{q_alpha}} =frac{partialvec{r}}{partial q_eta}frac{partialdot{q_eta}}{partial dot{q_alpha}} =frac{partialvec{r}}{partial q_eta}delta_{alphaeta}=frac{partialvec{r}}{partial q_alpha}

再證引理2(對坐標的偏微分和對時間的微分可以互換順序):

frac{d}{dt}frac{partialvec{r}}{partial q_alpha}=frac{partial}{partial q_eta}frac{partialvec{r}}{partial q_alpha}dot{q_eta} =frac{partial}{partial q_alpha}frac{partialvec{r}}{partial q_eta}dot{q_eta} =frac{partialdot{vec{r}}}{partial q_alpha}

最後證明原命題:

h_alpha a_alpha=h_alpha vec{e_alpha}cdotvec{a} =frac{partialvec{r}}{partial q_alpha}cdotddot{vec{r}} =frac{d}{dt}left( dot{vec{r}}cdotfrac{partialvec{r}}{partial q_alpha} 
ight)-dot{vec{r}}cdotfrac{d}{dt}frac{partialvec{r}}{partial q_alpha}\

對第一項逆用引理1,對第二項使用引理2:

h_alpha a_alpha=frac{d}{dt}left( dot{vec{r}}cdotfrac{partialdot{vec{r}}}{partialdot{ q_alpha}} 
ight)-dot{vec{r}}cdotfrac{partialdot{vec{r}}}{partial q_alpha}\ =frac{d}{dt}frac{partialfrac{1}{2}dot{vec{r}}^2}{partial q_alpha}-frac{d}{dt}frac{partialfrac{1}{2}dot{vec{r}}^2}{partial q_alpha}

a_alpha=frac{1}{h_alpha}left[ frac{d}{dt}frac{partial}{partialdot{q_alpha}}left( frac{v^2}{2}
ight)-frac{partial}{partial{q_alpha}}left( frac{v^2}{2}
ight) 
ight]

證畢。

3、自然坐標

對於軌道已知的運動:

弧坐標:沿軌道自參考點 O_s 至質點位置的有向弧段的數值。

位矢 vec{r}=vec{r}(s)

速度 vec{v}=frac{dvec{r}}{dt}=frac{dvec{r}}{ds}frac{ds}{dt}

定義切向單位矢 vec{	au}=frac{dvec{r}}{ds}

其模長的平方 vec{	au}cdotvec{	au}=left| frac{dvec{r}}{ds} 
ight|^2=1

速度可改寫為 vec{v}=frac{ds}{dt}vec{	au}

由於切向單位矢模長為常數,即:

frac{d(vec{	au}cdotvec{	au})}{dt}=2frac{dvec{	au}}{dt}cdotvec{	au}=0

可見 frac{dvec{	au}}{dt}otvec{	au} ,於是定義 left|frac{dvec{	au}}{dt}
ight|=left|frac{d	heta}{ds}
ight|=frac{1}{R} 為曲率, R 為曲率半徑。

可以看出上述矢量模長為 frac{1}{R} ,故可以定義主法向單位矢 vec{n}=Rfrac{dvec{	au}}{dt}

從而定義出副法向單位矢 vec{b}=vec{	au}	imesvec{n}

定義切向單位矢與主法向單位矢張成的平面為密切面;

定義切向單位矢與副法向單位矢張成的平面為從切面;

定義主副法向單位矢張成的平面為法面。

同時也可以定義切線、主法線、副法線。

此時速度可以表示為 vec{v}=vvec{	au}

加速度可以表示為 vec{a}=dot{vec{v}}=dot{v}vec{	au}+vdot{vec{	au}}=dot{v}vec{	au}+vdot{s}frac{dvec{	au}}{ds}=dot{v}vec{	au}+frac{v^2}{R}vec{n}

從中可以看出加速度始終位於密切面內,且 a_	au=dot{v}a_n=frac{v^2}{R}


最後舉一道例題:

例:一質點沿拋物線軌道 y=frac{x^2}{2p}A(-p,frac{p}{2})B(p,frac{p}{2}) 運動,任意時刻其切、法向加速度為一常數 C ,求質點在 AB 處速率之比。


解答:

加速度滿足:

a_	au=dot{v},a_n=frac{v^2}{R}

其比值:

frac{a_	au}{a_n}=frac{dv}{v}frac{R}{vdt}=frac{dv}{v}frac{R}{ds}=frac{dv}{v}frac{1}{d	heta}=C...(1)

有幾何關係:

	an{	heta}=frac{dy}{dx}=frac{x}{p}

(1) 式分離變數並積分: int_{v_A}^{v_B}frac{dv}{v}=int_{	heta_A}^{	heta_B}Cd	heta

即:

lnleft(frac{v_B}{v_A}
ight)=Carctanfrac{x}{p}{huge|}_{-p}^{p}=Cfrac{pi}{2}

於是

frac{v_B}{v_A}=e^{Cpi/2}

推薦閱讀:

能分析一下《全民炸金花》的勝率問題嗎?
《通天狄仁傑》中狄仁傑斷案用的博弈術是什麼?
輪椅上的不朽 ——紀念霍金先生
振動熱力學-19. Phonon-Phonon Anharmonicity
ABC KIDS《童趣大冒險》千年禁地中的聲音密碼究竟何方神聖?

TAG:物理學 | 力學 | 經典力學 |