【智能系统】采用智能算法参数整定的磁流变减摆自抗扰控制

团购合买资源类别：文库，文档格式：PDF，文档页数：9，文件大小：5.99MB

第17卷第5期智能系统学报 Vol.17 No.5 2022年9月 CAAI Transactions on Intelligent Systems Sep.2022 D0:10.11992/tis.202109012 网络出版地址：https:/kns.cnki.net/kcms/detail/23.1538.TP.20220617.1915.010.html 采用智能算法参数整定的磁流变减摆自抗扰控制董磊，陈增强2，孙明玮，孙青林，黄朝阳 (1.南开大学人工智能学院，天津300350：2.天津市智能机器人重点实验室，天津300350) 摘要：为了提高飞机起飞和降落过程中的稳定性及安全性，采用磁流变阻尼器进行飞机前轮摆振控制。本文描绘了改进的Bouc-wen模型磁流变阻尼器的力学特性，并将之用于一个三自由度的飞机前起落架模型中，而后为此设计了减摆自抗扰控制器(active disturbance rejection control,.ADRC):然后采用混沌分数阶天牛群算法整定优化自抗扰控制器的参数。最后，仿真结果显示经智能算法优化参数后的减摆自抗扰控制器能够较好地控制摆振现象，并且其性能具有一定的鲁棒性。关键词：飞机前起落架；摆振；磁流变阻尼器；改进的Bouc-wen模型；减摆；自抗扰控制；天牛群算法；混沌分数阶天牛群算法：参数优化：鲁棒性中图分类号：TP273 文献标志码：A文章编号：1673-4785(2022)05-0951-09 中文引用格式：董磊，陈增强，孙明玮，等.采用智能算法参数整定的磁流变减摆自抗扰控制[.智能系统学报，2022,17(5)： 951-959. 英文引用格式：DONG Lei,,CHEN Zengqiang,.SUN Mingwei,,etal.Active disturbance rejection control of the magneto--rheological shimmy damper based on the smart chaos fractional-order beetle swarm optimization algorithm CAAI transactions on intelli- gent systems,2022,17(⑤：951-959. Active disturbance rejection control of the magneto-rheological shimmy damper based on the smart chaos fractional-order beetle swarm optimization algorithm DONG Lei',CHEN Zengqiang",SUN Mingwei',SUN Qinglin',HUANG Zhaoyang' (1.College of Artificial Intelligence,Nankai University,Tianjin 300350,China;2.Key Laboratory of Intelligent Robots,Tianjin 300350,China) Abstract:To improve the stability and safety of an aircraft during the take-off and landing,a magneto-rheological damper is applied to the shimmy control of its front wheel.In this paper,the mechanical properties of a modified Bouc-Wen magneto-rheological damper are described,which is used in an aircraft landing gear model with three de- grees of freedom and for which a third-order active disturbance rejection controller(ADRC)is later designed;sub- sequently,the chaos fractional-order beetle swarm optimization is used to optimize the controller parameters.Finally, the simulation results indicate that the controller with optimized parameters exhibits good performance in shimmy redu- cing control with certain robustness. Keywords:aircraft front landing gear;shimmy;magneto-rheological damper;modified Bouc-wen model;shimmy re- duce;active disturbance rejection control;beetle swarm optimization,chaos fractional order beetle swarm optimization, parameters optimization;robustness 飞机现多采用前三点式起落架，这种结构方尤其是着陆过程，前起落架会受到机身大部分的式可以使飞机具有良好的地面运动稳定性、着陆负荷，加之跑道等因素的相互作用，前轮容易产易操纵性、飞行员视界良好等优点山。但是对前生基于中心的往复运动的摆振现象。这种现象若起落架提出了更高的要求，起飞和着陆过程中，不加以控制，会产生安全事故。常用的方式就是加入阻尼器来控制摆振。应用较多的是油液式阻收稿日期：2021-09-03.网络出版日期：2022-06-21 基金项目：国家自然科学基金项目(61973175,61573177.61973172). 尼器，它可以提供阻尼力来抑制摆振。近些年，通信作者：陈增强.E-mail:chenzq@nankai,.edu.cn. 磁流变液作为智能材料受到广泛的关注，而磁流

DOI: 10.11992/tis.202109012 网络出版地址: https://kns.cnki.net/kcms/detail/23.1538.TP.20220617.1915.010.html 采用智能算法参数整定的磁流变减摆自抗扰控制董磊1 ，陈增强1,2，孙明玮1 ，孙青林1 ，黄朝阳1 （1. 南开大学人工智能学院，天津 300350; 2. 天津市智能机器人重点实验室，天津 300350）摘要：为了提高飞机起飞和降落过程中的稳定性及安全性，采用磁流变阻尼器进行飞机前轮摆振控制。本文描绘了改进的 Bouc-wen 模型磁流变阻尼器的力学特性，并将之用于一个三自由度的飞机前起落架模型中，而后为此设计了减摆自抗扰控制器 (active disturbance rejection control, ADRC)；然后采用混沌分数阶天牛群算法整定优化自抗扰控制器的参数。最后，仿真结果显示经智能算法优化参数后的减摆自抗扰控制器能够较好地控制摆振现象，并且其性能具有一定的鲁棒性。关键词：飞机前起落架；摆振；磁流变阻尼器；改进的 Bouc-wen 模型；减摆；自抗扰控制；天牛群算法；混沌分数阶天牛群算法；参数优化；鲁棒性中图分类号：TP273 文献标志码：A 文章编号：1673−4785(2022)05−0951−09 中文引用格式：董磊, 陈增强, 孙明玮, 等. 采用智能算法参数整定的磁流变减摆自抗扰控制 [J]. 智能系统学报, 2022, 17(5): 951–959. 英文引用格式：DONG Lei, CHEN Zengqiang, SUN Mingwei, et al. Active disturbance rejection control of the magneto-rheological shimmy damper based on the smart chaos fractional-order beetle swarm optimization algorithm[J]. CAAI transactions on intelligent systems, 2022, 17(5): 951–959. Active disturbance rejection control of the magneto-rheological shimmy damper based on the smart chaos fractional-order beetle swarm optimization algorithm DONG Lei1 ，CHEN Zengqiang1,2 ，SUN Mingwei1 ，SUN Qinglin1 ，HUANG Zhaoyang1 (1. College of Artificial Intelligence, Nankai University, Tianjin 300350, China; 2. Key Laboratory of Intelligent Robots, Tianjin 300350, China) Abstract: To improve the stability and safety of an aircraft during the take-off and landing, a magneto –rheological damper is applied to the shimmy control of its front wheel. In this paper, the mechanical properties of a modified Bouc–Wen magneto–rheological damper are described, which is used in an aircraft landing gear model with three degrees of freedom and for which a third-order active disturbance rejection controller (ADRC) is later designed; subsequently, the chaos fractional-order beetle swarm optimization is used to optimize the controller parameters. Finally, the simulation results indicate that the controller with optimized parameters exhibits good performance in shimmy reducing control with certain robustness. Keywords: aircraft front landing gear; shimmy; magneto-rheological damper; modified Bouc-wen model; shimmy reduce; active disturbance rejection control; beetle swarm optimization; chaos fractional order beetle swarm optimization; parameters optimization; robustness 飞机现多采用前三点式起落架，这种结构方式可以使飞机具有良好的地面运动稳定性、着陆易操纵性、飞行员视界良好等优点[1]。但是对前起落架提出了更高的要求，起飞和着陆过程中，尤其是着陆过程，前起落架会受到机身大部分的负荷，加之跑道等因素的相互作用，前轮容易产生基于中心的往复运动的摆振现象。这种现象若不加以控制，会产生安全事故。常用的方式就是加入阻尼器来控制摆振。应用较多的是油液式阻尼器，它可以提供阻尼力来抑制摆振[2]。近些年，磁流变液作为智能材料受到广泛的关注，而磁流收稿日期：2021−09−03. 网络出版日期：2022−06−21. 基金项目：国家自然科学基金项目（61973175，61573177，61973172）. 通信作者：陈增强. E-mail：chenzq@nankai,edu.cn. 第 17 卷第 5 期智能系统学报 Vol.17 No.5 2022 年 9 月 CAAI Transactions on Intelligent Systems Sep. 2022

第17卷智能系统学报 ·952· 变阻尼器则是其重要应用之一。磁流变阻尼器的数，通过智能优化算法来寻找最优解。仿真结果表应用范围较广，文献[3-4]研究了将磁流变阻尼器明所设计控制器能够很好地控制摆振现象，而且应用于汽车悬架及座椅的减震（振）控制中；Dyke 面对参数的摄动，该控制器具有良好的鲁棒性。 SJ等将磁流变阻尼器应用于土木工程建筑物减轻地震伤害中；Nguyen团队I研究了一种磁流 1磁流变阻尼器变弹性隔离器，用于结构体的减振；文献[)中采磁流变阻尼器中的磁流变液具有良好的性能，用磁流变阻尼器进行直升机与地面共振的抑制研正常情况下呈现牛顿流体的特性，若外加磁场后，究。与传统的油液式阻尼器相比，磁流变阻尼器磁流变小颗粒在极短时间内后凝成“固体”，并且具有更快的响应速度，响应时间可以达到15ms 该过程是可逆的。这种特性，使得其具有良好的有的甚至可以达到6.5ms,且具有更大的阻尼力可控性。范围，而且其阻尼力的大小可以通过外加磁场的 1.1 磁流变阻尼器的力学模型方式进行调节，是一种半主动控制。磁流变阻尼为方便对磁流变阻尼器的描述与研究，研究器阻尼力可调的功能，使得一些控制方法能够应人员对其进行数学建模与分析。其中磁流变阻尼用于起落架的减摆控制。Kang等采用天棚阻器常用的模型包括：Bingham模型、Bouc-wen模尼控制对飞机起落架摆振进行控制；文献[9]采型、双曲线回归模型、Dahl模型等6-1。其中前两用自适应控制摆振问题；文献[10]将模糊PD控者是应用最广泛的模型。制应用于磁流变阻尼器控制；文献[11]用神经网本文采用文献[9]中改进的Bouc-wen模型，络逆模型通过阻尼力反解电流的方式实现阻尼器与Bouc-wen模型相比，物理概念更加明显，同时控制；文献[12]通过智能算法优化PD控制器的也能够很好地变现出低速时的滞回特性，模型原参数实现减摆控制。上述控制方法或控制设备较理图如图1所示。复杂或控制精度较差。天棚阻尼控制需要计算阻 Bouc-wen 尼系数非常依赖模型的精确性；自适应控制和神经网络逆模型设计比较复杂，难以实现，而且神经网络逆模型的准确度与训练数据息息相关； PID控制在控制效果上没有自抗扰控制精确1。因此，需要设计一种易于工程上实现又具有良好控制效果及鲁棒性的控制器。本文设计了自抗扰控制器(active disturbance rejection control,.ADRC) 来实现飞机前起落架的减摆控制。ADRC是由韩京清研究员提出的一种新型控制方法。它不依图1改进的Bouc-wen模型原理 Fig.1 Princple map of modified Bouc-wen model 赖对象模型，能够估计扰动并进行补偿，具有较高的控制精度和良好的鲁棒性。ADRC继承和发模型如式(1)所示：扬了PD控制和现代控制理论的优点，在某些领 f=c(i)文+a(i)z (1) 域表现出比PD控制更好的性能，有望能够取代式中：f表示阻尼力；x表示阻尼器的位移；表示 PID控制在工业控制中的地位。但是ADRC需滞回位移，其表达式如式(2)所示：要整定的参数比较多，比较依靠经验，所以这限之=-y闭zlm-1-B闭"+A元 (2) 制了它在实际工程中的推广应用。正因为ADRC c(①、a()分别为磁流变材料的粘性系数和滞参数整定的困难，一些智能优化算法被引入进行回力在阻尼器总阻尼力所占比重的调节参数，取参数整定。本文主要工作为：I)一改进的Bouc- 值为电流的函数。 Wen模型磁流变阻尼器应用到飞机前起落架进行 c()=13.5+2.083i+1.5632-1.302°N.s/mm 减摆控制。2)设计自抗扰控制器：将飞机起落架 a(0=900+2750i+18752-1563N.s/mm 和磁流变减摆器看作一个整体，将摆角作为控制式(2)中y、B、A为模型与滞回位移有关的参目标并进行误差反馈，将系统内部和外部的扰动数，通过调整这些参数可以控制阻尼力变化的平看作总扰动并用扩张状态观测器进行估计，然后稳度。y、B表示滞回环宽度和高度的调节系数，进行补偿。3)采用智能优化算法一混沌分数取值为150；A表示与最大阻尼力相关的参数，取阶天牛群算法整定控制器参数通过设计适应度函值为6：n表示滞回环圆滑系数，取值为2

变阻尼器则是其重要应用之一。磁流变阻尼器的应用范围较广，文献 [3-4] 研究了将磁流变阻尼器应用于汽车悬架及座椅的减震（振）控制中；Dyke S J 等 [5] 将磁流变阻尼器应用于土木工程建筑物减轻地震伤害中；Nguyen 团队[6] 研究了一种磁流变弹性隔离器，用于结构体的减振；文献 [7] 中采用磁流变阻尼器进行直升机与地面共振的抑制研究。与传统的油液式阻尼器相比，磁流变阻尼器具有更快的响应速度，响应时间可以达到 15 ms，有的甚至可以达到 6.5 ms，且具有更大的阻尼力范围，而且其阻尼力的大小可以通过外加磁场的方式进行调节，是一种半主动控制。磁流变阻尼器阻尼力可调的功能，使得一些控制方法能够应用于起落架的减摆控制。Kang 等 [8] 采用天棚阻尼控制对飞机起落架摆振进行控制；文献 [9] 采用自适应控制摆振问题；文献 [10] 将模糊 PID 控制应用于磁流变阻尼器控制；文献 [11] 用神经网络逆模型通过阻尼力反解电流的方式实现阻尼器控制；文献 [12] 通过智能算法优化 PID 控制器的参数实现减摆控制。上述控制方法或控制设备较复杂或控制精度较差。天棚阻尼控制需要计算阻尼系数非常依赖模型的精确性；自适应控制和神经网络逆模型设计比较复杂，难以实现，而且神经网络逆模型的准确度与训练数据息息相关； PID 控制在控制效果上没有自抗扰控制精确[13]。因此，需要设计一种易于工程上实现又具有良好控制效果及鲁棒性的控制器。本文设计了自抗扰控制器（active disturbance rejection control, ADRC）来实现飞机前起落架的减摆控制。ADRC 是由韩京清研究员提出的一种新型控制方法[14]。它不依赖对象模型，能够估计扰动并进行补偿，具有较高的控制精度和良好的鲁棒性。ADRC 继承和发扬了 PID 控制和现代控制理论的优点，在某些领域表现出比 PID 控制更好的性能，有望能够取代 PID 控制在工业控制中的地位[15]。但是 ADRC 需要整定的参数比较多，比较依靠经验，所以这限制了它在实际工程中的推广应用。正因为 ADRC 参数整定的困难，一些智能优化算法被引入进行参数整定。本文主要工作为：1) 一改进的 BoucWen 模型磁流变阻尼器应用到飞机前起落架进行减摆控制。2) 设计自抗扰控制器：将飞机起落架和磁流变减摆器看作一个整体，将摆角作为控制目标并进行误差反馈，将系统内部和外部的扰动看作总扰动并用扩张状态观测器进行估计，然后进行补偿。3) 采用智能优化算法−混沌分数阶天牛群算法整定控制器参数通过设计适应度函数，通过智能优化算法来寻找最优解。仿真结果表明所设计控制器能够很好地控制摆振现象，而且面对参数的摄动，该控制器具有良好的鲁棒性。 1 磁流变阻尼器磁流变阻尼器中的磁流变液具有良好的性能，正常情况下呈现牛顿流体的特性，若外加磁场后，磁流变小颗粒在极短时间内后凝成“固体”，并且该过程是可逆的。这种特性，使得其具有良好的可控性。 1.1 磁流变阻尼器的力学模型为方便对磁流变阻尼器的描述与研究，研究人员对其进行数学建模与分析。其中磁流变阻尼器常用的模型包括：Bingham 模型、Bouc-wen 模型、双曲线回归模型、Dahl 模型等[16-18]。其中前两者是应用最广泛的模型。本文采用文献 [9] 中改进的 Bouc-wen 模型，与 Bouc-wen 模型相比，物理概念更加明显，同时也能够很好地变现出低速时的滞回特性，模型原理图如图 1 所示。 Bouc-wen c (i) x f 图 1 改进的 Bouc-wen 模型原理 Fig. 1 Princple map of modified Bouc-wen model 模型如式 (1) 所示： f = c(i)· x˙ +α(i)·z (1) 式中： f 表示阻尼力；x表示阻尼器的位移； z 表示滞回位移，其表达式如式（2）所示： z˙ = −γ · |x˙| ·z· |z| n−1 −β · |x˙| · |z| n + A· x˙ (2) c(i)、α(i) 分别为磁流变材料的粘性系数和滞回力在阻尼器总阻尼力所占比重的调节参数，取值为电流的函数。 c(i) = 13.5+2.083i+1.563i 2 −1.302i 3 N·s/mm α(i) = 900+2750i+1875i 2 −1563i 3 N·s/mm γ β A γ β A n 式 (2) 中、、为模型与滞回位移有关的参数，通过调整这些参数可以控制阻尼力变化的平稳度。、表示滞回环宽度和高度的调节系数，取值为 150；表示与最大阻尼力相关的参数，取值为 6；表示滞回环圆滑系数，取值为 2。第 17 卷智能系统学报 ·952·

