logo

座標系での速度と加速度 📂古典力学

座標系での速度と加速度

球座標系での速度と加速度

v=r˙r^+rθ˙θ^+rϕ˙sin⁡θϕ^a=(r¨−rθ˙2−rϕ˙2sin⁡2θ)r^+(rθ¨+2r˙θ˙−rϕ˙2sin⁡θcos⁡θ)θ^+(rϕ¨sin⁡θ+2r˙ϕ˙sin⁡θ+2rθ˙ϕ˙cos⁡θ)ϕ^ \begin{align*} \mathbf{v} &=\dot{r} \hat {\mathbf{r}} +r \dot{\theta} \hat{ \boldsymbol{\theta}}+ r \dot{\phi} \sin{\theta} \hat{ \boldsymbol{\phi}} \\ \mathbf{a} &= (\ddot{r}-r\dot\theta^2-r\dot\phi^2\sin^2\theta)\hat{\mathbf{r}}+(r\ddot\theta+2\dot{r}\dot\theta-r\dot\phi^2\sin\theta\cos\theta)\hat{\boldsymbol{\theta}} \\ &\quad+(r\ddot\phi\sin\theta+2\dot{r}\dot\phi\sin\theta+2r\dot\theta\dot\phi\cos\theta)\hat{\boldsymbol{\phi}} \end{align*}

導出

596F6AA20.png

球座標系での単位ベクトルは以下の通りだ。

r^=cos⁡ϕsin⁡θx^+sin⁡ϕsin⁡θy^+cos⁡θz^θ^=cos⁡ϕcos⁡θx^+sin⁡ϕcos⁡θy^−sin⁡θz^ϕ^=−sin⁡ϕx^+cos⁡ϕy^ \begin{align*} \hat{\mathbf{r}} &= \cos \phi \sin \theta \hat{\mathbf{x}} + \sin \phi \sin \theta \hat{\mathbf{y}} + \cos\theta\hat{\mathbf{z}} \\ \hat{\boldsymbol{\theta}} &= \cos\phi \cos\theta \hat{\mathbf{x}} + \sin\phi \cos\theta \hat{\mathbf{y}} - \sin\theta\hat{\mathbf{z}} \\ \hat{\boldsymbol{\phi}} &= -\sin\phi \hat{\mathbf{x}} + \cos\phi \hat{\mathbf{y}} \end{align*}

これから球座標系で速度と加速度を順に求めよう。速度は位置を時間で微分して、加速度は速度を時間で微分して求めることができる。ちなみに r˙\dot{r} は「アル ドット」と読む。物理では文字の上の点は時間に対する微分を意味する。

r˙=drdt \dot{r}=\frac{dr}{dt}

速度と加速度を求める前に、まず単位ベクトルの微分を先に求めておこう。r^˙\dot{\hat{ \mathbf{r}}}、θ^˙\dot {\hat{ \boldsymbol{\theta}}}、ϕ^˙\dot {\hat{ \boldsymbol{\phi}}}を計算するとそれぞれ以下のようになる。

