正交曲线坐标系

                     

贡献者: addis; 零穹

  • 本文存在未完成的内容。
预备知识 柱坐标系,球坐标与直角坐标的转换

  1如果 $u, v, w$ 是三维空间中某曲线坐标系的三个坐标,空间任意一点的位置矢量 $ \boldsymbol{\mathbf{r}} $ 都是它们的函数 $ \boldsymbol{\mathbf{r}} (u, v, w)$。那么定义任意一点处三个单位矢量为

\begin{equation} \hat{\boldsymbol{\mathbf{u}}} = \frac{ \partial \boldsymbol{\mathbf{r}} /\partial u }{ \left\lvert \partial \boldsymbol{\mathbf{r}} /\partial u \right\rvert }\qquad \hat{\boldsymbol{\mathbf{v}}} = \frac{ \partial \boldsymbol{\mathbf{r}} /\partial v }{ \left\lvert \partial \boldsymbol{\mathbf{r}} /\partial v \right\rvert }\qquad \hat{\boldsymbol{\mathbf{w}}} = \frac{ \partial \boldsymbol{\mathbf{r}} /\partial w }{ \left\lvert \partial \boldsymbol{\mathbf{r}} /\partial w \right\rvert }~. \end{equation}
注意一般来说,这三个矢量会随着 $ \boldsymbol{\mathbf{r}} $ 改变。形象地说:当我们分别只把 $u, v, w$ 增加一点时,$ \boldsymbol{\mathbf{r}} $ 会分别沿 $ \hat{\boldsymbol{\mathbf{u}}} , \hat{\boldsymbol{\mathbf{v}}} , \hat{\boldsymbol{\mathbf{w}}} $ 方向移动(请以球坐标系和柱坐标系为例思考)。

   若对空间中任意一点,式 1 中的三个矢量都两两正交,那么这个曲线坐标系就是正交曲线坐标系(orthogonal curvilinear coordinate system)。常见的例子除了球坐标系柱坐标系 还有抛物线坐标系椭圆坐标系。直角坐标系也可以看作是正交曲线坐标系的一个特例。

例 1 球坐标系

   以下我们用式 1 重新推导球坐标系中的三个单位矢量(式 3 ),并证明它是一个正交曲线坐标系。

   球坐标系中,位置矢量可以表示为

