体積と回転体

同義語:Cavalieriの原理カヴァリエリの原理Pappus–Guldinusの定理Cavalieri's principle

概要

立体の体積は、$x$ 軸に垂直な平面による断面の面積 $S(x)$ の積分 $\int_a^bS(x)\,dx$ に等しい。体積を直方体による内外の近似で定めると、公式は立体が体積をもち各断面の面積 $S(x)$ が定まって積分可能なら成り立ち、断面積が各位置で等しい立体は体積も等しいという Cavalieri の原理が従う。連続な $0\le g\le f$ の間の図形を $x$ 軸の周りに回した体積は $\pi\int_a^b(f^2-g^2)\,dx$(円板法)、$0\le a$ のとき $0\le y\le f(x)$ の図形を $y$ 軸の周りに回した体積は $2\pi\int_a^bxf(x)\,dx$(殻の方法)で、どちらも「面積 × 重心が 1 周で動く距離」に等しい(Pappus–Guldinus の定理)。軸を含む平面による断面積を角度で積分したものは体積とは限らない。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 定積分, 区分求積と積分の定義, Riemann積分, 回転体, 重心

高校での出発点

高校では、立体の体積を断面積の積分で求める。立体が平面 $x=a$ と $x=b$ の間にあり、$x$ 軸に垂直な平面で切った断面の面積が $S(x)$ なら
$$ V=\int_a^bS(x)\,dx $$
である。特に、曲線 $y=f(x)$($f\ge0$)と $x$ 軸の間の部分を $x$ 軸の周りに回転させた立体では、断面が半径 $f(x)$ の円板なので $V=\pi\int_a^bf(x)^2\,dx$ となる。根拠としては、「$a$ から $x$ までの部分の体積を $V(x)$ とすると、厚さ $\Delta x$ の薄い板の体積はほぼ $S(x)\Delta x$ なので $V'(x)=S(x)$」という説明がされる。

断面が正三角形の立体

底面が円板 $x^2+y^2\le r^2$ で、$x$ 軸に垂直な平面で切った断面が、底面上の弦を 1 辺とする正三角形になる立体を考える。位置 $x$ の弦の長さは $2\sqrt{r^2-x^2}$ なので、断面積は $\frac{\sqrt3}4\bigl(2\sqrt{r^2-x^2}\bigr)^2=\sqrt3(r^2-x^2)$ であり、
$$ V=\int_{-r}^r\sqrt3(r^2-x^2)\,dx=\frac{4\sqrt3}3r^3 $$
である。同じ計算で、半径 $r$ の球(断面は半径 $\sqrt{r^2-x^2}$ の円板)の体積は $\int_{-r}^r\pi(r^2-x^2)\,dx=\frac43\pi r^3$ である。

この方法は強力だが、次の点が残っている。

  1. 体積とは何か。「$V'(x)=S(x)$」の議論は、体積 $V(x)$ がすでに定まっていることを前提にしている。(→ def-vcav-volume)
  2. 「薄い板の体積はほぼ $S(x)\Delta x$」はいつ正しいか。断面が $x$ とともに横へずれていく立体では、板は柱ではない。どういう切り方なら断面積を積分してよく、どういう切り方では誤るのか。(→ thm-vcav-slice、ex-vcav-fan)
  3. 回転体の公式は 1 つではない。$y$ 軸の周りの回転では、円板で切る方法と、同心の円筒(殻)に分ける方法がある。両者は一致するのか。回転体の体積と図形の重心の関係(Pappus–Guldinus の定理)はなぜ成り立つのか。(→ prop-vcav-agree、thm-vcav-pappus)
    本記事では、体積を直方体による内側と外側からの近似で定義し(積分の上限和・下限和と同じ考え)、断面積の公式と Cavalieri の原理(断面積が等しい 2 つの立体は体積が等しい)を証明する。そのうえで回転体の体積の 2 つの方法、Pappus–Guldinus の定理を証明し、断面の取り方を誤る反例を調べる。1 変数の積分については、連続関数や単調関数が積分可能であること、連続関数では分割の幅を $0$ に近づけると代表点によらず Riemann 和が積分に収束すること、微分積分学の基本定理を前提とする(証明は 区分求積と積分の定義 の記事や Leb26 Lemma 5.2.7・Proposition 5.2.11・Exercise 5.1.13・Theorem 5.3.1・5.3.3 にある)。
    高校の計算大学の概念ボックス
    立体の体積直方体による内側と外側の近似が一致したときの値def-vcav-volume
    $V=\int_a^bS(x)\,dx$断面積の積分(Fubini の定理の特別な場合)thm-vcav-slice
    断面積が等しい立体を比べるCavalieri の原理cor-vcav-cavalieri
    $\pi\int_a^bf(x)^2\,dx$断面が円環の立体への thm-vcav-slice の適用thm-vcav-disk
    $2\pi\int_a^bxf(x)\,dx$殻による挟み撃ち(円板法を部分積分したもの)thm-vcav-shell、prop-vcav-agree
    面積 × 重心の動く距離Pappus–Guldinus の定理thm-vcav-pappus

