球殻の慣性モーメント
📂古典力学球殻の慣性モーメント
公式
半径がa、質量がmの球殻の慣性モーメントは以下のようになる。
I=32ma2
導出
半径がaで質量がmの均一な球殻を考えてみよう。球の慣性モーメントを求める時と同じアイデアを使う。しかし、少し違いがある。球の慣性モーメントを計算する時のように、多数の円筒殻の合計として球殻を考えよう。

でもここで球の場合と同じように計算すると問題が生じる。球の場合、小さな円盤の体積を積分すると球の体積が出る。
∫−rrπ(r2−z2)dz=∫−rrπr2dz−∫−rrπz2dz=πr2[z]−rr−π[31z3]−rr=2πr3−32πr3=43πr3
しかし、小さな円盤の表面積の和の極限が球の表面積と同じではない。問題はまさにここで生じる。実際に以下の計算でこれを確認できる。
∫−rr2πxdz=∫−rr2πr2−z2dz
この時z≡rsinθに置き換えると
∫−rr⟹∫−2π2π&dz=rsinθdθ
∫−2π2π2πr2−r2sin2θrcosθdθ=2πr∫−2π2πr21−sin2θcosθdθ=2πr2∫−2π2πcos2dθ=πr2∫−2π2π(1+cos2θ)dθ=πr2[θ+21sin2θ]2π2π=πr2(π)=π2r2<4πr2
なので球の表面積に足りない。だから、球の表面積を円盤の表面積の和で近似することはできない。
次に、小さな円筒の表面積を少し異なる方法で近似してみよう。円筒の高さを上面と下面の垂直距離ではなく、側面の距離として置いて解こう。このように置いて積分すると、円筒の表面積が出る。興味があれば、自分でやってみて。

質量が一様にmで、半径がrの円筒殻の慣性モーメントはI=mr2なので、球殻の慣性モーメントは以下のようになる。
Isphere shell=∫dl=∫x2dm
上の図から得られたdmの値を代入すると、以下を得る。
∫x2ρ2πxadθ
この時、θによる積分をzによる積分に変えるために、以下の関係式を使って青色部分を上手く変えると以下のようになる。
z=asinθ⟹dz=acosθdθ,x=acosθ
わざわざzについて積分する理由はθについて積分するよりも簡単だからだ。
∫(a2−z2)ρ2πacosθadθ=ρ2πa∫−aa(a2−z2)dz=ρ2πa[a2z−31z3]−aa=ρ2πa(2a3−32a3)=ρ38πa4
そして、球殻の質量がm=ρ4πa2でρ=4πa2mなので、以下を得る。
Iz=ρ38πa4=4πa2m38πa4=32ma2
球と同じように、球殻もまたすべての方向から対称なので、以下が真だ。
Ix=Iy=Iz=32ma2
■