\begin{equation} \boldsymbol{\mathbf{r}} = r \hat{\boldsymbol{\mathbf{r}}} = r\sin\theta\cos\phi\, \hat{\boldsymbol{\mathbf{x}}} + r\sin\theta\sin\phi\, \hat{\boldsymbol{\mathbf{y}}} + r\cos\theta \hat{\boldsymbol{\mathbf{z}}} ~. \end{equation}
同样,球坐标系的三个单位矢量由三个坐标增加的方向确定
\begin{equation} \left\{\begin{aligned} \frac{\partial \boldsymbol{\mathbf{r}} }{\partial r} &= \sin\theta\cos\phi\, \hat{\boldsymbol{\mathbf{x}}} + \sin\theta\sin\phi\, \hat{\boldsymbol{\mathbf{y}}} + \cos\theta\, \hat{\boldsymbol{\mathbf{z}}} \\ \frac{\partial \boldsymbol{\mathbf{r}} }{\partial \theta} &= r\cos\theta\cos\phi\, \hat{\boldsymbol{\mathbf{x}}} + r\cos\theta\sin\phi\, \hat{\boldsymbol{\mathbf{y}}} - r\sin\theta\, \hat{\boldsymbol{\mathbf{z}}} \\ \frac{\partial \boldsymbol{\mathbf{r}} }{\partial \phi} &= -r\sin\theta\sin\phi\, \hat{\boldsymbol{\mathbf{x}}} + r\sin\theta\cos\phi\, \hat{\boldsymbol{\mathbf{y}}} \end{aligned}\right. ~\end{equation}
归一化得三个单位矢量为
\begin{equation} \begin{cases} \hat{\boldsymbol{\mathbf{r}}} = \sin\theta\cos\phi\, \hat{\boldsymbol{\mathbf{x}}} + \sin\theta\sin\phi\, \hat{\boldsymbol{\mathbf{y}}} + \cos\theta\, \hat{\boldsymbol{\mathbf{z}}} \\ \hat{\boldsymbol{\mathbf{\theta}}} = \cos\theta\cos\phi\, \hat{\boldsymbol{\mathbf{x}}} + \cos\theta\sin\phi\, \hat{\boldsymbol{\mathbf{y}}} - \sin\theta\, \hat{\boldsymbol{\mathbf{z}}} \\ \hat{\boldsymbol{\mathbf{\phi}}} = -\sin\phi\, \hat{\boldsymbol{\mathbf{x}}} + \cos\phi\, \hat{\boldsymbol{\mathbf{y}}} \end{cases}~ \end{equation}
不难验证对除原点外的任意坐标点,这三个单位矢量两两间内积为零,即两两垂直,所以球坐标系属于正交曲线坐标系。

1. 矢量场

  

未完成:举例:$r \hat{\boldsymbol{\mathbf{r}}} $,$r \hat{\boldsymbol{\mathbf{\phi}}} $ 等。以及和直角坐标的转换。

2. 坐标线,坐标面

   正如开头所说,若空间任意一点 $M$ 的坐标可用三个新的变量 $u,v,w$ 来表示的话,则 $(u,v,w)$ 称为点 $M$ 的曲线坐标。给定新变量以常数值 $A,B,C$,就得到三族曲坐标面。将 $ \boldsymbol{\mathbf{r}} (u,v,w)$ 写成分量的形式

\begin{equation} x=\varphi_1(u,v,w),\quad y=\psi_1(u,v,w),\quad z=\omega_1(u,v,w)~. \end{equation}
或写成 $u,v,w$ 的形式
\begin{equation} u=\varphi(x,y,z),\quad v=\psi(x,y,z),\quad w=\omega(x,y,z)~. \end{equation}
则这些新的曲坐标面在坐标系 $x,y,z$ 中的方程为
\begin{equation} \varphi(x,y,z)=A,\quad \psi(x,y,z)=B,\quad \omega(x,y,z)=C~. \end{equation}

   由不同族中取出任意两个曲坐标面,它们相交于一线

\begin{equation} \psi(x,y,z)=B_0,\quad \omega(x,y,z)=C_0~. \end{equation}
其中,$B_0,C_0$ 为确定的常数,这条线只有 $u$ 变化,这样的线叫作坐标线 $u$。同理还可得到坐标线 $v,w$.

   我们知道,对于含参曲线 $ \boldsymbol{\mathbf{r}} (t)$,$ \mathrm{d}{ \boldsymbol{\mathbf{r}} }/\mathrm{d}{t} $ 是其切线方向向量,而 $( \mathrm{d}{ \boldsymbol{\mathbf{r}} }/\mathrm{d}{t} )/ \left\lvert \mathrm{d}{ \boldsymbol{\mathbf{r}} }/\mathrm{d}{t} \right\rvert $ 是其切线单位向量。而对于曲坐标线 $u,v,w$,$u,v,w$ 分别是其参数,所以曲坐标线 $u,v,w$ 的单位向量便是式 1 中的 $ \hat{\boldsymbol{\mathbf{u}}} , \hat{\boldsymbol{\mathbf{v}}} , \hat{\boldsymbol{\mathbf{w}}} $。所以这三个向量两两正交,便意味着三族曲坐标线的切线两两正交,所以称其为正交曲线坐标系。

   三条坐标线的切线之间两两正交,这意味着其中任一条坐标线的法向量必与另两个曲坐标在这点构成的曲坐标面相切,而任一曲坐标面上一点的法线必与过该点的另一条坐标线相切。

3. 线元

   在曲线坐标系中,线元平方为

\begin{equation} \begin{aligned} & \,\mathrm{d}{s} ^2= \,\mathrm{d}{x} ^2+ \,\mathrm{d}{y} ^2+ \,\mathrm{d}{z} ^2= \left( \frac{\partial \varphi_1}{\partial u} \,\mathrm{d}{u} + \frac{\partial \varphi_1}{\partial v} \,\mathrm{d}{v} + \frac{\partial \varphi_1}{\partial w} \,\mathrm{d}{w} \right) ^2+\\ & \left( \frac{\partial \psi_1}{\partial u} \,\mathrm{d}{u} + \frac{\partial \psi_1}{\partial v} \,\mathrm{d}{v} + \frac{\partial \psi_1}{\partial w} \,\mathrm{d}{w} \right) ^2+ \left( \frac{\partial \omega_1}{\partial u} \,\mathrm{d}{u} + \frac{\partial \omega_1}{\partial v} \,\mathrm{d}{v} + \frac{\partial \omega_1}{\partial w} \,\mathrm{d}{w} \right) ^2~. \end{aligned} \end{equation}
此式为 $ \,\mathrm{d}{u} , \,\mathrm{d}{v} , \,\mathrm{d}{w} $ 的二次齐次多项式,其含有乘积 $ \,\mathrm{d}{u} \,\mathrm{d}{v} $ 的项的系数为
\begin{equation} 2 \left( \frac{\partial \varphi_1}{\partial u} \cdot \frac{\partial \varphi_1}{\partial v} + \frac{\partial \psi_1}{\partial u} \cdot \frac{\partial \psi_1}{\partial v} + \frac{\partial \omega_1}{\partial u} \cdot \frac{\partial \omega_1}{\partial v} \right) ~. \end{equation}
同样的,含有 $ \,\mathrm{d}{u} \,\mathrm{d}{w} $ 和 $ \,\mathrm{d}{v} \,\mathrm{d}{w} $ 的项系数分别为
\begin{equation} \begin{aligned} & 2 \left( \frac{\partial \varphi_1}{\partial u} \cdot \frac{\partial \varphi_1}{\partial w} + \frac{\partial \psi_1}{\partial u} \cdot \frac{\partial \psi_1}{\partial w} + \frac{\partial \omega_1}{\partial u} \cdot \frac{\partial \omega_1}{\partial w} \right) ~,\\ &2 \left( \frac{\partial \varphi_1}{\partial v} \cdot \frac{\partial \varphi_1}{\partial w} + \frac{\partial \psi_1}{\partial v} \cdot \frac{\partial \psi_1}{\partial w} + \frac{\partial \omega_1}{\partial v} \cdot \frac{\partial \omega_1}{\partial w} \right) ~. \end{aligned} \end{equation}

   显然,这些项的系数是各坐标线切线方向向量的点乘,而在正交曲线坐标系中,我们已经知道它们为 0,即正交曲线坐标系中线元平方式 9 可写为

\begin{equation} \,\mathrm{d}{s} ^2=H_u^2 \,\mathrm{d}{u} ^2+H_v^2 \,\mathrm{d}{v} ^2+H_w^2 \,\mathrm{d}{w} ^2~. \end{equation}
其中
\begin{equation} \begin{aligned} &H_u^2= \left( \frac{\partial \varphi_1}{\partial u} \right) ^2+ \left( \frac{\partial \psi_1}{\partial u} \right) ^2+ \left( \frac{\partial \omega_1}{\partial u} \right) ^2= \left\lvert \frac{\partial \boldsymbol{\mathbf{r}} }{\partial u} \right\rvert ^2~,\\ &H_v^2= \left( \frac{\partial \varphi_1}{\partial v} \right) ^2+ \left( \frac{\partial \psi_1}{\partial v} \right) ^2+ \left( \frac{\partial \omega_1}{\partial v} \right) ^2= \left\lvert \frac{\partial \boldsymbol{\mathbf{r}} }{\partial v} \right\rvert ^2~,\\ &H_w^2= \left( \frac{\partial \varphi_1}{\partial w} \right) ^2+ \left( \frac{\partial \psi_1}{\partial w} \right) ^2+ \left( \frac{\partial \omega_1}{\partial w} \right) ^2= \left\lvert \frac{\partial \boldsymbol{\mathbf{r}} }{\partial w} \right\rvert ^2~. \end{aligned} \end{equation}

   由于每条坐标线上仅有对应的坐标变化,线元 $ \,\mathrm{d}{ \boldsymbol{\mathbf{s}} } $ 可写成

\begin{equation} \,\mathrm{d}{ \boldsymbol{\mathbf{s}} } =H_u \,\mathrm{d}{u} \hat{\boldsymbol{\mathbf{u}}} +H_v \,\mathrm{d}{v} \hat{\boldsymbol{\mathbf{v}}} +H_w \,\mathrm{d}{w} \hat{\boldsymbol{\mathbf{w}}} ~. \end{equation}
其中,$H_u,H_v,H_w$ 称为拉梅系数(Lamé coefficients)
\begin{equation} H_u= \left\lvert \frac{\partial \boldsymbol{\mathbf{r}} }{\partial u} \right\rvert ~,\qquad H_v= \left\lvert \frac{\partial \boldsymbol{\mathbf{r}} }{\partial v} \right\rvert ~,\qquad H_w= \left\lvert \frac{\partial \boldsymbol{\mathbf{r}} }{\partial w} \right\rvert ~. \end{equation}
或者,线元可直接求得
\begin{equation} \begin{aligned} \,\mathrm{d}{ \boldsymbol{\mathbf{s}} } &= \,\mathrm{d}{ \boldsymbol{\mathbf{r}} } = \frac{\partial \boldsymbol{\mathbf{r}} }{\partial u} \,\mathrm{d}{u} + \frac{\partial \boldsymbol{\mathbf{r}} }{\partial v} \,\mathrm{d}{v} + \frac{\partial \boldsymbol{\mathbf{r}} }{\partial w} \,\mathrm{d}{w} \\ &=H_u \,\mathrm{d}{u} \hat{\boldsymbol{\mathbf{u}}} +H_v \,\mathrm{d}{v} \hat{\boldsymbol{\mathbf{v}}} +H_w \,\mathrm{d}{w} \hat{\boldsymbol{\mathbf{w}}} ~. \end{aligned} \end{equation}
注意,这里 $ \partial \boldsymbol{\mathbf{r}} /\partial u , \partial \boldsymbol{\mathbf{r}} /\partial v , \partial \boldsymbol{\mathbf{r}} /\partial w $ 是各曲坐标线的切向向量,并且总是分别指向 $u,v,w$ 增大的方向。 令 $ \,\mathrm{d}{s} _u=H_u \,\mathrm{d}{u} , \,\mathrm{d}{s} _v=H_v \,\mathrm{d}{v} , \,\mathrm{d}{s} _w=H_w \,\mathrm{d}{w} $,则式 14 式 16 可写成更简洁的形式
\begin{equation} \,\mathrm{d}{ \boldsymbol{\mathbf{s}} } = \,\mathrm{d}{s} _u \hat{\boldsymbol{\mathbf{u}}} + \,\mathrm{d}{s} _v \hat{\boldsymbol{\mathbf{v}}} + \,\mathrm{d}{s} _w \hat{\boldsymbol{\mathbf{w}}} ~. \end{equation}

4. 面积元、体积元

   面积元和体积元可用于 “正交曲线坐标系中的重积分”。

   求得正交曲线坐标系中的线元式 16 后,其面积元 $ \,\mathrm{d}{ \boldsymbol{\mathbf{S}} } $ 和体积元 $ \,\mathrm{d}{V} $ 非常任意求得,即是

\begin{equation} \begin{aligned} \,\mathrm{d}{ \boldsymbol{\mathbf{S}} } &= \,\mathrm{d}{s} _v \,\mathrm{d}{s} _w \hat{\boldsymbol{\mathbf{u}}} + \,\mathrm{d}{s} _u \,\mathrm{d}{s} _w \hat{\boldsymbol{\mathbf{v}}} + \,\mathrm{d}{s} _u \,\mathrm{d}{s} _v \hat{\boldsymbol{\mathbf{w}}} \\ &=H_vH_w \,\mathrm{d}{v} \,\mathrm{d}{w} \hat{\boldsymbol{\mathbf{u}}} +H_uH_w \,\mathrm{d}{u} \,\mathrm{d}{w} \hat{\boldsymbol{\mathbf{v}}} +H_uH_v \,\mathrm{d}{u} \,\mathrm{d}{v} \hat{\boldsymbol{\mathbf{w}}} ~, \end{aligned} \end{equation}
\begin{equation} \,\mathrm{d}{V} = \,\mathrm{d}{s} _u \,\mathrm{d}{s} _v \,\mathrm{d}{s} _w=H_uH_vH_w \,\mathrm{d}{u} \,\mathrm{d}{v} \,\mathrm{d}{w} ~. \end{equation}


1. ^ 本文参考 Wikipedia 相关页面

                     

© 小时科技 保留一切权利