回転運動の物理を深く理解しようとするとき、1軸まわりの回転だけでなく、3次元空間での任意の方向への回転を扱う必要が生じます。
そのとき、単純なスカラー量としての慣性モーメントでは記述が不十分となり、慣性テンソルという概念が必要になります。
慣性テンソルは3×3の行列として表現され、物体の質量分布の3次元的な非対称性をすべて含んだ情報を持ちます。
主軸・対角化・慣性楕円体など、線形代数と物理学が交差するこの分野は、剛体の3次元回転を正確に記述するための重要な基盤です。
本記事では、慣性テンソルの定義から行列表現・主軸の意味・対角化の方法・剛体の回転運動への応用まで、体系的にわかりやすく解説します。
慣性モーメントとテンソルの関係:結論と全体像
それではまず、慣性モーメントと慣性テンソルの関係の全体像について解説していきます。
1次元的な回転では「慣性モーメント」はスカラー値(単一の数)として扱われます。
しかし3次元空間での回転を考えると、角速度ベクトル ω と角運動量ベクトル L の関係が単純な比例関係にならないことがあります。
慣性テンソルの本質的な意味
一般の3次元回転では:L = I ω
ここで I は慣性テンソル(3×3行列)であり、スカラーではない。
この関係式が成り立つのは、物体の形状が回転軸に対して非対称な場合、角速度ベクトルと角運動量ベクトルの方向が一致しないためである。
つまり、慣性テンソルは「どの方向に回転させようとすると、角運動量がどの方向にどれだけ生じるか」を完全に記述する量です。
3×3の行列として表現される慣性テンソルの各要素は、慣性モーメント(対角成分)と慣性積(非対角成分)から構成されます。
慣性テンソルを適切な軸(主軸)まわりに表現すると対角行列になり、この状態が最も計算しやすい形です。
主軸まわりでは角速度と角運動量が平行になり、回転運動の解析が大幅に簡単になります。
慣性テンソルの定義と行列表現
続いては、慣性テンソルの具体的な定義と行列表現について確認していきます。
慣性テンソルの各成分の定義
慣性テンソルは次の3×3対称行列として定義されます。
慣性テンソル I の行列表現:
| Ixx -Ixy -Ixz |
| -Iyx Iyy -Iyz |
| -Izx -Izy Izz |
対角成分(慣性モーメント):
Ixx = ∫(y² + z²)dm
Iyy = ∫(x² + z²)dm
Izz = ∫(x² + y²)dm
非対角成分(慣性積):
Ixy = Iyx = ∫xy dm
Iyz = Izy = ∫yz dm
Ixz = Izx = ∫xz dm
対角成分の Ixx・Iyy・Izz は、それぞれ x 軸・y 軸・z 軸まわりの慣性モーメントです。
これらは常に正の値を持ちます。
非対角成分の慣性積(慣性乗積とも呼ぶ)は、質量分布の非対称性を表し、正負いずれの値も取りえます。
物体が座標軸に関して対称な質量分布を持つ場合、慣性積はゼロになります。
慣性テンソルの対称性と固有値
慣性テンソルが対称行列であることは、Ixy = Iyx などの関係から明らかです。
対称行列には重要な性質があり、必ず実数の固有値を持ち、異なる固有値に対応する固有ベクトルは直交します。
この性質が主軸の概念と直接結びついています。
慣性テンソルの固有値 I₁, I₂, I₃ を主慣性モーメントと呼び、対応する固有ベクトルの方向を主軸と呼びます。
主慣性モーメントは常に正の値(I₁, I₂, I₃ > 0)であり、互いに直交する3方向に関連します。
慣性テンソルの座標変換
慣性テンソルの値は、座標系の取り方によって変化します。
座標系を回転行列 R で変換したとき、慣性テンソルは次のように変換されます。
I’ = R I Rᵀ
ここで Rᵀ は R の転置行列
この変換則はテンソルの定義そのものであり、慣性テンソルが真のテンソル量であることを示している。
この変換は、異なる座標系で記述された慣性テンソルを統一的に扱うための基礎となります。
主軸と慣性テンソルの対角化
続いては、主軸の物理的意味と慣性テンソルの対角化について確認していきます。
主軸の物理的意味と求め方
主軸は、慣性テンソルを対角行列に変換する座標系の軸です。
主軸まわりに物体を回転させると、角速度ベクトルと角運動量ベクトルが平行になります。
これは物理的には「主軸まわりの回転では、軸受けに横方向の力がかからない」ことを意味します。
主軸は固定点(または重心)に対して物体に固定された特別な軸であり、多くの場合、物体の幾何学的対称軸と一致します。
主軸を求めるには、慣性テンソル行列の固有値問題を解きます。
主軸を求める固有値方程式:
I n = λ n
(I – λE)n = 0
det(I – λE) = 0 (特性方程式)
ここで n は主軸方向の単位ベクトル(固有ベクトル)
λ は主慣性モーメント(固有値)
E は単位行列
3次の行列式を展開すると3次方程式が得られ、その3つの実数解が主慣性モーメント I₁, I₂, I₃ になります。
対角化の手順と計算例
慣性テンソルの対角化を具体的な手順で確認しましょう。
対角化の手順:
①慣性テンソルの行列を計算する
②特性方程式 det(I – λE) = 0 を解いて固有値(主慣性モーメント)λ₁, λ₂, λ₃ を求める
③各固有値に対して固有ベクトル(主軸方向)n₁, n₂, n₃ を求める
④固有ベクトルを列に並べた直交行列 P を構成する
⑤対角化された慣性テンソル:PᵀIP = diag(I₁, I₂, I₃)
対角化後の慣性テンソルは対角成分に主慣性モーメントのみを持ち、慣性積はすべてゼロになります。
この主軸座標系が、剛体の回転運動を解析するために最も便利な座標系です。
対称性と主軸の関係
物体が持つ幾何学的対称性は、主軸の方向を決定するのに役立ちます。
物体が対称面を持つ場合、その面に垂直な方向が必ず主軸の一つになります。
物体が回転対称軸を持つ場合、その軸が主軸の一つです。
球対称な物体(均一な球など)では、任意の方向が主軸となり、3つの主慣性モーメントはすべて等しくなります。
| 物体の対称性 | 主軸の特徴 | 例 |
|---|---|---|
| 球対称 | 任意の方向が主軸、I₁=I₂=I₃ | 均一な球 |
| 軸対称(回転体) | 回転軸が主軸、I₁=I₂≠I₃ | 円柱・円錐 |
| 3つの対称面 | 各対称面の法線が主軸 | 直方体 |
| 対称性なし | 計算で求める必要あり | 不規則形状 |
剛体の3次元回転とオイラー方程式
続いては、慣性テンソルを使った剛体の3次元回転の解析とオイラー方程式について確認していきます。
オイラーの運動方程式
剛体の回転運動の方程式は、主軸座標系(体固定座標系)で記述するとオイラーの運動方程式として表されます。
オイラーの運動方程式(主軸座標系):
I₁ω̇₁ – (I₂ – I₃)ω₂ω₃ = N₁
I₂ω̇₂ – (I₃ – I₁)ω₃ω₁ = N₂
I₃ω̇₃ – (I₁ – I₂)ω₁ω₂ = N₃
ここで I₁, I₂, I₃:主慣性モーメント
ω₁, ω₂, ω₃:各主軸まわりの角速度成分
N₁, N₂, N₃:各主軸まわりのトルク成分
このオイラー方程式は一般的に非線形連立常微分方程式であり、解析的な一般解を持たないことが多いです。
主軸を使って記述することで、慣性積が消えて方程式が最もシンプルな形になります。
トルクなし自由回転とオイラーの運動方程式
外部トルクがゼロ(N₁ = N₂ = N₃ = 0)の場合、オイラー方程式は次のように簡単になります。
自由回転のオイラー方程式:
I₁ω̇₁ = (I₂ – I₃)ω₂ω₃
I₂ω̇₂ = (I₃ – I₁)ω₃ω₁
I₃ω̇₃ = (I₁ – I₂)ω₁ω₂
この方程式系には、エネルギー保存と角運動量保存の2つの保存則から2つの保存量が存在します。
主慣性モーメントがすべて異なる場合(I₁ ≠ I₂ ≠ I₃)の自由回転は「非対称コマ」と呼ばれ、その運動は楕円積分で記述されます。
2つの主慣性モーメントが等しい場合(I₁ = I₂ ≠ I₃)の自由回転は「対称コマ」と呼ばれ、首振り(歳差)運動が生じます。
人工衛星とジャイロスコープへの応用
慣性テンソルと主軸の概念は、人工衛星の姿勢制御やジャイロスコープの挙動解析に直接応用されます。
人工衛星の安定した姿勢は、主慣性モーメントの最大軸または最小軸まわりの回転では達成できますが、中間軸まわりの回転は不安定です。
これを「中間軸の定理」または「テニスラケットの定理」と呼び、日常的に観察できる物理現象として有名です。
コマの歳差運動・章動・自転の3種類の角運動はすべてオイラー角で記述され、慣性テンソルの主軸分析が不可欠です。
慣性楕円体と幾何学的表現
続いては、慣性テンソルの幾何学的な可視化ツールである慣性楕円体について確認していきます。
慣性楕円体の定義と性質
慣性楕円体は、任意の方向 n に対する慣性モーメントを幾何学的に表現するツールです。
方向 n に対する慣性モーメントを I(n) とすると、1/√I(n) を半径方向の距離として描いた楕円体が慣性楕円体です。
慣性楕円体の方程式(主軸座標系):
I₁x² + I₂y² + I₃z² = 1
半軸の長さ:1/√I₁, 1/√I₂, 1/√I₃
主軸がそのまま楕円体の軸と一致する。
慣性楕円体の主軸方向が物体の主軸と一致し、主軸の長さが主慣性モーメントの逆数の平方根になります。
ポアンソーの幾何学的描像
フランスの数学者ポアンソーは、自由回転する剛体の運動を純粋に幾何学的に描写する方法を考案しました。
この描像では、剛体に固定された慣性楕円体が、不変平面(角運動量に垂直な平面)の上を滑らずに転がる運動として自由回転が表現されます。
ポアンソーの描像は、複雑な3次元回転運動を直感的に理解するための強力な幾何学的ツールです。
エネルギーの等高線と角運動量の保存条件がポロダルとエルポロダルという2つの曲線として現れ、それらの接触点の軌跡がコマの運動を記述します。
数値計算・シミュレーションへの応用
現代の工学計算では、剛体の3次元回転を数値シミュレーションで解析することが一般的です。
有限要素法(FEM)ソフトウェアや多体動力学シミュレーターでは、慣性テンソルの計算と主軸の決定が自動的に行われます。
四元数(クォータニオン)を使った姿勢表現と慣性テンソルを組み合わせることで、特異点(ジンバルロック)を回避しながら正確な3次元回転シミュレーションが可能になります。
| 応用分野 | 使用される概念 | 具体的な応用例 |
|---|---|---|
| 航空宇宙工学 | 慣性テンソル・オイラー方程式 | 人工衛星の姿勢制御設計 |
| ロボット工学 | 主軸・慣性テンソル変換 | 多関節ロボットの動力学計算 |
| ゲーム物理エンジン | 主慣性モーメント・対角化 | リアルな剛体シミュレーション |
| 医療工学 | 慣性テンソル解析 | 人体運動の動力学モデリング |
| 車両工学 | 慣性テンソル・主軸 | 車両ロールオーバー解析 |
まとめ
本記事では、慣性モーメントとテンソルの関係から始まり、慣性テンソルの定義・行列表現・主軸・対角化・剛体の3次元回転への応用まで解説しました。
慣性テンソルは3×3の対称行列であり、対角成分が慣性モーメント、非対角成分が慣性積を表します。
主軸を基準とした座標系では慣性テンソルが対角行列となり、回転運動の解析が大幅に簡単になるのがポイントです。
オイラーの運動方程式・慣性楕円体・ポアンソーの描像など、慣性テンソルに関連する概念は、現代の工学シミュレーションや宇宙工学にも直接応用されています。
線形代数と物理学が融合したこの分野をしっかり理解することで、3次元回転運動の解析力が大きく向上するでしょう。