可旋转翼式弹道修正组件滚转控制技术研究

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第46卷 第4期2024年4月系统工程与电子技术
SystemsEngineeringa
ndElectronicsVol.46 No.4
Ap
ril2024文章编号:1001 506X(2024)04 1412 10 网址:www.sy
s ele.com收稿日期:20230511;修回日期:20230905;网络优先出版日期:20231027。

网络优先出版地址:http:
∥kns.cnki.net/kcms/detail/11.2422.TN.20231027.1701.012.html基金项目:上海市2022年度“科技创新行动计划”启明星项目申报指南/启明星培育(扬帆专项)(22YF1447000)
资助课题 通讯作者.
引用格式:郑秋实,许伟春,赵明翰,等.可旋转翼式弹道修正组件滚转控制技术研究[J].系统工程与电子技术,2024,46(4):1412 1421.犚犲犳犲狉犲狀犮犲犳狅狉犿犪狋:ZHENGQS,XUWC,ZHAOMH,etal.Researchonrollcontroltechnologyoftrajectoryc
orrectionfusewithactive canards[J].SystemsEngineeringa
ndElectronics,2024,46(4):1412 1421.可旋转翼式弹道修正组件滚转控制技术研究
郑秋实1, ,许伟春2,赵明翰1,李乃星1,包旭馨1
(1.上海无线电设备研究所,上海201109;2.上海航天技术研究院,上海201109)
摘 要:针对可旋转翼式弹道修正组件滚转通道控制中存在的未建模摩擦干扰、
参数不确定性和外部随机干扰造成的复合扰动问题,提出一种基于扩张状态观测器(extendedstateobserver,ESO)的滑模控制方法。

首先建立弹道修正组件滚转通道模型,将动力学模型中存在的外部干扰、未建模摩擦干扰和参数摄动整合为复合干扰,然后设计ESO对修正组件滚转通道模型中难以直接测定的状态变量以及复合干扰进行估计,并基于估计值结合滑模控
制理论设计滚转通道控制器,
实现对滚转角指令的精确跟踪。

综合考虑ESO和滑模控制器构成的闭环控制系统,利用Lyap
onov稳定性理论证明了所设计的闭环控制系统的稳定性。

最后,通过仿真实验分析,证明所设计的修正组件滚转通道控制器,对滚转角指令的瞬态响应和稳态性能优异,同时可以有效抑制系统复合扰动,具备较强的鲁棒性。

关键词:可旋转翼式弹道修正组件;滚转控制技术;扩张状态观测器;滑模控制;扰动估计和抑制技术中图分类号:TJ413 文献标志码:A 犇犗犐:10.12305/j.
issn.1001 506X.2024.04.30犚犲狊犲犪狉犮犺狅狀狉狅犾犾犮狅狀狋狉狅犾狋犲犮犺狀狅犾狅犵狔狅犳狋狉犪犼犲犮狋狅狉狔犮
狅狉狉犲犮狋犻狅狀犳狌狊犲狑犻狋犺犪犮狋犻狏犲 犮犪狀犪狉犱狊
ZHENGQiushi1, ,XUWeichun2,ZHAOMinghan1,LINaixing1,BAOXu
xin1(1.犛犺犪狀犵犺犪犻犚犪犱犻狅犈狇狌犻狆犿犲
狀狋犚犲狊犲犪狉犮犺犐狀狊狋犻狋狌狋犲,犛犺犪狀犵犺犪犻201109,犆犺犻狀犪;2.犛犺犪狀犵犺犪犻犃犮犪犱犲犿狔狅犳犛狆犪犮犲犳犾犻犵犺狋犜犲犮犺狀狅犾狅犵狔,犛犺
犪狀犵犺犪犻201109,犆犺犻狀犪) 犃犫狊狋狉犪犮狋:Aimingatthecompositedisturbanceproblemcausedbyu
nmodeledfrictiondisturbance,parameteruncertainty,a
ndexternalrandomdisturbanceintherollingchannelcontrolforatrajectorycorrectionfusewithactive canards,aslidingmodecontrolmethodbasedontheextendedstateobserver(ESO)isproposed.Firstly,themodeloftherollingchannelofthetrajectoryc
orrectionfuseisestablished.Theexternaldisturbances,unmodeledfrictionforceandparameterperturbationsareinteg
ratedintoacompositedisturbance.ThentheESOisdesignedtoestimatethestatevariablesandcompositedisturbancesthataredifficulttodirectly
measureinrollchannelmodelofcorrectionfuse,andtherollchannelcontrollerisdesignedbasedontheestimatedvaluesandtheslidingmodecontroltheorytorealizetheaccuratetrackingoftherollanglecommand.ComprehensivelyconsideringtheclosedloopcontrolsystemcomposedofESOandslidingmodecontroller,Lyapunovstabilitytheoryisusedtoprovethestabilityoftheproposedloopcontrolsystem.Finally,thesimulationresultverifiesthattherollchannelhasgreatperformanceofinstantresponsetorollangleindicator,whichcaneffectivelysurpassthecompositedisturbance,andhasstrongrobustnessandgoodstability.
犓犲狔狑狅
狉犱狊:trajectorycorrectionfusewithactive canards;rollcontrol;extendedstateobserver(ESO);slidingmodecontrol;disturbanceestimationandcompensationtechniq
ues0 引 言为满足现代战争对低成本、高精度制导弹药的需求,简易制导弹药技术应运而生。

