近空间可变翼飞行器具有强非线性、快时变、强耦合以及不确定性等动态特性,给飞行器的控制系统设计带来了困难[1-2]。顾臣风等提出针对近空间飞行器爬升段的跟踪控制研究,采用反馈线性化的方式得到飞行器速度和高度的非线性方程[3]。吴雨珊等采用基于动态逆-PID的控制方式实现了对近空间飞行器的巡航段的姿态控制[4]。针对近空间可变翼飞行器不确定性的动态特点,焦鑫等提出二型模糊滑模控制的方法、自适应鲁棒控制的方法和基于灰色预测滑模控制方式,设计了鲁棒控制器,实现了对飞行器的稳定控制[5-7]。
滑动模态控制相对于传统控制方法具有快速响应、对应参数变化及扰动不灵敏,无需系统在线辨识等优点[8-9]。虽然滑模控制方法具有众多优点,但是传统滑模控制方法仍然存在收敛速度慢,抖振严重等问题。梅红等提出的双幂次趋近律滑模控制方法,解决了收敛速度慢的方法,但是仍存在抖振的问题[10]。神经网络具有大规模并行处理、分布式信息存储、良好的自组织自学习能力,能够充分逼近复杂的非线性关系,具有很好的鲁棒性和容错能力,可以增强控制中控制器的适应能力[11-12]。
针对近空间可变翼飞行器小翼伸缩过程中系统结构和参数存在不确定性问题,本文设计了快速双幂次趋近律滑模控制和神经网络结合的自适应控制器,保证了小翼伸缩切换过程稳定。
1 近空间可变翼飞行器数学模型近空间可变翼飞行器是一个极其复杂的系统,包括自身气动布局、发动机模型、空气动力学模型、运动学模型等。在进行飞行器建模的研究时,不仅要考虑飞行器自身的气动特性、动态结构、弹性形变和气动热等因素,也要考虑大气环境的影响等。这对于飞行器的精确建模来说,十分复杂,难以实现。飞行器纵向模型主要由以下几个方面组成。
1.1 气动力模型飞行器受到的升力和阻力计算公式如下
$ L = \frac{1}{2}\rho {V^2}s{C_L} $ | (1) |
$ D = \frac{1}{2}\rho {V^2}s{C_D} $ | (2) |
式中:ρ为大气密度,V为速度,s为机翼面积,CL为升力系数,CD为阻力系数。
1.2 推力模型近空间可变翼飞行器的发动机推力模型等效为二阶模型[13]:
$ \dot {\beta = - 2\zeta \omega \dot {\beta -}} {\omega ^2}\beta + {\omega ^2}{\beta _c} $ | (3) |
式中:β为发动机节流阀调定值,ζ为阻尼比,ω为固有频率。
飞行器推力计算公式如下
$ {C_T} = \left\{ \begin{array}{l} 0.025\beta ,&\beta < 1\\ 0.022 + 0.003\beta ,&\beta \ge 1 \end{array} \right. $ | (4) |
$ T = \frac{1}{2}\rho {V^2}s{C_T} $ | (5) |
式中:CT为发动机推力系数,T为发动机推力。
1.3 小翼伸缩模型近空间可变翼飞行器在验证仿真中小翼伸出和小翼收回的过程不是瞬变的,本文将其描述不同状态下机翼面积变化,小翼伸缩模型如下
$ \left\{ \begin{array}{l} s = 369,\;\;\;t \le 25\\ s = 369\exp \left( { - 2t + 50} \right) + \\ \;\;\;\;\;389\left( {1 - \exp \left( { - 2t + 50} \right)} \right),\;\;\;25 < t < 30\\ s = 389,\;\;\;\;t \ge 30 \end{array} \right. $ | (6) |
小翼在25 s时开始伸出,伸出前s为369 m2,在30 s时完成小翼的伸出,小翼完全伸出s为389 m2。小翼伸出过程中s变化影响飞行器升力系数CL、阻力系数CD和俯仰力矩系数CM发生变化。
在飞行器小翼伸缩过程中,由于机翼面积发生变化,使得飞行器的气动参数发生变化,存在参数的摄动问题,传统的控制器对这种变化的适应能力有限,针对此问题本文设计了小翼伸缩自适应滑模控制器。
1.4 动力学模型$ \dot V = \frac{{T\cos \alpha - D}}{m} - \frac{\mu }{{{r^2}}}\sin \gamma $ | (7) |
$ \dot \gamma = \frac{{L + T\sin \alpha }}{{mV}} - \frac{{\left( {\mu - {V^2}r} \right)\cos \gamma }}{{V{r^2}}} $ | (8) |
$ \dot q = \frac{{{M_y}}}{{{I_y}}} $ | (9) |
$ \dot \alpha = q - \dot \gamma $ | (10) |
$ \dot h = V\sin \gamma $ | (11) |
$ \dot m = - \frac{T}{{{I_{sp}}}} $ | (12) |
式中α、γ、q、h、m、μ、r、My、Iy分别为飞行器的迎角、俯仰角、俯仰角速率、高度、质量、地球引力常量、地球半径、俯仰力矩和转动惯量。
2 自适应滑模控制器设计针对近空间可变翼飞行器小翼伸缩过程设计的自适应滑模控制器由滑模控制器和更新滑模趋近律参数的BP神经网络两部分构成。自适应滑模控制器采用快速双幂次滑模趋近律。其结构如图 1所示。
近空间可变翼飞行器的高度非线性和严重耦合性导致其控制系统的设计存在较大难度。基于微分同胚理论、李导数和李括号等非线性系统的基本理论,对近空间飞行器模型的飞行速度V和飞行高度h进行输入输出精确反馈线性化处理,既可以保证模型的准确度又可以降低设计难度。
自适应滑模控制器设计的目标是通过设计合适的控制输入量δe、βc,使控制系统在小翼伸出变化过程中飞行器的速度V和飞行高度h能够良好的跟踪设定的指令信号Vd和hd。
2.1 飞行器纵向模型精确反馈线性化本文针对近空间可变翼飞行器的纵向非线性模型(7)~(10),对飞行速度V和飞行高度h进行精确反馈线性化处理,得到:
$ \left[ \begin{array}{l} {V'''}\\ {h^{\left( 4 \right)}} \end{array} \right] = \left[ \begin{array}{l} {{V'''}_0}\\ h_0^{\left( 4 \right)} \end{array} \right] + \left[ {\begin{array}{*{20}{c}} {{b_{11}}}&{{b_{12}}}\\ {{b_{21}}}&{{b_{22}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{\beta _c}}\\ {{\delta _e}} \end{array}} \right] $ | (13) |
其中
$ {b_{11}} = \frac{{\rho {V^2}S{c_\beta }\omega _n^2}}{{2m}}\cos \alpha $ | (14) |
$ {b_{12}} = - \frac{{{c_e}\rho {V^2}S{c_A}}}{{2m{I_y}}}\left( {T\sin \alpha + {D_\alpha }} \right) $ | (15) |
$ {b_{21}} = \frac{{\rho {V^2}S{c_\beta }\omega _n^2}}{{2m}}\sin \left( {\alpha + \gamma } \right) $ | (16) |
$ \begin{array}{l} {b_{22}} = \frac{{{c_e}\rho {V^2}S{c_A}}}{{2m{I_y}}}\left[ {T\cos \left( {\alpha + \gamma } \right) + {L_\alpha }\cos \gamma + } \right.\\ \;\;\;\;\;\;\;\;\left. {{T_\alpha }\sin \left( {\alpha + \gamma } \right) - {D_\alpha }\sin \gamma } \right] \end{array} $ | (17) |
式中:cA为平均气动弦长,ce和cβ为与δe和βc相关的启动系数,Iy为转动惯量。
2.2 快速双幂次趋近律滑模控制律设计飞行器实际飞行与指令信号的跟踪误差表示为
$ {e_V}\left( t \right) = V\left( t \right) - {V_d}\left( t \right) $ | (18) |
$ {e_h}\left( t \right) = h\left( t \right) - {h_d}\left( t \right) $ | (19) |
针对近空间可变翼飞行器纵向非线性模型选择适当的积分滑模面[4]:
$ \mathit{\boldsymbol{S}} = \left[ {\begin{array}{*{20}{c}} {{S_V}}\\ {{S_h}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{{\left( {{\rm{d}}/{\rm{d}}t + {\lambda _V}} \right)}^3}\int {{e_V}\left( \tau \right){\rm{d}}\tau } }\\ {{{\left( {{\rm{d}}/{\rm{d}}t + {\lambda _h}} \right)}^4}\int {{e_h}\left( \tau \right){\rm{d}}\tau } } \end{array}} \right] $ | (20) |
式中:λv、λh为待设计正常数,积分项用于抑制稳态误差。
设计的积分滑模面S对时间t求导:
$ \mathit{\boldsymbol{S'}} = \left[ {\begin{array}{*{20}{c}} {{{S'}_V}}\\ {{{S'}_h}} \end{array}} \right] = \left[ \begin{array}{l} {{e'''}_v} + 3{\lambda _V}{{e'}_V} + 3\lambda _V^2{{e'}_V} + \lambda _V^3{e_V}\\ e_h^{\left( 4 \right)} + 4{\lambda _h}{{e'''}_h} + 6\lambda _h^2{{e'}_h} + 4\lambda _h^3{{e'}_h} + \lambda _h^2{e_h} \end{array} \right] $ | (21) |
对指令跟踪误差进行求导:
$ {{e'''}_h}\left( t \right) = V'''\left( t \right) - {{V'''}_d}\left( t \right) $ | (22) |
$ {{e'''}_h}\left( t \right) = h{\left( t \right)^{\left( 4 \right)}} - h_d^{\left( 4 \right)}\left( t \right) $ | (23) |
并将其简化为
$ {\mathit{\boldsymbol{u}}_{eq}} = {\mathit{\boldsymbol{B}}^{\rm{T}}}{\left( {\mathit{\boldsymbol{B}}{\mathit{\boldsymbol{B}}^{\rm{T}}}} \right)^{ - 1}}\left( { - \left[ \begin{array}{l} {{V'''}_0}\\ h_0^{\left( 4 \right)} \end{array} \right] + \left[ \begin{array}{l} {{V'''}_d}\\ h_d^{\left( 4 \right)} \end{array} \right] - \left[ {\begin{array}{*{20}{c}} {{F_V}}\\ {{F_h}} \end{array}} \right]} \right) $ | (24) |
设计快速双幂次趋近律:
$ \begin{array}{l} \mathit{\boldsymbol{S'}} = \\ \left[ {\begin{array}{*{20}{c}} { - {k_{V1}}{{\left| {{S_V}} \right|}^{{\beta _V}}}{\mathop{\rm sgn}} \left( {{S_V}} \right) - {k_{V2}}{{\left| {{S_V}} \right|}^{{\alpha _V}}}{\mathop{\rm sgn}} \left( {{S_V}} \right) - {k_{V3}}{S_V}}\\ { - {k_{h1}}{{\left| {{S_h}} \right|}^{{\beta _h}}}{\mathop{\rm sgn}} \left( {{S_h}} \right) - {k_{h2}}{{\left| {{S_h}} \right|}^{{\alpha _h}}}{\mathop{\rm sgn}} \left( {{S_h}} \right) - {k_{h3}}{S_h}} \end{array}} \right] \end{array} $ | (25) |
式中:kh1、kh2、kh3、kV1、kV2、kV3>0,0<αV, αh<1, βV、βh>1, sgn(SV)、sgn(Sh)为符号函数。
滑模控制器由两部分构成,分别为等效控制项(24) 和快速双幂次趋近律控制项(25),可以得到快速双幂次趋近滑模控制律为
$ \begin{array}{l} \mathit{\boldsymbol{u}} = {\mathit{\boldsymbol{B}}^{\rm{T}}}{\left( {\mathit{\boldsymbol{B}}{\mathit{\boldsymbol{B}}^{\rm{T}}}} \right)^{ - 1}}\left( { - \left[ \begin{array}{l} {{V'''}_0}\\ h_0^{\left( 4 \right)} \end{array} \right] + \left[ \begin{array}{l} {{V'''}_d}\\ h_d^{\left( 4 \right)} \end{array} \right] - \left[ {\begin{array}{*{20}{c}} {{F_V}}\\ {{F_h}} \end{array}} \right]} \right) + \\ \left[ {\begin{array}{*{20}{c}} { - {k_{V1}}{{\left| {{S_V}} \right|}^{{\beta _V}}}{\mathop{\rm sgn}} \left( {{S_V}} \right) - {k_{V2}}{{\left| {{S_V}} \right|}^{{\alpha _V}}}{\mathop{\rm sgn}} \left( {{S_V}} \right) - {k_{V3}}{S_V}}\\ { - {k_{h1}}{{\left| {{S_h}} \right|}^{{\beta _h}}}{\mathop{\rm sgn}} \left( {{S_h}} \right) - {k_{h2}}{{\left| {{S_h}} \right|}^{{\alpha _h}}}{\mathop{\rm sgn}} \left( {{S_h}} \right) - {k_{V3}}{S_h}} \end{array}} \right] \end{array} $ | (26) |
本节对设计的快速双幂次趋近律滑模控制的稳定性进行验证,首先定义李亚普诺夫函数:
$ V = \frac{1}{2}{\mathit{\boldsymbol{S}}^{\rm{T}}}\mathit{\boldsymbol{S}} $ | (27) |
函数对时间t的导数为
$ \begin{array}{*{20}{c}} {V' = {S_V}{{{\dot S}'}_V} + {S_h}{{{\dot S}'}_h} = }\\ { - {k_{V1}}{{\left| {{S_V}} \right|}^{{\beta _V}}}{S_V}{\mathop{\rm sgn}} \left( {{S_V}} \right) - {k_{V2}}{{\left| {{S_V}} \right|}^{{\alpha _V}}}{S_V}{\mathop{\rm sgn}} \left( {{S_V}} \right) - }\\ {{k_{V3}}{S_V}{S_V} - {k_{h1}}{{\left| {{S_h}} \right|}^{{\beta _h}}}{S_h}{\mathop{\rm sgn}} \left( {{S_h}} \right) - }\\ {{k_{h2}}{{\left| {{S_h}} \right|}^{{\alpha _h}}}{S_h}{\mathop{\rm sgn}} \left( {{S_h}} \right) - {k_{V3}}{S_h}{S_h} = }\\ { - {k_{V1}}{{\left| {{S_V}} \right|}^{{\beta _{V + 1}}}} - {k_{V2}}{{\left| {{S_V}} \right|}^{{\alpha _{V + 1}}}} - {k_{V3}}S_V^2 - }\\ {{k_{h1}}{{\left| {{S_h}} \right|}^{{\beta _{h + 1}}}} - {k_{h2}}{{\left| {{S_h}} \right|}^{{\alpha _{h + 1}}}} - {k_{V3}}S_h^2} \end{array} $ | (28) |
当kV1、kV2、kV3、kh1、kh2、kh3为正常数时,V′≤0,证明系统状态够到达切换面,满足系统稳定条件。
2.4 神经网络设计近空间可变翼飞行器小翼伸出过程中,飞行器气动参数产生摄动。研究小翼伸缩过程的自适应控制,利用BP神经网络对小翼伸缩过程中所对应的滑模趋近律参数进行学习、训练。在小翼伸出过程中实时更新不同状态下的趋近律参数。
BP神经网络结构如图 2所示。
图 2中,X为BP神经网络的输入为近空间可变翼飞行器小翼伸缩过程中的机翼面积SW、迎角α、马赫数Ma,通过训练后网络的输出o为同步更新小翼伸缩过程中快速双幂次滑模趋近律的系数av、bv、ah、bh、kv1、kv2、kv3、kh1、kh2、kh3。
3 仿真验证为了验证设计的控制器能够稳控制近空间可变翼飞行器小翼伸缩过程,仿真条件及参数设置如下:飞行器模型在25 s时小翼由收回状态开始伸出,30 s时小翼完全伸出,机翼面积由收回状态下s=369变为s=389。小翼伸出过程中,气动参数发生变化,且小翼伸缩过程具有不确定性。
1) 小翼收回模态下速度和高度初值为
$ {V_{ro}} = 4590\;{\rm{m}}/{\rm{s}},{h_{ro}} = 33\;528\;{\rm{m}} $ |
指令信号为
$ {V_{rd}} = 4630\;{\rm{m}}/{\rm{s}},{h_{rd}} = 33\;588\;{\rm{m}} $ |
2) 小翼伸出模态下速度和高度初值为
$ {V_{so}} = 4630\;{\rm{m}}/{\rm{s}},{h_{so}} = 33\;588\;{\rm{m}} $ |
指令信号为
$ {V_{sd}} = 4590\;{\rm{m}}/{\rm{s}},{h_{rd}} = 33\;528\;{\rm{m}} $ |
仿真验证对比传统滑模控制和自适应滑模控制两种小翼伸出过程的控制效果。图 3、4给出了飞行器速度和高度响应曲线,可以看出自适应滑模控制具有较好的快速性和抑制抖动的能力,上升时间叫滑模控制提升50%。图 5、6给出了飞行器的状态航迹角和迎角以的变化曲线,从中可以看出自适应滑模控制下航迹角的变化较滑模控制减小30%,迎角和俯仰角变化率抖动小,表现出飞行器的多面抖动也小。图 7、8给出控制输入舵面偏转角以及油门开合度的变化曲线,自适应滑模控制的舵偏量为滑模控制的50%,且油门开合未达到饱和,滑模控制的油门开合饱和,说明自适应滑模控制速度和高度跟踪性能更优。
近空间可变翼飞行器在小翼伸出过程中具有参数不确定性,为了验证快速双幂次趋近律滑模控制和神经网络结合的自适应滑模控制器的鲁棒性,若机翼浸润面积Sw、大气密度ρ、俯仰力矩惯性积、俯仰力矩系数中的参数ce、翼弦长c存在摄动,表示为
$ \begin{array}{*{20}{c}} {\left| {\Delta {I_y}} \right|/{I_{y0}} = 0.05,\left| {\Delta c} \right|/{c_0} = 0.05}\\ {\left| {\Delta {S_w}} \right|/{S_{w0}} = 0.05,\left| {\Delta \rho } \right|/{\rho _0} = 0.05}\\ {\left| {\Delta {c_e}} \right|/{c_{e0}} = 0.05} \end{array} $ |
在数值仿真过程中加入不确定参数摄动,分析快速双幂次趋近律滑模控制和神经网络结合的自适应滑模控制器的鲁棒性。
与未加参数摄动的控制输入对比分析,结果如图 9~14所示,在加入参数摄动的小翼伸出过程中,速度和高度同未加参数摄动的情况基本相同,升降舵偏角发生微小的变化,但与整个偏转幅度相比可以忽略不计。因此,针对小翼伸缩过程设计的快速双幂次趋近律滑模控制和神经网络结合的自适应滑模控制器在小翼伸缩过程存在参数摄动的情况下具有较好的控制效果,鲁棒性较强。
1) 快速双幂次趋近律滑模控制和神经网络结合的自适应滑模控制器能更好的抑制升降舵偏角在小翼伸缩瞬间的跳变,具有更好的跟踪精度,并且具有较强的鲁棒性。
2) 使用快速双幂次滑模趋近律实现滑模面的快速稳定收敛。
3) 使用神经网络更新滑模控制律参数。未来的研究中可针对近空间飞行器模型存在的弹性摄动进行深入研究。
[1] |
吴宏鑫, 孟斌. 高超声速飞行器控制研究综述[J]. 力学进展, 2009, 39(6): 756-765. WU Hongxin, MENG Bin. Summary of control research on hypersonic vehicle[J]. Advances in mechanics, 2009, 39(6): 756-765. DOI:10.6052/1000-0992-2009-6-J2009-088 (0) |
[2] |
方存光, 孙勇, 王伟. 高超声速飞行器及其飞行状态控制[J]. 控制工程, 2008(S1): 21-23. FANG Cunguang, SUN Yong, WANG Wei. Hypersonic vehicle and its flight state control[J]. Control engineering, 2008(S1): 21-23. (0) |
[3] |
顾臣风, 江驹, 吴雨珊. 近空间飞行器爬升跟踪控制[J]. 哈尔滨工程大学学报, 2016, 36(11): 1526-1531. GU Chenfeng, JIANG Ju, WU Yushan. A near space vehicle ascent tracking control[J]. Journal of Harbin Engineering University, 2016, 36(11): 1526-1531. (0) |
[4] |
吴雨珊, 江驹, 甄子洋, 等. 基于动态逆-PID的高超声速飞行器巡航姿态控制[J]. 飞行力学, 2015(2): 145-149. WU Yushan, JIANG JU, ZHEN Ziyang, et al. hypersonic cruise attitude based on-PID[J]. Control of dynamic inverse flight dynamics, 2015(2): 145-149. (0) |
[5] | JIAO Xin, BARIS F, JIANG Ju, et al. Adaptive mode switching of hypersonic morphing aircraft based on type-2 TSK fuzzy sliding mode control[J]. Sciece China information sciences, 2015, 58(7): 1-15. (0) |
[6] |
焦鑫, 江驹. 非线性系统自适应鲁棒控制器设计[J]. 哈尔滨工程大学学报, 2016(3): 402-407. JIAO Xin, JIANG Zi. Design of adaptive robust controller for nonlinear system[J]. Journal of Harbin Engineering University, 2016(3): 402-407. (0) |
[7] |
焦鑫, 江驹, 孟宏鹏. 基于灰色预测滑模控制的鲁棒飞行控制器设计[J]. 飞行力学, 2015(6): 542-546. JIAO Xin, JIANG Zi, MENG Hongpeng. Robust flight controller design based on grey prediction sliding mode control[J]. Flight mechanics, 2015(6): 542-546. (0) |
[8] | SUN H, LI S, SUN C. Finite time integral sliding mode control of hypersonic vehicles[J]. Nonlinear dynamics, 2013, 73(1-2): 229-244. DOI:10.1007/s11071-013-0780-4 (0) |
[9] | YU S, YU X, SHIRINZADEH B, et al. Continuous finitetime control for robotic manipulators with terminal sliding mode[J]. Automatica, 2005, 41(11): 1957-1964. DOI:10.1016/j.automatica.2005.07.001 (0) |
[10] | MEI H, WANG Y. Fast convergent sliding mode variable structure control of robot[J]. Information and control, 2009, 38(5): 552-557. (0) |
[11] |
朱坚民, 沈正强, 李孝茹, 等. 基于神经网络反馈补偿控制的磁悬浮球位置控制[J]. 仪器仪表学报, 2014, 35(5): 976-986. ZHU Jianmin, SHEN Zhengqiang, LI Xiaoru, et al. Position control of magnetic levitation ball based on neural network feedback compensation control[J]. Journal of instrumentation, 2014, 35(5): 976-986. (0) |
[12] |
吴玉香, 王聪. 不确定机器人的自适应神经网络控制与学习[J]. 控制理论与应用, 2013, 30(8): 990-997. WU Yuxiang, WANG Cong. Adaptive neural network control and learning for uncertain robots.[J]. Control theory and application, 2013, 30(8): 990-997. (0) |
[13] |
丁洪波, 蔡洪. 基于特征参数的吸气式高超声速飞行器上升段轨迹优化[J]. 固体火箭技术, 2012, 35(3): 296-300. DING Hongbo, CAI Hong. Optimization of ascent trajectory of an air breathing hypersonic vehicle based on characteristic parameters[J]. Solid rocket technology, 2012, 35(3): 296-300. (0) |
[14] | XU H, MIRMIRANI M D, LOANNOU P A. Adaptive sliding mode control design for a hypersonic flight vehicle[J]. Journal of guidance, control, and dynamics, 2004, 27(5): 829-838. DOI:10.2514/1.12596 (0) |
[15] | WANG Q, STENGEL R F. Robust nonlinear control of a hypersonic aircraft[J]. Journal of guidance control & dynamics, 2000, 23(23): 577-585. (0) |