r^˙=ddt(cos⁡ϕsin⁡θx^+sin⁡ϕsin⁡θy^+cos⁡θz^)=dcos⁡ϕdtsin⁡θx^+cos⁡ϕdsin⁡θdtx^+dsin⁡ϕdtsin⁡θy^+sin⁡ϕdsin⁡θdty^+dcos⁡θdtz^=dcos⁡ϕdϕdϕdtsin⁡θx^+cos⁡ϕdsin⁡θdθdθdtx^+dsin⁡ϕdϕdϕdtsin⁡θy^+sin⁡ϕdsin⁡θdθdθdty^+dcos⁡θdθdθdtz^=−ϕ˙sin⁡ϕsin⁡θx^+θ˙cos⁡ϕcos⁡θx^+ϕ˙cos⁡ϕsin⁡θy^+θ˙sin⁡ϕcos⁡θy^−θ˙sin⁡θz^=θ˙(cos⁡ϕcos⁡θx^+sin⁡ϕcos⁡θy^−sin⁡θz^)+ϕ˙sin⁡θ(−sin⁡ϕx^+cos⁡ϕy^)=θ˙θ^+ϕ˙sin⁡θϕ^ \begin{align*} \dot{\hat{ \mathbf{r}}} &= \frac{d}{dt}(\cos\phi\sin\theta\hat{\mathbf{x}} +\sin\phi\sin\theta\hat{\mathbf{y}}+\cos\theta\hat{\mathbf{z}} ) \\ &= \frac{d\cos\phi}{dt}\sin\theta\hat{\mathbf{x}} + \cos\phi\frac{d \sin\theta}{dt}\hat{\mathbf{x}}+\frac{d \sin\phi}{dt}\sin\theta\hat{\mathbf{y}} + \sin\phi\frac{d\sin\theta}{dt}\hat{\mathbf{y}}+\frac{d \cos\theta}{dt}\hat{\mathbf{z}} \\ &= \frac{d\cos\phi}{d\phi}\frac{d\phi}{dt}\sin\theta\hat{\mathbf{x}} + \cos\phi\frac{d \sin\theta}{d\theta}\frac{d\theta}{dt}\hat{\mathbf{x}} \\ &\quad +\frac{d \sin\phi}{d\phi}\frac{d\phi}{dt}\sin\theta\hat{\mathbf{y}} + \sin\phi\frac{d\sin\theta}{d\theta}\frac{d\theta}{dt}\hat{\mathbf{y}}+ \frac{d \cos\theta}{d\theta}\frac{d\theta}{dt}\hat{\mathbf{z}} \\ &= -\dot\phi\sin\phi \sin\theta\hat{\mathbf{x}} + \dot\theta\cos\phi\cos\theta\hat{\mathbf{x}}+\dot\phi\cos\phi \sin\theta\hat{\mathbf{y}} + \dot\theta\sin\phi\cos\theta\hat{\mathbf{y}}-\dot\theta\sin\theta\hat{\mathbf{z}} \\ &= \dot{\theta}(\cos\phi\cos\theta\hat{\mathbf{x}} + \sin\phi\cos\theta\hat{\mathbf{y}} -\sin\theta \hat{\mathbf{z}}) + \dot{\phi}\sin\theta (-\sin\phi\hat{\mathbf{x}}+\cos\phi\hat{\mathbf{y}}) \\ &= \dot{\theta} \hat{ \boldsymbol{\theta} } + \dot{\phi} \sin\theta \hat{ \boldsymbol{\phi}} \end{align*}

θ^˙=ddt(cos⁡ϕcos⁡θx^+sin⁡ϕcos⁡θy^−sin⁡θz^)=−ϕ˙sin⁡ϕcos⁡θx^−θ˙cos⁡ϕsin⁡θx^+ϕ˙cos⁡ϕcos⁡θy^−θ˙sin⁡ϕsin⁡θy^−θ˙cos⁡θz^=−θ˙(cos⁡ϕsin⁡θx^+sin⁡ϕsin⁡θy^+cos⁡θz^)+ϕ˙cos⁡θ(−sin⁡ϕx^+cos⁡ϕy^)=−θ˙r^+ϕ˙cos⁡θϕ^ \begin{align*} \dot {\hat{ \boldsymbol{\theta}}} &= \frac{d}{dt}(\cos\phi\cos\theta\hat{\mathbf{x}} + \sin\phi\cos\theta\hat{\mathbf{y}} -\sin\theta \hat{\mathbf{z}}) \\ &= -\dot{\phi}\sin\phi\cos\theta\hat{\mathbf{x}}-\dot\theta\cos\phi\sin\theta\hat{\mathbf{x}}+\dot{\phi} \cos\phi \cos\theta \hat{\mathbf{y}}-\dot{\theta} \sin\phi \sin\theta \hat{\mathbf{y}}-\dot\theta\cos\theta\hat{\mathbf{z}} \\ &= -\dot\theta (\cos\phi\sin\theta\hat{\mathbf{x}}+\sin\phi\sin\theta \hat{\mathbf{y}} + \cos\theta\hat{\mathbf{z}}) + \dot{\phi}\cos\theta (-\sin\phi\hat{\mathbf{x}}+\cos\phi\hat{\mathbf{y}}) \\ &= - \dot{\theta} \hat{ \mathbf{r}} + \dot{\phi} \cos \theta \hat {\boldsymbol{\phi}} \end{align*}

