2. 飞行器动力学基础

飞行器在飞行过程中，质量并非恒定不变。无论是消耗燃料的化学火箭、投放载荷的运输机，还是执行空中加油的受油机，质量的变化都会直接改变飞行器的加速度响应、姿态动力学特性以及推力需求。理解刚体运动方程是分析这一切的基石，而将质量变化项显式纳入动力学模型，则是设计自适应控制系统、处理质量扰动问题的前提。本节将系统梳理刚体运动方程的建立过程，重点剖析质量变化对线动量与角动量的耦合影响，并给出推力与重力随质量变化的建模方法。

2.1 刚体运动方程：从常质量到变质量的扩展

2.1.1 核心原理：变质量系统的牛顿-欧拉方程

经典刚体动力学假设质量恒定，其线动量方程为 F = m·a，角动量方程为 M = I·α + ω×(I·ω)。然而，飞行器在飞行中持续消耗燃料或投放载荷，质量 m(t) 和惯性张量 I(t) 均为时间的函数。此时必须回到牛顿第二定律的原始形式：F = d(mV)/dt，展开后得到：

F = m(t)·dV/dt + ṁ(t)·V

其中 ṁ(t) 为质量变化率（负值表示质量减少）。右侧第二项 ṁ·V 是质量变化引起的附加力，其物理本质是：被排出或吸入的质量携带动量，对飞行器主体产生反作用力。对于火箭发动机，这一项正是推力的来源之一；对于投放载荷的飞机，这一项表现为突然的动量损失。

角动量方程同样需要扩展。设飞行器本体坐标系原点取在瞬时质心，惯性张量 I(t) 时变，角动量 H = I(t)·ω。对时间求导：

M = I(t)·dω/dt + dI(t)/dt·ω + ω×(I(t)·ω)

其中 dI/dt·ω 项反映了质量分布变化（如燃料消耗导致惯性矩减小）对角加速度的影响。若质心位置本身也在漂移（例如液体燃料消耗导致质心移动），还需引入质心偏移引起的附加力矩项。

2.1.2 判断方法：何时必须考虑质量变化效应

工程实践中并非所有飞行器都需要完整变质量模型。判断依据如下：

判断准则：若 ṁ·V 项与 m·a 项量级比超过 5%，或 dI/dt·ω 项与 I·dω/dt 项量级比超过 5%，则必须采用变质量模型。对于高精度姿态控制（如卫星入轨、导弹末制导），即使质量变化率较小，也需在模型中显式加入。

2.1.3 案例拆解：液体火箭助推段的动力学建模

以某型两级液体运载火箭为例，起飞质量 300 吨，燃料质量 240 吨，发动机总推力 4000 kN，比冲 300 s。助推段持续 150 秒，燃料消耗率约为 1.6 吨/秒。

初始现象：在飞行至 60 秒时，火箭轴向加速度从 1.2g 上升至 2.8g，姿态控制系统出现明显的俯仰通道振荡，增益裕度下降。

分析过程：计算质量变化率 ṁ = -1600 kg/s，飞行速度 V = 900 m/s，则 ṁ·V = -1.44×10⁶ N，相当于推力的 36%。这一附加力项在固定增益控制器设计时被忽略，导致实际被控对象增益比模型高 36%，相位裕度下降约 12°。同时，燃料消耗使惯性张量 I_yy 从 2.1×10⁷ kg·m² 降至 1.2×10⁷ kg·m²，dI/dt·ω 项在俯仰角速度 5°/s 时产生约 8.7×10⁵ N·m 的附加力矩，相当于喷管摆动力矩的 15%。

解决方案：在控制律中引入质量变化补偿项：

δ_cmd = δ_base + K_m·(ṁ/m)·V + K_I·(dI/dt)·ω

其中 K_m 和 K_I 为自适应增益，通过在线质量估算器（基于燃料流量计和加速度计融合）实时更新。修改后，俯仰通道振荡幅值降低 70%，增益裕度恢复至 8 dB 以上。

复盘经验：变质量效应在助推段后期尤为显著，控制器设计必须覆盖全飞行包线内的质量范围，而非仅设计点。

2.2 质量变化对动力学的影响：质心漂移与惯性张量时变

2.2.1 核心原理：质量分布变化的力学效应

质量变化不仅改变总质量，还改变质量分布。对于液体燃料飞行器，燃料消耗顺序（先消耗哪个贮箱）直接影响质心位置。质心偏移 Δr_cg 会产生附加力矩：

ΔM = Δr_cg × F_total

其中 F_total 为作用在飞行器上的合力（推力+气动力+重力）。若质心沿机体轴向后移动，则俯仰静稳定性增加（静稳定裕度增大），但操纵效率下降；若质心前移，则稳定性降低，可能进入静不稳定状态。

惯性张量的时变同样关键。I(t) 的变化率 dI/dt 与燃料消耗速率、燃料在贮箱内的晃动有关。对于细长体飞行器，I_yy 和 I_zz 的变化最为显著，而 I_xx（滚转轴）变化相对较小。惯性张量变化直接改变姿态动力学的自然频率：

ω_n = sqrt(M_α / I_yy)

当 I_yy 减小时，自然频率升高，若控制器带宽不变，则相位滞后增大，稳定性裕度下降。

2.2.2 操作流程：质量特性在线估算方法

