9.1 解析模型

9.1 解析模型

本节摘要：解析模型指能用闭式公式直接算出结果的模型——真空抛物线、平射近似、内弹道解析解是弹道学三大经典。它们精度有限，却承担着数值时代不可替代的三种角色：秒级的量级估算、数值结果的校验标尺、参数敏感性的直观来源。本节梳理解析模型的谱系与适用边界，并演示"闭式解当标尺"的正确用法。本章从计算链的最便宜一档开始。

位置感：第 4 章已经用过真空抛物线（解析模型的最著名代表），本节把它放进完整谱系，并给出"什么时候用公式、什么时候上积分器"的工程判断。

一页纸上的弹道

解析模型按复杂度排成谱系，每一档用一行代码的钱回答一类问题：

真空闭式解。4.1 节的三件套（R、T、H）是零成本基线：任何数值弹道计算的第一步都是拿真空解对表——积分器在 Cd=0 时必须复现闭式解到小数位，否则后面的一切都是空中楼阁。这是"解析模型当标尺"的用法：精度不够的模型，恰好适合检验更复杂模型的正确性。

平射近似。西亚切在 19 世纪给出的聪明简化：对小射角、弹道低伸的情形（枪械与直射火炮），把空气密度当常数、把射角的高次项丢掉，弹道方程能积出闭式或查表形式。它至今活在两处：轻武器弹道学的教学体系、以及"低伸弹道"直觉的来源。适用边界同样清楚：射角超过 15° 后误差迅速失控——高射界的弹必须回到数值积分。

内弹道解析解。第 2 章的膛压-时间曲线，在简化假设（等燃面、绝热、集中参数）下可以积分出峰值膛压与初速的闭式估计。工程上它用于方案筛选：十几种装药方案先用解析式扫一遍量级，留下两三种上数值程序细化。粗、快、方向正确，是解析模型在工程链条里的标准站位。

点质量闭式修正。介于纯解析与全数值之间的小聪明：在真空解上叠加一阶阻力修正项，得到"带阻力的近似射程"。4.1 节的数值实验显示真空解高估三倍，而一阶修正能把误差压到两成上下——对脑内估算，两成已经奢侈地够用。

动手：把闭式解当标尺用

解析模型最有教学价值的用法是校验数值积分器。下面这个真空对照测试是所有弹道代码的第一道单元测试：

import math

def vacuum_closed_form(v0, theta_deg, g=9.8):
    """真空闭式解：射程、飞行时间、顶点高"""
    th = math.radians(theta_deg)
    R = v0**2 * math.sin(2*th) / g
    T = 2 * v0 * math.sin(th) / g
    H = (v0*math.sin(th))**2 / (2*g)
    return R, T, H

def numeric_vacuum(v0, theta_deg, g=9.8, dt=0.001):
    """数值积分（无阻力）：应复现闭式解"""
    th = math.radians(theta_deg)
    x = y = 0.0
    vx, vy = v0*math.cos(th), v0*math.sin(th)
    while y >= 0:
        vx, vy = vx, vy - g*dt      # 唯一的力：重力
        x, y = x + vx*dt, y + vy*dt
    return x

R, T, H = vacuum_closed_form(850, 45)
print(f"闭式解：射程 {R:.1f} m，时间 {T:.2f} s，顶点 {H:.1f} m")
print(f"数值解：射程 {numeric_vacuum(850, 45):.1f} m")
# 两行输出应在 0.1% 内吻合——吻合才允许加阻力项继续开发
​

这个动作看似平凡，却体现了计算工程的纪律：每一层复杂性都必须在上一层验证通过后叠加。积分器加阻力前先对真空表，加密度分层前先对均匀大气表，加六自由度前先对点质量——每一级都有上一级当免费标尺。

解析模型的思想红利

除了当标尺，闭式解还赠送三种直觉，数值解给不了：

平方律直觉。R = v²sin2θ/g 一眼看出射程对初速平方敏感——初速提升 10%，真空射程涨 21%。带阻力后敏感度降到约 1.5-2 倍（4.2 节的放大器反向），但"速度是弹道学第一货币"的直觉由此扎根。

对称性直觉。sin2θ 的对称性告诉你 30° 与 60° 真空等距；真实弹道破坏对称后，"哪个射角更优"就变成空气分层结构的问题——解析解给出基准线，偏差本身携带信息。

极限直觉。让 v 趋于无穷，真空射程趋于无穷而真实射程趋于某个极限——空气税的"封顶效应"直接写在解析式的对比里。远射程火炮与其拼命提初速，不如优化弹道系数与弹道（爬升到稀薄大气）——这个战略判断用两行公式就能论证。

常见坑

• 拿解析模型做最终设计依据。闭式解的每一项简化（均匀大气、常数 Cd、忽略姿态）在正式设计里都要偿还——它的岗位是标尺与筛选，不是交付；
• 在解析解适用域外硬用。平射近似套在大射角榴弹上，误差能到五成——用之前先确认自己的弹道在不在"低伸"这条边界内。

一道"先解析后数值"的完整示范

把本节的方法论串成一个真实感的工程小案例：为 800 m 训练靶道预估一块挡弹土堤的最小厚度。

第一步 解析粗估（选弹道顶点）：
  低伸弹道在 800 m 的下坠 H ≈ ½g·t²，
  平均速度约 750 m/s → t ≈ 1.07 s → 下坠约 5.6 m
  弹道在 800 m 处低于瞄准线约 5-6 m → 土堤核心高度有数了

第二步 解析定风险角：
  落点弹道倾角 θ ≈ arctan(g·t / v) ≈ arctan(9.8×1.07/750) ≈ 0.8°
  ——入射极平，按 5.3 节的临界角规律，土堤坡面必须
    设计成让任何命中处的入射角都远在跳弹临界角以内

第三步 数值复核：
  用 9.2 节的点质量积分器复算 800 m 弹道高与倾角，
  与解析估值差 5% 以内 → 解析答案可用，土堤按此加余量
​

注意工作次序：解析解先给量级与设计约束（高度、角度），数值解只做最后复核。如果反过来先跑数值再想约束，省下的那十分钟往往要在设计返工时加倍还回去。

本节要点回顾

• 谱系五档：真空解、平射近似、内弹道解析、点质量数值、六自由度/有限元，成本与精度逐级上升；
• 标尺用法：解析解是数值代码的单元测试基准，每层复杂性在上层验证后叠加；
• 思想红利：平方律、对称性、极限封顶——闭式解免费赠送的工程直觉；
• 站位纪律：解析模型管筛选与量级，数值模型管设计与交付，岗位不可互换；
• 边界意识：平射近似限低伸弹道，真空解限粗估——用前先核对自己在哪条边界内。

标尺立好了，下一节把主力工具搬上来：从 RK4 积分器写到有限元，弹道计算的完整兵器谱。