Open Access
Issue
JNWPU
Volume 44, Number 1, February 2026
Page(s) 134 - 142
DOI https://doi.org/10.1051/jnwpu/20264410134
Published online 27 April 2026

© 2026 Journal of Northwestern Polytechnical University. All rights reserved.

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

随着航空技术的发展, 非常规布局诸如翼身融合、联翼布局、倾转旋翼等逐渐应用于民用航空器领域, 这些新布局普遍具有高冗余操纵效应器/操纵面的新颖特征。

如何在概念设计阶段就充分考虑安全性和适航性的要求成为非常规布局成功进入民用领域的关键问题之一。虽然当前有效的适航规章针对的是传统布局航空器, 但其对航空器操纵面设计提出约束的底层安全逻辑不会改变。多操纵面适航性评估的传统方法依赖大量时域仿真, 导致计算成本高、迭代周期长, 且难以在早期识别物理层设计缺陷(如操纵力矩冗余不足), 导致后期设计迭代的风险增加[1]。转矩可达集(attainable moment subset, AMS)[2-7]能直观展示操纵力矩, 在飞行控制领域得到广泛应用。然而, 现有研究尚未建立AMS与适航设计的联系, 特别是缺乏面向飞机概念设计阶段的快速评价方法。

针对上述问题, 本文提出一种基于AMS的多操纵面飞机机动性适航快速评价方法, 通过AMS进行飞机物理架构的筛选, 构建“设计-评价-迭代-仿真”的评价流程。该方法首先将适航条款对机动性的要求(滚转、俯仰、偏航)转化为适航约束AMS并作为基准; 对比初始设计的AMS, 若其覆盖适航约束则进一步评价典型故障模式(如舵面卡阻、丧失等)的剩余AMS覆盖性; 任一阶段不满足要求则对操纵面进行迭代设计, 通过评价则进行时域仿真。本文以多操纵面翼身融合飞机为算例, 仿真实验表明, 本方法较传统方法其优势体现在: 通过AMS进行正常与故障状态机动性适航评价, 从力矩角度快速评价当前设计是否满足适航要求, 提升了设计与评估的效率。研究成果为多控制效应器/多操纵面民用航空器在概念设计阶段的适航考虑提供了支撑。

1 转矩可达集

转矩可达集的概念最早由Durham在1993年首先提出[8-10]。最初是为了研究多操纵面设计下,针对控制分配问题如何确定控制边界,他将飞机的操纵面u(t)(umin(t)≤u(t)≤umax(t))在控制效率矩阵Be: RmRk(m>k)映射下形成的几何体Φ称为转矩可达集。具体而言, AMS是由滚转、俯仰和偏航力矩构成的凸多面体, 分配系统只能对AMS内的转矩向量进行有效分配和控制。

1.1 转矩可达集的求解方法

转矩可达集是由操纵面的物理约束和控制效率矩阵共同确定, 设bkuk的效率向量, 那么控制效率矩阵Be改写为: Mathematical equation。求解Φ时, 首先调整控制效率矩阵Be的任意一轴与一个面垂直, 这里引入变换矩阵T, 控制效率矩阵Be中的第i行和第j列为bibj

Mathematical equation(1)

Mathematical equation(2)

Mathematical equation(3)

求出T使其满足

Mathematical equation(4)

Mathematical equation(5)

求解该不定方程即可求出变换矩阵T, 令t13为任意非0值, 即可求解T的另外2个元素, 令t13=1, 则(5)式变为

Mathematical equation(6)

将(6)式转换为矩阵方形式

Mathematical equation(7)

若满足

Mathematical equation(8)

Mathematical equation(9)

Mathematical equation(10)

则说明控制效率矩阵Be不满足任意3×3矩阵满秩, 则(9)式的求解方法失效, 这种情况可在Be中添加小扰动量, 破坏向量的相关性, 从而使用(9)式求解。