1. 燃料消耗积分法：通过流量计积分获得已消耗燃料质量，结合初始质量计算当前质量。精度取决于流量计误差（通常 1~3%）。
2. 加速度计-推力模型法：利用已知推力模型和实测加速度，反推质量 m = (T - D) / a_x。该方法对气动阻力模型敏感，适合高空低动压段。
3. 扩展卡尔曼滤波（EKF）融合：将质量、质心位置、惯性张量作为状态变量，融合加速度计、陀螺仪、燃油流量计、作动器位置等多源信息进行在线估计。典型状态向量为 [m, r_cg_x, r_cg_y, r_cg_z, I_xx, I_yy, I_zz]。
4. 查表插值法：预先通过地面试验或 CFD 计算不同燃料余量下的质量特性表，飞行中根据燃料余量查表。简单可靠，但无法覆盖非标称工况。

2.2.3 案例拆解：空中加油机的质量特性变化与配平补偿

某型加油机在执行任务时，起飞质量 180 吨，其中燃油 90 吨。在 30 分钟内输出 40 吨燃油给受油机，同时自身消耗 5 吨燃油。加油过程中，质心位置从 32% MAC（平均气动弦）后移至 38% MAC，I_yy 从 4.2×10⁶ kg·m² 降至 3.1×10⁶ kg·m²。

初始现象：加油开始后，飞行员报告飞机有持续低头力矩，自动驾驶仪配平舵面偏转从 -2° 增至 -8°，接近作动器行程极限。

分析判断：质心后移导致静稳定裕度从 8% MAC 增至 14% MAC，升力产生的低头力矩增大。同时，燃油从机翼油箱泵入机身油箱，改变了质量分布，I_yy 减小使俯仰响应加快，自动驾驶仪原增益偏小，响应滞后。

实际操作：在飞控系统中增加质量特性前馈补偿模块：

δ_trim = δ_trim_base + K_cg·(r_cg - r_cg_ref) + K_I·(I_yy - I_yy_ref)

同时将俯仰增益调度为 I_yy 的函数：K_q = K_q0·(I_yy0 / I_yy)。修改后，配平舵面偏转恢复至 -3°，自动驾驶仪在加油全过程中保持平稳。

结果对比：补偿前，加油过程中俯仰姿态偏差最大达 ±2.5°；补偿后降至 ±0.4°。作动器饱和时间从 12 秒降至 0 秒。

2.3 推力与重力变化建模

2.3.1 核心原理：推力模型与重力场模型

推力模型：发动机推力是油门、飞行高度、马赫数、大气温度的函数。对于涡扇发动机，可用简化模型：

T = T_sl·δ·f(Ma) · (油门系数)

其中 T_sl 为海平面静推力，δ 为压力比（P/P_0），f(Ma) 为马赫数影响函数。对于火箭发动机，推力近似为：

T = ṁ_fuel·Isp·g0 + (P_e - P_a)·A_e

其中第一项为动量推力，第二项为压力推力。当飞行高度增加、大气压 P_a 降低时，压力推力增大，总推力上升。

重力模型：重力加速度随高度变化：

g(h) = g0·(R_e / (R_e + h))²

其中 R_e 为地球半径（6378 km），g0 = 9.80665 m/s²。在近地空间（h < 100 km），g 变化小于 3%，但高超声速飞行器或亚轨道飞行器必须考虑。重力方向始终指向地心，在飞行器本体坐标系中需通过姿态角投影。

2.3.2 判断方法：推力与重力建模的精度需求

2.3.3 案例拆解：高超声速飞行器爬升段推力-重力耦合建模

某型高超声速试验飞行器采用超燃冲压发动机，飞行包线为 Ma 5~7，高度 25~40 km。在爬升段，发动机推力随动压变化显著，而重力加速度随高度增加而减小，两者耦合导致纵向加速度出现非线性。

初始现象：在 30 km 高度、Ma 6 条件下，飞行器纵向加速度出现 0.15g 的波动，与地面仿真结果偏差 20%。

分析过程：检查发现，地面仿真中推力模型采用了固定高度下的推力曲线，未考虑高度变化引起的压力推力增量。实际飞行中，高度从 28 km 升至 32 km，大气压从 1600 Pa 降至 870 Pa，压力推力项 (P_e - P_a)·A_e 增大了约 18%。同时，重力加速度从 9.71 m/s² 降至 9.68 m/s²，虽然变化仅 0.3%，但叠加推力增加后，净加速度变化被放大。

解决方案：建立推力-高度-马赫数三维插值表，并在飞控模型中引入重力梯度补偿项：

a_z_cmd = (T·sinα + L - m·g(h)) / m + Δg(h)

其中 Δg(h) = g(h) - g0 作为前馈补偿。修改后，纵向加速度波动降至 0.03g 以内，仿真与飞行数据吻合度提升至 95% 以上。

2.3.4 常见误区与复盘

误区一：将质量变化率 ṁ 视为常数。实际上，燃料消耗率随油门位置和飞行状态变化，火箭的 ṁ 在级间分离时还会发生阶跃跳变。

误区二：忽略质心偏移对力矩的影响。即使总质量变化不大，质心移动 1% MAC 也可能产生相当于 10% 舵面效率的力矩变化。

误区三：推力模型中使用海平面推力而不修正高度。在 15 km 高度，涡扇发动机推力可能降至海平面的 40%，若不修正，控制器将严重高估可用推力。

误区四：重力模型使用常数 g0。对于高空长航时飞行器（HALE），高度 20 km 时 g 已减小 0.6%，累积速度误差在 10 分钟飞行中可达 3.5 m/s，对航迹精度影响显著。

复盘建议：在飞行器设计初期就建立完整的变质量六自由度仿真模型，将质量变化、质心漂移、推力-重力耦合纳入控制器设计基线。飞行试验前，通过蒙特卡洛仿真覆盖质量变化率的不确定性范围，确保控制器在全包线内具有足够的鲁棒性。