logo

円筒座標系における速度と加速度 📂古典力学

円筒座標系における速度と加速度

円柱座標系での速度と加速度

$$ \begin{align*} \mathbf{v}&=\dot{r} \hat{\mathbf{r}} + r \dot{\phi} \hat{\boldsymbol{\phi}}+\dot{z} \hat{\mathbf{z}} \\ \mathbf{a} &= (\ddot r -r\dot{\phi} ^2)\hat{\mathbf{r}} + (2\dot{r} \dot{\phi} + r\ddot{\phi})\hat{\boldsymbol{\phi}} + \ddot{z}\hat{\mathbf{z}} \end{align*} $$

導出

596DB6CE0.jpg

円柱座標系では単位ベクトルは以下の通りです。

$$ \begin{align*} \boldsymbol{\rho}&=x\hat{\mathbf{x}}+y \hat{\mathbf{y}} +z\hat{\mathbf{z}}=r\hat{\mathbf{r}} +z\hat{\mathbf{z}} \\ \hat{\mathbf{r}} &= \hat{\mathbf{r}}(\phi) = \cos\phi \hat{\mathbf{x}} + \sin\phi \hat{\mathbf{y}} \\ \hat{\boldsymbol{\phi}} &= \hat{\mathbf{r}}(\phi+\pi/2) = -\sin\phi \hat{\mathbf{x}} + \cos\phi \hat{\mathbf{y}} \\ \hat{\mathbf{z}} &= \hat{\mathbf{z}} \end{align*} $$ 速度は位置を時間に対して微分し、加速度は速度を時間に対して微分して求められます。参考までに、$\dot{r}$は「ドット」と読みます。物理学で文字の上に点は時間に対する微分を意味します。

$$ \dot{r} = \frac{dr}{dt} $$

速度

$\boldsymbol{\rho}$を$t$に対して微分すると、以下のようになります。

$$ \begin{align*} \mathbf{v}&=\frac{d \boldsymbol{\rho}}{dt} \\ &=\frac{d}{dt}(r \hat{\mathbf{r}}+z\hat{\mathbf{z}}) \\ &=\frac{d r}{dt}\hat{\mathbf{r}} + r\frac{d \hat{\mathbf{r}}}{dt} +\frac{dz}{dt}\hat{\mathbf{z}} +z\frac{d \hat{\mathbf{z}}}{dt} \\ &=\dot{r} \hat{\mathbf{r}} +r \dot{\hat {\mathbf{r}}}+\dot{z} \hat{\mathbf{z}} \end{align*} $$

この時点で、直交座標系の単位ベクトルは時間によって変わりません。つまり、以下のようです。

$$ \dot{\hat{\mathbf{x}}}=\dot{\hat{\mathbf{y}}}=\dot{\hat{\mathbf{z}}} = \mathbf{0} $$

$\dot{\hat{\mathbf{r}}}$を計算すると、以下のようになります。

$$ \begin{align*} \dot{\hat {\mathbf{r}}} = \frac{d}{dt}(\hat{\mathbf{r}}) &= \frac{d}{dt}(\cos\phi \hat{\mathbf{x}}) + \frac{d}{dt}(\sin\phi \hat{\mathbf{y}}) \\ &= \frac{d\cos\phi}{dt}\hat{\mathbf{x}} + \frac{d\sin\phi}{dt}\hat{\mathbf{y}} \\ &= \frac{d\cos\phi}{d \phi}\frac{d \phi}{dt}\hat{\mathbf{x}}+\frac{d\sin\phi}{d \phi}\frac{d \phi}{dt}\hat{\mathbf{y}} \\ &= -\sin\phi \frac{d \phi}{dt}\hat{\mathbf{x}}+\cos\phi \frac{d \phi}{dt}\hat{\mathbf{y}} \\ &= \frac{d \phi }{dt}(-\sin\phi \hat{\mathbf{x}}+\cos\phi \hat{\mathbf{y}}) \\ &= \dot{\phi} \hat{\boldsymbol{\phi}} \end{align*} $$

従って、速度は以下のようになります。

$$ \mathbf{v}=\dot{r} \hat{\mathbf{r}} + r \dot{\phi} \hat{\boldsymbol{\phi}}+\dot{z} \hat{\mathbf{z}} $$

加速度

$\mathbf{v}$を$t$に対して微分すると、以下のようになります。

$$ \begin{align*} \mathbf{a}=\frac{d \mathbf{v}}{dt} &= \frac{d}{dt}(\dot{r} \hat{\mathbf{r}} + r \dot{\phi} \hat{\boldsymbol{\phi}} +\dot{z} \hat{\mathbf{z}}) \\ &= \ddot r \hat{\mathbf{r}} +\dot{r} \dot{ \hat{\mathbf{r}}} + \dot{r} \dot{\phi} \hat{\boldsymbol{\phi}} + r \ddot{\phi} \hat{\boldsymbol{\phi}} + r \dot{\phi} \dot{ \hat{\boldsymbol{\phi}}} +\ddot{z}\hat{\mathbf{z}} \end{align*} $$

$\dot{\hat{\boldsymbol{\phi}}}$を計算してみると、以下のようになります。

$$ \begin{align*} \dot{ \hat{\boldsymbol{\phi}}} = \frac{d}{dt}(\hat{\boldsymbol{\phi}}) &= \frac{d}{dt}(-\sin\phi \hat{\mathbf{x}})+\frac{d}{dt}(\cos\phi \hat{\mathbf{y}}) \\ &= -\frac{d\sin\phi}{dt}\hat{\mathbf{x}} +\frac{d\cos\phi}{dt}\hat{\mathbf{y}} \\ &= -\frac{d\sin\phi}{d \phi}\frac{d \phi}{dt}\hat{\mathbf{x}}+\frac{d\cos\phi}{d \phi}\frac{d \phi}{dt}\hat{\mathbf{y}} \\ &= \dfrac{d\phi}{dt} (-\cos\phi \hat{\mathbf{x}}-\sin\phi \hat{\mathbf{y}}) \\ &= - \dot{\phi} \hat{\mathbf{r}} \end{align*} $$

従って、代入して整理すると、以下のようになります。

$$ \begin{align*} \mathbf{a} &= \ddot r \hat{\mathbf{r}} +\dot{r} \dot{ \hat{\mathbf{r}}} + \dot{r} \dot{\phi} \hat{\boldsymbol{\phi}} + r \ddot{\phi} \hat{\boldsymbol{\phi}} + r \dot{\phi} \dot{ \hat{\boldsymbol{\phi}}} +\ddot{z}\hat{\mathbf{z}} \\ &= \ddot r \hat{\mathbf{r}} +\dot{r} \dot{\phi}\hat{\boldsymbol{\phi}} + \dot{r} \dot{\phi} \hat{\boldsymbol{\phi}} + r \ddot{\phi} \hat{\boldsymbol{\phi}} -r \dot{\phi} \dot{\phi} \hat{\mathbf{r}} +\ddot{z}\hat{\mathbf{z}} \\ &= (\ddot r -r\dot{\phi} ^2)\hat{\mathbf{r}} + (2\dot{r} \dot{\phi} + r\ddot{\phi})\hat{\boldsymbol{\phi}} +\ddot{z}\hat{\mathbf{z}} \end{align*} $$


円柱座標系は極座標系に高さ$z$だけを加えたものであるため、速度と加速度の公式もそれぞれ極座標系では$\dot{z} \hat{\mathbf{z}}$、$\ddot{z}\hat{\mathbf{z}}$項を加えた形になります。

参照