25-第二十五章 曲面积分
依赖于
被以下题目直接调用
正文部分
25.1.1 例题
设 S 为 z=x2+y2 被 x2+y2=2ax 割下的部分,求
I=∬S(x2y2+y2z2+z2x2)dS.
解 1
见图 25.1,在直角坐标系中计算
zx=zx,zy=zy,
故
1+zx2+zy2=2,
从而
I=∬x2+y2≤2ax[x2y2+(x2+y2)2]2dxdy.
用极坐标变换求上述二重积分,则
I=2∫−π/2π/2dθ∫02acosθ(r4cos2θsin2θ+r4)rdr=2∫−π/2π/2(cos2θsin2θ+1)(61r602acosθ)dθ=62(2a)6∫−π/2π/2cos6θ(cos2θsin2θ+1)dθ=8292πa6.
解 2
用参数式计算。z=x2+y2 在球坐标系中的方程为
φ=4π,因此 S 的参数方程为
x=21rcosθ,y=21rsinθ,z=21r,(r,θ)∈D.
又 S 的边界线
z=x2+y2,x2+y2=2ax
的球坐标表示为
φ=4π,r2sin2φ=2arsinφcosθ.
于是
D={(r,θ)∣−2π≤θ≤2π, 0≤r≤22acosθ}.
计算得
E=2r2,F=0,G=1.
最后得到
I=∫−π/2π/2dθ∫022acosθ(41r4cos2θsin2θ+41r4)2rdr=8292πa6.
25.1.2 例题
求上半球面 z=a2−x2−y2 被 x2+y2=ax 截取部分的面积与质心坐标,其中 a>0。
解
见图 22.11,这就是 Viviani 体的上表面。由于
zx=a2−x2−y2−x,zy=a2−x2−y2−y,
故所求面积
S=∬x2+y2≤axa2−x2−y2adxdy=∫−π/2π/2dθ∫0acosθa2−r2ardr=(π−2)a2.
下面求质心坐标,由对称性知 y0=0,且
x0=(π−2)a21∬SxdS=(π−2)a22∬x2+y2≤axy≥0a2−x2−y2axdxdy=(π−2)a2∫0π/2cosθdθ∫0acosθa2−r2r2dr(r=acost)=(π−2)a2∫0π/2cosθdθ∫π/2θasinta2cos2t(−asint)dt=π−2a∫0π/2cosθ(2π−θ−sinθcosθ)dθ=3(π−2)2a.
又
z0=(π−2)a21∬SzdS=(π−2)a21∬x2+y2≤axa2−x2−y2aa2−x2−y2dxdy=4(π−2)πa.
25.1.3 练习题
求 ∬Sz2dS,其中
- S 为 z2=x2+y2 在 x2+y2+z2=R2 内的上半部分 (z≥0);
- S 为 x=rsinαcosθ, y=rsinαsinθ, z=rcosα, 0≤r≤a, 0≤θ≤2π。
解
- 令 r=x2+y2,则在锥面上 z=r,而由 r2+z2=R2 得
0≤r≤R/2。又
dS=1+zr2rdrdθ=2rdrdθ.
因而
∬Sz2dS=∫02π∫0R/2r2⋅2rdrdθ=8π2R4.
- 已给参数方程中
rr=(sinαcosθ,sinαsinθ,cosα),rθ=(−rsinαsinθ,rsinαcosθ,0),
故 ∣rr×rθ∣=rsinα,于是
∬Sz2dS=∫02π∫0ar2cos2α⋅rsinαdrdθ=2πa4cos2αsinα.
求 ∬S(x+y+z)dS,其中 S 为上半单位球面 x2+y2+z2=1 (z≥0)。
解
由对称性,∬SxdS=∬SydS=0。上半球面上 z=1−x2−y2,故
zdS=dxdy,于是
∬S(x+y+z)dS=∬SzdS=∬x2+y2≤1dxdy=π.
求 ∬S(x+y+z)2dS,其中 S 为单位球面 x2+y2+z2=1。
解
展开后得
(x+y+z)2=x2+y2+z2+2xy+2yz+2zx.
单位球面对坐标平面都对称,故三项交叉项积分都为 0;又在球面上
x2+y2+z2=1,所以
∬S(x+y+z)2dS=∬SdS=4π.
求 ∬S(x4−y4+y2z2−z2x2+1)dS,其中 S 是锥面 z2=x2+y2 被柱面 x2+y2=2x 割下的部分。
解
在锥面上 z2=x2+y2,故
x4−y4+y2z2−z2x2+1=x4−y4+y2(x2+y2)−(x2+y2)x2+1=1.
所求即为该曲面部分的面积。它由上下两片组成;在每一片上可写成
z=±x2+y2,故 dS=2dxdy。投影区域是
x2+y2≤2x,其面积为 π,所以
∬S(x4−y4+y2z2−z2x2+1)dS=22∬x2+y2≤2xdxdy=22π.
求 ∬S∣xyz∣dS,其中
- S 为 ∣x∣+∣y∣+∣z∣=1;
- S 是抛物面 z=x2+y2 被 z=1 割下的部分。
解
- 由对称性,只算第一卦限内的平面片 x+y+z=1 即可。在该片上
z=1−x−y,dS=3dxdy,故
∬S∣xyz∣dS=83∫01dx∫01−xxy(1−x−y)dy=83⋅1201=153.
- 取极坐标 x=rcosθ, y=rsinθ, z=r2。则
0≤r≤1,0≤θ≤2π,且
dS=1+4r2rdrdθ.
因而
∬S∣xyz∣dS=∫02π∫01r4∣cosθsinθ∣1+4r2rdrdθ=2∫01r51+4r2dr=4201255−1.
求 ∬S(x2+y2+z2)dS,S 为 ∣x∣+∣y∣+∣z∣=a。
解
仍由对称性,只算第一卦限内的平面片 x+y+z=a。在该片上
dS=3dxdy,故
∬S(x2+y2+z2)dS=83∫0adx∫0a−x[x2+y2+(a−x−y)2]dy=23a4.
求
F(t)=∬x2+y2+z2=t2f(x,y,z)dS,
其中
f(x,y,z)={x2+y2,0,z≥x2+y2,z<x2+y2.
解
用球坐标
x=tsinφcosθ, y=tsinφsinθ, z=tcosφ。
只有当 z≥x2+y2,即 0≤φ≤π/4 时被积函数非零,于是
F(t)=∫02π∫0π/4t2sin2φ⋅t2sinφdφdθ=2πt4∫0π/4sin3φdφ=6π(8−52)t4.
求
F(x,y,z,t)=∬Sf(ξ,η,ζ)dS,
其中 S 为
(x−ξ)2+(y−η)2+(z−ζ)2=t2,
t>0,(x,y,z) 为满足 x2+y2+z2>a>0 的定点,
f(ξ,η,ζ)={1,0,ξ2+η2+ζ2<a2,ξ2+η2+ζ2≥a2.
解
记 d=x2+y2+z2(>a),并取极轴沿着 OP。设球面上点
(ξ,η,ζ) 与 OP 的夹角为 φ,则
ξ2+η2+ζ2=d2+t2−2dtcosφ.
条件 ξ2+η2+ζ2<a2 等价于
cosφ>2dtd2+t2−a2.
因而只有在 ∣t−d∣<a 时有交,且交成一个球冠。其高为
h=t−2dd2+t2−a2=2da2−(t−d)2,
故所求面积
F(x,y,z,t)=2πth=dπt[a2−(t−d)2].
综上,
F(x,y,z,t)=⎩⎨⎧dπt[a2−(t−d)2],0,∣t−d∣<a,∣t−d∣≥a.
求 ∬S(1+x+y)2dS,其中 S 为立体
x+y+z≤1, x≥0, y≥0, z≥0 的表面。
解
该立体的表面由四个面组成:x=0, y=0, z=0, x+y+z=1。前三个面的贡献分别为
∬x=0(1+x+y)2dS=∫01dy∫01−y(1+y)2dz=1−ln2,
∬y=0(1+x+y)2dS=1−ln2,
∬z=0(1+x+y)2dS=∫01dx∫01−x(1+x+y)2dy=ln2−21.
在斜面 z=1−x−y 上,dS=3dxdy,故其贡献为
3∫01dx∫01−x(1+x+y)2dy=3(ln2−21).
相加即得
∬S(1+x+y)2dS=23−3+(3−1)ln2.
求 ∬Sz∣x∣dS,其中 S 是柱面
x2+y2=2ay 被锥面 z=x2+y2 和平面 z=2a 所截下的部分。
解
取参数方程
x=asin2θ,y=2asin2θ,z=z,0≤θ≤π.
由 x2+y2=2ay 知这正是所给柱面,且
∣rθ×rz∣=2a.
又由锥面与平面给出
2asinθ≤z≤2a。于是
∬Sz∣x∣dS=4a2∫0π∣sinθcosθ∣(∫2asinθ2azdz)dθ=8a2∫0π/2sinθcosθlnsinθ1dθ=2a2.
求 ∬S(x2+y2)dS,S 是锥面
z=x2+y2 被 z=1 割下的部分。
解
在锥面上令 z=r,0≤r≤1,则
dS=2rdrdθ,x2+y2=r2.
因而
∬S(x2+y2)dS=∫02π∫01r2⋅2rdrdθ=2π2.
25.2.1 例题
设 Σ 为上半单位球面 z=1−(x2+y2),取内侧,求
I=∬Σdydz+dzdx+dxdy.
解 1
用直角坐标系计算:
I=∬Σdydz+∬Σdzdx+∬Σdxdy=I1+I2+I3.
计算 I1:Σ=Σ1+Σ2,其中
Σ1:x=1−y2−z2,(y,z)∈Dyz={y2+z2≤1, z≥0},取后侧,
Σ2:x=−1−y2−z2,(y,z)∈Dyz,取前侧.
则
I1=∬Σ1dydz+∬Σ2dydz=−∬Dyzdydz+∬Dyzdydz=0.
同理有
I2=0,I3=−∬Dxydxdy=−π.
最后得到
I=−π.
解 2
用参数方程
Σ: x=sinφcosθ,y=sinφsinθ,z=cosφ,0≤φ≤2π,0≤θ≤2π.
计算行列式
A=∂(φ,θ)∂(y,z)=sin2φcosθ,B=∂(φ,θ)∂(z,x)=sin2φsinθ,
C=∂(φ,θ)∂(x,y)=sinφcosφ.
因为 (A,B,C) 的方向与上半球面 S 内侧的法线方向相反,故积分号前取 “−” 号,得到
I=−∬D(sin2φcosθ+sin2φsinθ+sinφcosφ)dφdθ=−∫0π/2dφ∫02π[sin2φ(cosθ+sinθ)+sinφcosφ]dθ=−2π∫0π/2sinφcosφdφ=−π.
求
I=∬Σ(z+x)dydz+(x+y)dzdx+(y+z)dxdy,
其中 Σ 是由 x2+y2=1, z=1 及三个坐标平面围成的立体在第一卦限的部分的表面,取外侧。
解
见图 25.2,因为 Σ 分块较多(需分 5 块),不便于用参数式,故应在直角坐标系中计算。记
Σ=Σ1+Σ2+Σ3+Σ4+Σ5,其中
Σ1 为圆柱面,Σ2 为下底面,Σ3 为上底面,
Σ4 为左侧面 y=0,Σ5 为右侧面 x=0。则有
I=∬Σ(z+x)dydz+∬Σ(x+y)dzdx+∬Σ(y+z)dxdy=I1+I2+I3.
计算 I1,有
I1=i=1∑5∬Σi(x+z)dydz.
因为 Σ2,Σ3,Σ4 在 yOz 平面上的投影面积为零,于是
∬Σi(x+z)dydz=0,i=2,3,4,
而
Σ1:x=1−y2,(y,z)∈Dyz={0≤y≤1, 0≤z≤1},取前侧,
Σ5:x=0,(y,z)∈Dyz,取后侧.
于是
I1=∬Σ1(x+z)dydz+∬Σ5(x+z)dydz=∬Dyz(1−y2+z)dydz−∬Dyzzdydz=∫01∫011−y2dydz=4π.
用类似的方法可求出 I2=4π。
然后求
I3=i=1∑5∬Σi(y+z)dxdy.
显然有
∬Σi(y+z)dxdy=0,i=1,4,5,
而
Σ2:z=0,(x,y)∈Dxy={x2+y2≤1, x≥0, y≥0},取下侧,
Σ3:z=1,(x,y)∈Dxy,取上侧.
于是
I3=∬Σ2(y+z)dxdy+∬Σ3(y+z)dxdy=−∬Dxyydxdy+∬Dxy(1+y)dxdy=∬Dxydxdy=4π.
最后得到
I=I1+I2+I3=43π.
注
我们已经知道,利用第二型曲线积分可以计算平面图形的面积。类似地,利用第二型曲面积分也可以计算空间立体的体积。
25.2.2 例题
求
I=∬Σxyz(y2z2+z2x2+x2y2)dS,
其中 Σ 为第一卦限中的球面 x2+y2+z2=a2 (x≥0, y≥0, z≥0)。
解
见图 25.3,不论用参数式或直角坐标式,直接计算均相当复杂。取 Σ 的上侧,则 (x,y,z) 处的单位外法向量为
(ax,ay,az).
利用公式 (25.4),
I=a∬Σ(y3z3ax+z3x3ay+x3y3az)dS=a∬Σy3z3dydz+z3x3dzdx+x3y3dxdy=3a∬Σx3y3dxdy=3a∬Dxyx3y3dxdy.
其中
Dxy={x2+y2≤a2, x≥0, y≥0}.
作极坐标变换得
I=3a∫0π/2dθ∫0ar7sin3θcos3θdr=643a9∫0π/2sin32θdθ=321a9.
求
I=∬Σ(y−z)dydz+(z−x)dzdx+(x−y)dxdy,
其中 Σ 是球面 x2+y2+z2=2Rx 被柱面 x2+y2=2rx (0<r<R) 截下的位于 z≥0 的部分,取外侧。
解
见图 25.4,改写球面方程为 (x−R)2+y2+z2=R2,其外侧的法向量为
n=(Rx−R,Ry,Rz).
由公式 (25.4),有
I=∬Σ(y−z,z−x,x−y)⋅ndS=R1∬Σ[(y−z)(x−R)+(z−x)y+(x−y)z]dS=∬Σ(z−y)dS.
由于 Σ 关于 xOz 平面对称,而函数 y 是奇函数,于是
I=∬ΣzdS=∬x2+y2≤2rx2Rx−x2−y22Rx−x2−y2Rdxdy=R⋅πr2=πRr2.
25.2.3 练习题
设曲面 Σ 的方程为 z=z(x,y), (x,y)∈D,且 z(x,y) 在 D 中连续可微,证明:
∬ΣP(x,y,z)dydz+Q(x,y,z)dzdx+R(x,y,z)dxdy
=±∬D{−Pzx−Qzy+R}z=z(x,y)dxdy.
当 Σ 取上侧时符号取 $+$'',当 $\Sigma$ 取下侧时符号取 −”。
解
由两类曲面积分的关系,
∬ΣPdydz+Qdzdx+Rdxdy=∬Σ(Pcosα+Qcosβ+Rcosγ)dS.
若 Σ 取上侧,写成 z=z(x,y),则
(cosα,cosβ,cosγ)=1+zx2+zy2(−zx,−zy,1),dS=1+zx2+zy2dxdy,
故
∬ΣPdydz+Qdzdx+Rdxdy=∬D(−Pzx−Qzy+R)z=z(x,y)dxdy.
若取下侧,法向量反向,故整体只差一个负号。
若 Σ 分块光滑,且关于 xOy 平面对称,f(x,y,z) 在 Σ 上连续,且满足
f(x,y,z)=−f(x,y,−z),问:
∬Σf(x,y,z)dS=2∬Σ1f(x,y,z)dS
还是等于 0?(其中 Σ1 是 Σ 在 xOy 平面以上的部分。)
解
设 Σ2 是 Σ1 关于 xOy 平面的对称部分。反射
(x,y,z)↦(x,y,−z) 不改变面积元 dS,而
f(x,y,−z)=−f(x,y,z),所以
∬Σ2f(x,y,z)dS=−∬Σ1f(x,y,z)dS.
故
∬Σf(x,y,z)dS=0.
若 Σ 分块光滑,且关于 xOy 平面对称,R(x,y,z) 在 Σ 上连续,满足
R(x,y,z)=−R(x,y,−z),问:
∬ΣR(x,y,z)dxdy=2∬Σ1R(x,y,z)dxdy
还是等于 0?(其中 Σ1 是 Σ 在 xOy 平面以上的部分,取侧与 Σ 取侧相一致。)
解
写成
∬ΣR(x,y,z)dxdy=∬ΣR(x,y,z)cos(n,z)dS.
反射到下半部分时,R(x,y,z) 变号,而 cos(n,z) 也随所取侧变号,二者乘积不变,故
∬ΣR(x,y,z)dxdy=2∬Σ1R(x,y,z)dxdy.
设 Σ 是平面 Π 内的一个有界区域,其面积为 S,Π 取上侧的法向量为 n,且
cos(n,z)=μ。证明:Σ 在 xOy 平面上的投影的面积为 μS,并利用这个结果重新计算例题 21.4.2。
解
对平面内的面积元有 dxdy=cos(n,z)dS=μdS。积分即得投影面积
Sxy=∬Σdxdy=μ∬ΣdS=μS.
例题 21.4.2 中,平面 2x+y+z=0 的上侧法向量可取
6(2,1,1),故 μ=1/6。把 z=−2x−y 代入椭球方程
3x2+y2+2z2=1,得投影椭圆
14x2+12xy+9y2≤6.
其面积为
Sxy=14⋅9−62π⋅6=π52.
因而原椭圆面积
S=μSxy=6π52=2π53.
求 I1=∬Σzdxdy, I2=∬Σz2dxdy,
Σ 是球面 x2+y2+z2=a2,取外侧。
解
在球面 x2+y2+z2=a2 的外侧,
cos(n,z)=az.
故
I1=∬Σzdxdy=a1∬Σz2dS,I2=∬Σz2dxdy=a1∬Σz3dS.
由对称性,
∬Σz2dS=31∬Σ(x2+y2+z2)dS=31a2⋅4πa2=34πa4,
而 ∬Σz3dS=0,故
I1=34πa3,I2=0.
计算
∬Sxzdydz+yxdzdx+yzdxdy,
S 是圆柱面 x2+y2=1 在 −1≤z≤1 及 x≥0 的部分,取前侧。
解
取参数方程
x=cosθ,y=sinθ,z=z,−2π≤θ≤2π, −1≤z≤1.
此时前侧法向量对应的面积向量为
ndS=(cosθ,sinθ,0)dθdz.
因而
∬Sxzdydz+yxdzdx+yzdxdy=∫−11dz∫−π/2π/2(zcos2θ+sin2θcosθ)dθ=2∫−π/2π/2sin2θcosθdθ=34.
求
∬Sx(z2−y2)dydz+y(x2−z2)dzdx+z(y2−x2)dxdy,
其中 S 是 y2+z2=1 被 x=0, x=1, z+y=0, z−y=0 截取的上方部分,取外侧。
解
取参数方程
x=x,y=cosθ,z=sinθ,0≤x≤1, 4π≤θ≤43π.
这里上方部分对应 z≥∣y∣。外侧面积向量为
ndS=(0,cosθ,sinθ)dxdθ.
所以
∬Sx(z2−y2)dydz+y(x2−z2)dzdx+z(y2−x2)dxdy=∫01dx∫π/43π/4x2cos2θdθ=31∫π/43π/4cos2θdθ=−31.
25.3.1 例题
求
I=∬Σ4xzdydz−2yzdzdx+(1−z2)dxdy,
其中 Σ 是曲线 z=ey (0≤y≤a) 绕 z 轴旋转生成的旋转面,取下侧。
解
Σ 的方程为
z=ex2+y2(x2+y2≤a2).
直接计算比较复杂,考虑用 Gauss 公式。由于 Σ 不闭,需要添加辅助面
Σ1:z=ea,x2+y2≤a2,
取上侧。见图 25.5,设 Σ 与 Σ1 围成的区域为 D。令
P=4xz,Q=−2yz,R=1−z2,
则
∂x∂P+∂y∂Q+∂z∂R=0.
由公式 (25.5),得
I=(∬Σ+∬Σ1−∬Σ1)(1−z2)dxdy=−∬Σ1(1−z2)dxdy=(e2a−1)∬x2+y2≤a2dxdy=(e2a−1)πa2.
计算曲面积分
I=∬S(ax2+by2+cz2)3/2xdydz+ydzdx+zdxdy,
其中 S 是球面 x2+y2+z2=1,取外侧 (a>0, b>0, c>0)。
解 1
记
P(x,y,z)=(ax2+by2+cz2)3/2x,Q(x,y,z)=(ax2+by2+cz2)3/2y,
R(x,y,z)=(ax2+by2+cz2)3/2z,
则在不包含原点的任何区域上
∂x∂P+∂y∂Q+∂z∂R=0.
为了利用 Gauss 公式,对充分小的 ε>0,作闭曲面
Sε={ax2+by2+cz2=ε2},
取外侧。由 Gauss 公式
I=∬Sε(ax2+by2+cz2)3/2xdydz+ydzdx+zdxdy=ε31∬Sεxdydz+ydzdx+zdxdy.
上述积分在 Sε 的外侧。再一次用 Gauss 公式,则
I=ε33∭ax2+by2+cz2≤ε2dxdydz=ε33⋅34π⋅abcε3=abc4π.
解 2(不用 Gauss 公式而直接计算)
利用单位球面的参数方程
x=sinφcosθ,y=sinφsinθ,z=cosφ,
0≤φ≤π,0≤θ≤2π.
计算得到
A=∂(φ,θ)∂(y,z)=sin2φcosθ,B=∂(φ,θ)∂(z,x)=sin2φsinθ,
C=∂(φ,θ)∂(x,y)=sinφcosφ.
容易看出,(A,B,C) 的方向与单位球面外侧法线方向相同,故积分号前取 “+” 号。由 (25.3) 得
I=8∫0π/2∫0π/2(asin2φcos2θ+bsin2φsin2θ+ccos2φ)3/2sinφdφdθ.
先计算对 φ 的积分,令 cosφ=t,则
∫0π/2(asin2φcos2θ+bsin2φsin2θ+ccos2φ)3/2sinφdφ=∫01[(acos2θ+bsin2θ)−(acos2θ+bsin2θ−c)t2]3/2dt=acos2θ+bsin2θ1[(acos2θ+bsin2θ)−(acos2θ+bsin2θ−c)t2]1/2tt=0t=1=c1⋅acos2θ+bsin2θ1.
最后得到
I=c8∫0π/2acos2θ+bsin2θdθ=c8∫0+∞a+bt2dt(t=tanθ)=c8⋅ab1arctan(abt)0+∞=abc4π.
注
利用 Gauss 公式来计算曲面积分在很多情况下是一种有效的手段,但要注意使用 Gauss 公式的条件,要弄清楚在什么情况下要“挖洞”(即用封闭曲面把 P,Q,R 无定义或不可微的点围住)以及选择什么曲面“挖洞”计算更简便。
利用 Gauss 公式可导出用曲面积分表示 R3 中具有逐片光滑边界的有界闭区域 Ω 的体积公式
∣Ω∣=∬∂Ωxdydz=∬∂Ωydzdx=∬∂Ωzdxdy=31∬∂Ωxdydz+ydzdx+zdxdy.
上述积分在 ∂Ω 的外侧进行。由两类曲面积分之间的关系,又有
∣Ω∣=31∬∂Ω(xcosα+ycosβ+zcosγ)dS,
其中 (cosα,cosβ,cosγ) 为 ∂Ω 的单位外法向量。如果 ∂Ω 有参数表示
x=x(u,v),y=y(u,v),z=z(u,v),(u,v)∈D,
则
V=31∬Dxxuxvyyuyvzzuzvdudv.
特别地,若一个立体的表面在球坐标系中由方程
r=r(φ,θ),(φ,θ)∈D
给出,则
V=31∬Dr3(φ,θ)sinφdφdθ.
25.3.2 练习题
利用 Gauss 公式计算积分:
- ∬Sy(x−z)dydz+z2dzdx+(y2+xz)dxdy,
其中 S 是正立方体 {(x,y,z)∣0≤x≤a, 0≤y≤a, 0≤z≤a} (a>0) 的表面,取内侧;
- ∬Σ(x3+x)dydz+(y2−zx)dzdx+(z3+z)dxdy,
其中 Σ 是球面 x2+y2+z2=2z,取外侧;
- 设 A1=x3−x2y+z3, A2=xy2+y3, A3=xz+z2,
Σ 是由 yOz 平面上的抛物线 z=1−y2 与 z=0 所围成的平面区域绕 z 轴旋转后所得的旋转体的表面,取外侧。试求
∬Σ(∂y∂A3−∂z∂A2)dydz+(∂z∂A1−∂x∂A3)dzdx+(∂x∂A2−∂y∂A1)dxdy.
解
- 取
F=(y(x−z),z2,y2+xz).
则
divF=∂x∂[y(x−z)]+∂y∂z2+∂z∂(y2+xz)=x+y.
表面取内侧,故所求积分等于外侧通量的相反数,从而
∬S⋯=−∭[0,a]3(x+y)dxdydz=−a4.
- 取
F=(x3+x, y2−zx, z3+z).
则
divF=3x2+3z2+2.
球体为 (x,y,z)=(u,v,w+1) 下的单位球 u2+v2+w2≤1,故
∭Ωx2dV=∭B1u2dV=154π,
∭Ωz2dV=∭B1(w+1)2dV=∭B1w2dV+∭B11dV=154π+34π.
因而
∬Σ⋯=3⋅154π+3(154π+34π)+2⋅34π=15124π.
- 记 A=(A1,A2,A3)。所给积分正是
∬Σ(∇×A)⋅ndS.
由 Gauss 公式及 div(∇×A)=0,
∬Σ(∇×A)⋅ndS=∭Ωdiv(∇×A)dV=0.
先添加辅助面,再用 Gauss 公式计算下列曲面积分:
- ∬Σ(x2cosα+y2cosβ+z2cosγ)dS,
其中 Σ 是锥面 z2=x2+y2 在 0≤z≤h 的一段,
(cosα,cosβ,cosγ) 为 Σ 上的单位法向量,其方向为下方;
- ∬Σx3dydz+y3dzdx+z3dxdy,
其中 Σ 为球面 x2+y2+z2=a2 之上半部分,取上侧;
∬Σ(a3x3+y3z3)dydz+(b3y3+z3x3)dzdx+(c3z3+x3y3)dxdy,
其中 Σ 为椭球面
a2x2+b2y2+c2z2=1, x≥0,取后侧。
解
- 令 F=(x2,y2,z2),并在上口添加圆盘
S0: z=h, x2+y2≤h2。锥面与圆盘围成的立体为
Ω={x2+y2≤z2, 0≤z≤h}。因锥面取下方,恰为该立体的外侧,故
∬Σ(x2cosα+y2cosβ+z2cosγ)dS=∬ΣF⋅ndS.
又
divF=2x+2y+2z,
其中 x,y 的体积分因对称性为零,故
∭ΩdivFdV=2∭ΩzdV=2∫0h∫02π∫0zzrdrdθdz=2πh4.
而
∬S0F⋅ndS=∬x2+y2≤h2h2dxdy=πh4.
因而
∬ΣF⋅ndS=2πh4−πh4=−2πh4.
- 令 F=(x3,y3,z3),并补上底圆盘
S0: z=0, x2+y2≤a2。因 S0 上通量为零,故
∬Σx3dydz+y3dzdx+z3dxdy=∭x2+y2+z2≤a2, z≥03(x2+y2+z2)dV.
用球坐标即得
∬Σ⋯=3∫02π∫0π/2∫0aρ2⋅ρ2sinφdρdφdθ=56πa5.
- 记
F=(a3x3+y3z3, b3y3+z3x3, c3z3+x3y3).
用圆盘 S0: x=0, b2y2+c2z2≤1 封闭之。由于
F⋅(−1,0,0)x=0=−y3z3,
其在 S0 上积分由奇偶性为零。于是闭曲面的外侧通量就是曲面
x≥0 这一半椭球面的外侧通量。又
divF=a33x2+b33y2+c33z2.
在半椭球内有
∭x2dV=152πa3bc,∭y2dV=152πab3c,∭z2dV=152πabc3,
故外侧通量为
52π(ab+ac+bc).
题中取后侧,与外侧相反,所以
∬Σ⋯=−52π(ab+ac+bc).
F(x,y,z) 是定义在 R3 上的光滑函数,且 F(x,y,z)=0 是一个以原点为顶点的锥面 Σ。
如果 Σ 与平面 Π:Ax+By+Cz=D 围成一个锥体,证明:此锥体的体积
V=31SH,
其中 S 为平面 Π 上锥底部分的面积,H 为顶点到锥底的高。
解
取向量场 F=(x,y,z),则 divF=3。对锥体应用 Gauss 公式得
3V=∬∂ΩF⋅ndS.
锥侧上的母线都过原点,故径向向量 (x,y,z) 与锥侧切平面平行,从而
F⋅n=0。在底面 Π:Ax+By+Cz=D 上,单位法向量取与外侧一致的方向,
则对底面上每一点都有
F⋅n=H,
这里 H 是原点到平面 Π 的距离。于是
3V=∬底面HdS=HS,
即
V=31SH.
求由曲面 (x2+y2+z2)2=a2xy 所围成的立体的体积。
解
用球坐标
x=ρsinφcosθ,y=ρsinφsinθ,z=ρcosφ.
曲面方程化为
ρ4=a2ρ2sin2φcosθsinθ=2a2ρ2sin2φsin2θ.
因而内部区域满足
0≤ρ≤2asinφsin2θ,
其中 sin2θ≥0,即
θ∈[0,π/2]∪[π,3π/2]。由对称性,
V=2∫0π/2∫0π∫02asinφsin2θρ2sinφdρdφdθ=32a3(∫0πsin4φdφ)(∫0π/2(sin2θ)3/2dθ)=162πa3B(45,21).
求
∬Σ(x3+y3)dydz+(x3+2x2y)dzdx−x2zdxdy,
其中 Σ 是单叶双曲面 x2+y2−z2=1 在 0≤z≤3 的部分,取外侧。
解
记
F=(x3+y3, x3+2x2y, −x2z).
取上、下底圆盘把曲面封闭。由
divF=4x2
得闭曲面的外侧通量为
∭Ω4x2dV=4∫03∫02π∫01+z2r2cos2θ⋅rdrdθdz=5243π.
上底圆盘 z=3, x2+y2≤4 上的通量为
∬−x2zdxdy=−43π,
下底 z=0 上通量为 0。故所求曲面积分为
5243π−(−43π)=5443π.
V={(x,y,z)∣x2+y2<z<1},S=∂V,求积分
∬Syzdzdx+(x2+y2)zdxdy,
积分沿外法线方向。
解
令
F=(0, yz, (x2+y2)z).
则
divF=z+x2+y2.
由 Gauss 公式,
∬Syzdzdx+(x2+y2)zdxdy=∭V(z+x2+y2)dV.
用柱坐标 0≤r≤1, r2≤z≤1,得
∭V(z+r2)dV=2π∫01∫r21(z+r2)rdzdr=2π.
求第二型曲面积分
∬Szdydz+cosydzdx+dxdy,
其中 S 为 x2+y2+z2=1 的外侧。
解
令 F=(z, cosy, 1),则
divF=−siny.
单位球关于平面 y=0 对称,而 siny 是奇函数,因此
∬Szdydz+cosydzdx+dxdy=∭x2+y2+z2≤1(−siny)dV=0.
25.3.3 例题
求
I=∮C(y2−z2)dx+(z2−x2)dy+(x2−y2)dz,
其中 C 是立方体 {(x,y,z)∣0≤x≤a, 0≤y≤a, 0≤z≤a} 的表面与平面
x+y+z=23a 的交线,取向从 z 轴正向看去是逆时针方向。
分析
见图 25.6,分六段积分的计算量很大,且 C 也不便于表示为一个统一的参数式。因 C 为闭曲线,且
P=y2−z2,Q=z2−x2,R=x2−y2
连续可微,故考虑用 Stokes 公式。
解
令 Σ 为 x+y+z=23a 被 C 所围的一块,取上侧,则 C 的取向与 Σ 的取侧相容。应用 Stokes 公式 (25.13),
I=∬Σ31∂x∂y2−z231∂y∂z2−x231∂z∂x2−y2dS=31∬Σ[−4(x+y+z)]dS=−34∬Σ23adS=−23a∣Σ∣=−23a⋅433a2=−29a3.
设 Σ 是分片光滑的闭曲面,n 为 Σ 上的单位外法向量,证明:
I=∬Σcos(n,x)∂x∂Pcos(n,y)∂y∂Qcos(n,z)∂z∂RdS=0,
其中分两种情形:(1)P,Q,R 在 Ω 上二阶连续可微,Ω 为 Σ 所围的立体;(2)P,Q,R 在 Σ 上二阶连续可微。
证
对情形 (1) 用 Gauss 公式:
I=∬Σ(∂y∂R−∂z∂Q)dydz+(∂z∂P−∂x∂R)dzdx+(∂x∂Q−∂y∂P)dxdy=∭Ω[∂x∂(∂y∂R−∂z∂Q)+∂y∂(∂z∂P−∂x∂R)+∂z∂(∂x∂Q−∂y∂P)]dxdydz=0.
情形 (2) 参见图 25.7。在 Σ 上任取一条逐段光滑的闭曲线 C,C 分 Σ 为两部分
Σ1 与 Σ2。在 Σ1,Σ2 上分别应用 Stokes 公式,则对于 i=1,2,
∬Σicos(n,x)∂x∂Pcos(n,y)∂y∂Qcos(n,z)∂z∂RdS=∮CiPdx+Qdy+Rdz.
因 Σ1,Σ2 分居 C 两侧,故 C1,C2 为同一条曲线 C,只是它们的定向相反。若记 C1 为 C+,则 C2 为 C−,从而
I=(∬Σ1+∬Σ2)cos(n,x)∂x∂Pcos(n,y)∂y∂Qcos(n,z)∂z∂RdS=(∮C++∮C−)Pdx+Qdy+Rdz=0.
试用 Stokes 公式计算
I=∮C(y2+z2)dx+(z2+x2)dy+(x2+y2)dz,
其中 C 为 x2+y2+z2=2Rx 与 x2+y2=2rx 的交线
(0<r<R, z>0),C 的定向使得 C 所包围的球面上较小区域保持在左边。
解
见图 25.4,设 S 为球面 x2+y2+z2=2Rx 被柱面
x2+y2=2rx 所截部分的外侧。由 Stokes 公式 (25.13),
I=∬SRx−R∂x∂y2+z2Ry∂y∂z2+x2Rz∂z∂x2+y2dS=R2∬S[(y−z)(x−R)+(z−x)y+(x−y)z]dS=2∬S(z−y)dS=2∬SzdS=2R∬x2+y2≤2rxdxdy=2πr2R.
25.3.4 练习题
设 C 是平面 xcosα+ycosβ+zcosγ−p=0 上逐段光滑的闭曲线,C 所界的面积为 S,C 的定向与
(cosα,cosβ,cosγ) 成右手系,试计算积分
∮Cdxcosαxdycosβydzcosγz.
解
记常向量
n=(cosα,cosβ,cosγ),r=(x,y,z).
则题中积分可写成
∮C(n×r)⋅dr.
设 Σ 是 C 在该平面内所围成的区域,其法向量就是 n。由 Stokes 公式,
∮C(n×r)⋅dr=∬Σcurl(n×r)⋅ndS.
因 n 为常向量,直接计算得
curl(n×r)=2n,故
∮Cdxcosαxdycosβydzcosγz=2∬ΣdS=2S.
求 ∮C(y−z)dx+(z−x)dy+(x−y)dz,
C 为 x2+y2=1 与 x+y+z=1 的交线,从 x 轴正向看是逆时针方向。
解
取平面片 Σ: x+y+z=1, x2+y2≤1,其法向量取
3(1,1,1),与题设定向相容。设
F=(y−z, z−x, x−y).
则
curlF=(−2,−2,−2).
该平面片在 xOy 面上的投影是单位圆盘,故
∣Σ∣=cos(n,z)π=3π.
由 Stokes 公式,
∮C(y−z)dx+(z−x)dy+(x−y)dz=(−2,−2,−2)⋅3(1,1,1)∣Σ∣=−6π.
求
∫C(x3+3xz2y)dx+(x3+3y2z)dy+(y3+3z2x)dz,
其中 C 是 z=a2−x2−y2 与 x=y 的交线,自
A(2a,2a,0)到B(−2a,−2a,0).
解
该曲线可参数化为
x=y=2acost,z=asint,0≤t≤π.
代入后
∫C⋯=∫0π[(x3+3xz2y)x′+(x3+3y2z)y′+(y3+3z2x)z′]dt=802a4(15π−16a).
用 Stokes 公式求
∫Cex+z{[(x+1)y2+1]dx+2xydy+xy2dz},
其中 C 是右半柱面 ∣x∣+∣y∣=a (y>0) 与平面 y=z 的交线上从 (−a,0,0) 到 (a,0,0) 的一段 (a>0)。
解
设
A=(−a,0,0),P=(0,a,a),B=(a,0,0).
曲线 C 正是折线 A→P→B。添上线段 BA 后得到闭曲线
Γ,它围成平面 y=z 内的三角形
Σ={(x,s,s)∣−a≤x≤a, 0≤s≤a−∣x∣}.
取参数方程 (x,s,s),其法向量可取
(0,1,−1),与 Γ 的方向 A→P→B→A 相容。设
F=ex+z([(x+1)y2+1], 2xy, xy2),
则
curlF=(0,ex+z,0).
由 Stokes 公式,
∮ΓF⋅dr=∬ΣcurlF⋅ndS=∫−aadx∫0a−∣x∣ex+sds=2(2a−1)ea+e−a.
而在线段 BA 上有 y=z=0,故
∫BAex+z{[(x+1)y2+1]dx+2xydy+xy2dz}=∫a−aexdx=e−a−ea.
因此
∫Cex+z{[(x+1)y2+1]dx+2xydy+xy2dz}=2(2a+1)ea−e−a.
设 C 是空间任一逐段光滑的简单闭曲线,f(x),g(x),h(x) 是任意连续函数,证明:
∮C[f(x)−yz]dx+[g(y)−zx]dy+[h(z)−xy]dz=0.
解
设
F(x)=∫0xf(t)dt,G(y)=∫0yg(t)dt,H(z)=∫0zh(t)dt.
则
d(F(x)+G(y)+H(z)−xyz)=[f(x)−yz]dx+[g(y)−zx]dy+[h(z)−xy]dz.
因而该微分式是全微分,沿任意闭曲线的积分都为 0。
求
∬Σcosα∂x∂x−zcosβ∂y∂x3−yzcosγ∂z∂−3xy2dS,
其中 Σ 是 x2+y2+z2=R2 在 z≥0 的部分,
(cosα,cosβ,cosγ) 是 Σ 下侧的单位法向量。
解
由 Stokes 公式,该曲面积分等于边界圆周 C: x2+y2=R2, z=0 上的线积分。由于 Σ 取下侧,故 C 的正向从 z 轴正向看去是顺时针方向。把 z=0, dz=0 代入得
∮C(x−z)dx+(x3−yz)dy−3xy2dz=∮Cxdx+x3dy=∮Cx3dy.
再由 Green 公式,
∮Cx3dy=−∬x2+y2≤R23x2dxdy=−43πR4.
求 ∮Cydx+zdy+xdz,其中 C 是
x2+y2+z2=a2 与 x+y+z=0 的交线,从 z 轴正向看是逆时针方向。
解
取平面圆盘 Σ: x+y+z=0, x2+y2+z2≤a2。由题设定向,
其单位法向量取
n=3(1,1,1),
因为它的 z 分量为正。设 F=(y,z,x),则
curlF=(−1,−1,−1).
圆盘半径为 a,面积为 πa2,故
∮Cydx+zdy+xdz=∬ΣcurlF⋅ndS=−3πa2.
25.3.5 例题
对于微分式
z(x2y1−x2+z21)dx+xy2zdy+(x2+z2x−xy1)dz,
判断原函数的存在性并求出之。
解 1
容易验证:
∂y∂P=∂x∂Q=−x2y2z,∂z∂Q=∂y∂R=xy21,
∂x∂R=∂z∂P=x2y1+(x2+z2)2z2−x2.
因此该微分式有原函数。根据微分式的特点,为计算简单起见取
z0=0, x0,y0>0,积分路径为
(x0,y0,0)→(x,y0,0)→(x,y,0)→(x,y,z),
则
φ(x,y,z)=∫0z(x2+z2x−xy1)dz+C=arctanxz−xyz+C.
解 2
求原函数时也可用下面求不定积分的方法:由
∂z∂φ=x2+z2x−xy1,
则
φ(x,y,z)=∫(x2+z2x−xy1)dz=arctanxz−xyz+ψ(x,y).
由此得
∂x∂φ=−x2+z2z+x2yz+∂x∂ψ,∂y∂φ=xy2z+∂y∂ψ.
由 ∂x∂φ=P, ∂y∂φ=Q 得
∂x∂ψ=∂y∂ψ=0.
即 ψ(x,y) 为常数,所以
φ(x,y,z)=arctanxz−xyz+C.
25.3.6 练习题
证明:下列微分式为全微分,并求出其原函数:
- (x2−2yz)dx+(y2−2zx)dy+(z2−2xy)dz;
[(x2−y2)2x−x1+2x2]dx+[y1−(x2−y2)2y+3y3]dy+5z3dz.
解
- 直接观察得
d(3x3+y3+z3−2xyz)=(x2−2yz)dx+(y2−2zx)dy+(z2−2xy)dz.
故它是全微分,其原函数可取
3x3+y3+z3−2xyz+C.
- 在每个连通分支上,
d[−2(x2−y2)1−ln∣x∣+ln∣y∣+32x3+43y4+45z4]
=[(x2−y2)2x−x1+2x2]dx+[y1−(x2−y2)2y+3y3]dy+5z3dz.
故它也是全微分,原函数如上。
求 ∫(1,2,3)(6,1,1)yzdx+xzdy+xydz。
解
因为
yzdx+xzdy+xydz=d(xyz),
故
∫(1,2,3)(6,1,1)yzdx+xzdy+xydz=xyz(1,2,3)(6,1,1)=6−6=0.
设 C 是由球面 x2+y2+z2=a2 上的任一点沿任一路径运动到球面
x2+y2+z2=b2 (b>a) 上的任一点的轨迹,C 分段光滑,证明:
∫Cr3(xdx+ydy+zdz)=51(b5−a5),
其中 r=x2+y2+z2。
解
由
dr=rxdx+ydy+zdz,
得
r3(xdx+ydy+zdz)=r4dr=d(5r5).
因而对任意分段光滑路径 C,
∫Cr3(xdx+ydy+zdz)=51r5起点终点=51(b5−a5).
25.4.4 例题
设 ω=dxi1∧dxi2∧⋯∧dxik,则 dω=0。
证
dω=d(1dxi1∧dxi2∧⋯∧dxik)=d1∧dxi1∧dxi2∧⋯∧dxik=0.
设 ω=Pdx+Qdy+Rdz 为 R3 上的 C1 类 1-形式,计算 dω。
解
dω=(dP)∧dx+(dQ)∧dy+(dR)∧dz=(∂x∂Pdx+∂y∂Pdy+∂z∂Pdz)∧dx+(∂x∂Qdx+∂y∂Qdy+∂z∂Qdz)∧dy+(∂x∂Rdx+∂y∂Rdy+∂z∂Rdz)∧dz=(∂y∂R−∂z∂Q)dy∧dz+(∂z∂P−∂x∂R)dz∧dx+(∂x∂Q−∂y∂P)dx∧dy.
设 f∈Λ0 是 C2 类的 0-形式,即 f 是二次连续可微的函数,则 d2f=0。
证
由于
∂xi∂xj∂2f=∂xj∂xi∂2f(i,j=1,2,⋯,n),
所以
d2f=d(df)=d(i=1∑n∂xi∂fdxi)=i=1∑nj=1∑n∂xj∂xi∂2fdxj∧dxi=i<j∑(∂xi∂xj∂2f−∂xj∂xi∂2f)dxi∧dxj=0.
25.4.6 例题
极坐标变换 x=rcosθ, y=rsinθ,有
dx∧dy=∂(r,θ)∂(x,y)dr∧dθ=rdr∧dθ.
证
由
x=rcosθ,y=rsinθ
得
dx=cosθdr−rsinθdθ,dy=sinθdr+rcosθdθ.
于是
dx∧dy=(cosθdr−rsinθdθ)∧(sinθdr+rcosθdθ)=r(cos2θ+sin2θ)dr∧dθ=rdr∧dθ.
25.5.1 习题课教案一例
1. 关于第一类曲面积分在正交变换下的不变性
第一类曲面积分在正交变换下的不变性,即
∬Sf(X)dS=∬Σf(ATU)dΣ,(25.18)
其中
X=xyz,U=uvw,
A 是正交矩阵,满足
U=AX,(25.19)
AT 为 A 的转置,曲面 Σ 是曲面 S 在正交变换 A 下的像。
讲解
公式 (25.18) 的证明如下:设 S 的参数方程为
S: x=x(s,t),y=y(s,t),z=z(s,t),(s,t)∈D.
通过 (25.19) 可得 Σ 的参数方程为
Σ: u=u(s,t),v=v(s,t),w=w(s,t),(s,t)∈D.
则
∬Sf(X)dS=∬Df(X(s,t))EG−F2dsdt,(25.20)
∬Σf(ATU)dΣ=∬Df(ATU(s,t))E1G1−F12dsdt,(25.21)
其中
E=Xs2,F=Xs⋅Xt,G=Xt2,
E1=Us2,F1=Us⋅Ut,G1=Ut2.
由 (25.19) 得
Us=AXs,Ut=AXt.
从而
E=E1,F=F1,G=G1.(25.22)
由 (25.20)—(25.22) 知 (25.18) 成立。
设 S 是球面 x2+y2+z2=1,f 是连续函数,证明:
∬Sf(ax+by+cz)dS=2π∫−11f(ua2+b2+c2)du,
其中 a,b,c 是常数。
证
不妨设 a2+b2+c2=0。设
uvw=Axyz,
其中矩阵 A 为正交矩阵,且 A 的第一行的元素为
a2+b2+c2a,a2+b2+c2b,a2+b2+c2c.
由 (25.18) 知
∬Sf(ax+by+cz)dS=∬u2+v2+w2=1f(ua2+b2+c2)dS=2∫−11∫−1−u21−u2f(ua2+b2+c2)1−u2−v21dudv=4∫−11f(ua2+b2+c2)du∫01−u21−u2−v2dv=4∫−11f(ua2+b2+c2)du∫0π/2dt(v=1−u2sint)=2π∫−11f(ua2+b2+c2)du.
2. 第二型曲面积分符号的确定(公式 (25.3) 的证明)
在计算第二型曲面积分时,要将第二型曲面积分化为二重积分,其中的困难之一是如何决定二重积分前的符号。在直角坐标系中计算时往往还比较容易,当用参数方程计算时就要麻烦一些。
讲解
下面是公式 (25.3) 的证明:首先由两类曲面积分之间的关系有
∬SPdydz+Qdzdx+Rdxdy=∬S(Pcosα+Qcosβ+Rcosγ)dS,
这里 (cosα,cosβ,cosγ) 为 S 所指定的一侧的单位法向量。我们知道 S 在参数方程下的法向量
n=±(A,B,C).
于是
cosαcosβcosγ=±A2+B2+C2AA2+B2+C2BA2+B2+C2C.
由此可以看出,其中 ± 号的选取应该是使 ±(A,B,C) 的方向与 S 所选的那一侧的法向量
(cosα,cosβ,cosγ) 的方向一致。另一方面,由
A2+B2+C2=EG−F2 得到
dS=A2+B2+C2dudv,
于是 (25.3) 式成立。
求第二类曲面积分
∬Sxdydz+ydzdx+zdxdy,
其中 S 是柱体 x2+y2≤1 被两个平面 z=0 和 z=4 所截部分的边界,积分沿边界的外侧。
解
设上底面为 S1,下底面为 S2,侧面为 S3,则
∬S1xdydz+ydzdx+zdxdy=∬S1zdxdy=4∬S1dxdy=4π.
在 S2 上,z=dz=0,于是在 S2 上的积分为零。下面计算在 S3 上的积分。令
x=cosθ,y=sinθ,z=z,
(θ,z)∈D={0≤θ≤2π, 0≤z≤4},
则
A=cosθ,B=sinθ,C=0.
于是
∬S3xdydz+ydzdx+zdxdy=∬D(cos2θ+sin2θ)dθdz=8π.
这里积分号前取正号是因为向量 (A,B,C) 与柱面外侧法方向一致,最后得到
∬Sxdydz+ydzdx+zdxdy=12π.
3. 球面上的曲面积分(证明关系式 (25.28))
若曲面 S 是以 M0=(x0,y0,z0) 为球心、r 为半径的球面
∂Br(M0),则其参数方程为
x=x0+rsinφcosθ,y=y0+rsinφsinθ,z=z0+rcosφ,
0≤θ≤2π,0≤φ≤π.
讲解
由此计算出
E=xθ2+yθ2+zθ2=r2sin2φ,
F=xθxφ+yθyφ+zθzφ=0,G=xφ2+yφ2+zφ2=r2.
于是
EG−F2=r2sinφ.
与直角坐标系中三重积分化为球坐标积分的 Jacobi 行列式一样,但这里的 r 是常数。于是
∬∂Br(M0)f(x,y,z)dSr=∫02π∫0πf(x0+rsinφcosθ, y0+rsinφsinθ,z0+rcosφ)r2sinφdθdφ.(25.23)
其中 dSr 是半径为 r 的球面的面积元。从上面的计算可以看出
dSr=r2sinφdθdφ,(25.24)
dS1=sinφdθdφ.(25.25)
于是就有
dSr=r2dS1.(25.26)
如果记
sinφcosθ=α1,sinφsinθ=α2,cosφ=α3,
注意这里 α1,α2,α3 都是 φ,θ 的函数,则 (25.23) 可改写为
∬∂Br(M0)f(x,y,z)dSr=∬∂Br(0)f(x0+rα1, y0+rα2, z0+rα3)dSr=r2∬∂B1(0)f(x0+rα1, y0+rα2, z0+rα3)dS1.(25.27)
另一方面,(α1,α2,α3) 恰是球面 ∂Br(0)(也是 ∂B1(0))上的单位法向量 n,方向指向球面的外部,因此上述表达式又可以简写为
∬∂Br(M0)f(x,y,z)dSr=∬∂Br(0)f(M0+rn)dSr=r2∬∂B1(0)f(M0+rn)dS1.(25.28)
设 Br(M0) 是以 M0=(x0,y0,z0) 为心,r 为半径的球,
∂Br(M0) 是以 M0=(x0,y0,z0) 为心,r 为半径的球面,证明:
(1)∭Br(M0)f(x,y,z)dxdydz=∫0r∬∂Bρ(M0)f(x,y,z)dSρdρ,
(2)drd∭Br(M0)f(x,y,z)dxdydz=∬∂Br(M0)f(x,y,z)dSr.
证
只证 (1),因为在 (1) 的两边对 r 求导便得到 (2)。作球坐标变换
x=x0+ρsinφcosθ,y=y0+ρsinφsinθ,z=z0+ρcosφ,
0≤θ≤2π,0≤φ≤π,0≤ρ≤r.
则
∭Br(M0)f(x,y,z)dxdydz=∫0r∫0π∫02πf(x0+ρsinφcosθ, y0+ρsinφsinθ,z0+ρcosφ)ρ2sinφdθdφdρ=∫0r∬∂Bρ(M0)f(x,y,z)dSρdρ.
最后一个等号是由 (25.23) 得出。
25.5.3 参考题
设 S 是平面 x+y+z=t 上被球面 x2+y2+z2=1 所割下的部分,
φ(x,y,z)={1−x2−y2−z2,0,x2+y2+z2≤1,x2+y2+z2>1.
证明:
∬Sφ(x,y,z)dS={18π(3−t2)2,0,∣t∣≤3,∣t∣>3.
解
作正交变换,使新坐标满足
u=3x+y+z.
则球面仍化为 u2+v2+w2=1,平面化为 u=t/3。当 ∣t∣>3 时该平面与单位球无交,积分自然为 0。当 ∣t∣≤3 时,截面 S 在 uvw 坐标中就是圆盘
v2+w2≤1−3t2,
且在该平面上
φ=1−u2−v2−w2=1−3t2−v2−w2.
因而
∬SφdS=∬v2+w2≤1−3t2(1−3t2−v2−w2)dvdw=2π∫01−t2/3(1−3t2−ρ2)ρdρ=18π(3−t2)2.
证明:
∫02πdx∫0πsinyesiny(cosx−sinx)dy=2π(e2−e−2).
解
把积分看成单位球面的第一类曲面积分。令
X=sinycosx,Y=sinysinx,Z=cosy,
则 X2+Y2+Z2=1,且 dS=sinydxdy,所以
I=∬X2+Y2+Z2=1eX−YdS.
应用例题 25.5.1,其中 (a,b,c)=(1,−1,0),得
I=2π∫−11e2udu=2π(e2−e−2).
因而题末常数应为 2π。
设 Ω 为空间第一卦限中的区域,函数 f(x,y,z) 在 Ω 上有连续一阶偏导数。S 为 Ω 中任一光滑闭曲面,试给出第二型曲面积分
∬Sf(x,y,z)(xdydz+ydzdx+zdxdy)=0
的充分必要条件,并证明之。
解
记
F=f(x,y,z)(x,y,z).
则题中积分就是 ∬SF⋅ndS。由 Gauss 公式,
∬SF⋅ndS=∭ΩdivFdV.
而
divF=∂x∂(xf)+∂y∂(yf)+∂z∂(zf)=xfx+yfy+zfz+3f.
因此该曲面积分对 Ω 中任意光滑闭曲面都为 0 的充分必要条件是
xfx+yfy+zfz+3f=0
在 Ω 内处处成立。
设 Σ 是光滑的闭曲面,围成的区域为 Ω,n 为 Σ 上单位外法向量,
(x0,y0,z0) 为 Ω 内固定一点,(x,y,z)∈Σ,
r=(x−x0,y−y0,z−z0),证明:
∬Σcos(n,r)dS=2∭Ω∣r∣dxdydz.
解
取向量场
F=∣r∣r=(∣r∣x−x0,∣r∣y−y0,∣r∣z−z0).
则
F⋅n=cos(n,r),
而直接计算得
divF=∣r∣2.
由 Gauss 公式,
∬Σcos(n,r)dS=∭ΩdivFdV=2∭Ω∣r∣dxdydz.
已给平面 Π:Ax+By+Cz=D,对于 Π 的任一定向,求
ω=P(x,y,z)dx+Q(x,y,z)dy+R(x,y,z)dz,使得沿 Π 上任意逐段光滑简单封闭曲线
Γ(Γ 的定向与 Π 的定向一致)恒有
∮Γω=S(Γ),
其中 S(Γ) 为 Γ 在 Π 上所围区域的面积。
解
设
N=A2+B2+C2,n=N(A,B,C)
是与题设定向一致的单位法向量。取
ω=21(n×r)⋅dr,
即
ω=2N1[(Bz−Cy)dx+(Cx−Az)dy+(Ay−Bx)dz].
由第 25.3.4 节第 1 题,
∮Γω=S(Γ).
设 Γ 是 R3 中逐段光滑简单封闭定向曲线,对于 (x,y,z)∈/Γ,定义
P(x,y,z)=∮Γr3(ζ−z)dη−(η−y)dζ,
Q(x,y,z)=∮Γr3(ξ−x)dζ−(ζ−z)dξ,
R(x,y,z)=∮Γr3(η−y)dξ−(ξ−x)dη,
其中 (ξ,η,ζ)∈Γ 是积分变元,
r=(ξ−x)2+(η−y)2+(ζ−z)2.
证明:
∂x∂Q=∂y∂P,∂y∂R=∂z∂Q,∂z∂P=∂x∂R.
解
令
ω=Pdx+Qdy+Rdz.
记
u=ξ−x,v=η−y,w=ζ−z,r=u2+v2+w2.
在对参数求导时可将偏导数移入积分号内,于是
∂y∂P−∂x∂Q=∮Γ[−∂v∂r3wdη−vdζ+∂u∂r3udζ−wdξ]=−∮Γd(r3w)=0.
同理
∂z∂Q−∂y∂R=−∮Γd(r3u)=0,
∂x∂R−∂z∂P=−∮Γd(r3v)=0.
故
∂x∂Q=∂y∂P,∂y∂R=∂z∂Q,∂z∂P=∂x∂R.
在上题中求函数 u(x,y,z),使得 du=Pdx+Qdy+Rdz。
解
由上题知 ω=Pdx+Qdy+Rdz 在 R3∖Γ 上是闭的。严格地说,除非再去掉一张以 Γ 为边界的定向光滑曲面,否则它一般没有全局单值原函数。若取这样的曲面 Σ,则在
R3∖Σ 的每个连通分支上可定义
u(M)=∫M0Mω,
这里积分路径限制在同一连通分支内。由于该分支中闭路积分为零,上式与路径无关,且显然满足
du=Pdx+Qdy+Rdz.
这个 u 就是曲线 Γ 对点 M 的有向立体角;跨过 Σ 时,两侧的函数值相差 4π。
利用 Gauss 公式证明 Archimedes 的流体静力学定律:物体在液体中所受的浮力等于物体排开液体的重量,方向垂直向上。
解
设浸没在液体中的物体表面为 Σ,所围体积为 Ω,液体密度为 ρ,z 轴竖直向上。静水压满足
∇p=(0,0,−ρg).
液体对物体的合力为
F=−∬ΣpndS.
对三个分量分别用 Gauss 公式,得
F=−∭Ω∇pdV=∭Ω(0,0,ρg)dV=(0,0,ρg∣Ω∣).
这正是排开液体的重量,方向垂直向上。
评论
支持 Markdown 和 LaTeX 数学公式。