此时已经求得T= [t11 t12 t13] 的值, 将TB相乘后, B(1, i)和B(1, j)中将有1个元素为0, 而B的第一行中其他元素的符号决定了其他操纵面的设置方式: 若该元素为正, 则对应的操纵面应设置为最大值; 若该元素为负, 则设置为最小值。通过这种方式, 可以生成一个转矩可达集Φ的边界面。反之, 若将正元素对应的操纵面设置为最小值, 负元素对应的操纵面设置为最大值, 则可生成另一个边界面。为了完整计算Φ的所有边界面, 需对每2个控制执行机构的组合重复上述计算。对于m个控制执行机构和3个力矩轴, 所需计算次数为

Mathematical equation(11)

式中, m是控制执行机构的数目[11]。

1.2 操纵面故障情况的转矩可达集求解

对于操纵面故障情况, 常见的故障模式包括操纵面卡阻、操纵面松浮和操纵面损伤[12]。

操纵面卡阻是指操纵面固定在某一位置, 无法产生预期的控制效果, 同时可能引入额外的附加力和力矩。在求解转矩可达集时, 操纵面卡阻的情况可视为控制效率不变, 若卡阻在δ位置, 则操纵面位置限制变为umin=umax=δ。

操纵面松浮是指操纵面失效, 也就是操纵面的控制效率为0, 可以简化看作松浮的操纵面不会对飞机产生任何气动力和气动力矩, 即飞机的控制输入产生的作用为0。求解转矩可达集时,假设第i片操纵面松浮,控制效率矩阵Be对应的故障位置变为

Mathematical equation(12)

操纵面损伤, 会导致飞机质量、质心、转动惯量以及飞机气动系数等参数的变化, 假设操纵面质量远远小于整个飞机的质量, 则可以视为上述参数的变化可忽略。此时, 操纵面产生的控制效率会降低, 假设第i片操纵面正常状态下控制效能为δ, 操纵面损伤d%时, 实际有效的控制效能为δ′=(1-d%)δ, 求解转矩可达集时, 控制效率矩阵Be对应的故障位置变为

Mathematical equation(13)

2 机动性快速评价方法

2.1 评价流程

根据以上内容, 建立下述评价流程:

1) 初始构型设计: 获取设计方案的相关数据, 如: 质量、气动外型、动力装置、操纵面设计等;

2) 适航约束建模: 将适航条款的机动性要求转化为力矩约束, 建立适航约束的AMS;

3) 正常状态评价: 计算飞机初始设计AMS, 评价其是否覆盖适航约束AMS;

4) 故障状态评价: 确定需要评价的故障清单, 分别计算故障状态剩余AMS并评价其是否覆盖适航约束AMS;

5) AMS机动性适航评价: 若步骤3)~4)中任一评价未通过, 则说明飞机设计方案不符合适航要求并修改设计方案; 全部通过则进入步骤6);

6) 时域仿真: 对通过AMS评价的设计方案进行时域仿真, 通过可得到符合适航要求的飞机设计方案。

评价流程图如图 1所示。

Thumbnail: 图1 Refer to the following caption and surrounding text. 图1

评价流程图

2.2 机动性适航要求

根据《运输类飞机适航标准》(CCAR-25-R4)[13]分析得出: 正常状态下, 飞机纵向能产生-1.0g~2.5g过载, 横向操纵性要求飞机在7 s内完成60°的滚转, 航向操纵性要求临界发动机失效后, 可以左/右航向改变15°; 操纵面故障状态下, 飞机纵向能产生0.8g~1.3g过载, 横向操纵性要求飞机11 s内完成60°的滚转, 并且能够30°坡度定常左转弯或右转弯。将上述适航要求转化为操纵力矩要求, 转化到AMS的坐标系中, 定量评价飞机的机动能力是否满足适航要求。

2.3 适航约束的转矩可达集

根据2.2节中适航要求, 将三轴机动性要求转化为AMS。计算适航要求规定动作所需的最小三轴力矩, 并通过力矩映射至AMS坐标系, 构建适航约束AMS作为基准, 实现适航要求与操纵能力的可视化对比分析。

