CN112859617A - 一种iPI无模型自适应全局非奇异快速终端滑模控制方法 - Google Patents
一种iPI无模型自适应全局非奇异快速终端滑模控制方法 Download PDFInfo
- Publication number
- CN112859617A CN112859617A CN202110178099.4A CN202110178099A CN112859617A CN 112859617 A CN112859617 A CN 112859617A CN 202110178099 A CN202110178099 A CN 202110178099A CN 112859617 A CN112859617 A CN 112859617A
- Authority
- CN
- China
- Prior art keywords
- model
- sliding mode
- terminal sliding
- ipi
- control method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 116
- 230000003044 adaptive effect Effects 0.000 title claims description 59
- 238000004088 simulation Methods 0.000 claims abstract description 51
- 238000013459 approach Methods 0.000 claims abstract description 20
- 238000004134 energy conservation Methods 0.000 claims abstract description 11
- 238000013178 mathematical model Methods 0.000 claims abstract description 9
- 230000003313 weakening effect Effects 0.000 claims abstract description 4
- 230000005855 radiation Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 11
- 230000002457 bidirectional effect Effects 0.000 claims description 5
- 230000004069 differentiation Effects 0.000 claims description 4
- 230000006978 adaptation Effects 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 3
- 230000000191 radiation effect Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims 4
- 239000002994 raw material Substances 0.000 claims 4
- 238000013461 design Methods 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 13
- 239000010453 quartz Substances 0.000 description 12
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 12
- 238000010438 heat treatment Methods 0.000 description 7
- 239000000463 material Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 229910018487 Ni—Cr Inorganic materials 0.000 description 1
- HMDDXIMCDZRSNE-UHFFFAOYSA-N [C].[Si] Chemical compound [C].[Si] HMDDXIMCDZRSNE-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- VNNRSPGTAMTISX-UHFFFAOYSA-N chromium nickel Chemical compound [Cr].[Ni] VNNRSPGTAMTISX-UHFFFAOYSA-N 0.000 description 1
- 229910002804 graphite Inorganic materials 0.000 description 1
- 239000010439 graphite Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种iPI无模型自适应全局非奇异快速终端滑模控制方法,包括,根据能量守恒定律,建立高超声速飞行器气动热地面模拟系统数学模型及无模型控制的超局部模型;利用非线性ESO观测器对所述高超声速飞行器气动热地面模拟系统进行未知扰动预测;基于全局非奇异快速终端滑模面削弱趋近抖振、收敛速度慢、奇异问题;根据滑模可达性条件定义等效控制率和自适应趋近律,得到滑模控制率,完成目标跟踪。本发明通过全局非奇异快速终端滑模面的设计,去除了趋近模态,从而减缓了抖振现象,加快了在滑动模态上的控制速度,对滑模面条件的限制使得控制器不存在奇异问题,结合了自适应的方法有效解决收敛停滞问题。
Description
技术领域
本发明涉及航空航天自动化的技术领域,尤其涉及一种iPI无模型自适应全局非奇异快速终端滑模控制方法。
背景技术
石英灯加热器作为地面热试验模拟的一种较为常见方式,相较于镍铬棒、碳硅棒、石墨棒作为加热元件的热辐射试验中,它具有使用寿命长、热效率高、升温速度快、结构紧凑、可控性能强等特点,因而被广泛应用到高超声速飞行器气动热地面模拟试验中。作为模拟飞行器热环境的石英灯加热器,要求其具有良好的稳定性、准确性、快速性和抗干扰的能力,这就需要设计控制方法来达到这些性能指标的要求。
传统的控制方法可以对石英灯加热器输出温度进行跟踪,但需要精确的石英灯加热器模型参数。而石英灯加热器由于外部的未知扰动和自身内部的不确定项的影响,使得传统基于模型的控制方式不能够应用在石英灯加热器上。因此提出了无模型的控制方式作用在石英灯加热器上,超局部模型是无模型控制的一种方式,但是仅适用于单输入单输出的系统中且要求转化后的被控对象表达式只能是一阶或二阶的。采用无模型控制中的超局部模型结合ESO扩展状态观测器将整个单输入单输出的石英灯加热器的数学模型进行线性化处理,所有的未知扰动包括确定的、不确定的、未知的、已知的都通过ESO扩展观测器来跟踪。
另外,传统的PID控制的优点是结构简单,易于调试,但是控制的精度不够高,收敛速度较慢,超调较大。而作为现代控制理论中的滑模控制不论是在响应速度、抗干扰能力还是控制精度上都优于传统的PID控制,并且其对模型精确性的要求不高,因此被广泛应用。其中,终端滑模除了具有传统滑模控制的特点,还能够在有限时间收敛。但因为终端滑模也有趋近和滑动两个阶段,存在趋近阶段抖动、收敛停滞、控制奇异等问题。
发明内容
本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊,而这种简化或省略不能用于限制本发明的范围。
鉴于上述现有存在的问题,提出了本发明。
因此,本发明提供了一种iPI无模型自适应全局非奇异快速终端滑模控制方法,能够解决无法在有限时间内跟踪期望温度、削弱控制过程中的抖振现象的问题。
为解决上述技术问题,本发明提供如下技术方案:包括,根据能量守恒定律,建立高超声速飞行器气动热地面模拟系统数学模型及无模型控制的超局部模型;利用非线性ESO观测器对所述高超声速飞行器气动热地面模拟系统进行未知扰动预测;基于全局非奇异快速终端滑模面削弱趋近抖振、收敛速度慢、奇异问题;根据滑模可达性条件定义等效控制率和自适应趋近律,得到滑模控制率,完成目标跟踪。
作为本发明所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的一种优选方案,其中:所述高超声速飞行器气动热地面模拟系统包括,非接触辐射加热器、电功率调节装置和量热传感器;根据所述能量守恒定律建立输入输出能量守恒等式,得到当前温度T1和双向晶闸管的导通角α之间的数学关系,如下,
其中,等式左边UI为输入电压即电源两端电压,R为非接触辐射加热器的电阻之和,α为双向晶闸管的导通角,等式右边分别为用于非接触辐射加热器自身消耗的内能、对流换热过程中损失的热能、热传导过程中损失的热能、热辐射效应输出的热能,c、m、T1、T0、A、ε、Δt分别为非接触辐射加热器的比热容、质量、当前温度、初始温度、表面积、黑度系数、工作时间,β、λ、σ、F分别为对流换热系数、导热系数、斯蒂芬-玻尔兹曼常数、角系数。
作为本发明所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的一种优选方案,其中:包括,当被控对象模型是单输入单输出系统时,则将所述被控对象模型转变成所述无模型控制的超局部模型,如下,
y(n)=G+χu(t)
其中,y(n)表示为输出量y对时间t的n阶导数,n一般取1或者2,u表示为输入量,G表示为所有未知扰动的集合,既包含了外界扰动和系统内部非线性扰动,χ表示为非物理意义的可调参数。
作为本发明所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的一种优选方案,其中:包括,根据所述无模型控制的超局部模型,将所述输入输出能量守恒等式两边除以Δt并进行移项,得到所述高超声速飞行器气动热地面模拟系统的数学模型,如下,
其中,为T1对时间Δt的导数,α分别对应所述无模型控制的超局部模型中的y(n)、u;而sin2α给系统带来的是周期性的震动,并没有对系统整体的收敛产生影响,含有sin2α的项可以看作输入扰动,AεσFT1 4可以看作系统的高阶输出扰动,因此可以看作既包含输入扰动又包含输出扰动的全部扰动之和,对应于超局部模型的G,G可以通过观测器来观测。
作为本发明所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的一种优选方案,其中:还包括,
T1=x3,G=x4
e1=z1-x3,e2=z2-x4
β1>0,β2>0
作为本发明所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的一种优选方案,其中:还包括,定义输出的跟踪误差表达式如下,
e(t)=y*-y
其中,e为跟踪误差,y*为输出目标;根据所述无模型控制的超局部模型,通过闭环控制得到无模型控制器,如下,
其中,uaux是依据全局非奇异快速终端滑模控制。
作为本发明所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的一种优选方案,其中:建立所述全局非奇异快速终端滑模面,包括,
η>0,ι>0
作为本发明所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的一种优选方案,其中:还包括,建立跟踪误差e与辅助控制器uaux之间的数学关系,如下,
作为本发明所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的一种优选方案,其中:包括,联立融合所述辅助控制器、所述等效控制器和所述趋近率,得到所述高超声速飞行器气动热地面模拟系统的iPI无模型自适应全局非奇异快速终端滑模控制器u(t),如下,
其中,uaux=ueq+ucor。
本发明的有益效果:本发明将全局非奇异快速终端滑模控制方法结合无模型控制超局部模型应用到高超声速飞行器气动热地面模拟系统中,将高超声速飞行器气动热地面模拟系统数学模型线性化处理,利用非线性ESO观测器对全部的未知扰动进行观测;同时,全局非奇异快速终端滑模面的设计,去除了趋近模态,使整个控制阶段都在滑动模态上进行,从而减缓了抖振现象,加快了在滑动模态上的控制速度,对滑模面条件的限制使得控制器不存在奇异问题,结合了自适应的方法有效解决收敛停滞问题。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。其中:
图1为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的流程示意图;
图2为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的高超声速飞行器气动热地面模拟系统工作流程示意图;
图3(a)为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的高超声速导弹三维结构示意图;
图3(b)为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的高超声速导弹二维尺寸示意图;
图4(a)为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的高超声速导弹气动热的有限元仿真示意图;
图4(b)为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的高超声速导弹气动热的壁面平均温度采样示意图;
图4(c)为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的高超声速导弹气动热的数据拟合示意图;
图5为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的控制框架示意图;
图6为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的高超声速飞行器气动热地面模拟系统在iPI无模型自适应全局非奇异快速终端滑模控制方法(1)、iPI无模型自适应全局终端滑模控制方法(2)、iPI控制方法(3)和传统PID方法(4)下的输出温度曲线图(a)和局部放大图(b);
图7为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的高超声速飞行器气动热地面模拟系统在iPI无模型自适应全局非奇异快速终端滑模控制方法(1)、iPI无模型自适应全局终端滑模控制方法(2)、iPI控制方法(3)和传统PID方法(4)下的跟踪误差曲线图(a)和局部放大图(b);
图8为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的外扰(a)下高超声速飞行器气动热地面模拟系统在iPI无模型自适应全局非奇异快速终端滑模控制方法(1)、iPI无模型自适应全局终端滑模控制方法(2)、iPI控制方法(3)和传统PID方法(4)下的输出温度曲线图(b)和局部放大图(c);
图9为本发明一个实施例所述的iPI无模型自适应全局非奇异快速终端滑模控制方法的外扰下高超声速飞行器气动热地面模拟系统在iPI无模型自适应全局非奇异快速终端滑模控制方法(1)、iPI无模型自适应全局终端滑模控制方法(2)、iPI控制方法(3)和传统PID方法(4)下的跟踪误差曲线图(a)和局部放大图(b)。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合说明书附图对本发明的具体实施方式做详细的说明,显然所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明的保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
其次,此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例,也不是单独的或选择性的与其他实施例互相排斥的实施例。
本发明结合示意图进行详细描述,在详述本发明实施例时,为便于说明,表示器件结构的剖面图会不依一般比例作局部放大,而且所述示意图只是示例,其在此不应限制本发明保护的范围。此外,在实际制作中应包含长度、宽度及深度的三维空间尺寸。
同时在本发明的描述中,需要说明的是,术语中的“上、下、内和外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一、第二或第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
本发明中除非另有明确的规定和限定,术语“安装、相连、连接”应做广义理解,例如:可以是固定连接、可拆卸连接或一体式连接;同样可以是机械连接、电连接或直接连接,也可以通过中间媒介间接相连,也可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1
参照图1~图5,为本发明的第一个实施例,提供了一种iPI无模型自适应全局非奇异快速终端滑模控制方法,本发明方法基于高超声速飞行器气动热地面模拟系统无模型控制的超局部模型,结合iPI、非线性ESO观测器、全局非奇异快速终端滑模面、等效控制率、自适应趋近律,设计控制器u(t),实现目标跟踪;参照图5,为本发明高超声速飞行器气动热地面模拟系统iPI无模型自适应全局非奇异快速终端滑模控制框图,具体包括:
S1:根据能量守恒定律,建立高超声速飞行器气动热地面模拟系统数学模型及无模型控制的超局部模型。其中需要说明的是,高超声速飞行器气动热地面模拟系统包括:
非接触辐射加热器、电功率调节装置和量热传感器;
根据能量守恒定律建立输入输出能量守恒等式,得到当前温度T1和双向晶闸管的导通角α之间的数学关系,如下,
其中,等式左边UI为输入电压即电源两端电压,R为非接触辐射加热器的电阻之和,α为双向晶闸管的导通角,等式右边分别为用于非接触辐射加热器自身消耗的内能、对流换热过程中损失的热能、热传导过程中损失的热能、热辐射效应输出的热能,c、m、T1、T0、A、ε、Δt分别为非接触辐射加热器的比热容、质量、当前温度、初始温度、表面积、黑度系数、工作时间,β、λ、σ、F分别为对流换热系数、导热系数、斯蒂芬-玻尔兹曼常数、角系数;
当被控对象模型是单输入单输出系统时,则将被控对象模型转变成无模型控制的超局部模型,如下,
y(n)=G+χu(t)
其中,y(n)表示为输出量y对时间t的n阶导数,n一般取1或者2,u表示为输入量,G表示为所有未知扰动的集合,既包含了外界扰动和系统内部非线性扰动,χ表示为非物理意义的可调参数。
S2:利用非线性ESO观测器对高超声速飞行器气动热地面模拟系统进行未知扰动预测。本步骤需要说明的是:
根据无模型控制的超局部模型,将输入输出能量守恒等式两边除以Δt并进行移项,得到高超声速飞行器气动热地面模拟系统的数学模型,如下,
其中,为T1对时间Δt的导数,α分别对应无模型控制的超局部模型中的y(n)、u;而sin2α给系统带来的是周期性的震动,并没有对系统整体的收敛产生影响,含有sin2α的项可以看作输入扰动,AεσFT1 4可以看作系统的高阶输出扰动,因此可以看作既包含输入扰动又包含输出扰动的全部扰动之和,对应于超局部模型的G,G可以通过观测器来观测;
T1=x3,G=x4
对高超声速飞行器气动热地面模拟系统进行未知扰动预测,建立非线性ESO观测器,如下,
e1=z1-x3,e2=z2-x4
β1>0,β2>0
定义输出的跟踪误差表达式如下,
e(t)=y*-y
其中,e为跟踪误差,y*为输出目标;
根据无模型控制的超局部模型,通过闭环控制得到无模型控制器,如下,
为了削弱观测扰动,在无模型控制器上添加辅助控制器uaux,如下,
其中,uaux是依据全局非奇异快速终端滑模控制。
S3:基于全局非奇异快速终端滑模面削弱趋近抖振、收敛速度慢、奇异问题。其中还需要说明的是,建立全局非奇异快速终端滑模面,包括:
η>0,ι>0
建立跟踪误差e与辅助控制器uaux之间的数学关系,如下,
S4:根据滑模可达性条件定义等效控制率和自适应趋近律,得到滑模控制率,完成目标跟踪。本步骤还需要说明的是:
联立融合辅助控制器、等效控制器和趋近率,得到高超声速飞行器气动热地面模拟系统的iPI无模型自适应全局非奇异快速终端滑模控制器u(t),如下,
其中,uaux=ueq+ucor;
进一步的,建立Lyapunov函数稳定性判据表达式,验证iPI无模型自适应全局非奇异快速终端滑模控制方法的收敛性,如下:
参照图2,高超声速飞行器气动热地面模拟系统工作流程主要包括以下步骤:
(1)高超声速飞行器气动热数据采集:对高超声速导弹,通过给定飞行环境和壁面材料型号对其进行有限元数值模拟;采集每一时刻导弹壁面的平均温度,对采样数据进行线性拟合,作为整个高超声速飞行器气动热地面模拟系统的期望输出值,即目标值,为了和实际控制器输出值进行对比。
(2)高超声速飞行器气动热地面模拟控制系统:设计控制器对石英灯加热系统进行控制;将目标值加载到控制板中,通过控制板改变双向晶闸管的导通角α进而改变输出电压U,不同的输出电压U值对应不同的石英灯加热系统电功率P,通过传感器获得石英灯加热器输出的实际温度T1,与目标值对比得到跟踪误差e,再通过闭环反馈给控制器进行调节双向晶闸管的导通角α,最终达到跟踪控制。
(3)地面模拟试验反馈:石英灯加热器对试件进行加热试验,再通过对试件的性能进行检测,分析材料的可行性,选定材料,若不可以就更换材料再次经过第一步操作,从而优化热防护系统设计。
参照图3,是有限元仿真所绘制的高超声速导弹,导弹的具体参数为:总长7600mm,弹体长4270mm,弹体直径1168.4mm,制导部分夹角7°,导引头半径30mm,夹角12.84°,飞行环境为高度32km,速度为6.0马赫数,攻角0°巡航。
参照图4,是导弹进行有限元仿真的仿真模拟图、壁面平均温度采样图和平均温度曲线拟合图,拟合曲线为:
y*=1.848*10-5t8-0.001497t7+0.04917t6-0.8402t5
+7.977t4-41.34t3+93.82t2+151t+287.4
参照图5,是高超声速飞行器气动热地面模拟系统iPI无模型自适应全局非奇异快速终端滑模控制框图,是对图2的第二步高超声速飞行器气动热地面模拟控制系统的进一步的说明,根据图5的示意,控制器u(t)由4部分组成:1、自适应全局非奇异快速终端滑模面提供等效控制率和自适应趋近律;2、无模型控制iPI闭环反馈控制率;3、期望目标的一阶微分;4、利用非线性ESO观测器对系统未知扰动G的观测值
优选的,本实施例还需要说明的是,与现有技术相比,本实施例公开了一种iPI无模型自适应全局非奇异快速终端滑模控制方法,旨在通过采用基于iPI无模型自适应全局非奇异快速终端滑模控制方法能够在有限时间内跟踪期望温度,并且削弱控制过程中的抖振现象,加快收敛速度,保证控制全过程的鲁棒性能;其中,采用无模型控制中的超局部模型结合非线性ESO扩展状态观测器,实现对系统线性化;采用全局非奇异快速终端滑模面、等效控制率、自适应趋近律解决了存在趋近阶段抖动、收敛停滞、控制奇异等问题。
实施例2
参照图6~图9,为本发明的第二个实施例,该实施例不同于第一个实施例的是,提供了一种iPI无模型自适应全局非奇异快速终端滑模控制方法的测试对比验证,具体包括:
本实施例中将采用高超声速飞行器气动热地面模拟系统在iPI无模型自适应全局非奇异快速终端滑模控制方法(1)、iPI无模型自适应全局终端滑模控制方法(2)、iPI控制方法(3)和传统PID方法(4)下分别对高超声速飞行器气动热地面模拟系统的输出温度和跟踪误差进行实时测量对比。
测试环境:参照图4(c),将高超声速飞行器气动热地面模拟系统运行在仿真平台模拟跟踪期望目标曲线,分别利用高超声速飞行器气动热地面模拟系统在iPI无模型自适应全局非奇异快速终端滑模控制方法(1)、iPI无模型自适应全局终端滑模控制方法(2)、iPI控制方法(3)和传统PID方法(4)下进行测试并获得测试结果数据;四种方法都将在开启自动化测试设备并运用MATLB软件编程实现对比方法的仿真测试,根据试验结果得到仿真数据,每种方法各测试4组数据,每组数据采样15s,计算获得每组数据输入温度和跟踪误差,与仿真模拟输入的期望目标温度进行对比计算误差。
参照图6、图7、图8和图9,是高超声速飞行器气动热地面模拟系统在iPI无模型自适应全局非奇异快速终端滑模控制方法(1)、iPI无模型自适应全局终端滑模控制方法(2)、iPI控制方法(3)和传统PID方法(4)下的输出温度曲线图和局部放大图、误差跟踪曲线对比图和局部放大图,以及外扰动下的误差跟踪曲线对比图及局部放大图。
参照图7(a),为外扰动为时变电阻R,如下:
R=3.08×(1+0.0045y*)
iPI无模型自适应全局终端滑模控制方法(2):
iPI控制方法(3):
具体实施例各项参数设置如下表所示:
表1:高超声速飞行器气动热地面模拟系统参数表。
表2:iPI无模型自适应全局非奇异快速终端滑模控制方法数据表。
表3:iPI无模型自适应全局终端滑模控制方法数据表。
表4:iPI和PID数据表。
参照图6,能够直观的看出,4种方法都可以有效地跟踪目标曲线,在图6(b)中,当时间为0~0.2s时,方法(3)和方法(4)出现了较大的超调量。
参照图7,能够直观的看出,方法(2)前期出现了一定的稳态误差,且抖振较为明显;方法(3)的快速性较差以及抖振较为明显;方法(4)出现了较大的超调量。
参照图8,能够直观的看出,在外扰的情况下,4种方法大体上可以跟踪目标曲线,在图8(c)中,当时间0~0.2s时,方法3和方法(4)前期未能有效的跟踪目标曲线,存在较大的超调量。
参照图9,能够直观的看出,在外扰的情况下,方法(2)前期的控制精度较差,存在一定的误差;方法(3)的快速性较差以及后期高频切换幅度较大;方法(4)的超调量大且快速性慢。
参照图6~图9,能够分析出,本发明的控制方法从稳态误差、快速性、超调量、控制精度4个方面均优于其它3种方法,得益于本发明一种iPI无模型自适应全局非奇异快速终端滑模控制方法削弱控制过程中的抖振现象,加快收敛速度,保证控制全过程的鲁棒性能。
应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (9)
1.一种iPI无模型自适应全局非奇异快速终端滑模控制方法,其特征在于:包括,
根据能量守恒定律,建立高超声速飞行器气动热地面模拟系统数学模型及无模型控制的超局部模型;
利用非线性ESO观测器对所述高超声速飞行器气动热地面模拟系统进行未知扰动预测;
基于全局非奇异快速终端滑模面削弱趋近抖振、收敛速度慢、奇异问题;
根据滑模可达性条件定义等效控制率和自适应趋近律,得到滑模控制率,完成目标跟踪。
2.根据权利要求1所述的iPI无模型自适应全局非奇异快速终端滑模控制方法,其特征在于:所述高超声速飞行器气动热地面模拟系统包括,非接触辐射加热器、电功率调节装置和量热传感器;
根据所述能量守恒定律建立输入输出能量守恒等式,得到当前温度T1和双向晶闸管的导通角α之间的数学关系,如下,
其中,等式左边UI为输入电压即电源两端电压,R为非接触辐射加热器的电阻之和,α为双向晶闸管的导通角,等式右边分别为用于非接触辐射加热器自身消耗的内能、对流换热过程中损失的热能、热传导过程中损失的热能、热辐射效应输出的热能,c、m、T1、T0、A、ε、Δt分别为非接触辐射加热器的比热容、质量、当前温度、初始温度、表面积、黑度系数、工作时间,β、λ、σ、F分别为对流换热系数、导热系数、斯蒂芬-玻尔兹曼常数、角系数。
3.根据权利要求1或2所述的iPI无模型自适应全局非奇异快速终端滑模控制方法,其特征在于:包括,
当被控对象模型是单输入单输出系统时,则将所述被控对象模型转变成所述无模型控制的超局部模型,如下,
y(n)=G+χu(t)
其中,y(n)表示为输出量y对时间t的n阶导数,n一般取1或者2,u表示为输入量,G表示为所有未知扰动的集合,既包含了外界扰动和系统内部非线性扰动,χ表示为非物理意义的可调参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110178099.4A CN112859617B (zh) | 2021-02-08 | 2021-02-08 | 一种iPI无模型自适应全局非奇异快速终端滑模控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110178099.4A CN112859617B (zh) | 2021-02-08 | 2021-02-08 | 一种iPI无模型自适应全局非奇异快速终端滑模控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112859617A true CN112859617A (zh) | 2021-05-28 |
CN112859617B CN112859617B (zh) | 2021-12-07 |
Family
ID=75989759
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110178099.4A Active CN112859617B (zh) | 2021-02-08 | 2021-02-08 | 一种iPI无模型自适应全局非奇异快速终端滑模控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112859617B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114355779A (zh) * | 2022-01-10 | 2022-04-15 | 南京工业大学 | 一种结构热试验非线性扩展状态观测器全局滑模控制方法 |
CN114371623A (zh) * | 2022-01-10 | 2022-04-19 | 南京工业大学 | 一种结构热试验神经网络观测器分数阶全局滑模控制方法 |
CN115598978A (zh) * | 2022-10-12 | 2023-01-13 | 北京临近空间飞行器系统工程研究所(Cn) | 一种高速飞行器全局快速非奇异终端滑模姿态控制方法 |
CN118244643A (zh) * | 2024-04-25 | 2024-06-25 | 电子科技大学 | 一种伺服位置跟踪系统的无模型自适应滑模控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105116725A (zh) * | 2015-07-17 | 2015-12-02 | 浙江工业大学 | 基于扩张状态观测器的伺服系统自适应滑模控制方法 |
CN109194219A (zh) * | 2018-09-25 | 2019-01-11 | 湖南工业大学 | 基于无模型非奇异终端滑模控制永磁同步电机方法及系统 |
CN109343549A (zh) * | 2018-10-09 | 2019-02-15 | 北京航空航天大学 | 一种飞行器姿态控制方法、系统、介质及设备 |
CN111459051A (zh) * | 2020-04-23 | 2020-07-28 | 河北工业大学 | 一种带扰动观测器的离散终端滑模无模型控制方法 |
CN111865169A (zh) * | 2020-07-21 | 2020-10-30 | 南京航空航天大学 | 一种超声电机伺服系统的无模型积分滑模控制方法 |
-
2021
- 2021-02-08 CN CN202110178099.4A patent/CN112859617B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105116725A (zh) * | 2015-07-17 | 2015-12-02 | 浙江工业大学 | 基于扩张状态观测器的伺服系统自适应滑模控制方法 |
CN109194219A (zh) * | 2018-09-25 | 2019-01-11 | 湖南工业大学 | 基于无模型非奇异终端滑模控制永磁同步电机方法及系统 |
CN109343549A (zh) * | 2018-10-09 | 2019-02-15 | 北京航空航天大学 | 一种飞行器姿态控制方法、系统、介质及设备 |
CN111459051A (zh) * | 2020-04-23 | 2020-07-28 | 河北工业大学 | 一种带扰动观测器的离散终端滑模无模型控制方法 |
CN111865169A (zh) * | 2020-07-21 | 2020-10-30 | 南京航空航天大学 | 一种超声电机伺服系统的无模型积分滑模控制方法 |
Non-Patent Citations (3)
Title |
---|
PENG GAO 等: "An Adaptive Super Twisting Nonlinear Fractional Order PID Sliding Mode Control of Permanent Magnet Synchronous Motor Speed Regulation System Based on Extended State Observer", 《IEEE ACCESS》 * |
ZHANG GUANGMING 等: "Disturbance Observer-Based Sliding Mode Control for DC-DC Boost Converter", 《2018 IEEE 3RD ADVANCED INFORMATION TECHNOLOGY, ELECTRONIC AND AUTOMATION CONTROL CONFERENCE (IAEAC)》 * |
侯明冬 等: "一种无模型自适应积分终端滑模控制方法", 《控制与决策》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114355779A (zh) * | 2022-01-10 | 2022-04-15 | 南京工业大学 | 一种结构热试验非线性扩展状态观测器全局滑模控制方法 |
CN114371623A (zh) * | 2022-01-10 | 2022-04-19 | 南京工业大学 | 一种结构热试验神经网络观测器分数阶全局滑模控制方法 |
CN114355779B (zh) * | 2022-01-10 | 2022-07-19 | 南京工业大学 | 一种结构热试验非线性扩展状态观测器全局滑模控制方法 |
CN115598978A (zh) * | 2022-10-12 | 2023-01-13 | 北京临近空间飞行器系统工程研究所(Cn) | 一种高速飞行器全局快速非奇异终端滑模姿态控制方法 |
CN118244643A (zh) * | 2024-04-25 | 2024-06-25 | 电子科技大学 | 一种伺服位置跟踪系统的无模型自适应滑模控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112859617B (zh) | 2021-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112859617B (zh) | 一种iPI无模型自适应全局非奇异快速终端滑模控制方法 | |
CN109883660B (zh) | 一种热模拟试验控制方法 | |
CN112965382B (zh) | 一种基于神经网络观测器非线性全局滑模无模型控制方法 | |
CN112987566B (zh) | 一种气动热的超螺旋非线性分数阶滑模无模型控制方法 | |
CN112987569B (zh) | 一种扰动上界自适应分数阶全局终端滑模无模型控制方法 | |
CN112965384B (zh) | 一种自适应神经网络全局非奇异快速终端滑模无模型方法 | |
CN112965364B (zh) | 一种高超声速飞行器气动热地面模拟系统建模及iPID方法 | |
CN105093931B (zh) | 一种航空发动机非线性系统控制器设计方法 | |
CN113051831B (zh) | 机床热误差自学习预测模型建模方法及热误差控制方法 | |
CN113625564B (zh) | 一种结构热试验基于系统模型的滑模控制方法 | |
CN108803336A (zh) | 一种航空发动机自适应lqg/ltr控制器设计方法 | |
CN107944137B (zh) | 高超声速飞行器弹道状态多场耦合的热气动弹性计算技术 | |
CN113495486B (zh) | 一种结构热试验基于扩展状态观测器的模型预测控制方法 | |
CN108828954B (zh) | 气候风洞自适应预测控制系统及其控制方法 | |
CN114371623B (zh) | 一种结构热试验神经网络观测器分数阶全局滑模控制方法 | |
CN114706295B (zh) | 用于空天飞机强度测试的热试验中快时变热载荷控制方法 | |
CN109298735A (zh) | 差示扫描量热仪恒速升温过程的前馈-反馈复合控制方法 | |
CN114355779B (zh) | 一种结构热试验非线性扩展状态观测器全局滑模控制方法 | |
CN115236975B (zh) | 一种导弹气动热地面有限元分析及ipd非线性滑模控制方法 | |
CN113900373B (zh) | 一种结构热试验的全局非奇异快速终端滑模控制方法 | |
CN117608197B (zh) | 一种基于自适应鲁棒滑模技术的低温风洞总温控制方法 | |
CN113065201A (zh) | 一种考虑滑移修正的辐射平衡温度计算方法 | |
Grigorie et al. | Design and experimental validation of a combined PI and Bi-positional laws controller for delaying the transition from laminar flow to turbulent flow over a morphing wing | |
Panishev et al. | Modernization of Chamber Furnace Electric Equipment | |
Song et al. | Research on Sapphire Furnace Temperature Control Strategy based on Active Disturbance Rejection Control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |