格林公式、高斯公式、斯托克斯公式的应用【范本模板】

的应用

摘要

格林（Green）公式，斯托克斯(Stokes）公式和高斯（Gauss）公式是多元函数积分学中的三个基本公式，它们分别建立了曲线积分与二重积分、曲面积分与三重积分、曲线积分和曲面积分的联系。它们建立了向量的散度与通量、旋度与环量之间的关系,除了在数学上应用于计算多元函数积分，在其他领域也有很多重要的应用。本文将主要从这三个公式与物理学之间的联系展开介绍它们的其他应用，其中包括应用于GPS面积测量仪，确定外部扰动重力场，应用于保守场以及推证阿基米德定律和高斯定理等，帮助人们加深对格林公式、斯托克斯公式和高斯公式的理解，从而能够更准确地应用此三个公式.

关键词：格林公式 斯托克斯公式 高斯公式 散度 旋度 应用

目录

一、引言 ......................................... 1 二、格林（Green)公式的应用 ....................... 1

（一） 格林公式的定义 ............................. 1 1、单连通区域的概念 ............................ 1 2、区域的边界曲线的正向规定 .................... 1 3、陈述 ........................................ 1 （二）格林公式的物理原型 .......................... 1 1、物理原型 .................................... 2 2、 计算方法 ................................... 2 （三）格林公式与GPS面积测量仪 .................... 3 1.应用曲线积分计算平面区域面积 ................. 3 2。GPS面积测量仪的数学原理 ..................... 4 3。实验结果 .................................... 5 4.进一步讨论 ................................... 5 (四）应用格林积分直接以地面边值确定外部扰动重力场 . 6 1。扰动重力位的地面边值问题 .................... 6 2。地面边值问题的格林公式表示 .................. 6

三、Stokes公式的应用 ............................. 8

（一）Stokes公式简介 ............................. 8 （二）环量与环量密度 .............................. 9 (三）环量的应用 ................................... 9

1。开尔文定理 .................................. 9 2.开尔文定理的推论 ............................ 10 3。升力 ....................................... 10 （四)旋度 ........................................ 11 （五） 旋度的应用 ................................ 12 1。 平面矢量场的旋度 .......................... 12 2.环流量是区域S内有无漩涡的量度 .............. 12 3。旋度是矢量场某点漩涡强度的量度 ............. 13 4.空间矢量场的旋度 ............................ 14

四、Gauss公式的应用 ............................. 16

1、 数学中的高斯公式 .......................... 16 2、 保守场的推导 .............................. 17 3、 高斯公式在电场中的运用 .................... 17 4、 高斯定理在万有引力场中的应用 .............. 19 5.高斯公式推证阿基米德浮力定律 ................ 21 6。高斯公式推证静电场中的高斯定理 ............. 22 7。高斯公式与散度 ............................. 24五、结语 ........................................ 25 六、参考文献 .................................... 26

一、引言

格林(Green）公式，斯托克斯（Stokes）公和高斯（Gauss)公式是多元函数积分学中的三个基本公式,它们分别建立了曲线积分与二重积分、曲面积分与三重积分、曲线积分和曲面积分的联系。它们有很强的物理意义即建立了向量的散度与通量、旋度与环量之间的关系，因此它们有许多重要的应用,在数学上它们主要用来简化某些多元函数积分的运算,而在其他各个专业领域它们也有很多重要的应用。接下来将一一介绍它们在不同专业中的应用。

二、格林（Green）公式的应用

（一）格林公式的定义

Green公式反映了第二型平面线积分与二重积分的联系。

1、单连通区域的概念

设D为平面区域,如果D内任一闭曲线所围的部分区域都属于D，则D称为平面单连通区域；否则称为复连通区域。

通俗地讲，单连通区域是不含\"洞\"（包括”点洞\"）与\"裂缝\"的区域。

2、区域的边界曲线的正向规定

设L是平面区域D的边界曲线，规定L的正向为：当观察者沿的这个方向行走时，平面区域（也就是上面的D)内位于他附近的那一部分总在他的左边. 简言之：区域的边界曲线的正向应符合条件：人沿曲线走，区域在左边，人走的方向就是曲线的正向。 3、陈述

设闭区域D由分段光滑的曲线L围成,函数Px,y及Qx,y在D上具有一阶连续偏导数，则有 PQ (1)xyPdxQdy LD 其中L是D的取正向的边界曲线。公式⑴叫做格林（green）公式。格林公式沟通了二重积分与对坐标的曲线积分之间的联系,因此其应用十分地广泛

［1]. (二)格林公式的物理原型

1