1) 俯仰力矩

假设在目标过载的保持期间, 飞机的状态量不再随时间快速变化, 并且是纯纵向平面的机动, 无侧滑, 气动力和力矩是线性的。所需俯仰力矩Cmreq

Mathematical equation(14)

Mathematical equation(15)

式中: Cm0为零升力俯仰力矩系数;Cma为纵向静稳定导数;αn为目标n过载下需要的迎角;Cmq为俯仰阻尼导数;qn为目标过载机动时的俯仰角速率;为平均几何弦长;V为速度;T为发动机推力;zT为发动机安装位置到飞机水平面的距离;S为机翼参考面积;Q为动压。

2) 滚转力矩

横向操纵要求将操纵时间和偏转角度转化为滚转力矩, 假设该过程为匀速, 滚转角加速度=0, 因此仅考虑滚转阻尼项, 则所需滚转力矩系数Clreq

Mathematical equation(16)

式中: p为滚转角速率;Clp为滚转阻尼导数。

3) 偏航力矩

偏航操纵需满足单发失效时平衡非对称推力偏航力矩的基本要求, 同时保持航向操纵能力, 则所需偏航力矩系数Cnreq

Mathematical equation(17)

式中: T为正常工作的发动机推力;yT为发动机到机体y轴的距离;b为翼展;C为偏航静稳定导数;C为侧力对侧滑角的导数;r为偏航角速率。

3 仿真实验与分析

3.1 飞机模型介绍

本文采用多操纵面翼身融合(blended wing body, BWB)飞机模型作为初始构型, 如图 2所示。

Thumbnail: 图2 Refer to the following caption and surrounding text. 图2

BWB飞机模型

BWB飞机由升降舵、机翼两侧2对升降副翼以及1对垂尾构成[14]。

3.2 正常状态机动性评价

选取飞机在3 000 m高度, Ma=0.3时的控制效率矩阵Be1

Mathematical equation(18)

操纵面位置限制如表 1所示。

表1

位置限制

根据1.1节的方法, 通过初始设计飞机的控制效率矩阵B e1和操纵面位置限制, 求解初始设计AMS, 采用2.3节的方法, 计算出适航要求下的转矩可达集。三轴力矩系数的区间结果为

Mathematical equation

Mathematical equation

将三轴操纵面力矩系数的约束区间绘制在转矩可达集坐标系中。在正常状态下, 判断飞机初始设计AMS是否覆盖正常适航约束AMS。如图 3~4所示。

Thumbnail: 图3 Refer to the following caption and surrounding text. 图3

初始构型AMS与适航约束覆盖性

Thumbnail: 图4 Refer to the following caption and surrounding text. 图4

初始构型AMS与适航约束覆盖性

图 3~4表明, 初始设计AMS在Cm轴未能覆盖适航约束AMS, 说明初始设计在纵向机动上不满足适航要求, 需要对操纵面进行修改。

3.3 操纵面设计修改

由3.2节分析可知, BWB飞机因该构型纵向操纵力矩臂较短, 导致初始构型的纵向机动能力未满足适航要求; 而滚转与偏航方向的操纵能力存在明显冗余。

在操纵面优化设计中, 现有升降舵布局已完全饱和, 但机翼后缘内侧区域仍存在未利用的气动面。为提高纵向操纵效能, 本研究参考主流BWB机型及传统民机的舵面布局[15-17], 提出改进方案:

1) 操纵面面积与数量修改: 增加3号升降副翼, 并延长原有1~2号升降副翼展向长度, 如图 5所示;

Thumbnail: 图5 Refer to the following caption and surrounding text. 图5

操纵面设计修改

2) 气动模型简化: 将舵效简化为与偏转角呈线性关系;

3) 耦合效应简化: 忽略操纵面数量和面积改变对飞行力学耦合的影响。

气动分析得到迭代设计的控制效率矩阵Be2

Mathematical equation(19)

