1. はじめに

1
2013 年 01 月 06 日

地球の重力について
新潟工科大学 情報電子工学科 竹野茂治

1

はじめに

本稿では、3 重積分の極座標変換の例として、地球の重力の計算を紹介する。そしてそれにより、ついでに次のようないくつかの疑問にも答えたい。

• 疑問 1: 地球のような星から受ける万有引力は、本来は星の各部分から受ける小さい万有引力の合力だが、その合力は、星の中心 (重心) 一点に質量が集中していると見た場合の引力と見てよいのかどうか。
• 疑問 2: 地球の外部ではなく、地球の内部では重力はどのようにかかるのか。すなわち、中心に進むにつれて重力は大きくなるのか小さくなるのか。
• 疑問 3: 地球の内部は実は空洞で、その中に地底人が住む、という SF があるが、彼らが受ける重力はどれくらいか。

123

図 1: 3 つの疑問

2

ある物体から受ける万有引力

まず、形と質量のある物体 V から、点 A に置いた質量 m の小さい物体が受ける万有引力について考えてみる。

2. ある物体から受ける万有引力

2V が、点 P に置かれた質量 M の小さい物体である場合、点 A はその物体から、距離の 2 乗に反比例し、それぞれの質量に比例して、互いに引き合う力である 万有引力−→mM G APF = −→ · −→|AP|2 |AP|−→ −→を受ける (G > 0: 万有引力定数)。なお、AP/|AP| は A から P に向かう単位ベクトルであることに注意する。
複数の点 Pj に質量 Mj の小さい物体がある場合は、その合力−−→APj
F =−−→ 2 · −−→|APj |j=1 |APj |n∑mMj Gを受けることになる。

3 次元領域 V の場合は、V 内の各点 P での密度を ψ(P) とすると、V を小片 ∆Vj に分解して考えれば、その小片の質量はほぼ ψ(Pj )∆Vj なので (Pj ∈ ∆Vj )、V から A の物体が受ける万有引力は、
n∑mGψ(Pj )∆VjF =j=1−−→|APj |2−−→APj· −−→|APj |にほぼ等しく、厳密にはその極限としての体積分 (3 重積分) で表されることになる (積−→分変数は x = OP):−→−→∫∫∫mGψ(P) AP
x − OAF =· −→ dv = mG
ψ(x)−→−→ dxdydz
V
V|AP|2|AP||x − OA|3∫(1)

この公式自体は、A が V の内部にあっても変わらないが、その場合は分母が 0 となりうるので、厳密には (1) は広義積分となる。
さて、V の重心点 Q の位置ベクトルは、n∑−−→ψ(Pj )∆Vj OPjj=1n∑
j=1ψ(Pj )∆Vj

3. 地球のような星の場合3の極限として得られるので、V の総質量を MV とすれば、(1 ∫−→−→OQ =ψ(P)OP dvMV V)∫MV =ψ(P)dvVと表される。よって、∫V−→ −→
ψ(P)(OP − OA) dv =
∫

=V

∫
V−→−→
ψ(P)OP dv − OA→−−→ψ(P) O dv − MV OA∫ψ(P) dv
Vより、もし重心 Q に質量 MV が集中していると考えると、それが A におよぼす万有引力は、(∫)

1−→−→−→mMV G
ψ(P)OP dv − OA
mMV G AQMV VF =¯¯∫−→ 2 · −→ =¯1−→−→¯¯3|AQ||AQ|¯
ψ(P)OPdv−OA¯¯MV V∫−→mG ψ(x)(x − OA) dv
V= ¯¯3∫¯ 1¯−→¯ψ(x)(x − OA) dv ¯¯¯MV Vとなるが、これは (1) とは明らかに異なり、一般に両者は等しくはない。
よって、重心に質量が集中していると考えてよいのは、特別な場合であることがわかる。

地球のような星の場合

3

地球のような星の場合、V は球で1 、さらに、その質量は中心から層状に対称に分布すると考えられるので、地球の中心を原点 O とすれば、ψ(x) は r = |x| のみの関数

ψ(x) = ρ(r)と見ることができる2 。この条件の元で (1) を極座標変換を用いて計算してみる。

12

厳密には地球は球ではなく少しゆがみがある。
このような層状の対称性を 球対称 と呼ぶことがある。

3. 地球のような星の場合

