2. 风场与大气扰动基础

无人飞行器在空中受到的“风”并不是一个单一量。对飞控与抗扰算法而言,真正需要区分的是不同时间尺度、空间尺度和物理成因的大气运动:它们作用在机体上的方式不同,建模时所需的自由度、频带和随机性描述也不同。若把风统一当成一个常值扰动,控制器在遇到湍流或下击暴流时就会严重失配。本节先建立“风场分层”的概念,再逐类说明作用机理与建模需求。

一、按时间与空间尺度给风场分层

工程上常把大气扰动按尺度分为三类。第一类是平均风,指在数十秒到数分钟内相对稳定的风速风向,可视为准常值,主要产生配平偏移与航迹漂移。第二类是阵风与湍流,时间尺度从零点几秒到十几秒,属于随机过程,主要激发机体的短周期响应与结构振动。第三类是风切变、热气流与下击暴流,属于强空间梯度或强垂直运动事件,可能在数秒内造成显著的高度与姿态突变,是安全边界问题而非单纯精度问题。

建模时首先要明确参照系约定:本文统一采用机体坐标系描述气动力,采用地坐标系(NED,北-东-地)描述风场与航迹。风速矢量在地坐标系下表示为 V_w = [u_w, v_w, w_w],其中 w_w 向下为正。相对气流速度为地速与风速之差,即 V_a = V_g - V_w,空速、迎角与侧滑角均由 V_a 计算。这一约定在后续所有图示与公式中保持一致。

二、平均风:配平偏移与航迹漂移

平均风对无人飞行器的作用可分解为两部分。其一是气动配平变化:水平风改变相对气流方向,使迎角与侧滑角偏离无风配平值,需要升降舵、副翼或方向舵提供附加偏转。其二是航迹漂移:若飞控只稳定姿态而不补偿风,机体将保持姿态但沿风矢量方向漂移,形成位置误差。以某小型固定翼为例,巡航空速 18 m/s,遭遇 6 m/s 正侧风,若仅做姿态保持,1 分钟后横向漂移约 360 m,远超常规航点容差。

平均风的建模需求相对简单:给定风速、风向与高度剖面即可。常用对数律或幂律描述近地风剖面,例如幂律 V(h) = V_ref · (h/h_ref)^α,其中 α 在开阔地取 0.14 左右,在城市或林地可达 0.25~0.4。建模时需注意:平均风不是“无扰动”,它决定了抗扰控制器的基准工作点。

三、阵风与湍流:随机过程与频带匹配

阵风是离散的、可辨识的风速突增,常用“1-cos”模型描述:u_g(t) = (U_g/2)·[1 - cos(πt/T_g)],其中 U_g 为阵风幅值,T_g 为持续时间,典型值 1~5 s。湍流则是连续随机场,工程上多用 Dryden 或 von Kármán 功率谱生成。Dryden 纵向谱为 Φ_u(ω) = σ_u²·(2L_u/πV)·1/[1+(L_uω/V)²],其中 σ_u 为湍流强度,L_u 为湍流尺度。

湍流建模的关键不是“加噪声”,而是频带匹配:湍流能量集中在 0.1~2 Hz,恰好覆盖多数小型无人机的短周期模态频率。若控制器带宽低于湍流主频带,扰动将被“穿透”到姿态与高度通道;若带宽过高,又会放大传感器噪声。因此湍流模型必须与飞控带宽、传感器采样率一起评估。下表给出典型参数示例(属教学假设值,非实测)。

扰动类型典型幅值时间/空间尺度主要影响通道
平均风3~10 m/s分钟级航迹、配平
阵风2~8 m/s1~5 s姿态、空速
湍流σ 0.5~3 m/s0.1~2 Hz姿态、角速率
风切变0.05~0.2 s⁻¹数十米空速、升力
热气流1~5 m/s 上升百米级高度、能量
下击暴流10~30 m/s 下沉数百米高度、姿态、空速

四、风切变、热气流与下击暴流:强梯度事件

风切变指风速或风向在空间上的急剧变化,典型形式是水平风切变与垂直风切变。垂直风切变直接改变相对气流:当无人机从顺风切变进入逆风切变时,空速骤增,升力上升,机体被“顶起”;反之空速骤降,可能触发失速告警。建模时通常用切变率 ∂V/∂h 或 ∂V/∂x 描述,微下击暴流的切变率可达 0.1 s⁻¹ 量级。

热气流是垂直上升气流,对滑翔类无人机是能量来源,对多旋翼则是高度扰动。其空间尺度通常为几十到几百米,中心上升速度 1~5 m/s。建模可用高斯型垂直速度场 w(x,y) = w_0·exp(-r²/2R²),其中 w_0 为中心强度,R 为半径。

下击暴流是最危险的低空风场事件:强下沉气流触地后向外扩散,形成“下沉—外流”结构。无人机穿越时先遭遇逆风增大(空速上升),随后进入下沉区(升力损失),再遭遇顺风切变(空速骤降),三者叠加可在数秒内造成掉高与姿态失控。建模需同时给出下沉速度剖面与外流水平风剖面,常用双涡环或工程经验剖面近似。

风场扰动分类与对无人机的作用通道 大气扰动源 平均风 准常值 · 分钟级 阵风 1-cos · 1~5 s 湍流 Dryden · 0.1~2 Hz 风切变 空间梯度 · 0.1 s⁻¹ 热气流 垂直上升 · 1~5 m/s 下击暴流 下沉+外流 · 10~30 m/s 作用机理 改变相对气流 V_a = V_g − V_w 迎角 α、侧滑角 β 偏移 气动力/力矩变化 升力、阻力、俯仰/滚转力矩 姿态与角速率响应 短周期模态被激发 空速与高度突变 升力损失、掉高、失速风险 航迹漂移与位置误差 航点容差、任务失败 建模需求 常值+剖面模型 幂律/对数律 离散阵风模型 1-cos 幅值/时长 随机湍流谱 Dryden/von Kármán 空间梯度场 切变率、涡环剖面 频带与带宽匹配 与飞控带宽联合评估 参照系约定:地坐标系 NED 描述风场,机体坐标系描述气动力;V_a = V_g − V_w

五、建模需求与工程取舍

把上述扰动落到仿真与控制器设计中,需要回答三个问题:用什么模型、放在哪个通道、以什么频带注入。平均风与风切变适合作为确定性输入,直接叠加到地速与空速计算中;阵风适合作为离散事件触发,用于测试瞬态响应;湍流适合作为连续随机过程,用于评估统计性能与鲁棒性;热气流与下击暴流则更适合作为场景化事件,用于安全边界验证。

一个常见的工程做法是分层建模:在配平层用平均风确定基准舵面偏度;在控制层用 Dryden 湍流谱生成 u_w, v_w, w_w 三通道随机风,采样率不低于 50 Hz;在安全层用下击暴流剖面做蒙特卡洛打靶,统计掉高与姿态超调。需要注意的是,湍流强度 σ 与尺度 L 必须与飞行高度、地形和气象条件匹配,不能直接套用固定值。

建模检查要点:① 风速矢量是否在地坐标系定义并与空速计算一致;② 湍流谱的频带是否覆盖机体短周期频率;③ 风切变与下击暴流是否作为独立安全场景而非普通扰动;④ 仿真采样率是否满足湍流最高频率的 5~10 倍。

理解这些扰动的物理成因与尺度差异,是后续设计抗扰控制器的前提:平均风决定工作点,阵风与湍流决定带宽与鲁棒性,风切变、热气流与下击暴流决定安全边界。只有把“风”拆成这几种可建模的对象,控制器才能有针对性地分配权限与增益。