根据Be2进行迭代设计后的机动性评价。

图 6表明, 迭代设计后纵向AMS体积较初始设计扩大1.52倍, 完整覆盖适航约束AMS, 且保持滚转与偏航操纵能力的冗余。

Thumbnail: 图6 Refer to the following caption and surrounding text. 图6

设计修改AMS与适航约束覆盖性

3.4 正常状态机动性评价

AMS评价正常状态符合适航要求, 进行时域仿真, 首先进行控制律与分配律设计。基本控制律如图 7a)~7b)所示, 其中各个增益参数设定依据飞行品质要求, 采用经典控制律方法设计。纵向为俯仰通道和发动机通道, 俯仰通道采取C*控制增稳设计, 发动机通道输入为油门开度; 滚转通道使用滚转角速率反馈降低滚转角振荡, 侧向过载反馈来消除侧滑角; 偏航通道采用偏航角速率以及侧滑角反馈改善荷兰滚模态的阻尼比、消除侧滑角, 并提高航向稳定性[19], 控制分配方法选用全操纵面伪逆法。

Thumbnail: 图7 Refer to the following caption and surrounding text. 图7

BWB飞机飞行控制律示意图

根据2.2节中的机动性适航要求, 正常状态纵向产生-1.0g~2.5g的法向过载, 横航向7 s滚转60°, 选定在高度3 000 m、Ma=0.3时相应的指令信号, 图 8a)~8b)分别为飞机拉-1g与2.5g过载的机动动作, 图 8c)为过载-1g时的操纵面响应。

Thumbnail: 图8 Refer to the following caption and surrounding text. 图8

纵向仿真结果图

图 9a)为飞机做7 s滚转60°机动动作的滚转角响应图, 图 9b)为操纵面响应图。

Thumbnail: 图9 Refer to the following caption and surrounding text. 图9

横向仿真结果图

由仿真结果可得, 该设计方案在正常状态下的AMS与时域仿真评价结果均符合适航要求。

3.5 操纵面故障状态机动性评价

采用2.3节的方法, 计算得到操纵面故障后适航条款所规定的三轴力矩系数范围为

Mathematical equation

1) 纵向评价

在高度3 000 m、Ma=0.3的条件下,对操纵面卡阻的适航规定: 根据0.8g~1.3g过载量的操纵杆机动位置对应的升降舵偏转最大位置,经仿真后选取升降舵急偏后的卡阻位置为下偏9°。

首先进行AMS评价,图 10a)为正常状态AMS与升降舵故障后剩余AMS的对比,剩余AMS缩小至48%,图 10b)表明升降舵故障剩余AMS覆盖适航约束AMS,符合要求并进行时域仿真。

Thumbnail: 图10 Refer to the following caption and surrounding text. 图10

升降舵故障适航评价

本文分别构建升降舵、升降副翼、方向舵故障情况,进行AMS评价与时域仿真。仿真条件设置为:①推杆机动至0.8g的注入故障,在正常平飞状态下,在15 s升降舵突然急偏,并朝下卡阻在8.9°位置处,之后在第20 s进行推杆机动;②拉杆机动至1.3g的注入故障流程与0.8g相同。

图 10c)~10f)分别为推杆机动0.8g,拉杆机动至1.3g时的过载变化及单侧的各个升降副翼的偏角响应。

2) 横向评价

在高度3 000 m、Ma=0.3的条件下,根据操纵面卡阻的适航规定,飞机的仿真评估飞行条件可选取1号升降副翼故障和单侧方向舵故障,经仿真可确定最外侧单侧升降副翼卡阻位置为向下(上)偏7°(方向舵正常),单片方向舵卡阻位置为左(右)偏8°(1号升降副翼正常)。

(1) 1号升降副翼下偏7°卡阻,首先进行AMS评价。由图 11a)为正常状态AMS与升降舵故障后剩余AMS的对比,剩余AMS缩小至60%,图 11b)表明升降舵故障剩余AMS覆盖适航约束AMS,符合要求并进行时域仿真。

