4. 固定翼六自由度运动方程

在固定翼气动系数辨识中，六自由度运动方程是连接“可测飞行状态”与“待辨识气动参数”的桥梁。没有这组方程，试飞数据只是时间序列；有了它，才能把角速度、加速度、舵面偏度与气流角组织成可回归的模型结构。本节从牛顿-欧拉方程出发，说明力与力矩平衡关系、刚体假设的边界，以及工程上常用的简化形式。

一、核心原理：从牛顿第二定律到六自由度

固定翼飞机在空间中的运动可分解为质心平动的三个自由度和绕质心转动的三个自由度。设机体坐标系为 \(O_b x_b y_b z_b\)，速度向量为 \(V=[u,v,w]^T\)，角速度向量为 \(\omega=[p,q,r]^T\)，外力合力为 \(F\)，外力矩合力为 \(M\)，则牛顿-欧拉方程可写为：

平动：m (dV/dt + ω × V) = F
转动：I dω/dt + ω × (I ω) = M

其中 \(m\) 为质量，\(I\) 为惯性矩阵。第一式描述质心运动，第二式描述绕质心的转动。两式中的叉乘项分别代表科氏加速度与陀螺力矩，是机体旋转带来的耦合项。若忽略这些耦合，在大机动或快速滚转状态下会产生显著建模误差。

展开到分量形式，平动方程可写为：

m( u̇ + q w - r v ) = X
m( v̇ + r u - p w ) = Y
m( ẇ + p v - q u ) = Z

转动方程可写为：

Ixx ṗ - Ixz ṙ + q r (Izz - Iyy) - Ixz p q = L
Iyy q̇ + r p (Ixx - Izz) + Ixz (p² - r²) = M
Izz ṙ - Ixz ṗ + p q (Iyy - Ixx) + Ixz q r = N

这里 \(X,Y,Z\) 为气动力与推力在机体轴上的合力，\(L,M,N\) 为对应力矩。气动系数辨识的任务，就是把 \(X,Y,Z,L,M,N\) 表达为 \(V,\alpha,\beta,p,q,r,\delta_e,\delta_a,\delta_r\) 的函数，再通过飞行数据估计其中的导数项。

二、力与力矩平衡：方程右边的物理来源

方程左边是惯性项，右边是外力与外力矩。对固定翼而言，右侧通常包括气动力、发动机推力与重力。重力在机体轴上的投影与姿态角有关：

X_g = -m g sinθ
Y_g =  m g sinφ cosθ
Z_g =  m g cosφ cosθ

气动力通常写在气流坐标系中，再投影到机体轴。升力 \(L\)、阻力 \(D\)、侧力 \(Y\) 与迎角 \(\alpha\)、侧滑角 \(\beta\) 的关系决定了投影矩阵。力矩则主要来自气动焦点偏移、舵面偏转和角速度阻尼。辨识时，常把力矩系数写成：

Cl = Cl0 + Clβ β + Clp (p b / 2V) + Clr (r b / 2V) + Clδa δa + Clδr δr
Cm = Cm0 + Cmα α + Cmq (q c / 2V) + Cmδe δe
Cn = Cn0 + Cnβ β + Cnp (p b / 2V) + Cnr (r b / 2V) + Cnδa δa + Cnδr δr

这些系数导数就是气动系数辨识的核心输出。方程中的力与力矩平衡，实际上要求我们同时匹配平动与转动两组测量：加速度计提供比力，陀螺提供角速度，气压与惯导提供空速和姿态。若只匹配其中一组，辨识结果往往不满足物理一致性。

三、刚体假设与工程简化形式

刚体假设包含三点：飞机质量与惯性矩在辨识窗口内不变；机体为刚性，忽略结构弹性模态；质量分布对称，通常取 \(I_{xy}=I_{yz}=0\)，仅保留 \(I_{xz}\)。对大多数固定翼飞机，这一假设在低频段足够准确；但若机体存在明显燃油晃动或大展弦比机翼柔性，则需引入等效弹性或时变惯量修正。

工程简化常见有三类。第一，小扰动线性化：在配平点附近将方程线性化，得到纵向与横航向解耦的线性模型，适合低阶等效系统辨识。第二，准定常气动假设：忽略非定常迟滞，把气动力矩写成瞬时状态的代数函数，适合大多数常规试飞。第三，忽略 \(I_{xz}\) 或忽略某些角速度耦合项，用于初步参数估计，但最终模型应逐步放开。

简化形式	保留项	忽略项	适用场景
全六自由度	全部惯性耦合、\(I_{xz}\)	无	大机动、快速滚转、高精度辨识
小扰动线性化	一阶导数、线性气动	高阶项、大迎角非线性	配平点附近、低阶等效系统
纵向/横航向解耦	各自通道主要项	交叉耦合	短周期、荷兰滚初步分析
准定常气动	瞬时状态代数关系	非定常迟滞、动态失速	常规试飞、参数初步估计

四、案例场景：从试飞数据到方程残差

假设一次纵向试飞中，采样得到 \(V,\alpha,q,\theta,a_x,a_z,\delta_e\)。用简化纵向方程可写出：

m( u̇ + q w ) = X_a + X_g + T
m( ẇ - q u ) = Z_a + Z_g
Iyy q̇ = M_a + M_t

把气动系数模型代入，得到关于 \(C_{L0},C_{L\alpha},C_{Lq},C_{L\delta e},C_{m0},C_{m\alpha},C_{mq},C_{m\delta e}\) 的线性回归问题。若方程残差在机动段突然增大，通常说明简化假设失效，例如迎角进入非线性区、舵面速率限制或推力模型不准。此时应回到全六自由度方程，检查耦合项与 \(I_{xz}\) 的影响。

五、常见误区与复盘

• 误区一：把机体轴方程直接当气流轴方程。两者相差迎角与侧滑角投影，混用会导致升阻力系数符号和量级错误。
• 误区二：忽略重力投影。在姿态变化较大的机动中，重力项与气动力量级相当，漏掉会严重污染 \(C_L\) 与 \(C_m\) 的估计。
• 误区三：过早解耦。横航向耦合在滚转机动中不可忽略，解耦模型可能给出看似拟合好但物理错误的参数。
• 误区四：把刚体假设当绝对真理。柔性机体、燃油晃动和传感器安装偏差都会破坏刚体假设，需在残差分析中识别。

复盘时，建议按“方程结构—气动模型—测量误差—简化假设”四层排查。先确认六自由度方程左右量纲与坐标系一致，再检查气动系数模型是否覆盖主要非线性，然后评估传感器偏差与延迟，最后判断简化形式是否匹配当前机动。只有方程、模型与数据三者自洽，辨识出的气动系数才具备工程可信度。