ϕ^˙=ddt(−sin⁡ϕx^+cos⁡ϕy^)=−ϕ˙cos⁡ϕx^−ϕ˙sin⁡ϕy^=−ϕ˙(cos⁡ϕx^+sin⁡ϕy^) \begin{equation} \begin{aligned} \dot {\hat{ \boldsymbol{\phi}}} &= \frac{d}{dt}(-\sin\phi\hat{\mathbf{x}}+\cos\phi\hat{\mathbf{y}}) \\ &= -\dot\phi\cos\phi\hat{\mathbf{x}}-\dot\phi\sin\phi\hat{\mathbf{y}} \\ &= -\dot\phi (\cos\phi\hat{\mathbf{x}} + \sin\phi\hat{\mathbf{y}}) \end{aligned} \end{equation}

ϕ^˙\dot {\hat{\boldsymbol{\phi}}}の結果は、他の成分と異なりr^\hat{\mathbf{r}}とθ^\hat{\boldsymbol{\theta}}に直接整理されない。よく見るとz^\hat{\mathbf{z}}の成分がないことが分かる。r^\hat{\mathbf{r}}にsin⁡θ\sin\thetaを掛けてθ^\hat{\boldsymbol{\theta}}にcos⁡θ\cos\thetaを掛けて両方を足すとz^\hat{\mathbf{z}}の成分がなくなる。これを利用しよう。

 sin⁡θr^+cos⁡θθ^=cos⁡ϕsin⁡2θx^+sin⁡ϕsin⁡2θy^+sin⁡θcos⁡θz^+cos⁡ϕcos⁡2θx^+sin⁡ϕcos⁡2θy^−sin⁡θcos⁡θz^=cos⁡ϕ(sin⁡2θ+cos⁡2θ)x^+sin⁡ϕ(sin⁡2θ+cos⁡2θ)y^=cos⁡ϕx^+sin⁡ϕy^ \begin{align*} \ \sin\theta\hat{\mathbf{r}}+\cos\theta \hat{\boldsymbol{\theta}} &= \cos\phi\sin^2\theta\hat{\mathbf{x}}+\sin\phi\sin^2\theta\hat{\mathbf{y}}+\sin\theta\cos\theta \hat{\mathbf{z}} \\ &\quad + \cos\phi\cos^2\theta\hat{\mathbf{x}} +\sin\phi\cos^2\theta\hat{\mathbf{y}} -\sin\theta\cos\theta \hat{\mathbf{z}} \\ &= \cos\phi (\sin^2\theta+\cos^2\theta)\hat{\mathbf{x}}+\sin\phi (\sin^2\theta+\cos^2\theta)\hat{\mathbf{y}} \\ &= \cos\phi\hat{\mathbf{x}}+\sin\phi\hat{\mathbf{y}} \end{align*}

上の式を(1)(1)に代入すると以下の通りだ。

ϕ^˙=−ϕ˙sin⁡θr^−ϕ˙cos⁡θθ^ \dot{\hat{\boldsymbol{\phi}}}=-\dot\phi\sin\theta\hat{\mathbf{r}} - \dot\phi\cos\theta \hat{\boldsymbol{\theta}}

速度

r\mathbf{r}をttで微分すると以下の通りだ。

v=drdt=ddt(rr^)=drdtr^+rdr^dt=r˙r^+rr^˙=r˙r^+rθ˙θ^+rϕ˙sin⁡θϕ^ \begin{align*} \mathbf{v} = \frac{d \mathbf{r}}{dt} = \frac{d}{dt}(r \hat{\mathbf{r}}) = \frac{d r}{dt}\hat{\mathbf{r}} + r\frac{d \hat{\mathbf{r}}}{dt} &= \dot{r} \hat{\mathbf{r}} + r \dot{\hat{\mathbf{r}}} \\ &= \dot{r} \hat{\mathbf{r}} + r\dot{\theta} \hat{ \boldsymbol{\theta} } + r\dot{\phi} \sin\theta \hat{ \boldsymbol{\phi}} \end{align*}

事前に計算しておいた単位ベクトルの微分を代入するだけで結果が得られる。