其中,弹道修正组件技术凭借其成本低、
结构简单、精度高等优点,受到越来越多的青睐,具有广阔的发展前景[13]。

弹道修正组件上安装有制导控
 第4期
郑秋实等:可旋转翼式弹道修正组件滚转控制技术研究·1413 · 制模块和执行机构,可以采用同等封装直接替换炮弹引信,将常规无控炮弹改进为具备制导控制能力的弹道修正弹。

在飞行过程中,修正组件可以实时获取弹道信息,发出指令,
控制执行机构对飞行弹道进行修正,从而减小弹着点偏差,
提高命中精度[46]。

弹道修正组件的执行机构通常采用鸭舵机构,根据舵片能否旋转,弹道修正组件可分为固定翼式和可旋转翼式[7]。

固定翼式修正组件的结构较为简单,其鸭舵偏转角是固定不变的,无法改变控制力大小和方向,因此很难在飞行弹道内进行全程修正。

此类固定翼式修正组件通常采用电磁执行机构控制组件的滚转运动。

程杰等[8]研究了不同发射条件下修正组件电磁执行机构的工况,为后续电磁执行机构的设计提供了初期设计指标,随后通过瞬态数值计算和动态风洞试验研究了双旋弹内回路的滚转特性,并通过飞行试验进行验证[9]。

殷婷婷等[10]基于电磁执行机构力矩输出特性测试试验,研究建立了修正组件控制响应模型,随后基于得出的模型,提出了一种基于扩张状态观测器(extendedstateobserver,ESO)的转速控制器[11],通过设计输出反馈型ESO对系统整体扰动进行估计,并以前馈的方式进行补偿,以减弱模型不确定性和扰动对跟踪性能的影响。

张鑫等[12]通过对电磁执行机构电气原理的分析,建立了更为准确的修正组件滚转通道动力学模型,并设计了两自由度H∞回路整形滚转通道控制器,实现了精确的滚转角控制。

相比于固定翼式,可旋转翼式修正组件的鸭舵偏转角是连续可控的,可以借鉴常规旋转弹的控制方式,对弹药的飞行弹道进行制导控制,因此能够获得更高的精度。

Ilg[13]提出的智能弹药控制方案采用了比例式舵机控制的一对鸭舵进行单通道正弦控制,以实现二维弹道修正。

朱大林[14]在其博士论文中采用了两对可旋转鸭舵进行双通道旋转控制,实现了双旋弹道修正弹的二维弹道修正,并对其飞行稳定性进行了分析[15]。

但是两者并没有对修正组件滚转通道进行控制,而是通过自然减旋运动使组件转速保持在低速旋转状态。

为获取更高的控制精度,Theodoulis等[16]基于H∞控制理论研究设计了双旋弹倾斜转弯(bank to turn,BTT)自驾仪,并指出对修正组件滚转角的快速精确控制是极其关键的。

针对修正组件滚转控制问题,Salman等[17]采用了经典比例积分微分(proportion integ
ral derivative,PID)控制方法和H∞控制理论设计了滚转通道控制器,Theodoulis等[18]采用固定结构H1控制方法设计了155mm鸭舵控制双旋弹的前体滚转通道自驾仪,实现了对双旋弹前体滚转角的快速精确控制。

综合前人研究发现,使用可旋转翼式弹道修正组件的弹道修正弹,可以采用常规旋转弹使用的BTT控制方式,以获取更高的制导控制精度,而修正组件滚转角的快速精确控制是实现BTT控制的关键。

然而,在实际的弹道修正组件系统中,受限于安装空间和低成本要求,无法安装过多的传感器,许多系统状态变量无法通过直接测量获取,因此一些复杂度较高的控制方法并不适用。

同时,由于弹道修正弹飞行跨域大,系统具有强非线性,并且存在气动参数摄动、摩擦等未建模干扰力矩以及外部干扰等问题,这些问题大大增加了修正组件滚转控制的难度。

目前针对可旋转翼式修正组件滚转控制的研究较少,已有的几篇文献(文献[1618]
)均基于线性化鲁棒控制理论,缺乏稳定性保证,并且设计难度大、耗时长,难以在实际工程中实现。

针对此类非线性不确定系统的一种有效方法是滑模控制法[1924],其具有响应速度快、对扰动和参数变化不敏感等优点,但是同时具有抖振问题、结构复杂等缺点,并且滑模控制需要的已知系统状态量较多,因此普通滑模控制法并不适用于弹道修正组件。

