本文追溯经典物理学从伽利略到麦克斯韦近300年的发展历程,展示力学与电磁学如何从分散的自然观察逐步成长为严密的理论体系,最终实现光的电磁理论统一。
经典物理学(Classical Physics)是指20世纪初相对论和量子力学诞生之前建立的物理理论体系。它以牛顿力学(Newtonian Mechanics)和麦克斯韦电磁学(Maxwellian Electromagnetism)为两大支柱,构成了近代科学的第一个完整范式。
经典物理学之所以重要,不仅是因为它奠定了现代工程技术的理论基础,更在于它建立了"用数学语言描述自然规律"的科学方法——这一方法论至今仍是所有自然科学的核心范式。
经典物理学的发展大致可分为三个阶段:
时期
时间跨度
核心成就
代表人物
奠基期
1543-1687
日心说、行星运动定律、力学三定律
哥白尼、开普勒、伽利略、牛顿
拓展期
1687-1800
分析力学、热学理论、电学开端
拉格朗日、拉普拉斯、卡诺、库仑
统一期
1800-1890
电磁理论统一、能量守恒、统计力学
法拉第、麦克斯韦、焦耳、玻尔兹曼
以具体数据说明这一体系的成就:1687年牛顿发表《自然哲学的数学原理》(Philosophiæ Naturalis Principia Mathematica ),用3条运动定律和1条万有引力定律,统一解释了地球上苹果下落和天空中行星运动的规律。而200年后,麦克斯韦用4个方程(麦克斯韦方程组)统一了电学、磁学和光学。这种"以最少原理解释最多现象"的追求,成为物理学最鲜明的精神气质。
伽利略·伽利莱(Galileo Galilei,1564-1642)被誉为"现代科学之父"。他的核心贡献不在于发现了某个具体定律,而在于确立了实验与数学相结合的科学方法 ——即先通过可控实验获取数据,再通过数学分析发现函数关系。
斜面实验:定量研究的典范
伽利略著名的斜面实验,是他挑战亚里士多德运动学说的关键证据。他让小球从不同倾斜角度的斜面滚下,通过巧妙设计的水钟测量时间。
设想一个斜面长度为 L = 2 m L = 2\ \text{m}L = 2 m ,高度为 h = 0.5 m h = 0.5\ \text{m}h = 0 . 5 m ,则倾角 θ \thetaθ 满足:
sin θ = h L = 0.5 2 = 0.25 \sin \theta = \frac{h}{L} = \frac{0.5}{2} = 0.25
sin θ = L h = 2 0 . 5 = 0 . 2 5
小球沿斜面方向的加速度 a = g sin θ a = g \sin \thetaa = g sin θ (其中 g ≈ 9.8 m/s 2 g \approx 9.8\ \text{m/s}^2g ≈ 9 . 8 m/s 2 ):
a = 9.8 × 0.25 = 2.45 m/s 2 a = 9.8 \times 0.25 = 2.45\ \text{m/s}^2
a = 9 . 8 × 0 . 2 5 = 2 . 4 5 m/s 2
由匀加速运动公式 s = 1 2 a t 2 s = \frac{1}{2}at^2s = 2 1 a t 2 ,球从静止滚下1米所需时间:
t = 2 s a = 2 × 1 2.45 ≈ 0.90 s t = \sqrt{\frac{2s}{a}} = \sqrt{\frac{2 \times 1}{2.45}} \approx 0.90\ \text{s}
t = a 2 s = 2 . 4 5 2 × 1 ≈ 0 . 9 0 s
伽利略通过大量这样的测量发现:
距离与时间的平方成正比 :s ∝ t 2 s \propto t^2s ∝ t 2
下落速度与质量无关 :轻重物体在同一高度释放,同时落地(排除空气阻力)
后者在著名的比萨斜塔传说中形象地体现出来。虽然这个故事可能并不真实,但伽利略确实通过思想实验证明:如果亚里士多德是对的(重物落得快),那么把轻物和重物绑在一起,会出现逻辑矛盾。
伽利略的惯性原理
伽利略通过另一个思想实验——让小球在理想光滑平面上滚动——推翻了亚里士多德"力是维持运动的原因"的观点。他论证:一个物体在无外力作用时,将保持匀速直线运动。这成为牛顿第一定律的直接前身。
在牛顿之前,约翰内斯·开普勒(Johannes Kepler,1571-1630)利用第谷·布拉赫(Tycho Brahe)长达20年的精确天文观测数据,总结出了行星运动的三条定律。
第谷数据与开普勒的计算
第谷的观测精度达到 ± 4 \pm 4± 4 角分(约 0.067 0.0670 . 0 6 7 度),是前无古人的。开普勒以火星为主要研究对象,因为它的轨道偏离圆形最明显。经过长达6年的计算——如果手算,开普勒进行了约70次反复尝试——他最终发现:
第一定律(椭圆轨道定律,1609年): 所有行星绕太阳运动的轨道都是椭圆,太阳位于椭圆的一个焦点上。
假设椭圆的半长轴为 a aa ,半短轴为 b bb ,偏心率为 e ee ,则满足:
e = 1 − b 2 a 2 e = \sqrt{1 - \frac{b^2}{a^2}}
e = 1 − a 2 b 2
例如地球轨道:a ≈ 1.496 × 1 0 11 m a \approx 1.496 \times 10^{11}\ \text{m}a ≈ 1 . 4 9 6 × 1 0 1 1 m ,e ≈ 0.0167 e \approx 0.0167e ≈ 0 . 0 1 6 7 (接近圆形);而水星轨道 e ≈ 0.2056 e \approx 0.2056e ≈ 0 . 2 0 5 6 ,明显偏离圆形。
第二定律(面积定律,1609年): 行星与太阳的连线在相等时间内扫过相等的面积。
用公式表示:设 d A / d t dA/dtd A / d t 为面积变化率,r rr 为日-星距离,v ⊥ v_\perpv ⊥ 为垂直于径向的速度分量,则:
d A d t = 1 2 r v ⊥ = 常量 \frac{dA}{dt} = \frac{1}{2} r v_\perp = \text{常量}
d t d A = 2 1 r v ⊥ = 常量
这意味着行星在近日点(离太阳最近)时运动最快,在远日点时最慢。以地球为例:
近日点(1月初):r peri ≈ 1.471 × 1 0 11 m r_{\text{peri}} \approx 1.471 \times 10^{11}\ \text{m}r peri ≈ 1 . 4 7 1 × 1 0 1 1 m ,轨道速度约 30.3 km/s 30.3\ \text{km/s}3 0 . 3 km/s
远日点(7月初):r apj ≈ 1.521 × 1 0 11 m r_{\text{apj}} \approx 1.521 \times 10^{11}\ \text{m}r apj ≈ 1 . 5 2 1 × 1 0 1 1 m ,轨道速度约 29.3 km/s 29.3\ \text{km/s}2 9 . 3 km/s
两者速度差约为 1 km/s 1\ \text{km/s}1 km/s (约3.3%),这就是为什么北半球冬季比夏季稍短(因为地球在近日点运动更快)。
第三定律(周期定律,1619年): 行星绕太阳公转周期的平方与轨道半长轴的立方成正比。
T 2 ∝ a 3 或 T 2 a 3 = 常数 T^2 \propto a^3 \quad \text{或} \quad \frac{T^2}{a^3} = \text{常数}
T 2 ∝ a 3 或 a 3 T 2 = 常数
下表列出太阳系行星的轨道参数,验证第三定律:
行星
半长轴 a aa (AU)
周期 T TT (年)
a 3 a^3a 3
T 2 T^2T 2
T 2 / a 3 T^2/a^3T 2 / a 3
水星
0.387
0.241
0.058
0.058
1.000
金星
0.723
0.615
0.378
0.378
1.000
地球
1.000
1.000
1.000
1.000
1.000
火星
1.524
1.881
3.540
3.538
0.999
木星
5.203
11.862
140.8
140.7
0.999
表中可见 T 2 / a 3 T^2/a^3T 2 / a 3 几乎精确恒定,误差不到 0.2%。这个精确的比例关系后来成为牛顿推导万有引力定律的关键线索。
伊萨克·牛顿(Isaac Newton,1643-1727)在1687年出版的《自然哲学的数学原理》中,集前人研究之大成,建立了完整的力学体系。
牛顿三大运动定律
牛顿第一定律(惯性定律):任何物体都保持静止或匀速直线运动状态,除非外力迫使它改变这种状态。
牛顿第二定律:物体所受合力等于其质量与加速度的乘积。
∑ F ⃗ = m a ⃗ \sum \vec{F} = m \vec{a}
∑ F = m a
这是一个矢量方程。以一维情况为例:一个质量为 m = 5 kg m = 5\ \text{kg}m = 5 kg 的物体,受到 F 1 = 20 N F_1 = 20\ \text{N}F 1 = 2 0 N 向右和 F 2 = 8 N F_2 = 8\ \text{N}F 2 = 8 N 向左的力,则合力 ∑ F = 20 − 8 = 12 N \sum F = 20 - 8 = 12\ \text{N}∑ F = 2 0 − 8 = 1 2 N ,加速度:
a = ∑ F m = 12 5 = 2.4 m/s 2 a = \frac{\sum F}{m} = \frac{12}{5} = 2.4\ \text{m/s}^2
a = m ∑ F = 5 1 2 = 2 . 4 m/s 2
方向向右。
牛顿第三定律(作用力与反作用力定律):两个物体之间的作用力和反作用力大小相等、方向相反,且作用在同一直线上。
F ⃗ 12 = − F ⃗ 21 \vec{F}_{12} = -\vec{F}_{21}
F 1 2 = − F 2 1
一个常见的误解是:既然作用力与反作用力大小相等,为什么物体还能运动?关键在于它们作用在不同物体上 。比如人推墙,墙也推人——但墙的推力作用在人身上,推动人向反方向运动(即人后退)。
万有引力定律
牛顿从开普勒第三定律和牛顿第二定律出发,推导出万有引力公式。他的推理大致如下:
考虑质量为 m mm 的行星绕质量为 M MM 的太阳做匀速圆周运动,由牛顿第二定律:
F = m a = m ⋅ v 2 r F = m a = m \cdot \frac{v^2}{r}
F = m a = m ⋅ r v 2
其中 v vv 为轨道速度,r rr 为轨道半径。周期 T = 2 π r / v T = 2\pi r / vT = 2 π r / v ,代入得:
F = m ⋅ 4 π 2 r T 2 F = m \cdot \frac{4\pi^2 r}{T^2}
F = m ⋅ T 2 4 π 2 r
根据开普勒第三定律 T 2 ∝ r 3 T^2 \propto r^3T 2 ∝ r 3 ,设 T 2 = k r 3 T^2 = k r^3T 2 = k r 3 (k kk 为常数),则:
F = m ⋅ 4 π 2 r k r 3 = 4 π 2 k ⋅ m r 2 F = m \cdot \frac{4\pi^2 r}{k r^3} = \frac{4\pi^2}{k} \cdot \frac{m}{r^2}
F = m ⋅ k r 3 4 π 2 r = k 4 π 2 ⋅ r 2 m
也就是说,引力与距离平方成反比,与质量成正比。再根据牛顿第三定律,引力是相互的,所以引力也应该与太阳质量 M MM 成正比。因此:
F = G M m r 2 F = G \frac{M m}{r^2}
F = G r 2 M m
其中 G \approx 6.674 \times 10^-11\ \textN·m^2/\textkg^2 为万有引力常数。
万有引力的验证:哈雷彗星
一个著名的验证案例是哈雷彗星(Halley's Comet)。埃德蒙·哈雷(Edmond Halley)在1705年利用牛顿的引力理论计算了24颗彗星的轨道,发现1531年、1607年和1682年出现的三颗彗星具有相似的轨道参数,从而预言它们实际上是同一颗彗星,每约76年回归一次。他大胆预言该彗星将于1758年再次回归——虽然哈雷本人未能活到那一天,但1758年圣诞节,彗星如约而至,天文学界一片沸腾。这是牛顿引力理论的第一次重大成功预测。
牛顿力学的意义远超物理本身,它建立了一种被后世称为"牛顿范式"的科学方法论:
公理化体系 :从少量基本定律(三大运动定律+万有引力定律)出发,通过严格数学推导,演绎出所有力学现象的解释
因果确定性 :给定初始条件(位置和速度),物体的未来运动轨迹被唯一确定——这就是决定论(Determinism)的思想根源
数学化描述 :使用微积分这一全新数学语言精确描述运动和力的关系
牛顿范式如此成功,以至于18世纪的科学家普遍相信"宇宙像时钟一样精确运行"——这正是拉普拉斯(Pierre-Simon Laplace)提出"拉普拉斯妖"(Daemon of Laplace)思想实验的背景:如果有一个智能体知道所有粒子的位置和速度,它就能预知宇宙的一切。
牛顿的矢量力学虽然完备,但在处理复杂约束系统时变得极为繁琐。18世纪的数学家们致力于寻找更优雅、更通用的力学表述。
拉格朗日力学
约瑟夫-路易·拉格朗日(Joseph-Louis Lagrange,1736-1813)在1788年出版的《分析力学》(Mécanique Analytique )中,完全用数学分析取代了几何方法。他引入广义坐标 (Generalized Coordinates)的概念——即用于完全描述系统状态的任意一组独立参数。
对于一个单摆(一个质量为 m mm 的球由长度为 l ll 的轻杆悬挂),如果用牛顿力学分析,需要处理张力 T TT 和重力 m g mgm g 两个力,并利用约束 x 2 + y 2 = l 2 x^2 + y^2 = l^2x 2 + y 2 = l 2 。而拉格朗日力学只需要一个广义坐标:摆角 θ \thetaθ 。
系统的动能 T TT 和势能 V VV 分别为:
T = 1 2 m ( l θ ˙ ) 2 = 1 2 m l 2 θ ˙ 2 T = \frac{1}{2} m (l \dot{\theta})^2 = \frac{1}{2} m l^2 \dot{\theta}^2
T = 2 1 m ( l θ ˙ ) 2 = 2 1 m l 2 θ ˙ 2
V = m g l ( 1 − cos θ ) V = m g l (1 - \cos \theta)
V = m g l ( 1 − cos θ )
拉格朗日量 L = T − V L = T - VL = T − V 。代入拉格朗日方程 (Lagrange's Equation):
d d t ( ∂ L ∂ θ ˙ ) − ∂ L ∂ θ = 0 \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{\theta}} \right) - \frac{\partial L}{\partial \theta} = 0
d t d ( ∂ θ ˙ ∂ L ) − ∂ θ ∂ L = 0
计算过程更简化——不出现约束力,不需要矢量分解。
数值示例:双摆系统
考虑一个更复杂的系统:双摆(Double Pendulum)——两个摆串联。用牛顿力学分析需要4个变量(两个位置的x,y坐标)和2个约束方程,但用拉格朗日力学只需要2个广义坐标(θ 1 , θ 2 \theta_1, \theta_2θ 1 , θ 2 )。下表对比两种方法的复杂度:
特征
牛顿力学
拉格朗日力学
变量数
4(x 1 , y 1 , x 2 , y 2 x_1, y_1, x_2, y_2x 1 , y 1 , x 2 , y 2 )
2(θ 1 , θ 2 \theta_1, \theta_2θ 1 , θ 2 )
约束方程
2(杆长恒定)
0(已通过广义坐标消去)
约束力计算
需要引入张力
自动消去
方程类型
矢量微分方程
标量微分方程
这正是拉格朗日力学的核心优势:由系统自由度的数目决定变量个数,而非由粒子个数决定。
威廉·罗恩·哈密顿(William Rowan Hamilton,1805-1865)进一步将力学改写为更对称的形式。他引入哈密顿量 H = T + V H = T + VH = T + V (系统总能量),以及广义动量 p i = ∂ L / ∂ q ˙ i p_i = \partial L / \partial \dot{q}_ip i = ∂ L / ∂ q ˙ i ,得到哈密顿方程:
q ˙ i = ∂ H ∂ p i , p ˙ i = − ∂ H ∂ q i \dot{q}_i = \frac{\partial H}{\partial p_i}, \quad \dot{p}_i = -\frac{\partial H}{\partial q_i}
q ˙ i = ∂ p i ∂ H , p ˙ i = − ∂ q i ∂ H
这一形式不仅优美对称,更重要的是揭示了力学与光学之间的深刻联系,为后来量子力学的建立奠定了基础。
一维谐振子示例
考虑一个质量为 m = 1 kg m = 1\ \text{kg}m = 1 kg ,劲度系数 k = 4 N/m k = 4\ \text{N/m}k = 4 N/m 的弹簧振子。
哈密顿量:
H = p 2 2 m + 1 2 k x 2 = p 2 2 + 2 x 2 H = \frac{p^2}{2m} + \frac{1}{2} k x^2 = \frac{p^2}{2} + 2 x^2
H = 2 m p 2 + 2 1 k x 2 = 2 p 2 + 2 x 2
代入哈密顿方程:
x ˙ = ∂ H ∂ p = p \dot{x} = \frac{\partial H}{\partial p} = p
x ˙ = ∂ p ∂ H = p
p ˙ = − ∂ H ∂ x = − 4 x \dot{p} = -\frac{\partial H}{\partial x} = -4x
p ˙ = − ∂ x ∂ H = − 4 x
将此方程改写为 x ¨ = p ˙ = − 4 x \ddot{x} = \dot{p} = -4xx ¨ = p ˙ = − 4 x ,这是一个简谐振动方程,角频率 ω = k / m = 2 rad/s \omega = \sqrt{k/m} = 2\ \text{rad/s}ω = k / m = 2 rad/s ,周期 T = 2 π / ω = π ≈ 3.14 s T = 2\pi/\omega = \pi \approx 3.14\ \text{s}T = 2 π / ω = π ≈ 3 . 1 4 s 。
给定初始条件 x ( 0 ) = 1 m x(0) = 1\ \text{m}x ( 0 ) = 1 m ,p ( 0 ) = 0 p(0) = 0p ( 0 ) = 0 (即从最大位移释放),系统在相空间(Phase Space)( x , p ) (x, p)( x , p ) 中的轨迹是一个椭圆:p 2 / 4 + x 2 = 1 p^2/4 + x^2 = 1p 2 / 4 + x 2 = 1 (因为总能量 E = H = 2 × 1 2 = 2 J E = H = 2 \times 1^2 = 2\ \text{J}E = H = 2 × 1 2 = 2 J ),半轴长分别为 a x = 1 m a_x = 1\ \text{m}a x = 1 m ,a_p = 2\ \textkg·m/s。
这种相空间轨迹 的几何图像,在统计力学和混沌理论中有极其重要的作用。
年份
贡献者
核心贡献
数学工具
1687
牛顿
三大定律、万有引力
几何+微积分(流数法)
1736
欧拉
刚体力学方程
微分方程
1743
达朗贝尔
达朗贝尔原理
虚功原理
1788
拉格朗日
分析力学、广义坐标
变分法
1834
哈密顿
正则方程、最小作用量原理
偏微分方程
18-19世纪,牛顿力学在天体力学中取得了惊人的成就。
海王星发现:纸面上的行星
最著名的案例是海王星的发现。1821年,法国天文学家布瓦尔(Alexis Bouvard)发现天王星的观测位置与牛顿力学预测的理论位置存在偏差——最大达到约 2 22 角分(约 0.033 0.0330 . 0 3 3 度),远大于观测误差。
1845年,英国数学家约翰·库奇·亚当斯(John Couch Adams)和法国数学家奥本·勒维耶(Urbain Le Verrier)独立地利用牛顿引力理论反推,计算出一个未知行星的轨道和位置。
以数据展示计算难度:假设天王星受到一颗未知行星的引力摄动,其轨道参数需要求解12个未知数(未知行星的6个轨道要素和它的质量)。勒维耶经过反复迭代计算,最终预测该行星位于宝瓶座黄经 32 6 ∘ 326^\circ3 2 6 ∘ 处。
1846年9月23日,柏林天文台的加勒(Johann Gottfried Galle)在勒维耶预测位置偏差仅 1 ∘ 1^\circ1 ∘ 的范围内发现了海王星——误差不到预测值的 0.3%。这是牛顿理论最辉煌的胜利之一。
理论值与观测值的对比:
参数
预测值(勒维耶1846)
实际测量值
误差
轨道半长轴 a aa
30.0 AU
30.07 AU
0.2%
轨道偏心率 e ee
0.1076
0.0086
较大(预测偏高)
轨道倾角 i ii
—
1.77°
—
黄经
326°
328°
0.6%
发现位置偏差
—
< 1°
—
牛顿力学为工业革命提供了理论支撑。简单的杠杆原理、滑轮系统、斜面等机械,都可以用牛顿力学精确分析。
一个实际工程计算:桥的受力分析
一座简支梁桥,跨度为 L = 10 m L = 10\ \text{m}L = 1 0 m ,承受均匀分布荷载 w = 5000 N/m w = 5000\ \text{N/m}w = 5 0 0 0 N/m (约相当于每米500公斤)。根据牛顿力学,桥梁的反力为:
R A = R B = w L 2 = 5000 × 10 2 = 25000 N R_A = R_B = \frac{wL}{2} = \frac{5000 \times 10}{2} = 25000\ \text{N}
R A = R B = 2 w L = 2 5 0 0 0 × 1 0 = 2 5 0 0 0 N
梁中点的最大弯矩:
M_\textmax = \fracwL^28 = \frac5000 \times 1008 = 62500\ \textN·m
如果梁是矩形截面,高 h = 0.5 m h = 0.5\ \text{m}h = 0 . 5 m ,宽 b = 0.3 m b = 0.3\ \text{m}b = 0 . 3 m ,则截面模量 W = b h 2 / 6 W = bh^2/6W = b h 2 / 6 ,最大弯曲应力:
σ max = M max W = 62500 × 6 0.3 × 0. 5 2 = 5 MPa \sigma_{\text{max}} = \frac{M_{\text{max}}}{W} = \frac{62500 \times 6}{0.3 \times 0.5^2} = 5\ \text{MPa}
σ max = W M max = 0 . 3 × 0 . 5 2 6 2 5 0 0 × 6 = 5 MPa
这个值远低于普通钢材的屈服强度(约 250 MPa 250\ \text{MPa}2 5 0 MPa ),安全系数高达50倍——这就是为什么19世纪的铁路桥梁至今仍能使用的力学原因之一。
牛顿力学同样孕育了流体力学。丹尼尔·伯努利(Daniel Bernoulli,1700-1782)在其1738年的著作《流体动力学》(Hydrodynamica )中提出了著名的伯努利方程。对于一个沿流线的不可压缩理想流体:
1 2 ρ v 2 + ρ g h + p = 常量 \frac{1}{2} \rho v^2 + \rho g h + p = \text{常量}
2 1 ρ v 2 + ρ g h + p = 常量
其中 ρ \rhoρ 为流体密度,v vv 为流速,g gg 为重力加速度,h hh 为高度,p pp 为压强。
数值示例:管流中的速度变化
设想水流在管道中流动,管道截面积从 A 1 = 1 m 2 A_1 = 1\ \text{m}^2A 1 = 1 m 2 缩小到 A 2 = 0.25 m 2 A_2 = 0.25\ \text{m}^2A 2 = 0 . 2 5 m 2 。假设在高度相同处,入口速度 v 1 = 2 m/s v_1 = 2\ \text{m/s}v 1 = 2 m/s ,水的密度 ρ = 1000 kg/m 3 \rho = 1000\ \text{kg/m}^3ρ = 1 0 0 0 kg/m 3 。
由连续性方程 A 1 v 1 = A 2 v 2 A_1 v_1 = A_2 v_2A 1 v 1 = A 2 v 2 :
v 2 = v 1 ⋅ A 1 A 2 = 2 × 1 0.25 = 8 m/s v_2 = v_1 \cdot \frac{A_1}{A_2} = 2 \times \frac{1}{0.25} = 8\ \text{m/s}
v 2 = v 1 ⋅ A 2 A 1 = 2 × 0 . 2 5 1 = 8 m/s
由伯努利方程(忽略高度差),入口压强 p 1 = 200000 Pa p_1 = 200000\ \text{Pa}p 1 = 2 0 0 0 0 0 Pa (约2个大气压):
p 2 = p 1 + 1 2 ρ ( v 1 2 − v 2 2 ) = 200000 + 500 × ( 4 − 64 ) = 200000 − 30000 = 170000 Pa p_2 = p_1 + \frac{1}{2} \rho (v_1^2 - v_2^2) = 200000 + 500 \times (4 - 64) = 200000 - 30000 = 170000\ \text{Pa}
p 2 = p 1 + 2 1 ρ ( v 1 2 − v 2 2 ) = 2 0 0 0 0 0 + 5 0 0 × ( 4 − 6 4 ) = 2 0 0 0 0 0 − 3 0 0 0 0 = 1 7 0 0 0 0 Pa
可见流速增大导致压强下降——这就是文丘里流量计(Venturi Meter)和飞机机翼产生升力的工作原理。
18世纪末到19世纪初,人们开始认识到"热"并非物质(燃素说),而是能量的一种形式。这一认识经历了曲折的实验和争论过程。
焦耳的热功当量实验
詹姆斯·焦耳(James Joule,1818-1889)通过一系列精密的实验确定了热和功之间的换算关系。他著名的实验装置是一个带搅拌器的隔热容器,通过悬挂重物下落驱动搅拌器转动,使水温度升高。
在一个典型实验中,质量为 m = 26.32 kg m = 26.32\ \text{kg}m = 2 6 . 3 2 kg 的重物从高度 h = 1.60 m h = 1.60\ \text{m}h = 1 . 6 0 m 处下落,使 M = 6.31 kg M = 6.31\ \text{kg}M = 6 . 3 1 kg 的水温度升高 Δ T = 0.155 ∘ C \Delta T = 0.155\ ^\circ\text{C}Δ T = 0 . 1 5 5 ∘ C 。
重物释放的势能(即做的功):
W = m g h = 26.32 × 9.8 × 1.60 ≈ 412.7 J W = m g h = 26.32 \times 9.8 \times 1.60 \approx 412.7\ \text{J}
W = m g h = 2 6 . 3 2 × 9 . 8 × 1 . 6 0 ≈ 4 1 2 . 7 J
水吸收的热量(根据比热容 c_\text水 = 4186\ \textJ/(kg·°C)):
Q = M c 水 Δ T = 6.31 × 4186 × 0.155 ≈ 4094 J Q = M c_{\text{水}} \Delta T = 6.31 \times 4186 \times 0.155 \approx 4094\ \text{J}
Q = M c 水 Δ T = 6 . 3 1 × 4 1 8 6 × 0 . 1 5 5 ≈ 4 0 9 4 J
注意:这个数据表示需要用大约412.7焦耳的机械功才能使水温度升高0.155°C。焦耳的贡献在于定量证明了:无论通过哪种方式做功(摩擦、搅拌、压缩),同样数量的功总是产生等量的热。
焦耳在1843-1878年间进行了大量实验,最终确定的热功当量为 1 cal = 4.184 J 1\ \text{cal} = 4.184\ \text{J}1 cal = 4 . 1 8 4 J ,与现代公认值一致。
能量守恒定律(热力学第一定律)
综合焦耳等人的工作,赫尔曼·冯·亥姆霍兹(Hermann von Helmholtz,1821-1894)在1847年系统阐述了能量守恒定律:在一个孤立系统中,各种形式的能量可以相互转化,但总能量保持不变。
数学表述为:
Δ U = Q − W \Delta U = Q - W
Δ U = Q − W
其中 Δ U \Delta UΔ U 是系统内能变化,Q QQ 是系统吸收的热量,W WW 是系统对外做的功。
一个具体计算:
1 kg 1\ \text{kg}1 kg 的 0 ∘ C 0^\circ\text{C}0 ∘ C 冰融化为 0 ∘ C 0^\circ\text{C}0 ∘ C 水,需要吸收多少热量?
冰的熔化潜热 λ 冰 = 3.34 × 1 0 5 J/kg \lambda_{\text{冰}} = 3.34 \times 10^5\ \text{J/kg}λ 冰 = 3 . 3 4 × 1 0 5 J/kg :
Q = m λ 冰 = 1 × 3.34 × 1 0 5 = 3.34 × 1 0 5 J Q = m \lambda_{\text{冰}} = 1 \times 3.34 \times 10^5 = 3.34 \times 10^5\ \text{J}
Q = m λ 冰 = 1 × 3 . 3 4 × 1 0 5 = 3 . 3 4 × 1 0 5 J
这个热量相当于将 3.34 × 1 0 5 / ( 1 × 9.8 ) ≈ 34000 kg 3.34 \times 10^5 / (1 \times 9.8) \approx 34000\ \text{kg}3 . 3 4 × 1 0 5 / ( 1 × 9 . 8 ) ≈ 3 4 0 0 0 kg 的重物提升1米所做的功。或者相当于一个 100 W 100\ \text{W}1 0 0 W 的灯泡点亮 334000 / 100 = 3340 s ≈ 56 min 334000/100 = 3340\ \text{s} \approx 56\ \text{min}3 3 4 0 0 0 / 1 0 0 = 3 3 4 0 s ≈ 5 6 min 所消耗的电能。
卡诺循环
尼古拉·卡诺(Nicolas Carnot,1796-1832)在1824年出版的《关于火的动力之思考》中,从理想热机模型出发,揭示了热机效率的基本极限。
卡诺热机由两个等温过程和两个绝热过程组成。其效率只取决于高温热源温度 T H T_HT H 和低温热源温度 T C T_CT C :
η Carnot = 1 − T C T H \eta_{\text{Carnot}} = 1 - \frac{T_C}{T_H}
η Carnot = 1 − T H T C
数值示例: 一个蒸汽机,锅炉温度为 T H = 500 K T_H = 500\ \text{K}T H = 5 0 0 K (约 22 7 ∘ C 227^\circ\text{C}2 2 7 ∘ C ),冷凝器温度为 T C = 300 K T_C = 300\ \text{K}T C = 3 0 0 K (约 2 7 ∘ C 27^\circ\text{C}2 7 ∘ C ),则最大可能效率:
η Carnot = 1 − 300 500 = 0.40 = 40 % \eta_{\text{Carnot}} = 1 - \frac{300}{500} = 0.40 = 40\%
η Carnot = 1 − 5 0 0 3 0 0 = 0 . 4 0 = 4 0 %
这说明了为什么即使实际的蒸汽机效率只有 15-20%,理论上限也仅为 40%。即使将锅炉温度提高到 1000 K 1000\ \text{K}1 0 0 0 K (72 7 ∘ C 727^\circ\text{C}7 2 7 ∘ C ),效率也只能达到:
η Carnot = 1 − 300 1000 = 0.70 = 70 % \eta_{\text{Carnot}} = 1 - \frac{300}{1000} = 0.70 = 70\%
η Carnot = 1 − 1 0 0 0 3 0 0 = 0 . 7 0 = 7 0 %
下表展示不同温差下的卡诺效率:
T H T_HT H (K)
T C T_CT C (K)
Δ T \Delta TΔ T (K)
效率 η \etaη
400
300
100
25.0%
500
300
200
40.0%
600
300
300
50.0%
800
300
500
62.5%
1000
300
700
70.0%
可见,要提高热机效率,提高热源温度远比降低冷源温度有效(因为 T C T_CT C 受环境温度制约,最低约为 300 K 300\ \text{K}3 0 0 K )。
熵的引入
鲁道夫·克劳修斯(Rudolf Clausius,1822-1888)在1865年定义了熵 (Entropy)的概念。对于可逆过程,系统的熵变为:
Δ S = ∫ d Q T \Delta S = \int \frac{dQ}{T}
Δ S = ∫ T d Q
热力学第二定律表述为:孤立系统的熵永不减少,即 Δ S ≥ 0 \Delta S \geq 0Δ S ≥ 0 。
一个具体数值计算:
将一块温度为 T hot = 350 K T_{\text{hot}} = 350\ \text{K}T hot = 3 5 0 K 、质量为 m = 1 kg m = 1\ \text{kg}m = 1 kg 的铁块(比热容 c_\textFe = 450\ \textJ/(kg·K)),放入温度为 T cool = 300 K T_{\text{cool}} = 300\ \text{K}T cool = 3 0 0 K 的大量水中。铁块最终温度变为 300 K 300\ \text{K}3 0 0 K 。
铁块的熵变(假设可逆冷却):
Δ S Fe = ∫ m c Fe d T T = m c Fe ln T final T initial \Delta S_{\text{Fe}} = \int \frac{m c_{\text{Fe}} dT}{T} = m c_{\text{Fe}} \ln \frac{T_{\text{final}}}{T_{\text{initial}}}
Δ S Fe = ∫ T m c Fe d T = m c Fe ln T initial T final
Δ S Fe = 1 × 450 × ln 300 350 = 450 × ln ( 0.857 ) ≈ 450 × ( − 0.154 ) = − 69.3 J/K \Delta S_{\text{Fe}} = 1 \times 450 \times \ln \frac{300}{350} = 450 \times \ln(0.857) \approx 450 \times (-0.154) = -69.3\ \text{J/K}
Δ S Fe = 1 × 4 5 0 × ln 3 5 0 3 0 0 = 4 5 0 × ln ( 0 . 8 5 7 ) ≈ 4 5 0 × ( − 0 . 1 5 4 ) = − 6 9 . 3 J/K
水吸收热量 Q = m c Fe ( 350 − 300 ) = 1 × 450 × 50 = 22500 J Q = m c_{\text{Fe}} (350-300) = 1 \times 450 \times 50 = 22500\ \text{J}Q = m c Fe ( 3 5 0 − 3 0 0 ) = 1 × 4 5 0 × 5 0 = 2 2 5 0 0 J ,在 300 K 300\ \text{K}3 0 0 K 下水的熵变为:
Δ S 水 = Q T = 22500 300 = 75.0 J/K \Delta S_{\text{水}} = \frac{Q}{T} = \frac{22500}{300} = 75.0\ \text{J/K}
Δ S 水 = T Q = 3 0 0 2 2 5 0 0 = 7 5 . 0 J/K
总熵变 Δ S 总 = − 69.3 + 75.0 = 5.7 J/K > 0 \Delta S_{\text{总}} = -69.3 + 75.0 = 5.7\ \text{J/K} > 0Δ S 总 = − 6 9 . 3 + 7 5 . 0 = 5 . 7 J/K > 0 ,符合热力学第二定律。
如果反过来,热量自发地从低温水传到高温铁块(使铁块升温),则总熵变为 − 5.7 J/K -5.7\ \text{J/K}− 5 . 7 J/K ,违反第二定律——这正是时间箭头的微观热力学基础。
路德维希·玻尔兹曼(Ludwig Boltzmann,1844-1906)将热力学现象还原为大量分子运动的统计行为。他的关键公式将熵与微观状态数联系起来:
S = k B ln Ω S = k_B \ln \Omega
S = k B ln Ω
其中 k B ≈ 1.38 × 1 0 − 23 J/K k_B \approx 1.38 \times 10^{-23}\ \text{J/K}k B ≈ 1 . 3 8 × 1 0 − 2 3 J/K 是玻尔兹曼常数,Ω \OmegaΩ 是系统的微观状态数。
这个公式的精髓在于:它将宏观量(熵)与微观量(状态的计数)联系起来了。玻尔兹曼公式被刻在他的墓碑上,作为对他伟大贡献的永久纪念。
一个简单计算:
一个系统有 N = 100 N = 100N = 1 0 0 个独立自旋 1 / 2 1/21 / 2 粒子,每个粒子有两种可能的状态(上或下)。如果总自旋向上数为 n nn ,则微观状态数(组合数):
Ω = N ! n ! ( N − n ) ! \Omega = \frac{N!}{n!(N-n)!}
Ω = n ! ( N − n ) ! N !
当 n = 50 n = 50n = 5 0 时:
Ω 50 = 100 ! 50 ! × 50 ! ≈ 1.01 × 1 0 29 \Omega_{50} = \frac{100!}{50! \times 50!} \approx 1.01 \times 10^{29}
Ω 5 0 = 5 0 ! × 5 0 ! 1 0 0 ! ≈ 1 . 0 1 × 1 0 2 9
对应的熵:
S 50 = k B ln ( 1.01 × 1 0 29 ) ≈ 1.38 × 1 0 − 23 × 66.8 ≈ 9.22 × 1 0 − 22 J/K S_{50} = k_B \ln(1.01 \times 10^{29}) \approx 1.38 \times 10^{-23} \times 66.8 \approx 9.22 \times 10^{-22}\ \text{J/K}
S 5 0 = k B ln ( 1 . 0 1 × 1 0 2 9 ) ≈ 1 . 3 8 × 1 0 − 2 3 × 6 6 . 8 ≈ 9 . 2 2 × 1 0 − 2 2 J/K
当 n = 100 n = 100n = 1 0 0 (全部向上)时,Ω = 1 \Omega = 1Ω = 1 ,S = 0 S = 0S = 0 。
当 n = 60 n = 60n = 6 0 时:
Ω 60 = 100 ! 60 ! × 40 ! ≈ 1.37 × 1 0 28 \Omega_{60} = \frac{100!}{60! \times 40!} \approx 1.37 \times 10^{28}
Ω 6 0 = 6 0 ! × 4 0 ! 1 0 0 ! ≈ 1 . 3 7 × 1 0 2 8
可见 n = 50 n=50n = 5 0 的状态数远超 n = 60 n=60n = 6 0 ——这正是为什么系统会趋向于均匀分布(最大熵状态)的统计力学解释。
吉尔伯特的开创
电磁学的科学开端可以追溯到威廉·吉尔伯特(William Gilbert,1544-1603)。他在1600年出版的《论磁》(De Magnete )中,通过系统的实验区分了电和磁现象。他发现摩擦琥珀能吸引轻小物体(这就是"电"一词的来源——希腊语 ελεκτρον,意为琥珀),而天然磁石能吸引铁。
库仑定律
查尔斯-奥古斯丁·库仑(Charles-Augustin de Coulomb,1736-1806)在1785年用扭秤实验精确测定了静电力与距离的关系。他的扭秤装置如下:一根细银丝悬挂一个水平绝缘杆,杆的一端带有一个带电小球,另一端是平衡锤。通过测量银丝的扭转角来测定静电力。
库仑发现两个点电荷之间的静电力与它们电量的乘积成正比,与距离的平方成反比:
F ⃗ 12 = k q 1 q 2 r 2 r ^ 12 \vec{F}_{12} = k \frac{q_1 q_2}{r^2} \hat{r}_{12}
F 1 2 = k r 2 q 1 q 2 r ^ 1 2
其中 k \approx 8.99 \times 10^9\ \textN·m^2/\textC^2(更常用的是真空介电常数 ε 0 = 1 / ( 4 π k ) ≈ 8.854 × 1 0 − 12 F/m \varepsilon_0 = 1/(4\pi k) \approx 8.854 \times 10^{-12}\ \text{F/m}ε 0 = 1 / ( 4 π k ) ≈ 8 . 8 5 4 × 1 0 − 1 2 F/m ,使得 k = 1 / ( 4 π ε 0 ) k = 1/(4\pi \varepsilon_0)k = 1 / ( 4 π ε 0 ) )。
数值示例: 两个电量分别为 q 1 = + 2 μ C q_1 = +2\ \mu\text{C}q 1 = + 2 μ C 和 q 2 = + 3 μ C q_2 = +3\ \mu\text{C}q 2 = + 3 μ C 的小球,相距 r = 0.5 m r = 0.5\ \text{m}r = 0 . 5 m :
F = 8.99 × 1 0 9 × 2 × 1 0 − 6 × 3 × 1 0 − 6 0. 5 2 F = 8.99 \times 10^9 \times \frac{2 \times 10^{-6} \times 3 \times 10^{-6}}{0.5^2}
F = 8 . 9 9 × 1 0 9 × 0 . 5 2 2 × 1 0 − 6 × 3 × 1 0 − 6
F = 8.99 × 1 0 9 × 6 × 1 0 − 12 0.25 = 8.99 × 1 0 9 × 2.4 × 1 0 − 11 F = 8.99 \times 10^9 \times \frac{6 \times 10^{-12}}{0.25} = 8.99 \times 10^9 \times 2.4 \times 10^{-11}
F = 8 . 9 9 × 1 0 9 × 0 . 2 5 6 × 1 0 − 1 2 = 8 . 9 9 × 1 0 9 × 2 . 4 × 1 0 − 1 1
F = 0.216 N F = 0.216\ \text{N}
F = 0 . 2 1 6 N
这个力约相当于 22 g 22\ \text{g}2 2 g 物体的重力——对于静电力来说相当可观。
如果两个电荷异号,库仑力变为吸引力;如果同号,则为排斥力。注意,库仑定律的形式与万有引力定律完全相同(都与距离平方成反比),只是系数和参与物理量不同(电荷 vs. 质量)。
静磁学
静磁学的发展比静电学更早,因为天然磁石的存在早已被人类认识。然而直到1819年汉斯·克里斯蒂安·奥斯特(Hans Christian Oersted,1777-1851)发现电流的磁效应后,电磁学的真正突破才开始。
奥斯特的实验(1820年)
奥斯特在哥本哈根大学的课堂上做了一个简单的实验:将一根导线平行于指南针放置,当给导线通电时,指南针发生了偏转。这个看似简单的现象,其革命性在于——它第一次证明了电和磁不是相互独立的,而是有内在联系的。
实验数据:当通过 1 A 1\ \text{A}1 A 的电流时,距离导线 r = 0.1 m r = 0.1\ \text{m}r = 0 . 1 m 处的磁感应强度为:
B = μ 0 I 2 π r = 4 π × 1 0 − 7 × 1 2 π × 0.1 = 2 × 1 0 − 6 T B = \frac{\mu_0 I}{2\pi r} = \frac{4\pi \times 10^{-7} \times 1}{2\pi \times 0.1} = 2 \times 10^{-6}\ \text{T}
B = 2 π r μ 0 I = 2 π × 0 . 1 4 π × 1 0 − 7 × 1 = 2 × 1 0 − 6 T
地磁场约为 5 × 1 0 − 5 T 5 \times 10^{-5}\ \text{T}5 × 1 0 − 5 T ,因此上述磁场约为地磁场的 2 / 50 = 4 % 2/50 = 4\%2 / 5 0 = 4 % ,足以使指南针发生可观测的偏转。
安培定律
安德烈-玛丽·安培(André-Marie Ampère,1775-1836)在奥斯特实验的启发下,迅速建立了两条平行载流导线之间相互作用力的数学理论。他发现,两条平行载流导线之间的单位长度作用力为:
F l = μ 0 I 1 I 2 2 π d \frac{F}{l} = \frac{\mu_0 I_1 I_2}{2\pi d}
l F = 2 π d μ 0 I 1 I 2
其中 d dd 是导线间距。
数值示例: 两条相距 d = 0.2 m d = 0.2\ \text{m}d = 0 . 2 m 的平行导线,分别通有 I 1 = 10 A I_1 = 10\ \text{A}I 1 = 1 0 A 和 I 2 = 15 A I_2 = 15\ \text{A}I 2 = 1 5 A 的同向电流,则单位长度的吸引力:
F l = 4 π × 1 0 − 7 × 10 × 15 2 π × 0.2 = 4 π × 1 0 − 7 × 150 2 π × 0.2 \frac{F}{l} = \frac{4\pi \times 10^{-7} \times 10 \times 15}{2\pi \times 0.2} = \frac{4\pi \times 10^{-7} \times 150}{2\pi \times 0.2}
l F = 2 π × 0 . 2 4 π × 1 0 − 7 × 1 0 × 1 5 = 2 π × 0 . 2 4 π × 1 0 − 7 × 1 5 0
= 4 × 1 0 − 7 × 150 2 × 0.2 = 6 × 1 0 − 5 0.4 = 1.5 × 1 0 − 4 N/m = \frac{4 \times 10^{-7} \times 150}{2 \times 0.2} = \frac{6 \times 10^{-5}}{0.4} = 1.5 \times 10^{-4}\ \text{N/m}
= 2 × 0 . 2 4 × 1 0 − 7 × 1 5 0 = 0 . 4 6 × 1 0 − 5 = 1 . 5 × 1 0 − 4 N/m
这个力相当微弱,但在工程中可以叠加——例如一个1000匝的线圈可以放大这个效应。
迈克尔·法拉第(Michael Faraday,1791-1867)是最伟大的实验物理学家之一。虽然他的数学知识有限,但他凭借天才的实验直觉,发现了电磁感应的核心规律,并引入了场 (Field)的概念——这可能是经典物理中最重要的理论创新之一。
电磁感应(1831年)
法拉第发现:当一个磁铁在线圈中运动时,线圈中会产生电流。更精确地说,通过改变穿过一个回路的磁通量,会在回路中产生感应电动势:
E = − d Φ B d t \mathcal{E} = -\frac{d\Phi_B}{dt}
E = − d t d Φ B
其中 Φ B = ∫ B ⃗ ⋅ d A ⃗ \Phi_B = \int \vec{B} \cdot d\vec{A}Φ B = ∫ B ⋅ d A 为磁通量。
数值示例: 一个圆形线圈(半径 r = 0.1 m r = 0.1\ \text{m}r = 0 . 1 m ,匝数 N = 50 N = 50N = 5 0 ),置于均匀磁场 B = 0.5 T B = 0.5\ \text{T}B = 0 . 5 T 中。在 Δ t = 0.2 s \Delta t = 0.2\ \text{s}Δ t = 0 . 2 s 内,磁场均匀减小到零。
初始磁通量(磁场垂直于线圈平面):Φ i = B × π r 2 = 0.5 × π × 0.01 ≈ 0.0157 Wb \Phi_i = B \times \pi r^2 = 0.5 \times \pi \times 0.01 \approx 0.0157\ \text{Wb}Φ i = B × π r 2 = 0 . 5 × π × 0 . 0 1 ≈ 0 . 0 1 5 7 Wb
末态磁通量:Φ f = 0 \Phi_f = 0Φ f = 0
感应电动势:
E = − N Δ Φ Δ t = − 50 × 0 − 0.0157 0.2 = 50 × 0.0785 = 3.93 V \mathcal{E} = -N \frac{\Delta \Phi}{\Delta t} = -50 \times \frac{0 - 0.0157}{0.2} = 50 \times 0.0785 = 3.93\ \text{V}
E = − N Δ t Δ Φ = − 5 0 × 0 . 2 0 − 0 . 0 1 5 7 = 5 0 × 0 . 0 7 8 5 = 3 . 9 3 V
如果线圈电阻为 R = 2 Ω R = 2\ \OmegaR = 2 Ω ,则感应电流:
I = E R = 3.93 2 = 1.97 A I = \frac{\mathcal{E}}{R} = \frac{3.93}{2} = 1.97\ \text{A}
I = R E = 2 3 . 9 3 = 1 . 9 7 A
法拉第电磁感应定律是发电机的理论基础——机械能转化为电能的核心机制。今天全球99%以上的电力生产都依赖法拉第的发现(通过涡轮机带动磁铁在线圈中旋转)。下表为不同磁场变化率下的感应电动势比较:
线圈匝数 N NN
面积 A AA (m 2 \text{m}^2m 2 )
磁场变化率 Δ B / Δ t \Delta B/\Delta tΔ B / Δ t (T/s)
感应电动势 E \mathcal{E}E (V)
50
0.0314
2.5
3.93
100
0.0314
2.5
7.86
50
0.0314
5.0
7.86
100
0.0314
5.0
15.7
力线(场)的概念
法拉第引入了"力线"(Lines of Force)这一直观概念来描述电场和磁场。他指出:
电荷周围有电力线,从正电荷出发,终止于负电荷
磁力线是闭合的,没有起点和终点
力线的密度表示场的强弱
尽管当时的许多数学物理学家对法拉第缺乏数学表述的做法不以为然,但他的"场"概念后来被麦克斯韦数学化,并成为整个物理学的核心范式——现代物理学中的规范场论直接继承了"场"的概念。
詹姆斯·克拉克·麦克斯韦(James Clerk Maxwell,1831-1879)被普遍认为是牛顿之后最伟大的理论物理学家。他在《电磁通论》(A Treatise on Electricity and Magnetism , 1873)中,将电磁学的一切规律总结为4个方程,并引入了位移电流 (Displacement Current)这一关键概念,完成了电磁理论的统一。
麦克斯韦方程组(积分形式)
∮ E ⃗ ⋅ d A ⃗ = Q ε 0 (高斯定律) \oint \vec{E} \cdot d\vec{A} = \frac{Q}{\varepsilon_0} \quad \text{(高斯定律)}
∮ E ⋅ d A = ε 0 Q ( 高斯定律 )
∮ B ⃗ ⋅ d A ⃗ = 0 (磁高斯定律) \oint \vec{B} \cdot d\vec{A} = 0 \quad \text{(磁高斯定律)}
∮ B ⋅ d A = 0 ( 磁高斯定律 )
∮ E ⃗ ⋅ d l ⃗ = − d Φ B d t (法拉第定律) \oint \vec{E} \cdot d\vec{l} = -\frac{d\Phi_B}{dt} \quad \text{(法拉第定律)}
∮ E ⋅ d l = − d t d Φ B ( 法拉第定律 )
∮ B ⃗ ⋅ d l ⃗ = μ 0 I + μ 0 ε 0 d Φ E d t (安培-麦克斯韦定律) \oint \vec{B} \cdot d\vec{l} = \mu_0 I + \mu_0 \varepsilon_0 \frac{d\Phi_E}{dt} \quad \text{(安培-麦克斯韦定律)}
∮ B ⋅ d l = μ 0 I + μ 0 ε 0 d t d Φ E ( 安培 - 麦克斯韦定律 )
麦克斯韦方程组(微分形式)
∇ ⋅ E ⃗ = ρ ε 0 \nabla \cdot \vec{E} = \frac{\rho}{\varepsilon_0}
∇ ⋅ E = ε 0 ρ
∇ ⋅ B ⃗ = 0 \nabla \cdot \vec{B} = 0
∇ ⋅ B = 0
∇ × E ⃗ = − ∂ B ⃗ ∂ t \nabla \times \vec{E} = -\frac{\partial \vec{B}}{\partial t}
∇ × E = − ∂ t ∂ B
∇ × B ⃗ = μ 0 J ⃗ + μ 0 ε 0 ∂ E ⃗ ∂ t \nabla \times \vec{B} = \mu_0 \vec{J} + \mu_0 \varepsilon_0 \frac{\partial \vec{E}}{\partial t}
∇ × B = μ 0 J + μ 0 ε 0 ∂ t ∂ E
每个方程的物理含义如下表:
方程
含义
关键物理洞察
∇ ⋅ E ⃗ = ρ / ε 0 \nabla \cdot \vec{E} = \rho/\varepsilon_0∇ ⋅ E = ρ / ε 0
电荷是电场的源
电场线从正电荷出发
∇ ⋅ B ⃗ = 0 \nabla \cdot \vec{B} = 0∇ ⋅ B = 0
无磁单极子
磁力线始终闭合
∇ × E ⃗ = − ∂ B ⃗ / ∂ t \nabla \times \vec{E} = -\partial \vec{B}/\partial t∇ × E = − ∂ B / ∂ t
变化的磁场产生电场
电磁感应的核心
∇ × B ⃗ = μ 0 J ⃗ + μ 0 ε 0 ∂ E ⃗ / ∂ t \nabla \times \vec{B} = \mu_0 \vec{J} + \mu_0 \varepsilon_0 \partial \vec{E}/\partial t∇ × B = μ 0 J + μ 0 ε 0 ∂ E / ∂ t
电流和变化电场产生磁场
位移电流是关键修正
位移电流的预测性贡献
麦克斯韦在第四方程中增加 μ 0 ε 0 ∂ E ⃗ / ∂ t \mu_0 \varepsilon_0 \partial \vec{E}/\partial tμ 0 ε 0 ∂ E / ∂ t 项,最初完全出于理论完备性的考虑——否则在电容器充电过程中,安培环路定律会出现矛盾(电流在电容器极板间中断)。然而这个修正带来了一个极其重要的预测:电磁波的存在。
从麦克斯韦方程组可以推导出电磁波方程。在真空中(无电荷和电流,即 ρ = 0 , J ⃗ = 0 \rho = 0, \vec{J} = 0ρ = 0 , J = 0 ),电场 E ⃗ \vec{E}E 满足:
∇ 2 E ⃗ = μ 0 ε 0 ∂ 2 E ⃗ ∂ t 2 \nabla^2 \vec{E} = \mu_0 \varepsilon_0 \frac{\partial^2 \vec{E}}{\partial t^2}
∇ 2 E = μ 0 ε 0 ∂ t 2 ∂ 2 E
这是一个经典的波动方程,波速为:
c = 1 μ 0 ε 0 c = \frac{1}{\sqrt{\mu_0 \varepsilon_0}}
c = μ 0 ε 0 1
将 μ 0 = 4 π × 1 0 − 7 H/m \mu_0 = 4\pi \times 10^{-7}\ \text{H/m}μ 0 = 4 π × 1 0 − 7 H/m 和 ε 0 ≈ 8.854 × 1 0 − 12 F/m \varepsilon_0 \approx 8.854 \times 10^{-12}\ \text{F/m}ε 0 ≈ 8 . 8 5 4 × 1 0 − 1 2 F/m 代入:
c = 1 4 π × 1 0 − 7 × 8.854 × 1 0 − 12 ≈ 2.998 × 1 0 8 m/s c = \frac{1}{\sqrt{4\pi \times 10^{-7} \times 8.854 \times 10^{-12}}} \approx 2.998 \times 10^8\ \text{m/s}
c = 4 π × 1 0 − 7 × 8 . 8 5 4 × 1 0 − 1 2 1 ≈ 2 . 9 9 8 × 1 0 8 m/s
巧合的是,这个数值与当时已知的光速(斐索1849年测量值 3.15 × 1 0 8 m/s 3.15 \times 10^8\ \text{m/s}3 . 1 5 × 1 0 8 m/s )惊人地一致。麦克斯韦由此做出了物理学史上最伟大的理论预测之一:光就是电磁波 。
赫兹实验(1887-1888)
海因里希·赫兹(Heinrich Hertz,1857-1894)在卡尔斯鲁厄大学设计了精密的电磁波产生和探测装置,以实验验证麦克斯韦的预言。
赫兹的发射器是一个偶极天线——两根金属杆,末端各有一个金属球,中间有微小间隙。当给天线加上高压时,间隙中产生火花放电(电振荡),从而辐射电磁波。接收器是一个类似的环形天线,在间隙处产生微小火花作为探测信号。
赫兹的关键实验数据:
发射电磁波频率:约 5 × 1 0 7 Hz 5 \times 10^7\ \text{Hz}5 × 1 0 7 Hz (50 MHz)
波长:约 6 m 6\ \text{m}6 m
测量波速:v = f λ = 5 × 1 0 7 × 6 = 3 × 1 0 8 m/s v = f \lambda = 5 \times 10^7 \times 6 = 3 \times 10^8\ \text{m/s}v = f λ = 5 × 1 0 7 × 6 = 3 × 1 0 8 m/s
该值与麦克斯韦预测的光速一致
赫兹还证明了电磁波具有与光相同的性质:
反射 :用金属板反射电磁波,检测到驻波
折射 :用沥青棱镜使电磁波发生偏折
偏振 :用金属线栅使电磁波偏振
衍射 :电磁波能绕过障碍物传播
麦克斯韦预言的实验验证时间线:
年份
贡献者
验证内容
1865
麦克斯韦
理论预言电磁波存在
1887
赫兹
产生并探测到电磁波
1888
赫兹
证明电磁波具有反射、折射等光学特性
1895
马可尼
远距离无线电报通信(超过1英里)
1901
马可尼
跨大西洋无线电报(英国-纽芬兰,约3500公里)
赫兹实验不仅验证了麦克斯韦理论,更开启了无线通信时代。当被问及其发现的实际用途时,赫兹回答:"它毫无用处——这只是一个实验,证明麦克斯韦大师是对的。我们只是拥有了神秘波,用肉眼看不见,但它确实存在。"——然而不到10年,马可尼就将电磁波用于商业通信。
经典力学最深远的思想遗产之一是拉普拉斯的决定论。它预设宇宙是严格因果律支配的:只要知道宇宙中所有粒子在某一时刻的位置和速度,就能推算出过去和未来的全部演化。
这种世界观在19世纪达到了巅峰。皮埃尔-西蒙·拉普拉斯(Pierre-Simon Laplace)在1814年写道:
"我们可以把宇宙的现状视为其过去的结果和未来的原因。如果有一个智能体,它知道所有使自然界运动的力和所有组成它的物体的位置……那么对它来说,没有任何东西是不确定的,未来就像过去一样呈现在它面前。"
然而20世纪的发展表明,经典决定论存在两个根本局限:
量子力学的不确定性原理 :无法同时精确知道粒子的位置和动量
混沌理论对初值的敏感依赖性 :即使经典力学本身,非线性系统的长时间演化也无法准确预测
今天我们知道,经典物理学的有效性受到两个边界的约束:
区域
适用理论
适用条件
宏观低速(v ≪ c v \ll cv ≪ c )
经典力学/牛顿力学
日常世界
宏观高速(v ∼ c v \sim cv ∼ c )
狭义/广义相对论
粒子加速器、天文现象
微观低速
量子力学
原子、分子
微观高速
量子场论
高能物理、粒子物理
经典力学在以下情况仍然精确有效:
物体速度小于光速的 1 % 1\%1 % (即 v < 3 × 1 0 6 m/s v < 3 \times 10^6\ \text{m/s}v < 3 × 1 0 6 m/s )
物体尺度大于原子尺度(r > 1 0 − 10 m r > 10^{-10}\ \text{m}r > 1 0 − 1 0 m )
引力场较弱(牛顿引力近似)
一个典型对比:
计算一个 m = 1 kg m = 1\ \text{kg}m = 1 kg 的物体在 F = 2 N F = 2\ \text{N}F = 2 N 恒力作用下加速 t = 10 s t = 10\ \text{s}t = 1 0 s 后的速度。牛顿力学给出:v = F t / m = 2 × 10 / 1 = 20 m/s v = Ft/m = 2 \times 10/1 = 20\ \text{m/s}v = F t / m = 2 × 1 0 / 1 = 2 0 m/s 。相对论修正因子 γ = 1 / 1 − ( v / c ) 2 ≈ 1 / 1 − ( 20 / 3 × 1 0 8 ) 2 ≈ 1 + 2.22 × 1 0 − 15 \gamma = 1/\sqrt{1 - (v/c)^2} \approx 1/\sqrt{1 - (20/3\times 10^8)^2} \approx 1 + 2.22 \times 10^{-15}γ = 1 / 1 − ( v / c ) 2 ≈ 1 / 1 − ( 2 0 / 3 × 1 0 8 ) 2 ≈ 1 + 2 . 2 2 × 1 0 − 1 5 ,修正量约为 1 0 − 15 10^{-15}1 0 − 1 5 量级,完全可忽略。
人物
生卒年份
核心贡献
名言/轶事
伽利略
1564-1642
实验科学方法、惯性原理、自由落体
"科学的真理不应建立在权威之上"
开普勒
1571-1630
行星运动三定律
"我宁愿成为一个尖锐的批评家"
牛顿
1643-1727
三大定律、万有引力、微积分
"如果我看得更远,是因为我站在巨人的肩膀上"
拉格朗日
1736-1813
分析力学、拉格朗日量
"拉格朗日是数学诗人在力学中的化身"
卡诺
1796-1832
热机理论、卡诺循环
享年36岁,其著作被长期忽视
法拉第
1791-1867
电磁感应、场的概念
仅读过小学,但被爱因斯坦称为物理学史上最具独创性的人物之一
焦耳
1818-1889
热功当量测定、能量守恒
啤酒酿造商出身,业余科学家
麦克斯韦
1831-1879
麦克斯韦方程组、电磁波预言
"电磁理论之父",享年仅48岁
玻尔兹曼
1844-1906
统计力学、玻尔兹曼熵公式
因学术争论导致抑郁,最终自杀
赫兹
1857-1894
电磁波实验验证
享年仅36岁
经典物理学300年的发展可以用一条清晰的逻辑链概括:
伽利略 建立实验方法 → 发现自由落体规律
开普勒 精确描述行星轨道 → 三定律
牛顿 统一天地运动 → 三大定律 + 万有引力
拉格朗日/哈密顿 实现力学数学化 → 分析力学
卡诺/焦耳/克劳修斯 建立热力学 → 能量守恒 + 熵
玻尔兹曼 从微观解释宏观 → 统计力学
库仑/奥斯特/法拉第 积累电磁知识 → 电学、磁学、电磁感应
麦克斯韦 统一电、磁、光 → 电磁波理论
赫兹 实验验证 → 经典物理范式完成
这九大步骤环环相扣,每一步都建立在前人的基础上。经典物理学的伟大在于它不仅仅是知识的总和,更是一套范式的确立:用数学精确描述自然,用实验严格验证理论 ——这套范式至今仍是所有自然科学的方法论基石。
正如爱因斯坦所说:"纯粹的逻辑思维不能给我们任何关于经验世界的知识;所有关于实在的知识,都从经验开始,又终于经验。"
Newton, I. (1687). Philosophiæ Naturalis Principia Mathematica . 中译本:《自然哲学的数学原理》,商务印书馆。
Maxwell, J. C. (1873). A Treatise on Electricity and Magnetism . Oxford: Clarendon Press.
Feynman, R. P., Leighton, R. B., & Sands, M. (1963). The Feynman Lectures on Physics . Addison-Wesley.
Kuhn, T. S. (1962). The Structure of Scientific Revolutions . University of Chicago Press. 中译本:《科学革命的结构》,北京大学出版社。
Goldstein, H., Poole, C., & Safko, J. (2001). Classical Mechanics (3rd ed.). Addison-Wesley.
Jackson, J. D. (1998). Classical Electrodynamics (3rd ed.). Wiley.
吴大猷. (1975). 《理论物理》第一册:古典力学. 联经出版.
赵凯华, 陈熙谋. (2003). 《电磁学》(第三版). 高等教育出版社.