3. 空气动力学基础
在无人机动力学仿真中，气动力与力矩是连接“飞行状态”与“运动响应”的核心桥梁。若气动模型失真，再精确的刚体六自由度方程也无法复现实飞表现。本节从工程仿真视角梳理升力、阻力、侧力与三轴力矩的产生机理，并给出螺旋桨与机翼的常用气动模型及简化假设，为后续配平、线性化与实飞对标奠定基础。
一、核心原理：气动力与力矩从何而来
气动力本质上来自空气与机体表面相对运动产生的压力分布与剪切应力分布。对无人机而言，可将其分解为沿气流坐标系的三个分量：升力 L（垂直于来流、位于机体对称面内）、阻力 D（平行于来流、反向）、侧力 Y（垂直于对称面）。三者通常写成无量纲系数形式：
L = 0.5 * rho * V^2 * S * C_L
D = 0.5 * rho * V^2 * S * C_D
Y = 0.5 * rho * V^2 * S * C_Y
M = 0.5 * rho * V^2 * S * c * C_M其中 rho 为空气密度，V 为空速，S 为参考面积，c 为参考弦长。力矩则来自气动力作用点与重心不重合：俯仰力矩 M 主要由升力作用点（气动中心）与重心的纵向偏移产生，滚转与偏航力矩由左右机翼/螺旋桨不对称力及侧力产生。理解“力作用于何处”比记住公式更重要——这正是仿真中必须区分气动中心与重心的原因。
二、螺旋桨气动模型：从动量理论到叶素简化
螺旋桨是多数多旋翼与固定翼无人机的核心推力源。仿真中常用两类模型：
- 动量理论模型：将桨盘视为一维流管，推力 T 与功率 P 由诱导速度决定，适合悬停与低速前飞的整体性能估算。
- 叶素-动量组合模型（BEMT）：沿桨叶半径积分升力与阻力，能反映桨距、转速与来流角的影响，精度高但计算量大。
工程仿真中更常用简化参数模型：以转速 n 为输入，用推力系数 C_T 与功率系数 C_P 描述：
T = C_T * rho * n^2 * D^4
P = C_P * rho * n^3 * D^5
Q = P / (2 * pi * n)其中 D 为桨径，Q 为反扭矩。该模型假设 C_T、C_P 在悬停附近近似常数，忽略桨叶挥舞、压缩性与桨间干扰。对多旋翼而言，反扭矩 Q 直接产生偏航力矩，是偏航通道建模的关键。
三、机翼气动模型与简化假设
固定翼无人机的升力与阻力通常用线性升力线模型描述。在小迎角范围内，升力系数随迎角线性增长：
C_L = C_L0 + C_L_alpha * alpha
C_D = C_D0 + k * C_L^2
C_M = C_M0 + C_M_alpha * alpha仿真中常见的简化假设包括：
- 准定常假设：忽略非定常分离与动态失速，气动力仅由当前状态决定。
- 小迎角线性假设：C_L 与 alpha 线性，失速后需切换非线性模型。
- 忽略压缩性：马赫数较低时密度与温度效应可忽略。
- 对称假设：左右机翼气动特性一致，侧力与滚转力矩由不对称舵面或侧滑产生。
- 舵面效率线性：升降舵、副翼、方向舵的操纵导数近似常数。
四、案例与场景：悬停与巡航的气动差异
以四旋翼为例，悬停时螺旋桨主要克服重力，来流速度低，C_T 近似常数；前飞时桨盘受到斜向来流，前进桨与后退桨的局部迎角不同，产生桨盘不对称升力，进而引发滚转与俯仰耦合。若仿真仍使用悬停 C_T，则前飞速度超过 8 m/s 后姿态误差会明显增大。
固定翼巡航场景中，升力需平衡重力，阻力由巡航速度决定。若忽略 C_D0 随雷诺数的变化，低速巡航阻力会被低估，导致仿真航程偏大。实飞对标时，应优先校准 C_L_alpha、C_D0 与 k 三个参数。
| 气动对象 | 主要输出 | 关键参数 | 常用简化 |
|---|---|---|---|
| 螺旋桨 | 推力 T、反扭矩 Q | C_T、C_P、D、n | 系数常数、忽略挥舞 |
| 机翼 | 升力 L、阻力 D、俯仰力矩 M | C_L0、C_L_alpha、C_D0、k | 线性升力、准定常 |
| 机身 | 阻力、侧力 | C_D_body、C_Y_beta | 忽略干扰 |
| 舵面 | 操纵力矩 | C_L_delta、C_M_delta | 效率线性 |
五、误区与复盘
复盘建议：每次实飞对标后，按“先配平、再线性化、后非线性”的顺序回查气动参数。若悬停推力偏差大，优先检查 C_T 与桨径；若巡航阻力偏差大，优先检查 C_D0 与雷诺数修正；若姿态耦合明显，检查桨盘不对称与侧滑侧力模型。只有让气动模型与飞行状态匹配，动力学仿真才能真正成为实飞的可信预演工具。