4、经典控制理论回顾:传递函数、根轨迹、频域分析(Bode/Nyquist)在飞控中的应用

飞行控制律设计的本质,是在“被控对象(飞机本体)—控制器(飞控计算机)—传感器—执行机构”构成的闭环回路中,通过调节开环增益与动态补偿结构,使系统同时满足稳定性、机动性与鲁棒性要求。经典控制理论虽然诞生于上世纪中叶,但其物理概念清晰、工程验证手段成熟,至今仍是飞控律初步设计与故障排查的第一工具。本节将系统回顾传递函数建模、根轨迹整形、Bode图与Nyquist判据四大核心模块,并重点说明它们在飞控纵向/横航向通道中的具体落点。

经典控制理论在飞控律设计中的工程映射 被控对象 飞机本体(刚体+气动) 传递函数 G(s) 短周期/长周期/滚转/荷兰滚 控制器 K(s) PID/超前滞后/陷波器 根轨迹分析 增益变化→极点移动轨迹 频域分析 Bode图/Nyquist曲线 稳定性判据 幅值裕度/相位裕度/Nyquist 闭环反馈(传感器测量→飞控计算机→舵面) 飞控律设计中的典型应用场景 • 纵向短周期阻尼不足 → 俯仰角速率反馈(增益+超前) • 荷兰滚模态不满足规范 → 偏航阻尼器(洗出网络+侧向加速度反馈) • 结构模态耦合 → 陷波器(频域整形) • 舵机延迟 → 相位裕度补偿 • 传感器噪声 → 低通滤波(高频衰减) • 增益调度 → 不同飞行状态切换

上图展示了经典控制理论在飞控律设计中的整体映射关系。从飞机本体建模得到传递函数,通过根轨迹观察闭环极点随增益的变化趋势,再借助Bode图与Nyquist曲线评估系统的频域特性与稳定裕度,最终形成完整的控制律设计闭环。下面逐一展开各模块的核心原理与工程应用。

4.1 传递函数:飞控律设计的“物理语言”

核心原理

传递函数将线性时不变系统的微分方程转换为复域代数表达式,其标准形式为:

G(s) = (b_m·s^m + b_(m-1)·s^(m-1) + ... + b_0) / (a_n·s^n + a_(n-1)·s^(n-1) + ... + a_0)

在飞控领域,传递函数的核心价值在于:将气动导数、质量特性、惯性矩等物理参数映射为可计算的极点/零点位置。例如,纵向短周期传递函数可简化为:

Δθ(s)/Δδe(s) = (M_δe·(s + 1/T_θ)) / (s² + 2·ζ_sp·ω_sp·s + ω_sp²)

其中,ω_sp为短周期固有频率,ζ_sp为短周期阻尼比,M_δe为升降舵操纵导数。这些参数直接决定了飞机的俯仰响应品质。

飞控中的典型传递函数

通道典型传递函数关键参数物理意义
纵向短周期Δθ/Δδeω_sp, ζ_sp俯仰姿态对升降舵的响应
纵向长周期ΔV/Δδtω_ph, ζ_ph速度对油门/阻力的响应
横航向滚转Δp/Δδaτ_r(滚转时间常数)滚转角速率对副翼的响应
横航向荷兰滚Δβ/Δδrω_d, ζ_d侧滑角对方向舵的响应

工程应用要点

  • 模型降阶:飞控设计通常忽略非主导极点(如执行机构高频动态),保留2~4阶主导模态。
  • 参数不确定性:气动导数随马赫数、迎角变化,需在包线内选取多个设计点。
  • 零点影响:传递函数零点会改变根轨迹形状,需特别关注右半平面零点(非最小相位系统)。
经验提醒:在建立飞控传递函数时,务必先进行量纲一致性检查。例如,俯仰角速率反馈增益的单位应为“度/度/秒”或“rad/(rad/s)”,避免因单位混乱导致增益计算错误。

4.2 根轨迹法:增益调参的“导航图”

核心原理

根轨迹描述的是当开环增益K从0变化到∞时,闭环特征方程根的移动轨迹。其绘制规则基于幅值条件和相角条件:

