2. 空气动力学基础

在固定翼气动系数辨识工程中,我们最终要估计的对象并不是“升力多少牛”“阻力多少牛”,而是无量纲气动系数。原因很直接:同一架飞机在不同高度、不同速度、不同载荷下,气动力绝对值差异巨大,但系数在相似条件下具有可比性和可建模性。因此,理解六分量气动力与力矩的物理来源,以及它们如何被无量纲化,是后续建模、试验设计和参数辨识的基础。

一、六分量:三个力与三个力矩

固定翼飞机在气流中受到的合力和合力矩,通常分解为体轴系下的六个分量:

  • 升力 L:垂直于来流方向、位于飞机对称面内的气动力分量,主要来自机翼上下表面压力差。
  • 阻力 D:平行于来流方向、阻碍飞机前进的气动力分量,包括摩擦阻力、压差阻力、诱导阻力和干扰阻力。
  • 侧力 Y:垂直于对称面、沿机体横轴方向的气动力分量,主要来自侧滑角、方向舵偏转和不对称气流。
  • 俯仰力矩 M:绕机体横轴、使机头抬头或低头的力矩,主要来自升力作用点与重心之间的力臂,以及平尾配平作用。
  • 滚转力矩 Lr:绕机体纵轴、使飞机左右倾斜的力矩,主要来自副翼偏转、侧滑和机翼不对称升力分布。
  • 偏航力矩 N:绕机体立轴、使机头左右偏转的力矩,主要来自方向舵、侧滑和垂尾侧力。

这六个分量共同决定了飞机的运动状态。辨识工作的核心,就是建立它们与飞行状态量(迎角、侧滑角、角速率、舵面偏度、马赫数等)之间的数学关系。

二、无量纲系数定义

为了消除尺度、速度和高度的影响,工程上采用动压 q̄ = ½ρV² 和参考面积 S、参考弦长 c̄、参考展长 b 进行无量纲化。常用定义如下:

分量 物理量 无量纲系数 典型表达式
升力 L CL L / (q̄S)
阻力 D CD D / (q̄S)
侧力 Y CY Y / (q̄S)
俯仰力矩 M Cm M / (q̄S c̄)
滚转力矩 Lr Cl Lr / (q̄S b)
偏航力矩 N Cn N / (q̄S b)

其中,力矩系数使用参考长度进行二次无量纲化:俯仰力矩用平均气动弦长 c̄,滚转和偏航力矩用翼展 b。这样得到的系数在几何相似的飞机之间具有可比性,也便于在辨识模型中统一量纲。

三、物理来源与状态依赖关系

每个系数都不是常数,而是多个状态量的函数。以纵向为例,升力系数和俯仰力矩系数通常写为:

C_L = C_L0 + C_Lα·α + C_Lq·(q c̄ / 2V) + C_Lδe·δe
C_m = C_m0 + C_mα·α + C_mq·(q c̄ / 2V) + C_mδe·δe

这里 α 是迎角,q 是俯仰角速率,δe 是升降舵偏度。CLα 反映迎角对升力的贡献,Cmq 反映俯仰阻尼,Cmδe 反映升降舵操纵效率。横航向类似,侧力、滚转力矩和偏航力矩主要依赖侧滑角 β、滚转角速率 p、偏航角速率 r 以及副翼和方向舵偏度。

关键理解:气动系数辨识不是“测力”,而是“从飞行数据中反推这些导数项”。每一个导数项都对应一个明确的物理效应,例如 CLα 对应机翼迎角升力线斜率,Cnβ 对应航向静稳定性,Clδa 对应副翼滚转操纵能力。

四、典型场景:从试飞数据到系数模型

假设一次纵向试飞中,飞机以固定油门做小振幅俯仰激励。飞控记录下迎角、俯仰角速率、升降舵偏度和法向过载。通过气动力模型,我们可以把法向过载换算为升力系数,把俯仰角加速度换算为俯仰力矩系数。然后利用最小二乘或极大似然估计,拟合出 CLα、CLq、Cmα、Cmq、Cmδe 等参数。

这个过程的关键在于:模型结构必须与物理来源一致。如果遗漏了俯仰阻尼项 Cmq,辨识结果就会把阻尼效应错误地归入 Cmα,导致短周期特性预测失真。同样,如果忽略侧滑对滚转力矩的贡献 Clβ,横航向耦合模态就会明显偏离真实飞机。

五、常见误区与复盘

误区一:把系数当成常数。气动系数随迎角、马赫数、雷诺数和舵面偏度变化,尤其在大迎角区会出现非线性甚至分离。辨识时必须限定飞行包线,或引入非线性项。

误区二:参考量不统一。不同文献可能用不同参考面积或参考弦长,导致系数数值不可比。工程中必须明确 S、c̄、b 的定义,并与风洞或 CFD 数据对齐。

误区三:忽略力矩参考点。俯仰力矩系数对参考点位置敏感。重心位置变化会直接改变 Cm0 和 Cmα,辨识前必须确认力矩参考点。

复盘时,建议把每个系数与物理效应一一对应,形成“系数—导数—物理意义—典型量级”的清单。这样在辨识结果异常时,可以快速判断是数据质量问题、模型结构问题,还是参数不可辨识问题。

固定翼六分量气动力/力矩与无量纲系数关系 机体轴系 升力 L → C_L 阻力 D → C_D 侧力 Y → C_Y 俯仰 M → C_m 滚转 L_r → C_l 偏航 N → C_n 无量纲化公式与状态依赖 C_L = L / (q̄S) = C_L0 + C_Lα·α + C_Lq·(q c̄ / 2V) + C_Lδe·δe C_D = D / (q̄S) = C_D0 + K·C_L² + C_Dδe·δe C_Y = Y / (q̄S) = C_Yβ·β + C_Yδr·δr C_m = M / (q̄S c̄) = C_m0 + C_mα·α + C_mq·(q c̄ / 2V) + C_mδe·δe C_l = L_r / (q̄S b) = C_lβ·β + C_lp·(p b / 2V) + C_lr·(r b / 2V) + C_lδa·δa + C_lδr·δr C_n = N / (q̄S b) = C_nβ·β + C_np·(p b / 2V) + C_nr·(r b / 2V) + C_nδa·δa + C_nδr·δr

掌握六分量的物理来源和无量纲定义,是后续建立气动模型、设计激励输入、选择辨识算法和评估参数可信度的前提。只有把“力从哪来、力矩绕哪转、系数怎么归一化”这三件事弄清楚,辨识结果才具有工程意义。