韩京清等[2526]针对一类存在不确定性的系统,提出一种ESO对系统的不确定性进行估计。

ESO可以估计系统的状态以及不确定性和干扰的影响,从而
实现干扰抑制或补偿。

文献[2729]均利用ESO对系统中的扰动和变量进行估计,从而达到良好的控制精度和抗干扰能力。

Wang等[30]对ESO、高增益观测器(highgainobserver,HGO)和滑模观测器(slidingmodeobserver,SMO)这三种观测器进行了比较研究,结果表明ESO在处理不确定性、干扰和传感器噪声方面总体性能更优。

基于以上研究,本文提出一种基于ESO的滑模控制方法,首先针对可旋转翼式弹道修正组件的特点,介绍了修正组件的控制原理,接着通过对修正组件所受力矩进行分析,建立了修正组件滚转通道的动力学模型;然后将系统存在的气动参数不确定性、未建模摩擦干扰和外部干扰整合为一个总的复合干扰,设计ESO对该复合干扰以及实际工程中无法直接测量的状态变量进行估计,随后结合滑模控制理论,设计具有足够鲁棒性的滚转通道控制器,并进行了稳定性分析;最后利用Simulink仿真实验进行验证,并对仿真结果进行分析及得出结论。

本文设计的修正组件
滚转控制器只需要一个滚转角反馈信号,既能够对速度、加速度和复合干扰信号进行观测,同时又能够对系统中未建模摩擦力、随机外部干扰和参数不确定性造成的复合扰动进行抑制,保证滚转控制器具有较高的滚转角跟踪精度和鲁棒性。

1 可旋转翼式修正组件滚转控制原理可旋转翼式修正组件通常采用电动舵机作为执行机构,
控制两对鸭舵按照指令连续偏转。

如图1所示,修正组件上安装有两对鸭舵,其中舵片1和舵片3构成一对差动舵,可为修正组件提供导转力矩,舵片2和舵片4构成一对同向舵,用于产生侧向控制力和控制力矩。

修正组件与弹体之间通过轴承连接,两者之间可以相对自由旋转。

当需要产生某方向控制作用时,滚转通道控制器通过控制差动舵偏转,控制修正组件旋转至预定角度,然后同向舵偏转,产生一定方向和大小的控制力,进而控制弹体飞行方向。

 ·1
414 ·系统工程与电子技术
第46卷
图1 Fig.1 Schematicdiagramofthetrajectoryc
orrectionfusewithactive canards
2 修正组件滚转通道建模
2.1 修正组件滚转力矩分析
首先对弹道修正组件在飞行过程中所受滚转力矩进行分析,如图2所示,修正组件受到的滚转力矩包括:差动舵1~差动舵3产生的控制力矩犜犆,由组件滚转引起的气动阻尼
力矩犜狆犳,组件与弹体之间的滚转阻尼力矩犜犳。

图中狆犳为组件滚转速度,狆犪为弹体滚转速度, 犳为组件滚转角度。

设定组件与弹体旋转方向相同(从弹头看向弹尾,逆时针旋转),并设定此方向为滚转速度和滚转力矩的正方向。

图2 修正组件所受力矩分析
Fig.2 Momentsanaly
sisofthecorrectionfuse差动舵控制力矩的表达式为
犜犆=12ρ
犛犔犞2犆犕δ狓δ狓(1)气动阻尼力矩表达式为
犜狆犳=12ρ
犛犔犞(2犔)
犞犆lpf狆犳(2)式中:ρ为大气密度;犛为参考面积;犔为参考长度;犞为弹丸飞行速度;δ狓为差动舵偏角;
犆犕δ狓和犆lpf分别为差动舵控制力矩系数导数和组件气动滚转阻尼力矩系数导数,可通过计算流体动力学(computationalfluiddynamics,CFD)气动仿真或者风洞试验获得。

组件与弹体之间的滚转阻尼力矩由轴承的粘性阻尼力矩和滚转摩擦力矩组成,粘性阻尼力矩与组件和弹体相对转速有关,滚转摩擦力与轴承所承受的法向力有关,具体表达式为
犜犳=犆犞(狆犪-狆犳)+犆犚|犉犖|sgn(狆犪-狆犳)(3)式中:犆犞为粘性阻尼力矩系数;犆犚为滚转摩擦力矩系数;
犉犖为作用在轴承上的法向力;sg
n(·)为符号函数。

其中,粘性阻尼力矩系数和滚转摩擦力矩系数的理论值较为复杂,一般可通过实验方法进行估计。

综上所述,修正组件滚转通道所受的总滚转力矩为
犜狓=犜犆-犜pf+犜犳(4)
完成对组件所受力矩的分析后,即可开始进行滚转通道建模。

2.2 修正组件滚转通道建模
修正组件滚转通道模型包括如下的姿态角运动学方程和转速动力学方程:
·犳=狆犳犐狓犳狆·
犳=犜烅烄烆
狓(5)式中:犐狓犳为修正组件的轴向转动惯量。

