柱坐标下的三重积分计算

为什么需要柱坐标

当积分区域是旋转体(如圆柱、圆锥、旋转抛物面等),或被积函数含有 x2+y2x^2 + y^2 时,柱坐标可以极大简化积分限的确定和积分计算。柱坐标本质上是”极坐标 + 直角坐标 zz“的混合坐标系。

柱坐标变换

柱坐标变换公式

直角坐标与柱坐标的关系: {x=ρcosθy=ρsinθz=z\begin{cases} x = \rho \cos\theta \\ y = \rho \sin\theta \\ z = z \end{cases}

其中 ρ0\rho \geq 0 是点到 zz 轴的距离(极径),θ[0,2π)\theta \in [0, 2\pi) 是极角,zRz \in \mathbb{R} 是竖坐标。

雅可比行列式: J=(x,y,z)(ρ,θ,z)=cosθρsinθ0sinθρcosθ0001=ρJ = \frac{\partial(x,y,z)}{\partial(\rho,\theta,z)} = \begin{vmatrix} \cos\theta & -\rho\sin\theta & 0 \\ \sin\theta & \rho\cos\theta & 0 \\ 0 & 0 & 1 \end{vmatrix} = \rho

因此体积元素:dv=dxdydz=ρdρdθdzdv = dx\,dy\,dz = \rho \, d\rho \, d\theta \, dz

柱坐标下的三重积分公式

Ωf(x,y,z)dv=Ωf(ρcosθ,ρsinθ,z)ρdρdθdz\iiint_\Omega f(x,y,z) \, dv = \iiint_{\Omega'} f(\rho\cos\theta, \rho\sin\theta, z) \cdot \rho \, d\rho \, d\theta \, dz

常见柱坐标区域

1. 圆柱体

圆柱体

Ω:x2+y2R2,azb\Omega: x^2 + y^2 \leq R^2, \quad a \leq z \leq b

柱坐标下: Ω:0ρR,0θ2π,azb\Omega': 0 \leq \rho \leq R, \quad 0 \leq \theta \leq 2\pi, \quad a \leq z \leq b

Ωfdv=abdz02πdθ0Rf(ρcosθ,ρsinθ,z)ρdρ\iiint_\Omega f \, dv = \int_a^b dz \int_0^{2\pi} d\theta \int_0^R f(\rho\cos\theta, \rho\sin\theta, z) \rho \, d\rho

2. 圆锥体

圆锥体

顶点在原点、开口向上的圆锥 z=x2+y2z = \sqrt{x^2+y^2} 与平面 z=hz = h 围成的区域:

柱坐标下:zzρ\rhohh,或 ρ\rho00zz

方法一(先对 zz 积分): 0ρh,0θ2π,ρzh0 \leq \rho \leq h, \quad 0 \leq \theta \leq 2\pi, \quad \rho \leq z \leq h

方法二(先对 ρ\rho 积分,先二后一): 0zh,0θ2π,0ρz0 \leq z \leq h, \quad 0 \leq \theta \leq 2\pi, \quad 0 \leq \rho \leq z

3. 旋转抛物面

旋转抛物面

z=x2+y2z = x^2 + y^2z=hz = h 围成的区域:

柱坐标下:z=ρ2z = \rho^2 0ρh,0θ2π,ρ2zh0 \leq \rho \leq \sqrt{h}, \quad 0 \leq \theta \leq 2\pi, \quad \rho^2 \leq z \leq h

适用场景总结

典型例题

例题1:圆柱区域

计算 Ω(x2+y2)dv\iiint_\Omega (x^2 + y^2) \, dv,其中 Ω:x2+y24,0z3\Omega: x^2 + y^2 \leq 4, \, 0 \leq z \leq 3

参考答案(4 个标签)
三重积分柱坐标圆柱体旋转体

解题思路: 圆柱体 + 被积函数含 x2+y2x^2+y^2,柱坐标是最佳选择。

详细步骤

  1. 柱坐标变换:x=ρcosθ,y=ρsinθ,z=zx = \rho\cos\theta, \, y = \rho\sin\theta, \, z = z,体积元素 ρdρdθdz\rho \, d\rho \, d\theta \, dz
  2. 积分区域:0ρ2,0θ2π,0z30 \leq \rho \leq 2, \, 0 \leq \theta \leq 2\pi, \, 0 \leq z \leq 3
  3. 被积函数:x2+y2=ρ2x^2+y^2 = \rho^2
  4. Ω(x2+y2)dv=03dz02πdθ02ρ2ρdρ=32π02ρ3dρ=6π[ρ44]02=6π4=24π\iiint_\Omega (x^2+y^2) dv = \int_0^3 dz \int_0^{2\pi} d\theta \int_0^2 \rho^2 \cdot \rho \, d\rho = 3 \cdot 2\pi \cdot \int_0^2 \rho^3 d\rho = 6\pi \cdot \left[\frac{\rho^4}{4}\right]_0^2 = 6\pi \cdot 4 = 24\pi

答案24π24\pi

例题2:圆锥区域

计算 Ωzdv\iiint_\Omega z \, dv,其中 Ω\Omega 由圆锥 z=x2+y2z = \sqrt{x^2+y^2} 和平面 z=1z = 1 围成。

参考答案(4 个标签)
三重积分柱坐标圆锥体先二后一

解题思路: 圆锥区域,被积函数仅含 zz,可以用柱坐标,也可以用先二后一法。

详细步骤

方法一:柱坐标(先对 zz 积分)

  1. 柱坐标下:0ρ1,0θ2π,ρz10 \leq \rho \leq 1, \, 0 \leq \theta \leq 2\pi, \, \rho \leq z \leq 1
  2. Ωzdv=02πdθ01ρdρρ1zdz=2π01ρ[z22]ρ1dρ=2π01ρ(12ρ22)dρ=π01(ρρ3)dρ=π(1214)=π4\iiint_\Omega z \, dv = \int_0^{2\pi} d\theta \int_0^1 \rho \, d\rho \int_\rho^1 z \, dz = 2\pi \int_0^1 \rho \cdot \left[\frac{z^2}{2}\right]_\rho^1 d\rho = 2\pi \int_0^1 \rho\left(\frac{1}{2} - \frac{\rho^2}{2}\right) d\rho = \pi \int_0^1 (\rho - \rho^3) d\rho = \pi\left(\frac{1}{2} - \frac{1}{4}\right) = \frac{\pi}{4}

方法二:先二后一

  1. 0z10 \leq z \leq 1,截面 Dz:x2+y2z2D_z: x^2+y^2 \leq z^2,面积 S(z)=πz2S(z) = \pi z^2
  2. Ωzdv=01zπz2dz=π01z3dz=π4\iiint_\Omega z \, dv = \int_0^1 z \cdot \pi z^2 dz = \pi \int_0^1 z^3 dz = \frac{\pi}{4}

答案π4\frac{\pi}{4}

例题3:旋转抛物面与球面围成

计算 Ωzdv\iiint_\Omega z \, dv,其中 Ω\Omegaz=2x2y2z = \sqrt{2-x^2-y^2}z=x2+y2z = x^2+y^2 围成。

参考答案(5 个标签)
三重积分柱坐标旋转抛物面球面交线

解题思路: 先求交线确定投影区域,再用柱坐标计算。

详细步骤

  1. 交线:2x2y2=x2+y2\sqrt{2-x^2-y^2} = x^2+y^2,令 r2=x2+y2r^2 = x^2+y^2,则 2r2=r2\sqrt{2-r^2} = r^2 平方得 2r2=r4r4+r22=0(r2+2)(r21)=0r=12-r^2 = r^4 \Rightarrow r^4+r^2-2=0 \Rightarrow (r^2+2)(r^2-1)=0 \Rightarrow r=1
  2. 交线在 z=1z=1 处,投影区域 D:x2+y21D: x^2+y^2 \leq 1
  3. 柱坐标下:0ρ1,0θ2π,ρ2z2ρ20 \leq \rho \leq 1, \, 0 \leq \theta \leq 2\pi, \, \rho^2 \leq z \leq \sqrt{2-\rho^2}
  4. Ωzdv=02πdθ01ρdρρ22ρ2zdz=2π01ρ[z22]ρ22ρ2dρ\iiint_\Omega z \, dv = \int_0^{2\pi} d\theta \int_0^1 \rho \, d\rho \int_{\rho^2}^{\sqrt{2-\rho^2}} z \, dz = 2\pi \int_0^1 \rho \cdot \left[\frac{z^2}{2}\right]_{\rho^2}^{\sqrt{2-\rho^2}} d\rho =π01ρ[(2ρ2)ρ4]dρ=π01(2ρρ3ρ5)dρ=π[11416]=π712=7π12= \pi \int_0^1 \rho[(2-\rho^2) - \rho^4] d\rho = \pi \int_0^1 (2\rho - \rho^3 - \rho^5) d\rho = \pi\left[1 - \frac{1}{4} - \frac{1}{6}\right] = \pi \cdot \frac{7}{12} = \frac{7\pi}{12}

答案7π12\frac{7\pi}{12}


练习题

练习1

计算 Ωx2+y2dv\iiint_\Omega \sqrt{x^2+y^2} \, dv,其中 Ω:x2+y29,1z2\Omega: x^2+y^2 \leq 9, \, -1 \leq z \leq 2

参考答案(3 个标签)
三重积分柱坐标圆柱体

解题思路: 圆柱体 + 被积函数含 x2+y2\sqrt{x^2+y^2},用柱坐标。

详细步骤

  1. 柱坐标:0ρ3,0θ2π,1z20 \leq \rho \leq 3, \, 0 \leq \theta \leq 2\pi, \, -1 \leq z \leq 2
  2. Ωx2+y2dv=12dz02πdθ03ρρdρ=32π03ρ2dρ=6π9=54π\iiint_\Omega \sqrt{x^2+y^2} dv = \int_{-1}^2 dz \int_0^{2\pi} d\theta \int_0^3 \rho \cdot \rho \, d\rho = 3 \cdot 2\pi \cdot \int_0^3 \rho^2 d\rho = 6\pi \cdot 9 = 54\pi

答案54π54\pi

练习2

计算 Ω(x2+y2)dv\iiint_\Omega (x^2+y^2) \, dv,其中 Ω\Omegax2+y2=2zx^2+y^2 = 2zz=2z = 2 围成。

参考答案(3 个标签)
三重积分柱坐标旋转抛物面

解题思路: 旋转抛物面区域,柱坐标计算。

详细步骤

  1. 交线:x2+y2=4x^2+y^2 = 4,投影区域 D:x2+y24D: x^2+y^2 \leq 4
  2. 柱坐标:0ρ2,0θ2π,ρ22z20 \leq \rho \leq 2, \, 0 \leq \theta \leq 2\pi, \, \frac{\rho^2}{2} \leq z \leq 2
  3. Ω(x2+y2)dv=02πdθ02ρdρρ2/22ρ2dz=2π02ρ3(2ρ22)dρ=2π02(2ρ3ρ52)dρ=2π[ρ42ρ612]02=2π(8163)=2π83=16π3\iiint_\Omega (x^2+y^2) dv = \int_0^{2\pi} d\theta \int_0^2 \rho \, d\rho \int_{\rho^2/2}^2 \rho^2 \, dz = 2\pi \int_0^2 \rho^3\left(2 - \frac{\rho^2}{2}\right) d\rho = 2\pi \int_0^2 (2\rho^3 - \frac{\rho^5}{2}) d\rho = 2\pi\left[\frac{\rho^4}{2} - \frac{\rho^6}{12}\right]_0^2 = 2\pi\left(8 - \frac{16}{3}\right) = 2\pi \cdot \frac{8}{3} = \frac{16\pi}{3}

答案16π3\frac{16\pi}{3}


总结

本文出现的符号

符号类型读音/说明在本文中的含义
ρ\rho极径rho点到 z 轴的距离,ρ0\rho \geq 0
θ\theta极角theta在 xOy 平面上与 x 轴正方向的夹角
zz竖坐标z coordinate与直角坐标的 z 相同
ρdρdθdz\rho \, d\rho \, d\theta \, dz柱坐标体积元素cylindrical volume element柱坐标系下的体积微元

中英对照

中文术语英文术语音标说明
柱坐标cylindrical coordinates/sɪˈlɪndrɪkəl koʊˈɔːrdɪneɪts/用极径、极角和竖坐标表示点的坐标系
旋转体solid of revolution/ˈsɒlɪd əv ˌrɛvəˈluːʃən/由平面曲线绕轴旋转而成的立体
圆柱体cylinder/ˈsɪləndər/圆柱面与两平面围成的立体
圆锥体cone/koʊn/圆锥面与平面围成的立体
旋转抛物面paraboloid of revolution/pəˈræbəlɔɪd əv ˌrɛvəˈluːʃən/抛物线绕轴旋转而成的曲面
雅可比行列式Jacobian determinant/dʒəˈkoʊbiən dɪˈtɜːrmɪnənt/坐标变换的偏导数行列式