マッカラーの公式のソースを表示

{{翻訳直後|[[:fr:Formule de MacCullagh]] 29 décembre 2016 à 15:13 (CEST)|date=2019年6月}}

[[File:James MacCullagh.png|thumb|right|150px|[[ジェームズ・マッカラー]] （[[1809年]] - [[1847年]]）]] 
'''マッカラーの公式'''（マッカラーのこうしき、{{Lang-fr-short|formule de MacCullagh}}）は、球状でない物体が、物体から十分に離れてはいるが、物体の径と比べれば極端には離れていない点で作る[[重力場|万有引力の場]]を記述するのに重要な公式である。万有引力ポテンシャルの基本的な定義が基になっているが、ポテンシャルの一般的なラプラス展開をするのではなく、2次近似までで切ることで得られる。[[測地学]]の理論や[[地球物理学]]で広く用いられる。

== 万有引力ポテンシャルの一般表示 == 
[[File:MacCullagh 1.png|thumb|300px]]

[[球面座標系]]で、ポテンシャルが作られる点 <math>P</math>（位置 <math>\mathbf{r}</math>）、ポテンシャルを作る点 <math>P'</math>（位置 <math>\mathbf{r'}</math>）の動径成分をそれぞれ <math>r</math> 、<math>r'</math> とする。座標の原点 <math>O</math> は当面は任意としておく。ベクトル <math>\mathbf{r}</math> と <math>\mathbf{r'}</math> が成す角を <math>\psi</math> 、<math>P'</math> での体積要素を <math>\mathrm dM(\mathbf{r'})</math> 、物体Bの全質量を <math>M</math> と書く。このとき、物体Bによって作られる万有引力ポテンシャルは一般的に

:<math>V(P) = G \int_{M} \frac{1}{\|\mathbf{r'} - \mathbf{r}\|}\, \mathrm dM(\mathbf{r'})</math>

と書ける<ref>ここでは <math>V</math> を仕事として定義しているが、物理学（の他分野）では、万有引力（重力）ポテンシャルは <math>\Phi = -V</math> を満たす相互作用エネルギー <math>\Phi</math> と定義されるのが慣習である。</ref>。ここで  