幅值条件:|K·G(s)·H(s)| = 1
相角条件:∠G(s)·H(s) = (2k+1)·180°

在飞控中,根轨迹法主要用于确定反馈增益的合理范围,以及观察极点移动趋势。例如,当增加俯仰角速率反馈增益时,短周期阻尼比会增大,但过大的增益可能导致高频模态失稳。

飞控根轨迹分析流程

  1. 建立开环传递函数:包含飞机本体、舵机、传感器模型。
  2. 绘制根轨迹:使用MATLAB/Simulink或Python控制库。
  3. 标注设计约束:在s平面上画出阻尼比等值线(如ζ=0.35~1.3)和自然频率范围。
  4. 选择增益:使闭环极点落在满足飞行品质规范(如MIL-F-8785C)的区域。
  5. 验证鲁棒性:在气动参数±20%摄动下,检查极点是否仍处于稳定区域。

案例拆解:纵向短周期阻尼器设计

背景与约束:某轻型公务机在巡航状态(高度8000m,马赫0.6)下,短周期阻尼比仅为0.25,低于适航要求的最小值0.35。飞控计算机采样率为100Hz,舵机带宽为30rad/s。

初始现象:飞行员反馈俯仰振荡明显,操纵后飞机出现持续3~4个周期的摆动。

关键输入数据:短周期传递函数为 Δθ/Δδe = 8.5(s+1.2)/(s²+1.8s+9.0),即ω_sp=3.0 rad/s,ζ_sp=0.3。

分析过程:

  1. 绘制开环根轨迹,发现增加俯仰角速率反馈增益Kq时,短周期极点向左半平面移动,阻尼比增大。
  2. 当Kq=0.15时,闭环极点位于s=-2.1±j2.8,对应ζ=0.6,满足要求。
  3. 但继续增大Kq至0.3时,出现高频复极点(s=-12±j15),可能激发舵机模态。

实际操作:选择Kq=0.18,并加入一阶低通滤波器(截止频率25rad/s)抑制高频噪声。

结果对比:闭环阻尼比提升至0.58,超调量从35%降至12%,调节时间从6.2s缩短至2.8s。

风险与代价:低通滤波器引入约2.3°的相位滞后,需通过前置超前补偿修正。

复盘经验:根轨迹法能直观展示增益变化对模态的影响,但必须结合频域分析验证相位裕度,避免“根轨迹看似稳定,实际频域裕度不足”的情况。

风险提示:根轨迹法只适用于线性时不变系统。当飞机迎角超过失速迎角或进入深失速状态时,气动导数发生剧烈非线性变化,根轨迹分析结果将不再有效,必须切换至非线性仿真或增益调度策略。

4.3 频域分析:Bode图与Nyquist判据的工程应用

核心原理

频域分析通过系统对正弦输入的稳态响应来评估性能。Bode图包含幅频特性与相频特性两条曲线,Nyquist曲线则是在复平面上绘制开环频率响应的极坐标图。飞控设计中,频域分析的核心指标包括:

  • 幅值裕度(GM):系统在相位穿越频率处,幅值低于0dB的余量,一般要求≥6dB。
  • 相位裕度(PM):系统在增益穿越频率处,相位高于-180°的余量,一般要求≥45°。
  • 增益穿越频率(ω_c):决定系统响应速度,飞控中通常取2~10 rad/s。

Bode图在飞控中的典型应用

应用场景Bode图特征设计对策
传感器噪声抑制高频段幅值过高加入低通滤波器,降低高频增益
结构模态耦合在结构频率处出现谐振峰设计陷波器,在特定频率处衰减
舵机延迟补偿相位曲线在穿越频率处下降过快加入超前校正,提升相位裕度
阵风扰动抑制低频段增益不足增加积分环节或提高低频增益

Nyquist判据的飞控应用

Nyquist判据通过观察开环频率响应曲线围绕(-1, j0)点的次数来判断闭环稳定性。在飞控中,Nyquist图特别适用于:

  1. 多回路系统分析:当存在内回路(如阻尼回路)和外回路(如姿态保持回路)时,Nyquist图能清晰展示各回路的稳定裕度。
  2. 时滞系统稳定性:舵机延迟、传感器采样延迟等纯时滞环节在Nyquist图上表现为螺旋线,便于评估延迟对稳定性的影响。
  3. 鲁棒性验证:在气动参数摄动下,Nyquist曲线会形成“曲线带”,检查曲线带是否包围(-1, j0)点。

