固定された回転軸に対する回転慣性のスカラー値
慣性モーメント フライホイールは 、回転運動の速度の変化を滑らかにするために、大きな慣性モーメントを備えています。 一般的な記号
私 SI単位 kg⋅m 2 その他のユニット
lbf·ft·s 2 他の量 からの導出
私
=
ら
ω
{\displaystyle I={\frac {L}{\omega }}}
寸法 M L 2
綱渡り師は、 長い棒の慣性モーメントを利用してバランスを取りながら綱を渡りながら歩きます。1890 年に ナイアガラ川を渡るサミュエル ディクソン。
戦闘機は機動性を向上させるために慣性モーメントを最小限に抑えるように設計されていますが、民間航空機ではそうではない場合がよくあります。
慣性モーメントは 、 質量慣性モーメント 、 角/回転質量 、 質量の二次モーメント 、または最も正確には 回転慣性 とも呼ばれ 、回転 軸に対して定義されます。これは、その軸の周りで適用される トルク と結果として生じる 角加速度 の比率です。これは、直線運動における 質量 と同じ役割を果たします 。特定の軸の周りの物体の慣性モーメントは、質量と軸に対する質量の分布の両方に依存し、質量と軸からの距離とともに増加します。
これは拡張可能な (加法的)特性 です。 質点 の場合、慣性モーメントは単に質量に回転軸に対する 垂直距離 の 2 乗を掛けた値です。剛体複合システムの慣性モーメントは、その構成要素であるサブシステムの慣性モーメントの合計です(すべて同じ軸を中心にとった値です)。最も単純な定義は、 軸 からの距離に対する質量の2 次 モーメント です。
平面内で回転するように制約された物体の場合、平面に垂直な軸の周りの慣性モーメント ( スカラー値) のみが重要です。3 次元内で自由に回転する物体の場合、モーメントは 対称 の 3 行 3 列のマトリックス で記述できます。このマトリックスは、互いに垂直な主軸のセットに対して 対角線 になり、軸の周りのトルクは互いに独立して作用します。
機械工学 では、単に「慣性」は「 慣性質量 」または「慣性モーメント」 を指すことが多い。 [1]
導入
物体が軸の周りを自由に回転する場合、 角運動量を変更するには トルク を適用する必要があります。任意の 角加速度( 角速度 の変化率) を 発生させるために必要なトルクの量は、 物体の慣性モーメントに比例します。慣性モーメントは、 SI単位ではキログラム メートルの 2 乗 (kg·m 2 ) 、ヤードポンド 法 または 米国 単位ではポンド フィート秒の 2 乗 (lbf·ft·s 2 ) で表されます 。
慣性モーメントは、線形運動学における質量 (慣性)と同じ役割を回転 運動学で果たします。どちらも、物体の運動の変化に対する抵抗を特徴付けます。慣性モーメントは、回転軸の周囲に質量がどのように分布しているかによって決まり、選択した軸によって異なります。点状の質量の場合、ある軸の周りの慣性モーメントは で与えられます。 ここで、 は軸からの点の距離、 は 質量です。拡張された剛体の場合、慣性モーメントは、すべての小さな質量片の合計に、回転軸からの距離の 2 乗を掛けた値になります。規則的な形状と均一な密度の拡張された物体の場合、この合計から、物体の寸法、形状、総質量に応じた簡単な式が生成されることがあります。
m
r
2
{\displaystyle mr^{2}}
r
{\displaystyle r}
m
{\displaystyle m}
1673年、 クリスティアーン・ホイヘンスは、 軸から吊るされた物体の振動に関する研究でこのパラメータを導入しました 。この物体は 複合振り子として知られています。 [2] 慣性モーメント( ラテン語 で「momentum inertiae」 ) という用語は、 レオンハルト・オイラー が 1765年 に著書 「固体の運動は剛体である」で導入し、 [2] [3] オイラーの第二法則 に組み込まれています 。
複合振り子の固有振動数は、重力によって振り子の質量に課されるトルクと、慣性モーメントによって定義される加速度に対する抵抗の比から得られる。この固有振動数を単一の質量点からなる単純振り子の固有振動数と比較すると、拡張された物体の慣性モーメントの数学的定式化が得られる。 [4] [5]
慣性モーメントは、 運動量 、 運動エネルギー 、および ニュートンの剛体運動の法則 にも現れ、その形状と質量を組み合わせた物理的パラメータとして表されます。慣性モーメントが平面運動と空間運動で現れる方法には興味深い違いがあります。平面運動では慣性モーメントを定義する単一のスカラーがありますが、空間運動では同じ計算で慣性モーメントの 3 × 3 行列が得られ、これは慣性行列または慣性テンソルと呼ばれます。 [6] [7]
回転するフライホイール の慣性モーメント は、機械内で適用されるトルクの変化に抵抗して回転出力を滑らかにするために使用されます。飛行機の縦軸、水平軸、垂直軸周りの慣性モーメントは、翼、昇降舵、方向舵の制御面にかかる操舵力が、飛行機のロール、ピッチ、ヨーの動きにどのような影響を与えるかを決定します。
意味 慣性モーメントは、
断面 の質量と、基準軸と 断面の 重心間の距離の二乗との積として定義されます。
回転するフィギュアスケーターは、腕を引くことで慣性モーメントを減らし、 角運動量保存則 によりより速く回転することができます。
回転椅子の実験のビデオ。慣性モーメントを示しています。回転する教授が腕を引くと、慣性モーメントが減少し、角運動量を保存するため、角速度が増加します。
慣性モーメント Iは 、系の正味 角運動量 L と主軸周りの 角速度 ω の比としても定義され、 [8] [9] 、
I = L ω . {\displaystyle I={\frac {L}{\omega }}.}
システムの角運動量が一定であれば、慣性モーメントが小さくなるにつれて、角速度は増加しなければなりません。これは、回転する フィギュアスケート選手が 伸ばした腕を引き締めたり、 ダイバーが 飛び込み中に体を丸めて 回転 を速くしたりするときに起こります。 [8] [9] [10] [11] [12] [13] [14]
物体の形状が変化しない場合、その慣性モーメントは、 ニュートンの運動の法則 において、物体に 加えられたトルク τ と主軸の周りの
角加速度 α の比として現れます。つまり、
τ = I α です。{\displaystyle \tau =I\alpha .}
単純な振り子 の場合、この定義により、振り子の質量 m と支点からの
距離 r に関する 慣性モーメント Iの式は次のようになります。
I
=
m
r
2
.
{\displaystyle I=mr^{2}.}
したがって、振り子の慣性モーメントは、物体の 質量 m と、回転軸までの
距離 rによって定義されるその形状、つまり幾何学的形状の両方に依存します。
この簡単な公式を一般化すると、任意の形状の物体の慣性モーメントは、すべての基本質量点dm の合計に、 軸 kに対する垂直距離 r の 2 乗を乗じたものとして定義されます 。したがって、任意の物体の慣性モーメントは、その質量の空間分布によって異なります。
一般に、質量m の物体がある場合、 特定の回転軸に依存して、軸周りの慣性モーメントが k となるような値で有効
半径 k を定義できます。
ここで、 kは軸周り の回転半径 と呼ばれます 。
I
=
m
k
2
,
{\displaystyle I=mk^{2},}
例
単振り子
数学的には、単振り子の慣性モーメントは、振り子の支点周りの重力によるトルクと、その支点周りの角加速度の比です。単振り子の場合、これは粒子の質量 と支点からの距離の2乗の積であることが分かっています 。
m
{\displaystyle m}
r
{\displaystyle r}
I
=
m
r
2
.
{\displaystyle I=mr^{2}.}
これは次のように表すことができます。単純な振り子の質量にかかる重力は、 振り子の運動面に対して垂直な軸の周りにトルクを生成します。ここで、 はトルク軸から振り子の質量中心までの距離ベクトルであり、 は 質量にかかる正味の力です。このトルクに関連して、弦と質量のこの軸周りの 角加速度 、が あります。質量は円に制限されているため、質量の接線加速度は です 。 トルク方程式は次のようになります。
τ
=
r
×
F
{\displaystyle {\boldsymbol {\tau }}=\mathbf {r} \times \mathbf {F} }
r
{\displaystyle \mathbf {r} }
F
{\displaystyle \mathbf {F} }
α
{\displaystyle {\boldsymbol {\alpha }}}
a
=
α
×
r
{\displaystyle \mathbf {a} ={\boldsymbol {\alpha }}\times \mathbf {r} }
F
=
m
a
{\displaystyle \mathbf {F} =m\mathbf {a} }
τ
=
r
×
F
=
r
×
(
m
α
×
r
)
=
m
(
(
r
⋅
r
)
α
−
(
r
⋅
α
)
r
)
=
m
r
2
α
=
I
α
k
^
,
{\displaystyle {\begin{aligned}{\boldsymbol {\tau }}&=\mathbf {r} \times \mathbf {F} =\mathbf {r} \times (m{\boldsymbol {\alpha }}\times \mathbf {r} )\\&=m\left(\left(\mathbf {r} \cdot \mathbf {r} \right){\boldsymbol {\alpha }}-\left(\mathbf {r} \cdot {\boldsymbol {\alpha }}\right)\mathbf {r} \right)\\&=mr^{2}{\boldsymbol {\alpha }}=I\alpha \mathbf {\hat {k}} ,\end{aligned}}}
ここで、 は 振り子の平面に垂直な単位ベクトルです。 (最後から2番目のステップでは、 と の垂直性を持つ ベクトル三重積展開を 使用します。) この量は、ピボットポイントの周りのこの単一の質量の 慣性モーメント です 。
k
^
{\displaystyle \mathbf {\hat {k}} }
α
{\displaystyle {\boldsymbol {\alpha }}}
r
{\displaystyle \mathbf {r} }
I
=
m
r
2
{\displaystyle I=mr^{2}}
この量は、 単純な振り子の 角運動量 にも現れます。これは、 振り子の質量のピボット周りの速度から計算されます。ここで、は、質量のピボットポイント周りの 角速度 です。この角運動量は
、前の式と同様の導出を使用して
求められます
。
I
=
m
r
2
{\displaystyle I=mr^{2}}
v
=
ω
×
r
{\displaystyle \mathbf {v} ={\boldsymbol {\omega }}\times \mathbf {r} }
ω
{\displaystyle {\boldsymbol {\omega }}}
L
=
r
×
p
=
r
×
(
m
ω
×
r
)
=
m
(
(
r
⋅
r
)
ω
−
(
r
⋅
ω
)
r
)
=
m
r
2
ω
=
I
ω
k
^
,
{\displaystyle {\begin{aligned}\mathbf {L} &=\mathbf {r} \times \mathbf {p} =\mathbf {r} \times \left(m{\boldsymbol {\omega }}\times \mathbf {r} \right)\\&=m\left(\left(\mathbf {r} \cdot \mathbf {r} \right){\boldsymbol {\omega }}-\left(\mathbf {r} \cdot {\boldsymbol {\omega }}\right)\mathbf {r} \right)\\&=mr^{2}{\boldsymbol {\omega }}=I\omega \mathbf {\hat {k}} ,\end{aligned}}}
同様に、振り子の質量の運動エネルギーは、振り子が軸の周りを回る速度によって定義され、
E
K
=
1
2
m
v
⋅
v
=
1
2
(
m
r
2
)
ω
2
=
1
2
I
ω
2
.
{\displaystyle E_{\text{K}}={\frac {1}{2}}m\mathbf {v} \cdot \mathbf {v} ={\frac {1}{2}}\left(mr^{2}\right)\omega ^{2}={\frac {1}{2}}I\omega ^{2}.}
これは、量が 、質量と物体の形状との組み合わせによって回転慣性を定義する方法であることを示しています。任意の形状の物体の慣性モーメントは、物体内のすべての質量要素の値の合計です 。
I
=
m
r
2
{\displaystyle I=mr^{2}}
m
r
2
{\displaystyle mr^{2}}
複合振り子
メンデンホール重力 計で使用された振り子 。1897 年の科学雑誌より。トーマス C. メンデンホールが 1890 年に開発した携帯型重力計は、地球の局所的な重力場の相対的な測定値を最も正確に提供しました。
複合 振り子は 、連続した形状の粒子の集合体から形成され、軸を中心に剛体回転する物体である。その慣性モーメントは、それを構成する各粒子の慣性モーメントの合計である。 [15] [16] : 395–396 [17] : 51–53 複合振り子の 固有 振動数 ( )は、その慣性モーメント に依存し 、
ここで は 物体の質量、 は局所的な重力加速度、 は 軸から物体の質量中心までの距離である。この振動周波数を小さな角変位にわたって測定することは、物体の慣性モーメントを測定する効果的な方法である。 [18] : 516–517
ω
n
{\displaystyle \omega _{\text{n}}}
I
P
{\displaystyle I_{P}}
ω
n
=
m
g
r
I
P
,
{\displaystyle \omega _{\text{n}}={\sqrt {\frac {mgr}{I_{P}}}},}
m
{\displaystyle m}
g
{\displaystyle g}
r
{\displaystyle r}
したがって、物体の慣性モーメントを決定するには、物体を都合のよい支点から吊り下げて、 目的の慣性モーメントの方向に対して垂直な平面内で自由に振動するようにし、その固有振動数または振動周期 ( ) を測定して、 を取得します。
ここで、 は振動周期 (持続時間) です (通常は複数の周期にわたって平均化されます)。
P
{\displaystyle P}
t
{\displaystyle t}
I
P
=
m
g
r
ω
n
2
=
m
g
r
t
2
4
π
2
,
{\displaystyle I_{P}={\frac {mgr}{\omega _{\text{n}}^{2}}}={\frac {mgrt^{2}}{4\pi ^{2}}},}
t
{\displaystyle t}
振動の中心
複合振り子と同じ固有振動数を持つ単振り子は、 支点から複合振り子の 振動中心と呼ばれる点までの長さを定義します。この点は、 打撃中心 とも一致します 。長さは、
次の
式から決定されます。
L
{\displaystyle L}
L
{\displaystyle L}
ω
n
=
g
L
=
m
g
r
I
P
,
{\displaystyle \omega _{\text{n}}={\sqrt {\frac {g}{L}}}={\sqrt {\frac {mgr}{I_{P}}}},}
L
=
g
ω
n
2
=
I
P
m
r
.
{\displaystyle L={\frac {g}{\omega _{\text{n}}^{2}}}={\frac {I_{P}}{mr}}.}
時計の「チクタク」と「ドク」の音を出す秒振り子は 、 左右に振れるのに1秒かかります。これは2秒の周期、つまり 振り子の固有振動数です。この場合、振動の中心までの距離は 次のように計算できます。
π
r
a
d
/
s
{\displaystyle \pi \ \mathrm {rad/s} }
L
{\displaystyle L}
L
=
g
ω
n
2
≈
9.81
m
/
s
2
(
3.14
r
a
d
/
s
)
2
≈
0.99
m
.
{\displaystyle L={\frac {g}{\omega _{\text{n}}^{2}}}\approx {\frac {9.81\ \mathrm {m/s^{2}} }{(3.14\ \mathrm {rad/s} )^{2}}}\approx 0.99\ \mathrm {m} .}
秒振り子の振動中心までの距離は、局所的な重力加速度の異なる値に対応するために調整する必要があることに注意してください。 ケーターの振り子は 、この特性を利用して局所的な重力加速度を測定する複合振り子であり、 重力計 と呼ばれます。
慣性モーメントの測定
車両や飛行機などの複雑なシステムの垂直軸周りの慣性モーメントは、システムを3点から吊り下げて三本線 振り子を 形成することで測定できます。三本線振り子は、垂直の重心軸を中心にねじれ振動するように設計された3本のワイヤーで支えられたプラットフォームです。 [19] 三本線振り子の振動周期から、システムの慣性モーメントが得られます。 [20]
断面2次モーメント
断面二次モーメントは 断面二次モーメント とも呼ばれ、その物理的な意味は質量慣性モーメントとはまったく異なります。これらの計算は、梁や柱の構造設計のために土木工学でよく使用されます。断面積は、x軸の垂直モーメント とy軸の水平モーメントについて計算されます 。
高さ( h )と幅( b )は、円を除いて線形の測定値であり、円は実質的に半幅から導出されます。
I
x
x
{\displaystyle I_{xx}}
I
y
y
{\displaystyle I_{yy}}
r
{\displaystyle r}
断面積モーメントは次のように計算される [21]
四角:
I
x
x
=
I
y
y
=
b
4
12
{\displaystyle I_{xx}=I_{yy}={\frac {b^{4}}{12}}}
長方形: そして;
I
x
x
=
b
h
3
12
{\displaystyle I_{xx}={\frac {bh^{3}}{12}}}
I
y
y
=
h
b
3
12
{\displaystyle I_{yy}={\frac {hb^{3}}{12}}}
三角:
I
x
x
=
b
h
3
36
{\displaystyle I_{xx}={\frac {bh^{3}}{36}}}
円形:
I
x
x
=
I
y
y
=
1
4
π
r
4
=
1
64
π
d
4
{\displaystyle I_{xx}=I_{yy}={\frac {1}{4}}{\pi }r^{4}={\frac {1}{64}}{\pi }d^{4}}
固定面内の動き
点質点
質量と半径が同じ 4 つの物体が、滑ることなく転がりながら平面を滑り降りていきます。
後ろから前へ: 球殻、 固体球、 円筒形のリング、および 固体シリンダー。
各物体がゴールラインに到達する時間は、その慣性モーメントによって決まります。( OGV バージョン )
物体の軸周りの慣性モーメントは、 物体内のすべての粒子についての合計で計算されます。ここで、は 指定された軸に対する垂直距離です。拡張された物体の運動の研究で慣性モーメントがどのように発生するかを確認するには、点質点の剛体集合体を考慮すると便利です。(この式は、慣性モーメントを完全には説明できないことを理解すれば、主軸以外の軸にも使用できます。 [22] )
m
r
2
{\displaystyle mr^{2}}
r
{\displaystyle r}
回転軸上の最も近い点である ピボットポイントからの 距離にある 質量体 の集合体の運動エネルギーを考えてみましょう。これは個々の質量体の運動エネルギーの合計です。 [18] : 516–517 [23] : 1084–1085 [23] : 1296–1300
N
{\displaystyle N}
m
i
{\displaystyle m_{i}}
r
i
{\displaystyle r_{i}}
P
{\displaystyle P}
E
K
=
∑
i
=
1
N
1
2
m
i
v
i
⋅
v
i
=
∑
i
=
1
N
1
2
m
i
(
ω
r
i
)
2
=
1
2
ω
2
∑
i
=
1
N
m
i
r
i
2
.
{\displaystyle E_{\text{K}}=\sum _{i=1}^{N}{\frac {1}{2}}\,m_{i}\mathbf {v} _{i}\cdot \mathbf {v} _{i}=\sum _{i=1}^{N}{\frac {1}{2}}\,m_{i}\left(\omega r_{i}\right)^{2}={\frac {1}{2}}\,\omega ^{2}\sum _{i=1}^{N}m_{i}r_{i}^{2}.}
これは、物体の慣性モーメントが各項の合計であることを示しています 。
m
r
2
{\displaystyle mr^{2}}
I
P
=
∑
i
=
1
N
m
i
r
i
2
.
{\displaystyle I_{P}=\sum _{i=1}^{N}m_{i}r_{i}^{2}.}
したがって、慣性モーメントは、回転軸の周りの粒子の質量と分布を組み合わせた物理的特性です。同じ物体の異なる軸を中心に回転すると、異なる慣性モーメントが生じることに注意してください。
指定された軸の周りを回転する連続体の慣性モーメントは、点粒子が無限にあることを除いて、同じ方法で計算されます。したがって、合計の制限はなくなり、合計は次のように表されます。
I
P
=
∑
i
m
i
r
i
2
{\displaystyle I_{P}=\sum _{i}m_{i}r_{i}^{2}}
別の式では、合計を 積分 に置き換えます。
I
P
=
∭
Q
ρ
(
x
,
y
,
z
)
‖
r
‖
2
d
V
{\displaystyle I_{P}=\iiint _{Q}\rho (x,y,z)\left\|\mathbf {r} \right\|^{2}dV}
ここで、 関数 は 各点 における質量密度を与え 、 は回転軸に垂直で回転軸上の点から 固体内の点まで伸びるベクトルであり、積分は 物体 の体積にわたって評価されます 。平面の慣性モーメントも同様で、質量密度が面積質量密度に置き換えられ、積分はその面積にわたって評価されます。
ρ
{\displaystyle \rho }
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
r
{\displaystyle \mathbf {r} }
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
V
{\displaystyle V}
Q
{\displaystyle Q}
断面二次モーメントに関する注記 :平面内を移動する物体の慣性モーメントと梁の 断面の二次モーメント は混同されることが多い。断面の形状を持つ物体の慣性モーメントは、 断面に垂直な - 軸周りのこの領域の二次モーメントに密度を加重した値である。これは 断面の極モーメント とも呼ばれ、 - 軸 と - 軸周りの二次モーメントの合計である。 [24] 梁 の応力は、荷重に応じて - 軸または- 軸 周りの断面積の二次モーメントを使用して計算される 。
z
{\displaystyle z}
x
{\displaystyle x}
y
{\displaystyle y}
x
{\displaystyle x}
y
{\displaystyle y}
例
細い棒の端に薄い円盤が取り付けられ、棒のもう一方の端にある支点を中心に振動する複合振り子 の慣性 モーメントは、細い棒と薄い円盤のそれぞれの質量中心の慣性モーメントの計算から始まります。 [23]
一定の断面積 と密度を持ち 、長さが 質量中心を通る垂直軸の周りの 細い棒 の慣性モーメントは積分によって決まります。 [23] :1301 軸を棒に 合わせ 、原点を棒の中心に置き、 棒の質量は どこにあるか。
s
{\displaystyle s}
ρ
{\displaystyle \rho }
ℓ
{\displaystyle \ell }
x
{\displaystyle x}
I
C
,
rod
=
∭
Q
ρ
x
2
d
V
=
∫
−
ℓ
2
ℓ
2
ρ
x
2
s
d
x
=
ρ
s
x
3
3
|
−
ℓ
2
ℓ
2
=
ρ
s
3
(
ℓ
3
8
+
ℓ
3
8
)
=
m
ℓ
2
12
,
{\displaystyle I_{C,{\text{rod}}}=\iiint _{Q}\rho \,x^{2}\,dV=\int _{-{\frac {\ell }{2}}}^{\frac {\ell }{2}}\rho \,x^{2}s\,dx=\left.\rho s{\frac {x^{3}}{3}}\right|_{-{\frac {\ell }{2}}}^{\frac {\ell }{2}}={\frac {\rho s}{3}}\left({\frac {\ell ^{3}}{8}}+{\frac {\ell ^{3}}{8}}\right)={\frac {m\ell ^{2}}{12}},}
m
=
ρ
s
ℓ
{\displaystyle m=\rho s\ell }
一定の厚さ 、半径 、密度の 薄い円板 の、中心を通り面に対して垂直な軸( 回転対称 軸 に平行)の周りの慣性モーメントは 積分によって決定されます。 [23] : 1301 [ 検証失敗 ] 軸を円板の軸に 合わせ、体積要素を と定義すると 、 は その 質量になります。
s
{\displaystyle s}
R
{\displaystyle R}
ρ
{\displaystyle \rho }
z
{\displaystyle z}
d
V
=
s
r
d
r
d
θ
{\displaystyle dV=sr\,dr\,d\theta }
I
C
,
disc
=
∭
Q
ρ
r
2
d
V
=
∫
0
2
π
∫
0
R
ρ
r
2
s
r
d
r
d
θ
=
2
π
ρ
s
R
4
4
=
1
2
m
R
2
,
{\displaystyle I_{C,{\text{disc}}}=\iiint _{Q}\rho \,r^{2}\,dV=\int _{0}^{2\pi }\int _{0}^{R}\rho r^{2}sr\,dr\,d\theta =2\pi \rho s{\frac {R^{4}}{4}}={\frac {1}{2}}mR^{2},}
m
=
π
R
2
ρ
s
{\displaystyle m=\pi R^{2}\rho s}
複合振り子の慣性モーメントは、棒と円盤の慣性モーメントをピボット ポイントの周りに加えることで 、次のように得られます。 ここで、は 振り子の長さです。慣性モーメントを質量中心から振り子のピボット ポイントに移動するために、平行軸定理が使用されていることに注意してください。
P
{\displaystyle P}
I
P
=
I
C
,
rod
+
M
rod
(
L
2
)
2
+
I
C
,
disc
+
M
disc
(
L
+
R
)
2
,
{\displaystyle I_{P}=I_{C,{\text{rod}}}+M_{\text{rod}}\left({\frac {L}{2}}\right)^{2}+I_{C,{\text{disc}}}+M_{\text{disc}}(L+R)^{2},}
L
{\displaystyle L}
標準ボディ形状の慣性モーメントの公式 のリストは、 より単純な形状のボディのアセンブリとして複雑なボディの慣性モーメントを取得する方法を提供します。 平行軸定理は、 個々のボディの基準点をアセンブリの基準点にシフトするために使用されます。
もう1つの例として、質量中心を通る軸の周りの一定密度の固体球の慣性モーメントを考えてみましょう。これは、考慮する軸に沿って中心を持つ球を形成できる薄い円板の慣性モーメントを合計することによって決定されます。球の表面が式 [23]で定義されている場合 :1301
x
2
+
y
2
+
z
2
=
R
2
,
{\displaystyle x^{2}+y^{2}+z^{2}=R^{2},}
すると、-軸に沿った 断面における円板の 半径の2乗 は
r
{\displaystyle r}
z
{\displaystyle z}
z
{\displaystyle z}
r
(
z
)
2
=
x
2
+
y
2
=
R
2
−
z
2
.
{\displaystyle r(z)^{2}=x^{2}+y^{2}=R^{2}-z^{2}.}
したがって、球の慣性モーメントは、軸に沿った円板の慣性モーメントの合計です 。
ここで、は 球の質量です。
z
{\displaystyle z}
I
C
,
sphere
=
∫
−
R
R
1
2
π
ρ
r
(
z
)
4
d
z
=
∫
−
R
R
1
2
π
ρ
(
R
2
−
z
2
)
2
d
z
=
1
2
π
ρ
[
R
4
z
−
2
3
R
2
z
3
+
1
5
z
5
]
−
R
R
=
π
ρ
(
1
−
2
3
+
1
5
)
R
5
=
2
5
m
R
2
,
{\displaystyle {\begin{aligned}I_{C,{\text{sphere}}}&=\int _{-R}^{R}{\tfrac {1}{2}}\pi \rho r(z)^{4}\,dz=\int _{-R}^{R}{\tfrac {1}{2}}\pi \rho \left(R^{2}-z^{2}\right)^{2}\,dz\\[1ex]&={\tfrac {1}{2}}\pi \rho \left[R^{4}z-{\tfrac {2}{3}}R^{2}z^{3}+{\tfrac {1}{5}}z^{5}\right]_{-R}^{R}\\[1ex]&=\pi \rho \left(1-{\tfrac {2}{3}}+{\tfrac {1}{5}}\right)R^{5}\\[1ex]&={\tfrac {2}{5}}mR^{2},\end{aligned}}}
m
=
4
3
π
R
3
ρ
{\textstyle m={\frac {4}{3}}\pi R^{3}\rho }
剛体
慣性モーメントが大きいシリンダーは、より多くの位置エネルギーを回転運動エネルギーに変換する必要があるため、より小さな加速度で斜面を転がり落ちます。
機械システム が固定された平面と平行に動くように制約されている場合 、システム内の物体の回転は この平面と平行な軸の周りで発生します。この場合、このシステム内の質量の慣性モーメントは、 極慣性モーメント と呼ばれるスカラーです。極慣性モーメントの定義は、運動量、運動エネルギー、および粒子の剛体の平面運動に関するニュートンの法則を考慮することで得られます。 [15] [18] [25] [26]
k
^
{\displaystyle \mathbf {\hat {k}} }
粒子 のシステム が 剛体に組み立てられる場合、システムの運動量は、基準点 に対する位置 と絶対速度で表すことができます 。
ここで、 はシステムの角速度、 は の速度です 。
n
{\displaystyle n}
P
i
,
i
=
1
,
…
,
n
{\displaystyle P_{i},i=1,\dots ,n}
R
{\displaystyle \mathbf {R} }
v
i
{\displaystyle \mathbf {v} _{i}}
Δ
r
i
=
r
i
−
R
,
v
i
=
ω
×
(
r
i
−
R
)
+
V
=
ω
×
Δ
r
i
+
V
,
{\displaystyle {\begin{aligned}\Delta \mathbf {r} _{i}&=\mathbf {r} _{i}-\mathbf {R} ,\\\mathbf {v} _{i}&={\boldsymbol {\omega }}\times \left(\mathbf {r} _{i}-\mathbf {R} \right)+\mathbf {V} ={\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}+\mathbf {V} ,\end{aligned}}}
ω
{\displaystyle {\boldsymbol {\omega }}}
V
{\displaystyle \mathbf {V} }
R
{\displaystyle \mathbf {R} }
平面運動の場合、角速度ベクトルは 運動平面に垂直な単位ベクトルに沿って向きます。 基準点から 点までの単位ベクトル と単位ベクトルを導入すると 、
k
{\displaystyle \mathbf {k} }
e
i
{\displaystyle \mathbf {e} _{i}}
R
{\displaystyle \mathbf {R} }
r
i
{\displaystyle \mathbf {r} _{i}}
t
^
i
=
k
^
×
e
^
i
{\displaystyle \mathbf {\hat {t}} _{i}=\mathbf {\hat {k}} \times \mathbf {\hat {e}} _{i}}
e
^
i
=
Δ
r
i
Δ
r
i
,
k
^
=
ω
ω
,
t
^
i
=
k
^
×
e
^
i
,
v
i
=
ω
×
Δ
r
i
+
V
=
ω
k
^
×
Δ
r
i
e
^
i
+
V
=
ω
Δ
r
i
t
^
i
+
V
{\displaystyle {\begin{aligned}\mathbf {\hat {e}} _{i}&={\frac {\Delta \mathbf {r} _{i}}{\Delta r_{i}}},\quad \mathbf {\hat {k}} ={\frac {\boldsymbol {\omega }}{\omega }},\quad \mathbf {\hat {t}} _{i}=\mathbf {\hat {k}} \times \mathbf {\hat {e}} _{i},\\\mathbf {v} _{i}&={\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}+\mathbf {V} =\omega \mathbf {\hat {k}} \times \Delta r_{i}\mathbf {\hat {e}} _{i}+\mathbf {V} =\omega \,\Delta r_{i}\mathbf {\hat {t}} _{i}+\mathbf {V} \end{aligned}}}
これは、平面内を移動する粒子の剛体システムの相対位置ベクトルと速度ベクトルを定義します。
外積に関する注意 : 物体が地面と平行に動く場合、物体のすべての点の軌道はこの地面と平行な平面上にあります。つまり、物体が回転する場合は、必ずこの平面に垂直な軸を中心に回転します。平面運動は、回転軸が点として現れるように、この地面に投影されて表現されることがよくあります。この場合、物体の角速度と角加速度はスカラーであり、回転軸に沿ったベクトルであるという事実は無視されます。これは通常、このトピックの導入で好まれます。しかし、慣性モーメントの場合、質量と形状の組み合わせは外積の幾何学的特性から恩恵を受けます。このため、平面運動に関するこのセクションでは、物体の角速度と角加速度は地面に垂直なベクトルであり、外積の演算は空間剛体運動の研究で使用されるものと同じです。
角運動量
剛体粒子系の平面運動の角運動量ベクトルは [15] [18]で与えられる。
L
=
∑
i
=
1
n
m
i
Δ
r
i
×
v
i
=
∑
i
=
1
n
m
i
Δ
r
i
e
^
i
×
(
ω
Δ
r
i
t
^
i
+
V
)
=
(
∑
i
=
1
n
m
i
Δ
r
i
2
)
ω
k
^
+
(
∑
i
=
1
n
m
i
Δ
r
i
e
^
i
)
×
V
.
{\displaystyle {\begin{aligned}\mathbf {L} &=\sum _{i=1}^{n}m_{i}\Delta \mathbf {r} _{i}\times \mathbf {v} _{i}\\&=\sum _{i=1}^{n}m_{i}\,\Delta r_{i}\mathbf {\hat {e}} _{i}\times \left(\omega \,\Delta r_{i}\mathbf {\hat {t}} _{i}+\mathbf {V} \right)\\&=\left(\sum _{i=1}^{n}m_{i}\,\Delta r_{i}^{2}\right)\omega \mathbf {\hat {k}} +\left(\sum _{i=1}^{n}m_{i}\,\Delta r_{i}\mathbf {\hat {e}} _{i}\right)\times \mathbf {V} .\end{aligned}}}
質量の中心 を 基準点として
C
{\displaystyle \mathbf {C} }
Δ
r
i
e
^
i
=
r
i
−
C
,
∑
i
=
1
n
m
i
Δ
r
i
e
^
i
=
0
,
{\displaystyle {\begin{aligned}\Delta r_{i}\mathbf {\hat {e}} _{i}&=\mathbf {r} _{i}-\mathbf {C} ,\\\sum _{i=1}^{n}m_{i}\,\Delta r_{i}\mathbf {\hat {e}} _{i}&=0,\end{aligned}}}
そして質量中心に対する慣性モーメントを 次のように
定義する。
I
C
{\displaystyle I_{\mathbf {C} }}
I
C
=
∑
i
m
i
Δ
r
i
2
,
{\displaystyle I_{\mathbf {C} }=\sum _{i}m_{i}\,\Delta r_{i}^{2},}
すると角運動量の式は [23]のように簡略化される :1028
L
=
I
C
ω
k
^
.
{\displaystyle \mathbf {L} =I_{\mathbf {C} }\omega \mathbf {\hat {k}} .}
剛体システムの動きに垂直で質量の中心を通る軸の周り の慣性モーメントは、 極慣性モーメント として知られています。具体的には、軸(または極)からの直交距離に関する
質量の 2 番目のモーメント です。
I
C
{\displaystyle I_{\mathbf {C} }}
角運動量が一定であれば、慣性モーメントが減少すると角速度が増加します。フィギュア スケーターは腕を引き締めることで慣性モーメントを変更できます。したがって、腕を伸ばしたスケーターが達成した角速度は、腕を引き締めると慣性モーメントが減少するため、角速度が大きくなります。ただし、フィギュア スケーターは剛体ではありません。
運動エネルギー
この 1906 年型回転式剪断機は、2 つのフライホイールの慣性モーメントを利用して運動エネルギーを蓄え、そのエネルギーを放出して金属材料を切断します (国際技術図書館、1906 年)。
平面内を運動する剛体粒子系の運動エネルギーは [15] [18]で与えられる。
E
K
=
1
2
∑
i
=
1
n
m
i
v
i
⋅
v
i
,
=
1
2
∑
i
=
1
n
m
i
(
ω
Δ
r
i
t
^
i
+
V
)
⋅
(
ω
Δ
r
i
t
^
i
+
V
)
,
=
1
2
ω
2
(
∑
i
=
1
n
m
i
Δ
r
i
2
t
^
i
⋅
t
^
i
)
+
ω
V
⋅
(
∑
i
=
1
n
m
i
Δ
r
i
t
^
i
)
+
1
2
(
∑
i
=
1
n
m
i
)
V
⋅
V
.
{\displaystyle {\begin{aligned}E_{\text{K}}&={\frac {1}{2}}\sum _{i=1}^{n}m_{i}\mathbf {v} _{i}\cdot \mathbf {v} _{i},\\&={\frac {1}{2}}\sum _{i=1}^{n}m_{i}\left(\omega \,\Delta r_{i}\mathbf {\hat {t}} _{i}+\mathbf {V} \right)\cdot \left(\omega \,\Delta r_{i}\mathbf {\hat {t}} _{i}+\mathbf {V} \right),\\&={\frac {1}{2}}\omega ^{2}\left(\sum _{i=1}^{n}m_{i}\,\Delta r_{i}^{2}\mathbf {\hat {t}} _{i}\cdot \mathbf {\hat {t}} _{i}\right)+\omega \mathbf {V} \cdot \left(\sum _{i=1}^{n}m_{i}\,\Delta r_{i}\mathbf {\hat {t}} _{i}\right)+{\frac {1}{2}}\left(\sum _{i=1}^{n}m_{i}\right)\mathbf {V} \cdot \mathbf {V} .\end{aligned}}}
参照点をシステムの質量中心とする と、第2項はゼロになり、慣性モーメントを導入する と、運動エネルギーは [23]で与えられる :1084
C
{\displaystyle \mathbf {C} }
I
C
{\displaystyle I_{\mathbf {C} }}
E
K
=
1
2
I
C
ω
2
+
1
2
M
V
⋅
V
.
{\displaystyle E_{\text{K}}={\frac {1}{2}}I_{\mathbf {C} }\omega ^{2}+{\frac {1}{2}}M\mathbf {V} \cdot \mathbf {V} .}
慣性モーメントとは、 物体の
慣性極モーメントのこと です。
I
C
{\displaystyle I_{\mathbf {C} }}
ニュートンの法則
1920 年代の John Deere トラクター。エンジンに スポーク フライホイールが付いています。フライホイールの大きな慣性モーメントにより、トラクターの動作がスムーズになります。
ニュートンの剛体 粒子系に対する法則は、 基準点における 合力 とトルクで表すことができ、 [15] [18] となる。
ここで、は 各粒子の軌道を表す。
n
{\displaystyle n}
P
i
,
i
=
1
,
…
,
n
{\displaystyle P_{i},i=1,\dots ,n}
R
{\displaystyle \mathbf {R} }
F
=
∑
i
=
1
n
m
i
A
i
,
τ
=
∑
i
=
1
n
Δ
r
i
×
m
i
A
i
,
{\displaystyle {\begin{aligned}\mathbf {F} &=\sum _{i=1}^{n}m_{i}\mathbf {A} _{i},\\{\boldsymbol {\tau }}&=\sum _{i=1}^{n}\Delta \mathbf {r} _{i}\times m_{i}\mathbf {A} _{i},\end{aligned}}}
r
i
{\displaystyle \mathbf {r} _{i}}
剛体の 運動学では 、基準粒子の 位置と 加速度、 および剛体粒子 系の
角速度ベクトル と角加速度ベクトルに基づいて、粒子の加速度の式が次のように表されます。
P
i
{\displaystyle P_{i}}
R
{\displaystyle \mathbf {R} }
A
{\displaystyle \mathbf {A} }
ω
{\displaystyle {\boldsymbol {\omega }}}
α
{\displaystyle {\boldsymbol {\alpha }}}
A
i
=
α
×
Δ
r
i
+
ω
×
ω
×
Δ
r
i
+
A
.
{\displaystyle \mathbf {A} _{i}={\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i}+{\boldsymbol {\omega }}\times {\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}+\mathbf {A} .}
平面運動に制限されているシステムの場合、角速度ベクトルと角加速度ベクトルは運動面に対して垂直に向いており、この加速度方程式は単純化されます。この場合、加速度ベクトルは、 基準点から 点への 単位ベクトル と単位ベクトルを導入することで単純化できます 。
k
^
{\displaystyle \mathbf {\hat {k}} }
e
^
i
{\displaystyle \mathbf {\hat {e}} _{i}}
R
{\displaystyle \mathbf {R} }
r
i
{\displaystyle \mathbf {r} _{i}}
t
^
i
=
k
^
×
e
^
i
{\displaystyle \mathbf {\hat {t}} _{i}=\mathbf {\hat {k}} \times \mathbf {\hat {e}} _{i}}
A
i
=
α
k
^
×
Δ
r
i
e
^
i
−
ω
k
^
×
ω
k
^
×
Δ
r
i
e
^
i
+
A
=
α
Δ
r
i
t
^
i
−
ω
2
Δ
r
i
e
^
i
+
A
.
{\displaystyle {\begin{aligned}\mathbf {A} _{i}&=\alpha \mathbf {\hat {k}} \times \Delta r_{i}\mathbf {\hat {e}} _{i}-\omega \mathbf {\hat {k}} \times \omega \mathbf {\hat {k}} \times \Delta r_{i}\mathbf {\hat {e}} _{i}+\mathbf {A} \\&=\alpha \Delta r_{i}\mathbf {\hat {t}} _{i}-\omega ^{2}\Delta r_{i}\mathbf {\hat {e}} _{i}+\mathbf {A} .\end{aligned}}}
これにより、システムにかかるトルクは次のようになります。
τ
=
∑
i
=
1
n
m
i
Δ
r
i
e
^
i
×
(
α
Δ
r
i
t
^
i
−
ω
2
Δ
r
i
e
^
i
+
A
)
=
(
∑
i
=
1
n
m
i
Δ
r
i
2
)
α
k
^
+
(
∑
i
=
1
n
m
i
Δ
r
i
e
^
i
)
×
A
,
{\displaystyle {\begin{aligned}{\boldsymbol {\tau }}&=\sum _{i=1}^{n}m_{i}\,\Delta r_{i}\mathbf {\hat {e}} _{i}\times \left(\alpha \Delta r_{i}\mathbf {\hat {t}} _{i}-\omega ^{2}\Delta r_{i}\mathbf {\hat {e}} _{i}+\mathbf {A} \right)\\&=\left(\sum _{i=1}^{n}m_{i}\,\Delta r_{i}^{2}\right)\alpha \mathbf {\hat {k}} +\left(\sum _{i=1}^{n}m_{i}\,\Delta r_{i}\mathbf {\hat {e}} _{i}\right)\times \mathbf {A} ,\end{aligned}}}
ここで 、 は すべての粒子の平面に垂直な単位ベクトルです 。
e
^
i
×
e
^
i
=
0
{\displaystyle \mathbf {\hat {e}} _{i}\times \mathbf {\hat {e}} _{i}=\mathbf {0} }
e
^
i
×
t
^
i
=
k
^
{\displaystyle \mathbf {\hat {e}} _{i}\times \mathbf {\hat {t}} _{i}=\mathbf {\hat {k}} }
P
i
{\displaystyle P_{i}}
質量中心を 基準 点として、質量中心に対する慣性モーメントを定義すると 、結果として生じるトルクの式は次のように簡略化される [23] :1029
C
{\displaystyle \mathbf {C} }
I
C
{\displaystyle I_{\mathbf {C} }}
τ
=
I
C
α
k
^
.
{\displaystyle {\boldsymbol {\tau }}=I_{\mathbf {C} }\alpha \mathbf {\hat {k}} .}
剛体の空間運動と慣性行列
スカラー慣性モーメントは、粒子系が3次元空間で運動する剛体に組み立てられるときに、行列の要素として現れる。この慣性行列は、粒子の剛体系の角運動量、運動エネルギー、および結果として生じるトルクの計算に現れる。 [4] [5] [6] [7] [27]
粒子 のシステム が、 固定された基準フレームに対する 速度 で 座標 に配置されているとします。(移動する可能性のある) 基準点 の場合 、相対位置は で
、(絶対) 速度は です
。
ここで はシステムの角速度、 は の速度です 。
n
{\displaystyle n}
P
i
,
i
=
1
,
…
,
n
{\displaystyle P_{i},i=1,\dots ,n}
r
i
{\displaystyle \mathbf {r} _{i}}
v
i
{\displaystyle \mathbf {v} _{i}}
R
{\displaystyle \mathbf {R} }
Δ
r
i
=
r
i
−
R
{\displaystyle \Delta \mathbf {r} _{i}=\mathbf {r} _{i}-\mathbf {R} }
v
i
=
ω
×
Δ
r
i
+
V
R
{\displaystyle \mathbf {v} _{i}={\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}+\mathbf {V} _{\mathbf {R} }}
ω
{\displaystyle {\boldsymbol {\omega }}}
V
R
{\displaystyle \mathbf {V_{R}} }
R
{\displaystyle \mathbf {R} }
角運動量
外積は、最初のオペランドと演算子をの要素から構成される 歪対称行列 に結合することによって、 行列の乗算として同等に記述できること に注意します 。
[
b
]
{\displaystyle \left[\mathbf {b} \right]}
b
=
(
b
x
,
b
y
,
b
z
)
{\displaystyle \mathbf {b} =(b_{x},b_{y},b_{z})}
b
×
y
≡
[
b
]
y
[
b
]
≡
[
0
−
b
z
b
y
b
z
0
−
b
x
−
b
y
b
x
0
]
.
{\displaystyle {\begin{aligned}\mathbf {b} \times \mathbf {y} &\equiv \left[\mathbf {b} \right]\mathbf {y} \\\left[\mathbf {b} \right]&\equiv {\begin{bmatrix}0&-b_{z}&b_{y}\\b_{z}&0&-b_{x}\\-b_{y}&b_{x}&0\end{bmatrix}}.\end{aligned}}}
慣性行列は角運動量を考慮して構築され、 物体の基準点は質量中心に選ばれます 。 [4] [7] ここで、 ( )
を含む項の合計は 質量中心の 定義によりゼロになります 。
R
{\displaystyle \mathbf {R} }
C
{\displaystyle \mathbf {C} }
L
=
∑
i
=
1
n
m
i
Δ
r
i
×
v
i
=
∑
i
=
1
n
m
i
Δ
r
i
×
(
ω
×
Δ
r
i
+
V
R
)
=
(
−
∑
i
=
1
n
m
i
Δ
r
i
×
(
Δ
r
i
×
ω
)
)
+
(
∑
i
=
1
n
m
i
Δ
r
i
×
V
R
)
,
{\displaystyle {\begin{aligned}\mathbf {L} &=\sum _{i=1}^{n}m_{i}\,\Delta \mathbf {r} _{i}\times \mathbf {v} _{i}\\&=\sum _{i=1}^{n}m_{i}\,\Delta \mathbf {r} _{i}\times \left({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}+\mathbf {V} _{\mathbf {R} }\right)\\&=\left(-\sum _{i=1}^{n}m_{i}\,\Delta \mathbf {r} _{i}\times \left(\Delta \mathbf {r} _{i}\times {\boldsymbol {\omega }}\right)\right)+\left(\sum _{i=1}^{n}m_{i}\,\Delta \mathbf {r} _{i}\times \mathbf {V} _{\mathbf {R} }\right),\end{aligned}}}
V
R
{\displaystyle \mathbf {V_{R}} }
=
C
{\displaystyle =\mathbf {C} }
次に、相対位置ベクトル から得られる 歪対称行列 を使用して を定義できます。
ここで 、 によって定義される は
、質量中心 を基準として測定された粒子の剛体系の対称慣性行列です 。
[
Δ
r
i
]
{\displaystyle [\Delta \mathbf {r} _{i}]}
Δ
r
i
=
r
i
−
C
{\displaystyle \Delta \mathbf {r} _{i}=\mathbf {r} _{i}-\mathbf {C} }
L
=
(
−
∑
i
=
1
n
m
i
[
Δ
r
i
]
2
)
ω
=
I
C
ω
,
{\displaystyle \mathbf {L} =\left(-\sum _{i=1}^{n}m_{i}\left[\Delta \mathbf {r} _{i}\right]^{2}\right){\boldsymbol {\omega }}=\mathbf {I} _{\mathbf {C} }{\boldsymbol {\omega }},}
I
C
{\displaystyle \mathbf {I_{C}} }
I
C
=
−
∑
i
=
1
n
m
i
[
Δ
r
i
]
2
,
{\displaystyle \mathbf {I} _{\mathbf {C} }=-\sum _{i=1}^{n}m_{i}\left[\Delta \mathbf {r} _{i}\right]^{2},}
C
{\displaystyle \mathbf {C} }
運動エネルギー
剛体粒子系の運動エネルギーは、系の 質量中心 と慣性モーメントの行列で表すことができます。 粒子系を 速度の 座標に配置すると 、運動エネルギーは [4] [7] となり、
ここで は 質量中心に対する粒子の位置ベクトルです。
n
{\displaystyle n}
P
i
,
i
=
1
,
…
,
n
{\displaystyle P_{i},i=1,\dots ,n}
r
i
{\displaystyle \mathbf {r} _{i}}
v
i
{\displaystyle \mathbf {v} _{i}}
E
K
=
1
2
∑
i
=
1
n
m
i
v
i
⋅
v
i
=
1
2
∑
i
=
1
n
m
i
(
ω
×
Δ
r
i
+
V
C
)
⋅
(
ω
×
Δ
r
i
+
V
C
)
,
{\displaystyle E_{\text{K}}={\frac {1}{2}}\sum _{i=1}^{n}m_{i}\mathbf {v} _{i}\cdot \mathbf {v} _{i}={\frac {1}{2}}\sum _{i=1}^{n}m_{i}\left({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}+\mathbf {V} _{\mathbf {C} }\right)\cdot \left({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}+\mathbf {V} _{\mathbf {C} }\right),}
Δ
r
i
=
r
i
−
C
{\displaystyle \Delta \mathbf {r} _{i}=\mathbf {r} _{i}-\mathbf {C} }
この式を展開すると3つの項が得られる。
E
K
=
1
2
(
∑
i
=
1
n
m
i
(
ω
×
Δ
r
i
)
⋅
(
ω
×
Δ
r
i
)
)
+
(
∑
i
=
1
n
m
i
V
C
⋅
(
ω
×
Δ
r
i
)
)
+
1
2
(
∑
i
=
1
n
m
i
V
C
⋅
V
C
)
.
{\displaystyle E_{\text{K}}={\frac {1}{2}}\left(\sum _{i=1}^{n}m_{i}\left({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}\right)\cdot \left({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}\right)\right)+\left(\sum _{i=1}^{n}m_{i}\mathbf {V} _{\mathbf {C} }\cdot \left({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}\right)\right)+{\frac {1}{2}}\left(\sum _{i=1}^{n}m_{i}\mathbf {V} _{\mathbf {C} }\cdot \mathbf {V} _{\mathbf {C} }\right).}
質量の中心は で定義されるので
、この式の2番目の項はゼロです。歪対称行列を導入すると 、運動エネルギーは
∑
i
=
1
n
m
i
Δ
r
i
=
0
{\displaystyle \sum _{i=1}^{n}m_{i}\Delta \mathbf {r} _{i}=0}
[
Δ
r
i
]
{\displaystyle [\Delta \mathbf {r} _{i}]}
E
K
=
1
2
(
∑
i
=
1
n
m
i
(
[
Δ
r
i
]
ω
)
⋅
(
[
Δ
r
i
]
ω
)
)
+
1
2
(
∑
i
=
1
n
m
i
)
V
C
⋅
V
C
=
1
2
(
∑
i
=
1
n
m
i
(
ω
T
[
Δ
r
i
]
T
[
Δ
r
i
]
ω
)
)
+
1
2
(
∑
i
=
1
n
m
i
)
V
C
⋅
V
C
=
1
2
ω
⋅
(
−
∑
i
=
1
n
m
i
[
Δ
r
i
]
2
)
ω
+
1
2
(
∑
i
=
1
n
m
i
)
V
C
⋅
V
C
.
{\displaystyle {\begin{aligned}E_{\text{K}}&={\frac {1}{2}}\left(\sum _{i=1}^{n}m_{i}\left(\left[\Delta \mathbf {r} _{i}\right]{\boldsymbol {\omega }}\right)\cdot \left(\left[\Delta \mathbf {r} _{i}\right]{\boldsymbol {\omega }}\right)\right)+{\frac {1}{2}}\left(\sum _{i=1}^{n}m_{i}\right)\mathbf {V} _{\mathbf {C} }\cdot \mathbf {V} _{\mathbf {C} }\\&={\frac {1}{2}}\left(\sum _{i=1}^{n}m_{i}\left({\boldsymbol {\omega }}^{\mathsf {T}}\left[\Delta \mathbf {r} _{i}\right]^{\mathsf {T}}\left[\Delta \mathbf {r} _{i}\right]{\boldsymbol {\omega }}\right)\right)+{\frac {1}{2}}\left(\sum _{i=1}^{n}m_{i}\right)\mathbf {V} _{\mathbf {C} }\cdot \mathbf {V} _{\mathbf {C} }\\&={\frac {1}{2}}{\boldsymbol {\omega }}\cdot \left(-\sum _{i=1}^{n}m_{i}\left[\Delta \mathbf {r} _{i}\right]^{2}\right){\boldsymbol {\omega }}+{\frac {1}{2}}\left(\sum _{i=1}^{n}m_{i}\right)\mathbf {V} _{\mathbf {C} }\cdot \mathbf {V} _{\mathbf {C} }.\end{aligned}}}
したがって、粒子の剛体系の運動エネルギーは、
が質量中心に対する慣性行列であり、 が
総 質量である、 によって与えられます。
E
K
=
1
2
ω
⋅
I
C
ω
+
1
2
M
V
C
2
.
{\displaystyle E_{\text{K}}={\frac {1}{2}}{\boldsymbol {\omega }}\cdot \mathbf {I} _{\mathbf {C} }{\boldsymbol {\omega }}+{\frac {1}{2}}M\mathbf {V} _{\mathbf {C} }^{2}.}
I
C
{\displaystyle \mathbf {I_{C}} }
M
{\displaystyle M}
結果トルク
慣性行列は、ニュートンの第2法則を剛体粒子集合体に適用する際に現れる。このシステムで生じるトルクは、 [4] [7]
で、
ここで は 粒子の加速度である 。 剛体の 運動学 では、基準点の位置と加速度 、および剛体システムの角速度ベクトル と角加速度ベクトル によって、粒子の加速度の式 が次のように得られる。
τ
=
∑
i
=
1
n
(
r
i
−
R
)
×
m
i
a
i
,
{\displaystyle {\boldsymbol {\tau }}=\sum _{i=1}^{n}\left(\mathbf {r_{i}} -\mathbf {R} \right)\times m_{i}\mathbf {a} _{i},}
a
i
{\displaystyle \mathbf {a} _{i}}
P
i
{\displaystyle P_{i}}
P
i
{\displaystyle P_{i}}
R
{\displaystyle \mathbf {R} }
A
R
{\displaystyle \mathbf {A} _{\mathbf {R} }}
ω
{\displaystyle {\boldsymbol {\omega }}}
α
{\displaystyle {\boldsymbol {\alpha }}}
a
i
=
α
×
(
r
i
−
R
)
+
ω
×
(
ω
×
(
r
i
−
R
)
)
+
A
R
.
{\displaystyle \mathbf {a} _{i}={\boldsymbol {\alpha }}\times \left(\mathbf {r} _{i}-\mathbf {R} \right)+{\boldsymbol {\omega }}\times \left({\boldsymbol {\omega }}\times \left(\mathbf {r} _{i}-\mathbf {R} \right)\right)+\mathbf {A} _{\mathbf {R} }.}
質量中心を 基準点として、 交差積を表す歪対称行列を導入する と、
C
{\displaystyle \mathbf {C} }
[
Δ
r
i
]
=
[
r
i
−
C
]
{\displaystyle \left[\Delta \mathbf {r} _{i}\right]=\left[\mathbf {r} _{i}-\mathbf {C} \right]}
(
r
i
−
C
)
×
{\displaystyle (\mathbf {r} _{i}-\mathbf {C} )\times }
τ
=
(
−
∑
i
=
1
n
m
i
[
Δ
r
i
]
2
)
α
+
ω
×
(
−
∑
i
=
1
n
m
i
[
Δ
r
i
]
2
)
ω
{\displaystyle {\boldsymbol {\tau }}=\left(-\sum _{i=1}^{n}m_{i}\left[\Delta \mathbf {r} _{i}\right]^{2}\right){\boldsymbol {\alpha }}+{\boldsymbol {\omega }}\times \left(-\sum _{i=1}^{n}m_{i}\left[\Delta \mathbf {r} _{i}\right]^{2}\right){\boldsymbol {\omega }}}
計算では、
以下の証明に示すように、
三重 クロス積の ヤコビ恒等式
から得られる恒等式を使用します。
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
+
ω
×
(
(
ω
×
Δ
r
i
)
×
Δ
r
i
)
=
0
,
{\displaystyle \Delta \mathbf {r} _{i}\times \left({\boldsymbol {\omega }}\times \left({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}\right)\right)+{\boldsymbol {\omega }}\times \left(\left({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}\right)\times \Delta \mathbf {r} _{i}\right)=0,}
証拠
τ
=
∑
i
=
1
n
(
r
i
−
R
)
×
(
m
i
a
i
)
=
∑
i
=
1
n
Δ
r
i
×
(
m
i
a
i
)
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
a
i
]
…
cross-product scalar multiplication
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
a
tangential
,
i
+
a
centripetal
,
i
+
A
R
)
]
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
a
tangential
,
i
+
a
centripetal
,
i
+
0
)
]
{\displaystyle {\begin{aligned}{\boldsymbol {\tau }}&=\sum _{i=1}^{n}(\mathbf {r_{i}} -\mathbf {R} )\times (m_{i}\mathbf {a} _{i})\\&=\sum _{i=1}^{n}\Delta \mathbf {r} _{i}\times (m_{i}\mathbf {a} _{i})\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times \mathbf {a} _{i}]\;\ldots {\text{ cross-product scalar multiplication}}\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times (\mathbf {a} _{{\text{tangential}},i}+\mathbf {a} _{{\text{centripetal}},i}+\mathbf {A} _{\mathbf {R} })]\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times (\mathbf {a} _{{\text{tangential}},i}+\mathbf {a} _{{\text{centripetal}},i}+0)]\\\end{aligned}}}
最後の文では、 は静止しているか、加速されていない一定の速度で動いているか、固定された(世界)座標参照系の原点が質量の中心に置かれている ためです 。そして、積を合計に分配すると、次の式が得られます。
A
R
=
0
{\displaystyle \mathbf {A} _{\mathbf {R} }=0}
R
{\displaystyle \mathbf {R} }
C
{\displaystyle \mathbf {C} }
τ
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
a
tangential
,
i
+
Δ
r
i
×
a
centripetal
,
i
]
τ
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
Δ
r
i
×
(
ω
×
v
tangential
,
i
)
]
τ
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
]
{\displaystyle {\begin{aligned}{\boldsymbol {\tau }}&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times \mathbf {a} _{{\text{tangential}},i}+\Delta \mathbf {r} _{i}\times \mathbf {a} _{{\text{centripetal}},i}]\\{\boldsymbol {\tau }}&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times \mathbf {v} _{{\text{tangential}},i})]\\{\boldsymbol {\tau }}&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))]\end{aligned}}}
そして、最後の項には
次の ヤコビ恒等式が使用されます。
0
=
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
+
ω
×
(
(
ω
×
Δ
r
i
)
×
Δ
r
i
)
+
(
ω
×
Δ
r
i
)
×
(
Δ
r
i
×
ω
)
=
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
+
ω
×
(
(
ω
×
Δ
r
i
)
×
Δ
r
i
)
+
(
ω
×
Δ
r
i
)
×
−
(
ω
×
Δ
r
i
)
…
cross-product anticommutativity
=
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
+
ω
×
(
(
ω
×
Δ
r
i
)
×
Δ
r
i
)
+
−
[
(
ω
×
Δ
r
i
)
×
(
ω
×
Δ
r
i
)
]
…
cross-product scalar multiplication
=
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
+
ω
×
(
(
ω
×
Δ
r
i
)
×
Δ
r
i
)
+
−
[
0
]
…
self cross-product
0
=
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
+
ω
×
(
(
ω
×
Δ
r
i
)
×
Δ
r
i
)
{\displaystyle {\begin{aligned}0&=\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))+{\boldsymbol {\omega }}\times (({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\times \Delta \mathbf {r} _{i})+({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\times (\Delta \mathbf {r} _{i}\times {\boldsymbol {\omega }})\\&=\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))+{\boldsymbol {\omega }}\times (({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\times \Delta \mathbf {r} _{i})+({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\times -({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\;\ldots {\text{ cross-product anticommutativity}}\\&=\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))+{\boldsymbol {\omega }}\times (({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\times \Delta \mathbf {r} _{i})+-[({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})]\;\ldots {\text{ cross-product scalar multiplication}}\\&=\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))+{\boldsymbol {\omega }}\times (({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\times \Delta \mathbf {r} _{i})+-[0]\;\ldots {\text{ self cross-product}}\\0&=\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))+{\boldsymbol {\omega }}\times (({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\times \Delta \mathbf {r} _{i})\end{aligned}}}
ヤコビ恒等式 を適用した結果は 、次のように続きます。
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
=
−
[
ω
×
(
(
ω
×
Δ
r
i
)
×
Δ
r
i
)
]
=
−
[
(
ω
×
Δ
r
i
)
(
ω
⋅
Δ
r
i
)
−
Δ
r
i
(
ω
⋅
(
ω
×
Δ
r
i
)
)
]
…
vector triple product
=
−
[
(
ω
×
Δ
r
i
)
(
ω
⋅
Δ
r
i
)
−
Δ
r
i
(
Δ
r
i
⋅
(
ω
×
ω
)
)
]
…
scalar triple product
=
−
[
(
ω
×
Δ
r
i
)
(
ω
⋅
Δ
r
i
)
−
Δ
r
i
(
Δ
r
i
⋅
(
0
)
)
]
…
self cross-product
=
−
[
(
ω
×
Δ
r
i
)
(
ω
⋅
Δ
r
i
)
]
=
−
[
ω
×
(
Δ
r
i
(
ω
⋅
Δ
r
i
)
)
]
…
cross-product scalar multiplication
=
ω
×
−
(
Δ
r
i
(
ω
⋅
Δ
r
i
)
)
…
cross-product scalar multiplication
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
=
ω
×
−
(
Δ
r
i
(
Δ
r
i
⋅
ω
)
)
…
dot-product commutativity
{\displaystyle {\begin{aligned}\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))&=-[{\boldsymbol {\omega }}\times (({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\times \Delta \mathbf {r} _{i})]\\&=-[({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})({\boldsymbol {\omega }}\cdot \Delta \mathbf {r} _{i})-\Delta \mathbf {r} _{i}({\boldsymbol {\omega }}\cdot ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))]\;\ldots {\text{ vector triple product}}\\&=-[({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})({\boldsymbol {\omega }}\cdot \Delta \mathbf {r} _{i})-\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot ({\boldsymbol {\omega }}\times {\boldsymbol {\omega }}))]\;\ldots {\text{ scalar triple product}}\\&=-[({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})({\boldsymbol {\omega }}\cdot \Delta \mathbf {r} _{i})-\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot (0))]\;\ldots {\text{ self cross-product}}\\&=-[({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})({\boldsymbol {\omega }}\cdot \Delta \mathbf {r} _{i})]\\&=-[{\boldsymbol {\omega }}\times (\Delta \mathbf {r} _{i}({\boldsymbol {\omega }}\cdot \Delta \mathbf {r} _{i}))]\;\ldots {\text{ cross-product scalar multiplication}}\\&={\boldsymbol {\omega }}\times -(\Delta \mathbf {r} _{i}({\boldsymbol {\omega }}\cdot \Delta \mathbf {r} _{i}))\;\ldots {\text{ cross-product scalar multiplication}}\\\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))&={\boldsymbol {\omega }}\times -(\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot {\boldsymbol {\omega }}))\;\ldots {\text{ dot-product commutativity}}\\\end{aligned}}}
最終結果は次のようにメインの証明に置き換えることができます。
τ
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
Δ
r
i
×
(
ω
×
(
ω
×
Δ
r
i
)
)
]
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
ω
×
−
(
Δ
r
i
(
Δ
r
i
⋅
ω
)
)
]
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
ω
×
{
0
−
Δ
r
i
(
Δ
r
i
⋅
ω
)
}
]
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
ω
×
{
[
ω
(
Δ
r
i
⋅
Δ
r
i
)
−
ω
(
Δ
r
i
⋅
Δ
r
i
)
]
−
Δ
r
i
(
Δ
r
i
⋅
ω
)
}
]
…
ω
(
Δ
r
i
⋅
Δ
r
i
)
−
ω
(
Δ
r
i
⋅
Δ
r
i
)
=
0
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
ω
×
{
[
ω
(
Δ
r
i
⋅
Δ
r
i
)
−
Δ
r
i
(
Δ
r
i
⋅
ω
)
]
−
ω
(
Δ
r
i
⋅
Δ
r
i
)
}
]
…
addition associativity
{\displaystyle {\begin{aligned}{\boldsymbol {\tau }}&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i}))]\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+{\boldsymbol {\omega }}\times -(\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot {\boldsymbol {\omega }}))]\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+{\boldsymbol {\omega }}\times \{0-\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot {\boldsymbol {\omega }})\}]\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+{\boldsymbol {\omega }}\times \{[{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})-{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})]-\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot {\boldsymbol {\omega }})\}]\;\ldots \;{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})-{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})=0\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+{\boldsymbol {\omega }}\times \{[{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})-\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot {\boldsymbol {\omega }})]-{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})\}]\;\ldots {\text{ addition associativity}}\\\end{aligned}}}
τ
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
ω
×
{
ω
(
Δ
r
i
⋅
Δ
r
i
)
−
Δ
r
i
(
Δ
r
i
⋅
ω
)
}
−
ω
×
ω
(
Δ
r
i
⋅
Δ
r
i
)
]
…
cross-product distributivity over addition
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
ω
×
{
ω
(
Δ
r
i
⋅
Δ
r
i
)
−
Δ
r
i
(
Δ
r
i
⋅
ω
)
}
−
(
Δ
r
i
⋅
Δ
r
i
)
(
ω
×
ω
)
]
…
cross-product scalar multiplication
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
ω
×
{
ω
(
Δ
r
i
⋅
Δ
r
i
)
−
Δ
r
i
(
Δ
r
i
⋅
ω
)
}
−
(
Δ
r
i
⋅
Δ
r
i
)
(
0
)
]
…
self cross-product
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
ω
×
{
ω
(
Δ
r
i
⋅
Δ
r
i
)
−
Δ
r
i
(
Δ
r
i
⋅
ω
)
}
]
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
α
×
Δ
r
i
)
+
ω
×
{
Δ
r
i
×
(
ω
×
Δ
r
i
)
}
]
…
vector triple product
=
∑
i
=
1
n
m
i
[
Δ
r
i
×
−
(
Δ
r
i
×
α
)
+
ω
×
{
Δ
r
i
×
−
(
Δ
r
i
×
ω
)
}
]
…
cross-product anticommutativity
=
−
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
Δ
r
i
×
α
)
+
ω
×
{
Δ
r
i
×
(
Δ
r
i
×
ω
)
}
]
…
cross-product scalar multiplication
=
−
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
Δ
r
i
×
α
)
]
+
−
∑
i
=
1
n
m
i
[
ω
×
{
Δ
r
i
×
(
Δ
r
i
×
ω
)
}
]
…
summation distributivity
τ
=
−
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
Δ
r
i
×
α
)
]
+
ω
×
−
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
Δ
r
i
×
ω
)
]
…
ω
is not characteristic of particle
P
i
{\displaystyle {\begin{aligned}{\boldsymbol {\tau }}&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+{\boldsymbol {\omega }}\times \{{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})-\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot {\boldsymbol {\omega }})\}-{\boldsymbol {\omega }}\times {\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})]\;\ldots {\text{ cross-product distributivity over addition}}\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+{\boldsymbol {\omega }}\times \{{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})-\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot {\boldsymbol {\omega }})\}-(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})({\boldsymbol {\omega }}\times {\boldsymbol {\omega }})]\;\ldots {\text{ cross-product scalar multiplication}}\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+{\boldsymbol {\omega }}\times \{{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})-\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot {\boldsymbol {\omega }})\}-(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})(0)]\;\ldots {\text{ self cross-product}}\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+{\boldsymbol {\omega }}\times \{{\boldsymbol {\omega }}(\Delta \mathbf {r} _{i}\cdot \Delta \mathbf {r} _{i})-\Delta \mathbf {r} _{i}(\Delta \mathbf {r} _{i}\cdot {\boldsymbol {\omega }})\}]\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times ({\boldsymbol {\alpha }}\times \Delta \mathbf {r} _{i})+{\boldsymbol {\omega }}\times \{\Delta \mathbf {r} _{i}\times ({\boldsymbol {\omega }}\times \Delta \mathbf {r} _{i})\}]\;\ldots {\text{ vector triple product}}\\&=\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times -(\Delta \mathbf {r} _{i}\times {\boldsymbol {\alpha }})+{\boldsymbol {\omega }}\times \{\Delta \mathbf {r} _{i}\times -(\Delta \mathbf {r} _{i}\times {\boldsymbol {\omega }})\}]\;\ldots {\text{ cross-product anticommutativity}}\\&=-\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times {\boldsymbol {\alpha }})+{\boldsymbol {\omega }}\times \{\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times {\boldsymbol {\omega }})\}]\;\ldots {\text{ cross-product scalar multiplication}}\\&=-\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times {\boldsymbol {\alpha }})]+-\sum _{i=1}^{n}m_{i}[{\boldsymbol {\omega }}\times \{\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times {\boldsymbol {\omega }})\}]\;\ldots {\text{ summation distributivity}}\\{\boldsymbol {\tau }}&=-\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times {\boldsymbol {\alpha }})]+{\boldsymbol {\omega }}\times -\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times {\boldsymbol {\omega }})]\;\ldots \;{\boldsymbol {\omega }}{\text{ is not characteristic of particle }}P_{i}\end{aligned}}}
任意のベクトルに対して 、次の条件が成り立つことに注意してください。
u
{\displaystyle \mathbf {u} }
−
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
Δ
r
i
×
u
)
]
=
−
∑
i
=
1
n
m
i
(
[
0
−
Δ
r
3
,
i
Δ
r
2
,
i
Δ
r
3
,
i
0
−
Δ
r
1
,
i
−
Δ
r
2
,
i
Δ
r
1
,
i
0
]
(
[
0
−
Δ
r
3
,
i
Δ
r
2
,
i
Δ
r
3
,
i
0
−
Δ
r
1
,
i
−
Δ
r
2
,
i
Δ
r
1
,
i
0
]
[
u
1
u
2
u
3
]
)
)
…
cross-product as matrix multiplication
=
−
∑
i
=
1
n
m
i
(
[
0
−
Δ
r
3
,
i
Δ
r
2
,
i
Δ
r
3
,
i
0
−
Δ
r
1
,
i
−
Δ
r
2
,
i
Δ
r
1
,
i
0
]
[
−
Δ
r
3
,
i
u
2
+
Δ
r
2
,
i
u
3
+
Δ
r
3
,
i
u
1
−
Δ
r
1
,
i
u
3
−
Δ
r
2
,
i
u
1
+
Δ
r
1
,
i
u
2
]
)
=
−
∑
i
=
1
n
m
i
[
−
Δ
r
3
,
i
(
+
Δ
r
3
,
i
u
1
−
Δ
r
1
,
i
u
3
)
+
Δ
r
2
,
i
(
−
Δ
r
2
,
i
u
1
+
Δ
r
1
,
i
u
2
)
+
Δ
r
3
,
i
(
−
Δ
r
3
,
i
u
2
+
Δ
r
2
,
i
u
3
)
−
Δ
r
1
,
i
(
−
Δ
r
2
,
i
u
1
+
Δ
r
1
,
i
u
2
)
−
Δ
r
2
,
i
(
−
Δ
r
3
,
i
u
2
+
Δ
r
2
,
i
u
3
)
+
Δ
r
1
,
i
(
+
Δ
r
3
,
i
u
1
−
Δ
r
1
,
i
u
3
)
]
=
−
∑
i
=
1
n
m
i
[
−
Δ
r
3
,
i
2
u
1
+
Δ
r
1
,
i
Δ
r
3
,
i
u
3
−
Δ
r
2
,
i
2
u
1
+
Δ
r
1
,
i
Δ
r
2
,
i
u
2
−
Δ
r
3
,
i
2
u
2
+
Δ
r
2
,
i
Δ
r
3
,
i
u
3
+
Δ
r
2
,
i
Δ
r
1
,
i
u
1
−
Δ
r
1
,
i
2
u
2
+
Δ
r
3
,
i
Δ
r
2
,
i
u
2
−
Δ
r
2
,
i
2
u
3
+
Δ
r
3
,
i
Δ
r
1
,
i
u
1
−
Δ
r
1
,
i
2
u
3
]
=
−
∑
i
=
1
n
m
i
[
−
(
Δ
r
2
,
i
2
+
Δ
r
3
,
i
2
)
u
1
+
Δ
r
1
,
i
Δ
r
2
,
i
u
2
+
Δ
r
1
,
i
Δ
r
3
,
i
u
3
+
Δ
r
2
,
i
Δ
r
1
,
i
u
1
−
(
Δ
r
1
,
i
2
+
Δ
r
3
,
i
2
)
u
2
+
Δ
r
2
,
i
Δ
r
3
,
i
u
3
+
Δ
r
3
,
i
Δ
r
1
,
i
u
1
+
Δ
r
3
,
i
Δ
r
2
,
i
u
2
−
(
Δ
r
1
,
i
2
+
Δ
r
2
,
i
2
)
u
3
]
=
−
∑
i
=
1
n
m
i
[
−
(
Δ
r
2
,
i
2
+
Δ
r
3
,
i
2
)
Δ
r
1
,
i
Δ
r
2
,
i
Δ
r
1
,
i
Δ
r
3
,
i
Δ
r
2
,
i
Δ
r
1
,
i
−
(
Δ
r
1
,
i
2
+
Δ
r
3
,
i
2
)
Δ
r
2
,
i
Δ
r
3
,
i
Δ
r
3
,
i
Δ
r
1
,
i
Δ
r
3
,
i
Δ
r
2
,
i
−
(
Δ
r
1
,
i
2
+
Δ
r
2
,
i
2
)
]
[
u
1
u
2
u
3
]
=
−
∑
i
=
1
n
m
i
[
Δ
r
i
]
2
u
−
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
Δ
r
i
×
u
)
]
=
(
−
∑
i
=
1
n
m
i
[
Δ
r
i
]
2
)
u
…
u
is not characteristic of
P
i
{\displaystyle {\begin{aligned}-\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times \mathbf {u} )]&=-\sum _{i=1}^{n}m_{i}\left({\begin{bmatrix}0&-\Delta r_{3,i}&\Delta r_{2,i}\\\Delta r_{3,i}&0&-\Delta r_{1,i}\\-\Delta r_{2,i}&\Delta r_{1,i}&0\end{bmatrix}}\left({\begin{bmatrix}0&-\Delta r_{3,i}&\Delta r_{2,i}\\\Delta r_{3,i}&0&-\Delta r_{1,i}\\-\Delta r_{2,i}&\Delta r_{1,i}&0\end{bmatrix}}{\begin{bmatrix}u_{1}\\u_{2}\\u_{3}\end{bmatrix}}\right)\right)\;\ldots {\text{ cross-product as matrix multiplication}}\\[6pt]&=-\sum _{i=1}^{n}m_{i}\left({\begin{bmatrix}0&-\Delta r_{3,i}&\Delta r_{2,i}\\\Delta r_{3,i}&0&-\Delta r_{1,i}\\-\Delta r_{2,i}&\Delta r_{1,i}&0\end{bmatrix}}{\begin{bmatrix}-\Delta r_{3,i}\,u_{2}+\Delta r_{2,i}\,u_{3}\\+\Delta r_{3,i}\,u_{1}-\Delta r_{1,i}\,u_{3}\\-\Delta r_{2,i}\,u_{1}+\Delta r_{1,i}\,u_{2}\end{bmatrix}}\right)\\[6pt]&=-\sum _{i=1}^{n}m_{i}{\begin{bmatrix}-\Delta r_{3,i}(+\Delta r_{3,i}\,u_{1}-\Delta r_{1,i}\,u_{3})+\Delta r_{2,i}(-\Delta r_{2,i}\,u_{1}+\Delta r_{1,i}\,u_{2})\\+\Delta r_{3,i}(-\Delta r_{3,i}\,u_{2}+\Delta r_{2,i}\,u_{3})-\Delta r_{1,i}(-\Delta r_{2,i}\,u_{1}+\Delta r_{1,i}\,u_{2})\\-\Delta r_{2,i}(-\Delta r_{3,i}\,u_{2}+\Delta r_{2,i}\,u_{3})+\Delta r_{1,i}(+\Delta r_{3,i}\,u_{1}-\Delta r_{1,i}\,u_{3})\end{bmatrix}}\\[6pt]&=-\sum _{i=1}^{n}m_{i}{\begin{bmatrix}-\Delta r_{3,i}^{2}\,u_{1}+\Delta r_{1,i}\Delta r_{3,i}\,u_{3}-\Delta r_{2,i}^{2}\,u_{1}+\Delta r_{1,i}\Delta r_{2,i}\,u_{2}\\-\Delta r_{3,i}^{2}\,u_{2}+\Delta r_{2,i}\Delta r_{3,i}\,u_{3}+\Delta r_{2,i}\Delta r_{1,i}\,u_{1}-\Delta r_{1,i}^{2}\,u_{2}\\+\Delta r_{3,i}\Delta r_{2,i}\,u_{2}-\Delta r_{2,i}^{2}\,u_{3}+\Delta r_{3,i}\Delta r_{1,i}\,u_{1}-\Delta r_{1,i}^{2}\,u_{3}\end{bmatrix}}\\[6pt]&=-\sum _{i=1}^{n}m_{i}{\begin{bmatrix}-(\Delta r_{2,i}^{2}+\Delta r_{3,i}^{2})\,u_{1}+\Delta r_{1,i}\Delta r_{2,i}\,u_{2}+\Delta r_{1,i}\Delta r_{3,i}\,u_{3}\\+\Delta r_{2,i}\Delta r_{1,i}\,u_{1}-(\Delta r_{1,i}^{2}+\Delta r_{3,i}^{2})\,u_{2}+\Delta r_{2,i}\Delta r_{3,i}\,u_{3}\\+\Delta r_{3,i}\Delta r_{1,i}\,u_{1}+\Delta r_{3,i}\Delta r_{2,i}\,u_{2}-(\Delta r_{1,i}^{2}+\Delta r_{2,i}^{2})\,u_{3}\end{bmatrix}}\\[6pt]&=-\sum _{i=1}^{n}m_{i}{\begin{bmatrix}-(\Delta r_{2,i}^{2}+\Delta r_{3,i}^{2})&\Delta r_{1,i}\Delta r_{2,i}&\Delta r_{1,i}\Delta r_{3,i}\\\Delta r_{2,i}\Delta r_{1,i}&-(\Delta r_{1,i}^{2}+\Delta r_{3,i}^{2})&\Delta r_{2,i}\Delta r_{3,i}\\\Delta r_{3,i}\Delta r_{1,i}&\Delta r_{3,i}\Delta r_{2,i}&-(\Delta r_{1,i}^{2}+\Delta r_{2,i}^{2})\end{bmatrix}}{\begin{bmatrix}u_{1}\\u_{2}\\u_{3}\end{bmatrix}}\\&=-\sum _{i=1}^{n}m_{i}[\Delta r_{i}]^{2}\mathbf {u} \\[6pt]-\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times \mathbf {u} )]&=\left(-\sum _{i=1}^{n}m_{i}[\Delta r_{i}]^{2}\right)\mathbf {u} \;\ldots \;\mathbf {u} {\text{ is not characteristic of }}P_{i}\end{aligned}}}
最後に、結果を使用して次のように主な証明を完了します。
τ
=
−
∑
i
=
1
n
m
i
[
Δ
r
i
×
(
Δ
r
i
×
α
)
]
+
ω
×
−
∑
i
=
1
n
m
i
Δ
r
i
×
(
Δ
r
i
×
ω
)
]
=
(
−
∑
i
=
1
n
m
i
[
Δ
r
i
]
2
)
α
+
ω
×
(
−
∑
i
=
1
n
m
i
[
Δ
r
i
]
2
)
ω
{\displaystyle {\begin{aligned}{\boldsymbol {\tau }}&=-\sum _{i=1}^{n}m_{i}[\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times {\boldsymbol {\alpha }})]+{\boldsymbol {\omega }}\times -\sum _{i=1}^{n}m_{i}\Delta \mathbf {r} _{i}\times (\Delta \mathbf {r} _{i}\times {\boldsymbol {\omega }})]\\&=\left(-\sum _{i=1}^{n}m_{i}[\Delta r_{i}]^{2}\right){\boldsymbol {\alpha }}+{\boldsymbol {\omega }}\times \left(-\sum _{i=1}^{n}m_{i}[\Delta r_{i}]^{2}\right){\boldsymbol {\omega }}\end{aligned}}}
したがって、粒子の剛体システムにかかる合成トルクは、
質量中心に対する慣性行列
によって与えられます。
τ
=
I
C
α
+
ω
×
I
C
ω
,
{\displaystyle {\boldsymbol {\tau }}=\mathbf {I} _{\mathbf {C} }{\boldsymbol {\alpha }}+{\boldsymbol {\omega }}\times \mathbf {I} _{\mathbf {C} }{\boldsymbol {\omega }},}
I
C
{\displaystyle \mathbf {I_{C}} }
平行軸定理
物体の慣性行列は基準点の選択に依存する。質量中心に対する慣性行列 と他の点に対する慣性行列の間には有用な関係がある 。この関係は平行軸定理と呼ばれる。 [4] [7]
C
{\displaystyle \mathbf {C} }
R
{\displaystyle \mathbf {R} }
基準点を基準として測定された剛体粒子系に対して得られる 慣性行列を考える 。
I
R
{\displaystyle \mathbf {I_{R}} }
R
{\displaystyle \mathbf {R} }
I
R
=
−
∑
i
=
1
n
m
i
[
r
i
−
R
]
2
.
{\displaystyle \mathbf {I} _{\mathbf {R} }=-\sum _{i=1}^{n}m_{i}\left[\mathbf {r} _{i}-\mathbf {R} \right]^{2}.}
を剛体の質量の中心と
すると、 は 質量の中心から 基準点へのベクトル
です 。この式を使用して慣性行列を計算します。
C
{\displaystyle \mathbf {C} }
R
=
(
R
−
C
)
+
C
=
d
+
C
,
{\displaystyle \mathbf {R} =(\mathbf {R} -\mathbf {C} )+\mathbf {C} =\mathbf {d} +\mathbf {C} ,}
d
{\displaystyle \mathbf {d} }
C
{\displaystyle \mathbf {C} }
R
{\displaystyle \mathbf {R} }
I
R
=
−
∑
i
=
1
n
m
i
[
r
i
−
(
C
+
d
)
]
2
=
−
∑
i
=
1
n
m
i
[
(
r
i
−
C
)
−
d
]
2
.
{\displaystyle \mathbf {I} _{\mathbf {R} }=-\sum _{i=1}^{n}m_{i}[\mathbf {r} _{i}-\left(\mathbf {C} +\mathbf {d} \right)]^{2}=-\sum _{i=1}^{n}m_{i}[\left(\mathbf {r} _{i}-\mathbf {C} \right)-\mathbf {d} ]^{2}.}
外積に分配して
I
R
=
−
(
∑
i
=
1
n
m
i
[
r
i
−
C
]
2
)
+
(
∑
i
=
1
n
m
i
[
r
i
−
C
]
)
[
d
]
+
[
d
]
(
∑
i
=
1
n
m
i
[
r
i
−
C
]
)
−
(
∑
i
=
1
n
m
i
)
[
d
]
2
.
{\displaystyle \mathbf {I} _{\mathbf {R} }=-\left(\sum _{i=1}^{n}m_{i}[\mathbf {r} _{i}-\mathbf {C} ]^{2}\right)+\left(\sum _{i=1}^{n}m_{i}[\mathbf {r} _{i}-\mathbf {C} ]\right)[\mathbf {d} ]+[\mathbf {d} ]\left(\sum _{i=1}^{n}m_{i}[\mathbf {r} _{i}-\mathbf {C} ]\right)-\left(\sum _{i=1}^{n}m_{i}\right)[\mathbf {d} ]^{2}.}
最初の項は、質量中心に対する 慣性行列です。2 番目と 3 番目の項は、質量中心の定義によりゼロになります。最後の項は、システムの合計質量に、 から構築された 歪対称行列の 2 乗を乗じたものです 。
I
C
{\displaystyle \mathbf {I_{C}} }
C
{\displaystyle \mathbf {C} }
[
d
]
{\displaystyle [\mathbf {d} ]}
d
{\displaystyle \mathbf {d} }
結果は平行軸定理であり、
は 質量中心から 基準点へのベクトルです 。
I
R
=
I
C
−
M
[
d
]
2
,
{\displaystyle \mathbf {I} _{\mathbf {R} }=\mathbf {I} _{\mathbf {C} }-M[\mathbf {d} ]^{2},}
d
{\displaystyle \mathbf {d} }
C
{\displaystyle \mathbf {C} }
R
{\displaystyle \mathbf {R} }
マイナス符号に関する注意 : 基準点に対する位置ベクトルの歪対称行列を使用すると、各粒子の慣性行列は という形式になります。 これは、平面運動で現れる に似ています 。ただし、これを正しく機能させるには、マイナス符号が必要です。このマイナス符号は、 の歪対称性特性を使用して、必要に応じて項 に吸収できます 。
−
m
[
r
]
2
{\displaystyle -m\left[\mathbf {r} \right]^{2}}
m
r
2
{\displaystyle mr^{2}}
m
[
r
]
T
[
r
]
{\displaystyle m\left[\mathbf {r} \right]^{\mathsf {T}}\left[\mathbf {r} \right]}
[
r
]
{\displaystyle [\mathbf {r} ]}
平面におけるスカラー慣性モーメント
方向が単位ベクトルで指定され 、点で物体を通過する 指定軸の周りの物体のスカラー慣性モーメント は 、次の式で表されます。 [7]
ここで、 は基準点 に対するシステムの慣性モーメント行列 、 は ベクトル から得られる歪対称行列です 。
I
L
{\displaystyle I_{L}}
k
^
{\displaystyle \mathbf {\hat {k}} }
R
{\displaystyle \mathbf {R} }
I
L
=
k
^
⋅
(
−
∑
i
=
1
N
m
i
[
Δ
r
i
]
2
)
k
^
=
k
^
⋅
I
R
k
^
=
k
^
T
I
R
k
^
,
{\displaystyle I_{L}=\mathbf {\hat {k}} \cdot \left(-\sum _{i=1}^{N}m_{i}\left[\Delta \mathbf {r} _{i}\right]^{2}\right)\mathbf {\hat {k}} =\mathbf {\hat {k}} \cdot \mathbf {I} _{\mathbf {R} }\mathbf {\hat {k}} =\mathbf {\hat {k}} ^{\mathsf {T}}\mathbf {I} _{\mathbf {R} }\mathbf {\hat {k}} ,}
I
R
{\displaystyle \mathbf {I_{R}} }
R
{\displaystyle \mathbf {R} }
[
Δ
r
i
]
{\displaystyle [\Delta \mathbf {r} _{i}]}
Δ
r
i
=
r
i
−
R
{\displaystyle \Delta \mathbf {r} _{i}=\mathbf {r} _{i}-\mathbf {R} }
これは次のようにして導かれます。粒子 の剛体集合 が 座標 を持っているとします 。 を参照点 として を選択し、 参照点 を通る 単位ベクトルによって定義される直線 L の周りの慣性モーメントを計算します 。 この直線から粒子への垂直ベクトルは、 に投影される成分を除去することによって から得られます 。
ここで は 慣性行列との混同を避けるため単位行列であり、 は 直線 に沿った 単位ベクトルから形成される外積行列です 。
n
{\displaystyle n}
P
i
,
i
=
1
,
…
,
n
{\displaystyle P_{i},i=1,\dots ,n}
r
i
{\displaystyle \mathbf {r} _{i}}
R
{\displaystyle \mathbf {R} }
k
^
{\displaystyle \mathbf {\hat {k}} }
R
{\displaystyle \mathbf {R} }
L
(
t
)
=
R
+
t
k
^
{\displaystyle \mathbf {L} (t)=\mathbf {R} +t\mathbf {\hat {k}} }
P
i
{\displaystyle P_{i}}
Δ
r
i
{\displaystyle \Delta \mathbf {r} _{i}}
k
^
{\displaystyle \mathbf {\hat {k}} }
Δ
r
i
⊥
=
Δ
r
i
−
(
k
^
⋅
Δ
r
i
)
k
^
=
(
E
−
k
^
k
^
T
)
Δ
r
i
,
{\displaystyle \Delta \mathbf {r} _{i}^{\perp }=\Delta \mathbf {r} _{i}-\left(\mathbf {\hat {k}} \cdot \Delta \mathbf {r} _{i}\right)\mathbf {\hat {k}} =\left(\mathbf {E} -\mathbf {\hat {k}} \mathbf {\hat {k}} ^{\mathsf {T}}\right)\Delta \mathbf {r} _{i},}
E
{\displaystyle \mathbf {E} }
k
^
k
^
T
{\displaystyle \mathbf {\hat {k}} \mathbf {\hat {k}} ^{\mathsf {T}}}
k
^
{\displaystyle \mathbf {\hat {k}} }
L
{\displaystyle L}
このスカラー慣性モーメントを物体の慣性行列に関連付けるには、 となる歪対称行列を導入します。すると、 が単位ベクトルである
ことに注意して 恒等式が得られます
。
[
k
^
]
{\displaystyle \left[\mathbf {\hat {k}} \right]}
[
k
^
]
y
=
k
^
×
y
{\displaystyle \left[\mathbf {\hat {k}} \right]\mathbf {y} =\mathbf {\hat {k}} \times \mathbf {y} }
−
[
k
^
]
2
≡
|
k
^
|
2
(
E
−
k
^
k
^
T
)
=
E
−
k
^
k
^
T
,
{\displaystyle -\left[\mathbf {\hat {k}} \right]^{2}\equiv \left|\mathbf {\hat {k}} \right|^{2}\left(\mathbf {E} -\mathbf {\hat {k}} \mathbf {\hat {k}} ^{\mathsf {T}}\right)=\mathbf {E} -\mathbf {\hat {k}} \mathbf {\hat {k}} ^{\mathsf {T}},}
k
^
{\displaystyle \mathbf {\hat {k}} }
垂直ベクトルの大きさの2乗は
|
Δ
r
i
⊥
|
2
=
(
−
[
k
^
]
2
Δ
r
i
)
⋅
(
−
[
k
^
]
2
Δ
r
i
)
=
(
k
^
×
(
k
^
×
Δ
r
i
)
)
⋅
(
k
^
×
(
k
^
×
Δ
r
i
)
)
{\displaystyle {\begin{aligned}\left|\Delta \mathbf {r} _{i}^{\perp }\right|^{2}&=\left(-\left[\mathbf {\hat {k}} \right]^{2}\Delta \mathbf {r} _{i}\right)\cdot \left(-\left[\mathbf {\hat {k}} \right]^{2}\Delta \mathbf {r} _{i}\right)\\&=\left(\mathbf {\hat {k}} \times \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\right)\cdot \left(\mathbf {\hat {k}} \times \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\right)\end{aligned}}}
この式の簡略化には、点積と外積を入れ替えた 三重スカラー積の恒等式を使用します
。積を交換し、 と が 直交していることに注目して簡略化します。
(
k
^
×
(
k
^
×
Δ
r
i
)
)
⋅
(
k
^
×
(
k
^
×
Δ
r
i
)
)
≡
(
(
k
^
×
(
k
^
×
Δ
r
i
)
)
×
k
^
)
⋅
(
k
^
×
Δ
r
i
)
,
{\displaystyle \left(\mathbf {\hat {k}} \times \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\right)\cdot \left(\mathbf {\hat {k}} \times \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\right)\equiv \left(\left(\mathbf {\hat {k}} \times \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\right)\times \mathbf {\hat {k}} \right)\cdot \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right),}
Δ
r
i
{\displaystyle \Delta \mathbf {r} _{i}}
k
^
{\displaystyle \mathbf {\hat {k}} }
(
k
^
×
(
k
^
×
Δ
r
i
)
)
⋅
(
k
^
×
(
k
^
×
Δ
r
i
)
)
=
(
(
k
^
×
(
k
^
×
Δ
r
i
)
)
×
k
^
)
⋅
(
k
^
×
Δ
r
i
)
=
(
k
^
×
Δ
r
i
)
⋅
(
−
Δ
r
i
×
k
^
)
=
−
k
^
⋅
(
Δ
r
i
×
Δ
r
i
×
k
^
)
=
−
k
^
⋅
[
Δ
r
i
]
2
k
^
.
{\displaystyle {\begin{aligned}&\left(\mathbf {\hat {k}} \times \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\right)\cdot \left(\mathbf {\hat {k}} \times \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\right)\\={}&\left(\left(\mathbf {\hat {k}} \times \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\right)\times \mathbf {\hat {k}} \right)\cdot \left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\\={}&\left(\mathbf {\hat {k}} \times \Delta \mathbf {r} _{i}\right)\cdot \left(-\Delta \mathbf {r} _{i}\times \mathbf {\hat {k}} \right)\\={}&-\mathbf {\hat {k}} \cdot \left(\Delta \mathbf {r} _{i}\times \Delta \mathbf {r} _{i}\times \mathbf {\hat {k}} \right)\\={}&-\mathbf {\hat {k}} \cdot \left[\Delta \mathbf {r} _{i}\right]^{2}\mathbf {\hat {k}} .\end{aligned}}}
したがって、方向を 通る 直線の周りの慣性モーメントは 、計算から得られます
。
ここで、は 基準点に対するシステムの慣性モーメント行列です 。
L
{\displaystyle L}
R
{\displaystyle \mathbf {R} }
k
^
{\displaystyle \mathbf {\hat {k}} }
I
L
=
∑
i
=
1
N
m
i
|
Δ
r
i
⊥
|
2
=
−
∑
i
=
1
N
m
i
k
^
⋅
[
Δ
r
i
]
2
k
^
=
k
^
⋅
(
−
∑
i
=
1
N
m
i
[
Δ
r
i
]
2
)
k
^
=
k
^
⋅
I
R
k
^
=
k
^
T
I
R
k
^
,
{\displaystyle {\begin{aligned}I_{L}&=\sum _{i=1}^{N}m_{i}\left|\Delta \mathbf {r} _{i}^{\perp }\right|^{2}\\&=-\sum _{i=1}^{N}m_{i}\mathbf {\hat {k}} \cdot \left[\Delta \mathbf {r} _{i}\right]^{2}\mathbf {\hat {k}} =\mathbf {\hat {k}} \cdot \left(-\sum _{i=1}^{N}m_{i}\left[\Delta \mathbf {r} _{i}\right]^{2}\right)\mathbf {\hat {k}} \\&=\mathbf {\hat {k}} \cdot \mathbf {I} _{\mathbf {R} }\mathbf {\hat {k}} =\mathbf {\hat {k}} ^{\mathsf {T}}\mathbf {I} _{\mathbf {R} }\mathbf {\hat {k}} ,\end{aligned}}}
I
R
{\displaystyle \mathbf {I_{R}} }
R
{\displaystyle \mathbf {R} }
これは、慣性行列を使用して、物体内の任意の指定された回転軸の周りの物体の慣性モーメントを計算できることを示しています。
慣性テンソル
同じ物体でも、回転軸が異なれば、その軸の周りの慣性モーメントも異なります。一般に、物体がすべての軸に対して対称でない限り、慣性モーメントは等しくなりません。 慣性モーメント テンソルは 、物体のすべての慣性モーメントを 1 つの量でまとめる便利な方法です。慣性モーメントは空間内の任意の点を基準に計算できますが、実用上は質量中心が最も一般的に使用されます。
意味
質点 からなる剛体の場合 、慣性モーメント テンソルは 次のように表される。
N
{\displaystyle N}
m
k
{\displaystyle m_{k}}
I
=
[
I
11
I
12
I
13
I
21
I
22
I
23
I
31
I
32
I
33
]
.
{\displaystyle \mathbf {I} ={\begin{bmatrix}I_{11}&I_{12}&I_{13}\\I_{21}&I_{22}&I_{23}\\I_{31}&I_{32}&I_{33}\end{bmatrix}}.}
その構成要素は次のように定義される。
I
i
j
=
d
e
f
∑
k
=
1
N
m
k
(
‖
r
k
‖
2
δ
i
j
−
x
i
(
k
)
x
j
(
k
)
)
{\displaystyle I_{ij}\ {\stackrel {\mathrm {def} }{=}}\ \sum _{k=1}^{N}m_{k}\left(\left\|\mathbf {r} _{k}\right\|^{2}\delta _{ij}-x_{i}^{(k)}x_{j}^{(k)}\right)}
どこ
i
{\displaystyle i}
、、およびの場合 、それぞれ 1、2 、または 3 に等しい 。
j
{\displaystyle j}
x
{\displaystyle x}
y
{\displaystyle y}
z
{\displaystyle z}
r
k
=
(
x
1
(
k
)
,
x
2
(
k
)
,
x
3
(
k
)
)
{\displaystyle \mathbf {r} _{k}=\left(x_{1}^{(k)},x_{2}^{(k)},x_{3}^{(k)}\right)}
テンソルが計算される点から 質点へのベクトルであり、
m
k
{\displaystyle m_{k}}
δ
i
j
{\displaystyle \delta _{ij}}
クロネッカーのデルタ です 。
定義により、は 対称テンソル であることに注意してください 。
I
{\displaystyle \mathbf {I} }
対角要素はより簡潔に次のように記述される。
I
x
x
=
d
e
f
∑
k
=
1
N
m
k
(
y
k
2
+
z
k
2
)
,
I
y
y
=
d
e
f
∑
k
=
1
N
m
k
(
x
k
2
+
z
k
2
)
,
I
z
z
=
d
e
f
∑
k
=
1
N
m
k
(
x
k
2
+
y
k
2
)
,
{\displaystyle {\begin{aligned}I_{xx}\ &{\stackrel {\mathrm {def} }{=}}\ \sum _{k=1}^{N}m_{k}\left(y_{k}^{2}+z_{k}^{2}\right),\\I_{yy}\ &{\stackrel {\mathrm {def} }{=}}\ \sum _{k=1}^{N}m_{k}\left(x_{k}^{2}+z_{k}^{2}\right),\\I_{zz}\ &{\stackrel {\mathrm {def} }{=}}\ \sum _{k=1}^{N}m_{k}\left(x_{k}^{2}+y_{k}^{2}\right),\end{aligned}}}
一方、非対角要素は 慣性積 とも呼ばれ、
I
x
y
=
I
y
x
=
d
e
f
−
∑
k
=
1
N
m
k
x
k
y
k
,
I
x
z
=
I
z
x
=
d
e
f
−
∑
k
=
1
N
m
k
x
k
z
k
,
I
y
z
=
I
z
y
=
d
e
f
−
∑
k
=
1
N
m
k
y
k
z
k
.
{\displaystyle {\begin{aligned}I_{xy}=I_{yx}\ &{\stackrel {\mathrm {def} }{=}}\ -\sum _{k=1}^{N}m_{k}x_{k}y_{k},\\I_{xz}=I_{zx}\ &{\stackrel {\mathrm {def} }{=}}\ -\sum _{k=1}^{N}m_{k}x_{k}z_{k},\\I_{yz}=I_{zy}\ &{\stackrel {\mathrm {def} }{=}}\ -\sum _{k=1}^{N}m_{k}y_{k}z_{k}.\end{aligned}}}
ここで、は 物体が x 軸を中心に回転するときの - 軸周り の慣性モーメントを表し、は 物体が - 軸を中心に回転するときの - 軸 周りの慣性モーメントを表します 。
I
x
x
{\displaystyle I_{xx}}
x
{\displaystyle x}
I
x
y
{\displaystyle I_{xy}}
y
{\displaystyle y}
x
{\displaystyle x}
これらの量は、スカラー慣性モーメントと同様に、質量密度関数で記述される分散質量を持つ物体に一般化することができる。
I
=
∭
V
ρ
(
x
,
y
,
z
)
(
‖
r
‖
2
E
3
−
r
⊗
r
)
d
x
d
y
d
z
,
{\displaystyle \mathbf {I} =\iiint _{V}\rho (x,y,z)\left(\|\mathbf {r} \|^{2}\mathbf {E} _{3}-\mathbf {r} \otimes \mathbf {r} \right)\,dx\,dy\,dz,}
ここで 、 はそれらの 外積 、 E 3 は 3×3 単位行列 、 V は オブジェクトを完全に含む空間領域です。
r
⊗
r
{\displaystyle \mathbf {r} \otimes \mathbf {r} }
あるいは、角運動量演算子を 使って次のように書くこともできます 。
[
r
]
x
=
r
×
x
{\displaystyle [\mathbf {r} ]\mathbf {x} =\mathbf {r} \times \mathbf {x} }
I
=
∭
V
ρ
(
r
)
[
r
]
T
[
r
]
d
V
=
−
∭
Q
ρ
(
r
)
[
r
]
2
d
V
{\displaystyle \mathbf {I} =\iiint _{V}\rho (\mathbf {r} )[\mathbf {r} ]^{\textsf {T}}[\mathbf {r} ]\,dV=-\iiint _{Q}\rho (\mathbf {r} )[\mathbf {r} ]^{2}\,dV}
慣性テンソルは慣性行列と同じように、任意の軸の周りのスカラー慣性モーメントを計算するために使用できます 。
n
{\displaystyle \mathbf {n} }
I
n
=
n
⋅
I
⋅
n
,
{\displaystyle I_{n}=\mathbf {n} \cdot \mathbf {I} \cdot \mathbf {n} ,}
ここで、 ドット積は 、成分テンソルの対応する要素で取られます。 計算によって、慣性項の積、たとえば が得られ
、これは、 物体が - 軸の周りを回転するときの - 軸
の周りの慣性モーメントとして解釈できます 。
I
12
{\displaystyle I_{12}}
I
12
=
e
1
⋅
I
⋅
e
2
,
{\displaystyle I_{12}=\mathbf {e} _{1}\cdot \mathbf {I} \cdot \mathbf {e} _{2},}
x
{\displaystyle x}
y
{\displaystyle y}
2次のテンソルの成分は行列にまとめることができる。慣性テンソルの場合、この行列は次のように表される。
I
=
[
I
11
I
12
I
13
I
21
I
22
I
23
I
31
I
32
I
33
]
=
[
I
x
x
I
x
y
I
x
z
I
y
x
I
y
y
I
y
z
I
z
x
I
z
y
I
z
z
]
=
[
∑
k
=
1
N
m
k
(
y
k
2
+
z
k
2
)
−
∑
k
=
1
N
m
k
x
k
y
k
−
∑
k
=
1
N
m
k
x
k
z
k
−
∑
k
=
1
N
m
k
x
k
y
k
∑
k
=
1
N
m
k
(
x
k
2
+
z
k
2
)
−
∑
k
=
1
N
m
k
y
k
z
k
−
∑
k
=
1
N
m
k
x
k
z
k
−
∑
k
=
1
N
m
k
y
k
z
k
∑
k
=
1
N
m
k
(
x
k
2
+
y
k
2
)
]
.
{\displaystyle \mathbf {I} ={\begin{bmatrix}I_{11}&I_{12}&I_{13}\\[1.8ex]I_{21}&I_{22}&I_{23}\\[1.8ex]I_{31}&I_{32}&I_{33}\end{bmatrix}}={\begin{bmatrix}I_{xx}&I_{xy}&I_{xz}\\[1.8ex]I_{yx}&I_{yy}&I_{yz}\\[1.8ex]I_{zx}&I_{zy}&I_{zz}\end{bmatrix}}={\begin{bmatrix}\sum _{k=1}^{N}m_{k}\left(y_{k}^{2}+z_{k}^{2}\right)&-\sum _{k=1}^{N}m_{k}x_{k}y_{k}&-\sum _{k=1}^{N}m_{k}x_{k}z_{k}\\[1ex]-\sum _{k=1}^{N}m_{k}x_{k}y_{k}&\sum _{k=1}^{N}m_{k}\left(x_{k}^{2}+z_{k}^{2}\right)&-\sum _{k=1}^{N}m_{k}y_{k}z_{k}\\[1ex]-\sum _{k=1}^{N}m_{k}x_{k}z_{k}&-\sum _{k=1}^{N}m_{k}y_{k}z_{k}&\sum _{k=1}^{N}m_{k}\left(x_{k}^{2}+y_{k}^{2}\right)\end{bmatrix}}.}
剛体力学では、 慣性テンソルの成分として、
やなどの 、 軸、 軸を明示的に識別する表記法を使用するのが一般的です。
x
{\displaystyle x}
y
{\displaystyle y}
z
{\displaystyle z}
I
x
x
{\displaystyle I_{xx}}
I
x
y
{\displaystyle I_{xy}}
代替慣性慣性
SolidWorks、Unigraphics NX/Siemens NX、MSC Adams などの一部の CAD および CAE アプリケーションでは、慣性積に別の規則が使用されています。この規則によると、慣性積の式からマイナス記号が削除され、代わりに慣性マトリックスに挿入されます。
I
x
y
=
I
y
x
=
d
e
f
∑
k
=
1
N
m
k
x
k
y
k
,
I
x
z
=
I
z
x
=
d
e
f
∑
k
=
1
N
m
k
x
k
z
k
,
I
y
z
=
I
z
y
=
d
e
f
∑
k
=
1
N
m
k
y
k
z
k
,
I
=
[
I
11
I
12
I
13
I
21
I
22
I
23
I
31
I
32
I
33
]
=
[
I
x
x
−
I
x
y
−
I
x
z
−
I
y
x
I
y
y
−
I
y
z
−
I
z
x
−
I
z
y
I
z
z
]
=
[
∑
k
=
1
N
m
k
(
y
k
2
+
z
k
2
)
−
∑
k
=
1
N
m
k
x
k
y
k
−
∑
k
=
1
N
m
k
x
k
z
k
−
∑
k
=
1
N
m
k
x
k
y
k
∑
k
=
1
N
m
k
(
x
k
2
+
z
k
2
)
−
∑
k
=
1
N
m
k
y
k
z
k
−
∑
k
=
1
N
m
k
x
k
z
k
−
∑
k
=
1
N
m
k
y
k
z
k
∑
k
=
1
N
m
k
(
x
k
2
+
y
k
2
)
]
.
{\displaystyle {\begin{aligned}I_{xy}=I_{yx}\ &{\stackrel {\mathrm {def} }{=}}\ \sum _{k=1}^{N}m_{k}x_{k}y_{k},\\I_{xz}=I_{zx}\ &{\stackrel {\mathrm {def} }{=}}\ \sum _{k=1}^{N}m_{k}x_{k}z_{k},\\I_{yz}=I_{zy}\ &{\stackrel {\mathrm {def} }{=}}\ \sum _{k=1}^{N}m_{k}y_{k}z_{k},\\[3pt]\mathbf {I} ={\begin{bmatrix}I_{11}&I_{12}&I_{13}\\[1.8ex]I_{21}&I_{22}&I_{23}\\[1.8ex]I_{31}&I_{32}&I_{33}\end{bmatrix}}&={\begin{bmatrix}I_{xx}&-I_{xy}&-I_{xz}\\[1.8ex]-I_{yx}&I_{yy}&-I_{yz}\\[1.8ex]-I_{zx}&-I_{zy}&I_{zz}\end{bmatrix}}={\begin{bmatrix}\sum _{k=1}^{N}m_{k}\left(y_{k}^{2}+z_{k}^{2}\right)&-\sum _{k=1}^{N}m_{k}x_{k}y_{k}&-\sum _{k=1}^{N}m_{k}x_{k}z_{k}\\[1ex]-\sum _{k=1}^{N}m_{k}x_{k}y_{k}&\sum _{k=1}^{N}m_{k}\left(x_{k}^{2}+z_{k}^{2}\right)&-\sum _{k=1}^{N}m_{k}y_{k}z_{k}\\[1ex]-\sum _{k=1}^{N}m_{k}x_{k}z_{k}&-\sum _{k=1}^{N}m_{k}y_{k}z_{k}&\sum _{k=1}^{N}m_{k}\left(x_{k}^{2}+y_{k}^{2}\right)\end{bmatrix}}.\end{aligned}}}
慣性慣性を決定する(主軸法)
慣性データがあっても、 どの慣性規則が使用されているかがわからなければ、主軸も持っているかどうかを判断できます。主軸法では、次の 2 つの仮定から慣性行列を作成します。
(
I
x
x
,
I
y
y
,
I
z
z
,
I
x
y
,
I
x
z
,
I
y
z
)
{\displaystyle (I_{xx},I_{yy},I_{zz},I_{xy},I_{xz},I_{yz})}
標準慣性規則が使用されています 。
(
I
12
=
I
x
y
,
I
13
=
I
x
z
,
I
23
=
I
y
z
)
{\displaystyle (I_{12}=I_{xy},I_{13}=I_{xz},I_{23}=I_{yz})}
代替慣性規則が使用されています 。
(
I
12
=
−
I
x
y
,
I
13
=
−
I
x
z
,
I
23
=
−
I
y
z
)
{\displaystyle (I_{12}=-I_{xy},I_{13}=-I_{xz},I_{23}=-I_{yz})}
次に、2 つの行列の固有ベクトルを計算します。固有ベクトルが主軸に平行な行列は、使用されている慣性規則に対応します。
テンソル成分の導出
原点を通る回転軸から方向への 粒子の 距離 は であり 、 は 単位ベクトルである。軸上の慣性モーメントは
r
{\displaystyle r}
x
{\displaystyle \mathbf {x} }
n
^
{\displaystyle \mathbf {\hat {n}} }
|
x
−
(
x
⋅
n
^
)
n
^
|
{\displaystyle \left|\mathbf {x} -\left(\mathbf {x} \cdot \mathbf {\hat {n}} \right)\mathbf {\hat {n}} \right|}
n
^
{\displaystyle \mathbf {\hat {n}} }
I
=
m
r
2
=
m
(
x
−
(
x
⋅
n
^
)
n
^
)
⋅
(
x
−
(
x
⋅
n
^
)
n
^
)
=
m
(
x
2
−
2
x
(
x
⋅
n
^
)
n
^
+
(
x
⋅
n
^
)
2
n
^
2
)
=
m
(
x
2
−
(
x
⋅
n
^
)
2
)
.
{\displaystyle I=mr^{2}=m\left(\mathbf {x} -\left(\mathbf {x} \cdot \mathbf {\hat {n}} \right)\mathbf {\hat {n}} \right)\cdot \left(\mathbf {x} -\left(\mathbf {x} \cdot \mathbf {\hat {n}} \right)\mathbf {\hat {n}} \right)=m\left(\mathbf {x} ^{2}-2\mathbf {x} \left(\mathbf {x} \cdot \mathbf {\hat {n}} \right)\mathbf {\hat {n}} +\left(\mathbf {x} \cdot \mathbf {\hat {n}} \right)^{2}\mathbf {\hat {n}} ^{2}\right)=m\left(\mathbf {x} ^{2}-\left(\mathbf {x} \cdot \mathbf {\hat {n}} \right)^{2}\right).}
行列転置を 使用して方程式を書き直します 。
I
=
m
(
x
T
x
−
n
^
T
x
x
T
n
^
)
=
m
⋅
n
^
T
(
x
T
x
⋅
E
3
−
x
x
T
)
n
^
,
{\displaystyle I=m\left(\mathbf {x} ^{\textsf {T}}\mathbf {x} -\mathbf {\hat {n}} ^{\textsf {T}}\mathbf {x} \mathbf {x} ^{\textsf {T}}\mathbf {\hat {n}} \right)=m\cdot \mathbf {\hat {n}} ^{\textsf {T}}\left(\mathbf {x} ^{\textsf {T}}\mathbf {x} \cdot \mathbf {E_{3}} -\mathbf {x} \mathbf {x} ^{\textsf {T}}\right)\mathbf {\hat {n}} ,}
ここで、 E 3 は 3×3の 単位行列 です。
これは慣性モーメントのテンソル式につながる。
I
=
m
[
n
1
n
2
n
3
]
[
y
2
+
z
2
−
x
y
−
x
z
−
y
x
x
2
+
z
2
−
y
z
−
z
x
−
z
y
x
2
+
y
2
]
[
n
1
n
2
n
3
]
.
{\displaystyle I=m{\begin{bmatrix}n_{1}&n_{2}&n_{3}\end{bmatrix}}{\begin{bmatrix}y^{2}+z^{2}&-xy&-xz\\[0.5ex]-yx&x^{2}+z^{2}&-yz\\[0.5ex]-zx&-zy&x^{2}+y^{2}\end{bmatrix}}{\begin{bmatrix}n_{1}\\[0.7ex]n_{2}\\[0.7ex]n_{3}\end{bmatrix}}.}
複数の粒子の場合、この式が正しいことを確認するには、慣性モーメントが加算されるということを思い出すだけで十分です。
並進の慣性テンソル
を物体 の質量中心 で計算した慣性テンソルとし 、を 物体の変位ベクトルとします。移動した物体の元の質量中心に対する慣性テンソルは次のように求められます。
ここで 、 は物体の質量、 E 3 は 3 × 3 単位行列、 は 外積 です 。
I
0
{\displaystyle \mathbf {I} _{0}}
R
{\displaystyle \mathbf {R} }
I
=
I
0
+
m
[
(
R
⋅
R
)
E
3
−
R
⊗
R
]
{\displaystyle \mathbf {I} =\mathbf {I} _{0}+m[(\mathbf {R} \cdot \mathbf {R} )\mathbf {E} _{3}-\mathbf {R} \otimes \mathbf {R} ]}
m
{\displaystyle m}
⊗
{\displaystyle \otimes }
回転の慣性テンソル
物体の回転を表す 行列を と する 。回転する物体の慣性テンソルは次のように与えられる: [28]
R
{\displaystyle \mathbf {R} }
I
=
R
I
0
R
T
{\displaystyle \mathbf {I} =\mathbf {R} \mathbf {I_{0}} \mathbf {R} ^{\textsf {T}}}
異なる参照フレームにおける慣性行列
ニュートンの第2法則における慣性行列の使用は、その成分が慣性フレームに平行な軸を基準として計算され、物体に固定された参照フレームを基準として計算されないことを前提としています。 [7] [25] これは、物体が移動すると、慣性行列の成分が時間とともに変化することを意味します。対照的に、物体に固定されたフレームで測定された慣性行列の成分は一定です。
ボディフレーム
質量中心に対するボディフレームの慣性行列を と表し 、慣性フレームに対するボディフレームの向きを回転行列 で定義します 。ここで、
ボディ固定座標フレームの
ベクトルは 慣性フレームの座標を持ちます。すると、慣性フレームで測定されたボディの慣性行列は次のように表されます。
I
C
B
{\displaystyle \mathbf {I} _{\mathbf {C} }^{B}}
A
{\displaystyle \mathbf {A} }
x
=
A
y
,
{\displaystyle \mathbf {x} =\mathbf {A} \mathbf {y} ,}
y
{\displaystyle \mathbf {y} }
x
{\displaystyle \mathbf {x} }
I
C
=
A
I
C
B
A
T
.
{\displaystyle \mathbf {I} _{\mathbf {C} }=\mathbf {A} \mathbf {I} _{\mathbf {C} }^{B}\mathbf {A} ^{\mathsf {T}}.}
体が動くと が変化するのに対し、 は一定のままであること に注目してください 。
A
{\displaystyle \mathbf {A} }
I
C
B
{\displaystyle \mathbf {I} _{\mathbf {C} }^{B}}
主軸
物体座標系で測定すると、慣性行列は定数の実対称行列となる。実対称行列は 回転行列 と対角行列の積に 固有値分解 され、
次
のように表される。
Q
{\displaystyle \mathbf {Q} }
Λ
{\displaystyle {\boldsymbol {\Lambda }}}
I
C
B
=
Q
Λ
Q
T
,
{\displaystyle \mathbf {I} _{\mathbf {C} }^{B}=\mathbf {Q} {\boldsymbol {\Lambda }}\mathbf {Q} ^{\mathsf {T}},}
Λ
=
[
I
1
0
0
0
I
2
0
0
0
I
3
]
.
{\displaystyle {\boldsymbol {\Lambda }}={\begin{bmatrix}I_{1}&0&0\\0&I_{2}&0\\0&0&I_{3}\end{bmatrix}}.}
回転行列の列は 物体の主軸の方向を定義し、定数 、、 は 主慣性モーメント と呼ばれます。この結果は 、JJシルベスター(1852) によって最初に示され、 シルベスターの慣性法則 の一種です 。 [29] [30] 慣性モーメントが最も大きい主軸は、 図形軸 または 図形の軸 と呼ばれることもあります。
Q
{\displaystyle \mathbf {Q} }
I
1
{\displaystyle I_{1}}
I
2
{\displaystyle I_{2}}
I
3
{\displaystyle I_{3}}
おもちゃの コマは 回転する剛体の一例で、 コマ という言葉は剛体の種類の名前に使われます。すべての主慣性モーメントが異なる場合、 質量の中心 を通る主軸は一意に指定され、剛体は 非対称コマ と呼ばれます。2 つの主モーメントが同じ場合、剛体は 対称コマ と呼ばれ、対応する 2 つの主軸は一意に選択できません。3 つの主モーメントがすべて同じ場合、剛体は 球状コマ と呼ばれ(球状である必要はありません)、どの軸も主軸と見なすことができます。つまり、どの軸についても慣性モーメントは同じです。
主軸は、多くの場合、オブジェクトの対称軸と一致します。剛体が の順序の対称軸を持つ場合、つまり、 指定された軸を中心に 360° / m の回転に対して対称である場合、その軸は主軸です。 のとき 、剛体は対称トップです。剛体が互いに平行でも垂直でもない少なくとも 2 つの対称軸を持つ場合、それは球状トップ、たとえば立方体やその他の プラトン立体 です。
m
{\displaystyle m}
m
>
2
{\displaystyle m>2}
車両 の 動き は、 ヨー、ピッチ、ロール の観点から説明されることが多く、 これらは通常、3 つの主軸の周りの回転にほぼ相当します。車両が左右対称である場合、主軸の 1 つは、横方向 (ピッチ) 軸と正確に対応します。
この数学的現象の実際の例としては、自動車のタイヤ のバランスをとるという日常的な作業が挙げられます 。これは基本的に、車輪の質量の分布を調整して、車輪の慣性主軸が車軸と一直線になるようにし、車輪がぐらつかないようにすることを意味します。
回転分子は 非対称、対称、球状トップに分類され、 回転スペクトル の構造はタイプごとに異なります。
楕円
半主直径が 、 、 と表示された楕円体 。
a
{\displaystyle a}
b
{\displaystyle b}
c
{\displaystyle c}
物体座標系における慣性モーメント行列は、 ポアンソの楕円体 と呼ばれる物体の表面を定義する2次形式です。 [31] 主軸に沿った質量中心に対する慣性行列をとれば、表面
または
は
物体座標系における 楕円体
を定義します 。この方程式を の形式で書き表すと、
この楕円体の半主直径は次のように表されることがわかります。
Λ
{\displaystyle {\boldsymbol {\Lambda }}}
x
T
Λ
x
=
1
,
{\displaystyle \mathbf {x} ^{\mathsf {T}}{\boldsymbol {\Lambda }}\mathbf {x} =1,}
I
1
x
2
+
I
2
y
2
+
I
3
z
2
=
1
,
{\displaystyle I_{1}x^{2}+I_{2}y^{2}+I_{3}z^{2}=1,}
(
x
1
/
I
1
)
2
+
(
y
1
/
I
2
)
2
+
(
z
1
/
I
3
)
2
=
1
,
{\displaystyle \left({\frac {x}{1/{\sqrt {I_{1}}}}}\right)^{2}+\left({\frac {y}{1/{\sqrt {I_{2}}}}}\right)^{2}+\left({\frac {z}{1/{\sqrt {I_{3}}}}}\right)^{2}=1,}
a
=
1
I
1
,
b
=
1
I
2
,
c
=
1
I
3
.
{\displaystyle a={\frac {1}{\sqrt {I_{1}}}},\quad b={\frac {1}{\sqrt {I_{2}}}},\quad c={\frac {1}{\sqrt {I_{3}}}}.}
この楕円体上の点を その大きさと方向 で定義するとします。 ここで は単位ベクトルです。すると、慣性行列と 方向の軸の周りの スカラー慣性モーメントの間には、上に示した関係から次の式 が得られます
。
x
{\displaystyle \mathbf {x} }
x
=
‖
x
‖
n
{\displaystyle \mathbf {x} =\|\mathbf {x} \|\mathbf {n} }
n
{\displaystyle \mathbf {n} }
I
n
{\displaystyle I_{\mathbf {n} }}
n
{\displaystyle \mathbf {n} }
x
T
Λ
x
=
‖
x
‖
2
n
T
Λ
n
=
‖
x
‖
2
I
n
=
1.
{\displaystyle \mathbf {x} ^{\mathsf {T}}{\boldsymbol {\Lambda }}\mathbf {x} =\|\mathbf {x} \|^{2}\mathbf {n} ^{\mathsf {T}}{\boldsymbol {\Lambda }}\mathbf {n} =\|\mathbf {x} \|^{2}I_{\mathbf {n} }=1.}
したがって、慣性楕円体上の
方向の 点の大きさは
x
{\displaystyle \mathbf {x} }
n
{\displaystyle \mathbf {n} }
‖
x
‖
=
1
I
n
.
{\displaystyle \|\mathbf {x} \|={\frac {1}{\sqrt {I_{\mathbf {n} }}}}.}
参照
参考文献
^ エスクディエ、マルセル、アトキンス、トニー(2019)。機械工学辞典(第2版)。オックスフォード大学出版局。doi : 10.1093 / acref/9780198832102.001.0001。ISBN 978-0-19-883210-2 。
^ ab Mach, Ernst (1919). The Science of Mechanics. pp. 173–187 . 2014年 11月21日 閲覧 。
^ オイラー、レオンハルト (1765)。 Theoria motus corporum Solidorum seu pigidorum: Ex primis nostrae cognitionis principiis stackilita et ad omnes motus, qui in huiusmodi corpora cadere possunt, accommodata [固体または剛体の運動理論: 私たちの知識の第一原理に基づいて確立され、次のようなすべての運動に適しています。そのような体では発生する可能性があります。](ラテン語)。ロストックとグライフスヴァルト (ドイツ): AF Röse。 p. 166.ISBN 978-1-4297-4281-8 。 166 ページより: 「定義 7. 422. 慣性体運動量は、eujuspiam 軸を中心として、合計オムニウム製品、クェ・オリウントゥル、四個体距離当たりの身体要素を単数形にする。」 (定義 7. 422。任意の軸に関する物体の慣性モーメントは、物体の個々の要素に軸からの距離の 2 乗を乗算した場合に生じる積のすべての合計です。)
^ abcdef マリオン、JB; ソーントン、ST (1995)。 粒子とシステムの古典的ダイナミクス (第4版)。トムソン 。ISBN 0-03-097302-3 。
^ ab Symon, KR (1971). 力学 (第3版). Addison-Wesley. ISBN 0-201-07392-7 。
^ ab Tenenbaum, RA (2004). 応用ダイナミクスの基礎 . Springer. ISBN 0-387-00887-X 。
^ abcdefgh Kane, TR; Levinson, DA (1985). ダイナミクス、理論およびアプリケーション 。ニューヨーク:McGraw-Hill。
^ ab Winn, Will (2010). 理解できる物理学入門: 第1巻 - 力学. AuthorHouse. p. 10.10. ISBN 978-1449063337 。
^ abフラートン、ダン (2011)。優等生向け物理学要点。シリービーグルプロダクション 。pp . 142–143。ISBN 978-0983563334 。
^ Wolfram, Stephen (2014). 「Spinning Ice Skater」. Wolfram Demonstrations Project . Mathematica, Inc. 2014年 9月30日 閲覧 。
^ Hokin, Samuel (2014). 「フィギュアスケートのスピン」。 日常の物理学。 2014年 9月30日 閲覧 。
^ ブライトハウプト、ジム (2000)。『物理学の理解:上級者向け』ネルソン・トーマス、p. 64。ISBN 0748743146 。
^ クロウェル、ベンジャミン (2003)。 保存則 。光と物質。pp. 107。ISBN 0970467028 アイス スケーターの角運動量保存則。
^ Tipler, Paul A. (1999). 科学者とエンジニアのための物理学、第1巻:力学、振動と波、熱力学。マクミラン。p. 304。ISBN 1572594918 。
^ abcde ポール、バートン(1979年6月)。 平面機械の運動学と動力学 。プレンティスホール 。ISBN
978-0135160626 。
^ ハリデイ、デイビッド、レスニック、ロバート、ウォーカー、ジェール (2005)。 物理学の基礎 (第 7 版)。ホーボーケン、ニュージャージー: ワイリー 。ISBN
9780471216438 。
^ French, AP (1971). 振動と波 . ボカラトン、フロリダ州: CRC プレス. ISBN
9780748744473 。
^ abcdef ウイッカー、ジョン・J.; ペノック、ゴードン・R.; シグリー、ジョセフ・E. (2010). 機械とメカニズムの理論 (第4版). オックスフォード大学出版局. ISBN
978-0195371239 。
^ C. Couch および J. Mayes、「MOI のための三線振り子」、Happresearch.com、2016 年。
^ グレーシー、ウィリアム、簡略化された複合振り子法による飛行機の慣性モーメントの実験的測定、NACA 技術ノート No. 1629、1948
^ Morrow, HW; Kokernak, Robert (2011). 静力学と材料強度 (第7版)。ニュージャージー:Prentice Hall。pp. 192–196。ISBN 978-0135034521 。
^ この場合、この慣性モーメントは、その軸に沿って加えられたトルクがその軸を中心に回転を引き起こすことのみを説明します。ただし、主軸に沿っていないトルクは、他の軸を中心に回転を引き起こします。
^ abcdefghi Ferdinand P. Beer、E. Russell Johnston, Jr.、Phillip J. Cornwell (2010)。 エンジニアのためのベクトル力学:ダイナミクス (第9版)。ボストン:McGraw- Hill。ISBN 978-0077295493 。
^ Walter D. Pilkey、「弾性梁の解析と設計:計算方法」、John Wiley、2002年。
^ ab ゴールドスタイン、H. (1980)。 古典力学 (第2版)。アディソン・ウェズレー 。ISBN 0-201-02918-9 。
^ LD Landau および EM Lifshitz、「力学」第 1 巻、第 2 版、Pergamon Press、1969 年。
^ LW Tsai、「ロボット分析:シリアルおよびパラレルマニピュレータのメカニズム」、John-Wiley、NY、1999年。
^ David, Baraff. 「物理ベースモデリング - 剛体シミュレーション」 (PDF) 。Pixar Graphics Technologies 。
^ Sylvester, JJ (1852). 「すべての同次二次多項式は、実直交代入によって正負の平方の和の形に簡約可能であるという定理の証明」 (PDF) 。 Philosophical Magazine 。第 4 シリーズ。4 ( 23): 138–142。 doi :10.1080/14786445208647087 。 2008 年 6 月 27 日 閲覧 。
^ ノーマン、CW (1986)。 学部 代数学 。 オックスフォード大学出版局 。pp.360-361。ISBN 0-19-853248-2 。
^ メイソン、マシュー T. (2001)。ロボット操作の力学。MIT プレス 。ISBN 978-0-262-13396-8 . 2014年 11月21日 閲覧 。
外部リンク
ウィキメディア・コモンズには、慣性モーメント に関連するメディアがあります 。
2次元および3次元における角運動量と剛体回転
剛体の回転と慣性モーメントに関する講義ノート
慣性モーメントテンソル
慣性モーメントの入門レッスン: 垂直のポールが倒れないようにする (Java シミュレーション)
慣性モーメントの求め方に関するチュートリアル。さまざまな基本図形の問題と解答付き。
操作の力学に関するメモ: 角慣性テンソル