在工科的“高等数学”教材中,格林公式这部分都是先给出定理,然后加以证明、应用。讲这部分内容时，总有学生询问同一问题,即人们怎样想到这个公式，怎样想到曲线积分与重积分会有这样的数值上的联系?能否将格林公式的来源即物理原型加入教材呢?在教学中,试着加入这部分内容，并对公式作了简单的符号记法，简化了公式,降底了出错率,并对应用总结了几个类型.多年的实践证明，效果是很好的，下面就将加入的内容介绍如下[2］:

1、物理原型

在流体物理学中,称满足下述三个条件的“流速场”为“平面稳定流动”。 （1）场中每一点的速度都不随时间改变,只是位t的函数即

VP(x,y)iQ(x,y)j

(2）所论流体介于两个互相平行的平面之间(为方便，不妨设平面间距离为l 个单位)其中之一称为底面(往往底面即为xoy坐标面）。 （3）垂宜于底面的直线上的各点流速相等, 并平行于底面。

在这种“ 平面稳定流动” 中,我们来计算单位时间内流过曲线C的流体体积即流t 密度（ 其实是流过以C 为准线、高为l 的柱体的流体体积； 简单用面积表示) 其中C 是平面上一个闭的、无重点, 光滑曲线。无重点， 是指曲线

X(t),Y(t)，当t1t2时，点((t1),(t1))与((t2),(t2))总是相异的。

2、计算方法

（1）在C上任取一小段弧线△S,在△t时间内流过△S的流体面积,近似于一个

平行四边形的面积，它的一个边长是△S另一个相邻边长是流程Vt

因此面积为sVtcos(Vn)s(Vn)t

其中n是C的单位法向量

单位时间内流体面积为：Vns

ncos,cos 由曲线积分定义有总的流体面Vnds;l为C的全长，设l0则PcosQcosds

0l设l为点（x,y)处的切线，与x轴夹角cossin,coscos

PsinQcosdsPdyQdx

0l（2）的计算可以从另一个角度来计算,那就是先算出流过场内每一个微dxdy

2

在单位时间内散发出去的流体的面积,然后求其总和。 设上述曲线C所围平面区城为G，在G内任取一个微元dxdy

显然在单位时间内从左边流进（x轴方向）这个微元的流体面积近似于Pdy ,而从右边流出的面积近似于PP'xdxdy,（P'xdx为偏增量的近似)。因此这个微元在单位时间内沿x方向（净）散发出去流体面积近似于

PP'xdxdyPdyP'xdxdy。同理沿y方向（净)散出去的流体面积近似于

Q'ydxdy，所以总的和为P'xQ'ydxdy

PQ由重积分的定义得：dxdy cxyPQ有（1)、（2）可得：PdyQdxdxdy cxy这是场论中最根本的公式，即格林公式的原型。

(三）格林公式与GPS面积测量仪

格林公式作为多元微积分中联系平面曲线积分与二重积分的一个重要公式,不仅给出了一个有效计算平面曲线积分的方法，而且给出了一种已知边界曲线方程的平面区域面积的计算方法．在这部分的教学内容中，传统应用主要局限于纯几何与物理问题的解决，很少应用于生活实际问题的讨论．本文在基于微元法的基础上，讨论了GPS面积测量仪测量平面区域面积的数学原理，并在教学实践中，将其以引入性问题和课程探索性实验的形式作为曲线积分教学内容的扩充,实现了抽象的数学理论与方法和生活实际的有效结合[3]．

1。应用曲线积分计算平面区域面积

设D为xOy平面上的闭区域,其边界∂D由光滑或分段光滑闭曲线组成，函数P(x，y）和Q(x,y）在D上有连续的一阶偏导数，则有

LP(x,y)dxQ(x,y)dy(DQP)dxdy（1） xy

其中D的方向为关于区域D的正方向．曲线正方向的确定使用“左手法则”，即当一个人沿着该方向行走的时候，区域位于左手一侧．式（1)对于平面单连通区域或多连通区域都成立．

当式（1)中的二重积分被积函数为常数时,可以使用左端关于坐标的曲线积

3

分计算封闭曲线围成的平面区域的面积．即若则有

QPA xySD1ALP(x,y)dxQ(x,y)dy

（2）

因此，只要构造合适的P（x, y)和Q(x，y)，就可以通过封闭曲线∂D上的第二类曲线积分计算其围成的平面区域D的面积．则

SD

12LydxxdyLydxLxdy

(3）

2。GPS面积测量仪的数学原理

利用格林公式或二重积分方法计算平面区域的面积时，一般需要知道其边界曲线方程，而在实际生活中，这样的边界曲线方程是很难知道的，因此无法直接使用它们来完成对面积的精确计算．GPS面积测量仪则给出了比较好的平面区域面积的近似计算方法．只要手持测量仪绕行测量区域一周．仪器就可以通过自动记录行进路线的坐标，计算所围绕区域的近似面积．

设由边界曲线3D围成的区域和使用GPS测量仪记录的平面坐标为

P,2...,n) i(xi,yi)(i1

图1 目标区域与记录点位置

由式(2)可知，在闭曲线方程已知的情况下,对其围成的封闭区域面积的计算可以转换为曲线积分计算．假设闭曲线方程未知，则根据积分的存在性,借助于微元法思想，封闭曲线可以近似为由有向线段 的并，其中 其中Pn1P1，即

(4) DPP12P2P3...PnP1

从而有

4

1n（5) SDP(xi,yi)(xi1xi)Q(xi,yi)(yi1yi)Ai1

其中xn1x1,yn1y1。

3。实验结果 下面以参数方程

x=4sint—sin4t （6） y=4cost-cos4t

确定的封闭曲线为例，在Mathematica中进行数学实验验证。

