三重积分的对称性与应用

对称性在三重积分中的应用

关于坐标面的对称性

定理1

设积分区域 Ω\Omega 关于 xOyxOy 坐标面对称(即若 (x,y,z)Ω(x,y,z) \in \Omega,则 (x,y,z)Ω(x,y,-z) \in \Omega),Ω1\Omega_1Ω\OmegaxOyxOy 面上方的部分(z0z \geq 0)。

  • f(x,y,z)=f(x,y,z)f(x, y, -z) = -f(x, y, z)(即 ff 关于 zz奇函数),则 Ωf(x,y,z)dv=0\iiint_\Omega f(x,y,z) \, dv = 0

  • f(x,y,z)=f(x,y,z)f(x, y, -z) = f(x, y, z)(即 ff 关于 zz偶函数),则 Ωf(x,y,z)dv=2Ω1f(x,y,z)dv\iiint_\Omega f(x,y,z) \, dv = 2\iiint_{\Omega_1} f(x,y,z) \, dv

几何解释
推论
证明
符号说明

类似地,有关于 yOzyOz 面(对应 xx 的奇偶性)和 xOzxOz 面(对应 yy 的奇偶性)的对称性结论。

关于原点的对称性

Ω\Omega 关于原点对称(即若 (x,y,z)Ω(x,y,z) \in \Omega,则 (x,y,z)Ω(-x,-y,-z) \in \Omega):

  • f(x,y,z)=f(x,y,z)f(-x,-y,-z) = -f(x,y,z),则 Ωfdv=0\iiint_\Omega f \, dv = 0
  • f(x,y,z)=f(x,y,z)f(-x,-y,-z) = f(x,y,z),则 Ωfdv=2Ω1fdv\iiint_\Omega f \, dv = 2\iiint_{\Omega_1} f \, dv

轮换对称性

若积分区域 Ω\Omega 关于 x,y,zx, y, z 具有轮换对称性(即交换任意两个变量后区域不变,如球体 x2+y2+z2R2x^2+y^2+z^2 \leq R^2),则:

Ωf(x)dv=Ωf(y)dv=Ωf(z)dv\iiint_\Omega f(x) \, dv = \iiint_\Omega f(y) \, dv = \iiint_\Omega f(z) \, dv

因此: Ω[f(x)+f(y)+f(z)]dv=3Ωf(x)dv\iiint_\Omega [f(x) + f(y) + f(z)] \, dv = 3\iiint_\Omega f(x) \, dv

物理应用

1. 空间物体的质量

质量

设空间物体占有区域 Ω\Omega,体密度为 ρ(x,y,z)\rho(x,y,z),则物体的质量为: M=Ωρ(x,y,z)dvM = \iiint_\Omega \rho(x,y,z) \, dv

2. 质心

定理2

空间物体的质心坐标 (xˉ,yˉ,zˉ)(\bar{x}, \bar{y}, \bar{z})

xˉ=1MΩxρ(x,y,z)dv,yˉ=1MΩyρ(x,y,z)dv,zˉ=1MΩzρ(x,y,z)dv\bar{x} = \frac{1}{M}\iiint_\Omega x\rho(x,y,z) \, dv, \quad \bar{y} = \frac{1}{M}\iiint_\Omega y\rho(x,y,z) \, dv, \quad \bar{z} = \frac{1}{M}\iiint_\Omega z\rho(x,y,z) \, dv

其中 M=Ωρ(x,y,z)dvM = \iiint_\Omega \rho(x,y,z) \, dv 是总质量。

若物体均匀(ρ\rho 为常数),则质心称为形心xˉ=1VΩxdv,yˉ=1VΩydv,zˉ=1VΩzdv\bar{x} = \frac{1}{V}\iiint_\Omega x \, dv, \quad \bar{y} = \frac{1}{V}\iiint_\Omega y \, dv, \quad \bar{z} = \frac{1}{V}\iiint_\Omega z \, dv 其中 VV 是区域 Ω\Omega 的体积。

几何解释
推论
证明
符号说明

3. 转动惯量

定理3

空间物体关于各坐标轴的转动惯量:

  • 关于 xx 轴:Ix=Ω(y2+z2)ρ(x,y,z)dvI_x = \iiint_\Omega (y^2 + z^2)\rho(x,y,z) \, dv
  • 关于 yy 轴:Iy=Ω(x2+z2)ρ(x,y,z)dvI_y = \iiint_\Omega (x^2 + z^2)\rho(x,y,z) \, dv
  • 关于 zz 轴:Iz=Ω(x2+y2)ρ(x,y,z)dvI_z = \iiint_\Omega (x^2 + y^2)\rho(x,y,z) \, dv
  • 关于原点:IO=Ω(x2+y2+z2)ρ(x,y,z)dv=12(Ix+Iy+Iz)I_O = \iiint_\Omega (x^2 + y^2 + z^2)\rho(x,y,z) \, dv = \frac{1}{2}(I_x + I_y + I_z)
几何解释
推论
证明
符号说明

4. 引力

空间物体对位于点 (x0,y0,z0)(x_0, y_0, z_0) 处质量为 mm 的质点的引力 F=(Fx,Fy,Fz)\vec{F} = (F_x, F_y, F_z)