案例拆解:横航向荷兰滚抑制

背景与约束:某战斗机在跨声速飞行时,荷兰滚模态阻尼比仅为0.08,严重偏离规范要求(≥0.4)。飞控系统采用偏航阻尼器,通过方向舵反馈偏航角速率。

初始现象:飞行员报告飞机在侧风条件下出现明显的“蛇形”摆动,摆动周期约2.5s。

关键输入数据:荷兰滚模态传递函数为 Δβ/Δδr = 12.5(s+0.5)/(s²+0.8s+6.25),即ω_d=2.5 rad/s,ζ_d=0.16。

分析过程:

  1. 绘制Bode图发现,在ω=2.5 rad/s处,相位滞后约-160°,相位裕度仅12°。
  2. Nyquist曲线显示,开环频率响应在第三象限接近(-1, j0)点,距离仅0.15,表明稳定裕度严重不足。
  3. 设计偏航阻尼器增益K_r=0.8,并加入洗出网络(高通滤波器,截止频率0.5 rad/s)以消除稳态侧滑。

实际操作:在Simulink中搭建非线性模型,加入阻尼器后进行线性化,验证闭环极点位置。

结果对比:荷兰滚阻尼比提升至0.45,相位裕度增至52°,幅值裕度达到8.5dB。

风险与代价:洗出网络在低频段引入相位超前,可能影响协调转弯性能,需通过侧向加速度反馈进行补偿。

复盘经验:频域分析能同时评估多个频率点的性能,特别适合处理多模态耦合问题。但Bode图无法直接显示时域响应,需结合时域仿真验证。

工程经验:在飞控律设计初期,建议先使用Bode图确定回路带宽和稳定裕度,再用根轨迹法微调增益,最后通过非线性仿真验证。这种“频域定边界、根轨迹选增益、时域做验证”的流程能显著提高设计效率。

4.4 常见误区与总结复盘

误区一:过度依赖根轨迹而忽视频域裕度

根轨迹法关注极点位置,但无法直接反映系统的相位裕度和幅值裕度。例如,某设计虽然闭环极点位于左半平面,但开环Bode图在穿越频率附近相位下降过快,导致实际相位裕度不足20°,系统对参数摄动极为敏感。

误区二:忽略执行机构与传感器的动态特性

飞控设计中的传递函数往往只包含飞机本体,但舵机的饱和、速率限制、延迟以及传感器的噪声、量化误差都会显著影响系统性能。建议在分析中加入以下环节:

舵机模型:G_act(s) = ω_act² / (s² + 2·ζ_act·ω_act·s + ω_act²)
传感器模型:G_sensor(s) = e^(-τ_s·s) / (T_s·s + 1)

误区三:线性分析结果直接用于非线性系统

飞机在大迎角、高滚转速率等非线性区域,气动导数变化剧烈。线性分析只能作为初步设计依据,必须通过六自由度非线性仿真进行最终验证。

总结复盘

工具优势局限飞控中的最佳用途
传递函数物理意义清晰,便于参数化仅适用于线性时不变系统建立飞控设计模型
根轨迹直观显示增益变化对极点的影响无法直接评估频域裕度确定反馈增益范围
Bode图同时展示幅频/相频特性无法直接显示时域响应评估稳定裕度、设计滤波器
Nyquist曲线适用于时滞系统与多回路分析绘制复杂,不如Bode图直观验证时滞系统的稳定性

经典控制理论在飞控律设计中的核心价值在于提供了一套从物理模型到控制参数的完整映射方法。传递函数将气动数据转化为可计算的数学模型,根轨迹法指导增益选择,Bode图与Nyquist判据确保系统具有足够的稳定裕度。在实际工程中,这四种工具往往需要组合使用,并辅以非线性仿真验证,才能设计出满足飞行品质规范的控制律。