4V は半径 R の球であるとし、極座標を
x = (r cos φ cos θ)e1 + (r sin φ cos θ)e2 + (r sin θ)e3(0 ≤ r ≤ R, 0 ≤ φ < 2π, −π/2 ≤ θ ≤ π/2)(2)のように導入する。ここで、e1 , e2 , e3 は互いに直交する単位ベクトルとするが、軸方−→向の基本ベクトルと取る必要はないので、計算を簡単にするため、e3 を OA の方向に取り、e1 , e2 はそれに応じて取り直すこととする。ただし、A = O の場合は後で別に−→−→−→考えることとし、まずは |OA| =6 0 の場合を考える。この場合は、OA = |OA|e3 となることに注意する。
まずはこの変数変換のヤコビアン D(x)/D(r, φ, θ) であるが、これはD(x)D(r, φ, θ)¯¯¯ cos φ cos θ −r sin φ cos θ −r cos φ sin θ ¯¯¯¯¯¯= ¯ sin φ cos θ r cos φ cos θ −r sin φ sin θ ¯¯¯¯¯¯sin θ0r cos θ¯¯¯ cos φ cos θ − sin φ − cos φ sin θ ¯¯¯¯¯2¯= r cos θ ¯ sin φ cos θ cos φ − sin φ sin θ ¯¯¯¯¯¯sin θ0cos θ= r2 cos θ(cos2 φ cos2 θ + sin2 φ sin2 θ + cos2 φ sin2 θ + sin2 φ cos2 θ)= r2 cos θ (≥ 0)となる。また、−→−→
x − OA = (r cos φ cos θ)e1 + (r sin φ cos θ)e2 + (r sin θ − |OA|)e3
より、

−→
−→
|x − OA|2 = r2 cos2 φ cos2 θ + r2 sin2 φ cos2 θ + (r sin θ − |OA|)2−→−→= r2 − 2r|OA| sin θ + |OA|2となるので、(1) は極座標 (2) により以下のようになる:−→(r cos φ cos θ)e1 + (r sin φ cos θ)e2 + (r sin θ − |OA|)e3
F = mG
ρ(r)−→−→W(r2 − 2r|OA| sin θ + |OA|2 )3/2×r2 cos θdrdφdθ(3)∫∫∫3. 地球のような星の場合

5

ここで、W = {(r, φ, θ); 0 ≤ r ≤ R, 0 ≤ φ < 2π, −π/2 ≤ θ ≤ π/2} であり、
∫ 2π

∫ 2πcos φ dφ =sin φ dφ = 000なので、F の e1 , e2 方向成分は 0 となり、よって e3 成分のみ残ることになるが、これは φ にはよらないので、∫ R

F = 2πmGe30−→
ρ(r)r2 (r sin θ − |OA|) cos θdrdθ−→−→−π/2 (r 2 − 2r|OA| sin θ + |OA|2 )3/2
∫ π/2(4)となる。今、f (t) (t > 0) を、∫ π/2f (t) =(t sin θ − 1) cos θdθ−π/2 (t2 − 2t sin θ + 1)3/2(5)とすると、(fr−→|OA|)−→= |OA|2−→(r sin θ − |OA|) cos θdθ−→−→−π/2 (r 2 − 2r|OA| sin θ + |OA|2 )3/2∫ π/2なので、(4) は f (t) を用いて2πmG ∫ R
F = −→ e3
ρ(r)r2 f20|OA|()r−→ dr|OA|(6)と書けることになる。

f (t) を積分するために u = t2 − 2t sin θ + 1 と置換すると、t sin θ − 1 =t2 + 1 − ut2 − 1 − u−1=,22cos θdθ = −1
du2tより、t 6= 1 のとき f (t) はf (t) =∫ t2 −2t+1

t2 − 1 − u
2u3/2
t2 +2t+1[(−12t])du =21 ∫ (t+1)2 2{(t − 1)u−3/2 − u−1/2 }du
24t (t−1)[]2

u=(t+1)√ u=(t+1)
2
1
1 t2 − 1 + u√−(t2 − 1) √ − 2 u==−4tu2tu
u=(t−1)2u=(t−1)2

3. 地球のような星の場合