Thumbnail: 图11 Refer to the following caption and surrounding text. 图11

1号升降副翼故障适航评价

在5 s时飞机左侧1号升降副翼下偏7°,然后在20 s时飞机建立左滚30°的坡度转弯。图 11c)~11d)为1号升降副翼故障的飞机滚转角和另一侧正常的升降副翼响应。

(2) 单片方向舵左偏8°卡阻,首先进行AMS评价。由图 12a)为正常状态AMS与升降舵故障后剩余AMS的对比,剩余AMS缩小至61%,图 12b)表明方向舵故障剩余AMS覆盖适航约束AMS,符合要求并进行时域仿真。

Thumbnail: 图12 Refer to the following caption and surrounding text. 图12

方向舵故障适航评价

单侧方向舵的故障,其故障注入与升降副翼一致。图 12c)~12d)为单侧方向舵故障时的飞机滚转角和另一侧正常的升降副翼和方向舵偏角响应。

3.6 实验结果分析

正常状态下,AMS与时域仿真结果均符合适航要求,故障状态下,分别验证了3种操纵面卡阻情况,分析如下:

图 10c)~10f)可知,升降舵急偏后卡阻后,飞机过载和舵偏量均较小,此时飞机为了保持稳定飞行,利用其他升降副翼进行配平,20 s推拉杆进行机动完成过载指令,符合适航要求。

图 11c)可知,左侧1号升降副翼卡阻后,滚转角大约发生了3.5°的变化,20 s给定30°坡度定常转弯指令后完成滚转指令,符合适航要求。

图 12c)可知,单侧方向舵卡阻后,滚转角大约发生了2.6°的变化,20 s给定30°坡度定常转弯指令后完成滚转指令,符合适航要求。

通过AMS评价操纵面正常与故障情况和时域仿真结果均符合适航条件,说明改进后的设计方案满足机动性适航要求。

4 结论

本研究将转矩可达集理论与适航要求相结合,研究了飞机在正常和故障状态下的机动性适航快速评价方法,飞行仿真表明了本文提出的方法对于机动性适航评价的有效性。得出以下结论:①本方法改进了传统依赖大量仿真的试错模式,减少冗余的仿真实验;②通过AMS直观展示并量化飞机的三轴操纵能力和操纵面故障后的剩余飞行能力,为设计优化提供依据;③建立了“设计-评价-迭代-仿真”的闭环评价流程,在飞机概念设计阶段为多操纵面飞机的机动性适航评价提供了新的思路,具有工程应用价值。