体積の定義

座標空間の直方体(各辺が座標軸に平行なもの)$[a_1,b_1]\times[a_2,b_2]\times[a_3,b_3]$ の体積は $(b_1-a_1)(b_2-a_2)(b_3-a_3)$ とする。有界な立体 $K$ を含む直方体 $R$ をとり、$R$ の各辺の区間を分割すると、$R$ は小さな直方体(セル)に分かれる。この分け方を $R$ の分割という。

体積と面積

有界な集合 $K\subset\mathbb{R}^3$ と、$K$ を含む直方体 $R$ の分割 $P$ について、
$$ v_*(K;P):=(\text{$K$ に含まれるセルの体積の和}),\qquad v^*(K;P):=(\text{$K$ と共通点をもつセルの体積の和}) $$
とおく。$\sup_Pv_*(K;P)=\inf_Pv^*(K;P)$($P$ は $R$ のすべての分割を動く)のとき、$K$ は体積をもつといい、この共通の値を $K$ の体積 $\mathrm{vol}(K)$ とする。平面の有界な集合 $E\subset\mathbb{R}^2$ についても、長方形のセルで同じように $a_*(E;Q)$、$a^*(E;Q)$ を定め、両者の上限と下限が一致するとき $E$ は面積をもつといい、共通の値を面積 $\mathrm{area}(E)$ とする。

$v_*(K;P)$ と $v^*(K;P)$ は、$K$ 上で $1$、それ以外で $0$ をとる関数 $\chi_K$ の、分割 $P$ についての下限和と上限和である。セルが $K$ に含まれるときに限り $\chi_K$ のセル上の下限が $1$ になり、セルが $K$ と交わるときに限り上限が $1$ になるからである。したがって $K$ が体積をもつことは $\chi_K$ が $R$ 上で積分可能であることと同じで、$\mathrm{vol}(K)=\int_R\chi_K$ である。これは Leb26b §10.5 の Jordan 可測集合とその体積の定義と一致し、$R$ の選び方によらない。定義からすぐ分かるように、体積は単調である($K\subset L$ なら、どの $P$ でも $v_*(K;P)\le v_*(L;P)$ なので $\mathrm{vol}(K)\le\mathrm{vol}(L)$)。また、平行移動しても、$\lambda>0$ 倍に拡大しても(分割も同じように移せばよいので)、体積はそれぞれ変わらない・$\lambda^3$ 倍になり、面積は変わらない・$\lambda^2$ 倍になる。

直角二等辺三角形の面積を正方形で挟む

$E=\{(x,y)\mid0\le y\le x\le1\}$ を含む正方形 $[0,1]^2$ を、1 辺 $\frac1n$ の $n^2$ 個の正方形のセル $\bigl[\frac{i-1}n,\frac in\bigr]\times\bigl[\frac{j-1}n,\frac jn\bigr]$ に分ける。セルが $E$ に含まれるのは $j\le i-1$ のとき、$E$ と共通点をもつのは $j\le i+1$ のときなので、
$$ a_*=\frac1{n^2}\cdot\frac{n(n-1)}2=\frac{n-1}{2n},\qquad a^*=\frac1{n^2}\cdot\frac{n^2+3n-2}2 $$
である($n=10$ で $0.45$ と $0.64$)。$n\to\infty$ で両者は $\frac12$ に近づくので、$E$ は面積をもち、$\mathrm{area}(E)=\frac12$ である。

次の 2 つは証明しない。

  • 加法性:体積をもつ $K_1,\dots,K_n$ の、どの 2 つの共通部分も体積 $0$ ならば、$\mathrm{vol}(K_1\cup\cdots\cup K_n)=\sum_j\mathrm{vol}(K_j)$ である(指示関数の等式 $\chi_{K\cup L}=\chi_K+\chi_L-\chi_{K\cap L}$ と積分の線形性 Leb26b Proposition 10.1.10 から従う)。平面や円柱面の有界な部分、有界閉集合の上で連続な関数のグラフ(例えば閉区間 $[a,b]$ 上の連続関数 $g$ による、$a\le\sqrt{x^2+z^2}\le b$ 上の曲面 $y=g(\sqrt{x^2+z^2})$)は体積 $0$ である。
  • 存在:本記事で扱う立体・図形(連続関数のグラフで囲まれた図形、球、円柱、錐、回転体など)は体積・面積をもつ。一般に、境界が測度 $0$ の有界集合は体積をもつ(Leb26b Proposition 10.5.1)。また、半径 $r$ の円板の面積は $\pi r^2$、内径 $p$・外径 $q$ の円環の面積は $\pi(q^2-p^2)$ である。