将各力矩表达式代
入式(5),可得: ·犳=狆犳犐狓犳狆·犳=-12ρ犛犔犞(2犔)
犞犆lpf狆犳+12ρ犛犔犞2犆犕δ狓δ狓+犜烅烄烆
犳(6)
式中:滚转阻尼力矩犜犳的理论值及计算值较为复杂,并且在实际工况中,还会存在弹道初始扰动、随机风等各种扰动的影响。

为简化分析,本文将滚转阻尼力矩和各种其他干扰力矩统一视为外部干扰犱ext进行处理。

同时,考虑将差动舵动力学模型视为如下一阶环节:
δ狓(狊)δ狓犮(
狊)=1犜狊狊+1(7)式中:δ狓(狊)为差动舵输出;δ狓犮(狊)为差动舵指令;犜狊为差动舵系统时间常数。

由式(6)和式(7)可得修正组件滚转通道模型框图如图3所示。

Fig.3 Blockdiag
ramoftherollchannelmodel图3中,犫1=12犐狓犳ρ犛犔犞2犔()
犞犆lpf,犫2=12犐狓犳ρ犛犔犞2犆犕δ狓。

忽略外部干扰犱ext的影响,可得修正组件滚转通道中由差动舵指令δ狓犮(狊)到滚转角 犳(狊)
的传递函数为 犳(狊)δ狓犮(狊)=犫2(犜狊
狊+1)(狊+犫1)狊=犫狊3+犪1狊2+犪2狊+犪3(8)式中:犪1=(犜-1狊+犫1);犪2=犜-1狊犫1;犪3=0;犫=犫2犜-1狊。

由此得到的是理想被控对象的传函,而在实际系统中还存在参数不确定性和外部干扰的影响。

本文将参数不确定性体现为式(8)中参数犪犻(犻=1,2,3)和犫的摄动,定义参数犪犻=犪犻0+Δ犪犻和犫=犫0+Δ犫。

其中,犪犻0和犫0代表各参数的标称值,Δ
犪犻和Δ犫代表各参数的摄动值。

当考虑参数不确定性和外部干扰时,将式(8)代表的被控对象改写为如下的空间状态方程:
狓·=犃狓+犅狌+犇犱eq狔=烅烄烆犆狓(9)式中:狓=[
狓1狓2狓3]T为系统状态变量,分别对应修正组件滚转运动中的滚转角、滚转角速度和滚转角加速度变
 第4期
郑秋实等:可旋转翼式弹道修正组件滚转控制技术研究·1415 · 量,狔= 犳为输出滚转角,狌=δ狓犮为输入的舵偏指令,系统矩阵犃,犅,犆和扰动输入矩阵犇的表达式如下:犃=010001-犪30-犪20-犪熿燀燄燅10犅=[0 0 犫0]T犆=[1 0 0]犇=[0 0 1]熿燀T(10)犱eq代表施加于系统的复合干扰,包含了未建模摩擦干扰、外部干扰和参数不确定性,其表达式为犱 -Δ犪3狓1-Δ犪2狓2-Δ犪1狓3+Δ犫狌+ω(11)式中:ω为外部干扰。

由此,得到考虑参数不确定性和外部干扰影响的滚转通道控制模型。

下面将针对此模型设计能够理想地跟踪输入信号并且对参数不确定性和外部干扰具有足够鲁棒性的滚转控制器。

3 基于犈犛犗的滑模控制器设计本节将ESO与滑模控制理论相结合,针对修正组件滚
转通道设计出滚转角控制器。

控制系统框图如图4所示。

图4 修正组件滚转通道控制系统框图
Fig.4 Blockdiagramofrollchannelcontrolsystemforcorrectionfuse首先,设计ESO对系统中无法直接测量的状态变量和由未建模摩擦干扰、外部干扰和参数不确定性等引起的复合
干扰量犱进行估计;然后利用滑模控制理论,设计滚转角控制器,并对控制量进行复合干扰补偿;最后通过Lyaponov稳定性理论验证闭环控制系统的稳定性。

3.1 犈犛犗设计不同于传统的观测器,ESO将作用在系统中的不确定性和扰动估计为原始系统的扩张状态,由此被称为ESO。

ESO是相对独立于系统的数学模型,
结构内部固有鲁棒性,性能更好并且易于实现。

考虑式(9)表示的系统模型,将复合干扰犱设计为系统的一个扩张状态狓4,则式(9)中的系统可改写为如下的状
态空间形式:
狓·1=狓2狓·2=狓3狓·3=狓4-犪30狓1-犪20狓2-犪10狓3+犫0狌狓·4=犺
狔=狓烅烄烆
1(12)式中:犺=犱·为复合干扰的变化率,该值为未知但是有界。

