4.1 自由飞行动力学：真空抛物线与真实轨迹
4.1 自由飞行动力学：真空抛物线与真实轨迹
本节摘要：自由飞行研究的起点是真空抛物线：射程 R = v²·sin2θ/g，45° 射角最远。这条公式在 850 m/s 初速下给出 70 多公里的"虚假繁荣"，而真实射程往往只有它的四分之一——差额全部记在空气阻力账上。本节先算清真空账，再把阻力加回来，给出点质量弹道的数值解法与真实轨迹的三段结构。这是第 3 章初始条件的第一次全链条使用。
位置感：内弹道给了初速，中间弹道给了初始姿态，本节开始消费这些输入。我们先故意忽略空气——不是为了凑一个好看公式，而是因为真空解是理解一切真实弹道的参照系：它是轨迹的"骨架"，空气只负责在上面做减法。
先把问题简化到真空
真空里弹丸只受重力，两个方向解耦：水平匀速，竖直匀加速（9.8 m/s²）。落点方程一行搞定：
真空抛物线三件套（初速 v0，射角 θ）： 飞行时间 T = 2·v0·sinθ / g 最大射程 R = v0²·sin(2θ) / g 弹道顶点 H = v0²·sin²θ / (2g) 例：v0 = 850 m/s，θ = 45° R = 850² × 1 / 9.8 ≈ 73 700 m ≈ 73.7 km T = 2 × 850 × 0.707 / 9.8 ≈ 122.7 s
73.7 公里、两分多钟的滞空——如果这成立，炮兵就不需要任何射表修正了。现实中 155 mm 榴弹在同样初速下的最大射程约 20-30 公里，滞空不过几十秒。真空公式高估射程约三倍，这就是空气的"税"。顺便注意公式里射角以 2θ 出现：θ=45° 最远；θ=30° 与 θ=60° 射程相同（sin60°=sin120°）——同一门炮的"高射界"与"低射界"弹道在真空里等价，实战中因空气分布不同而分道扬镳。
图：真空抛物线与真实弹道的对照
真实轨迹的三段结构
加回空气后，抛物线变成一条左右不对称的曲线，三段各有名堂：
上升段。弹丸刚出膛速度最高，阻力最大——注意阻力按速度平方走，速度掉一半，阻力掉四分之三。上升段是动能流失最快的阶段（步枪弹头几百米内就丢掉三成速度），但随着弹丸爬升、空气变稀、速度变慢，阻力迅速收敛。高射界射击（大射角）时弹丸会进入 10 公里以上的稀薄大气，阻力税明显减轻——这就是远程火炮偏爱高射界弹道的原因之一。
顶点。弹道倾角过零，速度到达全程最低。对旋转稳定弹丸这是最险的一段：速度低意味着陀螺"资本"相对充裕，但顶点附近的姿态扰动恢复得慢，散布圆往往在顶点附近被撕大。
下降段。重力重新掌舵，弹丸越落越陡，落点倾角比射角更陡（真空里两者应相等）。存速远低于初速：步枪弹在 800 米外存速常只剩初速的六成，穿甲能力随之下滑——外弹道算完的存速，正是终点弹道的开场数字。
动手算：把阻力加回来
点质量弹道的全部方程只有三行：两个方向的速度变化率加位置积分。用一个可运行的数值积分把真实弹道算出来：
import math def point_mass_trajectory(v0=850.0, theta_deg=45.0, dt=0.001): """点质量弹道积分（含平方阻力），返回射程与飞行时间。 Cd 取 0.25 的保守常数（真实弹丸随马赫数变化，见 4.2 节）""" g, rho = 9.8, 1.225 # 重力加速度与海平面空气密度 d, m = 0.155, 45.0 # 155 mm 弹径与弹重 A = math.pi * (d/2) ** 2 th = math.radians(theta_deg) x, y = 0.0, 0.0 vx, vy = v0 * math.cos(th), v0 * math.sin(th) t = 0.0 while y >= 0: v = math.hypot(vx, vy) Fd = 0.5 * rho * v * v * 0.25 * A # F = ½ρv²·Cd·A ax = -Fd / m * vx / v # 阻力方向与速度反向 ay = -Fd / m * vy / v - g vx += ax * dt; vy += ay * dt x += vx * dt; y += vy * dt t += dt return x, t R, T = point_mass_trajectory() print(f"带阻力的射程约 {R/1000:.1f} km，飞行时间约 {T:.0f} s") # 输出量级：射程 20 km 上下、时间数十秒——与真实 155 mm 榴弹同量级
十几行代码，把真空与现实的差距从 50 公里压缩到公里级。剩余误差来自三处：Cd 实际随马赫数变化（4.2 节）、空气密度随高度变化（第 7 章）、以及我们完全忽略的姿态效应（4.3 节）。工程建模的哲学在此显形：先抓主导项，再按需补细节。
常见坑
- 用真空公式估枪械弹道。手枪弹 350 m/s 的"真空射程"超过 12 公里，实际不到 2 公里——低速弹的阻力税更高（相对税基），真空公式失真更严重。
- 忽略射角与存速的关系。同一目标，用高射界打弹丸滞空长、受风时间长；用低射界打存速高、终点效应好。射角选择是弹道学与战术的双重决策。
低伸弹道的一次速算
枪械射击多为低伸弹道（射角几度以内），此时真空公式可以直接用来做"弹道高"速算——瞄具为什么要抬高、归零距离怎么选，全在这笔账里：
低伸近似：弹道下坠量 H ≈ ½ g t²，t = 距离 ÷ 平均速度 例：初速 940 m/s，飞到 300 m 用时约 0.36 s（带阻力） 下坠 ≈ 0.5 × 9.8 × 0.36² ≈ 0.63 m 含义：枪管必须上抬约 0.63 m/300 m ≈ 2.1 mrad 才能命中 ——这就是"归零"在做的几何补偿
再走一步：若归零在 300 m，目标在 150 m 处呢？飞行时间约 0.17 s，下坠仅 0.14 m，而瞄具仍按 0.63 m 抬着——弹着点会高出瞄准点约半米。这就是"近距反而打高"的经典现象，也是每一个归零距离只对应一条弹道高曲线的几何原因。会了这一步速算，看任何弹道高表都能心算验证。
本节要点回顾
- 真空三件套：R = v²sin2θ/g、T = 2v·sinθ/g、H = v²sin²θ/2g，45° 最远；
- 空气的税：850 m/s 的真空 73.7 km 到实际 20-30 km，差额按速度平方累积；
- 真实轨迹不对称：上升段丢速度最快，顶点姿态最险，下降段越落越陡、存速越低；
- 点质量可算：三行微分方程加数值积分，一个下午就能复现真实射程量级；
- 射角双解：30° 与 60° 真空等距，实战中因空气分层而分家——高射界走得远，低射界落得狠。
骨架有了，下一节把"空气税"的税率表摊开：阻力系数怎么随速度变、弹道系数怎么比、升力又在悄悄干什么。