由于该封闭曲线方程已知，所以由公式（3)，利用第二类曲线积分的直接计算方法,可得所围封闭区域面积为20π≈62．832．取参数增量分别为

t6,12,24,48,480,依次在曲线上取点，计算得到的结果分别为

53．196，60．086,62．122，62．653，62．830．若取P(x， y）=—y，Q(x， y)=0,或者P（x， y）=0，Q（x, y)=x虽然在近似计算形式上看似有所差别,但是在

Mathematica中以默认精度进行计算时，每个结果可以保持在小数点后13位一直相同，并且随着分割的细化，结果逼近直接计算得到的精确结果。

4。进一步讨论

使用边界点坐标方法计算区域的面积还有借助于微元法思想和辛普森公式容易验证的公式．对任一个平面凸区域D(即过该区域能做一组与区域边界曲线交点不多于两个的平行直线的区域），设正好夹住平面区域的两平行直线的距离为b．在两平行直线之间做n-2（偶数)条距离为b/n，平行于这两条直线的一组直线，各条直线夹在闭曲线∂D围成的区域D范围内的线段长度记作 （i=1，2，⋯,n—1)。

图2 平面凸区域面积近似方法

5

通过坐标系旋转或者存在有一组平行于Y轴的直线，b即为区域在z轴上投影区间的长度，这样实际上也就是由微元法构造定积分模型的形式．该方法思想简单，在实际计算中相对来说约束较多。除了以上借助于曲线上点坐标近似计算平面区域面积之外，另外也可以通过已知点列坐标,利用插值、拟合的方法获取近似边界分段曲线方程,然后利用二重积分或者第二类曲线积分计算面积．同时,这种近似计算的思想也适用于求曲线的弧长,比如椭圆周长的近似计算与一些不可积函数的积分计算.

(四)应用格林积分直接以地面边值确定外部扰动重力场

1.扰动重力位的地面边值问题

确定地球外部重力场和大地水准面是大地测量学的主要任务之一。确定地球外部重力场和大地水准面的斯托克斯理论需要将地面观测的重力异常归算至大地水准面,再采用调和函数球面边值的解式(如Stokes 公式)求得大地水准面及其外部的扰动重力位。归算将涉及对大地水准面.至地面的质量迁移，对大地水准面产生间接影响，而且由于归算对质量进行了调整，改变了外部扰动重力场，因此归算到大地水准面上的重力异常用以确定外部扰动重力位会导致结果的歪曲。直接以地面重力异常为边值的Molodensky问题从理论上避免了归算的困难，成为近代外部重力场研究的理论基石。然而，由于地球表面的复杂性，给这一问题的求解带来极大难度。

Molodensky 基于基本积分方程的小参数解法得到地面扰动位的级数解式。。提出将地面重力异常解析地延拓到一点的水准面上，再采用球面的Stokes 积分得到地面扰动位，其结果同样是级数的形式。也研究得到类似的级数解..则提出将地面重力异常调和地延拓到一个内部球面上,再由球面边值问题解逼近外部扰动位，其调和延拓需要求解Poisson积分方程.尽管这些理论解的途径有所不同，但在一定前提下它们是等价的。虽然经过线性化和地球椭球作球近似后的所谓简单Molodensky 问题的研究已得到几近完美的理论结果，但它们的实现仍具有相当大的困难.由于受到数据和高阶项计算稳定性的限制，目前在实际上通常只能考虑到一阶项.对于确定地球外部扰动重力场问题,上述解在应用上受到一定的限制。像Molodenky 解通常应用于地面，Moritz 的解析延拓解和Bjeharmmar 解虽然可以拓展到外部空中，但边值的延拓仍是一个较复杂的过程。本文侧重于应用的需要，讨论直接由地面边值确定外部扰动位的方法。

2.地面边值问题的格林公式表示

确定地球外部扰动重力位T 归结为下面的边值问题。

T0 在地面的外部

6

[4]

BT|f 在地面上



222式中222为Laplace 算子，B 代表某一泛函算子，f为已知泛函.

xyz由位理论知，T作为调和函数可以由格林第三恒等式表示为

Tp141T1[Tlnn(l)]d.式中，l是计算点p至地面Σ上面元dΣ的空

间距离。n是相对于调和空间的边界面外法线方向。取局部北东天坐标系，求法线方向导数得