Fx=GmΩρ(x,y,z)(xx0)[(xx0)2+(yy0)2+(zz0)2]3/2dvF_x = Gm\iiint_\Omega \frac{\rho(x,y,z)(x-x_0)}{[(x-x_0)^2+(y-y_0)^2+(z-z_0)^2]^{3/2}} \, dv Fy=GmΩρ(x,y,z)(yy0)[(xx0)2+(yy0)2+(zz0)2]3/2dvF_y = Gm\iiint_\Omega \frac{\rho(x,y,z)(y-y_0)}{[(x-x_0)^2+(y-y_0)^2+(z-z_0)^2]^{3/2}} \, dv Fz=GmΩρ(x,y,z)(zz0)[(xx0)2+(yy0)2+(zz0)2]3/2dvF_z = Gm\iiint_\Omega \frac{\rho(x,y,z)(z-z_0)}{[(x-x_0)^2+(y-y_0)^2+(z-z_0)^2]^{3/2}} \, dv

其中 GG 为引力常数。

典型例题

例题1:利用对称性

计算 Ω(x+y+z)2dv\iiint_\Omega (x + y + z)^2 \, dv,其中 Ω:x2+y2+z2R2\Omega: x^2+y^2+z^2 \leq R^2

参考答案(4 个标签)
三重积分对称性轮换对称性球体

解题思路: 展开后利用对称性和轮换对称性简化。

详细步骤

  1. (x+y+z)2=x2+y2+z2+2xy+2yz+2xz(x+y+z)^2 = x^2+y^2+z^2 + 2xy + 2yz + 2xz
  2. xyxy 关于 xx 是奇函数,区域关于 yOzyOz 面对称 xydv=0\Rightarrow \iiint xy \, dv = 0
  3. 同理 yzdv=0\iiint yz \, dv = 0xzdv=0\iiint xz \, dv = 0
  4. 由轮换对称性:x2dv=y2dv=z2dv\iiint x^2 dv = \iiint y^2 dv = \iiint z^2 dv
  5. (x2+y2+z2)dv=3x2dv\iiint (x^2+y^2+z^2) dv = 3\iiint x^2 dv,而由球坐标知 (x2+y2+z2)dv=4πR55\iiint (x^2+y^2+z^2) dv = \frac{4\pi R^5}{5}
  6. 因此 (x+y+z)2dv=4πR55\iiint (x+y+z)^2 dv = \frac{4\pi R^5}{5}

答案4πR55\frac{4\pi R^5}{5}

例题2:求均匀球体的转动惯量

求半径为 RR、密度为常数 ρ\rho 的均匀球体关于其直径的转动惯量。

参考答案(4 个标签)
三重积分转动惯量球体球坐标

解题思路: 取直径为 zz 轴,计算 Iz=Ω(x2+y2)ρdvI_z = \iiint_\Omega (x^2+y^2)\rho \, dv

详细步骤

  1. 球坐标下:x2+y2=r2sin2ϕx^2+y^2 = r^2\sin^2\phi,体积元素 r2sinϕdrdϕdθr^2\sin\phi \, dr \, d\phi \, d\theta
  2. Iz=ρ02πdθ0πdϕ0Rr2sin2ϕr2sinϕdrI_z = \rho \int_0^{2\pi} d\theta \int_0^\pi d\phi \int_0^R r^2\sin^2\phi \cdot r^2\sin\phi \, dr =ρ2π0πsin3ϕdϕ0Rr4dr= \rho \cdot 2\pi \int_0^\pi \sin^3\phi \, d\phi \int_0^R r^4 dr
  3. 0πsin3ϕdϕ=20π/2sin3ϕdϕ=223=43\int_0^\pi \sin^3\phi \, d\phi = 2\int_0^{\pi/2} \sin^3\phi \, d\phi = 2 \cdot \frac{2}{3} = \frac{4}{3}
  4. 0Rr4dr=R55\int_0^R r^4 dr = \frac{R^5}{5}
  5. Iz=ρ2π43R55=8πρR515I_z = \rho \cdot 2\pi \cdot \frac{4}{3} \cdot \frac{R^5}{5} = \frac{8\pi\rho R^5}{15}
  6. 总质量 M=ρ43πR3M = \rho \cdot \frac{4}{3}\pi R^3,故 Iz=25MR2I_z = \frac{2}{5}MR^2

答案Iz=25MR2I_z = \frac{2}{5}MR^2(其中 MM 为球体总质量)

例题3:求圆锥体的形心

求高为 hh、底面半径为 RR 的均匀正圆锥体的形心。

参考答案(5 个标签)
三重积分形心圆锥体柱坐标对称性

解题思路: 圆锥体关于 zz 轴对称,形心在 zz 轴上,只需计算 zˉ\bar{z}