■

加速度

v\mathbf{v}をttで微分すると以下の通りだ。

a=dvdt=ddt(r˙r^+rθ˙θ^+rϕ˙sin⁡θϕ^)=(r¨r^+r˙r^˙)+(r˙θ˙θ^+rθ¨θ^+rθ˙θ^˙)+(r˙ϕ˙sin⁡θϕ^+rϕ¨sin⁡θϕ^+rϕ˙θ˙cos⁡θϕ^+rϕ˙sin⁡θϕ^˙) \begin{align*} \mathbf{a} = \frac{d \mathbf{v}}{dt} &= \frac{d}{dt}(\dot{r} \hat{\mathbf{r}} +r\dot\theta\hat{\boldsymbol{\theta}}+ r\dot\phi\sin\theta \hat{\boldsymbol{\phi}}) \\ &= (\ddot r \hat{\mathbf{r}} + \dot{r} \dot{\hat{\mathbf{r}}}) + (\dot{r}\dot\theta\hat{\boldsymbol{\theta}} + r\ddot\theta\hat{\boldsymbol{\theta}} + r\dot\theta\dot{\hat{\boldsymbol{\theta}}}) \\ &\quad + (\dot{r} \dot\phi\sin\theta \hat{\boldsymbol{\phi}}+ r\ddot\phi\sin\theta\hat{\boldsymbol{\phi}}+ r\dot\phi\dot\theta\cos\theta\hat{\boldsymbol{\phi}}+ r\dot\phi\sin\theta\dot{\hat{\boldsymbol{\phi}}}) \end{align*}

すごく長い。じっくり解いてみよう。単位ベクトルの微分は既に求めてあるので、代入して整理すれば終わりだ。

a=[r¨r^+r˙(θ˙θ^+ϕ˙sin⁡θϕ^)]+[r˙θ˙θ^+rθ¨θ^+rθ˙(−θ˙r^+ϕ˙cos⁡θϕ^)]+[r˙ϕ˙sin⁡θϕ^+rϕ¨sin⁡θϕ^+rϕ˙θ˙cos⁡θϕ^+rϕ˙sin⁡θ(−ϕ˙sin⁡θr^−ϕ˙cos⁡θθ^)]=(r¨−rθ˙2−rϕ˙2sin⁡2θ)r^+(rθ¨+2r˙θ˙−rϕ˙2sin⁡θcos⁡θ)θ^+(rϕ¨sin⁡θ+2r˙ϕ˙sin⁡θ+2rθ˙ϕ˙cos⁡θ)ϕ^ \begin{align*} \mathbf{a} &= \left[ \ddot r \hat{\mathbf{r}} + \dot{r} (\dot{\theta} \hat{\boldsymbol{\theta}} + \dot{\phi} \sin\theta \hat{\boldsymbol{\phi}}) \right] + \left[ \dot{r}\dot\theta\hat{\boldsymbol{\theta}}+ r\ddot\theta\hat{\boldsymbol{\theta}}+r\dot\theta ( -\dot{\theta} \hat{\mathbf{r}} + \dot{\phi} \cos \theta \hat{\boldsymbol{\phi}}) \right] \\ &\quad + \left[ \dot{r} \dot\phi\sin\theta\hat{\boldsymbol{\phi}}+ r\ddot\phi\sin\theta\hat{\boldsymbol{\phi}}+ r\dot\phi\dot\theta\cos\theta\hat{\boldsymbol{\phi}}+ r\dot\phi\sin\theta (-\dot\phi\sin\theta\hat{\mathbf{r}} - \dot\phi\cos\theta\hat{\boldsymbol{\theta}}) \right] \\ &= (\ddot{r}-r\dot\theta^2-r\dot\phi^2\sin^2\theta)\hat{\mathbf{r}}+(r\ddot\theta+2\dot{r}\dot\theta-r\dot\phi^2\sin\theta\cos\theta)\hat{\boldsymbol{\theta}} \\ &\quad +(r\ddot\phi\sin\theta+2\dot{r}\dot\phi\sin\theta+2r\dot\theta\dot\phi\cos\theta)\hat{\boldsymbol{\phi}} \end{align*}

■

関連項目