·953· 董磊，等：采用智能算法参数整定的磁流变诚摆自抗扰控制第5期 1.2改进的Bouc-wen模型的力学特性进行该模型的激励测试，输入激励信号幅值 ant 为3mm,频率为5Hz,输入电流一次为0、0.2、0.4、 0.6、0.8、1.0A,得到位移与阻尼力的关系如图1 所示，速度与阻尼力的关系如图2所示。 ESO 2500 2000 图4自抗扰控制器结构 1500 Fig.4 Principle map of ADRC 1000 ESO是自抗扰控制的核心部分，它可以将系 500 统内扰外扰，确定性和不确定性的扰动当作总扰 0 -500 动一同进行估计。 -1000 NLSEF将误差的反馈形式表示为非线性，和 -1500 扰动估计补偿共同作用形成控制信号。 -2000 -2500 2.2自抗扰控制器的设计 3 - -10 1 x/mm 2.2.1飞机前起落架模型飞机前起落架采用一个三自由度的非线性简图2位移与阻尼力的关系 Fig.2 Displacement vs.damping force 易模型2,101，如式(3)(4)所示： I,=M1(p)+M2(p)+M3(a)+M(p/m)+M5(3) 从图2中可以看出，位移与阻尼力之间的关系具有良好的对称性，随着电流的增大，阻尼力增大。 +y=w+e-04 (4) 而且在低速运动中，阻尼力呈现了滞回特性，说式中：M1是轮胎上方扭矩连杆和转向杆一起提供明该模型能够较好地刻画磁流变阻尼器的特性。的线性弹簧扭矩；M2是油气式缓冲器轴承之间的 2 自抗扰控制粘性摩擦和减摆器的组合阻尼力矩；M3、M,分别为由侧滑和横摆引起的轮胎力矩和轮胎阻尼力 2.1自抗扰控制基本原理矩；M是磁流变阻尼器力矩。y是飞机机轮侧向位移；自抗扰的思想最早可以追溯到我国古代指南 v是飞机滑行速度；是机轮摆角。车的设计思想，是一种主动抗扰的方法。它主要 M1~Ms的表达式如式(5)(13)所示，M,是围绕包括跟踪-微分器(tracking differential,.TD)、扩张轮胎中心的扭转力矩；F,是以稳定距为杠杆作用状态观测器(extended state observer,,ESO)、非线性的转弯力矩，其他变量或参数以表1的形式给出。误差反馈率(nonlinear stste error feedback,NLSEF) 表1起落架模型参数及扰动估计补偿等部分，其结构图如图3所示。 Table 1 the parameters of the landing gear model 图中y是系统输出，是控制信号，r是给定信号；参数数值 b是模型参数。惯性矩I/kgm) 2500 垂直负荷FN 4200 2000 1500 线性扭矩系数c/(N.m)-rad -15000 1000 组合阻尼系数k/Nms)rad 500 -35 0 稳定距elm 0.12 -500 轮胎力矩导数c/(mrad -1000 -2 08A -1500 1.0A 侧向力导数CFarad 20 -2000 轮胎力矩极值角a() 10 -2500 100 -60 -2020 60100 轮胎力极值角() v/mm-s) 接触半长alm 0.1 图3速度与阻尼力的关系 Fig.3 Velocity vs.damping force 胎面宽度常数kNm2)rad -126 松弛长度σm 0.3 TD可以最快地跟踪给定信号的方法来提取微分信号，除此之外它还有安排过渡过程、配置滑行速度v(ms) 0-50 系统零极点、数字整流等应用。诚摆力臂dm 0.1

