青で示された部分が球冠の一例である。
球冠 (英語ではspherical cap, spherical domeやspherical segment of one baseという)とは、平面 により切断された球 の一部のこと。平面が球の中心を通り、球冠の高さが球体の半径 と等しいときには半球となる。
球冠の体積 と曲面の面積は、次の値を組み合わせることで計算できる。
球の半径
r
{\displaystyle r}
球冠の底の半径
a
{\displaystyle a}
球冠の高さ
h
{\displaystyle h}
球の中心から球冠の頂点(極)までの線と球冠の底を形作る円板 の端との間の極角
θ
{\displaystyle \theta }
r
{\displaystyle r}
と
h
{\displaystyle h}
を用いる
a
{\displaystyle a}
と
h
{\displaystyle h}
を用いる
r
{\displaystyle r}
と
θ
{\displaystyle \theta }
を用いる
体積
V
=
π
h
2
3
(
3
r
−
h
)
{\displaystyle V={\frac {\pi h^{2}}{3}}(3r-h)}
[1]
V
=
1
6
π
h
(
3
a
2
+
h
2
)
{\displaystyle V={\frac {1}{6}}\pi h(3a^{2}+h^{2})}
V
=
π
3
r
3
(
2
+
cos
θ
)
(
1
−
cos
θ
)
2
{\displaystyle V={\frac {\pi }{3}}r^{3}(2+\cos \theta )(1-\cos \theta )^{2}}
表面積
A
=
2
π
r
h
{\displaystyle A=2\pi rh}
[1]
A
=
π
(
a
2
+
h
2
)
{\displaystyle A=\pi (a^{2}+h^{2})}
A
=
2
π
r
2
(
1
−
cos
θ
)
{\displaystyle A=2\pi r^{2}(1-\cos \theta )}
ϕ
{\displaystyle \phi }
が地理座標 における緯度 を示す場合、
θ
+
ϕ
=
π
/
2
=
90
∘
{\displaystyle \theta +\phi =\pi /2=90^{\circ }\,}
である。
h
{\displaystyle h}
と
r
{\displaystyle r}
の関係は
0
≤
h
≤
2
r
{\displaystyle 0\leq h\leq 2r}
であれば問題ない。例えば、図の赤い部分は
h
>
r
{\displaystyle h>r}
の球冠である。
r
{\displaystyle r}
と
h
{\displaystyle h}
を用いる式は、ピタゴラスの定理 を用いて
r
{\displaystyle r}
の代わりに球冠の底面の半径
a
{\displaystyle a}
を用いる式に書き換えることができる。
r
2
=
(
r
−
h
)
2
+
a
2
=
r
2
+
h
2
−
2
r
h
+
a
2
,
{\displaystyle r^{2}=(r-h)^{2}+a^{2}=r^{2}+h^{2}-2rh+a^{2}\,,}
つまり
r
=
a
2
+
h
2
2
h
.
{\displaystyle r={\frac {a^{2}+h^{2}}{2h}}\,.}
となる。
これを式に代入すると
V
=
π
h
2
3
(
3
a
2
+
3
h
2
2
h
−
h
)
=
1
6
π
h
(
3
a
2
+
h
2
)
,
{\displaystyle V={\frac {\pi h^{2}}{3}}\left({\frac {3a^{2}+3h^{2}}{2h}}-h\right)={\frac {1}{6}}\pi h(3a^{2}+h^{2})\,,}
A
=
2
π
(
a
2
+
h
2
)
2
h
h
=
π
(
a
2
+
h
2
)
.
{\displaystyle A=2\pi {\frac {(a^{2}+h^{2})}{2h}}h=\pi (a^{2}+h^{2})\,.}
となる。
球面扇形の体積から直感的に表面積を導出する[ 編集 ]
以下の議論ベースの計算とは別であるが、球冠の表面積は直感的な議論により球面扇形の体積
V
s
e
c
{\displaystyle V_{sec}}
から導出することができる[2] 。
A
=
3
r
V
s
e
c
=
3
r
2
π
r
2
h
3
=
2
π
r
h
.
{\displaystyle A={\frac {3}{r}}V_{sec}={\frac {3}{r}}{\frac {2\pi r^{2}h}{3}}=2\pi rh\,.}
直感的な議論は、総体積を無限小の三角錐 の総体積を合計することに基づいている。角錐 の体積の式
V
=
1
3
b
h
′
{\displaystyle V={\frac {1}{3}}bh'}
(
b
{\displaystyle b}
は各角錐の底面(球体の表面にある)の無限面積 で、
h
′
{\displaystyle h'}
は底面から頂点(球の中心)までの各角錐の高さ)を用いる。極限をとったときの各
h
′
{\displaystyle h'}
は定数であり、球の半径
r
{\displaystyle r}
と等しいため、無限小の角錐の底面積の合計は球冠の表面積に等しくなる。
V
s
e
c
=
∑
V
=
∑
1
3
b
h
′
=
∑
1
3
b
r
=
r
3
∑
b
=
r
3
A
{\displaystyle V_{sec}=\sum {V}=\sum {\frac {1}{3}}bh'=\sum {\frac {1}{3}}br={\frac {r}{3}}\sum b={\frac {r}{3}}A}
緑の領域を回転させると、高さ
h
{\displaystyle h}
で球の半径
r
{\displaystyle r}
の球冠を作ることができる
体積と表面積の式は次の関数
f
(
x
)
=
r
2
−
(
x
−
r
)
2
=
2
r
x
−
x
2
{\displaystyle f(x)={\sqrt {r^{2}-(x-r)^{2}}}={\sqrt {2rx-x^{2}}}}
(
x
∈
[
0
,
h
]
{\displaystyle x\in [0,h]}
)
を調べ、表面積に対しては回転面 の式、体積に対しては回転体 の式を用いることで導出される。
表面積は
A
=
2
π
∫
0
h
f
(
x
)
1
+
f
′
(
x
)
2
d
x
{\displaystyle A=2\pi \int _{0}^{h}f(x){\sqrt {1+f'(x)^{2}}}\,dx}
である。
f
{\displaystyle f}
の導関数は
f
′
(
x
)
=
r
−
x
2
r
x
−
x
2
{\displaystyle f'(x)={\frac {r-x}{\sqrt {2rx-x^{2}}}}}
であるから
1
+
f
′
(
x
)
2
=
r
2
2
r
x
−
x
2
{\displaystyle 1+f'(x)^{2}={\frac {r^{2}}{2rx-x^{2}}}}
となる。よって表面積の式は
A
=
2
π
∫
0
h
2
r
x
−
x
2
r
2
2
r
x
−
x
2
d
x
=
2
π
∫
0
h
r
d
x
=
2
π
r
[
x
]
0
h
=
2
π
r
h
{\displaystyle A=2\pi \int _{0}^{h}{\sqrt {2rx-x^{2}}}{\sqrt {\frac {r^{2}}{2rx-x^{2}}}}\,dx=2\pi \int _{0}^{h}r\,dx=2\pi r\left[x\right]_{0}^{h}=2\pi rh}
となる。体積は
V
=
π
∫
0
h
f
(
x
)
2
d
x
=
π
∫
0
h
(
2
r
x
−
x
2
)
d
x
=
π
[
r
x
2
−
1
3
x
3
]
0
h
=
π
h
2
3
(
3
r
−
h
)
{\displaystyle V=\pi \int _{0}^{h}f(x)^{2}\,dx=\pi \int _{0}^{h}(2rx-x^{2})\,dx=\pi \left[rx^{2}-{\frac {1}{3}}x^{3}\right]_{0}^{h}={\frac {\pi h^{2}}{3}}(3r-h)}
となる。
半径
r
1
{\displaystyle r_{1}}
と
r
2
{\displaystyle r_{2}}
の球が交差したときの和集合 の体積は[3] 、
V
=
V
(
1
)
−
V
(
2
)
,
{\displaystyle V=V^{(1)}-V^{(2)}\,,}
ここで
V
(
1
)
=
4
π
3
r
1
3
+
4
π
3
r
2
3
{\displaystyle V^{(1)}={\frac {4\pi }{3}}r_{1}^{3}+{\frac {4\pi }{3}}r_{2}^{3}}
は2つの独立した球の体積の合計で
V
(
2
)
=
π
h
1
2
3
(
3
r
1
−
h
1
)
+
π
h
2
2
3
(
3
r
2
−
h
2
)
{\displaystyle V^{(2)}={\frac {\pi h_{1}^{2}}{3}}(3r_{1}-h_{1})+{\frac {\pi h_{2}^{2}}{3}}(3r_{2}-h_{2})}
は交差を作り出す2つの球冠の体積の合計。
d
≤
r
1
+
r
2
{\displaystyle d\leq r_{1}+r_{2}}
が2つの球の中心間の距離であるとき、変数
h
1
{\displaystyle h_{1}}
と
h
2
{\displaystyle h_{2}}
を消すと[4] [5]
V
(
2
)
=
π
12
d
(
r
1
+
r
2
−
d
)
2
(
d
2
+
2
d
(
r
1
+
r
2
)
−
3
(
r
1
−
r
2
)
2
)
.
{\displaystyle V^{(2)}={\frac {\pi }{12d}}(r_{1}+r_{2}-d)^{2}\left(d^{2}+2d(r_{1}+r_{2})-3(r_{1}-r_{2})^{2}\right)\,.}
となる。
2つの平行な平板で囲まれた球台の表面積は、それぞれの球冠の表面積の差である。半径
r
{\displaystyle r}
で高さが
h
1
{\displaystyle h_{1}}
と
h
2
{\displaystyle h_{2}}
の球冠の場合、表面積は
A
=
2
π
r
|
h
1
−
h
2
|
,
{\displaystyle A=2\pi r|h_{1}-h_{2}|\,,}
であり、地理座標である緯度
ϕ
1
{\displaystyle \phi _{1}}
と
ϕ
2
{\displaystyle \phi _{2}}
を用いると[6]
A
=
2
π
r
2
|
sin
ϕ
1
−
sin
ϕ
2
|
,
{\displaystyle A=2\pi r^{2}|\sin \phi _{1}-\sin \phi _{2}|\,,}
となる。例えば、地球を半径6371 kmの球と仮定すると、北極(2016年8月現在、北極圏である緯度66.56°より北[7] )の表面積は、2π ·63712 |sin 90° − sin 66.56°| = 2104万km2 で、地球の総表面積の0.5·|sin 90° − sin 66.56°| = 4.125%である。
この式を用いることで、地球の表面積の半分が南緯30°と北緯30°の間にあることを示すことができる。この範囲は熱帯 を包含する。
回転楕円体のドーム(spheroidal dome )は、ドームが円対称(回転軸を持つ)になるように回転楕円体 の一部を切り取ることで得られる。楕円体 から楕円体のドームも同様に得られる。
一般的に、
n
{\displaystyle n}
次元ユークリッド空間における高さ
h
{\displaystyle h}
で半径
r
{\displaystyle r}
の超球冠の
n
{\displaystyle n}
次元の体積は[8]
V
=
π
n
−
1
2
r
n
Γ
(
n
+
1
2
)
∫
0
arccos
(
r
−
h
r
)
sin
n
(
t
)
d
t
{\displaystyle V={\frac {\pi ^{\frac {n-1}{2}}\,r^{n}}{\,\Gamma \left({\frac {n+1}{2}}\right)}}\int \limits _{0}^{\arccos \left({\frac {r-h}{r}}\right)}\sin ^{n}(t)\,\mathrm {d} t}
で与えられる。ここで
Γ
{\displaystyle \Gamma }
(ガンマ関数 )は
Γ
(
z
)
=
∫
0
∞
t
z
−
1
e
−
t
d
t
{\displaystyle \Gamma (z)=\int _{0}^{\infty }t^{z-1}\mathrm {e} ^{-t}\,\mathrm {d} t}
で与えられる。
V
{\displaystyle V}
の式は、n次元球体 単位の体積
C
n
=
π
n
/
2
/
Γ
[
1
+
n
2
]
{\displaystyle C_{n}={\scriptstyle \pi ^{n/2}/\Gamma [1+{\frac {n}{2}}]}}
や超幾何関数
2
F
1
{\displaystyle {}_{2}F_{1}}
、正規化不完全ベータ関数
I
x
(
a
,
b
)
{\displaystyle I_{x}(a,b)}
を用いて
V
=
C
n
r
n
(
1
2
−
r
−
h
r
Γ
[
1
+
n
2
]
π
Γ
[
n
+
1
2
]
2
F
1
(
1
2
,
1
−
n
2
;
3
2
;
(
r
−
h
r
)
2
)
)
=
1
2
C
n
r
n
I
(
2
r
h
−
h
2
)
/
r
2
(
n
+
1
2
,
1
2
)
{\displaystyle V=C_{n}\,r^{n}\left({\frac {1}{2}}\,-\,{\frac {r-h}{r}}\,{\frac {\Gamma [1+{\frac {n}{2}}]}{{\sqrt {\pi }}\,\Gamma [{\frac {n+1}{2}}]}}{\,\,}_{2}F_{1}\left({\tfrac {1}{2}},{\tfrac {1-n}{2}};{\tfrac {3}{2}};\left({\tfrac {r-h}{r}}\right)^{2}\right)\right)={\frac {1}{2}}C_{n}\,r^{n}I_{(2rh-h^{2})/r^{2}}\left({\frac {n+1}{2}},{\frac {1}{2}}\right)}
と表すことができ、
表面積の式
A
{\displaystyle A}
は、n次元球体単位の表面積
A
n
=
2
π
n
/
2
/
Γ
[
n
2
]
{\displaystyle A_{n}={\scriptstyle 2\pi ^{n/2}/\Gamma [{\frac {n}{2}}]}}
を用いて
A
=
1
2
A
n
r
n
−
1
I
(
2
r
h
−
h
2
)
/
r
2
(
n
−
1
2
,
1
2
)
{\displaystyle A={\frac {1}{2}}A_{n}\,r^{n-1}I_{(2rh-h^{2})/r^{2}}\left({\frac {n-1}{2}},{\frac {1}{2}}\right)}
と表すことができる。ここで
0
≤
h
≤
r
{\displaystyle 0\leq h\leq r}
である。
それより前に[9] (1986, USSR Academ. Press) 次の式が導出されていた。
A
=
A
n
p
n
−
2
(
q
)
,
V
=
C
n
p
n
(
q
)
{\displaystyle A=A_{n}p_{n-2}(q),V=C_{n}p_{n}(q)}
, where
q
=
1
−
h
/
r
(
0
≤
q
≤
1
)
,
p
n
(
q
)
=
(
1
−
G
n
(
q
)
/
G
n
(
1
)
)
/
2
{\displaystyle q=1-h/r(0\leq q\leq 1),p_{n}(q)=(1-G_{n}(q)/G_{n}(1))/2}
,
G
n
(
q
)
=
∫
0
q
(
1
−
t
2
)
(
n
−
1
)
/
2
d
t
{\displaystyle G_{n}(q)=\int \limits _{0}^{q}(1-t^{2})^{(n-1)/2}dt}
.
奇数
n
=
2
k
+
1
:
{\displaystyle n=2k+1:}
に対しては
G
n
(
q
)
=
∑
i
=
0
k
(
−
1
)
i
(
k
i
)
q
2
i
+
1
2
i
+
1
{\displaystyle G_{n}(q)=\sum _{i=0}^{k}(-1)^{i}{\binom {k}{i}}{\frac {q^{2i+1}}{2i+1}}}
.
n
→
∞
{\displaystyle n\to \infty }
および
q
n
=
const.
{\displaystyle q{\sqrt {n}}={\text{const.}}}
である場合、
p
n
(
q
)
→
1
−
F
(
q
n
)
{\displaystyle p_{n}(q)\to 1-F({q{\sqrt {n}}})}
(
F
(
)
{\displaystyle F()}
は標準正規分布 の積分)であることが示されている[10] 。
より定量的な方法でこれを書くと[11] 、境界
A
/
A
n
=
n
Θ
(
1
)
⋅
[
(
2
−
h
/
r
)
h
/
r
]
n
/
2
{\displaystyle A/A_{n}=n^{\Theta (1)}\cdot [(2-h/r)h/r]^{n/2}}
が与えられる。
大きな球冠(
n
→
∞
{\displaystyle n\to \infty }
で
(
1
−
h
/
r
)
4
⋅
n
=
O
(
1
)
{\displaystyle (1-h/r)^{4}\cdot n=O(1)}
)の場合、境界は
n
Θ
(
1
)
⋅
e
−
(
1
−
h
/
r
)
2
n
/
2
{\displaystyle n^{\Theta (1)}\cdot e^{-(1-h/r)^{2}n/2}}
と簡単にすることができる。
^ a b Polyanin, Andrei D; Manzhirov, Alexander V. (2006), Handbook of Mathematics for Engineers and Scientists , CRC Press, p. 69, ISBN 9781584885023 , https://books.google.com/books?id=ge6nk9W0BCcC&pg=PA69 .
^ “Unizor - Geometry3D - Spherical Sectors ”. YouTube . Zor Shekhtman. 31 Dec 2018 閲覧。
^ Connolly, Michael L. (1985). “Computation of molecular volume”. Journal of the American Chemical Society 107 (5): 1118–1124. doi :10.1021/ja00291a006 .
^ Pavani, R.; Ranghino, G. (1982). “A method to compute the volume of a molecule”. Computers & Chemistry 6 (3): 133–135. doi :10.1016/0097-8485(82)80006-5 .
^ Bondi, A. (1964). “Van der Waals volumes and radii”. The Journal of Physical Chemistry 68 (3): 441–451. doi :10.1021/j100785a001 .
^ Scott E. Donaldson, Stanley G. Siegel (2001). Successful Software Development . ISBN 9780130868268 . https://books.google.com/books?id=lrix5MNRiu4C&pg=PA354 29 August 2016 閲覧。
^ “Obliquity of the Ecliptic (Eps Mean) ”. Neoprogrammics.com. 2014年5月13日 閲覧。
^ Li, S (2011). “Concise Formulas for the Area and Volume of a Hyperspherical Cap”. Asian Journal of Mathematics & Statistics 4 (1): 66–70. doi :10.3923/ajms.2011.66.70 .
^ Chudnov, Alexander M. (1986). “On minimax signal generation and reception algorithms (rus.)” . Problems of Information Transmission 22 (4): 49–54. http://mi.mathnet.ru/ppi958 .
^ Chudnov, Alexander M (1991). “Game-theoretical problems of synthesis of signal generation and reception algorithms (rus.)” . Problems of Information Transmission 27 (3): 57–65. http://mi.mathnet.ru/ppi570 .
^ Anja Becker, Léo Ducas, Nicolas Gama, and Thijs Laarhoven. 2016. New directions in nearest neighbor searching with applications to lattice sieving. In Proceedings of the twenty-seventh annual ACM-SIAM symposium on Discrete algorithms (SODA '16), Robert Kraughgamer (Ed.). Society for Industrial and Applied Mathematics, Philadelphia, PA, USA, 10-24.
Richmond, Timothy J. (1984). “Solvent accessible surface area and excluded volume in proteins: Analytical equation for overlapping spheres and implications for the hydrophobic effect”. Journal of Molecular Biology 178 (1): 63–89. doi :10.1016/0022-2836(84)90231-6 . PMID 6548264 .
Lustig, Rolf (1986). “Geometry of four hard fused spheres in an arbitrary spatial configuration”. Molecular Physics 59 (2): 195–207. Bibcode : 1986MolPh..59..195L . doi :10.1080/00268978600102011 .
Gibson, K. D.; Scheraga, Harold A. (1987). “Volume of the intersection of three spheres of unequal size: a simplified formula”. The Journal of Physical Chemistry 91 (15): 4121–4122. doi :10.1021/j100299a035 .
Gibson, K. D.; Scheraga, Harold A. (1987). “Exact calculation of the volume and surface area of fused hard-sphere molecules with unequal atomic radii”. Molecular Physics 62 (5): 1247–1265. Bibcode : 1987MolPh..62.1247G . doi :10.1080/00268978700102951 .
Petitjean, Michel (1994). “On the analytical calculation of van der Waals surfaces and volumes: some numerical aspects”. Journal of Computational Chemistry 15 (5): 507–523. doi :10.1002/jcc.540150504 .
Grant, J. A.; Pickup, B. T. (1995). “A Gaussian description of molecular shape”. The Journal of Physical Chemistry 99 (11): 3503–3510. doi :10.1021/j100011a016 .
Busa, Jan; Dzurina, Jozef; Hayryan, Edik; Hayryan, Shura (2005). “ARVO: A fortran package for computing the solvent accessible surface area and the excluded volume of overlapping spheres via analytic equations”. Computer Physics Communications 165 (1): 59–96. Bibcode : 2005CoPhC.165...59B . doi :10.1016/j.cpc.2004.08.002 .
ウィキメディア・コモンズには、
球冠 に関連するカテゴリがあります。