:<math>\|\mathbf{r'} - \mathbf{r}\| = r \sqrt{1 - 2\frac{r'}{r}\cos \psi + \left(\frac{r'}{r}\right)^{2}} </math>

で、<math>G</math> は[[アイザック・ニュートン|ニュートン]]の[[万有引力定数]]である。

<math>-1 \le x \le +1</math> のとき <math>(1 + x)^{-1/2} = 1 - x/2 + 3x^2/8 + o(x^3)</math> だから<ref><math>\textstyle\frac{1}{\sqrt{1 - 2 \frac{r'}{r} \cos \psi + \left(\frac{r'}{r}\right)^2}} = 1 - \frac12\left[- 2 \frac{r'}{r} \cos \psi + \left(\frac{r'}{r}\right)^2\right] + \frac{3}{8}\left[- 2 \frac{r'}{r} \cos \psi + \left(\frac{r'}{r}\right)^2\right]^2 + \dots</math> <math>\textstyle = 1 + \frac{r'}{r} \cos \psi - \frac12 \left(\frac{r'}{r}\right)^2 + \frac38 \left[4 \left(\frac{r'}{r}\right)^2 \cos^2\psi - 4 \left(\frac{r'}{r}\right)^3 \cos \psi + \left(\frac{r'}{r}\right)^4\right] + \dots </math><math>\textstyle = 1 + \frac{r'}{r} \cos \psi - \frac{1}{2} \left(\frac{r'}{r}\right)^2 + \frac{3}{2} \left(\frac{r'}{r}\right)^2 \cos^2\psi + o\left[\left(\frac{r'}{r}\right)^3\right]</math></ref>、<math>x = - 2(r'/r) \cos \psi + (r'/r)^2</math> とおくと 

:<math>\frac{1}{\sqrt{1 - 2 \frac{r'}{r} \cos \psi + \left(\frac{r'}{r}\right)^2}} = 1 + \frac{r'}{r} \cos \psi - \frac{1}{2} \left(\frac{r'}{r}\right)^2 + \frac{3}{2} \left(\frac{r'}{r}\right)^2 \cos^2\psi + o\left[\left(\frac{r'}{r}\right)^3\right]</math>

よって、

:<math>\frac{r}{\|\mathbf{r'} - \mathbf{r}\|} = 1 + \frac{r'}{r} \cos \psi + \left(\frac{r'}{r}\right)^2 - \frac32 \left(\frac{r'}{r}\right)^2 \sin^2\psi + \dots</math>

が得られる<ref><math>\textstyle\frac{r}{\|\mathbf{r'} - \mathbf{r}\|} = 1 + \frac{r'}{r} \cos \psi + \frac{1}{2} \left(\frac{r'}{r}\right)^2 (3 \cos^2\psi - 1) + \dots = 1 + \frac{r'}{r} \cos \psi + \frac{1}{2} \left(\frac{r'}{r}\right)^2 (3 - 3 \sin^2\psi - 1) + \dots </math> <math>\textstyle = 1 + \frac{r'}{r} \cos \psi + \left(\frac{r'}{r}\right)^2 - \frac{3}{2} \left(\frac{r'}{r}\right)^2 \sin^2\psi + \dots</math> .</ref> 
。これより <math>P</math> での万有引力ポテンシャルは次のようになる。

:<math>V(P) = \frac Gr \int_M \mathrm dM + \frac G{r^2} \int_M r' \cos \psi \mathrm dM + \frac G{r^3} \int_M r'^2 \mathrm dM - \frac{3G}{2r^3} \int_M r'^2 \sin^2\psi \mathrm dM + \dots</math>

== マッカラーの定理 == 

このように万有引力ポテンシャルの2次までの近似 <math>V(P)</math> は、<math>V_0(P)</math>, <math>V_1(P)</math>, <math>V_2(P)</math>, <math>V_3(P)</math> の4つの項に分解できることがわかる。つまり

:<math>V(P) = V_0(P) + V_1(P) + V_2(P) + V_3(P) .</math>

第1項

:<math>V_0 = G r^{-1} \int_M\, \mathrm dM = GM r^{-1}</math>

は球対称な質量分布のポテンシャルに対応する。これは物体の全質量が1点に集められている場合と同じである。減少の度合 1/r が他の項と比べ緩やかであるため、非常に遠距離ではこの項が支配的となる。これは[[ケプラーの法則|惑星のケプラー軌道]]を確立するのに用いられた、古典的なニュートンポテンシャルである。'''単極子'''（monopôle）、または単極子的質量分布（distribution de masse monopolaire）と呼ばれる。

[[File:MacCullagh 2.png|thumb|right|300px]] 
第2項

:<math>V_1 = G r^{-2} \int_M r' \cos \psi\, \mathrm dM</math>

は'''双極子'''（dipôle）、または双極子的質量分布（distribution de masse dipolaire）に対応する。座標の原点を物体の質量中心に一致させると、この項は消える<ref>このように、（近似を続けつつ）万有引力ポテンシャルの双極子項が消えるように座標軸の取り方を決めることで、地球の質量中心を正確に決定できる。</ref>。実際、右の図を参照すると次のように書くことができる。

:<math>V_1 = G r^{-2} \int_M r' \cos \psi\, \mathrm dM = G r^{-2} \int_M x\, \mathrm dM = G M r^{-2} x_0</math>

ここで <math>x_0</math> は物体の質量中心の Ox-軸成分。

第3項と第4項、<math>V_2</math>, <math>V_3</math> は'''四極子'''（quadrupôle）、または四極子的質量分布（distribution de masse quadrupolaire）である。<math>V_2</math> の計算を続けると   

:<math>
\begin{align}
V_2 &= G r^{-3} \int_M r'^2 \, \mathrm dM = G r^{-3} \int_M (x^2+y^2+z^2) \mathrm dM \\
&= \frac12 G r^{-3} \int_M [(y^2+z^2)+(x^2+z^2)+(x^2+y^2)] \, \mathrm dM = \frac12 G r^{-3} (A+B+C)
\end{align}
</math>

となる。ここで <math>A</math>, <math>B</math>, <math>C</math> はそれぞれOx, Oy, Oz-軸のまわりの[[慣性モーメント]]を表す。慣性モーメントの平均 <math>\bar I</math> を 

:<math>\bar I = \frac13 (A+B+C)</math>,</center>

とおくと、

:<math>V_2 = \frac32 G r^{-3} \bar I</math>

となる。<math>V_3</math> に対しては、

:<math>V_3 = -\frac32 G r^{-3} \int_M r'^2 \sin^2 \psi\, \mathrm dM = -\frac32 G r^{-3} \int_M r'^2 \cos^2 \alpha\, \mathrm dM = -\frac32 G r^{-3} I</math> 

ここで <math>I</math> は <math>\overline{OP}</math>-軸のまわりの慣性モーメントである。これより、次のマッカラーの公式（マッカラーの定理とも）が得られる。

:<math>V(P) = G M r^{-1} - \frac32 G r^{-3} (I - \bar I)</math>

この公式は、地球や他の惑星等の球状に近い物体に適用するのに非常に便利である。物体に何らかの対称性があれば、近似は十分遠距離でも有効である。明らかに、2次までの近似であることから有効な適用範囲には限界があり、外部重力場をさらに正確に表すためにはより高次の調和関数が必要になる。

== 対称性のある物体に対するマッカラーの公式 == 

質量中心を通る軸 <math>\overline{OP}</math> のまわりの慣性モーメントは、慣性テンソルを表す 3×3行列を[[対角化]]することで得られる主慣性モーメント <math>A</math>, <math>B</math>, <math>C</math> を使って

:<math>I = n_x^2 A + n_y^2 B + n_z^2 C</math> 

と書くことができる。ここで <math>n_x</math>, <math>n_y</math>, <math>n_z</math> は <math>\overline{OP}</math> （単位ベクトル）の慣性主軸 Ox, Oy, Oz への方向余弦を表す。

物体が Oz-軸を中心に回転する（Oz について軸対称である）と仮定する。このとき <math>A = B</math> である。<math>\phi</math> を <math>\overline{OP}</math> が Oxy-平面となす角とする<ref>地球物理学の文脈では、 <math>\phi</math> は一般に地点の緯度（もしくは緯度にほぼ一致する角）を表す。</ref>。このとき <math>n_z = \sin \phi</math> であり、<math>n_x^2 + n_y^2 + n_z^2 = 1</math> だから方向余弦は
 
:<math>\sin^2 \phi = 1 - n_x^2 - n_y^2</math>

である。<math>A = B</math> としているので、マッカラーの公式は

:<math>V(P) = GM r^{-1} - \frac12 G r^{-3} (3 I - 2 A - C)</math> 

となる。慣性モーメント I は

:<math>I = A (n_x^2+n_y^2) + C \sin^2 \phi = A (1 - \sin^2 \phi) + C \sin^2 \phi</math> 

となるから、回転体に対するマッカラーの公式は最終的に次の形になる
<ref><math>V(P) = GM r^{-1} - \frac12 G r^{-3} (3 I - 2 A - C)</math>
<math>= GM r^{-1} - \frac12 G r^{-3} (3 A - 3 A \sin^2 \phi + 3 C \sin^2 \phi - 2 A - C) </math>
<math>= GM r^{-1} - \frac12 G r^{-3} [A - C - 3 (A - C) \sin^2 \phi]</math>
<math>= GM r^{-1} + \frac12 G r^{-3} (C - A) (1 - 3 \sin^2 \phi)
</math></ref>。

:<math>V(P) = GM r^{-1} + \frac12 G r^{-3} (C - A) (1 - 3 \sin^2 \phi)</math>

== 地球の力学的形状係数 == 

軸対称性があるとき、万有引力ポテンシャルの多重極展開を2次まで行ったものは常に次の一般式で書ける。 

:<math>
\begin{align}
V(P) &= GM r^{-1} - J_2 GM a^2 r^{-3} P_2 (\cos \theta) \\
&= GM r^{-1} - \frac12 J_2 GM a^2 r^{-3} (1 - 3 \sin^2 \phi)
\end{align}
</math>

ここで 

:<math>P_2 (\cos \theta) = \frac32 \sin^2 \phi - \frac12 ,\quad (\phi = \frac{\pi}{2} -  \theta) </math>

は2次の[[ルジャンドル多項式]]である。地球の場合、2次の帯調和重力係数 <math>J_2</math> はしばしば地球の力学的形状係数（facteur de forme géodynamique）と呼ばれる。この式をマッカラーの公式と等しいと置くと、

:<math>J_2 = \frac{C - A}{Ma^2}</math>

この最後の関係式もマッカラーの公式と呼ばれることがある。

== 注釈 == 
<references/>

{{Earth-sci-stub}}
{{Physics-stub}}

{{DEFAULTSORT:まつからあのこうしき}}
[[Category:地球物理学]]
[[Category:ポテンシャルエネルギー]]
[[Category:人名を冠した数式]]