運動座標系のシステマティックな導出(1)

運動座標系への座標変換や，速度・加速度の記述について，線素および基底ベクトルを基本情報としてシステマティックに導出する手順を整理してみたい。 $\quad$

2次元極座標

目的とするところがよく見えるように，まず簡単な平面極座標について考えよう。デカルト座標および，極座標によるベクトルの表現をそれぞれの座標系の基底（単位ベクトル）を用いて，

$\boldsymbol{A} = A_x\boldsymbol{e}_x + A_y\boldsymbol{e}_y = A_r\boldsymbol{e}_r + A_\phi\boldsymbol{e}_\phi$

と書く。平面極座標は，座標変数を $(r,\theta)$ とすることも多いが，3次元極座標との関連を考えると方位角は $\phi$ を用いるのがよいと思う。

さて，まずベクトルの成分 $(A_x,A_y)$ と $(A_r,A_\phi)$ の間の変換だが，

$A_r = \boldsymbol{A}\cdot\boldsymbol{e}_r = A_x\boldsymbol{e}_x\cdot\boldsymbol{e}_r + A_y\boldsymbol{e}_y\cdot\boldsymbol{e}_r$

であるから，

$\left(\begin{matrix}A_r\\A_\phi\end{matrix}\right)=\left(\begin{matrix}\boldsymbol{e}_x\cdot\boldsymbol{e}_r\qquad\boldsymbol{e}_y\cdot\boldsymbol{e}_r\\ \boldsymbol{e}_x\cdot\boldsymbol{e}_\phi\qquad\boldsymbol{e}_y\cdot\boldsymbol{e}_\phi\end{matrix}\right)\left(\begin{matrix}A_x\\A_y\end{matrix}\right)$

と書ける。基底どうしの内積すなわち座標軸どうしの方向余弦を成分として，変換行列を書き下ろせばよい。しかし，この作業はさほど難しくはないがやや煩雑である。そこで，ベクトル成分の変換は，微小変位の変換に等しいことから，座標変換の微分をとってベクトル成分の変換を導出するのが簡便である。
f:id:yokkun831:20181026165038j:plain

座標変換

$x = r\cos\phi\quad,\qquad y = r\sin\phi$

の微分をとると，

$dx = dr\cos\phi - rd\phi\sin\phi \quad,\qquad dy = dr\sin\phi + rd\phi\cos\phi$

すなわち，

$\left(\begin{matrix}dx\\dy\end{matrix}\right)=\left(\begin{matrix}\cos\phi \quad -\sin\phi\\\sin\phi\qquad\cos\phi\end{matrix}\right)\left(\begin{matrix}dr\\rd\phi\end{matrix}\right)$

となるが，極座標の線素により

$d\boldsymbol{r} = dr\boldsymbol{e}_r + rd\phi\boldsymbol{e}_\phi$

と書けるから，微小変位の極座標成分は，

$(d\boldsymbol{r})_r = dr\quad,\qquad (d\boldsymbol{r})_\phi = rd\phi$

である。したがって，ベクトル成分の変換は

$\left(\begin{matrix}A_x\\A_y\end{matrix}\right)=\left(\begin{matrix}\cos\phi \quad -\sin\phi\\\sin\phi\qquad\cos\phi\end{matrix}\right)\left(\begin{matrix}A_r\\A_\phi\end{matrix}\right)$

となる。逆変換は，転置行列となるから

$\left(\begin{matrix}A_r\\A_\phi\end{matrix}\right)=\left(\begin{matrix}\quad\cos\phi \quad \sin\phi\\-\sin\phi\quad\cos\phi\end{matrix}\right)\left(\begin{matrix}A_x\\A_y\end{matrix}\right)$

である。デカルト座標と極座標との間の変換は，座標系の回転に他ならない。

次に速度，加速度を導出する。運動座標系では，基底（単位ベクトル）そのものが変化するから，その微小変化を考えるのが簡明である。方位角の微小変位 $d\phi$ にともなう基底の変化は，

$d\boldsymbol{e}_r = d\phi\boldsymbol{e}_\phi \quad,\quad d\boldsymbol{e}_\phi = -d\phi\boldsymbol{e}_r$

であるから，基底の時間変化率は

$\dot{\boldsymbol{e}_r}=\dot\phi\boldsymbol{e}_\phi \quad , \quad \dot{\boldsymbol{e}_\phi}=-\dot\phi\boldsymbol{e}_r$

となる。
f:id:yokkun831:20181026165141j:plain

なお，デカルト座標の基底にいったんもどして直接微分してもよい。基底の変換は，

$\left(\begin{matrix}\boldsymbol{e}_r\\\boldsymbol{e}_\phi\end{matrix}\right) =\left(\begin{matrix}\quad\cos\phi \quad \sin\phi\\ -\sin\phi\quad\cos\phi\end{matrix}\right)\left(\begin{matrix}\boldsymbol{e}_x\\\boldsymbol{e}_y\end{matrix}\right)$

これを微分して，

$\left(\begin{matrix}\boldsymbol{\dot{e}}_r\\\boldsymbol{\dot{e}}_\phi\end{matrix}\right) = \dot{\phi}\left(\begin{matrix}-\sin\phi \qquad \cos\phi\\ -\cos\phi\quad-\sin\phi\end{matrix}\right)\left(\begin{matrix}\boldsymbol{e}_x\\\boldsymbol{e}_y\end{matrix}\right)\\ 　　　 = \dot{\phi}\left(\begin{matrix}-\sin\phi \qquad \cos\phi\\ -\cos\phi\quad-\sin\phi\end{matrix}\right)\left(\begin{matrix}\cos\phi \quad -\sin\phi\\ \sin\phi\qquad\cos\phi\end{matrix}\right)\left(\begin{matrix}\boldsymbol{e}_r\\\boldsymbol{e}_\phi\end{matrix}\right)\\ 　　　 = \dot{\phi}\left(\begin{matrix}\,\,0\quad 1\\ -1 \quad 0\end{matrix}\right)\left(\begin{matrix}\boldsymbol{e}_r\\\boldsymbol{e}_\phi\end{matrix}\right)$

を得る。

$\boldsymbol{r}=r\boldsymbol{e}_r$ であるから，速度は

$\dot{\boldsymbol{r}} = \displaystyle\frac{d}{dt}(r\boldsymbol{e}_r) = \dot r \boldsymbol{e}_r + r\dot{\boldsymbol{e}_r} = \dot r \boldsymbol{e}_r + r\dot\phi\boldsymbol{e}_\phi$

したがって，加速度は

$\ddot{\boldsymbol{r}} = \displaystyle\frac{d}{dt}(\dot{\boldsymbol{r}}) = \ddot{r}\boldsymbol{e}_r + \dot{r}\dot{\boldsymbol{e}_r}+\dot{r}\dot{\phi}\boldsymbol{e}_\phi + r\ddot{\phi}\boldsymbol{e}_\phi+r\dot{\phi}\dot{\boldsymbol{e}_\phi} \\= \ddot{r}\boldsymbol{e}_r+2\dot{r}\dot{\phi}\boldsymbol{e}_\phi+r\ddot{\phi}\boldsymbol{e}_\phi-r\dot{\phi}^2\boldsymbol{e}_r\\ = (\ddot{r}-r\dot{\phi}^2)\boldsymbol{e}_r + \displaystyle\frac{1}{r}\frac{d}{dt}(r^2\dot{\phi})\boldsymbol{e}_\phi$

となる。

運動座標系のシステマティックな導出(2) - 科学のおもちゃ箱@Hatena

（初稿：2010/04/17）