针对式(12)
代表的扩张系统,建立如下的状态观测器:狓^·1=狓^2+β1(狓1-狓^1)狓^·
2=狓^3+β2(狓1-狓^1)狓^·3=狓^4+β3(狓1-狓^1)-犪30狓^1-犪20狓^2-犪10狓^3+犫0狌狓^·
4=β4(狓1-狓^1烅烄烆)(13)式中:狓^犻(犻=1,2,3,4)为ESO对系统中各状态变量的估计
值;β犻(犻=1,2,3,4)为观测器增益,取正实数。

定义狏=[狏1 狏2 狏3 狏4]T为ESO的观测误差向量,其中狏犻=狓犻-狓^犻(犻=1,2,3,4)为ESO估计误差,则可得ESO观测误差状态方程为
狏·=犃-狏+犎(14)
式中:犃-
=-β1100-β2010-β3-犪30-犪20-犪101-β4熿燀燄燅000;犎=000熿燀燄燅犺。

可以通过β犻(犻=1,2,3,4)使犃-为Hurwitz,
则存在任意给定的对称正定阵犙和犘满足Lyapunov方程:犃-T犘+犘犃-
+犙=0,即存在ε使得狏1≤ε成立。

3.2 基于犈犛犗的滑模控制器设计本节将利用第3.1节中ESO对系统状态和总干扰的估计值,结合滑模控制理论进行控制器设计。

针对式(9)表示的修正组件滚转通道系统模型,定义滚转角跟踪误差为犲=狓1-狓1犱(15)
式中:狓1犱为滚转角指令,设计滑模函数为狊=犮1犲+犮2犲·+犲¨=犮1犲+犮2(狓2-狓·1犱)+(狓3-狓¨1犱)(16)式中:犮1和犮2为常数,且满足Hurwitz稳定条件,对式(16)求导可得:狊·=犮1犲·+犮2犲¨+犲-=犮1(狓2-狓·1犱)+犮2(狓3-狓¨1犱)+(-犪30狓1-犪20狓2-犪10狓3+犫0狌+犱-狓…1犱)(17)在实际工程应用中,系统状态狓2和狓3难以直接测定,因此控制器设计使用的是第
3.1节所设计的ESO的估计值狓^2和狓^3,
同时利用ESO对滑模控制器复合干扰犱的估计值狓^4进行干扰补偿,则实际控制器滑模函数为狊^=犮1犲+犮2(狓^2-狓·1犱)+(狓^3-狓¨1犱)(18)采用指数趋近项狊·=-犽狊,由式(17)和式(18)可得基于ESO的滑模控制器:狌=-1犫0
[犮1(狓^2-狓·1犱)+犮2(狓^3-狓¨1犱)-犪30狓1-犪20狓^2-犪10狓^3+狓^4-狓…1犱+犽狊^](19)4 闭环控制系统稳定性首先定义滑模控制器的Lyapunov函数为
 ·1
416 ·系统工程与电子技术
第46卷
犞狊=12狊
2(20)对犞狊求导可得:犞·狊=狊狊·=犮1(狓2-狓·1犱)+犮2(狓3-狓¨1犱)+(-犪30狓1-犪20狓2-犪10狓3+犫0狌+犱-狓…1犱)(21)将式(19)代入式(21)并化简后可得:犞·狊=-犽狊2+狊[(犮1-犪20)狏~2+(犮2-犪10)狏~3+狏~4](
22)取|(犮1-犪20)狏~2+(犮2-犪10)狏~3+狏~4|≤Δmax,则有:犞·狊≤-犽狊2+狊Δmax≤-犽狊2+12
(狊2+Δ2max)=(-犽-)
12狊2+12Δ2max=-(2犽-1)犞狊+12Δ2max(23)引理1[31] 针对犞:[0,∞)∈犚,不等式方程犞·
≤-α犞+
犳,
狋≥狋0≥0的解为犞(狋)≤e-α(狋-狋0)犞(狋0)+∫
狋狋0e-α(狋-τ)犳(τ)dτ(
24)式中:α为任意常数。

根据引理1,令α=2犽-1,犳=12Δ2max,则式(
23)的解为犞狊(狋)≤e-α(狋-狋0)犞狊(狋0)+12Δ2max∫
狋狋0e-α(狋-τ)dτ=
e-(2犽-1)(狋-狋0)犞狊(狋0)-12(2犽-1
)Δ2max(1-e-α(狋-狋0))(25)当取犽>12时,
有:lim狋→∞犞狊(狋)≤12(2
犽-1)Δ2max(26)即狋→∞时,犞狊(狋)收敛于12(2犽-1)Δ2max,且收敛速度取
决于控制增益犽和观测器参数。