References

  1. Zhang Shuai, Xia Ming, Zhong Bowen. Evolution and technical factors influencing civil aircraft aerodynamic configuration[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(1): 30–44 (in Chinese) [Google Scholar]
  2. Shen Xiaoming, Lü Xinbo, Huang Zhenwei. Research on the control allocation for aircrafts with multi-control surfaces[J]. Aeronautical Science & Technology, 2016, 27(9): 8–11 (in Chinese) [Google Scholar]
  3. Fan Weiqiang, Li Jian, Zhang Chong, et al. Overall preliminary design of new configuration quadtiltrotor aircraft[J]. Flight Dynamics, 2024, 42(3): 26–32 (in Chinese) [Google Scholar]
  4. Cang Yu. Integrated design and optimization of theaerodynamic shape for aircraft under hybrid-wing-body configuration[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2018. (in Chinese) [Google Scholar]
  5. Durham W C. Attainable moments for the constrained control allocation problem[J]. Journal of Guidance, Control, and Dyna-mics, 1994, 17(6): 1371–1373. [Article] [Google Scholar]
  6. Durham W C. Constrained control allocation-three-moment problem[J]. Journal of Guidance, Control, and Dynamics, 1994, 17(2): 330–336. [Article] [Google Scholar]
  7. Bordingnon K A, Durham W C. Closed-form solutions to constrained control allocation problem[J]. Journal of Guidance, Control, and Dynamics, 1995, 18(5): 1000–1007. [Article] [Google Scholar]
  8. Shi Jingping, Qu Xiaobo, Bi Kejun. Multiple effector control allocation method based on attainable moment subset[J]. Journal of Northwestern Polytechnical University, 2014, 32(1): 123–130. [Article] (in Chinese) [Google Scholar]
  9. Shi Jingping, Qu Xiaobo, Li Guangwen. An effective design of tailored generalizedinverse control allocation method[J]. Journal of Northwestern Polytechnical University, 2013, 31(1): 8–13. [Article] (in Chinese) [Google Scholar]
  10. Shi Jingping, Zhang Weiguo. Efficient optimization method of daisy chain control an allocation based on genetic algorithm[J]. Journal of Northwestern Polytechnical University, 2012, 30(4): 582–588. [Article] (in Chinese) [Google Scholar]
  11. Durham W, Bordignon K A, Beck R. Aircraft control allocation[M]. Chichester, UK: John Wiley & Sons, 2017 [Google Scholar]
  12. Shi Jingping, Qu Xiaobo. Control allocation theory and its application for aircraft with multiple control surfaces[M]. Beijing: National Defense Industry Press, 2017 (in Chinese) [Google Scholar]
  13. 中国民用航空局. 中国民用航空规章-第25部-运输类飞机适航标准[S]. CCAR-25-R4, 2016. [Google Scholar]
  14. Yang Jianzhong, Yin Zhuyu, Yang Shibin. Airworthiness assessment of lateral-directional stability characteristics for blended wing body aircraft[J]. Civil Aircraft Design & Research, 2024(2): 19–26 (in Chinese) [Google Scholar]
  15. Garmendia D C, Chakraborty I, Trawick D R, et al. Assessment of electrically actuated redundant control surface layouts for a hybrid wing body concept[C]//14th AIAA Aviation Technology, Integration, and Operations Conference, 2014 [Google Scholar]
  16. Hileman J I, Spakovszky Z S, Drela M, et al. Airframe design for silent fuel-efficient aircraft[J]. Journal of Aircraft, 2010, 47(3): 956–969. [Article] [Google Scholar]
  17. Risch T, Cosentino G, Regan C, et al. X-48B flight test progress overview[C]//47th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, 2009 [Google Scholar]
  18. Yin Zhuyu. Control allocation design of BWB aircraftand airworthiness stability and handling characteristics evaluation[D]. Tianjin: Civil Aviation University of China, 2024. (in Chinese) [Google Scholar]

All Tables

表1

位置限制

All Figures

Thumbnail: 图1 Refer to the following caption and surrounding text. 图1

评价流程图

In the text
Thumbnail: 图2 Refer to the following caption and surrounding text. 图2

BWB飞机模型

In the text
Thumbnail: 图3 Refer to the following caption and surrounding text. 图3

初始构型AMS与适航约束覆盖性

In the text
Thumbnail: 图4 Refer to the following caption and surrounding text. 图4

初始构型AMS与适航约束覆盖性

In the text
Thumbnail: 图5 Refer to the following caption and surrounding text. 图5

操纵面设计修改

In the text
Thumbnail: 图6 Refer to the following caption and surrounding text. 图6

设计修改AMS与适航约束覆盖性

In the text
Thumbnail: 图7 Refer to the following caption and surrounding text. 图7

BWB飞机飞行控制律示意图

In the text
Thumbnail: 图8 Refer to the following caption and surrounding text. 图8

纵向仿真结果图

In the text
Thumbnail: 图9 Refer to the following caption and surrounding text. 图9

横向仿真结果图

In the text
Thumbnail: 图10 Refer to the following caption and surrounding text. 图10

升降舵故障适航评价

In the text
Thumbnail: 图11 Refer to the following caption and surrounding text. 图11

1号升降副翼故障适航评价

In the text
Thumbnail: 图12 Refer to the following caption and surrounding text. 图12

方向舵故障适航评价

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.