1.2 改进的 Bouc-wen 模型的力学特性进行该模型的激励测试，输入激励信号幅值为 3 mm，频率为 5 Hz，输入电流一次为 0、0.2、0.4、 0.6、0.8、1.0 A，得到位移与阻尼力的关系如图 1 所示，速度与阻尼力的关系如图 2 所示。 −3 −2 −1 0 1 2 3 x/mm −2 500 −2 000 −1 500 −1 000 −500 0 500 1 000 1 500 2 000 2 500 f/N 0 A 0.2 A 0.4 A 0.6 A 0.8 A 1.0 A 图 2 位移与阻尼力的关系 Fig. 2 Displacement vs. damping force 从图 2 中可以看出，位移与阻尼力之间的关系具有良好的对称性，随着电流的增大，阻尼力增大。而且在低速运动中，阻尼力呈现了滞回特性，说明该模型能够较好地刻画磁流变阻尼器的特性。 2 自抗扰控制 2.1 自抗扰控制基本原理 y u r b0 自抗扰的思想最早可以追溯到我国古代指南车的设计思想，是一种主动抗扰的方法。它主要包括跟踪–微分器 (tracking differential, TD)、扩张状态观测器 (extended state observer,ESO)、非线性误差反馈率 (nonlinear stste error feedback, NLSEF) 及扰动估计补偿等部分，其结构图如图 3 所示。图中是系统输出，是控制信号，是给定信号；是模型参数。 v/(mm·s−1) −2 500 −2 000 −1 500 −1 000 −500 0 500 1 000 1 500 2 000 2 500 f/N 0 A 0.2 A 0.4 A 0.6 A 0.8 A 1.0 A −100 −60 −20 20 60 100 图 3 速度与阻尼力的关系 Fig. 3 Velocity vs. damping force TD 可以最快地跟踪给定信号的方法来提取微分信号，除此之外它还有安排过渡过程、配置系统零极点、数字整流等应用。 NLSEF Plant ESO 1/b0 b0 z1 z2 z3 u0 u y r 图 4 自抗扰控制器结构 Fig. 4 Principle map of ADRC ESO 是自抗扰控制的核心部分，它可以将系统内扰外扰，确定性和不确定性的扰动当作总扰动一同进行估计。 NLSEF 将误差的反馈形式表示为非线性，和扰动估计补偿共同作用形成控制信号。 2.2 自抗扰控制器的设计 2.2.1 飞机前起落架模型飞机前起落架采用一个三自由度的非线性简易模型[2,10,19] ，如式 (3)(4) 所示： Izφ¨ = M1(φ)+ M2( ˙φ)+ M3(α)+ M4( ˙φ/v)+ M5 (3) y˙ + v σ y = vφ+(e−a) ˙φ (4) M1 M2 M3 M4 M5 y v φ 式中：是轮胎上方扭矩连杆和转向杆一起提供的线性弹簧扭矩；是油气式缓冲器轴承之间的粘性摩擦和减摆器的组合阻尼力矩；、分别为由侧滑和横摆引起的轮胎力矩和轮胎阻尼力矩；是磁流变阻尼器力矩。是飞机机轮侧向位移；是飞机滑行速度；是机轮摆角。 M1 M5 Mz Fy ~ 的表达式如式 (5)~(13) 所示，是围绕轮胎中心的扭转力矩；是以稳定距为杠杆作用的转弯力矩，其他变量或参数以表 1 的形式给出。表 1 起落架模型参数 Table 1 the parameters of the landing gear model 参数数值惯性矩Iz /(kg·m2 ) 1 垂直负荷Fz /N 4200 线性扭矩系数c/(N·m)·rad−1 −15000 组合阻尼系数k/(N·m·s)·rad−1 −35 稳定距e/m 0.12 轮胎力矩导数cMα/(m·rad−1) −2 侧向力导数cFα/rad−1 20 轮胎力矩极值角αg /(°) 10 轮胎力极值角δ/(°) 5 接触半长a/m 0.1 κ (N·m2 )·rad 胎面宽度常数 −1 / −126 松弛长度σ/ m 0.3 v (m·s −1 滑行速度 / ) 0~50 减摆力臂d/m 0.1 ·953· 董磊，等：采用智能算法参数整定的磁流变减摆自抗扰控制第 5 期