然后定义ESO的Lyap
onov函数为犞狅=狏T犘狏(27)
对式(27)求导可得:犞·狅=狏·T犘狏+狏T犘狏·=(犃-狏+犎)T犘狏+狏T犘(犃-狏+犎)=狏T犃-T犘狏+犎T犘狏+狏T犘犃-狏+狏T犘犎=狏T(犃-T犘+犘犃-
)狏+
犎T犘狏+狏T犘犎=-狏T犙狏+犎T犘狏+狏T犘犎≤
-λmin(犙)狏2+2Φ犘狏
(28)式中:λmin(犙)为犙的最小特征值,由式(28)可知ESO的收敛条件为
狏≤2Φ犘λmin(犙
)(29)综合考虑ESO和滑模控制器构成的闭环系统,定义Lyap
unov函数为犞^=犞狊+犞狅=12
狊2+狏T犘狏(30)对式(30)求导,则有:犞^·=犞·狊+犞·狅≤-(2犽-1)犞狊+12Δ2max-
λmin(犙)狏2+2Φ犘狏
(31)通过调节ESO和滑模控制器的参数,可保证犞^·
≤0,即可证明本文所设计的基于ESO的滑模控制器的闭环系统
是稳定的。

5 仿真结果与分析
在实际飞行过程中,弹道修正弹在不同时刻的飞行速度、转速、所受气动力等均不同,而修正组件的转速和受力也不同。

因此,为验证本文所设计的滚转通道控制器的有效性,对所建立的滚转控制系统进行仿真分析,从而得出最终结论。

假设1 弹道修正弹初始发射速度为犞=500m/s,出
炮口后弹体转速为1445rad/s,修正组件转速接近30rad/s。

图5为弹丸在全弹道内的飞行速度曲线,从图中可以看出,弹丸发射后,飞行速度先减小,30s左右弹丸飞至弹道顶点,随后速度开始增大。

图5 飞行速度曲线
Fig.5 Flightsp
eedcurce图6为出炮口后弹体和修正组件的转速曲线,从图中可以看出,弹丸发射后,弹体转速随时间衰减,但是仍然维持高速旋转;修正组件在炮管内运动时受到摩擦力作用而转速增大,在出炮口后同时受到气动阻力力矩影响,在全弹
道内维持约30rad/s左右的低速旋转运动。

图6 修正组件转速曲线Fig.6 Sp
inratecurveofthecorrectionfuse 
修正组件主要物理参数如表1所示。

 第4期
郑秋实等:可旋转翼式弹道修正组件滚转控制技术研究
·1417 ·
 表1 物理参数
犜犪犫犾犲1 犘犺狔狊犻犮犪犾狆犪狉犪犿犲狋犲狉狊参数取值参考面积犛/m20.01887
参考长度犔/m0.155
组件滚转惯量犐狓犳/(kg·m2)0.012舵机时间常数犜犛
0.0083
修正组件主要气动参数包括组件气动滚转阻尼力矩系数导数犆lpf和差动舵控制力矩系数导数犆犕δ狓。

在飞行过程中,这两个系数不是固定不变的,而是与弹丸的飞行马赫数有关,表2给出了这两个气动参数随飞行马赫数的变化情况,图7给出了在参考弹道中这两个气动参数随时间的变化曲线。

表2 气动参数犜犪犫犾犲2 犃犲狉狅犱狔
狀犪犿犻犮狆犪狉犪犿犲狋犲狉狊参数飞行马赫数0.60.91.11.52.0
犆lpf-0.6050-1.4090-1.4210-0.55800.3730犆犕δ狓-0.0016-0.0038-0.0042-0.0015-0.0010
图7 气动参数曲线
Fig.7 Aerody
namicparametercurves5.1 基于犈犛犗的滑模控制器性能仿真
弹道修正弹在弹道前半段,其外部环境和弹体本身运动状态变化较为剧烈,因此一般采取弹道顶点作为制导控制的起控点。

本文取飞行时间30s作为仿真起始点,仿真起始时刻弹体速度为210m/s,修正组件滚转角为50°,转速为37rad/s。

起控后,修正组件滚转角首先回到零位,然后开始跟踪幅值为±40°的阶跃指令。

在仿真开始后第18s,假设修正组件受到外部干扰力矩作用,产生的转动加速度干扰量
为犱ext=300rad/s
2,将其按照下式转换到输入端,即可得到式(11)
中的等效复合扰动值:犱eq=犱ext·犜-1狊+犱·
ext
(32)由于施加的是常值外部干扰,有犱·
ext=0,
式(32)经计算后可得犱eq=3.6×104rad/s
3。

为保证ESO估计信号的快速性和对噪声的低敏感性,
ESO的带宽折中取为ωE=80Hz,则经计算可得ESO的增益:
β1=4×2π×ωE=2.010×103β2=6×2π×ωE×ωE=1.515×106β3=4×2π×ωE×ωE×ωE=5.080×108β4=2π×ωE×ωE×ωE×ωE=6.383×10烅烄烆10(33)综合考虑控制系统稳态和动态性能,保证滚转控制系
统调节时间小于0.5s,超调量小于5%,稳态误差小于0.1°,滑模控制器参数经调试后配置为犮1=90,犮2=17,犽=200。

仿真结果如图8~图1
2所示。

