万有引力与航天
本课从开普勒三定律出发，推导万有引力定律，并深入探讨天体运动、卫星轨道、宇宙速度等航天知识。通过矢量分析与几何图示，攻克双星系统、变轨能量等难点，培养建模与计算能力。
1 课堂导入：从苹果到行星
牛顿坐在苹果树下，被掉落的苹果启发了万有引力的存在。但更令人惊叹的是，他意识到月亮也如同苹果一般，被地球“拉向”地面——只是因为月亮的水平速度足够大，才一直坠落而不落地，形成了绕地球的圆周运动。正是这种“天上”与“人间”统一的引力思想，奠定了近代天体力学的基础。今天我们就沿着牛顿的足迹，系统学习万有引力与航天知识。
2 开普勒三定律——行星运动的几何规律
在牛顿之前，第谷·布拉赫积累了大量精确的行星观测数据，他的助手开普勒经过多年分析，总结出三条定律。
2.1 第一定律（椭圆定律）
所有行星绕太阳运动的轨道都是椭圆，太阳处于椭圆的一个焦点上。
2.2 第二定律（面积定律）
对任意一个行星来说，它与太阳的连线在相等的时间内扫过相等的面积。这意味着行星在近日点速度大，远日点速度小。
2.3 第三定律（周期定律）
所有行星轨道的半长轴的三次方与公转周期的平方的比值都相等，即 $\dfrac{a^3}{T^2} = k$，其中 $k$ 与中心天体有关。
对于太阳系，$k = \dfrac{GM_\text{太阳}}{4\pi^2}$。该定律不仅适用于椭圆轨道，也适用于圆轨道（$a$ 即半径）。
3 万有引力定律
3.1 定律内容与公式
自然界中任何两个物体都相互吸引，引力的方向在它们的连线上，大小与它们的质量乘积成正比，与它们距离的平方成反比：
$$F = G\frac{m_1 m_2}{r^2}$$其中 $G = 6.67 \times 10^{-11} \ \mathrm{N \cdot m^2 / kg^2}$，称为引力常量，由卡文迪许通过扭秤实验测得。
3.2 万有引力与重力的区别（困难点）
在地球表面，物体所受万有引力 $F_\text{引}$ 有两个效果：一是提供随地球自转所需的向心力 $F_\text{向}$，二是产生重力 $mg$。注意：$g$ 是重力加速度，并非引力加速度。
关系：$G\dfrac{Mm}{R^2} = mg + m\omega^2 R\cos^2\varphi$ （矢量式，向心力 $m\omega^2 r_{\perp}$）。
在两级处，向心力为零，$g_{\text{极}} = \dfrac{GM}{R^2}$；在赤道处，$g_{\text{赤}} = \dfrac{GM}{R^2} - \omega^2 R$。
4 天体运动的核心应用
4.1 计算中心天体质量
对于绕中心天体做圆周运动的卫星，由 $G\dfrac{Mm}{r^2} = m\dfrac{v^2}{r} = m\omega^2 r = m\dfrac{4\pi^2}{T^2}r$，可得
$$M = \frac{v^2 r}{G} = \frac{\omega^2 r^3}{G} = \frac{4\pi^2 r^3}{GT^2}$$其中 $r$ 为卫星轨道半径，$T$ 为周期。如果不考虑卫星质量，$M$ 为中心天体质量。
4.2 双星系统（困难）
两个质量分别为 $m_1$、$m_2$ 的星体，绕它们连线上的一点（质心）做圆周运动，距离恒定 $L$。
关键关系：
- $m_1 r_1 = m_2 r_2$（质心位置）
- $\dfrac{Gm_1 m_2}{L^2} = m_1 \omega^2 r_1 = m_2 \omega^2 r_2$
- $\omega = \sqrt{\dfrac{G(m_1+m_2)}{L^3}}$，周期 $T = 2\pi \sqrt{\dfrac{L^3}{G(m_1+m_2)}}$
5 宇宙速度与卫星轨道
5.1 三个宇宙速度
- 第一宇宙速度（环绕速度）：$v_1 = \sqrt{\dfrac{GM}{R}} \approx 7.9 \ \text{km/s}$，是物体在地面附近绕地球做匀速圆周运动的速度。
- 第二宇宙速度（脱离速度）：$v_2 = \sqrt{\dfrac{2GM}{R}} \approx 11.2 \ \text{km/s}$，使物体脱离地球引力束缚的最小速度。
- 第三宇宙速度（逃逸太阳系）：$v_3 \approx 16.7 \ \text{km/s}$（需借助地球公转速度）。
5.2 卫星变轨与能量变化（困难）
从低轨道到高轨道需要加速（增加动能），卫星进入椭圆转移轨道，在远地点再次加速进入高圆轨道。变轨过程中引力做负功，总机械能增加。注意：高轨道运行速度反而更小（$v = \sqrt{GM/r}$）。
6 典型例题
例题1（中等难度）
已知地球半径为 R，地球表面重力加速度为 g，不考虑地球自转。求第一宇宙速度的表达式（用 R、g 表示）。
思路：第一宇宙速度即近地卫星的环绕速度，万有引力提供向心力，而地面附近 $mg \approx G\dfrac{Mm}{R^2}$，可消去 GM。
步骤：
- 由重力近似等于万有引力：$mg = G\dfrac{Mm}{R^2} \Rightarrow GM = gR^2$
- 对近地卫星：$G\dfrac{Mm}{R^2} = m\dfrac{v^2}{R}$
- 代入得 $\dfrac{gR^2}{R} = v^2$，所以 $v_1 = \sqrt{gR}$
易错点：不要直接认为 $v_1 = \sqrt{\dfrac{GM}{R}}$ 想当然，要善于利用已知条件转化。
例题2（较难）
一双星系统，质量分别为 m 和 2m，相距 L。求 (1) 它们的轨道半径之比；(2) 公转周期。
思路：利用质心位置和向心力公式。
解答：
- 设质心到质量为 m 的星的距离为 r₁，到 2m 的为 r₂，则 $r_1 + r_2 = L$，且 $m r_1 = 2m r_2$，得 $r_1 : r_2 = 2:1$，故 $r_1 = \dfrac{2}{3}L$，$r_2 = \dfrac{1}{3}L$。
- 对于质量为 m 的星，$G\dfrac{m \cdot 2m}{L^2} = m \omega^2 r_1$，代入 r₁ 得 $\dfrac{2Gm}{L^2} = \omega^2 \cdot \dfrac{2}{3}L$，故 $\omega = \sqrt{\dfrac{3Gm}{L^3}}$，周期 $T = \dfrac{2\pi}{\omega} = 2\pi \sqrt{\dfrac{L^3}{3Gm}}$。
注意：解题时务必写清对哪个星体列方程，且引力距离是 L 而不是轨道半径。
例题3（困难）
某宇航员登上一未知星球，在星球表面用弹簧测力计测得质量为 m 的物体重力为 F。他乘坐探测器在距该星球表面 h 高处绕星做圆周运动，周期为 T。星球半径为 R，且已知引力常量 G。求该星球的质量和平均密度。
思路：利用两个条件：表面重力与引力关系；高空轨道周期与半径关系。
步骤：
- 表面：$F = G\dfrac{Mm}{R^2}$ ①
- 轨道：$G\dfrac{Mm}{(R+h)^2} = m\left(\dfrac{2\pi}{T}\right)^2 (R+h)$ ②
- 由①得 $M = \dfrac{F R^2}{G m}$，代入②可检验，但题目要求求出 M，故直接代入。
- 计算密度：$\rho = \dfrac{M}{\frac{4}{3}\pi R^3} = \dfrac{F R^2}{G m} \cdot \dfrac{3}{4\pi R^3} = \dfrac{3F}{4\pi G m R}$
易漏：平均密度公式中体积用 $ rac{4}{3}\pi R^3$，而非轨道半径。(R+h) 只用于求周期关系检验，最终密度表达式中仅含 R。
7 常见误区
- 混淆引力常量 G 与重力加速度 g：G 是普适常量，g 因星球而异。
- 忽略地球自转：计算地面物体重力时，若题目明确“不考虑自转”，才可近似 $mg = GMm/R^2$；否则要考虑向心力。
- 卫星变轨速度判断失误：从低轨到高轨需要两次加速，但最终在高轨运行速度反而变小。
- 双星问题误判半径：两星轨道半径和等于间距，各自半径与质量成反比。
8 学习建议
- 画图分析：将天体运动画成几何图，标注力、半径、速度、周期，把抽象问题可视化。
- 记熟“黄金代换”：$GM = gR^2$，解决很多地表问题。
- 善用比较：对比卫星在不同轨道上的 $v$、$\omega$、$a$、$T$、机械能。
- 强化能量观点：卫星轨道变化由机械能变化引起，可以结合动能定理、引力势能公式 $E_p = -\dfrac{GMm}{r}$ 来理解。
9 知识链接
- 牛顿力学：万有引力定律是牛顿第二定律在天体上的延伸。
- 圆周运动：向心力 $F = m\omega^2 r$ 是分析轨道运动的基础。
- 功与能：卫星变轨伴随机械能变化，需要能量守恒与功能关系。
- 航天常识：了解我国“神舟”“天宫”系列任务中的轨道变化，增加学习兴趣。
10 习题自测
- （基础）设地球表面重力加速度为 $g_0$，地球半径为 $R$，求距地面 $2R$ 高处的重力加速度 $g$ 是 $g_0$ 的多少倍？
- （提高）某行星的半径是地球半径的 2 倍，质量是地球质量的 8 倍，则该行星的第一宇宙速度是地球第一宇宙速度的多少倍？
- （挑战）两个星球组成双星系统，它们之间的引力大小为 $F$，运行的角速度为 $\omega$，周期为 $T$。若两星球的质量和距离保持不变，而将它们间的引力突然变为 $4F$（假设系统稳定），则新的角速度变为多少？
参考答案：
- $\dfrac{1}{9}$（由 $g = \dfrac{GM}{(R+h)^2}$，$h=2R$ 时 $r=3R$，$g = \dfrac{GM}{9R^2} = \dfrac{g_0}{9}$）
- 2 倍（$v_1 = \sqrt{\dfrac{GM}{R}}$，代入 $M'=8M_{\oplus}$，$R'=2R_{\oplus}$，得 $v_1' = \sqrt{\dfrac{G\cdot 8M_{\oplus}}{2R_{\oplus}}} = 2\sqrt{\dfrac{GM_{\oplus}}{R_{\oplus}}}=2v_1$）
- $2\omega$（原引力 $F = G\dfrac{m_1 m_2}{L^2}$，引力变为 $4F$，但本质上引力由质量距离决定，所以不可能单独变；这里假设系统稳定即瞬间质量或距离改变，但题意模糊。更严谨：双星系统中 $F = m_1 \omega^2 r_1 = m_2 \omega^2 r_2$，且 $m_1 r_1 = m_2 r_2$，$F = \dfrac{Gm_1 m_2}{L^2}$。若 $F' = 4F$，则 $\omega'^2 r_1' = 4\omega^2 r_1$，但距离 L 不变时 $r_1$ 也不会变（质心不变），故 $\omega' = 2\omega$）