第17卷智能系统学报 ·954· M1=c9 (5) 长处，既有PSO算法搜索过程记录的个体历史最 M2=kφ (6) 优解及全局最优解，又加入了天牛须算法中的增 M3=M:-eFy (7) 量因子来开展最优解局部范围的搜索，改善了 M,={ 0,a叫l≥g PSO容易陷入局部最优的问题。文献[20]中将BSO cMaag/180-sin(180/ag.a), 其他 (8) F,=fCraaFs a56 算法与PSO算法、遗传算法(genetic algorithm, lcFa6F,sign(a),a≥6 (9) GA)在23组基准函数进行寻优测试，结果显示， Ma=K/v. (10) 除个别函数外，其他基准函数上，BSO算法具有 k=-0.15a2cr.F (11) 良好的搜索性能和效率。 Ms=f.d (12) 假设天牛群的种群规模为N,其中第个个体 a≈arctana=y1/o (13) 的位置更新公式为 2.2.2自抗扰控制 X+1=X+AV+(1-) (17) 针对上述磁流变阻尼器与飞机前起落架模型，式中：k为当前时刻的迭代次数；入为比例系数（取设计自抗扰控制器。因为模型的阶次为二阶，因一定值)；V代表第只天牛当前时刻的速度，速度此设计三阶的自抗扰控制器，其参考输入为0，将更新的公式为自抗扰简化为扩张状态观测器(extended state ob- V+1=wV:+cir(pbesui-X)+c2r2(gbest -X)(18) server,,ESO)、非线性反馈控制率(nonlinear state er-. 式中：为权重因子，其更新策略为线性递减，如 ror feedback,NLSEF)、扰动补偿三部分。式(19)所示；c1、c2为学习因子；n、n2为0~1之间的三阶ESO模型：随机数；Pe为第i只天牛的个体最优位置；em为 (e=z-y 此时全局最优位置。 i=z2-Bore 2=3-B2fale,0.5,)+bow (14) k =@max -(@max -Wmin)/K.k (19) 3=-Bos fal(e,0.25,6) 式中：oaax、ωmn分别表示权重的上下界；K为总送 NLSEF及扰动补偿：代次数。 e1=r-z1,e2=r-z2 除此之外，结为增量因子，其更新公式如式 uo=B1fal(e1.0.5.61)+B2fal(e2.0.25,62) u=4-3 (15) (20)所示： bo +=6V·sign(fX)-fX) (20) fal(·)为一非线性函数：式中：表示第k次的迭代步长，更新公式如下： e fal=-a,lel≤6 6l=eta·d (21) (16) el"sign(e).lel>6 式中：eta为一常数，一般取值0.95；sign为符号函将磁流变阻尼器的位移作为反馈量，产生的数；f为适应度函数，X、X分别为天牛的左须控制信号为电流，额定电流为1.2A,被控对象的和右须位置，其更新公式为观测对象为机轮摆角、机轮侧滑角、侧向位移。 X=X+·/2 (22) 可以看出，上述控制器中有B、B2、Bs、B1、B、b、 X1=X-·d'/2 6、61、62等多个参数，比较难整定，因此采用智能式中：*表示天牛左右两须之间的距离，更新公式为算法进行正定参数。 d*=6/G (23) 式中：c为一常数，一般取值为2。 3混沌分数阶天牛群算法 3.2混沌分数阶天牛群算法 3.2.1混沌化混沌分数阶天牛群算法(chaos fractional order 引入混沌来初始化种群。初始种群的分布， bettle swarm optimization,CFBSO)是将混沌思想、会影响算法的收敛速度。在天牛群算法中初始分数阶微分结构以及天牛群算法有机结合。种群一般是随机生成，种群的分布会出现不均匀 3.1天牛群算法的情况，这会减少多样性，限制搜索性能。因此，天牛群算法(bettle swarm optimization,BSO) 采用Tent混沌映射来进行种群的初始化，提高种是天牛须算法与粒子群算法(particle swarm optim- 群的多样性，如式(24)所示： ization,.PSO)的融合和改进2o。PS0中的个体用天牛来取代粒子，原来的每个粒子是一只天牛。 2-7-282021 (24) 在计算的迭代过程中位置更新，它结合了两者的而后通过式(25)映射到解空间：

M1 = cφ (5) M2 = kφ˙ (6) M3 = Mz −eFy (7) Mz = { 0, α| ⩾ αg cMααg/180 ·sin(180/αg ·α), 其他 (8) Fy = { cFααFz , α ⩽ δ cFαδFzsign(α), α ⩾ δ (9) M4 = κ/v ·φ˙ (10) κ = −0.15a 2 cFαFz (11) M5 = f · d (12) α ≈ arctanα = y1/σ (13) 2.2.2 自抗扰控制针对上述磁流变阻尼器与飞机前起落架模型，设计自抗扰控制器。因为模型的阶次为二阶，因此设计三阶的自抗扰控制器，其参考输入为 0，将自抗扰简化为扩张状态观测器 (extended state observer, ESO)、非线性反馈控制率 (nonlinear state error feedback, NLSEF)、扰动补偿三部分。三阶 ESO 模型：    e = z1 −y z˙1 = z2 −β01e z˙2 = z3 −β02 f al(e,0.5,δ)+b0u z˙3 = −β03 f al(e,0.25,δ) (14) NLSEF 及扰动补偿：    e1 = r −z1, e2 = r˙ −z2 u0 = β1 f al(e1,0.5,δ1)+β2 f al(e2,0.25,δ2) u = u0 −z3 b0 (15) fal(·) 为一非线性函数： fal =    e δ 1−α , |e| ⩽ δ |e| α sign(e), |e| > δ (16) β01 β02 β03 β1 β2 b0 δ δ1 δ2 将磁流变阻尼器的位移作为反馈量，产生的控制信号为电流，额定电流为 1.2 A，被控对象的观测对象为机轮摆角、机轮侧滑角、侧向位移。可以看出，上述控制器中有、、、、、、、、等多个参数，比较难整定，因此采用智能算法进行正定参数。 3 混沌分数阶天牛群算法混沌分数阶天牛群算法（chaos fractional order bettle swarm optimization , CFBSO）是将混沌思想、分数阶微分结构以及天牛群算法有机结合。 3.1 天牛群算法天牛群算法 (bettle swarm optimization, BSO) 是天牛须算法与粒子群算法 (particle swarm optimization, PSO) 的融合和改进[20]。PSO 中的个体用天牛来取代粒子，原来的每个粒子是一只天牛。在计算的迭代过程中位置更新，它结合了两者的长处，既有 PSO 算法搜索过程记录的个体历史最优解及全局最优解，又加入了天牛须算法中的增量因子来开展最优解局部范围的搜索，改善了 PSO 容易陷入局部最优的问题。文献 [20] 中将 BSO 算法与 PSO 算法、遗传算法 (genetic algorithm, GA) 在 23 组基准函数进行寻优测试，结果显示，除个别函数外，其他基准函数上，BSO 算法具有良好的搜索性能和效率。假设天牛群的种群规模为 N，其中第 i 个个体的位置更新公式为 X k+1 i = X k i +λV k i +(1−λ)ξ k i (17) k λ Vi i 式中：为当前时刻的迭代次数；为比例系数（取一定值）；代表第只天牛当前时刻的速度，速度更新的公式为 V k+1 i = ωkV k i +c1r1(pbest,i − X k i )+c2r2(gbest − X k i ) (18) ωk c1 c2 r1 r2 pbest,i i gbest 式中：为权重因子，其更新策略为线性递减，如式 (19) 所示；、为学习因子；、为 0~1 之间的随机数；为第只天牛的个体最优位置；为此时全局最优位置。 ωk = ωmax −(ωmax −ωmin)/K · k (19) 式中：ωmax、ωmin分别表示权重的上下界； K 为总迭代次数。 ξ k 除此之外， i 为增量因子，其更新公式如式 (20) 所示： ξ k+1 i = δ kV k i ·sign(f(X k ir)− f(X k il)) (20) δ k 式中：表示第 k 次的迭代步长，更新公式如下： δ k+1 = eta · δ k (21) eta sign f(·) X k il X k ir i 式中：为一常数，一般取值 0.95；为符号函数；为适应度函数，、分别为天牛的左须和右须位置，其更新公式为 { X k+1 il = X k i +V k i · d k /2 X k+1 ir = X k i −V k i · d k /2 (22) d 式中： k表示天牛左右两须之间的距离，更新公式为 d k = δ k /c (23) 式中：c为一常数，一般取值为 2。 3.2 混沌分数阶天牛群算法 3.2.1 混沌化引入混沌来初始化种群。初始种群的分布，会影响算法的收敛速度[21]。在天牛群算法中初始种群一般是随机生成，种群的分布会出现不均匀的情况，这会减少多样性，限制搜索性能。因此，采用 Tent 混沌映射来进行种群的初始化，提高种群的多样性，如式 (24) 所示： Zk+1 = { 2Zk , 0 ⩽ Zk ⩽ 0.5 2(1−Zk), 0.5 < Zk < 1 (24) 而后通过式 (25) 映射到解空间：第 17 卷智能系统学报 ·954·

·955· 董磊，等：采用智能算法参数整定的磁流变减摆自抗扰控制第5期 Xk=Zk(U-UD)+U (25) 其适应度值；式中U、U分别表示参数的上界和下界。 4)通过式(20)计算每个个体的增量因子，而若种群大小为N,需要优化的参数为m个，则后采用式(18)计算其速度；混沌初始化的过程为：首先产生m个0~1的随机 5)采用式(17更新每只天牛的位置：数，将之作为一个m维的向量；将此向量映射到 6)计算每只天牛的适应度值，更新个体最优解空间，记录其位置；通过式(24)产生新的m维解以及全局最优解；向量，再次映射，一共迭代Nm次；而后产生Nm(Wm>m 7)在全局最优解附近进行K次混沌搜索优化个m维的混沌序列，将混沌序列通过式(25)映射最优解；到解空间，从中选取适应度最好的N个作为初始 8)采用式(21)更新迭代步长：的天牛种群。 9)判断是否终止迭代，是则输出最优解，否在搜索最优解过程中，在全局最优解g✉的邻则返回第2)步。域内进行K次混沌搜索。过程为与上文方法一样算法流程如图5所示。得到K。个m维混沌序列，而后用式(26)得到K个开始计算每只天牛适应位于g邻域内的混沌搜索点，并将其适应度与度，更新个体最优全局最优适应度值进行比较来更新或保留全局最解与全局最优解初始化参数，并优解：初始化种群 Xx gbest RZk (26) 在全局最优解附近更新分数阶次混沌搜索以优化式中R为一常数，一般取值为0.3。更新两须距离 3.2.2分数阶更新迭代步长将分数阶引人到B$O中，可以使算法迭代过计算两须位置：程中天牛的运动具有分数阶微分的记忆特性四。计算适应度 D[xt】= 5(-1)T(a+1)x(t-hT) 更新增量因子最大迭代 Th+1)T(a-h+1) (27) 计算速度次数式(27)为α阶G-L分数阶导数的离散形式， T为采样周期，r为截断阶次。更新天牛位置结束式(19)中w取值为1，则图5混沌分数阶天牛群算法流程 Vet1-Vt cir(Pbeu-X)+car2(gbest -X) (28) Fig.5 Flowchart of CFBSO T=1时，可得 D[V]= (-1Y'T(a+1)x(t-h) 4仿真及结果与分析台rh+1ra-h+) (29) 4.1基准函数测试 cir(Poesti-X)+c2r2(gbes-X) 当r=4时，速度更新公式可以表示为将所提混沌分数阶天牛群算法与天牛群算法在23组基准函数20上进行性能测试，两种算法 V=cir(pbou-X)+czr2(ghea-X)+av+ 迭代次数均为1000次，测试结果如表2所示。其 1-awy1'+a1-u2-ay-+ (30) 中ave是算法运行30次的得到的最优值，std为 24a1-a2-a)3-)y 1 标准差。表2基准函数在不同算法下的优化结果比较式中：a采用式(31)更新： Table 2 Comparison of optimization results of benchmark a=0.8-0.5·k/K (31) function under different algorithms 式中K为迭代总次数。 CFBSO BSO 3.2.3算法流程函数 ave std ave std 加入混沌化与分数阶后，混沌分数阶天牛群 F1 11.762 15.173 48.287 60.287 算法流程如下： F2 3.3748 2.4714 6.1922 4.6076 1)初始化参数，混沌化初始种群： F3 168.47 225.87 287.26 637.92 2)如式(30)更新分数阶的阶次a,式(23)更 F4 0.4678 0.47890 3.9291 6.7811 新两须之间的距离； F5 84.525 30.012 516.25 3)用式(22)得到天牛两须的位置，然后计算 518.43