TTxTyTz nxnynznTTT,,)(x,y,z)为扰动重力矢量xyz根据位理论((xyzT,,)(cos,cos,cos)为法线的方向余弦。，因此可知n即为扰nnnn动重力在法线方向的分量,如图所示

图3 内法线方向示意图

对于所谓的简单Molodensky 问题，即将地球椭球近似为球面时,上述各元

2素的几何关系见下图。由距离公式l2rpr22rprcos式中l为P 点与dΣ 单

元处的距离；rp为P点的

球心距离；r为dΣ单元处的球心距离;ψ为P点到dΣ单元处的极距求导可知，

1112cos(l,n)，因此可得Tpnll4

1Tcos(l,n)[T]d。 lnl27

图4 球近似下各元素的几何关系

应用格林公式可以在实际地球表面上计算外部扰动位.其条件是需要同时具有地面上的扰动位和扰动重力矢量的观测值。这在实际应用中是有困难的。一方面,所需的边值条件很难满足.另一方面地球表面非常复杂，这就使得在地表起伏较大地区该式中的法线方向变化剧烈，其计算相当困难.尽管如此，格林公式提供了不需作任何边值的归算或延拓而以地球自然表面上的边值条件确定外部扰动重力场的唯一可能的解析形式。

三、Stokes公式的应用

（一）Stokes公式简介

Green公式给出了平面上沿闭曲线(C)的第二型线积分与(C）所围成平面区域上二重积分之间的关系.现在把它推广到空间，考察沿空间闭曲线(C)的第二型线积分与（C)上所张曲面的面积分之间的关系.

设区域(G)R3,P,Q,RC(1)((G)),（C）为(G)内一条分段光滑的有向简单闭合曲线，(S）是以（C）为边界且完全位于(G)内的任一分片光滑的有向曲面,（C)的方向与（S）的法向量符合右手螺旋法则，则

（C）PdxQdyRdz((S)RQPRQP)dydz()dzdx()dxdy yzzxxy称为Stokes公式。

设A=（P,Q,R），根据nabla算子的定义，Stokes公式可写成向量形式：

(C)Ads(A)dS(A)edS

n(S)(S)8

如果A为一平面场（P,Q）,（C)为一平面闭曲线，（C）所围成区域为()，这时，Stokes公式1就蜕化为Green公式，可见,Stokes公式是Green公式的推广[1].

(二）环量与环量密度

类似于平面向量场沿平面闭曲线(c)的环量，空间向量场沿空间闭曲线(c)的线积分

(C)A(x,y,z)dsP(x,y,z)dxQ(x,y,z)dyR(x,y,z)dz

(C)称为A(x，y,z)沿闭曲线(c)的环量，它同样表示了A绕(c)旋转趋势的大小。 以n为法向量，过点M作任一微小曲面(S)，它的边界曲线记为(C)。并选取(C)的正向使与n复合右手螺旋法则.当(S)很小时,A沿(C)的环量与小曲面(S)的面积之比，

ΔΓ1Ads ΔSΔS(c)近似的反映出A点在M点附近绕方向n的旋转趋势大小.让小曲面（）在保持n为其法向量的前提下任意缩向点M，若上述比值的极限存在,则称此极限值d为A在M点沿n方向的环量密度，记作，即

dSd1limlimds cAdS(Δs)MS(Δs)MS()（三）环量的应用

1。开尔文定理

流体动力学中的一个著名的定理。内容是：在无粘性、正压流体中（见正压流体），若外力有势，则沿由相同流体质点组成的封闭曲线的速度环量在随体运动过程中恒不变。

在流体力学中，沿封闭曲线的速度环量定义为线积分：

vdr

L式中Γ 为速度环量；v为速度矢量；dr为封闭曲线L的线段元矢量.速度环

9

量和涡通量（见涡旋)通过下列斯托克斯公式联系起来：

vdrQdS

Ls式中S是张在封闭曲线L上的曲面；Ω和dS分别为涡旋矢量和面积元矢量。

2.开尔文定理的推论

由开尔文定理可推出反映涡旋保持性的涡旋不生不灭定理：假设流体是无粘性和正压的，且外力有势，若初始时刻在某部分流体内无旋，则在此时刻以前或以后的任一时刻中，这部分流体皆无旋。反之，若初始时刻该部分流体有旋，则在以前或以后的任一时刻，这一部分流体皆有旋。因为若初始时刻某区域内的流体运动无旋,则根据斯托克斯公式（2),该区域内沿任一封闭曲线的速度环量为零。设过一时刻此区域内的流体运动到一新区域,从开尔文定理易见,在新区域内沿任一可能的封闭曲线的速度环量也为零.换言之，线积分与积分路径无关,它只是时间t以及变动点B的坐标r和固定点A的坐标r0的标量函数，可记为

故v，即存在速度势Φ(r,t)。由v()0，推出整个流动是无旋的。

vdr(r,t)(r,t)

A0B对于在重力场作用下的无粘性不可压缩均质流体,考察均匀来流定常绕流和从静止起动的流体运动。显然，两种情形都满足流体无粘性、正压和外力有势三个条件。流场中任一流体质点都来自无穷远处或初始的静止流体.因无穷远处均匀来流和静止流体都是无旋的,根据涡旋的不生不灭可以看出，整个流场都是无旋的。由此得到开尔文定理的一个重要推论：对于在工程实际中大量遇到的无粘性不可压缩均质流体在重力作用下的均匀来流定常绕流问题和静止起动问题，整个流体运动时时处处都是无旋的。由于无旋运动有些特殊性质，处理这类流动可作许多数学上的简化（见拉普拉斯无旋运动）. 3.升力

