3重積分(triple integrals)

**図 7.16:** 3重積分
$\begin{figure}\begin{center} \includegraphics[width=7cm]{CALCFIG/Fig7-6-1.eps} \end{center}\end{figure}$

2重積分と同様に3重積分も累次積分として表わすことができます．空間の有界閉領域をとし，の平面への正射影を $\Omega_{xy}$ とすると，は

$\displaystyle T = \{(x,y,z) : (x,y) \in \Omega_{xy}, \psi_{1}(x,y) \leq z \leq \psi_{2}(x,y) \}$

で表わせます．これより

$\displaystyle \iiint_{T}f(x,y,z)dxdydz = \iint_{\Omega_{xy}}[\int_{\psi_{1}(x,y)}^{\psi_{2}(x,y)}f(x,y,z)dz]dxdy$

が成り立ちます．さらに，

$\displaystyle \Omega = \{(x,y) : a \leq x \leq b, \phi_{1}(x) \leq y \leq \phi_{2}(x) \}$

ならば，

$\displaystyle \iiint_{T}f(x,y,z)dxdydz = \int_{a}^{b}dx \int_{\phi_{1}(x)}^{\phi_{2}(x)}dy\int_{\psi_{1}(x,y)}^{\psi_{2}(x,y)}f(x,y,z)dz$

が成り立ちます．

3重積分についても2重積分の場合と同様に次の定理が成り立ちます．

定理 7..6