断面積の積分と Cavalieri の原理

立体 $K$ と実数 $x$ について、$K_x:=\{(y,z)\mid(x,y,z)\in K\}$ を $K$ の位置 $x$ での断面という(平面 $x=\text{一定}$ で切った切り口を $yz$ 平面に写したもの)。

断面積の積分による体積

有界な立体 $K$ が体積をもち、$K$ は $a\le x\le b$ の範囲にあるとする。各 $x\in[a,b]$ で断面 $K_x$ が面積 $S(x)$ をもち、$S$ が $[a,b]$ で積分可能ならば、
$$ \mathrm{vol}(K)=\int_a^bS(x)\,dx $$
である。

セルを柱ごとにまとめる

$K$ を含む直方体を $R=[a,b]\times Q$($Q$ は $yz$ 平面の長方形)とする。$R$ の分割 $P$ は、$[a,b]$ の分割 $P_1=\{x_0,\dots,x_n\}$ と、$Q$ を長方形 $C$ に分ける分割 $Q'$ の組で、セルは $[x_{k-1},x_k]\times C$ の形である。
内側:$s_k:=\sum\{\mathrm{area}(C)\mid[x_{k-1},x_k]\times C\subset K\}$ とおくと、$v_*(K;P)=\sum_ks_k\,\Delta x_k$ である。$x\in[x_{k-1},x_k]$ なら、ここに現れる $C$ はどれも $K_x$ に含まれるので、$s_k\le a_*(K_x;Q')\le\mathrm{area}(K_x)=S(x)$ である(面積は内側の和の上限だから)。よって $s_k$ は $[x_{k-1},x_k]$ での $S$ の下限 $m_k$ 以下で、
$$ v_*(K;P)\le\sum_km_k\,\Delta x_k=L(S;P_1)\le\int_a^bS(x)\,dx $$
である($L(S;P_1)$ は $S$ の下限和)。
外側:$t_k:=\sum\{\mathrm{area}(C)\mid([x_{k-1},x_k]\times C)\cap K\ne\emptyset\}$ とおくと $v^*(K;P)=\sum_kt_k\,\Delta x_k$ である。$x\in[x_{k-1},x_k]$ で $C$ が $K_x$ と交われば、セル $[x_{k-1},x_k]\times C$ は $K$ と交わる。よって $S(x)\le a^*(K_x;Q')\le t_k$ であり、$S$ の上限和について $U(S;P_1)\le v^*(K;P)$、したがって $\int_a^bS\le v^*(K;P)$ である。
$P$ は任意なので、$\mathrm{vol}(K)=\sup_Pv_*(K;P)\le\int_a^bS\le\inf_Pv^*(K;P)=\mathrm{vol}(K)$ である。$\square$

証明は、1 変数の積分の上限和・下限和の考えを、セルを $x$ 方向の柱ごとにまとめて使っただけである。$x$ の代わりに $y$ や $z$ を使っても同じである。高校の「$V'(x)=S(x)$」の議論は、この定理の特別な場合を微分の言葉で述べたものといえる。

Cavalieriの原理

体積をもつ 2 つの立体 $K,L$ が $a\le x\le b$ の範囲にあり、すべての $x\in[a,b]$ で断面 $K_x$、$L_x$ が面積をもって $\mathrm{area}(K_x)=\mathrm{area}(L_x)$ であり、この面積が $x$ について積分可能ならば、$\mathrm{vol}(K)=\mathrm{vol}(L)$ である。

同じ積分で表される

thm-vcav-slice により、どちらの体積も同じ関数 $S(x)=\mathrm{area}(K_x)=\mathrm{area}(L_x)$ の積分に等しい。$\square$

断面の形は違ってよく、面積が各位置で等しければよい。Cavalieri の原理は、体積を計算せずに「形の分かっている立体」と比べるために使う。

斜めの柱と錐
  1. 面積 $B$ の平面図形を底面とし、底面を平行移動しながら高さ $h$ まで積み上げた柱(斜めの柱でもよい)は、どの高さの断面も底面の平行移動なので面積 $B$ であり、体積は $\int_0^hB\,dz=Bh$ である。直円柱と同じ体積になる。
  2. 面積 $B$ の底面と、底面から高さ $h$ の頂点をもつ錐(円錐・角錐など)を考える。頂点から距離 $t$ の平面での断面は、頂点を中心に底面を $\frac th$ 倍に縮小した図形なので、面積は $B\bigl(\frac th\bigr)^2$ である。よって
    $$ V=\int_0^hB\frac{t^2}{h^2}\,dt=\frac13Bh $$
    であり、底面の形によらず「底面積 × 高さ ÷ 3」となる。
半球と、円柱から円錐を除いた立体

半径 $r$ の半球 $\{x^2+y^2+z^2\le r^2,\ z\ge0\}$ と、半径 $r$・高さ $r$ の円柱から、上面を底面とし頂点を原点とする円錐をくり抜いた立体 $L$ を比べる。高さ $z$($0\le z\le r$)での半球の断面は半径 $\sqrt{r^2-z^2}$ の円板で、面積は $\pi(r^2-z^2)$ である。$L$ の断面は半径 $r$ の円板から半径 $z$ の円板を除いた円環で、面積は $\pi r^2-\pi z^2$ である。両者は等しいので、cor-vcav-cavalieri により
$$ \text{半球の体積}=\pi r^2\cdot r-\frac13\pi r^2\cdot r=\frac23\pi r^3 $$
であり、球の断面積を積分しなくても、円柱と円錐の体積(ex-vcav-prism-cone)から球の体積 $\frac43\pi r^3$ が得られる。

2 つの円柱の共通部分:切る向きを選ぶ

半径 $r$ の 2 つの円柱 $y^2+z^2\le r^2$($x$ 軸が中心軸)と $x^2+z^2\le r^2$($y$ 軸が中心軸)の共通部分 $K$ の体積を求める。
$x$ 軸に垂直な平面で切ると、断面は円板 $y^2+z^2\le r^2$ を帯 $|z|\le\sqrt{r^2-x^2}$ で切り取った図形で、面積は弓形の面積の公式を使う複雑な式になる。これに対し、$z$ 軸に垂直な平面で切ると、断面は $|x|\le\sqrt{r^2-z^2}$、$|y|\le\sqrt{r^2-z^2}$ の正方形で、面積は $4(r^2-z^2)$ である。よって
$$ \mathrm{vol}(K)=\int_{-r}^r4(r^2-z^2)\,dz=\frac{16}3r^3 $$
である。2 つの中心軸の両方に平行な平面で切ると、どちらの円柱も幅 $2\sqrt{r^2-z^2}$ の帯として現れ、断面が単純になる。

回転体の体積:円板法と殻の方法

以下、$f,g$ は $[a,b]$ で連続とする。

円板法

$0\le g\le f$ とする。図形 $D=\{(x,y)\mid a\le x\le b,\ g(x)\le y\le f(x)\}$ を $x$ 軸の周りに 1 回転させた立体の体積は
$$ V=\pi\int_a^b\bigl(f(x)^2-g(x)^2\bigr)\,dx $$
である。特に $g=0$ なら $V=\pi\int_a^bf(x)^2\,dx$ である。

断面は円環

位置 $x$ での断面は、軸からの距離が $g(x)$ 以上 $f(x)$ 以下の点の集合、すなわち内径 $g(x)$・外径 $f(x)$ の円環($g(x)=0$ なら円板)で、面積は $\pi(f(x)^2-g(x)^2)$ である。これは $x$ の連続関数なので積分可能であり、thm-vcav-slice から主張が従う。$\square$

回転放物体

$y=\sqrt x$($0\le x\le4$)と $x$ 軸の間の図形を $x$ 軸の周りに回転させると、thm-vcav-disk により
$$ V=\pi\int_0^4x\,dx=8\pi=25.13\ldots $$
である。これは外接する半径 $2$・長さ $4$ の円柱の体積 $16\pi$ のちょうど半分である。

$y$ 軸の周りの回転では、軸に垂直な平面($y=\text{一定}$)で切ると断面は円環の組合せになり、$y$ の関数として書くには $f$ の逆関数が要る。代わりに、回転体を同心の薄い円筒(殻)に分ける。

殻の方法

$0\le a< b$、$f\ge0$ とする。図形 $D=\{(x,y)\mid a\le x\le b,\ 0\le y\le f(x)\}$ を $y$ 軸の周りに 1 回転させた立体 $K$ の体積は
$$ V=2\pi\int_a^bx\,f(x)\,dx $$
である。

殻で内側と外側から挟む

点 $(x,y,z)$ の $y$ 軸からの距離を $\rho:=\sqrt{x^2+z^2}$ と書くと、$K=\{a\le\rho\le b,\ 0\le y\le f(\rho)\}$ である。$0\le p< q$、$c\ge0$ について、殻 $T(p,q,c):=\{p\le\rho\le q,\ 0\le y\le c\}$ の体積は $\pi(q^2-p^2)c$ である($c>0$ なら、殻は $0\le y\le c$ の範囲にあり、$y$ 軸に垂直な断面はどれも面積 $\pi(q^2-p^2)$ の円環($p=0$ なら円板)なので、thm-vcav-slice による。$c=0$ なら平面の一部で体積 $0$ である)。
$[a,b]$ の分割 $\{x_0,\dots,x_n\}$ をとり、$[x_{k-1},x_k]$ での $f$ の最小値・最大値を $m_k,M_k$ とする(連続関数は閉区間で最小値・最大値をとる。Leb26 Theorem 3.3.2)。$\rho\in[x_{k-1},x_k]$、$0\le y\le m_k$ なら $y\le f(\rho)$ なので $T(x_{k-1},x_k,m_k)\subset K$ であり、逆に $K\subset\bigcup_kT(x_{k-1},x_k,M_k)$ である。異なる $k$ の殻どうしは円柱面 $\rho=x_k$ でしか交わらない(共通部分の体積は $0$)ので、単調性と加法性から
$$ \sum_k\pi(x_k^2-x_{k-1}^2)\,m_k\le V\le\sum_k\pi(x_k^2-x_{k-1}^2)\,M_k $$
である。$c_k:=\frac{x_{k-1}+x_k}2$ とおくと $\pi(x_k^2-x_{k-1}^2)=2\pi c_k\,\Delta x_k$ である。
$\varphi(x):=2\pi xf(x)$ は連続である。$f(\xi_k)=m_k$ となる $\xi_k\in[x_{k-1},x_k]$ をとると、$|c_k-\xi_k|\le\Delta x_k$ なので、$M:=\max f$ として
$$ \Bigl|\sum_k2\pi c_km_k\,\Delta x_k-\sum_k\varphi(\xi_k)\,\Delta x_k\Bigr|\le2\pi\sum_k|c_k-\xi_k|\,m_k\,\Delta x_k\le2\pi M(b-a)\max_k\Delta x_k $$
である。分割の幅を $0$ に近づけると、右辺は $0$ に近づき、Riemann 和 $\sum_k\varphi(\xi_k)\Delta x_k$ は $\int_a^b\varphi$ に近づく。よって左の和は $\int_a^b\varphi$ に近づく。$M_k$ の側も同様である。$V$ は両者の間にあるので $V=\int_a^b2\pi xf(x)\,dx$ である。$\square$

殻の方法では、$f$ の逆関数が要らない。

正弦曲線の山を y 軸の周りに回す

$D=\{0\le x\le\pi,\ 0\le y\le\sin x\}$ を $y$ 軸の周りに回転させる。殻の方法では、部分積分により
$$ V=2\pi\int_0^\pi x\sin x\,dx=2\pi\bigl[-x\cos x+\sin x\bigr]_0^\pi=2\pi^2 $$
である。$y$ 軸に垂直に切る方法では、高さ $y$($0\le y\le1$)での断面が内径 $\arcsin y$・外径 $\pi-\arcsin y$ の円環なので、
$$ V=\pi\int_0^1\bigl((\pi-\arcsin y)^2-(\arcsin y)^2\bigr)\,dy=\pi\int_0^1\bigl(\pi^2-2\pi\arcsin y\bigr)\,dy=\pi\Bigl(\pi^2-2\pi\bigl(\tfrac\pi2-1\bigr)\Bigr)=2\pi^2 $$
となる($\int_0^1\arcsin y\,dy=\frac\pi2-1$)。値は一致するが、逆関数を 2 本使う分だけ手間がかかる。

2 つの方法が同じ立体の体積を表すことは、どちらも体積そのものに等しいことから分かる。体積を経由せずに、積分の計算だけで一致を確かめることもできる。

円板法と殻の方法の一致

$R>0$ とし、$f$ は $[0,R]$ で微分可能で $f'$ が連続、狭義減少で $f(R)=0$ とする。$H:=f(0)$、$g\colon[0,H]\to[0,R]$ を $f$ の逆関数とすると
$$ \pi\int_0^Hg(y)^2\,dy=2\pi\int_0^Rx\,f(x)\,dx $$
である。左辺は $D=\{0\le x\le R,\ 0\le y\le f(x)\}$ を $y$ 軸の周りに回転させた立体を $y$ 軸に垂直に切ったとき(断面は半径 $g(y)$ の円板)の体積、右辺は殻の方法による体積である。

置換積分と部分積分

左辺で $y=f(x)$ と置換する($x$ が $R$ から $0$ まで動くとき $y$ は $0$ から $H$ まで動き、$g(f(x))=x$)。置換積分の公式(Leb26 Theorem 5.3.5)により
$$ \int_0^Hg(y)^2\,dy=\int_R^0x^2f'(x)\,dx=-\int_0^Rx^2f'(x)\,dx $$
である。部分積分すると
$$ -\int_0^Rx^2f'(x)\,dx=-\bigl[x^2f(x)\bigr]_0^R+\int_0^R2xf(x)\,dx=\int_0^R2xf(x)\,dx $$
である($R^2f(R)=0$、$0^2f(0)=0$)。両辺に $\pi$ を掛ければよい。$\square$

円錐と回転放物面で確かめる

$f(x)=h\bigl(1-\frac xr\bigr)$($0\le x\le r$)なら、立体は底面の半径 $r$・高さ $h$ の円錐で、どちらの方法でも $\frac13\pi r^2h$ になる。$f(x)=H\bigl(1-\frac{x^2}{R^2}\bigr)$ なら回転放物面の内部で、$g(y)=R\sqrt{1-y/H}$ から $\pi\int_0^HR^2(1-\frac yH)\,dy=\frac12\pi R^2H$、殻の方法でも $2\pi\int_0^RxH(1-\frac{x^2}{R^2})\,dx=\frac12\pi R^2H$ である。

部分積分の境界項 $\bigl[x^2f(x)\bigr]_0^R$ が消えることが一致の鍵で、殻の方法は円板法を「部分積分したもの」と見ることができる。

Pappus–Guldinus の定理

回転体の体積は、回転させる図形の面積と重心の位置だけで決まる。まず、2 つのグラフの間の図形の重心を定める。

縦線図形の面積と重心

$g\le f$ を $[a,b]$ で連続な関数とし、$D:=\{(x,y)\mid a\le x\le b,\ g(x)\le y\le f(x)\}$ の面積 $A:=\int_a^b\bigl(f(x)-g(x)\bigr)dx$ が正であるとする(thm-vcav-slice と同じ証明で、この $A$ は def-vcav-volume の意味の $D$ の面積に等しい)。
$$ \bar x:=\frac1A\int_a^bx\bigl(f(x)-g(x)\bigr)dx,\qquad\bar y:=\frac1A\int_a^b\frac{f(x)^2-g(x)^2}2\,dx $$
とおき、点 $(\bar x,\bar y)$ を $D$ の重心という。

この定義は、$D$ を幅 $\Delta x$ の縦の細長い帯に分けたとき、位置 $x$ の帯がほぼ面積 $(f-g)\Delta x$ の長方形で、その重心が $\bigl(x,\frac{f+g}2\bigr)$ にあることから来ている。$\bar y$ の被積分関数 $\frac{f+g}2(f-g)$ は、帯の重心の高さと帯の面積の積である。厚さが一様な板の釣り合いの点が $(\bar x,\bar y)$ である。

Pappus–Guldinusの定理

def-vcav-centroid の図形 $D$ について、$g\ge0$ ならば、$D$ を $x$ 軸の周りに 1 回転させた立体の体積は $2\pi\bar y\cdot A$ である。また $a\ge0$ ならば、$D$ を $y$ 軸の周りに 1 回転させた立体の体積は $2\pi\bar x\cdot A$ である。すなわち、体積は「図形の面積 × 重心が 1 回転で動く距離」に等しい。

円板法と殻の方法に帰着する

$x$ 軸の周り:thm-vcav-disk により体積は $\pi\int_a^b(f^2-g^2)\,dx=2\pi\int_a^b\frac{f^2-g^2}2\,dx=2\pi\bar yA$ である。
$y$ 軸の周り:$D$ を $y$ 軸方向に平行移動すると、回転体も $y$ 軸方向に平行移動するだけなので体積は変わらず、$f-g$ と $\bar x$ も変わらない。そこで $f,g$ に同じ定数を足して、$g\ge0$ としてよい。$\{0\le y\le f(x)\}$ を回転させた立体は、$D$ を回転させた立体と $\{0\le y\le g(x)\}$ を回転させた立体の和で、両者は曲面 $y=g(\rho)$ の上でしか交わらない(共通部分の体積は $0$)。加法性と thm-vcav-shell により、$D$ を回転させた立体の体積は
$$ 2\pi\int_a^bxf(x)\,dx-2\pi\int_a^bxg(x)\,dx=2\pi\int_a^bx\bigl(f(x)-g(x)\bigr)dx=2\pi\bar xA $$
である。$\square$

対称な関数の積分で重み $x$ を区間の中点に置き換える技法は 対称性を用いた積分 で扱う。
本記事で証明したのは縦線図形(2 つのグラフの間の図形)の場合である。有限個の縦線図形に分けられる図形でも、体積と「面積 × 重心の座標」がどちらも分割について足し算になるので、同じ結論が成り立つ。一般の図形では重心を 2 重積分 $\bar y=\frac1A\iint_Dy\,dx\,dy$ で定め、同じ定理が成り立つが、ここでは証明しない。

トーラスの体積

$R\ge r>0$ とし、中心 $(0,R)$・半径 $r$ の円板 $D$ を $x$ 軸の周りに回転させた立体(トーラス、ドーナツ形)の体積を求める。$D$ は $f(x)=R+\sqrt{r^2-x^2}$、$g(x)=R-\sqrt{r^2-x^2}$($-r\le x\le r$)の間の縦線図形で、$g\ge0$ である。$A=\pi r^2$ であり、
$$ \frac{f^2-g^2}2=\frac{(f+g)(f-g)}2=R\cdot2\sqrt{r^2-x^2} $$
なので $A\bar y=R\int_{-r}^r2\sqrt{r^2-x^2}\,dx=R\cdot\pi r^2$、すなわち $\bar y=R$(円板の中心)である。thm-vcav-pappus により
$$ V=2\pi R\cdot\pi r^2=2\pi^2Rr^2 $$
である。$R=2$、$r=1$ なら $4\pi^2=39.47\ldots$ である。「円板を長さ $2\pi R$ の円柱に伸ばした体積」と同じになる。

逆向きに使う:半円板の重心

半径 $r$ の半円板 $\{x^2+y^2\le r^2,\ y\ge0\}$ を $x$ 軸の周りに回転させると半径 $r$ の球になる。thm-vcav-pappus から $\frac43\pi r^3=2\pi\bar y\cdot\frac12\pi r^2$ なので、$\bar y=\frac{4r}{3\pi}=0.424\ldots\times r$ である。同様に、頂点 $(0,0)$、$(h,0)$、$(h,r)$ の直角三角形($0\le y\le\frac rhx$)では $A=\frac12rh$、$\bar y=\frac r3$ で、$x$ 軸の周りの回転体(円錐)の体積 $2\pi\cdot\frac r3\cdot\frac12rh=\frac13\pi r^2h$ が再現される。

反例:断面の取り方を誤る

thm-vcav-slice は、平行な平面で切り、断面積を平面の間の距離で積分することを要求する。また thm-vcav-disk は、図形が回転軸の片側にあること($0\le g$)を要求する。これらを外すと誤る。

外した仮定崩れる主張ボックス
図形が回転軸の片側にある($0\le g$)円板法 $\pi\int(f^2-g^2)$ が体積になるex-vcav-straddle
平行な平面で切り、平面の間の距離で積分する断面積をパラメータで積分すれば体積になるex-vcav-fan
近似の誤差が $0$ に近づくこと(体積では成り立つ)円柱で近似した側面積 $2\pi\int f$ が正しいrem-vcav-surface
各断面が面積をもつ(体積をもつだけでは足りない)thm-vcav-slice の $S(x)$ がすべての $x$ で定まるex-vcav-no-section
反例:回転軸をまたぐ図形

$D=\{0\le x\le1,\ x-1\le y\le x\}$ を $x$ 軸の周りに回転させる。$D$ は $x$ 軸をまたいでいる($g(x)=x-1<0$)。円環の公式をそのまま使うと
$$ \pi\int_0^1\bigl(x^2-(x-1)^2\bigr)dx=\pi\int_0^1(2x-1)\,dx=0 $$
となり、明らかに誤りである。位置 $x$ の断面は、$y$ が $x-1$ から $x$ までの線分($0$ を含む)を回転させたものなので、半径 $\max(x,1-x)$ の円板である。正しい体積は
$$ \pi\int_0^1\max(x,1-x)^2\,dx=2\pi\int_{1/2}^1x^2\,dx=\frac7{12}\pi $$
である。$\pi\int(f^2+g^2)$ としても $\frac23\pi$ で誤る。軸の下側の部分は折り返すと上側の部分と重なり、重なった部分は二重に数えてはいけない。この反例は、$0\le g$ の仮定を外すと thm-vcav-disk の等式が成り立たないことを示す。

反例:軸を含む平面で扇形に切る

半径 $r$・高さ $h$ の直円柱を、中心軸を含む半平面で切ると、断面はどの向きでも $r\times h$ の長方形で、面積は $rh$ である。向きの角 $\theta$($0\le\theta<2\pi$)について断面積を積分すると $\int_0^{2\pi}rh\,d\theta=2\pi rh$ となるが、体積は $\pi r^2h$ であり、$r\ne2$ なら値が違う。
誤りの原因は、隣り合う 2 つの半平面の間の「板」の厚さが一定でないことである。軸から距離 $\rho$ の所での厚さは $\rho\,\Delta\theta$ で、軸から離れるほど厚い。正しくは、断面の各点に厚さの比 $\rho$ を掛けて $\int_0^{2\pi}\Bigl(\int_0^r\rho h\,d\rho\Bigr)d\theta=\pi r^2h$ となる(円柱座標への変数変換。本記事では証明しない)。内側の積分 $\int_0^r\rho h\,d\rho=\frac{r^2h}2$ は「断面の面積 $rh$ × 断面の重心の軸からの距離 $\frac r2$」であり、これは Pappus–Guldinus の定理と同じ形である。この反例は「どんな切り方でも、断面積をパラメータで積分すれば体積になる」という主張を破る。平行な平面で切り、平面の間の距離で積分するという条件は外せない。

同じ近似で表面積を求めると誤る

回転体の体積では、厚さ $\Delta x$ の部分を円柱で近似してよかった。同じ近似で回転体の側面積を求め、$2\pi\int_a^bf(x)\,dx$ としてはいけない。半径 $r$ の球では $2\pi\int_{-r}^r\sqrt{r^2-x^2}\,dx=\pi^2r^2=9.87\ldots\times r^2$ となり、正しい表面積 $4\pi r^2=12.57\ldots\times r^2$ と合わない。体積では、円柱で近似したときの誤差は各部分で $\pi(M_k^2-m_k^2)\Delta x_k$ 以下($m_k,M_k$ はその部分での $f$ の最小値・最大値)で、合計は $\pi f^2$ の上限和と下限和の差なので $0$ に近づく。側面積では、斜めの帯を縦の帯で置き換えると、傾き $f'$ の分だけ面積の比が $\sqrt{1+f'(x)^2}$ ずれ、分割を細かくしてもこのずれは消えない。正しい公式は $2\pi\int_a^bf(x)\sqrt{1+f'(x)^2}\,dx$ である(本記事では証明しない)。

Fubini の定理(紹介)

この節は紹介であり、証明しない。thm-vcav-slice は、多変数の積分の順序に関する Fubini の定理の特別な場合である。直方体 $R=[a,b]\times Q$ 上の関数 $F(x,y,z)$ の積分 $\iiint_RF$ も、上限和と下限和で 1 変数の場合と同じように定義される。Fubini の定理は、$F$ が $R$ 上で積分可能ならば
$$ \iiint_RF=\int_a^b\Bigl(\iint_QF(x,y,z)\,dy\,dz\Bigr)dx $$
が成り立つ、という定理である(Leb26b Theorem 10.2.2。内側の積分が存在しない $x$ があってもよいように、内側は下積分または上積分で述べられている)。$F=\chi_K$ とすると、左辺は $\mathrm{vol}(K)$、内側の積分は断面積 $S(x)$ であり、thm-vcav-slice になる。重心の 2 重積分による定義と def-vcav-centroid が一致することも、この定理から従う。
内側の積分が存在しない例がある。

反例:断面が面積をもたない立体

立方体 $[0,1]^3$ に、平面 $x=2$ 上の点の集合 $\{2\}\times\bigl((\mathbb{Q}\cap[0,1])\times(\mathbb{Q}\cap[0,1])\bigr)$ を付け加えた立体 $K$ を考える。付け加えた集合は正方形 $\{2\}\times[0,1]^2$ に含まれるので体積 $0$ で、$K$ は体積 $1$ をもつ。しかし位置 $x=2$ の断面は、座標が両方とも有理数の点の集合で、内側の和はつねに $0$、外側の和はつねに $1$ 以上なので面積をもたない。この反例は「体積をもつ立体の断面はどれも面積をもつ」という主張を破る。このため thm-vcav-slice では「各断面が面積をもつ」を仮定に入れた。

断面積が定まらない $x$ が有限個しかないなら、そこでの値を変えても積分は変わらないので、実際の計算で困ることは少ない。
Lebesgue 積分の理論では、Fubini の定理はより広い関数について成り立ち、「断面積が定まらない $x$ の集合は測度 $0$」という形で上の例を扱う。

判断の目安

立体の形切り方・方法注意
断面が相似な図形(錐)頂点からの距離 $t$ で切る。断面積は $t^2$ に比例(ex-vcav-prism-cone)底面の形によらず $\frac13Bh$
$x$ 軸の周りの回転体、$y=f(x)$ が与えられている円板法 $\pi\int f^2$、図形が 2 曲線の間なら $\pi\int(f^2-g^2)$図形が軸をまたぐなら不可(ex-vcav-straddle)。遠い方の境界で円板にする
$y$ 軸の周りの回転体、$y=f(x)$ が与えられている殻の方法 $2\pi\int xf$(thm-vcav-shell)$x\ge0$ の側にあること。逆関数で円板法にもできるが手間がかかる(ex-vcav-sine)
図形の面積と重心が分かっているPappus–Guldinus(thm-vcav-pappus)図形が軸の片側にあること
2 つの円柱の共通部分など断面が単純になる向きを探す(ex-vcav-steinmetz)どの向きでも体積は同じ
形の分かった立体と断面積が等しいCavalieri の原理(ex-vcav-archimedes)同じ向きの平行な平面で比べる
軸を含む平面で扇形に切りたくなるそのままでは不可(ex-vcav-fan)各点に軸からの距離を掛ける

さらに先へ

  • $n$ 次元の球の体積も Cavalieri の原理で求まる。半径 $1$ の $n$ 次元球の体積を $\omega_n$ とすると、断面は半径 $\sqrt{1-x^2}$ の $(n-1)$ 次元球なので $\omega_n=\omega_{n-1}\int_{-1}^1(1-x^2)^{(n-1)/2}\,dx$ である。
  • 極座標や円柱座標への変数変換では、ex-vcav-fan で現れた「厚さの比 $\rho$」が Jacobi 行列式として現れる。

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する