图8 滚转角指令、跟踪和观测值曲线
Fig.8 Curvesofrollanglecommand,tracking,a
ndobservedvalue
图10 系统状态狓3曲线
Fig.10 Sy
stemstatecurvesof狓3
 ·1418 ·系统工程与电子技术第46卷
图11 干扰量及其观测值曲线
Fig.11 Curvesofdisturbanceanditsobservedvalue
 第4期郑秋实等:可旋转翼式弹道修正组件滚转控制技术研究·1419 ·
图16 系统状态狓2曲线
Fig.16 Systemstatecurvesof狓2
图17 系统状态狓3曲线
Fig.17 Systemstatecurvesof狓3
图18 干扰量及其观测值曲线
Fig.18 Curvesofdisturbanceanditsobservedvalue
图19 控制量曲线
Fig.19 Controlquantitycurve
5.3 与其他控制器性能对比
为验证本文设计的ESO滑模控制器的性能,对系统模
型分别采用普通滑模控制和传统PID控制,并将控制结果
与本文控制器做对比研究,仿真结果如图20所示。

 ·1
420 ·系统工程与电子技术
第46卷
参考文献
[1]REGANFJ,SMITHJ.Aeroballisticsofaterminallycorrectedspinningprojectile(TCSP)[J].JournalofSp
acecraftandRockets,1975,12(12):733738.
[2]王中原,史金光,常思江,等.弹道修正弹技术发展综述[J].
弹道学报,2021,33(2):112.
WANGZY,SHIJG,CHANGSJ,etal.Reviewondevelop
mentoftechnologyoftrajectorycorrectionproj
ectile[J].JournalofBallistics,2021,33(2):112.
[3]邢炳楠,杜忠华,杜成鑫.二维弹道修正弹及其制导控制技术综述[J].国防科技大学学报,2021,43(4):5368.XINGBN,DUZH,DUCX.Reviewontwo dimensionaltra jectorycorrectionprojectileanditsguidanceandcontroltechnology[
J].JournalofNationalUniversityofDefenseTechnology,
2021,43(4):5368.
[4]PARKW,RYOOCK,KIMY,etal.Anewpracticalguidancelawforaguidedprojectile[C]∥Proc.oftheAIAAGuidance,Navigation,andControlConference,2011.[5]TIPANS,THEODOULISS,THAIS,etal.Nonlineardynamicinversionflightcontroldesignforguidedproj
ectiles[J].JournalofGuidance,Control,andDynamics,2020,43(12)
:16.[6]GAGNONE,LAUZONM.Coursecorrectionfuzeconceptanalysisforin service155mmspin stabilizedgunneryprojectiles[C]∥Proc.oftheAIAAGuidance,Navigation&ControlConference&Exhibit,2013.[7]柯知非,宋卫东.二维弹道修正组件发展现状及关键技术[J].飞航导弹,2018(5):8185.KEZF,SONGWD.Developmentstatusandkeytechnologiesoftwo dimensionaltrajectorycorrectionprojectiles[J].Aerody
namicMissileJournal,2
018(5):8185.[8]程杰,于纪言,王晓鸣,等.隔转鸭舵式弹道修正弹电磁执行机构工况研究[J].兵工学报,2014,35(12):20102015.CHENGJ,YUJY,WANGXM,etal.Researchonworkingconditionofelectromagneticactuatoroftrajectorycorrectionproj
ectilewithdecoup
ledcanards[J].ActaArmamentarii,2004,35(12):20102015.[9]程杰,王晓鸣,于纪言,等.隔转鸭舵式弹道修正弹双旋通道参
数辨识[J].兵工学报,2016,37(10):18121819.
CHENGJ,WANGXM,YUJY,etal.Parameterestimationofaxialdual spinsysteminatrajectorycorrectionproj
ectilewithdecou p
ledcanards[J].ActaArmamentarii,2016,37(10):18121819.[10]殷婷婷,贾方秀,于纪言,等.双旋弹丸弹道修正组件控制响应模型研究[J].北京理工大学学报,2018,38(2):183189.YINTT,JIAFX,YUJY,etal.Rollcontrolresponsemodelofthedual spintrajectoryc
orrectioncomponent[J].TransactionsofBeijingInstituteofTechnology,2
018,38(2):183189.[11]殷婷婷,贾方秀,于纪言,等.基于扩张状态观测器的双旋弹
丸舵翼转速预测控制[J].北京理工大学学报,2019,39(10)
:10571062,1068.YINTT,JIAFX,YUJY,etal.Directmodelpredictiverota
tingspeedcontrolforrudderindual spinprojectilewithoutp
ut feedbackextendedstateobserver[J].TransactionsofBeijingI
nsti tuteofTechnology,2
019,39(10):10571062,1068.[12]张鑫,姚晓先.固定翼双旋弹修正组件滚转控制研究[J].北京理工大学学报,2020,40(4):386395.
ZHANGX,YAOXX.Ro
llcontrolofcoursecorrectionfuzefordual spinproj
ectilewithfixed canards[J].TransactionsofBeijingInstituteofTechnology,2
020,40(4):386395.[13]ILGMD.Guidance,navigation,andcontrolformunitions[D].Philadelphia:DrexelUniversity,2008.[14]朱大林.双旋弹飞行特性与制导控制方法研究[D].北京:北京理工大学,2015.ZHUDL.Researchonflightcharacteristics,guidance,andcontrolforadual spinprojectile[D].Beijing:BeijingI
nstituteofTechnology,2
015.[15]ZHUDL,TANGSJ,GUOJ,etal.Flightstabilityo
fadual spinprojectilewithcanards[J].Proceeding
softheInstitutionofMechanicalEngineers,PartG:JournalofAerospaceEng
i neering,2
015,229(4):703716.[16]THEODOULISS,GASSMANV,BRUNNERT,etal.Robustbank to turnautopilotdesignforaclassof155mmspin stabi lizedcanard guidedprojectiles[C]∥Proc.oftheAIAAAtmos phericFlightMechanicsConference,2013.[17]SALMANMU,CHANGB,ILGMD.Controlofspinningflyingvehiclesutilizingpartialbodydespintechniques[C]∥Proc.ofthe7thAsianControlConference,2009.[18]THEODOULISS,GASSMANV,BRUNNERT,etal.Fixedstructurerobustcontroldesignforthe155mmcanard guidedprojectileroll channelautopilot[C]∥Proc.ofthe21stMediter raneanConferenceonControl&AutomationPlatanias Chania,2013.[19]刘金琨.滑模变结构控制MATLAB仿真[M].北京:清华大学出版社,2019.
LIUJK.Slidingmodecontroldesig
nandMATLABsimulation[M].Beijing:Ts
inghuaUniversityPress,2019.[20]刘金琨,孙富春.滑模变结构控制理论及其算法研究进展[J].控制理论与应用,2007,24(3):407418.LIUJK,SUNFC.Researchanddevelopmentontheorya
ndalgorithmsofslidingmodecontrol[J].ControlTheory&Ap p
lications,2007,24(3):407418.[21]UTKINVI.Variablestructuresystemswithslidingmo
des[J].IEEETrans.onAutomaticControl,1977,22(2):212222.
[22]YUX,KAYNAKO.Sliding modecontrolwithsoftcomputing:a
survey[
J].IEEETrans.onIndustrialElectronics,2009,56(9):32753285.[23]LIUL,DINGSH,YUXH.Second orderslidingmodecon troldesignsubjecttoanasymmetricoutputconstraint[J].IEEETrans.onCircuitsandSystemsII ExpressBriefs,2021,68(4):12781282.[24]LAGHROUCHES,HARMOUCHEM,CHITOURY,etal.Barrierfunction basedadaptivehigherorderslidingmo
decon trollers[
J].Automatica,2021,123:109355.[25]韩京清.一类不确定对象的扩张状态观测器[J].控制与决策,
 第4期
