=======================
角速度と速度との関係
=======================
この記事は、質点の角速度と速度との関係を幾何学的に説明するものです。2.に計算の説明を載せてありますが、
これは補助的なものです。内容がどうしても理解できない場合は、2.の方から読むのも良いかと思います。
1.角速度と速度の関係について
-----------------------------------
.. image:: okome-angularVelocity1.png
質点が、時刻 $t$ に 位置 ${\bm{x}(t)}$ にあるとします。またこの説明では ${\bm{x}(t)}$ は、座標3の ${z_{3}}$ 軸
上にあるものとしています(図1参照)。2.の方に詳しく書いていますが、座標3とは座標1を ${z_{1}}$ 軸周りに $\phi$ だけ回転し、
次に ${y_{1}}$ 軸周りに $\theta$ だけ回転させたものです。このときの位置 ${\bm{x}(t)}$ を座標3で表すと
\bm{x}_{3}(t)
&=
\begin{pmatrix}
x_{3}\\
y_{3}\\
z_{3}
\end{pmatrix}
=
\begin{pmatrix}
0\\
0\\
r
\end{pmatrix} \tag{1}
になります。 $\bm{x}_{3}$ とは $\bm{x}$ の座標3での成分表示のことを意味します。以下 ${\bm{x}_{n} (n= 1,2,3,4,5)}$ についても同じ意味で使います。
なお、速度と角速度についても同じ規約に従って書いるので、そのつもりで読んでいただきたいと思います。
次に座標3の原点から見たときの速度から、動径成分を取り除いたベクトルを ${\bm{v}_{3(rot)}}$ と書くとします。
すると ${\bm{v}_{3(rot)}}$ と速度 ${\bm{v}_{3}}$ の関係は次の式で書かれます。
&\bm{v}_{3}
=
\begin{pmatrix}
r\dot{\theta} \\
r\dot{\phi} \sin \theta \\
\dot{r}
\end{pmatrix}
=
\begin{pmatrix}
0 \\
0 \\
\dot{r}
\end{pmatrix}
+ \bm{v}_{3(rot)}
この式から ${\bm{v}_{3(rot)}}$ の成分が
\bm{v}_{3(rot)}
&=
r\begin{pmatrix}
\dot{\theta} \\
\dot{\phi} \sin \theta \\
0
\end{pmatrix} \tag{2}
になる事が直ちに分かります。これは座標3の原点周りの回転を表す速度ベクトルで、この値は座標原点の位置に依存する量です。
しかし今回の説明に限っては、全ての座標番号で原点は共通です。
.. image:: okome-angularVelocity2.png
ところで,ここでは角速度と速度を同時に説明するので、座標3のような軸の取り方が便利です.
図2を見てもらうと、 ${\bm{v}_{3(rot)}}$ と角速度の関係はすぐに分かると思います。その具体的な関係とは、 ${\bm{v}_{3(rot)}}$ の ${x_{3}}$ 成分
が ${y_{3}}$ 軸周りの回転、 ${y_{3}}$ 成分が ${-x_{3}}$ 軸回りの回転を表す成分になっているというものです。 ${z_{3}}$ 軸回りの回転は自転を意味するので、
大きさの無い質点の場合は定義することができません。以上のことから角速度の定義に従えば、角速度 $\bm{\omega_{3}}$ は次のようになることが分かります。
\bm{\omega}_{3} &=
\begin{pmatrix}
-\dot{\phi} \sin \theta\\
\dot{\theta} \\
0
\end{pmatrix} \tag{3}
ここで座標1について、少し補足しておきます。座標1とは時刻 $t$ における質点の位置 ${\bm{x}(t)}$ を次のような成分で表示する座標
のことで、数式で書くと次のようになります。
x_{1} &= r \sin\theta \cos\phi \tag{4}\\
y_{1} &= r \sin\theta \sin\phi \tag{5}\\
z_{1} &= r \cos\theta \tag{6}
図1を見れば分かる事ですが、これはまさに任意の位置を直交座標成分で書き表したものです。
次に座標1上での角速度ベクトル ${\bm{\omega}_{1}}$ も書いておきます。
\bm{\omega}_{1} =
\begin{pmatrix}
-\dot{\theta} \sin \phi - \dot{\phi} \sin \theta \cos \theta \cos \phi \\
\dot{\theta} \cos \phi - \dot{\phi} \sin \theta \cos \theta \sin \phi \\
\dot{\phi} \sin^{2} \theta
\end{pmatrix} \tag{7}
角速度の性格を眺めるには $\tag{3}$ の成分表示が役に立ちます。一方で、ある座標で任意の位置にある角速度の直交成分
を知るには $\tag{7}$ が必要です。これがここで説明したかったことの全てです。
理解できた方はさっそく下に載せてある練習問題を解いてみましょう。
2.速度を求めよう
------------------
.. image:: okome-angularVelocity3.png
ここでは具体的な計算を示しておきます。 ${\bm{v}_{3(rot)}}$ を求める方法は次の方法が考えられます。
計算の方針
^^^^^^^^^^^^
.. _A.:
A. 時刻 $t$ の質点の位置 ${\bm{x} (t)}$ が ${z_{3}}$ 軸上にのるような位置ベクトル ${\bm{x}_{3}}$ を、極座標変数 ${r,\theta,\phi}$ を用いて表していきます。
計算結果として、座標3上では動径方向以外の成分はゼロになることが分かるはずです。つまり、数式で書くと
\bm{x}_{3}
&=
\begin{pmatrix}
0\\
0\\
r
\end{pmatrix}
ということです。極座標変数で表すのは、動径成分と回転成分を分離して考えるのを容易にする狙いがあります。
.. _B.:
B. 時刻 ${t - \delta t}$ の質点の位置 ${\bm{x} ( t - \delta t )}$ が ${z_{5}}$ 軸上に乗るような
座標5上での位置ベクトル ${\bm{x}_{5}}$ をもとめていきます。
それができたら今度は、 ${\bm{x}_{3}}$ と ${\bm{x}_{5}}$ の関係式を求めることも考えます。計算結果は次のとおりです。
\bm{x}_{3} &= \bm{x}_{5}
-\begin{pmatrix}
r \delta \theta\\
r \delta \phi\\
0
\end{pmatrix}
この関係式を求める狙いは、座標3上での $\delta t$ 時間の微小変位量を求める事にあります。
.. _C.:
C. 質点の座標3上での速度成分を微分の定義に従って求めていきます。そのときB.で求めた式を利用して、次の計算結果が得られます。
\bm{v}_{3} &=
\begin{pmatrix}
0 \\
0 \\
\dot{r}
\end{pmatrix}
+ \begin{pmatrix}
r \dot{\theta} \\
r \dot{\phi} \sin \theta \\
0
\end{pmatrix}
計算
^^^^^
ここから「計算の方針」に従って計算を進めていきます。
`A.`_
まず時刻 $t$ における質点の位置 ${\bm{x}(t)}$ を座標1の成分表示を再掲しておきます。
x_{1} &= r \sin\theta \cos\phi \\
y_{1} &= r \sin\theta \sin\phi \\
z_{1} &= r \cos\theta
次に時刻 $t$ における質点の位置 ${\bm{x}(t)}$ が、 ${z_{3}}$ 軸上になるような座標変換を実行します。この変換は ${z_{1}}$ 軸周りに $\phi$ だけ回転し、
次に ${y_{1}}$ 軸周りに $\theta$ だけ回転させた変換です(図3参照)。結果が ${z_{3}}$ 成分が $r$ でそれ以外はゼロになる事は
計算しなくても分かります。しかし本当に考えたとおりになっているのかを調べておきます。この変換は後の計算にも利用するので、
できたら読んでいただいた方が分かりやすいかと思います。具体的な計算方法は次のとおりです。
[計算方法]
(ア)z軸周りに $\phi$ だけ回転
\begin{pmatrix}
x_{2}\\
y_{2}\\
z_{2}
\end{pmatrix}
&=
\begin{pmatrix}
\cos \phi & \sin \phi & 0\\
-\sin \phi & \cos \phi & 0\\
0 & 0 & 1
\end{pmatrix}
\begin{pmatrix}
x_{1}\\
y_{1}\\
z_{1}
\end{pmatrix} \tag{8}
(イ) ${y_{2}}$ 軸周りに $\theta$ だけ回転
\begin{pmatrix}
x_{3}\\
y_{3}\\
z_{3}
\end{pmatrix}
&= \begin{pmatrix}
\cos \theta & -\sin \theta & 0 \\
0 & 1 & 0 \\
\sin \theta & 0 & \cos \theta
\end{pmatrix}
\begin{pmatrix}
x_{2}\\
y_{2}\\
z_{2}
\end{pmatrix} \\
&=\begin{pmatrix}
\cos \theta & -\sin \theta & 0 \\
0 & 1 & 0 \\
\sin \theta & 0 & \cos \theta
\end{pmatrix}
\begin{pmatrix}
\cos \phi & \sin \phi & 0\\
-\sin \phi & \cos \phi & 0\\
0 & 0 & 1
\end{pmatrix}
\begin{pmatrix}
x_{1}\\
y_{1}\\
z_{1}
\end{pmatrix}\\
&=\begin{pmatrix}
\cos \theta \cos \phi & \cos \theta \sin \phi & -\sin \theta \\
-\sin \phi & \cos \phi & 0 \\
\sin \theta \cos \phi & \sin \theta \sin \phi & \cos \theta
\end{pmatrix}
\begin{pmatrix}
x_{1}\\
y_{1}\\
z_{1}
\end{pmatrix} \tag{9}
(ウ) $\tag{9}$ 式に $\tag{4}$ , $\tag{5}$ , $\tag{6}$ 式を代入
\begin{pmatrix}
x_{3}\\
y_{3}\\
z_{3}
\end{pmatrix}
&=
\begin{pmatrix}
0\\
0\\
r
\end{pmatrix}
以上の(ア),(イ),(ウ)の計算を実行すれば予測された結果である $\tag{1}$ 式が得られます。
これで座標3は、 ${z_{3}}$ 軸上に質点が時刻 $t$ にあるような座標になることが確かめられました。
ここでもう一つ気づかなければならない事があります。それは天下り的に与えた式 $\tag{4}$ , $\tag{5}$ , $\tag{6}$ 式は
(ア),(イ),(ウ)の変換を逆にたどっていけば得られるということです。
`B.`_
次に時刻 ${t - \delta t}$ に質点が、 ${z_{5}}$ 軸上に乗るような座標5上での位置ベクトル ${\bm{x}_{5}}$ を極座標変数を用いて書き表しておきます。
この座標変換は位置ベクトル ${\bm{x}_{1}}$ を $z_{1}$ を中心に ${\phi - \delta \phi}$ , ${y_{4}}$ 軸を中心に ${\theta - \delta \theta}$ だけ回転させた変換です。
ここに書いている ${\delta \phi}$ , ${\delta \theta}$ は微小な大きさの角です。
この変換は $\tag{9}$ 式を $\theta$ → $\theta$ $-$ $\delta$ $\theta$ , $\phi$ → $\phi$ $-$ $\delta$ $\phi$ , $\bm{x_{3}}$ → $\bm{x_{5}}$ と
書き換えた式を計算する事によって求まります。
\begin{pmatrix}
x_{5}\\
y_{5}\\
z_{5}
\end{pmatrix}
&=
\begin{pmatrix}
\cos (\theta - \delta \theta) \cos (\phi - \delta \phi) & \cos (\theta - \delta \theta) \sin (\phi - \delta \phi) & -\sin (\theta - \delta \theta) \\
-\sin (\phi - \delta \phi) & \cos (\phi - \delta \phi) & 0 \\
\sin (\theta - \delta \theta) \cos (\phi - \delta \phi) & \sin (\theta - \delta \theta) \sin (\phi - \delta \phi) & \cos (\theta + \delta \theta)
\end{pmatrix}
\begin{pmatrix}
x_{1}\\
y_{1}\\
z_{1}
\end{pmatrix}\tag{10}
ここで三角関数の
\sin (-\delta \psi) \simeq -\delta \psi,\cos (-\delta \psi) \simeq 1 (\delta \psi << 1) \tag{11}
という近似を用いると
\cos (\psi - \delta \psi) &= \cos \psi \cos (-\delta \psi) + \sin \psi \sin (-\delta \psi)\\
&= \cos \psi - \delta \psi \sin \psi \tag{12}\\
\sin (\psi - \delta \psi) &= \sin \psi \cos (-\delta \psi) - \sin (-\delta \psi) \cos \psi\\
&= \sin \psi + \delta \psi \cos \psi \tag{13}
が得られます。ここで近似式を用いておりますが、微分を行うときに、切り捨てた微小量は消えてしまいます。だから計算結果に影響する事はありません。 ${\tag{10}}$ 式に $\tag{12}$ , $\tag{13}$ に
それぞれ ${\theta - \delta \theta}$ , ${\phi - \delta \phi}$ を代入したものを放り込む事によって
\left(
\begin{array}{c}
x_{5}\\
y_{5}\\
z_{5}\\
\end{array}
\right)
&=
\left(
\begin{array}{ccc}
\cos \theta \cos \phi & \cos \theta\sin \phi & -\sin \theta \\
-\sin \phi & \cos \phi & 0 \\
\sin \theta \cos \phi & \sin \theta \sin \phi & \cos \theta \\
\end{array}
\right)
\left(
\begin{array}{c}
x_{1}\\
y_{1}\\
z_{1}\\
\end{array}
\right)
\\
&
+r\delta \theta
\left(
\begin{array}{ccc}
\sin \theta \cos \phi & \sin \theta \sin \phi & \cos \theta \\
0 & 0 & 0 \\
-\sin \theta \cos \phi & -\cos \theta \sin \phi & \sin \theta \\
\end{array}
\right)
\left(
\begin{array}{c}
x_{1}\\
y_{1}\\
z_{1}\\
\end{array}
\right)
\\
&+r\delta \phi
\left(
\begin{array}{ccc}
\cos \theta \sin \phi & -\cos \theta \cos \phi & 0 \\
\cos \phi & \sin \phi & 0 \\
\sin \theta \sin \phi & -\sin \theta \cos \phi & 0 \\
\end{array}
\right)
\left(
\begin{array}{c}
x_{1}\\
y_{1}\\
z_{1}\\
\end{array}
\right)
になることが分かります。更に $\tag{1}$ , $\tag{2}$ , $\tag{3}$ を代入すると
\begin{pmatrix}
x_{5}\\
y_{5}\\
z_{5}
\end{pmatrix}
&= \bm{x}_{3}
+\begin{pmatrix}
r \delta \theta\\
r \delta \phi\\
0
\end{pmatrix}\\
_{\cdot}{^{\cdot}}_{\cdot} \bm{x}_{3} &= \bm{x}_{5}
-\begin{pmatrix}
r \delta \theta\\
r \delta \phi\\
0
\end{pmatrix} \tag{14}
こうして ${\bm{x}_{3}}$ と ${\bm{x}_{5}}$ の関係を得ることができました。
`C.`_
時刻 $t$ に質点は ${z_{3}}$ 上あり、時刻 ${t- \delta t}$ に ${z_{5}}$ 上にあります。
この位置の変位を起こす間に ${\delta t}$ だけの時間が経過いるので、
座標3から見たときの速度 ${\bm{v}_{3}}$ は微分の定義式から計算する事ができます。 $\tag{1}$ 式
と $\tag{14}$ 式を微分の定義式に代入してやると
\bm{v}_{3} &= \lim_{\delta t \to 0} \frac{{\bm{x}_{3}}(t) - \bm{x}_{3} (t - \delta t )}{\delta t}
&=
\begin{pmatrix}
0 \\
0 \\
\dot{r}
\end{pmatrix}
+ \begin{pmatrix}
r \dot{\theta} \\
r \dot{\phi} \sin \theta \\
0
\end{pmatrix} \tag{15}
という計算結果が得られます。時刻 $t$ に質点は ${z_{3}}$ 軸方向上にあるので右辺第一項は動径方向を表しています。従って動径方向成分を除いた成分が、
座標原点から見たときの回転を表す速度成分になるわけです。座標3上での ${\bm{v}_{3(rot)}}$ を極座標変数を用いると、 ${\tag{2}}$ 式に書かれたとおりの結果に
なります。この速度 ${\bm{v}_{3(rot)}}$ は、動径成分においてのみ速さ $\dot{r}$ で運動する座標系から見た場合、その座標系での質点の速度になります。
読み返してもらえば分かりますが、この回転速度 ${\bm{v}_{3(rot)}}$ から角速度 ${\bm{\omega}_{3}}$ は直ちに分かります。
練習問題
-----------
練習問題は、かなりの分量があります。しかし今回の説明で、これだけの事が分かるのだと思うこともできると思います。
丁寧に取り組んでもらえれば、説明だけでは得られないものを感じる事ができるかもしれません。
練習問題1.
^^^^^^^^^^^
ここの説明で得られた ${\bm{\omega}_{3}}$ から、座標1での角速度 ${\bm{\omega}_{1}}$ を、座標変換の方法で求めてください。
これはよく角運動量 ${\bm{l}_{1}}$ の定義式と同時に与えられる公式
\bm{l}_{1} \equiv \bm{x}_{1} \times \bm{p}_{1} = mr^2 \bm{\omega}_{1}
で登場する角速度ベクトルと一致するはずです。この事実も同時に確認して下さい。
ヒントと答え( $\tag{7}$ 式)は以下に記しておきます。
(答え)
\bm{\omega}_{1} &=
\begin{pmatrix}
\omega_{x_{1}}\\
\omega_{y_{1}}\\
\omega_{z_{1}}
\end{pmatrix}
&=
\begin{pmatrix}
-\dot{\theta} \sin \phi - \dot{\phi} \sin \theta \cos \theta \cos \phi \\
\dot{\theta} \cos \phi - \dot{\phi} \sin \theta \cos \theta \sin \phi \\
\dot{\phi} \sin^{2} \theta
\end{pmatrix}
(ヒント)
オイラー角を利用した位置ベクトル ${\bm{x}_{3}}$ から ${\bm{x}_{1}}$ への変換と同じ操作をすれば求まります。
練習問題2.
^^^^^^^^^^^
座標1と3で、次の関係式を満たしている事を代入する事によって確かめてください。
\bm{v} = \bm{\omega} \times \bm{x} + \dot{r} \frac{\bm{x}}{r} \tag{16}
この関係式は速度と角速度を結ぶ一般的な関係式です。
練習問題3.
^^^^^^^^^^^
前問で確かめた式 $\tag{16}$ を ベクトルの回転_ の無限小回転を参考に成り立つ事を示してください。これではあまりにも不親切なので、ヒントも書いておきます。
(ヒント)
「ベクトルの回転」の"無限小回転"の $\tag{3}$ の法線ベクトル $\bm{n}$ は、今回の場合、常に角速度ベクトルと同じ向きを向く単位ベクトル
に等しくなります。つまり次のように書くことができるということです。
\bm{n} = \frac{\bm{\omega}}{|\bm{\omega}|} \tag{17}
すると位置ベクトルが動径方向に移動しないとし、 $\bm{n}$ を回転軸とする回転角を $\psi$ と書くとします。
すると"無限小回転"の $\tag{3}$ より、位置の微小変位は次のように書けます。
d \bm{x} = (\bm{\omega} \times \bm{x}) \frac{d \psi }{|\bm{\omega}|} \tag{18}
更に動径方向へも移動する一般の場合を考えます。すると位置の微小変化は動径方向への微小変位も加わり、次のようになります。
d \bm{x} &= (\bm{\omega} \times \bm{x}) \frac{d \psi}{|\bm{\omega}|} + dr \frac{\bm{x}}{r} \tag{19}
次にこれが微小時間 $dt$ の間に運動がなされたとしますと
\dot{\bm{x}} = (\bm{\omega} \times \bm{x}) \frac{\dot{\psi}}{|\bm{\omega}|} + \dot{r} \frac{\bm{x}}{r} \tag{20}
が成り立つことを示されます。最後に $\tag{20}$ 中の ${\dot{\psi}}$ が角速度の大きさ ${|\bm{\omega}|}$ に等しいことを示してください。
つまり次の式を示せばいいのです。
|\bm{\omega}| = \dot{\psi} \tag{21}
これは本文中の角速度 $\bm{\omega}$ の定義と ${\psi}$ の定義を比較することによって説明する事ができると思います。数式を使う必要はありません。
ここまでできますと、 $\tag{20}$ に $\tag{21}$ を代入することによって $\tag{16}$ 式が成り立つ事が示されます。以上の手順に従って計算をすれば、
この問題は解けると思います。
.. _ベクトルの回転 : http://hooktail.org/wiki/index.php?%BA%BA%C6%C9%2F%A5%D9%A5%AF%A5%C8%A5%EB%A4%CE%B2%F3%C5%BE%28Joh%C3%F8%29
練習問題4.
^^^^^^^^^^^^
座標3上での位置ベクトル、速度ベクトルをベクトル解析の方法で求めてください。実はこの方法の方が実践的です。
ここで、極座標 ${r, \theta , \phi}$ の単位ベクトルを与えておきます(この式の導出方法を忘れた方は後で確認しておく事を勧めます)。
&\bm{e}_{r} = \sin \theta \cos \theta \bm{e}_{x_{1}} + \sin \theta \sin \phi \bm{e}_{y_{1}} + cos \theta \bm{e}_{z_{1}} \tag{22}\\
&\bm{e}_{\theta} = \cos \theta \cos \theta \bm{e}_{x_{1}} + \cos \theta \sin \phi \bm{e}_{y_{1}} + cos \theta \bm{e}_{z_{1}} \tag{23}\\
&\bm{e}_{\phi} = = - \sin \theta \bm{e}_{x_{1}} + \cos \phi \bm{e}_{y_{1}} \tag{24}
ここで戸惑っている方もいるかもしれないので補足しておきますと、実は座標3は極座標(球面座標)と本質的には同じものです。
この問題はそれを確認するためのものです。最後に結果がセクション2.の `A.`_ での説明、つまり行列による計算の方法と一致する事を確かめておいてください。
練習問題5.
^^^^^^^^^^^
角速度をベクトル解析の方法で求めてください。これがいわゆる角速度ベクトルの定義式にあたるものです。これだけでは、不親切なので
またまたヒントを書いておきます。
(ヒント)
この問題は $\tag{16}$ とベクトルの公式を使えば求まります。具体的には $\bm{x}$ と $\tag{16}$ の両辺の外積を取ります。実際に計算してみると
次のようになります。
\bm{x} \times \bm{v} &= \bm{x} \times (\bm{\omega} \times \bm{x} + \dot{r} \frac{\bm{x}}{r})\\
&= \bm{x} \times (\bm{\omega} \times \bm{x}) + \frac{\dot{r}}{r} {\bm{x} \times \bm{x}}\\
&= \bm{\omega}(\bm{x} \cdot \bm{x}) - \bm{x} (\bm{x} \cdot \bm{\omega})\\
&= \bm{\omega} r^2
上の計算の3段目で、 $\bm{\omega}$ と $\bm{x}$ が垂直である事を用いています。このことはセクション1の説明を見れば、分かると思います。
すると結局、答えは次のとおりになります。
\bm{\omega} = \frac{\bm{x} \times \bm{v}}{r^2} \tag{25}
以上の手順をたどって答えを自分で求めてみてください。余談ですが、 $\tag{25}$ の両辺に慣性モーメント ${mr^{2}}$ をかけることによって角運動量が得られます。
慣性モーメントについては、また別の機会に説明をしたいと思います。少し長かったですがここまでたどりつけば、角速度についての一通りの理解は得られたと言えると思います。
お疲れ様でした。
@@author: ???@@
@@accept: ???@@