详细步骤

  1. 设圆锥顶点在原点,轴为 zz 轴,底面在 z=hz=h 处。圆锥面方程:z=hRx2+y2z = \frac{h}{R}\sqrt{x^2+y^2}
  2. 体积 V=13πR2hV = \frac{1}{3}\pi R^2 h
  3. 柱坐标下:0θ2π,0ρR,hRρzh0 \leq \theta \leq 2\pi, \, 0 \leq \rho \leq R, \, \frac{h}{R}\rho \leq z \leq h
  4. Ωzdv=02πdθ0Rρdρhρ/Rhzdz=2π0Rρ12[h2h2ρ2R2]dρ=πh20R(ρρ3R2)dρ=πh2(R22R24)=πh2R24\iiint_\Omega z \, dv = \int_0^{2\pi} d\theta \int_0^R \rho \, d\rho \int_{h\rho/R}^h z \, dz = 2\pi \int_0^R \rho \cdot \frac{1}{2}\left[h^2 - \frac{h^2\rho^2}{R^2}\right] d\rho = \pi h^2 \int_0^R \left(\rho - \frac{\rho^3}{R^2}\right) d\rho = \pi h^2\left(\frac{R^2}{2} - \frac{R^2}{4}\right) = \frac{\pi h^2 R^2}{4}
  5. zˉ=πh2R2/4πR2h/3=3h4\bar{z} = \frac{\pi h^2 R^2 / 4}{\pi R^2 h / 3} = \frac{3h}{4}

答案:形心在对称轴上,距顶点 3h4\frac{3h}{4} 处(即距底面 h4\frac{h}{4} 处)。


练习题

练习1

计算 Ωxyzdv\iiint_\Omega xyz \, dv,其中 Ω:x2+y2+z21,x0,y0,z0\Omega: x^2+y^2+z^2 \leq 1, \, x \geq 0, \, y \geq 0, \, z \geq 0(第一卦限的单位球)。

参考答案(4 个标签)
三重积分球坐标第一卦限乘积函数

解题思路: 第一卦限的单位球,球坐标下 ϕ\phiθ\theta 范围都是 [0,π/2][0, \pi/2]

详细步骤

  1. 球坐标:0r1,0ϕπ/2,0θπ/20 \leq r \leq 1, \, 0 \leq \phi \leq \pi/2, \, 0 \leq \theta \leq \pi/2
  2. xyz=r3sin2ϕcosϕsinθcosθxyz = r^3\sin^2\phi\cos\phi\sin\theta\cos\theta
  3. Ωxyzdv=0π/2sinθcosθdθ0π/2sin3ϕcosϕdϕ01r5dr\iiint_\Omega xyz \, dv = \int_0^{\pi/2} \sin\theta\cos\theta \, d\theta \int_0^{\pi/2} \sin^3\phi\cos\phi \, d\phi \int_0^1 r^5 dr =121416=148= \frac{1}{2} \cdot \frac{1}{4} \cdot \frac{1}{6} = \frac{1}{48}

答案148\frac{1}{48}

练习2

求均匀球体 x2+y2+z2R2x^2+y^2+z^2 \leq R^2 关于原点的转动惯量。

参考答案(4 个标签)
三重积分转动惯量球体球坐标

解题思路: 关于原点的转动惯量 IO=(x2+y2+z2)ρdvI_O = \iiint (x^2+y^2+z^2)\rho \, dv

详细步骤

  1. 球坐标下 x2+y2+z2=r2x^2+y^2+z^2 = r^2
  2. IO=ρ02πdθ0πsinϕdϕ0Rr2r2dr=ρ2π2R55=4πρR55I_O = \rho \int_0^{2\pi} d\theta \int_0^\pi \sin\phi \, d\phi \int_0^R r^2 \cdot r^2 dr = \rho \cdot 2\pi \cdot 2 \cdot \frac{R^5}{5} = \frac{4\pi\rho R^5}{5}
  3. 总质量 M=43πρR3M = \frac{4}{3}\pi\rho R^3,故 IO=35MR2I_O = \frac{3}{5}MR^2

答案IO=35MR2I_O = \frac{3}{5}MR^2


总结

本文出现的符号

符号类型读音/说明在本文中的含义
ρ(x,y,z)\rho(x,y,z)体密度volume density空间物体在点(x,y,z)处的体密度
MM质量mass空间物体的总质量
(xˉ,yˉ,zˉ)(\bar{x},\bar{y},\bar{z})质心坐标centroid coordinates质心的坐标
Ix,Iy,IzI_x, I_y, I_z转动惯量moment of inertia关于各坐标轴的转动惯量
VV体积volume空间区域的体积
GG引力常数gravitational constant万有引力常数

中英对照

中文术语英文术语音标说明
对称性symmetry/ˈsɪmətri/区域关于坐标面或原点的对称性质
轮换对称性cyclic symmetry/ˈsaɪklɪk ˈsɪmətri/交换变量后区域不变的性质
质心center of mass/ˈsɛntər əv mæs/质量中心
形心centroid/ˈsɛntrɔɪd/均匀物体的质心,即几何中心
转动惯量moment of inertia/ˈmoʊmənt əv ɪˈnɜːrʃə/物体转动惯性的量度
引力gravitational force/ˌɡrævɪˈteɪʃənəl fɔːrs/物体间的万有引力
体密度volume density/ˈvɒljuːm ˈdɛnsəti/单位体积的质量