郑秋实等:可旋转翼式弹道修正组件滚转控制技术研究
·1421 ·
 1995,10(1):8588.HANJQ.Extendedstateobserverforaclassofuncertain

lants[J].ControlandDecision,1995,10(1):8588.[26]HANJQ,HUANGY.Analysisanddesignforthesecondordernonlinearcontinuousextendedstatesobserver[J].Chinese
ScienceBulletin,2000,21(3):19381944.
[27]LIUCQ,LUOGZ,DUANXL,etal.AdaptiveLADRC baseddisturbancerej
ectionmethodforelectromechanicalservosystem[J].IEEETrans.onIndustryApp
lications,2020,56(1):876889.[28]LIUCQ,LUOGZ,CHENZ,etal.Measurementdelaycom pensatedLADRCbasedcurrentcontrollerdesignforPMSMdriveswithasimpleparametertuningmethod[J].ISATrans actions,2020,101:482492.[29]LUWQ,TANGB,JIKH,etal.Anewloadadaptiveidenti ficationmethodbasedonanimprovedslidingmodeobserverforPMSMpositionservosy
stem[J].IEEETrans.onPowerElec tronics,2021,36(3):32113223.
[30]WANGWW,GAOZQ.Acomparisonstudyofadvancedstateobserverdesigntechniques[C]∥Proc.oftheAmericanControlConference,2003.[31]IOANNOUPA,SUNJ.Robustadaptivecontrol[M].UpperSaddleRiver:PTRPrentice Hall,1996.作者简介
郑秋实(1993—),男,工程师,博士,主要研究方向为飞行器制导与控制。

许伟春(1964—)
,男,工程师,本科,主要研究方向为飞行器制导与控制。

赵明翰(1997—)
,男,工程师,硕士,主要研究方向为飞行器制导与控制。

李乃星(1993—),男,工程师,硕士,主要研究方向为伺服控制系统。

包旭馨(1995—),女,工程师,硕士,主要研究方向为惯性导航、伺服控制系统。

相关文档
最新文档