球冠

青で示された部分が球冠の一例である。

球冠(英語では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)}

となる。

応用

2つの交差する球の和集合と交差部の体積

半径 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}} と簡単にすることができる。

関連項目

ポータル 数学
ポータル 数学
  • 弓形 — 2次元の物体への類推
  • 立体角 — n次元球冠の式に含まれる
  • 球台
  • 半楕円体形
  • 球面扇形(英語版)
  • 球面楔形(英語版)

脚注

  1. ^ 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 .
  2. ^ “Unizor - Geometry3D - Spherical Sectors”. YouTube. Zor Shekhtman. 31 Dec 2018閲覧。
  3. ^ Connolly, Michael L. (1985). “Computation of molecular volume”. Journal of the American Chemical Society 107 (5): 1118–1124. doi:10.1021/ja00291a006. 
  4. ^ 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. 
  5. ^ Bondi, A. (1964). “Van der Waals volumes and radii”. The Journal of Physical Chemistry 68 (3): 441–451. doi:10.1021/j100785a001. 
  6. ^ Scott E. Donaldson, Stanley G. Siegel (2001). Successful Software Development. ISBN 9780130868268. https://books.google.com/books?id=lrix5MNRiu4C&pg=PA354 29 August 2016閲覧。 
  7. ^ “Obliquity of the Ecliptic (Eps Mean)”. Neoprogrammics.com. 2014年5月13日閲覧。
  8. ^ 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. 
  9. ^ 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. 
  10. ^ 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. 
  11. ^ 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. 

外部リンク

ウィキメディア・コモンズには、球冠に関連するカテゴリがあります。
  • Weisstein, Eric W. "Spherical cap". mathworld.wolfram.com (英語). Derivation and some additional formulas.
  • Online calculator for spherical cap volume and area.
  • Summary of spherical formulas.