変数変換 $\Phi : \left\{\begin{array}{c} x = \phi(u,v,w)\\ y = \psi(u,v,w)\\ z = \xi(u,v,w) \end{array}\right.$ によって，

空間の有界な閉領域

と

空間の有界な閉領域 $T^{\prime}$ が1対1に対応し， $\phi,\psi,\xi$ が

級で

$\displaystyle J = J(u,v,w) = \left\vert\begin{array}{ccc} x_{u}&x_{v}&x_{w}\\ y_{u}&y_{v}&y_{w}\\ z_{u}&z_{v}&z_{w} \end{array}\right \vert \neq 0$

とする．このとき，

が

上で連続ならば，

$\displaystyle \iiint_{T}f(x,y,z)dxdydz = \iiint_{T^{\prime}}f(\phi(u,v,w),\psi(u,v,w),\xi(u,v,w))\vert J\vert dudvdw$

が成り立つ．

ここで変数変換の中で特によく用いられる円柱座標変換と球面座標変換について学びます．

系 7..1

[円柱座標変換]

$\displaystyle \left\{\begin{array}{c} x = r\cos{\theta}\\ y = r\sin{\theta} ... ...eta}&0\\ \sin{\theta}& r\cos{\theta} & 0\\ 0&0&1 \end{array}\right \vert = r$

より

$\displaystyle \iiint_{T}f(x,y,z)dxdydz = \iiint_{T^{\prime}}f(r\cos{\theta},r\sin{\theta},z)dzrdrd\theta$

例題 7..11

$\displaystyle{T = \{(x,y,z) : x^2 + y^2 \leq 1, -\sqrt{1 - x^2 - y^2} \leq z \leq \sqrt{1 - x^2 - y^2}\}}$ のとき．次の値を求めてみましょう．

$\displaystyle I = \iiint_{T} \frac{1}{\sqrt{1 - x^2 - y^2}} dxdydz$

解円柱座標変換を用いると，は

$\displaystyle T' = \{(r,\theta, z) : 0 \leq r \leq 1, 0 \leq \theta \leq 2\pi, - \sqrt{1 - r^2} \leq z \leq \sqrt{1 - r^2} \}$

にうつるから，

$\displaystyle I$	$\displaystyle =$	$\displaystyle \iiint_{T'}\frac{1}{\sqrt{1 - r^2}} rdzdrd\theta$
	$\displaystyle =$	$\displaystyle \int_{0}^{2\pi}d\theta \int_{0}^{1} dr \int_{- \sqrt{1 - r^2}}^{\sqrt{1 - r^2}} \frac{r}{\sqrt{1 - r^2}}dz$
	$\displaystyle =$	$\displaystyle \int_{0}^{2\pi} \int_{0}^{1} 2r dr = 2\pi \cdot \left[r^{2} \right ]_{0}^{1} = 2\pi \ensuremath{ \blacksquare}$

**図 7.17:** 円柱座標変換
$\begin{figure}\begin{center} \includegraphics[width=5.6cm]{CALCFIG/Fig7-6-2-1.eps} \end{center}\end{figure}$

系 7..2

[球面座標変換]

$\displaystyle \left\{\begin{array}{c} x = \rho\sin{\phi}\cos{\theta}\\ y = \rho\sin{\phi}\sin{\theta}\\ z = \rho\cos{\phi} \end{array}\right .$

$\displaystyle \vert J(\rho,\phi,\theta)\vert = \left\vert\begin{array}{ccc} \si... ...\\ \cos{\phi}&-\rho\sin{\phi}&0 \end{array}\right \vert = \rho^{2} \sin{\phi}$

より

$\displaystyle \iiint_{T}f(x,y,z)dxdydz = \iiint_{T^{\prime}}f(\rho\sin{\phi}\co... ...o\sin{\phi}\sin{\theta},\rho\cos{\phi})\rho^{2} \sin{\phi} d\rho d\phi d\theta$

例題 7..12

$\displaystyle{T = \{(x,y,z) : x^2 + y^2 + z^2 \leq 1\}}$ のとき，次の値を求めてみましょう．

$\displaystyle I = \iiint_{T}\frac{dxdydz}{\sqrt{1 - x^2 - y^2 - z^2}}$

解球面座標変換を用いると，は

$\displaystyle T' = \{(\rho,\phi,\theta) : 0 \leq \rho \leq 1, 0 \leq \phi \leq \pi, 0 \leq \theta \leq 2\pi \}$

にうつるから，

$\displaystyle I$	$\displaystyle =$	$\displaystyle \iiint_{T'}\frac{\rho^{2} \sin{\phi}}{\sqrt{1 - \rho^2}} d\rho d\phi d\theta$
	$\displaystyle =$	$\displaystyle \int_{0}^{2\pi}d\theta \int_{0}^{\pi} \sin{\phi} d\phi \int_{0}^{1} \frac{\rho^2}{\sqrt{1 - \rho^2}} d \rho$
	$\displaystyle =$	$\displaystyle 2\pi \cdot 2 \cdot \left[\frac{1}{2}\left(- \rho \sqrt{1 - \rho^2} + \sin^{-1}{\rho} \right) \right ]_{0}^{1}$
	$\displaystyle =$	$\displaystyle 4\pi \cdot \frac{1}{2} \cdot \frac{\pi}{2} = \pi^2 \ensuremath{ \blacksquare}$

**図 7.18:** 球面座標変換
$\begin{figure}\begin{center} \includegraphics[width=5.6cm]{CALCFIG/Fig7-6-2-2.eps} \end{center}\end{figure}$

2重積分の応用でを面積として考えました．同じようにして，ここでは，を体積として考えることができます．

もし密度 $\rho$ がそれぞれの点で連続的に変化し， $\rho = f(x,y,z)$ ならば，立体の質量は

$\displaystyle M = \iiint_{T}f(x,y,z)dxdydz$

で表わせます．特に密度 $\rho = f(x,y,z) = 1$ のとき

$\displaystyle \iiint_{T}dxdydz$

は立体

の体積と考えることができます．

例題 7..13

次の図7.19のような三角錐

の体積を求めてみましょう．

**図 7.19:** 三角錐
$\begin{figure}\begin{center} \includegraphics[width=6cm]{CALCFIG/Fig7-6-3.eps} \end{center}\end{figure}$

解 3重積分を計算するためにを適当な座標軸面に正射影します．ここでは平面に正射影します．ここでV-simpleを用いると

$\displaystyle \Omega_{xy} = \{(x,y) : 0 \leq x \leq 1, 0 \leq y \leq 1-x \}$

これより

$\displaystyle T = \{(x,y,z) : 0 \leq x \leq 1, 0 \leq y \leq 1-x, 0 \leq z \leq 1-x-y \}$

よって

$\displaystyle V$	$\displaystyle =$	$\displaystyle \iint_{\Omega_{xy}}[\int_{0}^{1-x-y}dz]dxdy = \int_{0}^{1}\int_{0}^{1-x}(1-x-y)dydx$
	$\displaystyle =$	$\displaystyle \int_{0}^{1} \left[(1-x)y - \frac{y^2}{2}\right ]_{0}^{1-x} dx = \int_{0}^{1}\frac{(1-x)^{2}}{2}dx$
	$\displaystyle =$	$\displaystyle - \left[\frac{(1-x)^{3}}{6}\right ]_{0}^{1} = \frac{1}{6}$

となります． $\blacksquare$

実際，体積を求めるのが目的ならば，底面積 $\times$ 高さ $\displaystyle{\times \frac{1}{3}}$ よりすぐに求まります．そこで，ここでは，三角錐の各点における密度が異なる場合を考えてみます．

例題 7..14

この三角錐

の密度が $\rho(x,y,z) = xy$ のときの質量を求めてみましょう．

解

$\displaystyle M = \iiint_{T}xydxdydz = \int_{0}^{1}\int_{0}^{1-x}\int_{0}^{1-x-y}xy dzdydx$

ここで

$\displaystyle \int_{0}^{1-x-y}xydz = xy(1-x-y) = x(1-x)y - xy^2$

より

$\displaystyle \int_{0}^{1-x}\int_{0}^{1-x-y}xydzdy$	$\displaystyle =$	$\displaystyle \int_{0}^{1-x}(x(1-x)y - xy^2)dy$
	$\displaystyle =$	$\displaystyle \left[\frac{1}{2}x(1-x)y^2 - \frac{1}{3}xy^3\right ]_{0}^{1-x} = \frac{1}{6}x(1-x)^{3}$

よって

$\displaystyle M$	$\displaystyle =$	$\displaystyle \int_{0}^{1}\frac{1}{6}x(1-x)^{3}dx = \frac{1}{6}\int_{0}^{1}(x-3x^2 + 3x^3 - x^4)dx$
	$\displaystyle =$	$\displaystyle \frac{1}{6}\left[\frac{1}{2}x^2 - x^3 + \frac{3}{4}x^4 -\frac{1}{5}x^5\right ]_{0}^{1} = \frac{1}{120} \ensuremath{ \blacksquare}$

体積，質量と求めてきたので最後に重心を求めてみましょう．まず平面上に個の点 ${\rm P}_{1}, {\rm P}_{2}, \ldots, {\rm P}_{n}$ からなる質点系があって，点 ${\rm P}_{i}$ の座標は ${\rm P}_{i}(x_{i},y_{i})$ ．また ${\rm P}_{i}$ の質量が $M_{i}$ で与えられているとします．このときこれらの点がつりあうように軸に垂直な線 $x = \bar x$ を引くと，この線に対して反時計回りのモーメントと時計回りのモーメントはつりあっているはずです．よって

$\displaystyle \sum_{i=1}^{n}M_{i}(x_{i} - \bar x) = 0$

書き直すと

$\displaystyle \bar x = \frac{\sum_{i=1}^{n}M_{i}x_{i}}{\sum_{i=1}^{n}M_{i}}$

同様に

$\displaystyle \bar y = \frac{\sum_{i=1}^{n}M_{i}y_{i}}{\sum_{i=1}^{n}M_{i}}$

このようにして表わされた点 $(\bar x, \bar y)$ を 質点系の重心(centroid) といいます．

**図 7.20:** 質点系の重心
$\begin{figure}\begin{center} \includegraphics[width=8cm]{CALCFIG/Fig7-6-4.eps} \end{center}\end{figure}$

平面上の有界閉領域 $\Omega$ の各点 ${\rm P}(x,y)$ に密度 $\rho(x,y)$ が与えられているときは， $\Omega$ の分割 $\Delta : \Omega_{1},\Omega_{2},\ldots,\Omega_{n}$ をとり，それぞれの $\Omega_{i}$ の中から任意に点 ${\rm P}_{i}(x_{i},y_{i})$ を選び，質点系 ${\rm P}_{1}, {\rm P}_{2}, \ldots, {\rm P}_{n}$ の重心 $(x_{\Delta},y_{\Delta})$ を考えます． $\Omega_{i}$ の面積を $\Delta S_{i}$ とおくと

$\displaystyle x_{\Delta} = \frac{\sum_{i=1}^{n}\rho({\rm P}_{i})x_{i}\Delta S_{... ...ho({\rm P}_{i})y_{i}\Delta S_{i}}{\sum_{i=1}^{n}\rho({\rm P}_{i})\Delta S_{i}}$

となります．ここで $\vert\Delta\vert \rightarrow 0$ とすると $x_{\Delta},y_{\Delta}$ はそれぞれ次の $\bar x, \bar y$ に収束します．

$\displaystyle \bar x = \frac{1}{M}\iint_{\Omega}\rho(x,y)xdxdy, \bar y = \frac{1}{M}\iint_{\Omega}\rho(x,y)ydxdy$

ただし， $M = \iint_{\Omega}\rho(x,y)dxdy$ で質量(mass)を表わします．

空間の閉領域に，密度 $\rho = \rho(x,y,z)$ が与えられている場合にも，同様にして重心 $(\bar x, \bar y, \bar z)$ が次のように定まります．

$\displaystyle \bar x = \frac{1}{M}\iiint_{T}\rho(x,y,z)xdxdydz, \bar y = \frac{1}{M}\iiint_{T}\rho(x,y,z)ydxdydz$

$\displaystyle \bar z = \frac{1}{M}\iiint_{T}\rho(x,y,z)zdxdydz, M = \iiint_{T}\rho(x,y,z)dxdydz$

例題 7..15

さて，重心の復習が終わったところで上の三角錐

の重心を求めてみましょう．

解

$\displaystyle \bar x$	$\displaystyle =$	$\displaystyle \frac{1}{M}\iiint_{T}xyxdxdydz$
	$\displaystyle =$	$\displaystyle \frac{1}{M}\int_{0}^{1}\int_{0}^{1-x}\int_{0}^{1-x-y}x^2 ydzdydx$
	$\displaystyle =$	$\displaystyle \frac{1}{M}\int_{0}^{1}\int_{0}^{1-x}x^{2}y(1-x-y)dydx$
	$\displaystyle =$	$\displaystyle \frac{1}{M}\int_{0}^{1}\left[(\frac{x^{2}(1-x)y^{2}}{2} - \frac{x^{2}y^{3}}{3})\right ]_{0}^{1-x}dx$
	$\displaystyle =$	$\displaystyle \frac{1}{M}\int_{0}^{1}\frac{x^{2}(1-x)^{3}}{6} dx = \frac{1}{6M}\int_{0}^{1}(x^{2} - 3x^{3} + 3x^{4} - x^{5}) dx$
	$\displaystyle =$	$\displaystyle \frac{1}{6M}\left[\frac{x^{3}}{3} - \frac{3x^{4}}{4} + \frac{3x^{5}}{5} - \frac{x^{6}}{6} \right ]_{0}^{1} = \frac{1}{360M}$