升力， 也就是向上的力大于向下的力，其合力可以使物体上升. 这个力就是升力。升力的成因较复杂，因为要考虑实际流体的粘性、可压缩性等诸多条件.目前大多用的是库塔儒可夫斯基定理,它是工程师计算飞机升力最精确的方法。具体内容就是由绕翼环流导致升力，产生了上下压力差，这个压力差就是升力 (Y）,升力和向后的诱导阻力（d）合成为空气动力（R）.流过各个剖面升力总合就是机翼的升力.升力维持飞机在空中飞行. （1）升力的来源

升力来源于机翼上下表面气流的速度差导致的气压差.但机翼上下表面速度差的成因解释较为复杂，通常科普用的等时间论和流体连续性理论均不能完整解释速度差的成因。航空界常用二维机翼理论，主要依靠库塔条件、绕翼环量、库

10

塔—茹可夫斯基定理和伯努利定理来解释。 (2）库塔条件

在真实且可产生升力的机翼中，气流总是在后缘处交汇，否则在机翼后缘将会产生一个气流速度很大的点.这一条件被称为库塔条件，只有满足该条件，机翼才可能产生升力.

（3）库塔如茹可夫斯基方程式

由满足库塔条件所产生的绕翼环量导致了机翼上表面气流向后加速，由伯努利定理可推导出压力差并计算出升力，这一环量最终产生的升力大小亦可由库塔-茹可夫斯基方程计算（适用于不可压缩流体）：

物体单位长度上所受到的升力：

L（升力）v（气体密度流速环量值）

其中环量是流体的速度沿着一条闭曲线的路径积分.如果v是流体的速度,ds是沿着闭曲线C的单位向量，那么:

环量的量纲是长度的平方除以时间。这一方程同样可以计算马格努斯效应的

CVds气动力.不过以上理论仅适用于亚音速（更准确地说是Ma小于0.3），在超声速飞行时由于空气是可压缩的,伯努利定理不成立，此时无环流运动，升力主要靠机翼上下表面的激波所导致的压力差.当飞机以一定迎角在超声速流中飞行时上表面前端处与来流成一个凸面，形成膨胀波，气流流过膨胀波时压力下降，而下表面与来流形成一个凹面，导致激波，气流流过激波时压力增加。因此上表面压强小,下表面压强大，产生升力。

（四）旋度

若在场A(M）中一点M处存在这样一个向量，其方向为使A在点M环量密度最大的方向，其模等于环量密度的最大值，则称此向量为场A（M）在点M的旋度，记做rotA。则

dRQPRQP()cos()cos()cos dSyzzxxy为旋度的计算公式。

利用旋度，还可将Stokes公式写成下列形式

（C）AdsrotAdS（S）

11

（五）旋度的应用

1.平面矢量场的旋度

旋度最早是通过研究水流的涡旋建立起来的概念［5］。河水流动时，由于水有内摩擦力，因而靠两岸速度较小，河中间速度较大。故漂在水面上的救生圈一边顺流而下，一边还会旋转，这说明水中有涡旋，如下图所示.

图5 速度分布和涡旋特征

2.环流量是区域S内有无漩涡的量度

在平面流速场V(x,y)中作有向封闭曲线L，则流速场V沿L的环流(图3)

bcdaVdlVdlVdlVdlVdl

Labcd

V2dl0dl0acbd

V2abV1cd 0

在均匀流速场V(x,y)中，由于v1v2，所以V沿L的环流

12

Vdl0

L环流不等于零，在区域S内无涡旋。由特殊到一般，对于任一平面矢量场

A(x,y)如果

Adl0

L说明在区域S内有涡旋;如果

Adl0

L说明在区域S内无涡旋。因此环流是平面矢量场A在区域S内有无涡旋的量度.

3.旋度是矢量场某点漩涡强度的量度

环流的大小与封闭曲线L所包围的面积△s有关，所以不能用环流的大小

1Vdl来量来量度涡旋的强弱.而用环流与面积△S之比，即平均涡旋强度

SL1Vdl度△S区域内的涡旋强度.当S0，且收缩到P点时，用极限limS0SL来量度p点处的涡旋强度。此极限称为平面流速场V在p点的旋度，用V表示，即

1VlimVdlS0SL

 可见，旋度是环流对面积的变化率。

特殊到一般，任一平面矢量场直在p点的旋度

1AlimAdlS0SL

在直角坐标系中，平面矢量场A(x,y)在P(x,y)点的旋度

AxA（x,y）()k

xy13

Ay4。空间矢量场的旋度

例1 水池中的水漏掉时，会形成涡旋，如下图所示。

图6 水池中的涡旋

以p点为回心，作两个圆周L1和L2，两圆周面积相等,均为S,它们的法线n1,n2的方向与L1和L2的绕向符合右手螺旋法则.显然Vdl＞VdlL1L2

除以S,得

11Vdl＞VdlSL1SL2

可见，同一点p绕n1方向的平均涡旋强度大于绕n2方向的平均涡旋强度。为了量度场中任意一点的涡旋强度，必须把S取得很小,同时为了能比较不同点的涡旋强度，应使L的空间取向能得瓢最大的环流.在图6中，流速场V沿Ll的环流最大。现在，可以得到任一空间矢量场A旋度的定义：矢量场A在p点的旋度A是一个矢量，其大小为当面积S趋于零时单位面积上A的最大环流,其方向为当面积S的取向使得环流为最大时该面积S的法线方向（法线方向的单位矢量n，的方向与不的绕向符合右手螺旋法则)，即

1Alim(nAdl)

LS0S14

（x,y,z）（x,y,z）在直角坐标系中，矢量场A在P点的旋度

iA（x,y,z）xAxjyAykzAz

AzAyAxAzAyAx()i()i()kyzzxxy

例2绕定轴Z转动的刚体的角速度为。如图.求刚体上任一点P的线速度V的

旋度.

图7 绕定轴Z转动的刚体

解:点p的位置用位置矢量r来确定角速度

rxiyjzkk由力学知，点P的线速度



iVr0jk0yixjyz

xV的旋度

15

ijVxyyxk2k2 z0可见，速度场V的旋度V与刚体的旋转角速度之间有着密切的联系.

四、Gauss公式的应用

通过数学上高斯公式和保守场的定义，推导出物理学中的两个保守场:电场和万有引力场，把数学中的高斯公式运用到两个保守场中，结合类比思想，分别得出电场中的高斯定理和万有引力场中的高斯定理，并分别举例说明高斯定理在这两个场中的运用,万有引力场如同电场一样,是有源场，对于具有高对称性的质量西,万有引力高斯定理可简化万有引力场的求解过程。

数学是创立和发展物理学立论的主要工具，他的高度抽象性，使他能概括物理世界的基本结构[6］。数学中的高斯公式试分析大学物理中矢量场问题的重要工具[7］,电场中的高斯定理被广泛应用到求解具有高对称性的带电系周围的电场问题中[8]，而与电场同属于保守场的万有引力场中也能运用高斯定理,本文将把高斯公式推广到两个保守场中，得出各自场中的高斯定理，并举例说明其运用。

1、数学中的高斯公式

高斯公式是高等数学中曲面积分的一个重要公式，它可以把高斯面上的第二型曲面积分转化为所围体的三重积分[1］，描述为：

设空间区域V的边界曲面S是光滑的或分片光滑的，函数P（x,y，z）, Q(x，y，z), R（x，y,z)，在Q及S上具有已接连续篇倒数，S的方向为外发想，则

（

VsPQR)dxdydzxyz

PdydzQdzdxRdxdy （1）或

16

PQR（)dxdydzxyzV(PcosQcosRcos)dS (2)S

如果引入矢量函数,高斯公式又可以写为：

(3) dSadV SV数学意义为：矢量场通过闭合曲面S的流量等于此闭合曲面所围体V上每一点的在体V上的三重积分.

2、保守场的推导

数学中，保守场的电影以为:在矢量场中a汇总，若曲线积分与路径无关,只与起点和重点有关,这种常称为保守场。矢量场a为保守场的充要条件是： 0

大学物理中，电场力做功和万有引力做功都与路径无关，即：

(4) Edl0和Fdl0 LL根据斯托克斯定理得出：

(5) Adl(A)ds LS E

3、高斯公式在电场中的运用

0和F0 (6)

由此从数学推导,可以知道电场和万有引力场都是保守场［9］.

把高斯公式运用到大学物理的电场中，就是电场高斯定理.电场通量等于闭合曲面所围面V上每一点的散度在V上的三重积分（即

10dv,为电荷密度）

它可以表示数位：通过一个人已闭合曲面S的电通量r，等于该面所包围的所有电量的代数和

q除以0与闭合面外的电荷无关.即:

17

向四周辐射电场.

SEdS10qii内 （7）

在求解具有高对称性带电体系的电场分布时，电场高斯定理可以大大简化计算过程［5],还反映出电场的另一重要性质：静电场是有源场。正电荷是静电场的源头， 例:求一个均匀带正电球体内外的电场强度分布，设球体带电总量为Q，球体半径为R,如图

图8 带电量为Q的球体（虚线为高斯面）

如果用电场强度叠加原理,把带点体系分割成无穷个电荷元，每个电荷元视为点电荷，对无穷个电荷元产生的电场积分求和，势必牵扯到一个三重积分，仙人，这在数学计算中比较复杂。

如果应用高斯定理，由于电荷Q均匀分布在球面上，其电场分布具有球队成型,在任何与带电球面同心的球壳上的个点，电场强度大小均应相等，方向沿各自的矢径方向。

在球面外任取一点P,设想过P点做一个半径为r（r〉R）的球面，称之为高斯面,因球面上个点的法线方向与场强方向一致，所以通过该球面的电通量为:



此时该求面包为的电荷为:

根据高斯定理可得：

18

SEcosdSEdSE4πr2 （8）

Sqii内Q （9）

πr2 E4P点的电场强度为：

EQ0 （10）

1Q(rR) （ 11） 24π0r对于球内任一点p’，同样过p’作一半径为r(r〈R）的球形高斯面,通过该球面的电通量为：

r此时，该球面包围的电荷为：

2EcosdSEdSE4πr （12） ssQ43qr3πr3 (13) qi内34πR3Ri3由高斯定理可得：

Q3πrr （14） E4 30R2此时必有:

E

4、高斯定理在万有引力场中的应用

万有引力场和静电场同属于保守场，可以设想万有引力场如同电场一样是一种物理客观存在.质量为M的物体想四周辐射万有引力场,切万有引力定律运用类FGM1M2/r2与库仑定律 FQ1Q2/40r2相似,都服从平方反比定律。

（fEgdS)等于比思想,把高斯公式应用到万有引力场中，万有引力场通量

s1Qr （15）34π0R

闭合曲面所围体V上每一点的散度在体V上的三重积分.它可以表述为：通过一个任意闭合曲面S的万有引力场通量f等于该面所包围的所有质量的代数和

m 除以a，负号表明万有引力场为汇聚场，穿入闭合曲面S.

iii11EdS﹣m(a) （ 16） sg i内0a0i4G 电场高斯定理在求解具有高对称性带点体系的电场分布时有着十分重要的

19

地位,同理在求解具有高对称性质量体系的万有引力场分布时,用万有引力场中的高斯定理会大大简化求解过程.

例：求均质细棒中垂面上的收到的引力场强度，设棒长为2l，如图：

图9 质量均匀的细棒（圆柱体为高斯面）

如果用万有引力定律计算，选细棒中点O为原点，取坐标轴z沿细棒向上，由于细棒具有轴对称性，取纸平面作代表，细棒的中垂面与纸面的交线为中垂线OP。

整个细棒可以分割成一对一对的线元，其中每对线元dz和dz对于中垂线OP对称，这一对线质元在中垂线上任一点,从而合成P所产生的元引力场dEg和



dEg也对中垂线对称。他们在垂直于OP方向的分量互相抵消，从而合成矢量dEg

+dEg沿中垂线方向其大小为2dEkcos,其中：

dEg﹣G cosdzr2z2

rrzl22 (17)

r表示距离OP，表示线质量密度，则 dz是总质量.

2G Eg2dEcos﹣00lrdzr2z322﹣2Glrrl22 （18）

当细棒为无限长时,l这时的万有引力场强为:

2G （19） Eg﹣

20

1r 如果应用万有引力场中高斯定理，因该体系具有轴对称型，若以细棒为轴,在垂直于轴的平面内，同一圆周上的万有引力场大小处处香等,方向垂直于棒，向内汇聚,于是可选取以棒为轴，在垂直于轴的平面内，同一圆周上的万有引力场大小处处香等，方向垂直于棒，向内汇聚。于是可选取以棒为轴,半径为r长为l的封闭圆柱面为高斯面。圆柱的上、下两底面，Eg与dS的方向垂直；侧面的Eg与dS的方向相反,所以通过封闭圆柱体的万有引力通量为：

fEsgdS前面EgdS上下底面EgdS侧面EgdS0Eg2rl

(20)

根据高斯定理得：

Eg2rl﹣a10mi内﹣4Gl （21）

i得:

2G (22) Eg﹣公式（19）与公式（20）是相等的，可以看出，比起用万有引力定律来说，用万有引力场中的高斯定理求解问题比较简单，省去了很多积分过程，并且体现出万有引力场是保守场的性质。

5。高斯公式推证阿基米德浮力定律

在普通物理的教科书中，一般对阿基米德浮力定律都不作严格的数学证明，仅对它作一个说明。下面我们根据重力场中静止流体的压强分布，应用高斯公式给出一个证明.

一物体浮在液面上,液体表面的平面把浮体表面的封闭曲面S分为两部分

1rS1和S2，也把整个浮体分为两部分.其中浮在液面上的那部分为V1，浸没在液体中的那部分为V2。建立坐标系，取液体表面为x o y平面，Z轴的方向取为竖直向下。作用在曲面S1上的压强就是大气压P0，而作用在曲面S2上的压强则

为 PP0gz

式中P为液体的密度，z为曲面S2上某点处位于液面下的深度.作用在物

21

体上的浮力就是由于作用在物体下部的压强大于作用在物体上部的压强而产生的，我们来具体计算一下. 反，所以有：

F浮PdsPdsnPcosicosjxcoskdss因为作用在物体表面上任一面元上的压力总是与面元的法向矢量n方向相

siPdscosjPdscoskPdscosi(pi)dsj(pj)dsk(pk)dssssssss

（24)