6(1 t2 − 1 + (t + 1)2 t2 − 1 + (t − 1)2= −−2tt+1|t − 1|()22
1 2t + 2t 2t − 2tt−1= −−=−12tt+1|t − 1||t − 1|{=)0(t > 1 のとき),−2 (0 < t < 1 のとき)と計算される。よって、これを (6) に代入すると、結局4πmG ∫ R∧|OA|
ρ(r)r2 dr (a ∧ b = min{a, b})
F = − −→ e3
20|OA|−→(7)と書けることになる。
さらに、地球の中心から半径 s (> 0) までの部分の球の質量を M̂ (s) と書くことにすると、∫∫ sM̂ (s) =ρ(r)dv =0<r<s∫ π/2dr0−π/2∫ 2πdθρ(r)r2 cos θdφ = 4π0∫ sρ(r)r2 dr0となるので、これを用いれば (7) は、−→
mM̂ (R ∧ |OA|)G
F =−
e3−→
|OA|2(8)と書くこともできる。

−→−→
なお、A=O の場合は (r2 − 2r|OA| sin θ + |OA|2 )3/2 = r3 となるので、(3) は、∫∫∫

F = mGW∫ R= = mG0

ρ(r)(e1 cos φ cos θ + e2 sin φ cos θ + e3 sin θ) cos θdrdφdθρ(r)dr · (0e1 + 0e2 + 0e3 ) = 0となる。しかも、この計算は、A ∈ V の場合の広義積分が発散せずに収束することも示していることに注意する。A ∈ V の場合、A を中心とする小さな球 (A の近傍) の積分は、A を中心とする極座標で積分すれば、ψ(x) が ρ(r) ではなく r, φ, θ に依存する関数となるので 0 にはならないが、それ以外は上の計算とほぼ同じで分母の r は分子の r と丁度約分されてしまうので、積分は発散せずに収束することがわかる。

4. 考察47考察

本節では、公式 (8) を用いて、最初にあげたいくつかの疑問について考察してみる。

−→
まず、A が V の外にある場合は、|OA| > R なので−→M̂ (R ∧ |OA|) = M̂ (R) = MVとなり、よって引力はmMV G
F = − −→ e3|OA|2となるが、これは原点一点に質量 MV が集中している場合の引力と等しい。つまり、
V が球で、その密度分布が球対称な場合には、それが外部におよぼす引力は、質量が中心一点に集中していると考えてよいことになる (疑問 1 の答え)。

−→
−→
−→
次に A が V の内部にある場合は、|OA| < R なので M̂ (R ∧ |OA|) = M̂ (|OA|) より、−→
mM̂ (|OA|)G
F =−
e3−→
|OA|2となる。すなわちこの場合は、A がいるところより内側の部分の球による重力と同じになるので、A より外側の部分の引力の合力はちょうどつりあって 0 になっていることになる (図 2)。

AR

図 2: A が地球内部の場合
例えば密度が r によらず一定な場合を考えると、M̂ (s) =4π 3s ρ03(ρ(r) = ρ0 )4. 考察

8なので、重力はF =−4πρ0 mG −→|OA|e33−→となり、中心からの距離 |OA| に比例する。つまり、地球の中心から地表までは重力は中心 (原点) からの距離に比例し、地表から外は距離の 2 乗に反比例する、ということになり、よって、地球の重力は地表が一番強く、地球から離れても、地中に潜っても小さくなり、中心では無重力になる。
実際の地球は、中心と地表近くでは構成物質が異なるため密度は一定ではないから、地中での重力は中心からの距離には比例はしないが、似たような状況にはなっていて、地中での重力は地表より強くなるわけではなく、特に中心では重力は 0 になる (疑問 2の答え)。
最後に地球が空洞の場合を考えてみる。この場合、ρ(r) は 0 < r < T および r > R で−→−→0 とすればいい (図 3)。よって、|OA| < T では M̂ (|OA|) = 0 なので (8) より F = 0R
T

図 3: 空洞の場合となり、空洞部分では無重力となる (疑問 3 の答え)。
なお、そのような万有引力の影響がない場では、地球の自転による遠心力が無視できなくなり、その遠心力により自転軸の外側に「重力」を感じて、自転軸を上にして立つことになる。
古い SF で、地底の空洞世界にも小規模の太陽のようなものが地球の中心に浮かんでいて、それをエネルギーとして、「地表人」とは丁度逆向きに立って (地球の中心を天として) 生活する、といった図を見たような気がするが、実際は、赤道では地表人と地底人は逆向きに立つことになるが、極に近づくにつれ地底の重力方向は地面に対して斜めになっていき、遠心力も小さくなるので重力が小さくなっていき、極では 0 になる (図 4)。

5. 空洞世界の力の概算

9

図 4: 自転による遠心力
もし中心に小規模の太陽があるなら、それによる重力も無視できないかもしれないし、もしかすると、外にある月や太陽の引力もそれなりの大きさになってしまうのかもしれない。
その世界では、地上とは物理法則がだいぶ異なるので、それをちゃんと検討すれば面白い SF ネタになるかもしれないが3 、自転の遠心力は、最大の場所 (赤道) でも地球重力に比べて 2,3 桁位小さいので、いずれにせよそのような空洞内は、地上人の我々からすればほぼ無重力状態と言っていいだろう。

空洞世界の力の概算

5

最後に、前節に考察した、地中に空洞がある場合のその中でのいくつかの力の大きさを概算してみる。
まず、地底世界の最大の力である地球の自転の遠心力と、地表での重力との比較を行う。角速度 ω 、半径 r で円運動する質量 m の物体に働く遠心力は mrω 2 なので、rω 2と重力加速度 g = 9.8 [m/s2 ] を比較すればよい。地球の半径 R1 、角速度 ω1 は、地球の円周がほぼ 40000 [km] = 4.0 × 107 [m] なので、

R1 =34.0 × 107
2.0=× 107 [m],2ππω1 =2ππ=× 10−4 [rad/s]
24 × 60 × 60
4.32

すでにそのような考察を行った SF 小説があるかどうかは知らない。

5. 空洞世界の力の概算

10より、2.0 × πR1 ω12=× 10−1 = 0.00343
g4.322 × 9.8(9)なので、地表の最大自転遠心力 F2 は地表重力 F1 の千分の 3 程度ということになる。
地中の空洞内での最大遠心力はもちろん F2 以下であり、次はそれを打ち消すような小太陽の質量 M2 を考えてみる。空洞内の質量 m に働く引力の大きさは、mM2 G/R12 以上であるから、それが最大遠心力 mR1 ω12 を越えるとどの場所でも中心向きの力が勝ってしまうことになる。よって、その限界の質量 M2 では、

mM2 G= mR1 ω122R1となる。一方、地球の重力加速度 g はほぼ地球の万有引力によるものに等しいので、地球の質量 M1 (= 5.97 × 1027 [kg]) と (9) により、

mR1 ω12 = 0.00343 mg = 0.00343 ×mM1 G
mM2 G=
2R1R12から、結局

M2 = 0.00343 M1(10)となる。これよりも大きければ地底世界は全部小太陽に落ちていくことになる。もし
この小太陽の密度が太陽の密度 ρ3 = 1.41 × 103 [kg/m3 ] と同じであれば、その半径をR2 とすると、(10) は地球の密度 ρ1 = 5.52 × 103 [kg/m3 ] を用いて

M2 = ρ3 V2 =4π4π
ρ3 R23 = 0.00343 ×
ρ1 R1333と表され、よって√R2 = R1 3 0.00343 ×

ρ1= 0.238 R1ρ3となり、地球の半径の約 1/4 程度が限界であることがわかる。

5. 空洞世界の力の概算

11

最後に、月や太陽からの引力と地球の重力との比較を行う。まず太陽質量 M3 、および太陽と地球の距離 R3 は

M3 = 3.32 × 105 M1 ,R3 = 1.50 × 1011 [m] = 2.36 × 104 R1であり、よって太陽引力の大きさ F3 と地表面での重力の大きさ F1 との比は、

mM3 G/R32M3F3==2F1mM1 G/R1M1(R1R3)2= 3.32 × 105 × (1/2.36)2 × 10−8 = 5.96 × 10−4となり、地表重力の 1 万分の 6 程度となる。月の質量 M4 、月と地球の距離 R2 は

M4 = 0.0123 M1 ,R2 = 3.84 × 108 [m] = 6.03 × 101 R1

なので、月からの引力 F4 は、

M4F4 =M1(R1R4)2F1 = 0.0123 × (1/6.03)2 × 10−2 F1 = 3.38 × 10−6 F1となり、地表重力の 100 万分の 3 程度、ということになる。地表重力に代わり、最大遠心力 F0 と比較すれば、

F3 = 0.174 F0 ,F4 = 9.85 × 10−4 F0となる。つまり太陽引力 F4 は最大遠心力の 17% なので空洞世界ではかなり影響は強い。月の引力は最大遠心力の 1000 分の 1 程度なので太陽の引力ほどの影響はないが、しかしそれでも地表重力に対する太陽の引力の影響 (5.96 × 10−4 ) よりも桁が一つ大きいので、それなりの影響があることがわかる。
つまり、地底世界は、外の太陽も月も見えないが、それら配置の影響をかなり強く受けてしまうことになる。

参考文献[1] 国立天文台編、「理科年表 平成 2 年」(1990)、丸善.