CN112987566A - 一种气动热的超螺旋非线性分数阶滑模无模型控制方法 - Google Patents

一种气动热的超螺旋非线性分数阶滑模无模型控制方法 Download PDF

Info

Publication number
CN112987566A
CN112987566A CN202110170226.6A CN202110170226A CN112987566A CN 112987566 A CN112987566 A CN 112987566A CN 202110170226 A CN202110170226 A CN 202110170226A CN 112987566 A CN112987566 A CN 112987566A
Authority
CN
China
Prior art keywords
model
fractional order
nonlinear
sliding mode
supercoiled
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
Application number
CN202110170226.6A
Other languages
English (en)
Other versions
CN112987566B (zh
Inventor
张广明
柏志青
吕筱东
高鹏
王恒强
王月
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Tech University
Original Assignee
Nanjing Tech University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Tech University filed Critical Nanjing Tech University
Priority to CN202110170226.6A priority Critical patent/CN112987566B/zh
Publication of CN112987566A publication Critical patent/CN112987566A/zh
Application granted granted Critical
Publication of CN112987566B publication Critical patent/CN112987566B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/042Adaptive 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)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开了一种气动热的超螺旋非线性分数阶滑模无模型控制方法,包括,基于能量守恒定律,建立高超声速飞行器气动热地面模拟系统输入电能与输出温度之间的数学模型并转化为无模型控制的超局部模型;根据定义的输出跟踪误差、非线性函数和分数阶微积分,构建非线性分数阶滑膜面;结合所述非线性分数阶滑膜面、超螺旋趋近率、所述超局部模型及时间延时观测器,搭建超螺旋非线性分数阶滑膜的无模型控制器,抑制控制过程中的抖振。本发明对非线性分数阶滑模面的设计既保证了控制的稳定性、收敛速度又降低了稳态误差和饱和误差,而超螺旋趋近率的结合改善了滑模面的抖动问题。

Description