式中为n与三个坐标轴的夹角，应用在高斯公式，上式可化为体积分：

F浮i(pi)dvj(pj)dvk(pk)dvvvvidvjdvkdvxyzvvv(ijk)dvxyzvP0dvv

v2PdvPdvv1P0dv(P0gz)dv0gkdvv2gkv2v1v2 （25）

上式即为我们所熟知的阿基米德浮力定律的数学表达式,它表 明：浸在液体里的物体受到向上的浮力，浮力的大小等于物体所排开的液体的重量。

6。高斯公式推证静电场中的高斯定理

1静电场中的高斯定理:eEdss0(Sq

i内)在一般的普通物理教科书中,对高斯定理都不作严格的证明，而是利用电力线的概念加以说明，也有少数书中是采用引进“立体角”的概念来证明的.其

22

实，应用高斯公式可以很简洁地证明高斯定理，而且不需要引进“立体角”的概念。

我们先讨论单个点电荷的情况。在闭合曲面S内有一点电荷q，它所产生

qrdsqr的场强为E，则穿出S的电通量为eEds 3340sr40rs 注意,现在我们还不能利用高斯公式把上式中的面积分化为体积分。 因为

高斯公式成立的条件要求E在闭合曲面S所包围的区域V内是连续函数。但

qr显然E3在区域V内r=0处（即点电荷q处）不连续。我们可在闭

