“质面的位、球面的位、引力”这三个词放在一起乍看像是物理教材的目录碎片实际上是一条非常完整的逻辑链把质量从“点”解放到“面”研究一张质量面产生的引力位再把平面弯成球壳看“球面的位”如何在内外侧呈现完全不同的性格最后回到引力本身用位函数去解释那些我们凭直觉很难消化的事实。这篇笔记适合正在学经典力学或理论力学、被各种质量分布模型绕晕的人也适合想弄清楚“为什么球壳内部不受引力”的好奇型读者。我会尽量把推导过程讲清楚也会把我在计算中踩过的坑放在最后。1. 老问题的新角度质点模型为什么不够用1.1 质点模型的便利与局限万有引力定律的标准形式是对两个“质点”写的它们之间的距离 r 直接决定了相互作用。这个形式干净利落几乎所有入门课程都从这里切入。但真实存在的物体几乎都不是点你面前一个装满水的球壳、一块金属薄板、一条星系悬臂都是连续的质量堆叠。把每个质量元当成质点再做积分这听起来简单真正算起来却很容易被几何绕进去。质点模型的另一层局限是“位置”本身变得模糊。一个物体内部的某一点和物体外部的一点受到的引力贡献完全不同比如站在一个均匀球壳里面球壳各方向对你的吸引严格抵消站在外面则等效于全部质量集中在球心。这种差异用质点模型根本推不出来必须先引入连续分布的“面”“球”等形状。1.2 “位”到底是什么“位”就是引力势。它的正式定义是单位质量在引力场中具有的势能通常取无穷远处为零点记为 Φ(r)。引力场 g 与位的关系是 g -∇Φ也就是场强等于位函数下降最快的方向。这个关系的重要性在于直接算矢量场经常需要分解方向而先算标量场位只需要积分得到后再求梯度很多时候省掉一半麻烦。生活化理解位函数像一张地形的海拔图引力场则像当地“最陡的坡”。同一个海拔高度上你站得再高如果某处是平地的“势阱”那里的坡度依然可以是零我们关心引力往往也同时关心“高度”本身——比如从表面把物体搬到无穷远要克服多少能量这就是势差在做功。1.3 为什么先学“面”再学“球”后面第三部分你会看到匀质球壳的结论可以用很漂亮的几何解释但几何解释需要你先把“面”上的质量元如何对场贡献想清楚。先算无限大平面这个最简单的面分布能帮你建立微元分解和方向投影的直觉再到球壳时同一个方法加一个对称性就够了。2. 无限大质面的“位”一张平面能产生恒定场2.1 设定与微元分解取 x-y 平面为一个无限大的匀质“质面”面质量密度为 σ单位面积的质量kg/m²观察点 P 在 z 轴正方向距离平面高度为 h。由于平面无限任何水平偏移都不改变物理图景问题天然只有一维沿着 z 轴。把平面划分成以 P 到平面的垂足为圆心的圆环。半径在 ρ 到 ρdρ 之间的圆环面积是 2πρ dρ质量 dm σ·2πρ dρ。这个圆环上所有质量元到 P 的距离都为 r √(ρ²h²)且由于对称性水平方向合力为零只剩沿 z 方向的分量。每个质量元产生的引力沿连线方向与 z 轴的夹角为 θ垂直分量的比例是 cosθ h/r。这样直接写出贡献dF_z G·dm·cosθ / r² (G·σ·2πρ dρ·h) / (ρ²h²)^(3/2)这里唯一要小心的是r² 在分母上是距离平方而 cosθ 里的 h 和 r 也要同步代好。我最初就在这里把幂次算错过总是多写或少写一个 r导致积分发散或变成奇怪的表达式。2.2 积分结果与“与距离无关”的直观来源从 0 到无穷对 ρ 积分F_z 2πGσh ∫[0,∞] ρ/(ρ²h²)^(3/2) dρ令 u ρ²h²则 ρ dρ du/2分母变成 u^(3/2)从 uh² 积到无穷结果是 2/h。乘上前面的系数得到 F_z 2πGσ。注意这个结果里面 h 消失了。无限大匀质平面产生的引力场是垂直于平面的恒定场大小只由面密度 σ 决定与观测距离无关。我第一遍算出这个结果时完全不相信反复验算了几次才接受。直观解释也不复杂当你离平面更远时单个质量元对你的引力按 r² 衰减但同时你能“看到”的环带半径变大环带的面积从而质量按 ρ² 量级增长。衰减和增长刚好抵消于是总场恒定。这就好比在均匀起雾的天空里望出去远处的雾虽然更淡但看到的纵深更长整体亮度并不随距离明显变化。2.3 位函数与它的零点选择有了恒定场位函数立刻能写出来。取平面所在位置 z0 处位为零则在 z 点上方Φ(z) -∫[0,z] g·dz 2πGσ |z|因为场指向平面从 z 到 0 移动时势能降低但取绝对值后写法更方便。位函数呈 V 形在平面处最低往两侧线性上升。注意这里“位线性上升”和“场恒定”是互为印证的关系场是位的负梯度线性位的梯度是常数。这个 V 形位函数带来一个很有意思的现象如果你在质面上方释放一个物体它会以恒定加速度扑向平面但能量方面它每下降一段距离动能就线性增加相当于在一个均匀“斜坡”上加速。无限大平面因此能一直加速物体现实里当然没有无限大的引力平面但这并不妨碍它作为局部近似和教学模型的经典地位。2.4 有限盘和真实场景的衔接实际遇到更多的是有限半径的圆盘。离盘心很近、离盘面又近的区域有限盘的场和无限大平面几乎一致但一旦横向距离可与盘半径相比边缘效应就会出现场强开始下降。我在处理星系盘的简化模型时常常先用无限大平面估计中心区域的行为再在边缘处补一个修正。这个“先无限后有限”的思路在很多近似问题里都好用。3. 球面的“位”壳层内外的势与场如何两极分化3.1 球壳模型和三个结论现在把质量 M 均匀铺在一个半径 R 的球面上面密度 σ M/(4πR²)。这是一个匀质“球面”质面。它的引力结论可以浓缩成三句话球壳外部r R等效于全部质量集中在球心的质点场 g -GM/r²位 Φ -GM/r球壳内部r R场恒为零位恒为 Φ -GM/R球壳表面r R场在物理上取外部极限 -GM/R²位同为 -GM/R。这三句话是经典力学的“壳层定理”内容也是许多天文模型中反复使用的基石。先说外部结论从足够远处看球壳和点没有区别这是对称性的必然但内部结论就一点都不平庸了它告诉你一个匀质球壳内部处处不受引力。3.2 用微元投射法推导内部零场内部零场有一个漂亮的初等证明。把球壳看成无数对“对面元”观察点 P 在球内。任取一条过 P 的直线它与球壳交于 A、B 两个对称点方向相反。设 P 到 A 的距离为 a到 B 的距离为 b这两个小面元相对 P 的张角相同而面元面积正比于距离的平方A 处面元质量 dM_A 正比于 a²B 处面元质量 dM_B 正比于 b²。每个质量元对 P 的引力贡献正比于 dM/距离²于是 A 处贡献正比于 a²/a² 1B 处同样是 b²/b² 1。方向相反、大小相等自然抵消。把整个球壳按这种“对穿点”配完每一对都抵消合场强为零。这个证明不需要任何积分我第一次看到时觉得它比长篇积分漂亮得多。它也顺便解释了为什么球壳内部的结论如此干净不是某个分布凑巧而是每一对面元都精确相消。3.3 内部位为什么是常数而不为零“内部场为零”和“内部位是 -GM/R”这两件事并不矛盾场是位的负梯度某区域内梯度处处为零位就是常数常数是多少由边界处的连续性决定。球壳表面外侧位是 -GM/R内侧位也必须等于这个值否则跨越一个没有厚度的面时位函数会跳变这在物理上不允许位与做功路径有关不允许凭空多出能量差。于是球壳内部的物体虽然不受净力但它和无穷远之间的能量差依然存在从无穷远把物体搬进球壳内部总能量净减少 GMm/R。换句话说球壳是困住内部物体的“势阱”但阱底是一块完全失重的平地。我在科普里见过太多人把“不受力”误解成“没有势能”后面专门把这条列入避坑清单。3.4 从球壳到实心球场的分段表达球壳做完立刻可以叠加出均匀实心球总质量 M、半径 R的结果。把实心球想象成无数层同心的薄球壳叠加对外部而言每一层都等效为球心质点所以 r R 时整体等效为总质量 M 的质点对内部半径 r 处的观察点半径大于 r 的壳层全部不产生作用只有半径小于 r 的核心层在起作用。核心层质量 M_core M(r/R)³所以场强 g(r) -G M_core / r² -GM r / R³。这给出一个线性恢复力均匀实心球内部的引力场随半径线性增长在中心为零在表面达到 GM/R²。位函数也不难积分出来Φ(r) -3GM/(2R) GM r²/(2R³)内部Φ(r) -GM/r外部。注意这里中心位是 -3GM/(2R)比表面 -GM/R 还要低一半。也就是说一个物体从表面掉向球心势能一路走低动能越来越大到了中心虽然合力为零但速度却最大。这个“中心零力但速度最大”的细节我当年学的时候绕了很久。3.5 一个表格总结三场景下面这组表可以直接用来对照检查我的计算质量分布区域场强大小位函数无限大质面任一侧2πGσ恒定2πGσ|z|取平面为零点球壳外部r RGM/r²-GM/r球壳内部r R0-GM/R常数实心球外部r RGM/r²-GM/r实心球内部r RGM r/R³-3GM/(2R) GM r²/(2R³)看到内部场为零的同时还要看着位函数的数值——这两列经常被初学者混为一谈。4. 从位到引力梯度关系与实际计算中的桥接4.1 位函数为什么先算更省力直接算引力场是矢量的活三轴分量都要小心分解位函数是标量只需要一次积分。我处理“一个有限质面的位”这类问题时几乎总是先努力凑出 Φ再用 φ 求梯度。因为标量积分里对称性可以一次用完比如球对称直接只留径向分量而矢量积分要反复投影。以质面为例如果一开始就分别求 F_x、F_y、F_z三个方向都是复杂积分但先求出 V 形的 Φ再对 z 求导一步就得到 F_z水平方向因为对称性直接为零。这就是“位”的实用价值它替你把一场复杂的力场压缩成一个标量地形用的时候再展开。4.2 能量视角位函数如何直接参与轨道计算知道 Φ(r) 之后单位质量的机械能守恒可以写成½v² Φ(r) 常数。在一个球壳外部Φ -GM/r于是 ½v² - GM/r E。这个式子能推导轨道速度和逃逸速度。在球壳内部Φ -GM/R 是常数能量守恒直接保证内部物体自由运动时速率不变——不受力自然匀速直线运动与能量守恒完全自洽。这给零场区域一个动力学诠释内部物体确实在做匀速直线运动只是它的总能量比在无穷远静止时少 GM/R。这又是一个典型的“位不同力相同”的例子。4.3 泊松方程把面源和场跳跃连起来更统一地看待“质面”与“引力”可以用场论里的泊松方程∇²Φ 4πGρ。对无限大质面ρ 在平面处是一个 δ 函数式的面密度解出来的 Φ 是 V 形其二阶导数在平面处体现为面质量的“源”。跨过面层时场强法向分量发生跳变跳变量是 4πGσ。这个 4πGσ 的系数在解完无限大平面后其实已经摸到了上下两边的场都是 2πGσ方向相反跨面正负号翻转差值正好 4πGσ。记得这个关系后面处理多层球壳、两层介质交界时就不用每次重新推导。4.4 位函数连续但场可以不连续很多人在球壳表面会疑惑位函数在表面两边都是 -GM/R是连续的但场强呢外部是 GM/R²内部是 0跨过表面场强跳变了。这件事听上去违反直觉但完全合法位函数只要求连续场位的导数允许在面源处间断。就像地形高度不会在地图边界上突然跳变但坡度可以突然变化。这个“连续但导数间断”的点是所有面质量模型共有的特征也是用好位函数的关键。5. 我反复绕进去的几个坑势、场与边界5.1 坑一零场 零势球壳内部场为零很多人直接下结论“里面没有引力势能”。这是错的。势是能量角度的量场是受力度角度的量零力不等于零能量。拿最简单的例子说你先停在地面再站到屋顶两个位置的受力差别不大但势能差很大球壳内部就类似它是一个恒常的势能底不是没有势。5.2 坑二无限大质面的“无限”陷阱算无限大平面时场强与距离无关但位函数却线性增长。这意味着从平面延伸到无穷远位会变成无穷大所以有没有有限零点是个问题。实际应用里要么像我在 2.3 那样取平面处为零点要么只关心两点之间的势差绝不谈“平面上方的绝对势”。在模拟中通过截断半径把平面改成有限盘问题才会回到有限情形。5.3 坑三符号和零点约定不一致位函数里负号、零点是重灾区。有人习惯把无穷远定为零点球壳、球体场景有人习惯把平面定为零点质面场景只要混用梯度方向就会出错。我自己的习惯是凡是“有限质量分布在有限区域”的情况一律把无穷远设为零凡是“无限延展模型”明确写出零点并提醒自己只比较差。这个约定帮我少错了很多。5.4 坑四壳内零场不等于实心球内部零场球壳内部场为零实心球内部场不是零而是线性增长。这两个场景经常被混记。区别在于实心球内部有“体内层”的质量它会对内部偏外侧的观测点产生引力而球壳内部没有任何质量壳层位于你与球心之间都在你外侧。想清楚这句就再也不会混。5.5 一些适合动手做的练习如果想把这一套彻底消化我的建议是按顺序做三次手算第一次算无限大平面感受微元积分的对称性第二次算球壳用对穿点几何证明内部零场第三次把实心球当成无穷层球壳叠加推一遍内外部位函数。三次推完质面、球面、引力这三块就基本锁死在脑子里了。我自己的体会是位函数这东西第一眼很抽象但一旦习惯了“先标量后梯度”的思路它会变成你手里最顺手的工具。最后补一句我个人的感受连续质量分布的引力学真正困难的地方从来不是积分技巧而是把“位置”“势”“力”三个概念彻底分开、再重新焊到一起。我能给出的最好建议就是多画几张 V 形和 -GM/r 的图亲手算完一个模型再进入下一个。这篇笔记写到这里希望对正在啃经典力学的朋友有一点点帮助。
“质面的位、球面的位、引力”这三个词放在一起乍看像是物理教材的目录碎片实际上是一条非常完整的逻辑链把质量从“点”解放到“面”研究一张质量面产生的引力位再把平面弯成球壳看“球面的位”如何在内外侧呈现完全不同的性格最后回到引力本身用位函数去解释那些我们凭直觉很难消化的事实。这篇笔记适合正在学经典力学或理论力学、被各种质量分布模型绕晕的人也适合想弄清楚“为什么球壳内部不受引力”的好奇型读者。我会尽量把推导过程讲清楚也会把我在计算中踩过的坑放在最后。1. 老问题的新角度质点模型为什么不够用1.1 质点模型的便利与局限万有引力定律的标准形式是对两个“质点”写的它们之间的距离 r 直接决定了相互作用。这个形式干净利落几乎所有入门课程都从这里切入。但真实存在的物体几乎都不是点你面前一个装满水的球壳、一块金属薄板、一条星系悬臂都是连续的质量堆叠。把每个质量元当成质点再做积分这听起来简单真正算起来却很容易被几何绕进去。质点模型的另一层局限是“位置”本身变得模糊。一个物体内部的某一点和物体外部的一点受到的引力贡献完全不同比如站在一个均匀球壳里面球壳各方向对你的吸引严格抵消站在外面则等效于全部质量集中在球心。这种差异用质点模型根本推不出来必须先引入连续分布的“面”“球”等形状。1.2 “位”到底是什么“位”就是引力势。它的正式定义是单位质量在引力场中具有的势能通常取无穷远处为零点记为 Φ(r)。引力场 g 与位的关系是 g -∇Φ也就是场强等于位函数下降最快的方向。这个关系的重要性在于直接算矢量场经常需要分解方向而先算标量场位只需要积分得到后再求梯度很多时候省掉一半麻烦。生活化理解位函数像一张地形的海拔图引力场则像当地“最陡的坡”。同一个海拔高度上你站得再高如果某处是平地的“势阱”那里的坡度依然可以是零我们关心引力往往也同时关心“高度”本身——比如从表面把物体搬到无穷远要克服多少能量这就是势差在做功。1.3 为什么先学“面”再学“球”后面第三部分你会看到匀质球壳的结论可以用很漂亮的几何解释但几何解释需要你先把“面”上的质量元如何对场贡献想清楚。先算无限大平面这个最简单的面分布能帮你建立微元分解和方向投影的直觉再到球壳时同一个方法加一个对称性就够了。2. 无限大质面的“位”一张平面能产生恒定场2.1 设定与微元分解取 x-y 平面为一个无限大的匀质“质面”面质量密度为 σ单位面积的质量kg/m²观察点 P 在 z 轴正方向距离平面高度为 h。由于平面无限任何水平偏移都不改变物理图景问题天然只有一维沿着 z 轴。把平面划分成以 P 到平面的垂足为圆心的圆环。半径在 ρ 到 ρdρ 之间的圆环面积是 2πρ dρ质量 dm σ·2πρ dρ。这个圆环上所有质量元到 P 的距离都为 r √(ρ²h²)且由于对称性水平方向合力为零只剩沿 z 方向的分量。每个质量元产生的引力沿连线方向与 z 轴的夹角为 θ垂直分量的比例是 cosθ h/r。这样直接写出贡献dF_z G·dm·cosθ / r² (G·σ·2πρ dρ·h) / (ρ²h²)^(3/2)这里唯一要小心的是r² 在分母上是距离平方而 cosθ 里的 h 和 r 也要同步代好。我最初就在这里把幂次算错过总是多写或少写一个 r导致积分发散或变成奇怪的表达式。2.2 积分结果与“与距离无关”的直观来源从 0 到无穷对 ρ 积分F_z 2πGσh ∫[0,∞] ρ/(ρ²h²)^(3/2) dρ令 u ρ²h²则 ρ dρ du/2分母变成 u^(3/2)从 uh² 积到无穷结果是 2/h。乘上前面的系数得到 F_z 2πGσ。注意这个结果里面 h 消失了。无限大匀质平面产生的引力场是垂直于平面的恒定场大小只由面密度 σ 决定与观测距离无关。我第一遍算出这个结果时完全不相信反复验算了几次才接受。直观解释也不复杂当你离平面更远时单个质量元对你的引力按 r² 衰减但同时你能“看到”的环带半径变大环带的面积从而质量按 ρ² 量级增长。衰减和增长刚好抵消于是总场恒定。这就好比在均匀起雾的天空里望出去远处的雾虽然更淡但看到的纵深更长整体亮度并不随距离明显变化。2.3 位函数与它的零点选择有了恒定场位函数立刻能写出来。取平面所在位置 z0 处位为零则在 z 点上方Φ(z) -∫[0,z] g·dz 2πGσ |z|因为场指向平面从 z 到 0 移动时势能降低但取绝对值后写法更方便。位函数呈 V 形在平面处最低往两侧线性上升。注意这里“位线性上升”和“场恒定”是互为印证的关系场是位的负梯度线性位的梯度是常数。这个 V 形位函数带来一个很有意思的现象如果你在质面上方释放一个物体它会以恒定加速度扑向平面但能量方面它每下降一段距离动能就线性增加相当于在一个均匀“斜坡”上加速。无限大平面因此能一直加速物体现实里当然没有无限大的引力平面但这并不妨碍它作为局部近似和教学模型的经典地位。2.4 有限盘和真实场景的衔接实际遇到更多的是有限半径的圆盘。离盘心很近、离盘面又近的区域有限盘的场和无限大平面几乎一致但一旦横向距离可与盘半径相比边缘效应就会出现场强开始下降。我在处理星系盘的简化模型时常常先用无限大平面估计中心区域的行为再在边缘处补一个修正。这个“先无限后有限”的思路在很多近似问题里都好用。3. 球面的“位”壳层内外的势与场如何两极分化3.1 球壳模型和三个结论现在把质量 M 均匀铺在一个半径 R 的球面上面密度 σ M/(4πR²)。这是一个匀质“球面”质面。它的引力结论可以浓缩成三句话球壳外部r R等效于全部质量集中在球心的质点场 g -GM/r²位 Φ -GM/r球壳内部r R场恒为零位恒为 Φ -GM/R球壳表面r R场在物理上取外部极限 -GM/R²位同为 -GM/R。这三句话是经典力学的“壳层定理”内容也是许多天文模型中反复使用的基石。先说外部结论从足够远处看球壳和点没有区别这是对称性的必然但内部结论就一点都不平庸了它告诉你一个匀质球壳内部处处不受引力。3.2 用微元投射法推导内部零场内部零场有一个漂亮的初等证明。把球壳看成无数对“对面元”观察点 P 在球内。任取一条过 P 的直线它与球壳交于 A、B 两个对称点方向相反。设 P 到 A 的距离为 a到 B 的距离为 b这两个小面元相对 P 的张角相同而面元面积正比于距离的平方A 处面元质量 dM_A 正比于 a²B 处面元质量 dM_B 正比于 b²。每个质量元对 P 的引力贡献正比于 dM/距离²于是 A 处贡献正比于 a²/a² 1B 处同样是 b²/b² 1。方向相反、大小相等自然抵消。把整个球壳按这种“对穿点”配完每一对都抵消合场强为零。这个证明不需要任何积分我第一次看到时觉得它比长篇积分漂亮得多。它也顺便解释了为什么球壳内部的结论如此干净不是某个分布凑巧而是每一对面元都精确相消。3.3 内部位为什么是常数而不为零“内部场为零”和“内部位是 -GM/R”这两件事并不矛盾场是位的负梯度某区域内梯度处处为零位就是常数常数是多少由边界处的连续性决定。球壳表面外侧位是 -GM/R内侧位也必须等于这个值否则跨越一个没有厚度的面时位函数会跳变这在物理上不允许位与做功路径有关不允许凭空多出能量差。于是球壳内部的物体虽然不受净力但它和无穷远之间的能量差依然存在从无穷远把物体搬进球壳内部总能量净减少 GMm/R。换句话说球壳是困住内部物体的“势阱”但阱底是一块完全失重的平地。我在科普里见过太多人把“不受力”误解成“没有势能”后面专门把这条列入避坑清单。3.4 从球壳到实心球场的分段表达球壳做完立刻可以叠加出均匀实心球总质量 M、半径 R的结果。把实心球想象成无数层同心的薄球壳叠加对外部而言每一层都等效为球心质点所以 r R 时整体等效为总质量 M 的质点对内部半径 r 处的观察点半径大于 r 的壳层全部不产生作用只有半径小于 r 的核心层在起作用。核心层质量 M_core M(r/R)³所以场强 g(r) -G M_core / r² -GM r / R³。这给出一个线性恢复力均匀实心球内部的引力场随半径线性增长在中心为零在表面达到 GM/R²。位函数也不难积分出来Φ(r) -3GM/(2R) GM r²/(2R³)内部Φ(r) -GM/r外部。注意这里中心位是 -3GM/(2R)比表面 -GM/R 还要低一半。也就是说一个物体从表面掉向球心势能一路走低动能越来越大到了中心虽然合力为零但速度却最大。这个“中心零力但速度最大”的细节我当年学的时候绕了很久。3.5 一个表格总结三场景下面这组表可以直接用来对照检查我的计算质量分布区域场强大小位函数无限大质面任一侧2πGσ恒定2πGσ|z|取平面为零点球壳外部r RGM/r²-GM/r球壳内部r R0-GM/R常数实心球外部r RGM/r²-GM/r实心球内部r RGM r/R³-3GM/(2R) GM r²/(2R³)看到内部场为零的同时还要看着位函数的数值——这两列经常被初学者混为一谈。4. 从位到引力梯度关系与实际计算中的桥接4.1 位函数为什么先算更省力直接算引力场是矢量的活三轴分量都要小心分解位函数是标量只需要一次积分。我处理“一个有限质面的位”这类问题时几乎总是先努力凑出 Φ再用 φ 求梯度。因为标量积分里对称性可以一次用完比如球对称直接只留径向分量而矢量积分要反复投影。以质面为例如果一开始就分别求 F_x、F_y、F_z三个方向都是复杂积分但先求出 V 形的 Φ再对 z 求导一步就得到 F_z水平方向因为对称性直接为零。这就是“位”的实用价值它替你把一场复杂的力场压缩成一个标量地形用的时候再展开。4.2 能量视角位函数如何直接参与轨道计算知道 Φ(r) 之后单位质量的机械能守恒可以写成½v² Φ(r) 常数。在一个球壳外部Φ -GM/r于是 ½v² - GM/r E。这个式子能推导轨道速度和逃逸速度。在球壳内部Φ -GM/R 是常数能量守恒直接保证内部物体自由运动时速率不变——不受力自然匀速直线运动与能量守恒完全自洽。这给零场区域一个动力学诠释内部物体确实在做匀速直线运动只是它的总能量比在无穷远静止时少 GM/R。这又是一个典型的“位不同力相同”的例子。4.3 泊松方程把面源和场跳跃连起来更统一地看待“质面”与“引力”可以用场论里的泊松方程∇²Φ 4πGρ。对无限大质面ρ 在平面处是一个 δ 函数式的面密度解出来的 Φ 是 V 形其二阶导数在平面处体现为面质量的“源”。跨过面层时场强法向分量发生跳变跳变量是 4πGσ。这个 4πGσ 的系数在解完无限大平面后其实已经摸到了上下两边的场都是 2πGσ方向相反跨面正负号翻转差值正好 4πGσ。记得这个关系后面处理多层球壳、两层介质交界时就不用每次重新推导。4.4 位函数连续但场可以不连续很多人在球壳表面会疑惑位函数在表面两边都是 -GM/R是连续的但场强呢外部是 GM/R²内部是 0跨过表面场强跳变了。这件事听上去违反直觉但完全合法位函数只要求连续场位的导数允许在面源处间断。就像地形高度不会在地图边界上突然跳变但坡度可以突然变化。这个“连续但导数间断”的点是所有面质量模型共有的特征也是用好位函数的关键。5. 我反复绕进去的几个坑势、场与边界5.1 坑一零场 零势球壳内部场为零很多人直接下结论“里面没有引力势能”。这是错的。势是能量角度的量场是受力度角度的量零力不等于零能量。拿最简单的例子说你先停在地面再站到屋顶两个位置的受力差别不大但势能差很大球壳内部就类似它是一个恒常的势能底不是没有势。5.2 坑二无限大质面的“无限”陷阱算无限大平面时场强与距离无关但位函数却线性增长。这意味着从平面延伸到无穷远位会变成无穷大所以有没有有限零点是个问题。实际应用里要么像我在 2.3 那样取平面处为零点要么只关心两点之间的势差绝不谈“平面上方的绝对势”。在模拟中通过截断半径把平面改成有限盘问题才会回到有限情形。5.3 坑三符号和零点约定不一致位函数里负号、零点是重灾区。有人习惯把无穷远定为零点球壳、球体场景有人习惯把平面定为零点质面场景只要混用梯度方向就会出错。我自己的习惯是凡是“有限质量分布在有限区域”的情况一律把无穷远设为零凡是“无限延展模型”明确写出零点并提醒自己只比较差。这个约定帮我少错了很多。5.4 坑四壳内零场不等于实心球内部零场球壳内部场为零实心球内部场不是零而是线性增长。这两个场景经常被混记。区别在于实心球内部有“体内层”的质量它会对内部偏外侧的观测点产生引力而球壳内部没有任何质量壳层位于你与球心之间都在你外侧。想清楚这句就再也不会混。5.5 一些适合动手做的练习如果想把这一套彻底消化我的建议是按顺序做三次手算第一次算无限大平面感受微元积分的对称性第二次算球壳用对穿点几何证明内部零场第三次把实心球当成无穷层球壳叠加推一遍内外部位函数。三次推完质面、球面、引力这三块就基本锁死在脑子里了。我自己的体会是位函数这东西第一眼很抽象但一旦习惯了“先标量后梯度”的思路它会变成你手里最顺手的工具。最后补一句我个人的感受连续质量分布的引力学真正困难的地方从来不是积分技巧而是把“位置”“势”“力”三个概念彻底分开、再重新焊到一起。我能给出的最好建议就是多画几张 V 形和 -GM/r 的图亲手算完一个模型再进入下一个。这篇笔记写到这里希望对正在啃经典力学的朋友有一点点帮助。
平方复杂度无法被常规优化消除｜TaoToken 统一 Key 通道实测)