一种气动热的超螺旋非线性分数阶滑模无模型控制方法
技术领域
本发明涉及航空航天自动化的技术领域,尤其涉及一种气动热的超螺旋非线性分数阶滑模无模型控制方法。
背景技术
结构热实验常见的地面模拟热试验主要有风洞试验和辐射热试验,风洞试验成本高昂,环境条件要求严苛,在过高的飞行马赫数下不能精确的模拟飞行器的热环境,而辐射热试验以对流、传导、辐射等非接触方式将热能传递给试件,有效地模拟试件在飞行过程中的时序受热;其中,因为石英灯具有升温速度快、尺寸结构小、输出功率大、操作安全的特点,所以以钨丝为热辐射元件的石英灯加热器被广泛采纳;另外作为模拟飞行器热环境的实验设备首先要保证的是稳定性、控制的精度、收敛速度这三方面的要求,这就需要设计一个符合要求的控制方法来满足。
考虑到石英灯加热器在过程控制中存在一些未知的动态扰动,提出了无模型控制的超局部模型的方法,无模型控制的超局部模型要求被控对象是单输入单输出的系统,且被控对象的数学表达式只能是一阶、二阶的系统,滑模控制作为现代控制理论不可或缺的一条分支,相较于传统的控制方式具有快速响应、结构简单、对干扰不敏感的能力被学者深入研究,但由于对收敛速度和控制精度方面的严格要求,一些结合滑模的控制算法被广泛提出;分数阶滑模控制除了具有传统滑模控制的优点之外,还具有强的鲁棒性、灵活性,但现有的线性分数阶滑模具有不能减小抖动现象和控制精度不高的缺点,不能满足被控对象的高性能指标的要求。
发明内容
本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊,而这种简化或省略不能用于限制本发明的范围。
鉴于上述现有存在的问题,提出了本发明。
因此,本发明提供了一种气动热的超螺旋非线性分数阶滑模无模型控制方法,能够解决现有高超声速飞行器气动热地面模拟系统中进行温度跟踪存在误差、控制过程中的抖振问题。
为解决上述技术问题,本发明提供如下技术方案:包括,基于能量守恒定律,建立高超声速飞行器气动热地面模拟系统输入电能与输出温度之间的数学模型并转化为无模型控制的超局部模型;根据定义的输出跟踪误差、非线性函数和分数阶微积分,构建非线性分数阶滑模面;结合所述非线性分数阶滑模面、超螺旋趋近率、所述超局部模型及时间延时观测器,搭建超螺旋非线性分数阶滑模的无模型控制器,抑制控制过程中的抖振。
作为本发明所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的一种优选方案,其中:所述高超声速飞行器气动热地面模拟系统包括,非接触辐射加热器、电功率调节装置和量热传感器;所述非接触辐射加热器为石英灯加热器,所述电功率调节装置为双向晶闸管,所述量热传感器为热电偶传感器。
作为本发明所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的一种优选方案,其中:根据所述能量守恒定律构建输入输出能量守恒等式,包括,
Figure BDA0002938680540000021
其中,等式左边UI为输入电压即电源两端电压,R为非接触辐射加热器的电阻之和,α为双向晶闸管的导通角,等式右边分别为用于非接触辐射加热器自身消耗的内能、对流换热过程中损失的热能、热传导过程中损失的热能、热辐射效应输出的热能,c、m、T1、T0、A、ε、Δt分别为非接触辐射加热器的比热容、质量、当前温度、初始温度、表面积、黑度系数、工作时间,β、λ、σ、F分别为对流换热系数、导热系数、斯蒂芬-玻尔兹曼常数、角系数,从而得到当前温度T1和双向晶闸管的导通角α之间的数学关系。
作为本发明所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的一种优选方案,其中:包括,当被控对象模型是单输入单输出系统时,则将所述被控对象模型转化为无模型控制的超局部模型,如下,
y(n)=G+χu(t)
其中,y(n)表示为输出量y对时间t的n阶导数,n一般取1或者2,u表示为输入量,G表示为所有未知扰动的集合,既包含了外界扰动和系统内部非线性扰动,χ表示为非物理意义的可调参数。
作为本发明所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的一种优选方案,其中:还包括,根据所述超局部模型将所述输入输出能量守恒等式两边除以Δt并进行移项处理,获得无模型控制的超局部控制模型,如下,
Figure BDA0002938680540000031
其中,
Figure BDA0002938680540000032
为T1对时间Δt的导数,所述超局部控制模型中
Figure BDA0002938680540000033
α分别对应所述超局部模型的y(n)、u,而sin2α给系统带来的是周期性的震动,并没有对系统整体的收敛产生影响,含有sin2α的项可以看作输入扰动,
Figure BDA0002938680540000034
可以看作系统的高阶输出扰动,因此
Figure BDA0002938680540000035
可以看作既包含输入扰动又包含输出扰动的全部扰动之和,对应于超局部模型的G。
作为本发明所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的一种优选方案,其中:包括,定义输出的跟踪误差e表达式为:
e(t)=y*-y
其中,e为跟踪误差,y*为输出目标;非线性函数fal(e,γ,η)为:
Figure BDA0002938680540000036
其中0<γ<1,η>0,
Figure BDA0002938680540000037
非线性分数阶滑模面为:
Figure BDA0002938680540000038
其中,φna>0,φnb>0,φnc>0,用于调参,
Figure BDA0002938680540000039
为分数阶积分,
Figure BDA00029386805400000310
为分数阶微分,
Figure BDA00029386805400000311
ψ为可调参数。
作为本发明所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的一种优选方案,其中:包括,定义非线性分数阶滑模的所述超螺旋趋近率为:
Figure BDA00029386805400000312
其中,λ3>0,λ4>0,
Figure BDA0002938680540000041
将所述超局部模型得到,n取1,如下,
Figure BDA0002938680540000042
则所述时间延迟观测器为:
Figure BDA0002938680540000043
其中,
Figure BDA0002938680540000044
为扰动G的观测值,ν为延时的时间间隔,用于观测全部位置扰动G。
作为本发明所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的一种优选方案,其中:包括,非线性分数阶滑模面的一阶微分为:
Figure BDA0002938680540000045
输出的跟踪误差一阶微分为:
Figure BDA0002938680540000046
Figure BDA0002938680540000047
得到所述无模型控制器,如下,
Figure BDA0002938680540000048
所述无模型控制器利用所述超螺旋趋近率消除切换增益时的抖振,调用非线性函数和分数阶微积分完成快速收敛。
本发明的有益效果:本发明将超螺旋非线性分数阶滑模的无模型控制方法应用到高超声速飞行器气动热地面模拟系统上,一方面首次提出了石英灯加热器的数学模型表达式;另一方面所设计的控制方法提出了无模型的控制,将全部未知的扰动通过时间延时观测器来观测,将非线性的数学模型进行了线性化的处理;同时非线性分数阶滑模面的设计既保证了控制的稳定性、收敛速度又降低了稳态误差和饱和误差,而超螺旋趋近率的结合改善了滑模面的抖动问题。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。其中:
图1为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的流程示意图;
图2为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的高超声速飞行器气动热地面模拟系统工作流程示意图;
图3(a)为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的高超声速飞行器机翼三维结构示意图;
图3(b)为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的高超声速飞行器机翼二维尺寸示意图;
图4(a)为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的高超声速飞行器机翼气动热的有限元仿真示意图;
图4(b)为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的高超声速飞行器机翼气动热的壁面平均温度采样示意图;
图4(c)为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的高超声速飞行器机翼气动热的数据拟合示意图;
图5为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的无模型原理控制框架示意图;
图6为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的高超声速飞行器气动热地面模拟系统的跟踪拟合目标图(a)以及超螺旋非线性分数阶滑模的无模型控制方法(1)与传统线性分数阶控制方法(2)、iPID控制方法(3)、传统PID控制方法(4)对比的温度跟踪曲线(b)和局部放大图(c);
图7为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的在跟踪拟合目标下,高超声速飞行器气动热地面模拟系统的超螺旋非线性分数阶滑模的无模型控制方法(1)与传统线性分数阶控制方法(2)、iPID控制方法(3)、传统PID控制方法(4)对比的误差跟踪曲线(a)和局部放大图(b,c);
图8为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的在跟踪拟合目标下,高超声速飞行器气动热地面模拟系统的超螺旋非线性分数阶滑模的无模型控制方法(1),在外扰动下(a)与传统线性分数阶控制方法(2)、iPID控制方法(3)、传统PID控制方法(4)对比的温度跟踪曲线(b)和局部放大图(c);
图9为本发明一个实施例所述的气动热的超螺旋非线性分数阶滑模无模型控制方法的在跟踪拟合目标下,高超声速飞行器气动热地面模拟系统的超螺旋非线性分数阶滑模的无模型控制方法(1),在外扰动下与传统线性分数阶控制方法(2)、iPID控制方法(3)、传统PID控制方法(4)对比的误差跟踪曲线(a)和局部放大图(b,c)。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合说明书附图对本发明的具体实施方式做详细的说明,显然所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明的保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
其次,此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例,也不是单独的或选择性的与其他实施例互相排斥的实施例。
本发明结合示意图进行详细描述,在详述本发明实施例时,为便于说明,表示器件结构的剖面图会不依一般比例作局部放大,而且所述示意图只是示例,其在此不应限制本发明保护的范围。此外,在实际制作中应包含长度、宽度及深度的三维空间尺寸。
同时在本发明的描述中,需要说明的是,术语中的“上、下、内和外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一、第二或第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
本发明中除非另有明确的规定和限定,术语“安装、相连、连接”应做广义理解,例如:可以是固定连接、可拆卸连接或一体式连接;同样可以是机械连接、电连接或直接连接,也可以通过中间媒介间接相连,也可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1
参照图1~图5,为本发明的第一个实施例,提供了一种气动热的超螺旋非线性分数阶滑模无模型控制方法,本发明方法基于构建的高超声速飞行器气动热地面模拟系统无模型控制的超局部模型,将非线性分数阶滑模面、超螺旋趋近律、时间延时观测器结合进行控制器u(t)的设计,参照图5,则为本发明高超声速飞行器气动热地面模拟系统的超螺旋非线性分数阶滑模的无模型控制框图,具体实现步骤包括:
S1:基于能量守恒定律,建立高超声速飞行器气动热地面模拟系统输入电能与输出温度之间的数学模型并转化为无模型控制的超局部模型。其中需要说明的是,高超声速飞行器气动热地面模拟系统包括:
非接触辐射加热器、电功率调节装置和量热传感器;
非接触辐射加热器为石英灯加热器,电功率调节装置为双向晶闸管,量热传感器为热电偶传感器。
具体的,根据能量守恒定律构建输入输出能量守恒等式,包括:
Figure BDA0002938680540000071
其中,等式左边UI为输入电压即电源两端电压,R为非接触辐射加热器的电阻之和,α为双向晶闸管的导通角,等式右边分别为用于非接触辐射加热器自身消耗的内能、对流换热过程中损失的热能、热传导过程中损失的热能、热辐射效应输出的热能,c、m、T1、T0、A、ε、Δt分别为非接触辐射加热器的比热容、质量、当前温度、初始温度、表面积、黑度系数、工作时间,β、λ、σ、F分别为对流换热系数、导热系数、斯蒂芬-玻尔兹曼常数、角系数,从而得到当前温度T1和双向晶闸管的导通角α之间的数学关系。
进一步的,还包括:
当被控对象模型是单输入单输出系统时,则将被控对象模型转化为无模型控制的超局部模型,如下,
y(n)=G+χu(t)
其中,y(n)表示为输出量y对时间t的n阶导数,n一般取1或者2,u表示为输入量,G表示为所有未知扰动的集合,既包含了外界扰动和系统内部非线性扰动,χ表示为非物理意义的可调参数;
根据超局部模型将输入输出能量守恒等式两边除以Δt并进行移项处理,获得无模型控制的超局部控制模型,如下,
Figure BDA0002938680540000081
其中,
Figure BDA0002938680540000082
为T1对时间Δt的导数,超局部控制模型中
Figure BDA0002938680540000083
α分别对应超局部模型的y(n)、u,而sin2α给系统带来的是周期性的震动,并没有对系统整体的收敛产生影响,含有sin2α的项可以看作输入扰动,
Figure BDA0002938680540000084
可以看作系统的高阶输出扰动,因此
Figure BDA0002938680540000085
可以看作既包含输入扰动又包含输出扰动的全部扰动之和,对应于超局部模型的G。
S2:根据定义的输出跟踪误差、非线性函数和分数阶微积分,构建非线性分数阶滑模面。本步骤需要说明的是:
定义输出的跟踪误差e表达式,如下,
e(t)=y*-y
其中,e为跟踪误差,y*为输出目标;
非线性函数fal(e,γ,η)为:
Figure BDA0002938680540000091
其中0<γ<1,η>0,
Figure BDA0002938680540000092
非线性分数阶滑模面为:
Figure BDA0002938680540000093
其中,φna>0,φnb>0,φnc>0,用于调参,
Figure BDA0002938680540000094
为分数阶积分,
Figure BDA0002938680540000095
为分数阶微分,
Figure BDA0002938680540000096
ψ为可调参数;
定义非线性分数阶滑模的超螺旋趋近率为:
Figure BDA0002938680540000097
其中,λ3>0,λ4>0,
Figure BDA0002938680540000098
将超局部模型得到,n取1,如下,
Figure BDA0002938680540000099
则时间延迟观测器为:
Figure BDA00029386805400000910
其中,
Figure BDA00029386805400000911
为扰动G的观测值,ν为延时的时间间隔,用于观测全部位置扰动G。
S3:结合非线性分数阶滑模面、超螺旋趋近率、超局部模型及时间延时观测器,搭建超螺旋非线性分数阶滑模的无模型控制器,抑制控制过程中的抖振。其中还需要说明的是,包括:
非线性分数阶滑模面的一阶微分如下,
Figure BDA00029386805400000912
输出的跟踪误差一阶微分为:
Figure BDA00029386805400000913
Figure BDA0002938680540000101
得到无模型控制器,如下,
Figure BDA0002938680540000102
无模型控制器利用超螺旋趋近率消除切换增益时的抖振,调用非线性函数和分数阶微积分完成快速收敛;
构造Lyapunov分析稳定性,即构造Lyapunov稳定性判据,如下,
Figure BDA0002938680540000103
Figure BDA0002938680540000104
Figure BDA0002938680540000111
Figure BDA0002938680540000112
φnaγ|e(t)|γ-1>0、
Figure BDA0002938680540000113
s1(t)sign(s1(t))>0、s1(t)∫sign(s1(t))dt>0
Figure BDA0002938680540000114
只要λ3、λ4足够大,则满足
Figure BDA0002938680540000115
参照图2,整个高超声速飞行器气动热地面模拟系统工作流程主要包括一下步骤:
(1)高超声速飞行器气动热数据采集:对高超声速飞行器的机翼部分,通过给定飞行环境和壁面材料型号对其进行有限元数值模拟,采集每一时刻机翼壁面的平均温度,对采样数据进行线性拟合,作为整个高超声速飞行器气动热地面模拟系统的期望输出值,即目标值,为了和实际控制器输出值进行对比。
(2)高超声速飞行器气动热地面模拟控制系统:设计控制器对石英灯加热系统进行控制,将目标值加载到控制板中,通过控制板改变双向晶闸管的导通角α进而改变输出电压U,不同的输出电压U值对应不同的石英灯加热系统电功率P,通过传感器获得石英灯加热器输出的实际温度T1,与目标值对比得到跟踪误差e,再通过闭环反馈给控制器进行调节双向晶闸管的导通角α,最终达到跟踪控制。
(3)地面模拟试验反馈:石英灯加热器对试件进行加热试验,再通过对试件的性能进行检测,分析材料的可行性,选定材料,若不可以就更换材料再次经过第一步操作,从而优化热防护系统设计。
参照图3,为有限元仿真所绘制的机翼结构,机翼的具体参数为:翼根3550mm,翼展1250mm,前缘后掠角70°,后缘前掠角15°,板厚160mm,前缘半径40mm,材料为镍基高温合金GH1015,飞行环境为高度20km,速度为5.5马赫数,攻角10°巡航。
较佳的是,本发明为了能够保证高超声速飞行器气动热地面模拟系统,能够快速地进行温度跟踪并减小误差、抑制控制过程中的抖动现象,本实施例中提供了一种气动热的超螺旋非线性分数阶滑模无模型控制方法,其中,高超声速飞行器气动热地面模拟系统的无模型控制超局部模型,有效地避免非线性系统全局线性化的问题;在传统的滑模面上结合非线性函数fal(e,γ,η)和分数阶微积分
Figure BDA0002938680540000121
解决了快速收敛大超调的问题;采用超螺旋趋近律解决了切换增益时的抖振问题。
实施例2
参照图4~图9,为本发明的第二个实施例,该实施例不同于第一个实施例的是,提供了一种气动热的超螺旋非线性分数阶滑模无模型控制方法的测试验证,包括:
优选的,为了更好地对本发明方法中采用的技术效果加以验证说明,本实施例中选择以传统的线性分数阶控制方法与本发明方法进行对比测试,以科学论证的手段对比试验结果,验证本发明方法所具有的真实效果。
传统的线性分数阶控制方法不能减小抖动现象、控制精度不高,不能满足被控对象的高性能指标的要求,为验证本发明方法相对于传统方法具有较高的控制精度、抗干扰性、减小抖振,本实施例中将采用高超声速飞行器气动热地面模拟系统的超螺旋非线性分数阶滑模的无模型控制方法(1)在有无外扰动下,与传统线性分数阶控制方法(2)、iPID控制方法(3)、传统PID控制方法(4)分别对高超声速飞行器气动热地面模拟系统的输出温度和跟踪误差进行实时测量对比。
测试环境:参照图4(c),将高超声速飞行器气动热地面模拟系统运行在仿真平台模拟跟踪期望目标曲线,分别利用高超声速飞行器气动热地面模拟系统的超螺旋非线性分数阶滑模的无模型控制方法(1),在有无外扰动下,与传统线性分数阶控制方法(2)、iPID控制方法(3)、传统PID控制方法(4)进行测试并获得测试结果数据,两者都将在开启自动化测试设备并运用MATLB软件编程实现对比方法的仿真测试,根据实验结果得到仿真数据,每种方法各测试4组数据,每组数据采样18s,计算获得每组数据输入温度和跟踪误差,与仿真模拟输入的期望目标温度进行对比计算误差。
参照图4,是机翼进行有限元仿真的仿真模拟图、壁面平均温度采样图和平均温度曲线拟合图,拟合曲线为:
y*=7.224×10-6t6-0.001041×t5+0.05614×t4-1.353×t3+11.86×t2+43.25t+279.2
参照图5,本发明高超声速飞行器气动热地面模拟系统的超螺旋非线性分数阶滑模的无模型控制原理框图是对图2的高超声速飞行器地面模拟控制系统的进一步的说明图5的控制过程分为两个阶段:一、对闭环反馈控制的误差e做超螺旋非线性分数阶滑模控制,分为超螺旋的趋近阶段和非线性分数阶的滑动状态;二、是无模型控制超局部模型中,延时观测器对扰动G的在线观测。
则外扰动为时变电阻R:
R=3.08×(1+0.0045y*)
传统线性分数阶控制方法(2):
Figure BDA0002938680540000131
Figure DA00029386805465237031
iPID控制方法(3):
Figure BDA0002938680540000142
表1:各项参数设置表。
Figure BDA0002938680540000143
表2:本发明方法参数表。
Figure BDA0002938680540000144
表3:传统方法参数表。
Figure BDA0002938680540000145
表4:IPID和PID参数表。
Figure BDA0002938680540000146
参照图6~图9,为本发明在跟踪拟合下,高超声速飞行器气动热地面模拟系统的超螺旋非线性分数阶滑模的无模型控制方法(1),在有无外扰动下,与传统线性分数阶控制方法(2)、iPID控制方法(3)、传统PID控制方法(4)对比的温度跟踪曲线和局部放大图、误差跟踪曲线和局部放大图。
参照图6和图7,能够直观的看出本发明方法(1)的稳定性、控制精度、快速性都比较好,而方法(3)和(4)在图7中,当时间0~0.5s时,具有较大的超调量,方法(2)相比于方法(1)在图7中,出现大量的抖振现象。
参照图8和图9,能够直观的看出方法(2)的抖振明显,方法(3)快速性较差,方法(4)具有较大的超调量,从而说明对比的三种方法相较于方法(1)来说,抗干扰能力差。
应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (8)