40r曲面S内作一个以点电荷q为中心，以R为半径 的小球面S'的方向和S一样都

取为外法线方向.那么，在闭合曲面 S 和S之间的区域中，r0，E满足

'连续的条件,于是我们可利用高斯公式来计算通过区域的边界曲面的电通量。

qr eEdsEdv3dv40r （26）

r由直接计算可得，30 (r0)，所以e0.

r在面积分Eds，曲面的外法线方向在曲面S处与S的外法线方向

相同,而在S'曲面处则与S'的外法线方向相反。于是有

EdsEdse0 Edsss'所以EdsEds

ss'这说明通过包含点电荷的任意闭合曲面的电通量都与通过以 该点电荷为中心的任一球面的电通量相等。而通过球面的电通量很容 易算出：

qrqqq2'edsds4R 32240r40Rs'40R0s'所以eq0

当点电荷q不在闭合曲面S内时,则在S所包含的V区域内，r0，可直接利用高斯公式求出：

23

qeEdsEdv40svr3dv0 rv对于由一组点电荷q1q2qn所组成的带电体系来说，它们在空 间所产

生的总场强E是各点电荷单独存在时所产生的场强的迭加：

EE1E2En (27）

那么通过任意闭合曲面的电通量为:

EdsE1dsE2dsEndsssss

e1e2en式中e1、e2en是各个点电荷的电场通过闭合曲面S的电通量。由上述关于单个点电荷的结论可知： 当gi在S内时 eiqi0

1当gi在S外时ei0所以eei,0这样就证明了高斯定理.

7.高斯公式与散度

(S内)qi

设空间闭区域是由分片光滑的闭曲面所围成,若函数P（ x，y，z），Q（ x，y，z)与R（ x，y，z） 在上具有一阶连续偏导数,则有高斯公式（散度公式)：

