[論文レビュー] General solution to the Euler-Poisson equations of a free Lagrange top directly for the rotation matrix

Alexei A. Deriglazov|arXiv (Cornell University)|Mar 4, 2023

Aerospace Engineering and Control Systems被引用数 5

ひとこと要約

本稿では、オイラー・ポisson方程式に対する自由ラグランジュトップの一般解を、オイラー角などのパラメータ化を用いずに、回転行列 $ R_{ij}(t) $ の形で直接提示する。解は三角関数を用いて表現され、回転行列とポインソーの運動の幾何的描写を明示的に結びつける。これにより、運動量の軸まわりの回転と慣性系のz軸まわりの回転への分解が明らかになる。

ABSTRACT

The Euler-Poisson equations para determinar the rotation matrix of a rigid body can be solved without using of particular parameterization like the Euler angles. For the free Lagrange top, we obtain and discuss a general analytic solution, and compare it with the Poinsot's picture of motion.

研究の動機と目的

• 自由剛体のオイラー・ポisson方程式の一般解析解を、オイラー角のようなパラメータ化に依存せずに導出すること。
• 運動方程式から直接、回転行列 $ R_{ij}(t) $ の時間発展を、完全な直交構造を保ったまま表現すること。
• 解をポインソーの運動の幾何的描写（特に慣性楕円体が不変平面に滑らかに転がる様子）と結びつけること。
• 回転速度と1つの基底ベクトルの解から、回転行列全体を再構成できることを示すこと。
• 角度に基づくパrametrizationに内在する特異性を避ける、グローバルに有効な非局所的剛体回転の記述を提供すること。

提案手法

• ラグランジュ作用素を用いて剛体の力学を制約系として定式化し、$ R_{ij}(t) $ と $ \boldsymbol{\Omega}(t) $ のオイラー・ポisson方程式を導出する。
• 回転エネルギー（3）、運動量（4）、直交制約（5）といった保存量を用いて解の空間を制約する。
• 全系を最小限のセットに還元：基底ベクトル $ \boldsymbol{\gamma}(t) $ の1つのベクトル発展方程式（9）とオイラー方程式（8）を解く。
• 初期条件 $ (1,0,0), (0,1,0), (0,0,1) $ を満たす3つの正規直交基底ベクトル $ \mathbf{G}_1(t), \mathbf{G}_2(t), \mathbf{G}_3(t) $ を解いて、全回転行列 $ R(t) $ を構成する。
• 初期運動量成分と慣性モーメントに明示的な依存関係を持つ三角関数の形で解を表現する。
• 回転行列を $ R(t) = R_{\mathbf{m}}(t) \times R_{OZ}(t) $ の形に分解し、運動量の軸まわりの回転と慣性系z軸まわりの回転を順次表す。

実験結果

リサーチクエスチョン

• RQ1 自由剛体のオイラー・ポisson方程式は、オイラー角のようなパラメータ化を用いずに、回転行列 $ R_{ij}(t) $ に対して直接解けるか？
• RQ2 一般解 $ R_{ij}(t) $ はポインソーの剛体運動の幾何的記述とどのように関係するか？
• RQ3 三角関数と物理定数を用いて、ラグランジュトップの回転行列の明示的形は何か？
• RQ4 回転速度と1つの基底ベクトルの解からなる低次元系の解から、全回転行列をどのように再構成できるか？
• RQ5 回転行列の分解 $ R(t) = R_{\mathbf{m}}(t) \times R_{OZ}(t) $ の幾何的解釈は、円錐上での運動と慣性楕円体の転がりにどのように対応するか？

主な発見

• 回転行列 $ R_{ij}(t) $ の一般解が、時間全域で有効で特異性のない三角関数の形で明示的に導出された。
• 初期条件と $ R_{ki}R_{kj} = \delta_{ij} $ の保存量のおかげで、直交制約 $ R^T R = I $ が自動的に満たされる。
• 回転行列は $ R(t) = R_{\mathbf{m}}(t) \times R_{OZ}(t) $ の形に分解され、$ R_{\mathbf{m}}(t) $ は運動量ベクトルの軸まわりの回転、$ R_{OZ}(t) $ は周波数 $ \phi = \frac{(I_1 - I_3)m_3}{I_1 I_3} t $ でz軸まわりの回転を表す。
• 瞬間的な角速度 $ \boldsymbol{\omega}(t) $ が明示的に導出され、$ m_2, m_3, I_1, I_3 $ および時間に依存する余弦項に依存する。
• 角速度の大きさは $ \boldsymbol{\omega}^2 = k^2 + \frac{I_1 + I_3}{I_1 I_3} m_3 \phi $ を満たし、$ k = \frac{|\mathbf{m}|}{I_1} $ である。これにより複合回転周波数が確認される。
• 解はポインソーの描写と完全に一致する：物体固定軸の先端は円錐上を円を描き、慣性楕円体は不変平面に滑らかに転がる。