![[Cheat Physics]](image/cheat_logo.png "[Cheat Physics]"){kind=logo}
>>「力学」の平原
万有引力をクリア後に選択可能になります。
かなり早いうちに攻略可能となりますが、難易度が高めな上に、これを攻略しないと出現しないというクエストもほとんどないので後回しで構いません。
3 次元の極座標では 3 つのパラメータ \( (r,\theta,\phi) \) を使って位置を表します。
デカルト座標と極座標との関係は次のようになっています。
\[ \begin{align*} x \ &=\ r \, \sin\theta \, \cos\phi \\ y \ &=\ r \, \sin\theta \, \sin\phi \tag{1} \\ z \ &=\ r \, \cos\theta \end{align*} \]
ある点 \( P(r,\theta,\phi) \) を始点にして、それぞれのパラメータを別個に微小変化させてやると、それぞれ直交する方向に \( \diff r \) 、\( r \diff \theta \)、\( r \sin \theta \diff \phi \) だけ移動します。 この 3 辺で作られる微小な直方体の体積を次のように表すことができます。
\[ \diff V \ =\ r^2 \sin \theta \diff r \diff \theta \diff \phi \tag{2} \]
大きな球体はこの微小な直方体を使って埋め尽くすことが出来るので、 うまくパラメータの範囲を考えて微小体積を全て合計してやれば球全体の体積 \( V \) を求めることもできます。
\[ \begin{align*} V \ &=\ \int \diff V \\ &=\ \tint r^2 \sin \theta \diff r \diff \theta \diff \phi \\ &=\ \int^R_0 r^2 \diff r \int_0^\pi \sin \theta \diff \theta \int_0^{2\pi} \diff \phi \\ &=\ \frac{4}{3} \pi R^3 \tag{3} \end{align*} \]
天体の密度 \( \rho \) は内部まで一定ではありませんが、おそらくは球対称な質量分布になっているでしょうから \( \rho(r) \) と表せると仮定します。 微小体積 \( \diff V \) と密度 \( \rho(r) \) の積を作れば微小直方体の微小質量 \( \diff M \) が得られます。
\[ \diff M \ =\ \rho(r) \diff V \tag{4} \]
これを積分すれば、天体の全質量 \( M \) を求めることもできます。
\[ \begin{align*} M \ &=\ \int \diff M \\ &=\ \int \rho(r) \diff V \\ &=\ \tint \rho(r) \, r^2 \sin \theta \diff r \diff \theta \diff \phi \\ &=\ \int_0^R \rho(r)\, r^2 \diff r \ \int_0^\pi \sin \theta \diff \theta \ \int_0^{2\pi} \!\! \diff \phi \\ &=\ 4\pi \, \int_0^R \rho(r)\, r^2 \diff r \tag{5} \end{align*} \]
しかし密度 \( \rho(r) \) が具体的に定まっていないのでこれ以上は計算を進められません。 この (5) 式は後で使うことになります。
さて、この天体の中心から距離 \( s \) だけ離れたところに質量 \( m \) の質点を置いたときにどれだけの引力が働くかを計算してみましょう。 そのためには天体を構成する微小体積 \( \diff V \) が持つ質量 \( \diff M \) との間に働く力の大きさ \( \diff F \) を求めて、その全てを合計してやればいいわけです。
\[ \begin{align*} \diff F \ &=\ -G \frac{m \, \diff M}{d^2} \\[3pt] &=\ -G \frac{m \, \rho(r) \diff V}{d^2} \tag{6} \end{align*} \]
ここで使った \( d \) は微小質量 \( \diff M \) と質点 \( m \) との直線距離を表しています。 今からそれを具体的に求めてみましょう。 この質点 \( m \) をデカルト座標の \( (0,0,s) \) の点に置きます。 それぞれの微小質量 \( \diff M \) がデカルト座標 \( (x,y,z) \) の位置にあるとすると、距離は次のように表せます。
\[ \begin{align*} d^2 \ &=\ x^2 + y^2 + (z - s)^2 \\ &=\ x^2 + y^2 + z^2 - 2sz + s^2 \\ &=\ r^2 - 2sz + s^2 \tag{7} \end{align*} \]
ここに (1) 式を使えば (7) 式に含まれている \( z \) を極座標に書き換えることが出来て、距離 \( d \) は次のように全て極座標の変数で表せます
\[ d \ =\ \sqrt{r^2+s^2 - 2s r\,\cos\theta } \tag{8} \]
質点 \( m \) をわざわざ \( z \) 軸上に置いたのは、実はこの式に変数 \( \phi \) が混じらないための工夫なのでした。 なるべく単純な形にしておいたほうが後の計算が楽になるからです。
さて、これで微小な力の大きさ \( \diff F \) を求めることが出来るようになったのですが、このまま単純に合計するわけには行きません。 それぞれの微小質量からの力はそれぞれに向きが違うので、力どうしの打ち消し合いなどを考慮する必要があります。 そのためには成分ごとに分けて計算しなければなりません。 ところが今、質点は \( z \) 軸上に置いてあり、天体の形は \( z \) 軸に対称に分布しているので、\( x \) 軸方向の力も \( y \) 軸方向の力も完全に打ち消し合って 0 になることが、わざわざ計算をしなくても分かってしまいます。 合計の力はまっすぐ \( z \) 軸の上に乗るような方向を向いているはずです。 そこで、\( z \) 軸成分だけを求めて合計してやればいいわけです。
\( \diff M \) と \( m \) との直線距離が \( d \) で、その距離の \( z \) 成分が \( s - z \) ですから、この比率を考えて (6) 式に掛けてやれば \( \diff F \) の \( z \) 成分の大きさが計算できます。
\[ \begin{align*} \diff F_z \ &=\ -G \, \frac{m \, \rho(r) \diff V}{d^2} \cdot \frac{s-z}{d} \\[3pt] &=\ -G \, \frac{m \, \rho(r) \diff V}{r^2+s^2 - 2s r\,\cos\theta} \cdot \frac{s-r\,\cos\theta}{\sqrt{r^2+s^2 - 2s r\,\cos\theta}} \\[3pt] &=\ -G\,m \, \frac{\rho(r) \, (s-r\,\cos\theta)}{(r^2+s^2 - 2s r\,\cos\theta)^{3/2}} \diff V \\[3pt] &=\ -G\,m \, \frac{\rho(r) \, (s-r\,\cos\theta)}{(r^2+s^2 - 2s r\,\cos\theta)^{3/2}} \ r^2 \sin\theta \diff r \diff \theta \diff \phi \\ &=\ -G\,m \, \frac{\sin\theta \, (s-r\,\cos\theta)}{(r^2+s^2 - 2s r\,\cos\theta)^{3/2}} \ \rho(r) \, r^2 \diff r \diff \theta \diff \phi \tag{9} \end{align*} \]
あとはこれを積分してやるだけです。
\[ F_z \ =\ -G\,m \tint \frac{\sin\theta \, (s-r\,\cos\theta)}{(r^2+s^2 - 2s r\,\cos\theta)^{3/2}} \ \rho(r) \, r^2 \diff r \diff \theta \diff \phi \tag{10} \]
とは言ったものの簡単ではありません。 \( \phi \) についての積分は独立しているのですぐに計算できますが、\( r \) と \( \theta \) が入り混じっているため、 一方の変数で積分してからその結果を残りの変数で積分してやる必要があります。 \( \rho(r) \) が具体的に定まっていないので \( r \) での積分は後回しにして、このまま \( \theta \) での積分を先にやってしまいましょう。 式変形がややこしくなるので、\( \theta \) での積分に関係のある部分だけを取り出して計算を進めます。 この計算は部分積分を使えばうまく行くということを知っていないと難しいと思います。
\[ \begin{align*} &\int^\pi_0 \frac{\sin\theta}{(r^2+s^2 - 2s r\,\cos\theta)^{3/2}} (s-r\,\cos\theta) \diff \theta \\[3pt] =\ &\left[ -\frac{1}{sr} \, \frac{1}{\sqrt{r^2+s^2 - 2s r\,\cos\theta}} \, (s-r\,\cos\theta) \right]_0^\pi \ -\ \int_0^\pi \left( -\frac{1}{sr} \, \frac{1}{\sqrt{r^2+s^2 - 2s r\,\cos\theta}} \right)(r\,\sin\theta) \diff \theta \\[3pt] =\ &-\frac{1}{sr} \left[ \frac{s-r\,\cos\theta}{\sqrt{r^2+s^2 - 2s r\,\cos\theta}} \right]_0^\pi \ +\ \frac{1}{s} \int_0^\pi \frac{\sin\theta}{\sqrt{r^2+s^2 - 2s r\,\cos\theta}} \diff \theta \\[3pt] =\ &-\frac{1}{sr} \left[ \frac{s-r\,\cos\theta}{\sqrt{r^2+s^2 - 2s r\,\cos\theta}} \right]_0^\pi \ +\ \frac{1}{s} \left[ \frac{1}{sr} \sqrt{r^2+s^2 - 2s r\,\cos\theta} \right]_0^\pi \\[3pt] =\ &-\frac{1}{sr} \left( \frac{s+r}{\sqrt{r^2+s^2 + 2s r}} - \frac{s-r}{\sqrt{r^2+s^2 - 2s r}} \right) \ +\ \frac{1}{s^2 r} \left( \sqrt{r^2+s^2 + 2sr} - \sqrt{r^2+s^2 - 2sr}\right) \\[3pt] =\ &-\frac{1}{sr} \left( \frac{s+r}{\sqrt{(r+s)^2}} - \frac{s-r}{\sqrt{(r-s)^2}} \right) \ +\ \frac{1}{s^2 r} \left( \sqrt{(r+s)^2} - \sqrt{(r-s)^2} \right) \tag{11} \end{align*} \]
この根号を外すときに少し注意が必要です。 根号の中が 2 乗になっているので正の値であることは保証されています。 根号を外したときにも正の値になるようにしておかなくてはなりません。 今、質点 \( m \) を置いた地点は天体の外側にあるので、\( s>r \) という条件があります。 よって、\( \sqrt{(r-s)^2} \) と書かれている部分は \( (r-s) \) ではなく \( (s-r) \) のように書き直してやれば良いでしょう。
\[ \begin{align*} =\ &-\frac{1}{sr} \left( \frac{s+r}{s+r} - \frac{s-r}{s-r} \right) \ +\ \frac{1}{s^2 r} \Big( (s+r) - (s-r) \Big) \\ =\ &-\frac{1}{sr} \left( 1 - 1 \right) \ +\ \frac{1}{s^2 r} (2r) \\ =\ &\frac{2}{s^2} \tag{12} \end{align*} \]
この結果を (10) 式に戻してやりましょう。
\[ \begin{align*} F_z \ &=\ -G\,m \, \frac{2}{s^2} \int_0^R \rho(r)\, r^2 \diff r \int^{2\pi}_0 \diff \phi \\[3pt] &=\ -G\,m \, \frac{4\pi}{s^2} \int_0^R \rho(r)\, r^2 \diff r \tag{13} \end{align*} \]
ここで (5) 式を当てはめてやると、積分は天体の全質量 \( M \) で置き換えることができます。
\[ F_z \ =\ -G \frac{Mm}{s^2} \tag{14} \]
こうして、天体の質量分布が球対称であるような場合には、質量 \( m \) から距離 \( s \) だけ離れた地点に全質量が集中しているのと同じであると近似してよいということが証明できました。
ついでですから、質点 \( m \) が天体の内部にある場合についても考えてみましょう。 この場合、天体の微小質量 \( \diff M \) の位置によっては、\( r>s \) が成り立つ場合も出てきます。 原点を中心にした半径 \( s \) の球の外側に位置する微小質量 \( \diff M \) はどれもこの条件を満たします。 この場合には (11) 式の根号を外す時に正の値になるようにするために \( \sqrt{(r-s)^2} \) を \( (r-s) \) のように置き換えねばなりません。
\[ \begin{align*} (11)式 \ &=\ -\frac{1}{sr} \left( \frac{s+r}{s+r} - \frac{s-r}{r-s} \right) \ +\ \frac{1}{s^2 r} \left( (s+r) - (r-s) \right) \\ &=\ -\frac{1}{sr} \left( 1 - (-1) \right) \ +\ \frac{1}{s^2 r} (2s) \\ &=\ -\frac{2}{sr} \ +\ \frac{2}{sr} \\ &=\ 0 \tag{15} \end{align*} \]
きれいに 0 になってしまいました。 (10) 式に戻って計算する時に、\( r>s \) の条件を満たす部分だけは (15) 式の結果を使って、\( r<s \) の条件を満たす部分については (12) 式の結果を使うことになりますから、天体内部に置かれた質点 \( m \) は、 それより深いところに位置する球体領域からのみ引力を受けるのだと言えます。
これは質点 \( m \) が置かれた位置より外側の球殻からの引力は全てきれいに打ち消し合うことを意味しています。 もし天体の内部が中心部からきれいに球形にくり抜かれていれば、その巨大な空洞の内部ではどこであっても一切、天体からの引力を感じることがないわけです。
（要編集：ヒントが載っている魔導書ページへのリンク、記載されているページ数など）
calculus4d.png
calculus4c.png