(PQQ)dvxyzPdydzQdzdxRdxdy（是的整个边界曲面的外

侧）

在日常生活中，我们经常见到如图8用榔头钉钉子，图9灯泡或太阳向四周辐射光线,图10点燃的烟花向周围爆炸等现象.对这些现象进行对比观察，发现都具向四周散射的效果［10］。

我们不妨对钉子受力会被钉进桌面的现象进行受力分析,如图11，钉子受到榔头的力A=（P，Q,R） 作用，P，Q，R 分别是力A 在x,y，z 三个坐标轴方向上的分力，力A 的作用效果（变化率）等效于分力P，Q，R 在坐标轴方向

24

上的作用效果（变化率），从而我们可用度量A的（散射）作用效果，并称之为散度.因此，我们可以把高斯公式理解成:向量场（包含力场）A=（P，Q，R) 通过闭曲面流向外侧的通量（流量或辐射量）PdydzQdzdxRdxdy，等于A

的散度

PQR在闭曲面所围闭区域上的积分。 xyz

图10 用榔头钉钉子 图11 辐射光线

图12 点燃的烟花爆炸 图13 钉子的受力分析

五、结语

本文不仅详细剖析了多元函数积分学中三个基本公式，即格林（Green)公式，斯托克斯(Stokes）公式,高斯（Gauss)公式之间的关系，而且着重讨论了以上三个公式在在不同专业，尤其是在物理学上的应用。我们总结了诸如斯托克斯公式在环量、旋度上的原理和应用，高斯公式在保守场、散度的应用和阿基米德浮力定律的证明,以及格林公式应用于外部扰动动力场的计算和其在GPS面积测量仪中的应用。从中可以看出，这三大公式不仅在数学上应用广泛，而且还被用于跨学科，跨专业解决问题和公式定理的证明，并使之应用于实际中。以小见大，数

25

学学科不仅仅与物理学科关系密切，相辅相成，同时也是解决诸多实际问题的强力武器.

六、参考文献

[1］ 马知恩，王绵森,工科数学分析下册，北京：高等教育出版社 2006年 ［2] 包海臣，格林公式的物理原型及其它，呼和浩特职业学院 2005年

［3］ 刘雄伟，汪雄良，格林公式与GPS面积测量仪，高等数学研究，第18卷第一期：109-111页 2015年

[4] 田家磊，吴晓平，李姗姗，应用格林积分直接以地面边值确定外部扰动重力场，测绘学报，第44卷11期:1189-1195页 2015年

［5] 黄国良,矢量场散度和旋度的物理意义，西安矿业学院学报，第一期：71—77页 2003年

［6］陈修芳．从高斯公式到高斯定理［ J ]。科学之友。2015.5:11—1 3. [7］籍延坤．高斯定理的数学证明[ J ］．大连铁道学院学报．2004.25（3)13-16。 [8]温耐，王伟峰．高斯定理在静电场中的应用问题L J．物理通报，2010 （ 11 )：9.11． ［9］宋克慧，龙文海．关于保守场判据的一点注记[ J ］．蒙自师专学报：自然科学版．1993.1 0 ( 2 ):51-53。

［10］武秀荣,杨学工大学物理［M］北京：中国农业出版社.2009：97-102．

26