1.一种气动热的超螺旋非线性分数阶滑模无模型控制方法,其特征在于:包括,
基于能量守恒定律,建立高超声速飞行器气动热地面模拟系统输入电能与输出温度之间的数学模型并转化为无模型控制的超局部模型;
根据定义的输出跟踪误差、非线性函数和分数阶微积分,构建非线性分数阶滑模面;
结合所述非线性分数阶滑模面、超螺旋趋近率、所述超局部模型及时间延时观测器,搭建超螺旋非线性分数阶滑模的无模型控制器,抑制控制过程中的抖振。
2.根据权利要求1所述的气动热的超螺旋非线性分数阶滑模无模型控制方法,其特征在于:所述高超声速飞行器气动热地面模拟系统包括,非接触辐射加热器、电功率调节装置和量热传感器;
所述非接触辐射加热器为石英灯加热器,所述电功率调节装置为双向晶闸管,所述量热传感器为热电偶传感器。
3.根据权利要求1或2所述的气动热的超螺旋非线性分数阶滑模无模型控制方法,其特征在于:根据所述能量守恒定律构建输入输出能量守恒等式,包括,
Figure FDA0002938680530000011
其中,等式左边UI为输入电压即电源两端电压,R为非接触辐射加热器的电阻之和,α为双向晶闸管的导通角,等式右边分别为用于非接触辐射加热器自身消耗的内能、对流换热过程中损失的热能、热传导过程中损失的热能、热辐射效应输出的热能,c、m、T1、T0、A、ε、Δt分别为非接触辐射加热器的比热容、质量、当前温度、初始温度、表面积、黑度系数、工作时间,β、λ、σ、F分别为对流换热系数、导热系数、斯蒂芬-玻尔兹曼常数、角系数,从而得到当前温度T1和双向晶闸管的导通角α之间的数学关系。
4.根据权利要求3所述的气动热的超螺旋非线性分数阶滑模无模型控制方法,其特征在于:包括,
当被控对象模型是单输入单输出系统时,则将所述被控对象模型转化为无模型控制的超局部模型,如下,
y(n)=G+χu(t)
其中,y(n)表示为输出量y对时间t的n阶导数,n一般取1或者2,u表示为输入量,G表示为所有未知扰动的集合,既包含了外界扰动和系统内部非线性扰动,χ表示为非物理意义的可调参数。
5.根据权利要求4所述的气动热的超螺旋非线性分数阶滑模无模型控制方法,其特征在于:还包括,
根据所述超局部模型将所述输入输出能量守恒等式两边除以Δt并进行移项处理,获得无模型控制的超局部控制模型,如下,
Figure FDA0002938680530000021
其中,
Figure FDA0002938680530000022
为T1对时间Δt的导数,所述超局部控制模型中
Figure FDA0002938680530000023
α分别对应所述超局部模型的y(n)、u,而sin2α给系统带来的是周期性的震动,并没有对系统整体的收敛产生影响,含有sin2α的项可以看作输入扰动,AεσFT1 4可以看作系统的高阶输出扰动,因此
Figure FDA0002938680530000024
可以看作既包含输入扰动又包含输出扰动的全部扰动之和,对应于超局部模型的G。
6.根据权利要求5所述的气动热的超螺旋非线性分数阶滑模无模型控制方法,其特征在于:包括,
定义输出的跟踪误差e表达式为:
e(t)=y*-y
其中,e为跟踪误差,y*为输出目标;
非线性函数fal(e,γ,η)为:
Figure FDA0002938680530000025
其中0<γ<1,η>0,
Figure FDA0002938680530000026
非线性分数阶滑模面为:
Figure FDA0002938680530000031
其中,φna>0,φnb>0,φnc>0,用于调参,
Figure FDA0002938680530000032
为分数阶积分,
Figure FDA0002938680530000033
为分数阶微分,
Figure FDA0002938680530000034
ψ为可调参数。
7.根据权利要求6所述的气动热的超螺旋非线性分数阶滑模无模型控制方法,其特征在于:包括,
定义非线性分数阶滑模的所述超螺旋趋近率为:
Figure FDA0002938680530000035
其中,λ3>0,λ4>0,
Figure FDA0002938680530000036
将所述超局部模型得到,n取1,如下,
Figure FDA0002938680530000037
则所述时间延迟观测器为:
Figure FDA0002938680530000038
其中,
Figure FDA0002938680530000039
为扰动G的观测值,ν为延时的时间间隔,用于观测全部位置扰动G。
8.根据权利要求7所述的气动热的超螺旋非线性分数阶滑模无模型控制方法,其特征在于:包括,
非线性分数阶滑模面的一阶微分为:
Figure FDA00029386805300000310
输出的跟踪误差一阶微分为:
Figure FDA00029386805300000311
Figure FDA00029386805300000312
得到所述无模型控制器,如下,
Figure FDA0002938680530000041
所述无模型控制器利用所述超螺旋趋近率消除切换增益时的抖振,调用非线性函数和分数阶微积分完成快速收敛。
CN202110170226.6A 2021-02-08 2021-02-08 一种气动热的超螺旋非线性分数阶滑模无模型控制方法 Active CN112987566B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110170226.6A CN112987566B (zh) 2021-02-08 2021-02-08 一种气动热的超螺旋非线性分数阶滑模无模型控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110170226.6A CN112987566B (zh) 2021-02-08 2021-02-08 一种气动热的超螺旋非线性分数阶滑模无模型控制方法

