=======================
角速度1
=======================
角速度はある座標の原点から見たときの単位動径距離あたりの速度の回転成分と同じ大きさを持つ軸性ベクトルです。
つまり次元は時間分の1( $SI$ 単位系 $[$ $sec^{-1}$ $]$ )で速度の次元をもつわけではありません。また軸性ベクトル
は基本的に普通のベクトルとは違います。ここではまずオイラー角を用いて速度の回転成分を出すことについてのみ説明します。
角速度については次の項で説明します。次のところでもオイラー角を使って説明するのでここで少し回転操作に慣れておくのが良いと思います。
読むときは途中計算が省略されている部分が多々あるので手を動かす必要があるかと思います。
また図も挿入していないのでイメージがつかなければ絵を描く必要もあるかもしれません。
不親切に思われるかもしれませんが、手を動かさなければ理解できたと思われても自分で計算することはいつまでもできません。
なので頑張って手を動かしながら読んでいきましょう。
-----------------------------------------------
1.オイラー角を利用して速度の回転成分を求めよう
-----------------------------------------------
ここでは質点の速度の回転成分を求める事を考えます。その方法とは次の方法が考えられます。
[もとめる方法]
{1}時刻 $t$ の質点の位置 $\bm{x}$ $($ $t$ $)$ がある一つの軸上にのるような座標 $\bm{x"}$ を極座標変数 $r$ , $\theta$, $\phi$ を用いて表す。
これは動径方向以外はゼロになるはずです。極座標変数で表すのは動径成分と回転成分を分離させる狙いがあります。
{2}時刻 $t$ $+$ $\delta$ $t$ の質点の位置 $x$ $($ $t$ $+$ $\delta$ $t$ $)$ をある一つの軸上に乗るような
座標 ${x""}$ で時刻 $t$ $+$ $\delta$ $t$ の質点の位置と時刻 $t$ の質点の位置を極座標変数を用いて表す。
{3}{2}で求めた式を用いて時間微分を行い質点の座標 ${x""}$ 上での速度成分を求める。
これによって速度の回転成分と動径成分がそれぞれもとまります。
{1}
回転による座標変換にはオイラー角を利用するのが便利です。
まず時刻tにおける質点の位置 $x$ を直交直線座標で表したものを極座標変数で表現すると次のようになります。
x &= r \sin\theta \cos\phi \tag{1}\\
y &= r \sin\theta \sin\phi \tag{2}\\
z &= r \cos\theta \tag{3}
次に時刻tにおける質点の位置が ${z"}$ 軸上になるような座標変換を実行します。これはオイラー角を利用して $z$ 軸周りに $\phi$
だけ回転し、次に ${y'}$ 軸周りに $\theta$ だけ回転させてやれば良いわけです。結果は ${z"}$ 成分が $r$ でそれ以外はゼロになる事は
図を書けば計算しなくても分かります。しかし本当に考えたとおりになっているのかオイラー角の計算に慣れるためにも計算してみましょう。
計算する方法は次のとおりです。計算は各自でやってみてください。ここではその方法のみを示しておきます。
{計算方法}
[1]z軸周りに $\phi$ だけ回転
\begin{pmatrix}
x'\\
y'\\
z'
\end{pmatrix}
&=
\begin{pmatrix}
\cos \phi & \sin \phi & 0\\
-\sin \phi & \cos \phi & 0\\
0 & 0 & 1
\end{pmatrix}
\begin{pmatrix}
x\\
y\\
z
\end{pmatrix} \tag{4}
[2] ${y'}$ 軸周りに $\theta$ だけ回転
\begin{pmatrix}
x"\\
y"\\
z"
\end{pmatrix}
&= \begin{pmatrix}
\cos \theta & -\sin \theta & 0 \\
0 & 1 & 0 \\
\sin \theta & 0 & \cos \theta
\end{pmatrix}
\begin{pmatrix}
x'\\
y'\\
z'
\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\\
y\\
z
\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\\
y\\
z
\end{pmatrix} \tag{5}
[3] $\tag{5}$ 式に $\tag{1}$ , $\tag{2}$ , $\tag{3}$ 式を代入
\begin{pmatrix}
x"\\
y"\\
z"
\end{pmatrix}
&=
\begin{pmatrix}
0\\
0\\
r
\end{pmatrix} \tag{6}
[1],[2],[3]の計算を実行すれば予測された結果である $\tag{6}$ 式が得られます。
これで $\bm{x"}$ は ${z"}$ 軸上に質点が時刻 $t$ にあるような座標系になることが確かめられました。
{2}
次に時刻 ${t+ \delta t}$ に質点が ${z""}$ 軸上に乗るような座標 $\bm{x""}$ を極座標変数を用いて書き表します。
この座標変換は $x$ を $z$ を中心に ${\phi + \delta \phi}$ , ${y"'}$ 軸を中心に ${\theta + \delta \theta}$ だけ回転させることにあたります。
ここに書いている ${\delta \phi}$ , ${\delta \theta}$ は微小な大きさの角です。
この変換は $\tag{5}$ 式を $\theta$ → $\theta$ $+$ $\delta$ $\theta$ , $\phi$ → $\phi$ $+$ $\delta$ $\phi$ , $\bm{x"}$ → $\bm{x""}$ と
書き換えた式から求まります。その事は{1}と同じ理由から分かります。
\begin{pmatrix}
x""\\
y""\\
z""
\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\\
y\\
z
\end{pmatrix} \tag{5'}
ここで三角関数の
\sin \psi \simeq \psi,\cos \psi \simeq 1 (\psi << 1) \tag{7}
という近似を用いると
\cos (\psi + \delta \psi) &= \cos \psi \cos \delta \psi + \sin \psi \sin \delta \psi\\
&= \cos \psi + \delta \psi \sin \psi \tag{8}\\
\sin (\psi + \delta \psi) &= \sin \psi \sin \delta \psi - \sin \delta \psi \cos \psi\\
&= \sin \psi - \delta \psi \cos \psi \tag{9}
が得られます。 $\tag{5'}$ 式に $\tag{8}$ , $\tag{9}$ にそれぞれ ${\theta + \delta \theta}$ , ${\phi+ \delta \phi}$ を代入したものを放り込む事によって
\begin{pmatrix}
x""\\
y""\\
z""
\end{pmatrix}
&= \begin{pmatrix}
0\\
0\\
r
\end{pmatrix}
+ \delta \theta
\begin{pmatrix}
-\sin \theta \cos \phi & -\sin \theta \sin \phi & -\cos \theta \\
0 & 0 & 0 \\
\sin \theta \cos \phi & \cos \theta \sin \phi & -\sin \theta
\end{pmatrix}
\begin{pmatrix}
x\\
y\\
z
\end{pmatrix}
+ \delta \phi
\begin{pmatrix}
-\cos \theta \sin \phi & \cos \theta \cos \phi & 0 \\
-\cos \phi & -\sin \phi & 0 \\
-\sin \theta \sin \phi & \sin \theta \cos \phi & 0
\end{pmatrix}
\begin{pmatrix}
x\\
y\\
z
\end{pmatrix}
となることが分かります。更に(1),(2),(3)を代入すると
\begin{pmatrix}
x""\\
y""\\
z""
\end{pmatrix}
&=\begin{pmatrix}
-r \delta \theta\\
-r \delta \phi\\
r
\end{pmatrix} \tag{10}
を得ます。こうして時刻 $t$ の位置 $\bm{x}$ を座標 $\bm{x""}$ 上で極座標を用いて表す事ができました。
{3}
時刻tに質点は $z"$ 上あり、時刻 ${t+ \delta t}$ に $z""$ 上にあるので。このときの座標 $\bm{x""}$ から見た時の速度 $\bm{v}$ は微分の定義から
\bm{v} &= \lim_{\delta t \to 0} \frac{{\bm{x""}}(t + \delta t) - \bm{x""}(t)}{\delta t}\\
&=
\begin{pmatrix}
0 \\
0 \\
\dot{r}
\end{pmatrix}
+ \begin{pmatrix}
r \dot{\theta} \\
r \dot{\phi} \sin \theta \\
0
\end{pmatrix}
だと分かります。 $z""$ 軸方向上に質点があるので右辺第一項は動径方向を表しています。従って動径方向成分を除いた成分が
回転にかかわる速度になるわけです。座標 ${x""}$ 上では速度の回転成分を極座標変数を用いると
\bm{v}_{\omega}
&=
r\begin{pmatrix}
\dot{\theta} \\
\dot{\phi} \sin \theta \\
0
\end{pmatrix} \tag{12}
と表されます。
@@author: おこめ@@
@@accept: ?@@