- 追加された行はこの色です。
- 削除された行はこの色です。
#rst2hooktail_source
============================================================
剛体のオイラー角でのハミルトニアンを解く
============================================================
剛体の回転シリーズ番外編３です。
せっかく番外編２で剛体のハミルトニアンを求めたので、
剛体のハミルトニアンを解いて剛体の運動方程式を導いてみました。
復習
===================
まず、ハミルトニアンを確認します。
剛体のハミルトニアンを次のようなものでした。
<tex>
H &=\frac{1}{2 I_x \sin^2 \theta}\{ (p_\phi - \cos \theta p_\psi) \cos \psi - \sin \theta \sin \psi p_\theta \}^2 \\
&+ \frac{1}{2 I_y \sin^2 \theta}\{ (p_\phi - \cos \theta p_\psi)\sin \psi + \sin \theta \cos \psi p_\theta \}^2 \\
&+ \frac{p_\theta^2}{2 I_z} \tag{##}
</tex>
パラメータ $\lambda$ に対して、
<tex>
\dot{p}_\lambda = - \frac{\partial H}{\partial \lambda} \tag{##}
</tex>
<tex>
\dot{\lambda} = \frac{\partial H}{\partial p_\lambda} \tag{##}
</tex>
です。
それでは、さっそく式 $(2)$ を求めてみましょう。
<tex>
\dot{p}_\theta &= - \frac{\partial H}{\partial \theta} \\
&= \frac{\cos \theta}{I_x \sin^3 \theta}\{ (p_\phi - \cos \theta p_\psi) \cos \psi - \sin \theta \sin \psi p_\theta \} \\
&= \frac{\cos \theta}{I_x \sin^3 \theta}\{ (p_\phi - \cos \theta p_\psi) \cos \psi - \sin \theta \sin \psi p_\theta \}^2 \\
&- \frac{1}{I_x \sin^2 \theta}\{ (p_\phi - \cos \theta p_\psi) \cos \psi - \sin \theta \sin \psi p_\theta \}
(\cos \psi \sin \theta p_\psi - \cos \theta \sin \psi p_\theta) \\
&+ \frac{\cos \theta}{I_y \sin^3 \theta}\{ (p_\phi - \cos \theta p_\psi) \sin \psi + \sin \theta \cos \psi p_\theta \} \\
&+ \frac{\cos \theta}{I_y \sin^3 \theta}\{ (p_\phi - \cos \theta p_\psi) \sin \psi + \sin \theta \cos \psi p_\theta \}^2 \\
&- \frac{1}{I_y \sin^2 \theta}\{ (p_\phi - \cos \theta p_\psi) \cos \psi - \sin \theta \cos \psi p_\theta \}
(\sin \psi \sin \theta p_\psi + \cos \theta \cos \psi p_\theta) \tag{##}
</tex>
ここで、次のように $\alpha , \beta , \gamma$ を定義します。
<tex>
\alpha = \frac{\cos^2 \psi}{I_x}+\frac{\sin^2 \psi}{I_y} \tag{##}
</tex>
<tex>
\beta = \frac{1}{I_y}-\frac{1}{I_x} \tag{##}
</tex>
<tex>
\gamma = \frac{\sin^2 \psi}{I_x}+ \frac{\cos^2 \psi}{I_y} \tag{##}
</tex>
すると、式 $(4)$ の続きは、
<tex>
\dot{p}_\theta &= \frac{\cos \theta}{\sin^3 \theta}\alpha p_\phi^2 \\
&+ \frac{\sin \psi \cos \psi \cos \theta}{\sin^2 \theta}\beta p_\phi p_\theta \\
&- \frac{1 + \cos^2 \theta}{\sin^3 \theta} p_\phi p_\psi \\
&+ 0 \times p_\theta^2 \\
&- \frac{\sin \psi \cos \psi}{\sin^2 \theta} \beta p_\theta p_\psi \\
&+ \frac{\cos \theta}{\sin^3 \theta} \alpha p_\psi^2 \tag{##}
</tex>
@@author:クロメル@@
@@accept:2010-03-03@@
@@category:力学@@
@@id:equationOfRigidHamiltonian@@
:記事ソース/剛体のオイラー角でのハミルトニアンを解く のバックアップ差分(No.4)