Xk = ZK(Ub −Ul)+Ul (25) 式中 Ub、Ul分别表示参数的上界和下界。 Nm Nm Nm 若种群大小为 N，需要优化的参数为 m 个，则混沌初始化的过程为：首先产生 m 个 0~1 的随机数，将之作为一个 m 维的向量；将此向量映射到解空间，记录其位置；通过式 (24) 产生新的 m 维向量，再次映射，一共迭代次；而后产生 ( >m) 个 m 维的混沌序列，将混沌序列通过式 (25) 映射到解空间，从中选取适应度最好的 N 个作为初始的天牛种群。 gbest K0 K0 K0 gbest 在搜索最优解过程中，在全局最优解的邻域内进行次混沌搜索。过程为与上文方法一样得到个 m 维混沌序列，而后用式 (26) 得到个位于邻域内的混沌搜索点，并将其适应度与全局最优适应度值进行比较来更新或保留全局最优解： Xk = gbest +RZk (26) 式中 R 为一常数，一般取值为 0.3。 3.2.2 分数阶将分数阶引入到 BSO 中，可以使算法迭代过程中天牛的运动具有分数阶微分的记忆特性[22]。 D α [x(t)] = 1 T α   ∑r h=0 (−1)h Γ(α+1)x(t−hT) Γ(h+1)Γ(α−h+1)   (27) α T r 式 (27) 为阶 G-L 分数阶导数的离散形式，为采样周期，为截断阶次。式 (19) 中ωk取值为 1，则 V k+1 i −V k i = c1r1(pbest,i − X k i )+c2r2(gbest − X k i ) (28) T = 1 时，可得 D α [V k+1 i ] =   ∑r h=0 (−1)h Γ(α+1)x(t−h) Γ(h+1)Γ(α−h+1)   = c1r1(pbest,i − X k i )+c2r2(gbest − X k i ) (29) 当r = 4 时，速度更新公式可以表示为 V k+1 i = c1r1(pbest,i − X k i )+c2r2(gbest − X k i )+αV k i + 1 2 α(1−α)V k−1 i + 1 6 α(1−α)(2−α)V k−2 i + 1 24 α(1−α)(2−α)(3−α)V k−3 i (30) 式中：α采用式 (31) 更新： α = 0.8−0.5 · k/K (31) 式中 K 为迭代总次数。 3.2.3 算法流程加入混沌化与分数阶后，混沌分数阶天牛群算法流程如下： 1）初始化参数，混沌化初始种群； 2）如式 (30) 更新分数阶的阶次α，式 (23) 更新两须之间的距离； 3）用式 (22) 得到天牛两须的位置，然后计算其适应度值； 4）通过式 (20) 计算每个个体的增量因子，而后采用式 (18) 计算其速度； 5）采用式 (17) 更新每只天牛的位置； 6）计算每只天牛的适应度值，更新个体最优解以及全局最优解； 7）在全局最优解附近进行 K0次混沌搜索优化最优解； 8）采用式 (21) 更新迭代步长； 9）判断是否终止迭代，是则输出最优解，否则返回第 2) 步。算法流程如图 5 所示。开始初始化参数, 并初始化种群更新分数阶次, 更新两须距离计算两须位置; 计算适应度更新增量因子, 计算速度更新天牛位置结束最大迭代次数 Y N 更新迭代步长在全局最优解附近混沌搜索以优化计算每只天牛适应度, 更新个体最优解与全局最优解图 5 混沌分数阶天牛群算法流程 Fig. 5 Flowchart of CFBSO 4 仿真及结果与分析 4.1 基准函数测试将所提混沌分数阶天牛群算法与天牛群算法在 23 组基准函数[20] 上进行性能测试，两种算法迭代次数均为 1 000 次，测试结果如表 2 所示。其中 ave 是算法运行 30 次的得到的最优值，std 为标准差。表 2 基准函数在不同算法下的优化结果比较 Table 2 Comparison of optimization results of benchmark function under different algorithms 函数 CFBSO BSO ave std ave std F1 11.762 15.173 48.287 60.287 F2 3.374 8 2.4714 6.1922 4.607 6 F3 168.47 225.87 287.26 637.92 F4 0.467 8 0.47890 3.9291 6.781 1 F5 84.525 30.012 516.25 518.43 ·955· 董磊，等：采用智能算法参数整定的磁流变减摆自抗扰控制第 5 期