Publications (2)

Publication Number Publication Date
CN112987566A true CN112987566A (zh) 2021-06-18
CN112987566B CN112987566B (zh) 2021-11-02

Family

ID=76349182

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110170226.6A Active CN112987566B (zh) 2021-02-08 2021-02-08 一种气动热的超螺旋非线性分数阶滑模无模型控制方法

Country Status (1)

Country Link
CN (1) CN112987566B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113467245A (zh) * 2021-07-15 2021-10-01 北京信息科技大学 飞行器的分数阶滑模控制方法、装置及系统
CN113625564A (zh) * 2021-08-06 2021-11-09 南京工业大学 一种结构热试验基于系统模型的滑模控制方法
CN114355779A (zh) * 2022-01-10 2022-04-15 南京工业大学 一种结构热试验非线性扩展状态观测器全局滑模控制方法
CN114371623A (zh) * 2022-01-10 2022-04-19 南京工业大学 一种结构热试验神经网络观测器分数阶全局滑模控制方法
CN114721274A (zh) * 2022-04-24 2022-07-08 广东工业大学 一种基于改进fal函数的滑模控制器设计方法及系统
CN114967462A (zh) * 2022-06-01 2022-08-30 南京工业大学 一种塔式起重机输出抖动抑制的滑模控制方法
CN115586724A (zh) * 2022-10-27 2023-01-10 南京师范大学泰州学院 一种齿轮巡检机器人系统自适应分数阶全局滑模控制方法
CN116317775A (zh) * 2023-03-17 2023-06-23 铜陵学院 一种新型电动助力车驱动系统双分数阶积分滑模控制方法
CN117312778A (zh) * 2023-11-28 2023-12-29 中国航空工业集团公司沈阳空气动力研究所 风洞气动热试验数据滤波处理方法、电子设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103439133A (zh) * 2013-08-27 2013-12-11 北京航空航天大学 高超声速飞行器250℃/秒高速热冲击试验装置及方法
CN105548250A (zh) * 2016-01-13 2016-05-04 北京机电工程研究所 飞行器气动热试验的热流控制方法、装置及系统
US20180131021A1 (en) * 2016-03-16 2018-05-10 Huazhong University Of Science And Technology Thermoelectric cooperative control method for sofc system based on fractional order sliding mode variable structure
CN110635734A (zh) * 2019-09-24 2019-12-31 南京工业大学 一种新型永磁同步电机转速分数阶滑模控制方法
WO2020124938A1 (zh) * 2018-12-21 2020-06-25 南京埃斯顿机器人工程有限公司 柔性关节机械臂分数阶滑模优化控制方法
CN112747890A (zh) * 2020-12-29 2021-05-04 中国航天空气动力技术研究院 一种气动热力联合试验系统及试验方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103439133A (zh) * 2013-08-27 2013-12-11 北京航空航天大学 高超声速飞行器250℃/秒高速热冲击试验装置及方法
CN105548250A (zh) * 2016-01-13 2016-05-04 北京机电工程研究所 飞行器气动热试验的热流控制方法、装置及系统
US20180131021A1 (en) * 2016-03-16 2018-05-10 Huazhong University Of Science And Technology Thermoelectric cooperative control method for sofc system based on fractional order sliding mode variable structure
WO2020124938A1 (zh) * 2018-12-21 2020-06-25 南京埃斯顿机器人工程有限公司 柔性关节机械臂分数阶滑模优化控制方法
CN110635734A (zh) * 2019-09-24 2019-12-31 南京工业大学 一种新型永磁同步电机转速分数阶滑模控制方法
CN112747890A (zh) * 2020-12-29 2021-05-04 中国航天空气动力技术研究院 一种气动热力联合试验系统及试验方法

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
ALI SAFAEI,MUHAMMAD NASIRUDDIN MAHYUDDIN: "Adaptive Model-Free Control Based on an Ultra-Local Model With Model-Free Parameter Estimations for a Generic SISO System", 《IEEE ACCESS》 *
G. VENU; S. TARA KALYANI: "Design of Fractional Order based Super Twisting Algorithm for BLDC motor", 《2019 3RD INTERNATIONAL CONFERENCE ON TRENDS IN ELECTRONICS AND INFORMATICS (ICOEI)》 *
GAO,P;LV,XD;OUYANG,HM;MEI,L;ZHANG,GM: "A Novel Model-Free Intelligent Proportional-Integral Supertwisting Nonlinear Fractional-Order Sliding Mode Control of PMSM Speed Regulation System", 《COMPLEXITY》 *
PENG GAO,GUANGMING ZHANG,HUIMIN OUYANG,LEI MEI: "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》 *
WANGH.P.,GHAZALLY I.Y.MUSTAFA,Y.TIAN: "Model-free fractional-order sliding mode control for an active vehicle suspension system", 《ADVANCES IN ENGINEERING SOFTWARE》 *
兰天一,林辉,张克军: "石英灯辐射式气动热实验的分数阶迭代学习控制策略", 《西北工业大学学报》 *
吴大方,潘兵,高镇同,牟朦,朱林,王岳武: "超高温、大热流、非线性气动热环境试验模拟及测试技术研究", 《实验力学》 *
吴大方等: "模糊控制在导弹瞬态气动热流试验中的研究", 《北京航空航天大学学报》 *
吴大方等: "高超声速飞行器材料与结构气动热环境模拟方法及试验研究", 《航天器环境工程》 *
赵凯辉,殷童欢,张昌凡,李祥飞,陈跃,何静: "永磁同步电机无模型滑模控制方法研究", 《电子测量与仪器学报》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113467245A (zh) * 2021-07-15 2021-10-01 北京信息科技大学 飞行器的分数阶滑模控制方法、装置及系统
CN113467245B (zh) * 2021-07-15 2023-06-02 北京信息科技大学 飞行器的分数阶滑模控制方法、装置及系统
CN113625564B (zh) * 2021-08-06 2023-11-24 南京工业大学 一种结构热试验基于系统模型的滑模控制方法
CN113625564A (zh) * 2021-08-06 2021-11-09 南京工业大学 一种结构热试验基于系统模型的滑模控制方法
CN114355779A (zh) * 2022-01-10 2022-04-15 南京工业大学 一种结构热试验非线性扩展状态观测器全局滑模控制方法
CN114371623A (zh) * 2022-01-10 2022-04-19 南京工业大学 一种结构热试验神经网络观测器分数阶全局滑模控制方法
CN114355779B (zh) * 2022-01-10 2022-07-19 南京工业大学 一种结构热试验非线性扩展状态观测器全局滑模控制方法
CN114721274A (zh) * 2022-04-24 2022-07-08 广东工业大学 一种基于改进fal函数的滑模控制器设计方法及系统
CN114721274B (zh) * 2022-04-24 2022-11-01 广东工业大学 一种基于改进fal函数的滑模控制器设计方法及系统
CN114967462A (zh) * 2022-06-01 2022-08-30 南京工业大学 一种塔式起重机输出抖动抑制的滑模控制方法
CN114967462B (zh) * 2022-06-01 2023-05-30 南京工业大学 一种塔式起重机输出抖动抑制的滑模控制方法
CN115586724A (zh) * 2022-10-27 2023-01-10 南京师范大学泰州学院 一种齿轮巡检机器人系统自适应分数阶全局滑模控制方法
CN115586724B (zh) * 2022-10-27 2023-11-24 南京师范大学泰州学院 一种齿轮巡检机器人系统自适应分数阶全局滑模控制方法
CN116317775A (zh) * 2023-03-17 2023-06-23 铜陵学院 一种新型电动助力车驱动系统双分数阶积分滑模控制方法
CN116317775B (zh) * 2023-03-17 2023-12-08 铜陵学院 一种电动助力车驱动系统双分数阶积分滑模控制方法
CN117312778A (zh) * 2023-11-28 2023-12-29 中国航空工业集团公司沈阳空气动力研究所 风洞气动热试验数据滤波处理方法、电子设备及存储介质
CN117312778B (zh) * 2023-11-28 2024-01-30 中国航空工业集团公司沈阳空气动力研究所 风洞气动热试验数据滤波处理方法、电子设备及存储介质

