基于降阶模型的弹性飞行器QFT控制器设计与仿真
来源:保捱科技网
2010年12月 西北工业大学学报 Journal of Northwestern PolyteehnlcM University Dee.2010 第28卷第6期 Vo1.28 No.6 基于降阶模型的弹性飞行器Qrr控制器设计与仿真 李爱军,孟文杰,徐丽娜,刘世民 (西北工业大学自动化学院,陕西西安710072) 摘要:弹性飞行器具有很高的阶次,直接设计控制器会很困难,且不利于工程实现。为此,文章研究 了基于降阶模型的控制器设计方法。以某大型弹性飞行器12阶模型为研究对象,首先利用平衡截断 降阶方法对该模型降阶,所得到6阶的降阶模型能在较宽频段范围内近似全阶模型。然后,基于降阶 模型,使用Qrr方法设计飞行器的高度保持控制器,直接用于12阶模型的控制。给出了系统高度响 应、鲁棒性验证以及系统状态和四阶弹性变形模态的仿真结果曲线。仿真结果表明,无论是在标称状 态,还是当系统矩阵摄动的状态,高度响应都能快速跟踪指令信号,且无稳态误差,调节时间约为7 s, 系统的各阶结构弹性变形模态都得到了快速抑制,最终衰减为零。 关键词:弹性飞行器;模型降阶;平衡截断;定量反馈;飞行控制 中图分类号:V249 文献标识码:A 文章编号:1000-2758(2010)06-916-05 现代大型高速飞行器的结构弹性变形日趋严 在Nicho ̄图上形成模版和边界,以名义对象满足边 重,此类高速飞行器的设计要求结构质量轻,具有薄 界条件为要求进行控制器设计 J。该方法可以很 升力面、细长机身和低动压载荷系数等特点,这使得 好地实现对具有较大干扰和不确定性系统的控制。 飞行器结构柔性增加,结构变形增大。这种考虑了 定量反馈理论在飞行控制系统设计中已取得了一些 结构弹性变形的飞行器称为弹性飞行器。由于飞行 应用,但目前利用该方法设计弹性飞行器飞行控制 器的弹性运动与刚性运动之间会出现较强的耦合, 系统的研究还很少。本文以某弹性飞行器纵向模型 传统的针对刚性和弹性模态分别设计控制器的方法 为研究对象,首先利用平衡截断降阶方法进行模型 变得不再适用…。另外弹性飞行器在建模时需要 降阶,然后基于降阶后的模型,利用定量反馈理论方 引入大量的自由度(状态变量)来表示刚性模态和 法设计弹性飞行器的高度保持控制系统。仿真结果 弹性模态,因此飞行器模型往往具有很高的阶次,直 表明,利用该方法设计的控制器能很好地抑制飞行 接基于高阶模型进行控制器设计会很困难,同时设 器结构弹性变形模态的振荡,取得了满意的时域效 计所得的控制器阶次较高,不利于工程实现。因此 果,具有较强的鲁棒性。 研究适于弹性飞行器的模型降阶技术显得很有意 义。另外,由于建模误差、模型降阶以及干扰等不确 定性因素的影响,要求弹性飞行器控制器具有较好 1弹性飞行器模型降阶 的鲁棒性。 定量反馈理论(Quantitative Feedbaek Theory。 1.1平衡截断降阶算法 1981年,Moore[4 利用系统可控、可观的概念, Qrr)是以色列学者Horowitz提出的鲁棒控制系统 设计理论,同时也是一种针对具有高度不确定性反 最先提出平衡降阶方法,即所谓系统的内平衡实现 馈控制系统的工程设计理论 】。它将经典控制理 理论。考虑线性定常系统 =Ax+Bu 论频域校正思想推广到不确定系统的鲁棒控制器设 . 计中,将不确定范围和系统性能指标以定量的方式 y=Cx (1) I趺稿日期:2009-12-21 基金项目:航空科学基金(20080553018、20090753007)资助 作者简介:李爱军(1972一),西北工业大学副教授、博士,主要从事飞行器控制与仿真的研究。 第6期 李爱军等:基于降阶模型的弹性飞行器Qrr控制器设计与仿真 式中, ∈R , ∈R ,Y∈R ,假定系统可控、可观、 渐进稳定,则其可控、可观Grammian.阵是正定的 :自由度。以此可知,原模型降为6阶,降阶模型 与 全阶模型G的Bode频域响应比较如图1所示,由图 可知,平衡降阶所得的6阶模型能在较宽频带范围 很好的近似全阶模型。 I eA'BB eA dt =l e^ c c『e^‘dt (2) 且满足下面的李亚普诺夫方程 W +A +BB =0 WoA+A + C=0 (3) 如果 =wo,且为对角阵,则系统为平衡系统,通 常这个条件难以满足。Moore提出:当Ⅳc、W 正定 时,总存在一个非奇异线性变换阵 ,使变换后的实 现(TAT一,TB,CT )平衡。 文献[5]给出了r阵的简便计算方法: (1)对 、 矩阵分别进行Cholesky分解,即 w =L。L:,w =L。L: (2)求Cholesky分解因子L L。乘积的奇异值 分解,即 LT。L。= (3)求得非奇异线性变换阵 T=A一 UrLT。,T~=L VA一 对平衡实现后的系统可以做如下形式的分解 a o J( ̄i)+ ∽㈩ y(t)=【C c: (5) 则( 。,c )即为平衡降阶后所得的降阶 模型。 1.2弹性飞行器的模型降阶 以某高弹性飞行器纵向运动方程作为研究对 象。其基准飞行状态为高度h=1 500 nl,马赫数为 0.6的定常直线平飞。状态 =[11,,a, ,q,,7 , ],i =1,2,3,4, 为飞行器前向速度, 是攻角,0是俯 仰角,q是俯仰角速度, ,面 分别表示4阶广义位移 坐标及其速度,输人 是鸭翼偏转角 ,输出Y是驾 驶舱处的俯仰角速度qqDo系统状态方程形式如方程 (1)所示。 求得系统Hank ̄l奇异值为10.376,10.103 7, 1.977,1.898 2,0.942 4,0.738 4,0.545 4,0.527 9, 0.024 9,0.018,0.003 5,0.003 2,可见最后4个状态 的Hankel奇异值非常小,可以截断而不影响模型精 度,但此时降阶后的模型为8阶,阶次较高,降阶的 意义不大。比较实用的是保留奇异值较大的前6个 (a】频率/rad.s一 图1全阶模型与6阶降阶模型比较 2基于降阶模型的QFT控制器设计 2.1定量反馈理论的基本原理 Qrr通常研究的是一个二自由度的控制结构 (2-D结构),如图2所示。通常用单位反馈、控制器 和前置滤波器来减少被控对象因参数和扰动所造成 的变化。图中C和F为所需设计的控制器和前置滤 波器,d 和d2是外部干扰输人,屁和l,分别为系统的 输入和输出,P为被控对象,有着不确定的参数,且 不确定范围已知。 图2二自由度(2-D结构)控制系统结构框图 Qgr设计目标是:对于在一定范围内变化的P, 找到G和F,使得系统在指令输入R和干扰输人d 作用下闭环响应落在给定的指标范围内。下面 简单总结Qr'r的设计过程: 首先选择设计的频率点集和标称对象(选取那 些能够表征对象模版最大范围的频率特征点即 可),根据不确定对象在所选频率点集处的频率响 应来生成对象模版,根据模版和设计要求(通常包 括跟踪性能指标、抗干扰指标、鲁棒性能要求等,以 西北工业大学学报 第28卷 时域或频域的形式给出)在Nichols图上生成标称对 象的开环频率响应的各个设计边界,并取不同边界 的交集得到复合频域边界;接着进行回路整形过程, 即在Nichols图上通过加入零极点、增益来调整基准 求一致,所以二者可以设计。为保证要求的稳 定裕度,系统开环频率响应必须位于鲁棒稳定边界 之外,同时,为保证跟踪性能要求,在各个频率点上 的开环频率响应必须位于相应频率点处的跟踪边界 之上。在设计过程中,考虑到满足条件的情况下,系 统带宽应尽量小,所以G的增益不能太大,这就要 求开环响应曲线应尽可能的向边界靠拢。图3给出 了控制器作用下的开环频率响应。 对象的开环频率响应曲线的形状和位置,使其满足 设计边界要求,从而得到控制器G;然后根据闭环系 统频率响应要求进行前置滤波器F的设计。由于 所选择的频率点是有限的,因此最后必须加以性能 的分析验证。 2.2 QFT控制器设计 (1)计算出被控模型传函为 bss +b4s4+b3s +b2sz+bl +bo n,、 ,一s6+口5s5+口4s4+口3s3+口2s2+口1s+口o (6) (2)设定性能指标:QFT设计目标是设计控制 使得具有不确定性的系统稳定并且达到一定的性能 要求,闭环系统性能指标有: a)跟踪指标定义了由不确定性和干扰引起的 系统闭环跟踪响应的容许范围, (如)≤ ( ) ≤ ,式中 (s)为闭环传递函数, (s)和 (s) 分别为跟踪边界的上下限。对于本设计,跟踪指标 边界选为 )=筹 (7) (s) (8) b)稳定裕度指标定义为Ii_ 粉I≤p =1.3,它所对应的最小幅值裕度为 =1+l = 1.769,最小相位裕度为:咖村=180。一arccos(0.5 一1)=45.2。。 c)输入抗干扰指标定义为 II口1tJW,l T/・ 、ll ≤ :0’V 。. 5 (、 ,9) (3)控制器设计 由于QFT是一种图形设计法,所以控制器G的 设计过程叫做整形过程(Loop Shaping)。选择频率 集∞={0.01,1,6,7,8,10,25,5O,8O},然后得到这 些频率点的“对象模板”,并结合该模板和定义的性 能指标得到各性能边界,取它们的交集得到综合边 界。由Q171'理论知,控制器G的作用是减小被控对 象由于参数变化造成的对系统性能的影响,滤波器 F的作用是使系统的闭环响应曲线的形状和设计要 相位/(。) 图3控制器作用下的开环频率响应 如果仅考虑G的作用,系统闭环频率特性并不 能满足设计要求,因此必须设计前置滤波器,使闭环 频率特性曲线满足跟踪性能指标要求。前置滤波器 F设计方法与G频率整形过程基本相同。图4为该 滤波器作用下的闭环频率响应。 曼 墨 图4滤波器作用下闭环频率响应 对于本系统,设计得到的C(s)和F(s)如下 G(s)= 621.095(s+30.63)(s +0.01s+0.0001) (s +8.16e一5s+0.0001)(5 +146.2s+9361) (10) 舞 (11) 第6期 李爱军等:基于降阶 型 弹 壁 墨 量 塞 2 : 竺: 3仿真分析与鲁棒性验证 本文以基于降阶后的6阶模型设计的QFT控 制器直接对l2阶模型进行控制,设计其高度保持系 统。仿真结果如图5所示。由图可知,无论是在标 称状态,还是当A阵摄动±20%的状态,高度日都 一0 一 一2 O.1 ・0 一一0.1 0 一 O.1 O.1 O lO 2O t/s O.05 一 能快速跟踪指令信号,且无稳态误差,调节时间约为 0.05 7 s。表明利用QFT方法设计的控制器具有很好的 鲁棒性。 图5高度响应曲线及鲁棒性验证 图6、图7分别给出了标称系统的刚性、结构弹 性变形状态响应曲线。图7中由上至下依次为1— 4阶弹性变形模态,由图7可知,系统的各阶结构弹 性变形模态都得到了快速抑制,一 一最终衰减为零。 一 叠 图6刚性状态响应曲线 o 0 一 0 一 一O.O5 一0.05 0 10 20 t/s t/s (a)模态偏转位移 (b) 模态偏转速度 图7弹性状态响应曲线 4结论 由于考虑了弹性模态,弹性飞行器模型往往具 有高阶性,直接基于高阶模型进行分析与设计比较 困难,因此需要对模型进行降阶处理。平衡截断降 阶算法简单,所得降阶模型能在较宽的频段范围内 很好地近似全阶模型。本文利用该方法对某弹性飞 行器进行降阶,取得了很好的效果,证明该方法适于 弹性飞行器模型降阶。 QFT方法作为一种鲁棒控制系统设计方法,能 同时考虑对象的不确定性和系统的性能指标要求。 本文利用定量反馈理论方法对某弹性飞行器进行纵 向高度保持控制律设计。仿真结果表明,利用该方 法设计的控制器,鲁棒性、抗干扰性强,可以很好地 抑制飞行器的结构弹性变形模态的振荡,取得了很 好的时域控制效果。 I= 1J 1J 西1J 北工业1J 大学1J 学报 第28卷 参考文献: Schmidt D K,Raney D L.Modeling and Simulation of Flexible Flight Vehicles.Journal of Gmdance,Control and Dynamics, 2001,24(3):225~235 王增会,陈增强,孙青林,袁著祉.定量反馈理论发展综述.控制理论与应用,2006,23(3):403—408 Wang ̄nghui,Chen Zengqing,Sun Qianglin,Yuan Zhuzhi.Survey of the DeveIopment for Quantitative Feedback Theory. Control heTory&Applications,2006,23(3):403-408(in Chinese) 韦巍,吴树范,沈勇璋.定量反馈理论在鲁棒与容错飞行控制系中的应用.南京航空航天大学学报,2001,33(2):130 一l34 Wei Wei,Wu Shufn,ShenYoangzhang.Applicaiton ofQrrinRobust/FaultTolerant Fli出ControlSystemsDesign.Journalof Nanjing Uuniversity of Aeronautics&Astronautics,2001,33(2):130~134(in Chinese) Moore B C.Principal Component Analysis in Linear System:Contolrlability,Observability,and Model Reduction.IEEE Trans on Automatic Control,1981,26(1):17—3l Laub A J,Heath M T,Paige C C,Ward R C.Computation of Balnciang Transformaitons and Other Applieations ofSimultane- OIlS Diagonalization Algorihts.ImEEE Trans Automatic Control,1987,32(2):115—122 Applying Quantitative Feedback Theory(QFT)to Designing Controller for Flexible Flight Vehicle Li Aijun,Meng Wenjie,Xu Lina,Liu Shimin (Department ofAutomatic Control,Northwestern Polytechnicl aUniversiy,Xit all 710072,China) Abstract:Aim. I1le introduction of the full paper reviews Refs.1 through 3 and hen tproposes applying orr to de— signing controller.Sections 1 and 2 explain our controller desin metghod,which is illustrated by taking a certain type of lfexible flihtg vehicle as example.Their core consists 0f:(1)we employ the balanced truncation reduction metIId too reduce the number of orders of tl1e flexible flihtg vehicle model from 12 orders to 6;(2)the 6.order model approaches the performance of its full—order model in a wide frequency band range;(3)using the reduced. order model,we apply the Qrr to designing the longitudinal altitude—holding controller,which is diectrly used to control the 12一order model of lfexible flihtg vehicle.Section 3 simulates our design met}ld;tohe simulation resets。 presented in Figs.5,6 and 7,and their analysis show preliminarily that:(1)the coinnland sinalgs ofthe contro1. 1er hus desitned Cagn be tracked quickly by ist alittude response without steady-state error not only in normal state but also in system matrix perturbed satte;(2)the tuning time is about 7s;(3)the elasitc structure rtnsaform mo. dais of all orders Can be restrained quickly and become zero finally. Key words:design,control,flexible flight vehicle,order reduction,balanced truncation,quantiattive feedback hteory(Qrr),flight control