第17卷智能系统学报 ·956· 续表2 开始 CFBSO BSO 函数混沌初始化种群 ave std 赋值给ADRC参数 ave std F6 13.992 10.309 54.013 50.907 更新分数阶次运行减摆控制系统 F7 0.0056 0.0037 0.0044 0.0029 F8 -17934 3356.1 -12143 1097.8 更新两须距离输出性能指标 F9 52.816 26.916 23.600 15.471 更新两须位置计算适应度值 F10 0.34905 0.6387 0.6679 3.6334 F11 1.1906 0.1311 8.6136 17.232 更新个体最优解更新增量因子和全局最优解 F12 0.40605 0.2365 1.6958 1.1177 计算速度 F13 3.0101 0.8664 10.673 11.029 全局最优解附近混沌搜索 F14 0.998 0 0.998 0 更新位置 F15 0.0008 0.0004 0.0032 0.0069 更新迭代步长 F16 -1.0316 0 -1.0316 0 F17 0.398 0 0.398 0 最大迭代 F18 0 0 次数 F19 -3.3349 0.2811 -3.2191 0.6438 y F20 -1.2481 0.4047 -1.1193 0.6007 输出最优解 F21 -7.2404 0.3874 -0.3138 0.0676 F22 -6.8624 0.6383 -0.5227 0.8701 结束 F23 -7.6496 0.4815 -0.5798 0.9317 图6控制流程 F23 -7.6496 0.4815 -0.5798 0.9317 Fig.6 Flowchart of control process 4.3仿真结果及分析如表2所示，除函数7，CFBSO结果比BSO 将BSO算法和CFBSO算法都应用与控制器略差，及函数16、17、18结果一样，其余函数寻优结参数寻优。各运行5次，取适应度平均值，得到其果都优于BSO算法。所以所提算法，在处理单峰与迭代次数的关系如图7所示。在初始适应度及多峰函数问题上，比B$O算法具有显著优势。值上CFBSO小于BSO,说明混沌初始化可以改 4.2基于CFBSO的自抗扰控制善种群分布，达到最后的适应度收敛值，CFBSO 上文描述了自抗扰控制器和混沌分数阶天牛所用次数更少，且函数值更小，体现了算法的优群算法，现将两者结合，用该算法优化控制器参越性。数。控制器中的6取值为0.001,6、62都取值为0.01。 0.025 其余控制器参数用算法进行整定，它们的取值范 -CFBSO 围为B1∈0,1000]、Bm∈[0,5000]、B3∈[0,20000]、 0.020 ---BSO B1∈[0,50]B2∈[0,10]、bo∈[0,50]。设置适应度函数： J=(ωlel+wlel+w3lexI+-ωald-t0d 0.010 其中w,(i=1,2,3,4)为权重因子；e为反馈误差；ea、 0.005 e为扩张状态观测器所观测的y及的误差；为控 0 20 40 60 制量电流。权重因子的取值依次为0.59、0.2、迭代次数 0.2、0.01,因电流的量级与前三者的量级相差较图7进化曲线图大，所以控制信号电流的权重设置较小。 Fig.7 The diagram of evolution curve 加入CFBSO算法后，整个控制器的流程如取两种算法优化后取得的最优控制效果下的图6所示。初始参数设置y=0m,9=0.1rad。参数如表3所示

续表 2 函数 CFBSO BSO ave std ave std F6 13.992 10.309 54.013 50.907 F7 0.005 6 0.0037 0.004 4 0.0029 F8 −17 934 3356.1 −12 143 1097.8 F9 52.816 26.916 23.600 15.471 F10 0.349 05 0.6387 0.667 9 3.6334 F11 1.190 6 0.1311 8.613 6 17.232 F12 0.406 05 0.2365 1.695 8 1.1177 F13 3.010 1 0.8664 10.673 11.029 F14 0.998 0 0.998 0 F15 0.000 8 0.0004 0.003 2 0.0069 F16 −1.0316 0 −1.0316 0 F17 0.398 0 0.398 0 F18 3 0 3 0 F19 −3.3349 0.2811 −3.2191 0.6438 F20 −1.2481 0.4047 −1.1193 0.6007 F21 −7.2404 0.3874 −0.3138 0.0676 F22 −6.8624 0.6383 −0.5227 0.8701 F23 −7.6496 0.4815 −0.5798 0.9317 F23 −7.6496 0.4815 −0.5798 0.9317 如表 2 所示，除函数 7，CFBSO 结果比 BSO 略差，及函数 16、17、18 结果一样，其余函数寻优结果都优于 BSO 算法。所以所提算法，在处理单峰及多峰函数问题上，比 BSO 算法具有显著优势。 4.2 基于 CFBSO 的自抗扰控制 δ δ1 δ2 β01 ∈ [0,1 000] β02 ∈ [0,5 000] β03 ∈ [0,20 000] β1 ∈ [0,50] β2 ∈ [0,10] b0 ∈ [0,50] 上文描述了自抗扰控制器和混沌分数阶天牛群算法，现将两者结合，用该算法优化控制器参数。控制器中的取值为 0.001，、都取值为 0.01。其余控制器参数用算法进行整定，它们的取值范围为、、、、、。设置适应度函数： J = w ((ω1|e1|+ω2|ez1 |+ω3|ez2 |+ω4|u|)·t)dt ωi(i = 1,2,3,4) e1 ez1 ez2 y y˙ u 其中为权重因子；为反馈误差；、为扩张状态观测器所观测的及的误差；为控制量电流。权重因子的取值依次为 0.59、0.2、 0.2、0.01，因电流的量级与前三者的量级相差较大，所以控制信号电流的权重设置较小。 y = 0 φ = 0.1 加入 CFBSO 算法后，整个控制器的流程如图 6 所示。初始参数设置 m, rad。开始混沌初始化种群更新分数阶次更新两须距离更新两须位置更新增量因子计算速度更新位置输出最优解更新迭代步长全局最优解附近混沌搜索更新个体最优解和全局最优解计算适应度值输出性能指标运行减摆控制系统赋值给 ADRC 参数结束最大迭代次数 Y N 图 6 控制流程 Fig. 6 Flowchart of control process 4.3 仿真结果及分析将 BSO 算法和 CFBSO 算法都应用与控制器参数寻优。各运行 5 次，取适应度平均值，得到其与迭代次数的关系如图 7 所示。在初始适应度值上 CFBSO 小于 BSO，说明混沌初始化可以改善种群分布，达到最后的适应度收敛值，CFBSO 所用次数更少，且函数值更小，体现了算法的优越性。 0 20 40 60 迭代次数 0.005 0.010 0.015 0.020 0.025 适应度值 CFBSO BSO 图 7 进化曲线图 Fig. 7 The diagram of evolution curve 取两种算法优化后取得的最优控制效果下的参数如表 3 所示。第 17 卷智能系统学报 ·956·

表 3 优化后控制器参数 Table 3 The parameters of the controller 算法模块参数取值 CFBSO ESO β01 = 61.88, β02 = 4 999.82 β03 = 16 200.11,δ = 0.001 NLESF β1 = 45.65, β2 = 6.37,b0 = 35.94 δ1 = 0.01,δ2 = 0.01 BSO ESO β01 = 60.01, β02 = 4 998.66 β03 = 13 015.38,δ = 0.001 NLESF β1 = 0.15, β2 = 9.99,b0 = 42.53 δ1 = 0.01,δ2 = 0.01 因摆角、侧向位移以及滑移角变化趋势相同，本节仿真结果只给出摆角的变化曲线。首先，在没有加入磁流变阻尼器的情况下，摆角变化曲线如图 8 所示。 0 0.1 0.2 0.3 0.4 0.5 t/s (a) v=10 m/s −0.05 0 0.05 0.10 ψ/rad 0 0.1 0.2 0.3 0.4 0.5 t/s (b) v=30 m/s ψ/rad −0.2 −0.1 0 0.1 0.2 0.3 0 0.1 0.2 0.3 0.4 0.5 t/s (c) v=50 m/s ψ/rad −0.2 −0.1 0 0.1 0.2 图 8 无 MR，摆角变化曲线 Fig. 8 No MR, the response of yaw angle 从图 8 和图 9 中可以看到在该较小的速度下，两种方式下系统都可以达到稳态，但是加入磁流变阻尼器后稳定时间和振荡周期都明显的减少，两种算法下的控制效果几乎一致，稳定时间减至 0.13 s，说明加入磁流变阻尼器可以明显抑制摆振情况，然后将速度增大得到图 8(b)(c) 的仿真结果图。 0 0.1 0.2 0.3 0.4 0.5 t/s (a) v=10 m/s −0.05 0 0.05 0.10 ψ/rad CFBSO-ADRC BSO-ADRC 0 0.1 0.2 0.3 0.4 0.5 t/s (c) v=50 m/s −0.05 −0.10 0 0.05 0.10 ψ/rad CFBSO-ADRC BSO-ADRC 0 0.1 0.2 0.3 0.4 0.5 t/s (b) v=30 m/s −0.05 −0.10 0 0.05 0.10 ψ/rad CFBSO-ADRC BSO-ADRC 图 9 有 MR，摆角变化曲线 Fig. 9 With MR, the response of yaw angle 从图 8 中可以看出，当在一定范围内，速度增大，不加入磁流变阻尼器的情况下，前起落架的机轮的摆动会呈发散状，后会等幅振荡，造成系统的不稳定，在磁流变阻尼器的作用，摆动呈现逐渐衰减的趋势，并且可以在较短的时间内达到稳定状态。而且在速度分别为 30 m/s 和 50 m/s 时，从图 9 中可以看到，在后几个振荡周期内，CFBSO-AD- ·957· 董磊，等：采用智能算法参数整定的磁流变减摆自抗扰控制第 5 期

第17卷智能系统学报 ·958· RC的幅值和稳定时间更小。采用时间乘绝对误数的整定，与天牛群算法进行了对比，通过仿真差积分(integral of time and absolute error,,ITAE)指结果可知，控制器取得了良好的控制效果，并且标，如式(32)所示。表4为两种算法下减摆性能具有较好的响应时间与鲁棒性，体现了CFBSO算指标。法下减摆自抗扰控制的良好性能。 ITAE=le(t)ltdt (32) 在未来工作中，可以就文中所提算法进行改进，比如说将反向学习和搜索的策略引入考虑算表4减摆性能指标 Table 4 Performance index for shimmy reductiom 法中超参数的整定等。同时可以找寻更合适的智能优化算法来整定自抗扰控制器的参数。控制器速度/(m·s)稳定时间/s ITAE 10 0.130 1.161×10 参考文献： CFBSO-ADRC 30 0.278 3.703×10 50 0.333 4.365×10 [1]聂宏，魏小辉.飞机起落架动力学设计与分析M.西安：西北工业大学出版社，2013. 10 0.130 1.165x104 [2] 陈大伟，顾宏斌，吴东苏.基于磁流变阻尼器的起落架 BSO-ADRC 30 0.300 3.944×10 摆振半主动控制[J].中国机械工程，2010,21(12)： 50 0.353 4.706×10 1401-1405 如表4所示，所提算法在稳定时间和ITAE指 CHEN Dawei,GU Hongbin,WU Dongsu.Semi-active 标上都有更好的控制效果。 control of landing gear shimmy based on magneto-rheolo- 4.4鲁棒性分析 gical (MR)damper[J].China mechanical engineering. 上述仿真结果图说明当速度在一定范围变化 2010.21(12):1401-1405. 时，系统同样可以保持稳定。下面改变系统的一 [3] 武柏安，龙海洋，李耀刚，等.基于遗传算法的磁流变半些参数，观察控制系统的鲁棒性。表5为参数变主动悬架最优控制[J].机床与液压，2021,49(9) 化后CFBSO-ADRC控制器性能指标的变化。 109-114. WU Boan,LONG Haiyang,LI Yaogang,et al.Optimal 表5鲁棒性性能指标 Table 5 Performance index for robustness control of MR semi-active suspension based on genetic algorithm[J].Machine tool hydraulics,2021,49(9): 速度/ms) 变化参数变化范围% 稳定时间/s 5 109-114 0.143 F [4] ZHANG Hailong,WANG Enrong,ZHANG Ning,et al. +10 0.147 10 Semi-active sliding mode control of vehicle suspension +25 0.133 with magneto-rheological damper[J].Chinese journal of +50 0.170 mechanical engineering,2015,28(1):63-75. +5 0.351 [5] DYKE S J.SPENCER B F Jr.SAIN MK.et al.Model- +10 0.440 30 ing and control of magnetorheological dampers for seis- +25 0.381 mic response reduction[J].Smart materials and structures, +50 0.527 1996,5(5):565-575. +5 0.413 [6] NGUYEN X B.KOMATSUZAKI T.IWATA Y,et al +10 0.709 Modeling and semi-active fuzzy control of magnetorhe- 50 +25 0.486 ological elastomer-based isolator for seismic response re- +50 0.566 duction[J].Mechanical systems and signal processing, 2018.101:449-466 从表5中可以看到垂直负荷和摆角初始值在 [7] 王唯，夏品奇.采用磁流变阻尼器的直升机“地面共振定范围内增加时，3种速度的情况下系统仍能分析[.南京航空航天大学学报，2003,35(3)：264-267 保持稳定。结果显示，系统对摆角初始值变化的 WANG Wei,XIA Pingi.Ground resonance analysis of 容忍度更高一些。 helicopter with magnetorheological lag damper[J].Journ- 5结束语 al of Nanjing university of aeronautics astronautics. 2003,35(3):264-267 本文对飞机前起落架设计减摆自抗扰控制 [8] KANG B H,YOON J Y,KIM G W,et al.Landing effi- 器，并通过混沌分数阶天牛群算法进行控制器参 ciency control of a six degrees of freedom aircraft model

RC 的幅值和稳定时间更小。采用时间乘绝对误差积分（integral of time and absolute error, ITAE）指标，如式 (32) 所示。表 4 为两种算法下减摆性能指标。 ITAE = w |e(t)|tdt (32) 表 4 减摆性能指标 Table 4 Performance index for shimmy reductiom 控制器速度/(m·s−1) 稳定时间/s ITAE CFBSO-ADRC 10 0.130 1.161×10−4 30 0.278 3.703×10−4 50 0.333 4.365×10−4 BSO-ADRC 10 0.130 1.165×10−4 30 0.300 3.944×10−4 50 0.353 4.706×10−4 如表 4 所示，所提算法在稳定时间和 ITAE 指标上都有更好的控制效果。 4.4 鲁棒性分析上述仿真结果图说明当速度在一定范围变化时，系统同样可以保持稳定。下面改变系统的一些参数，观察控制系统的鲁棒性。表 5 为参数变化后 CFBSO-ADRC 控制器性能指标的变化。表 5 鲁棒性性能指标 Table 5 Performance index for robustness 速度/(m·s−1) 变化参数变化范围/% 稳定时间/s 10 Fz +5 0.143 +10 0.147 φ +25 0.133 +50 0.170 30 Fz +5 0.351 +10 0.440 φ +25 0.381 +50 0.527 50 Fz +5 0.413 +10 0.709 φ +25 0.486 +50 0.566 从表 5 中可以看到垂直负荷和摆角初始值在一定范围内增加时，3 种速度的情况下系统仍能保持稳定。结果显示，系统对摆角初始值变化的容忍度更高一些。 5 结束语本文对飞机前起落架设计减摆自抗扰控制器，并通过混沌分数阶天牛群算法进行控制器参数的整定，与天牛群算法进行了对比，通过仿真结果可知，控制器取得了良好的控制效果，并且具有较好的响应时间与鲁棒性，体现了 CFBSO 算法下减摆自抗扰控制的良好性能。在未来工作中，可以就文中所提算法进行改进，比如说将反向学习和搜索的策略引入考虑算法中超参数的整定等。同时可以找寻更合适的智能优化算法来整定自抗扰控制器的参数。参考文献：聂宏, 魏小辉. 飞机起落架动力学设计与分析 [M]. 西安: 西北工业大学出版社, 2013. [1] 陈大伟, 顾宏斌, 吴东苏. 基于磁流变阻尼器的起落架摆振半主动控制 [J]. 中国机械工程, 2010, 21(12): 1401−1405. CHEN Dawei, GU Hongbin, WU Dongsu. Semi-active control of landing gear shimmy based on magneto-rheological (MR) damper[J]. China mechanical engineering, 2010, 21(12): 1401−1405. [2] 武柏安, 龙海洋, 李耀刚, 等. 基于遗传算法的磁流变半主动悬架最优控制 [J]. 机床与液压, 2021, 49(9): 109−114. WU Boan, LONG Haiyang, LI Yaogang, et al. Optimal control of MR semi-active suspension based on genetic algorithm[J]. Machine tool & hydraulics, 2021, 49(9): 109−114. [3] ZHANG Hailong, WANG Enrong, ZHANG Ning, et al. Semi-active sliding mode control of vehicle suspension with magneto-rheological damper[J]. Chinese journal of mechanical engineering, 2015, 28(1): 63–75. [4] DYKE S J, SPENCER B F Jr, SAIN M K, et al. Modeling and control of magnetorheological dampers for seismic response reduction[J]. Smart materials and structures, 1996, 5(5): 565–575. [5] NGUYEN X B, KOMATSUZAKI T, IWATA Y, et al. Modeling and semi-active fuzzy control of magnetorheological elastomer-based isolator for seismic response reduction[J]. Mechanical systems and signal processing, 2018, 101: 449–466. [6] 王唯, 夏品奇. 采用磁流变阻尼器的直升机“地面共振” 分析 [J]. 南京航空航天大学学报, 2003, 35(3): 264−267. WANG Wei, XIA Pinqi. Ground resonance analysis of helicopter with magnetorheological lag damper[J]. Journal of Nanjing university of aeronautics & astronautics, 2003, 35(3): 264−267. [7] KANG B H, YOON J Y, KIM G W, et al. Landing efficiency control of a six degrees of freedom aircraft model [8] 第 17 卷智能系统学报 ·958·

·959· 董磊，等：采用智能算法参数整定的磁流变减摆自抗扰控制第5期 with magneto-rheological dampers:part 2-control simu- of-art research on dynamical,models of magnetorheolo- lation[J].Joural of intelligent material systems and struc- gical damper[J].Journal of vibration and shock,2006, tures,2021,32(12):1303-1315. 25(3:121-126.211. [9]FU Li.GUAN Wei.WANG Qi,et al.Modeling and ad- [18]SAVAIA G.PANZANI G,CORNO M,et al.Hammer- aptive control of magneto-rheological buffer system for stein-Wiener modelling of a magneto-rheological aircraft landing gear[J].Applied mathematical modelling, dampers considering the magnetization dynamics[J].Con- 2015,399y2509-2517. trol engineering practice,2021,112:104829. [10]李莹.飞机起落架飚流变减摆器控制系统的设计与实 [19]SOMIESKI G.Shimmy analysis of a simple aircraft nose 现D]天津：中国民航大学，2016：19-26， landing gear model using different mathematical meth- LI Ying.Design and implementation of the aircraft land- ods[J].Aerospace science and technology,1997,1(8): ing gear of MR shimmy damper control system[D]. 545-555. Tianjin:Civil Aviation University of China,2016:19-26. [20]WANG Tiantian,YANG Long,LIU Qiang.Beetle swarm [11]王恕浩.基于磁流变诚摆器的飞机前轮减摆控制研究 optimization algorithm:theory and application[J].Filo- [D1.天津：中国民航大学，2019：23-32 mat,2020,34(15)5121-5137 WANG Shuhao.Shimmy damping control research of [21]王坚浩，张亮，史超，等.基于混沌搜索策路的鲸鱼优化 aircraft's nose wheel based on magnetorheological 算法U.控制与决策，2019,34(9)：1893-1900 shimmy damper[D].Tianjin:Civil Aviation University of WANG Jianhao,ZHANG Liang,SHI Chao,et al.Whale China,2019:23-32. optimization algorithm based on chaotic search strategy[. [12]王璐琦.磁流变诚摆器控制器改进设计与算法研究 Control and decision,2019,34(9):1893-1900 D].天津：中国民航大学，2020：1-22 [22】]陈秋鑫.分数阶粒子群优化及其在聚类分析中的应用 WANG Luqi.Research on improved design and al- D].镇江：江苏科技大学，2019 gorithm of magnetorheological shimmy damper control- CHEN Qiuxin.Fractional-order particle swarm optimiza- ler[D].Tianjin:Civil Aviation University of China,2020: tion and its application in clustering analysis[D].Zhenji- 1-22 ang:Jiangsu University of Science and Technology,2019. [13]DONG Lei,CHEN Zengqiang,SUN Mingwei,et al. 作者简介： Study on magneto-rheological damper control of aircraft 董磊，硕士研究生，主要研究方向 landing gear based on LADRC[M]//Lecture Notes in 为智能预测控制、自抗扰控制。 Electrical Engineering.Singapore:Springer Singapore, 2021:830-840. [14]韩京清.自抗扰控制器及其应用.控制与决策，1998， 13(1):19-23 HAN Jingging.Auto-disturbances-rejection controller and its applications[J].Control and decision,1998,13(1): 陈增强，教授，博士生导师，中国自动化学会数据驱动控制、学习与优 19-23. 化专业委员会副主任，中国人工智能 [15]李杰，齐晓慧，万慧，等.自抗扰控制：研究成果总结与学会智能空天专业委员会副主任，天展望).控制理论与应用，2017,34(3)：281-295. 津市自动化学会理事。主要研究方向 LI Jie,QI Xiaohui,WAN Hui,et al.Active disturbance 为智能控制、预测控制、自抗扰控制。 rejection control:theoretical results summary and future 主持完成国家863项目和国家自然科 researches[J].Control theory applications,2017,34(3): 学基金项目6项，获得省部级科技奖5次，发表学术论文 300余篇。 281-295. [16]AHAMED R,FERDAUS MM,LI Yancheng.Advance- 孙明玮，教授，博士生导师，中国 ment in energy harvesting magneto-rheological fluid 仿真学会理事，中国自动化学会数据驱动控制、学习与优化专业委员会委 damper:a review[J].Korea-Australia rheology journal, 员，主要研究方向为飞行器制导与控 2016,28(4):355-379. 制、自抗扰控制。主持国防科技攻关 [17刀邓志党，高峰，刘献栋，等.磁流变阻尼器力学模型的研项目和国家自然科学基金项目4项，究现状[).振动与冲击，2006,25(3)：121-126,211. 获得国防科技进步奖3次，发表学术 DENG Zhidang,GAO Feng,LIU Xiandong,et al.State- 论文50余篇

with magneto-rheological dampers: part 2—control simulation[J]. Journal of intelligent material systems and structures, 2021, 32(12): 1303–1315. FU Li, GUAN Wei, WANG Qi, et al. Modeling and adaptive control of magneto-rheological buffer system for aircraft landing gear[J]. Applied mathematical modelling, 2015, 39(9): 2509–2517. [9] 李莹. 飞机起落架磁流变减摆器控制系统的设计与实现 [D]. 天津: 中国民航大学, 2016: 19−26. LI Ying. Design and implementation of the aircraft landing gear of MR shimmy damper control system[D]. Tianjin: Civil Aviation University of China, 2016: 19−26. [10] 王恕浩. 基于磁流变减摆器的飞机前轮减摆控制研究 [D]. 天津: 中国民航大学, 2019: 23−32. WANG Shuhao. Shimmy damping control research of aircraft's nose wheel based on magnetorheological shimmy damper[D]. Tianjin: Civil Aviation University of China, 2019: 23−32. [11] 王璐琦. 磁流变减摆器控制器改进设计与算法研究 [D]. 天津: 中国民航大学, 2020: 1−22. WANG Luqi. Research on improved design and algorithm of magnetorheological shimmy damper controller[D]. Tianjin: Civil Aviation University of China, 2020: 1−22. [12] DONG Lei, CHEN Zengqiang, SUN Mingwei, et al. Study on magneto-rheological damper control of aircraft landing gear based on LADRC[M]//Lecture Notes in Electrical Engineering. Singapore: Springer Singapore, 2021: 830−840. [13] 韩京清. 自抗扰控制器及其应用 [J]. 控制与决策, 1998, 13(1): 19−23. HAN Jingqing. Auto-disturbances-rejection controller and its applications[J]. Control and decision, 1998, 13(1): 19−23. [14] 李杰, 齐晓慧, 万慧, 等. 自抗扰控制: 研究成果总结与展望 [J]. 控制理论与应用, 2017, 34(3): 281−295. LI Jie, QI Xiaohui, WAN Hui, et al. Active disturbance rejection control: theoretical results summary and future researches[J]. Control theory & applications, 2017, 34(3): 281−295. [15] AHAMED R, FERDAUS M M, LI Yancheng. Advancement in energy harvesting magneto-rheological fluid damper: a review[J]. Korea-Australia rheology journal, 2016, 28(4): 355–379. [16] 邓志党, 高峰, 刘献栋, 等. 磁流变阻尼器力学模型的研究现状 [J]. 振动与冲击, 2006, 25(3): 121−126, 211. DENG Zhidang, GAO Feng, LIU Xiandong, et al. State- [17] of-art research on dynamical, models of magnetorheological damper[J]. Journal of vibration and shock, 2006, 25(3): 121−126, 211. SAVAIA G, PANZANI G, CORNO M, et al. Hammerstein-Wiener modelling of a magneto-rheological dampers considering the magnetization dynamics[J]. Control engineering practice, 2021, 112: 104829. [18] SOMIESKI G. Shimmy analysis of a simple aircraft nose landing gear model using different mathematical methods[J]. Aerospace science and technology, 1997, 1(8): 545–555. [19] WANG Tiantian, YANG Long, LIU Qiang. Beetle swarm optimization algorithm: theory and application[J]. Filomat, 2020, 34(15): 5121–5137. [20] 王坚浩, 张亮, 史超, 等. 基于混沌搜索策略的鲸鱼优化算法 [J]. 控制与决策, 2019, 34(9): 1893−1900. WANG Jianhao, ZHANG Liang, SHI Chao, et al. Whale optimization algorithm based on chaotic search strategy[J]. Control and decision, 2019, 34(9): 1893−1900. [21] 陈秋鑫. 分数阶粒子群优化及其在聚类分析中的应用 [D]. 镇江: 江苏科技大学, 2019. CHEN Qiuxin. Fractional-order particle swarm optimization and its application in clustering analysis[D]. Zhenjiang: Jiangsu University of Science and Technology, 2019. [22] 作者简介：董磊，硕士研究生，主要研究方向为智能预测控制、自抗扰控制。陈增强，教授，博士生导师，中国自动化学会数据驱动控制、学习与优化专业委员会副主任，中国人工智能学会智能空天专业委员会副主任，天津市自动化学会理事。主要研究方向为智能控制、预测控制、自抗扰控制。主持完成国家 863 项目和国家自然科学基金项目 6 项，获得省部级科技奖 5 次，发表学术论文 300 余篇。孙明玮，教授，博士生导师，中国仿真学会理事，中国自动化学会数据驱动控制、学习与优化专业委员会委员，主要研究方向为飞行器制导与控制、自抗扰控制。主持国防科技攻关项目和国家自然科学基金项目 4 项，获得国防科技进步奖 3 次，发表学术论文 50 余篇。 ·959· 董磊，等：采用智能算法参数整定的磁流变减摆自抗扰控制第 5 期

点击进入文档下载页（PDF格式）

已到末页，全文结束

点击下载（PDF格式）

浏览记录