Also Published As

Publication number Publication date
CN112987566B (zh) 2021-11-02

Similar Documents

Publication Publication Date Title
CN112987566B (zh) 一种气动热的超螺旋非线性分数阶滑模无模型控制方法
CN109883660B (zh) 一种热模拟试验控制方法
CN112965384B (zh) 一种自适应神经网络全局非奇异快速终端滑模无模型方法
CN112965382B (zh) 一种基于神经网络观测器非线性全局滑模无模型控制方法
CN112965364B (zh) 一种高超声速飞行器气动热地面模拟系统建模及iPID方法
CN112987569B (zh) 一种扰动上界自适应分数阶全局终端滑模无模型控制方法
CN113625564B (zh) 一种结构热试验基于系统模型的滑模控制方法
CN106528990B (zh) 一种基于泛函优化的高超声速尖锥外形热流密度建模方法
CN113495486A (zh) 一种结构热试验基于扩展状态观测器的模型预测控制方法
CN110096755B (zh) 固体蓄热炉内高温加热元件在线温度软测量方法及系统
CN112859617B (zh) 一种iPI无模型自适应全局非奇异快速终端滑模控制方法
CN114637277B (zh) 空天飞机测试实验的全方程热流密度测控系统及测控方法
CN107272412A (zh) 一种暂冲式风洞流场控制的辩识方法
CN115236975B (zh) 一种导弹气动热地面有限元分析及ipd非线性滑模控制方法
CN107092724A (zh) 一种考虑不确定性的热防护系统模型修正方法
CN113900373A (zh) 一种结构热试验的全局非奇异快速终端滑模控制方法
CN114371623B (zh) 一种结构热试验神经网络观测器分数阶全局滑模控制方法
CN114721450B (zh) 空天飞机结构强度测试中温度控制用热试验迭代控温方法
Ferrero et al. A numerical method for conjugate heat transfer problems in hypersonic flows
NOLL Aeroservoelasticity
CN208155557U (zh) 基于手自动双模式功率匹配回流风洞加热系统
Mouton Numerical Simulation of the Flow in the ONERA F1 Wind Tunnel
Poussot-Vassal et al. Aircraft feedback gust load control design and experimental wing bench validation
CN103245481A (zh) 基于变频技术的变负荷大型换热器气阻特性的检测方法
Chu et al. Dynamic Compensation Method for Temperature Measurement Error of High Temperature and High-speed Airflow

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