基于奇异摄动分解的弹性飞机乘坐品质控制.pdf

年 月第 卷第 期西 北 工 业 大 学 学 报 ：收稿日期：基金项目：国家自然科学基金（）资助作者简介：李钊星（），西北工业大学博士研究生，主要从事飞行器控制研究。通信作者：许斌（），西北工业大学教授，主要从事智能控制、先进导航研究。：基于奇异摄动分解的弹性飞机乘坐品质控制李钊星，王霞，王敏文，沈健，许斌（西北工业大学 自动化学院，陕西 西安；航空工业西安飞行自动控制研究所，陕西 西安）摘 要：考虑弹性飞机在风干扰作用下乘客乘坐舒适性降低，提出一种基于奇异摄动分解的乘坐品质控制方法。针对弹性客机动力学模型，采用奇异摄动理论将模型解耦为刚性慢变子系统和弹性快变子系统。考虑刚性子系统的附加时变扰动和气动不确定影响，设计干扰观测器估计扰动项并采用神经网络处理模型不确定性，将复合估计信息作为前馈补偿并结合俯仰角速率和法向过载跟踪误差反馈给出自适应鲁棒控制策略。针对弹性子系统设计非奇异终端滑模控制，实现振动主动抑制。综合刚弹子系统控制器设计形成弹性客机乘坐品质控制律，实现附加法向过载和弹性模态的快速抑制收敛。基于李雅普诺夫稳定性分析证明了系统的一致终值有界。仿真结果表明，在离散突风和大气紊流作用下所提出的控制方法能够降低客机关键位置处的附加法向过载，有效提升弹性客机乘坐品质。关 键 词：弹性飞机；奇异摄动分解；干扰观测器；非奇异终端滑模；乘坐品质控制中图分类号：文献标志码：文章编号：（）为实现节省燃油的经济性目标，现代客机呈现轻结构、大柔性和低阻尼等特点，气动弹性问题日益明显。客机飞行过程中不可避免遭遇不同类型、强度和方向的风场，造成机体弯曲振动、客机颠簸和晃动加剧，降低乘客乘坐舒适性，引发驾驶员操纵失误，对飞行安全构成较大隐患。为降低风场扰动对客机飞行安全与机上乘员的乘坐舒适性影响，早期研究由于工程技术限制，通常只能采取添加机械弹簧、阻尼器和惯性装置等结构被动抵抗风干扰。部分学者提出使用机械网络式减震器、流体制动器等部件实现客机乘坐品质改善。随着航空技术发展，借助前后缘襟翼、扰流板等非常规控制面提升乘坐舒适性成为当前工程实际中较为常见的一种方式。上述操纵装置通常还承担增加升力或减缓飞行速度等作用，响应较慢且不能根据风场特性自适应偏转。针对已有研究的局限性，部分学者基于常规控制舵面改善客机乘坐品质。文献借助非定常气动力有理函数拟合方法建立时域连续阵风响应方程，通过控制副翼和升降舵抑制飞机颠簸。文献基于状态空间涡格法建立飞行器阵风气动力模型，并通过鲁棒控制器设计副翼偏转，有效降低气动弹性系统对阵风扰动的敏感性。文献针对翼身融合布局客机操纵性差、抗干扰能力弱等问题，设计自抗扰控制器控制升降舵，改善乘客乘坐舒适性。考虑客机实际飞行过程风扰动的影响，部分研究基于传感器和观测器获得风扰动信息，并通过控制舵面自适应调节，提升乘坐品质。文献利用多普勒激光雷达传感器测量机头处的风场轮廓，据此设计 最优控制提升飞机乘坐品质。文献考虑风切变效应建立干扰动力学模型，基于法向过载反馈设计鲁棒控制器抑制风扰动。文献利用观测器估计风扰动产生的弯矩，设计 控制器减缓风扰动产生的附加机翼载荷。文献通过观测器估计阵风扰动信息并设计前馈控制器，与翼根载荷反馈控制器结合组成乘坐品质控制系统。部分学者考虑气动弹性影响，开展柔性客机乘坐品质研究。文献建立柔性机翼有限元模型并根据输出反馈和状态观测器进行控制器设计。文献基于阵风扰动下的柔性飞机状态空间方程，设西 北 工 业 大 学 学 报第 卷计预测控制律改善乘坐品质。文献考虑飞行过程中参数摄动及阵风扰动对弹性机翼气动特性的影响，建立参数不确定性的线性分式变换模型，并设计鲁棒控制器抑制机翼颤振。基于以上分析可知，已有研究忽略弹性模态在风扰动时带来的附加结构载荷，仅考虑刚体飞机进行乘坐品质控制器设计，部分研究考虑弹性振动影响建立刚弹一体化模型，仅依赖模型本身的稳定裕度使系统趋于稳定，未设计主动反馈控制抑制弹性模态。本文建立大型客机刚弹一体化模型，设计基于奇异摄动分解的乘坐品质控制方法，该方法将客机模型解耦为刚性慢变子系统和弹性快变子系统，并针对 个子系统分别设计控制器，通过二者共同作用改善弹性客机乘坐品质。针对刚性子系统，考虑风干扰引起的附加过载扰动和气动不确定，采用神经网络处理气动不确定性，设计扰动观测器估计风干扰，结合法向过载和俯仰角速率跟踪误差设计自适应鲁棒控制，针对弹性子系统设计非奇异终端滑模鲁棒控制抑制振动模态。仿真结果表明，在离散突风和大气紊流作用下所提方法能够降低飞机关键位置处的附加法向过载，实现振动主动抑制并有效提升弹性客机乘坐品质。问题描述 大型弹性客机动力学模型考虑如下大型弹性客机动力学模型 （）（）（）（）（）该模型包含 阶弹性模态 ，个状态 ，和 个控制输入 ，表示弹性模态广义坐标，和 分别表示飞行速度在平均体轴系中三轴方向上的投影，和 分别表示飞机在平均体轴系下绕三轴转动的角速率，和 分别表示飞机升降舵，副翼和方向舵，表示飞机质量，和分别表示飞机在平均体轴系下绕三轴的转动惯量和交叉惯性积，和 分别表示广义阻尼矩阵、广义刚度矩阵和弹性模态广义力，和（，）分别表示平均体轴系中三轴方向的合外力和力矩，具体表达式为（，）（，），式中：和（，）分别表示弹性模态阻尼比和自然频率；为飞机迎角；（）为动压；为空气密度；为飞机飞行速度，为机翼展长；为平均气动弦长；为机翼面积；，和 均为弹性模态广义坐标对应的非定常气动力影响系数矩阵；为气动力滞后根的状态变量对应的影响系数矩阵；为气动力滞后根对角阵；为弹性模态广义坐标导数系数矩阵；和（，）均表示气动参数；和（，）分别表示飞机平均体轴的气动力与力矩分量；和（，）分别表示飞机平均体轴的弹性模态广义力与力矩分量。动力学模型变换考虑弹性客机的法向过载为（）（）进一步可知第 期李钊星，等：基于奇异摄动分解的弹性飞机乘坐品质控制（）（）定义 ，并考虑弹性客机受到风场扰动后时变干扰 和，可知 ，（）式中：，其余变量表达式为 (（）)(（）（）)（）（）（）（）式中，和，分别表示 中，和对应的弹性模态广义力与力矩系数。奇异摄动分解弹性客机动力学模型定义 ，和，其中是的最小特征根。则（）式变换为 ，（）（，）（）令 ，则气动力滞后根的状态变量 迭代后，其值保持为。忽略弹性模态对刚性慢变量的影响，则（）式可近似写为 （）式中，。则刚性慢变子系统变换为 （）（）（）式可近似写为 （）式中 (（）)(（）（）)（）（）定义快变变量 和 ，新的时间尺度 ，快变子系统可写为 （）式中，控制目标基于奇异摄动分解后的刚性子系统和弹性子系西 北 工 业 大 学 学 报第 卷统，本文控制目标是对刚性慢变子系统设计自适应鲁棒控制，对弹性快变子系统设计非奇异终端滑模控制，综合两者设计得到乘坐品质控制律，降低风扰动下机头、机身和机尾位置处的附加法向过载，提升弹性客机乘坐品质。风扰动及乘坐品质评价准则典型风场扰动包括离散突风和大气紊流。采用国军标中的 离散突风和大气紊流模型模拟风干扰影响。离散突风模型为 （）式中：为风速；为风速的最大值；为离散突风的尺度；为突风空间变量。冯卡门大气紊流模型的谱密度函数为（）（）（）（）（）（）式中：为角频率；，和 为紊流尺度；，和为紊流速度的均方值。采用阵风减缓率、乘坐品质准则 指标以及乘坐不舒适指标作为乘坐品质评价准则。阵风减缓率定义为 （）（）式中：为不加入乘坐品质控制输出的过载响应均方根值；为加入乘坐品质控制后输出的过载响应均方根值。准则的表达式为 （）式中：为乘坐舒适性指数；和 为大气紊流引起的乘员所在位置处的法向过载均方根和侧向过载均方根。乘坐舒适性指数 的取值与乘坐品质的关系为：，乘坐品质为舒适；，乘坐品质为中等；，乘坐品质为不舒适；，乘坐品质为非常不舒适。国军标中乘坐不舒适性指标的表达式为（）（）（）（）式中：为乘坐不舒适指标（法向或横向）；（）为加速度加权函数（法向或横向）；（）为乘客位置处的传递特性；（）为冯卡门阵风强度功率谱密度；为客机弹性模态频率；为截断频率，当 高于截断频率时，在紊流中气动弹性响应不明显。乘坐不舒适指标的评价准则如表 所示。表 乘坐不舒适指标乘坐不舒适指标飞行阶段时间 均方根紊流强度超过概率长航时要求 以上短航时要求 以下 乘坐品质控制器设计针对奇异摄动分解得到的刚性子系统和弹性子系统分别进行控制器设计，并将控制输入叠加，共同组成乘坐品质控制的舵面偏转。刚性慢变子系统自适应鲁棒控制步骤 定义法向过载跟踪误差 ，为期望法向过载，利用神经网络逼近，可得法向过载误差动力学为 （）设计俯仰角速率虚拟控制 为 （）式中：为设计参数；（）为 的估计值；为最优权重向量 的估计值；为 的估计值。令虚拟控制量 通过一阶滤波器，可得 （）式中：（）（），为设计参数。定义俯仰角速率跟踪误差为 ，法向过载跟踪误差 的导数为第 期李钊星，等：基于奇异摄动分解的弹性飞机乘坐品质控制 （）式中：（），。为消除 的影响，设计补偿信号为 （）式中：（），在步骤 给出。定义补偿后的跟踪误差 和 ，则 的导数为 （）设计预测误差 为 （）式中，为设计参数。设计自适应更新律为（）（）（）式中，和 为设计参数。设计扰动观测器为（）（）（）式中：为设计参数。步骤 利用神经网络逼近，俯仰角速率误差动力学为 （）设计刚体慢变子系统升降舵偏转 为 （）式中：为设计参数；（）为的估计值；为最优权重向量 的估计值；为 的估计值。俯仰角速率跟踪误差 的导数为 （）式中：，。设计补偿信号为 （）式中，（）。补偿后的跟踪误差 的导数为 （）设计预测误差 为 （）式中：为设计参数。设计自适应更新律为（）（）（）式中，、和 为设计参数。设计扰动观测器为（）（）（）式中，为设计参数。弹性快变子系统非奇异终端滑模控制针对弹性快变子系统（）式，设计滑模面为 （）式中：为设计参数；和均为正奇数且满足 。滑模面误差动力学为 （）（）式中，。设计弹性快变子系统舵面偏转 为（）（）（）式中：为设计参数。基于刚性慢变子系统舵面偏转 和弹性快变子系统舵面偏转，得到乘坐品质改善控制系统的舵面偏转 。稳定性分析定理 考虑刚体慢变子系统（）式，设计升降舵偏转（）式，自适应更新律（）及（）式，扰动观测器（）及（）式，则李雅普诺夫函数（见（）式中的信号，是一致终值有界的。证明 选择李雅普诺夫函数 为 （）式中（）（）则 的导数为西 北 工 业 大 学 学 报第 卷 （）考虑以下不等式，式中：（）；（）；为常数；为常数。则（）式可进一步写为 ，（）（）式中 ，（）选择控制参数使得 和 ，可得 （）式中，。进一步可得 （）（）由（）式可知，李雅普诺夫函数（）中的信号是一致终值有界的。定理 考虑弹性快变子系统（）式，设计升降舵偏转（）式，则李雅普诺夫函数（见（）式）中的信号 是渐进稳定的。证明 选择李雅普诺夫函数 为（）对 求导可得（）（）（）（）（）（）当，由于 ，因此 ，故 是渐近稳定的，故 和 有界，进一步可知弹性模态 有界。仿真分析算例飞机为某大型双发民用运输客机，其空载质量为 ，绕质心俯仰转动惯量为 ，和 。以机头为参考坐标系原点，位置特性的主要参数如表 所示；各阶弹性模态的阻尼比为，自然频率如表 所示。在进行仿真测试时，算例飞机的初始状态设置为：速度 ，高度为 ，配平攻角为。舵机特性为（），舵面偏转速率限制为 ，位置限制为。表 算例飞机的位置特性主要参数机头驾驶舱坐标机身坐标机翼翼根坐标机尾坐标表 算例飞机的 阶弹性模态自然频率序号频率 模态说明机翼对称一弯机翼反对称一弯机身垂直一弯机身侧向一弯机翼对称二弯机翼面内反对称一弯机翼面内对称一弯平尾反对称一弯机翼面内对称二弯平尾对称一弯机翼面内反对称二弯高阶模态高阶模态参考国军标，设定 离散阵风的阵风尺度为 ，阵风强度为 ；设定大气紊流沿机体轴三轴的风干扰强度参数分别选取轻扰动和中等第 期李钊星，等：基于奇异摄动分解的弹性飞机乘坐品质控制扰动，轻扰动的风干扰强度为 ，紊流尺度为 ；中等扰动风干扰强度为 ，紊流尺度为。仿真步长和仿真周期分别设置为 和 。仿真过程中，飞机从初始时刻即遭遇风扰动。对于弹性快变子系统，控制器参数设置为 ，和。对于刚性慢变子系统，控制器参数设置为 和 ，一阶滤波器参数设置为，更新律参数设置为 ，和，扰动观测器参数设置为 和。为验证所提方法的有效性，本文所提基于奇异摄动分解的乘坐品质控制方法记为“”，而传统只反馈刚体状态，不反馈弹性模态的乘坐品质控制方法记为“”，不加入乘坐品质控制记为“”。客机处于阵风尺度为 的 离散阵风下，飞机机头、机身、机翼翼根和机尾处的过载响应结果如图 所示。图 客机关键位置处的法向过载（阵风）由仿真结果可知，种乘坐品质控制策略均能够减缓风扰动时客机的过载，针对客机各个位置的法向过载峰值，“”方法的减缓效果分别为，和。“”方法的减缓效果为、和，对阵风减缓效果更好。客机处于轻扰动和中等扰动大气紊流干扰下，客机机头、机身、机翼翼根和机尾处的过载响应结如图 和图 所示。图 客机关键位置处的法向过载（轻扰动大气紊流）西 北 工 业 大 学 学 报第 卷图 客机关键位置处的法向过载（中等扰动大气紊流）由仿真结果可知，“”方法相较于“”方法能够使客机关键位置处的过载更加平稳。轻扰动大气紊流引起的客机机头、机身、机翼翼根和机尾处的阵风减缓率如表 所示。由表 可知，当客机受到轻扰动大气紊流后，在客机机头、质心、机尾和机翼翼根处，采用“”方法的减缓率相较于“”方法均有所提升，在机尾位置处减缓效果提升最大，达到。表 大气紊流（轻扰动）减缓率控制器机头过载机身过载机尾过载机翼过载中等扰动大气紊流引起的客机机头、机身、机翼翼根和机尾处的阵风减缓率如表 所示。由表 可知，当客机受到中等扰动大气紊流后，在客机机头、质心、机尾和机翼翼根处，采用“”方法的减缓率相较于“”方法均有所提升，在机身位置处减缓率提升最大，达到。表 大气紊流（中等扰动）减缓率控制器机头过载机身过载机尾过载机翼过载为准确评价算例飞机的乘坐品质，分别列举无乘坐品质控制、传统“”方法以及本文提出的“”方法进行评估，处于中等大气紊流下客机关键位置处的乘坐品质 准则以及乘坐不舒适指标如表 和表 所示。表 大气紊流（中等扰动）准则控制器机头质心机尾机翼翼根 由表 可得，处于中等扰动大气紊流影响下，无乘坐品质控制时客机乘坐品质 准则为“中等”，采用乘坐品质控制的客机能够改善 准则数值，从“中等”改善为“舒适”，而本文所提的“”方法的 准则指标更小、改善效果更优。表 大气紊流（中等扰动）乘坐不舒适指标控制器机头质心机尾机翼翼根 由表 可得，客机受到中等扰动大气紊流影响时，采用乘坐品质控制的客机能改善乘坐不舒适指标。且本文所提的“”方法相较于“”方法，能将客机关键位置处的乘坐不舒适性指标降低至小于，满足长航时等级一标准。客机处于中等扰动大气紊流时，定义 为俯仰角速率估计值，为真实值，。定义 为 法 向 过 载 估 计 值，为 真 实 值，。第 期李钊星，等：基于奇异摄动分解的弹性飞机乘坐品质控制图 俯仰角速率、法向过载估计与估计误差 通过图 可知，在扰动观测器和神经网络作用下，自适应更新律实现了法向过载和俯仰角速率的准确估计。由图 可知，乘坐品质控制中升降舵偏转处于合理范围之内。图 为部分弹性模态示意图，在“”方法控制下相关弹性模态最终收敛到稳态值。图 升降舵偏转 图 部分弹性模态 综上所述，所提出的基于奇异摄动分解的乘坐品质控制策略能够在风扰动下有效降低弹性客机过载波动，提升客机机头、机身、机尾和机翼翼根位置处乘客的乘坐品质。结 论本文针对大型弹性客机乘坐品质提升研究基于奇异摄动分解的自适应鲁棒控制方法。主要研究内容总结如下：）采用奇异摄动方法对飞机进行刚弹模态解耦，忽略弹性模态广义速度和加速度的影响得到飞机刚性慢变子系统，引入新的时间尺度得到飞机弹性快变子系统。）针对刚性慢变子系统，借助扰动观测器估计附加扰动并采用神经网络处理模型不确定性，结合法向过载和俯仰角速率跟踪误差设计自适应鲁棒控制器，针对弹性快变子系统设计终端滑模控制器，结合 个子系统的控制输入组成乘坐品质控制律。）通过李雅普诺夫定理对系统稳定性进行了分析，仿真结果表明所提方法在客机受到 阵风或中等扰动大气紊流影响下，对关键位置处的阵风减缓效果为 和，弹性客机乘坐品质由“中等”改善为“舒适”，满足乘坐不舒适性指标中的长航时等级一标准。西 北 工 业 大 学 学 报第 卷参考文献：，：杨超，邱祈生，周宜涛，等 飞机阵风响应减缓技术综述 航空学报，（）：，（）：（），（）：，（）：，（）：，（）：吴志刚，陈磊，杨超，等 弹性飞机阵风响应建模与减缓方案设计中国科学：技术科学，（）：，（）：（）杨澜，安朝，谢长川，等 基于状态空间涡格法的阵风减缓分析 北京航空航天大学报，（）：，（）：（）骞恒浩，石鹏飞，王敏文，等 基于自抗扰的翼身融合客机控制器设计兵工自动化，（）：，（）：（），（）：，（）：，（）：，（）：王培涵，吴志刚，杨超，等 一种适用于弹性飞机飞行仿真的补丁方法 航空学报，（）：，（）：（），（）：，：陈洋，王正杰，郭士钧 多控制面柔性翼飞行器阵风减缓研究 北京理工大学学报，（）：，（）：（）刘祥，孙秦 一种弹性机翼的颤振主动抑制与阵风减缓方法西北工业大学学报，（）：，（）：（）许斌，王霞 基于时标分解的弹性高超声速飞行器智能控制 航空学报，（）：，（）：（）第 期李钊星，等：基于奇异摄动分解的弹性飞机乘坐品质控制 王立波，荆志伟，唐矗 波阵风中的弹性飞机动力学建模与仿真 航空学报（）：国防科学技术工业委员会 有人驾驶飞机（固定翼）飞行品质 ，孙权志 某大型飞机在大气紊流作用下乘座品质响应计算飞行力学，（）：，（）：（）国防科学技术工业委员会 有人驾驶飞机飞行控制系统通用规范 ，李玉忍，张宏宇，田梁波，等 无人机全电刹车系统级联快速终端滑模控制 西北工业大学学报，（）：，（）：（），；，：，：；引用格式：李钊星，王霞，王敏文，等 基于奇异摄动分解的弹性飞机乘坐品质控制 西北工业大学学报，（）：，（）：（）（：），