CN112363411A - 一种航空发动机动态矩阵控制器的设计方法 - Google Patents
一种航空发动机动态矩阵控制器的设计方法 Download PDFInfo
- Publication number
- CN112363411A CN112363411A CN202011282663.9A CN202011282663A CN112363411A CN 112363411 A CN112363411 A CN 112363411A CN 202011282663 A CN202011282663 A CN 202011282663A CN 112363411 A CN112363411 A CN 112363411A
- Authority
- CN
- China
- Prior art keywords
- prediction
- control
- output
- time
- model
- 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.)
- Withdrawn
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000013461 design Methods 0.000 title claims abstract description 14
- 238000005457 optimization Methods 0.000 claims description 14
- 230000004044 response Effects 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 9
- 238000005096 rolling process Methods 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 7
- 239000000446 fuel Substances 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 11
- 238000004422 calculation algorithm Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 3
- 239000000295 fuel oil Substances 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000007704 transition Effects 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
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
- F02C9/00—Controlling gas-turbine plants; Controlling fuel supply in air- breathing jet-propulsion plants
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明属于航空宇航推进理论与工程中的系统控制与仿真技术领域,提供了一种航空发动机动态矩阵控制器的设计方法,算法的特点是离线求解优化问题,将求解得到每个时刻的解序列存储在矩阵中,同时存储矩阵采用稀疏矩阵压缩存储方法,每个输出只取未来时刻的一个点进行优化,对每个输入只计算出一步控制作用,大幅度减少了计算量,为硬件在环和半实物仿真试验提供了理论基础。数值仿真表明,航空发动机可实现从起动到慢车、中间等状态的稳态和过渡态控制,同时满足执行机构和输出约束,设计的动态矩阵控制器具有一定的工程应用价值。
Description
技术领域
本发明提供了一种航空发动动态矩阵控制器的设计方法,属于航空宇航推进理论与工程中的系统控制与仿真领域。
背景技术
模型预测控制作为一种先进的控制策略在石油加工、电力、化工等辅助工业过程中得到广泛的应用,获得了很大的经济效益。多变量预测控制的基本目标是在满足各种约束条件的基础上,将被控变量尽可能在其设定值或容许值范围内。对于多变量系统,传统的模型预测控制算法是采用预测时域上所有预测值与设定值或参考轨迹偏差的平方和最小为优化目标函数来确定未来的控制作用,在每一采样时刻,需要求解维数较大的二次优化问题,因此模型预测控制执行周期通常为几分钟,适合处理时间响应较慢的过程,而在处理航空发动机控制和电机控制等快速动态问题时无法保证实时性。本发明提出一种实时性满足要求的动态矩阵控制算法,算法的特点是离线求解优化问题,将求解得到每个时刻的解序列存储在矩阵中,同时存储矩阵采用稀疏矩阵压缩存储方法,每个输出只取未来时刻的一个点进行优化,对每个输入只计算出一步控制作用,大幅度减少了计算量,为硬件在环和半实物仿真试验提供了理论基础。目前为止,没有专利公开基于航空发动机的动态矩阵控制器的设计方法。
发明内容
为了解决传统模型预测控制算法在实际应用中存在运算量大、实时性较差的问题,本发明提出了一种基于离线优化的动态矩阵控制器设计方法。
一种航空发动机动态矩阵控制器设计方法,包括以下步骤:
第一步,获得航空发动机辨识模型;
第二步,设计动态矩阵控制器;
第三步,数值仿真验证。
现具体阐述如下:
第一步,获得航空发动机辨识模型;
由航空发动机原理知,转速为常数的调节方案能保持发动机的机械载荷、热负荷状态近似不变,同时保证推力稳定;压比为常数和发动机推力和耗油率直接相关;本发明通过matlab自带的工具箱ident辨识得到航空发动机在标准状态点的线性模型,因此选择主燃油和压比为控制变量,即辨识输入数据,被控变量为高压转子转速和涡轮落压比,即辨识输出数据。
第二步,设计动态矩阵控制器;
动态矩阵控制器包括预测模型、滚动优化、反馈校正三个部分。传统的动态矩阵控制采用非线性模型获得预测模型,本发明采用辨识模型作为被控模型,离线获取响应数据,从而得到预测模型;滚动优化是在每一个采样时间点,算法只优化一个范围内的有限时域,得出所有优化控制序列后,只实施当前控制量给发动机,到下一采样时间点又再次进行优化预测;反馈校正是根据预测模型计算出下一采样时间点的输出,到下一采样时间点用实测输出与预测输出求误差,以此对未来的输出进行启发式修正。动态矩阵控制器通过将得到的二次规划问题离线求解,并将解序列存储在矩阵中,在线计算时只需要从矩阵中寻找即可,大大减少了计算量。
第三步,数值仿真验证。
为了验证设计的动态矩阵控制器的性能,基于某型军用小涵道比涡扇发动机非线性模型,在不同的工况点,手动设置油门杆角度使航空发动机在起动到慢车、中间状态完成稳态和过渡态控制,同时设置人为噪声和干扰验证控制器的鲁棒性。
本发明的有益效果:
传统模型预测控制算法需要在每个控制周期在线求解一个二次规划问题,无法满足航空发动机硬件在环仿真试验实时性的需求。本发明采用一种基于离线优化的动态矩阵控制器设计策略,大幅度减少了控制器计算量,为模型预测控制的嵌入式应用提供了参考。数值仿真表明,动态矩阵控制器在完成航空发动机稳态和过渡态控制的同时,还具有较强的鲁棒性和约束管理能力。本发明设计流程简单,可操作性强,为现代控制理论在FADEC上的应用提供了理论参考。
附图说明
图1是多变量控制架构图。
图2给定的油门杆角度图。
图3是航空发动机在H=0km,Ma=0工况下主燃油曲线。
图4是是航空发动机在H=0km,Ma=0工况下尾喷管面积曲线。
图5是航空发动机在H=0km,Ma=0工况下高压转子响应曲线。
图6是航空发动机在H=0km,Ma=0工况下压比响应曲线。
图7是航空发动机在H=5km,Ma=0.3工况下燃油曲线。
图8是是航空发动机在H=5km,Ma=0.3工况下尾喷管面积曲线。
图9是航空发动机在H=5km,Ma=0.3工况下高压转子响应曲线。
图10是航空发动机在H=5km,Ma=0.3工况下压比响应曲线。
具体实施方式
为了使本发明的目的、技术及优点更加清楚明白,以下结合附图及实例,对本发明进行进一步详细说明。
本发明是一种航空发动机动态矩阵控制器设计方法,具体设计步骤如下:
第一步,获得航空发动机辨识模型;
首先使发动机开环启动到慢车稳态点,分别给主燃油和尾喷管面积一个小扰动,获取主燃油、尾喷管面积、转速、压比的增量输入输出数据,利用ident工具箱即可航空发动机的辨识模型。
第二步,设计动态矩阵控制器;
步骤一,建立预测模型;
动态矩阵控制基于被控对象的阶跃响应来建立非参数模型,其适用于渐进稳定的线性对象。对被控对象的单位采样值:
ai=a(iT),i=1,2,…,N (1)
式中T为采样周期,N为建模时域,在时间t=NT之后,阶跃响应将趋于一个稳定值(稳态),即对象的动态信息可用一个有限区间:
a1,a2,a3,…,aN (2)
来表示,这个集合的参数就构成了DMC的模型参数,向量(2)为矩阵的模型向量。利用(2)就可以对未来时刻的输出值进行预测。在t时刻假设对未来N个时刻的初始预测值为:
当k时刻有一控制量Δu(t)时,可利用线性叠加原理算出其在未来预测时刻内的输出值:
同理,从k时刻开始,在M个连续的控制增量
Δu(k),…,Δu(k+M-1) (5)作用下可以得到未来各个时刻的输出值:
式(5)就是DMC所需要的预测模型。
步骤二,滚动优化;
滚动优化是在每一时刻t,都要再一次确定该时刻之后M个控制增量(5),使得预测输出(4)都能够不断的接近给定的目标值:
yr(t+i),i=1,…,P (7)
M,P分别代表控制步长和预测步长,一般来说M≤P≤N,保证可行性。
为了滚动计算每一时刻t的预测输出,必须确定最优的控制增量Δu(t),设计性能如下的性能指标函数:
qi,rj为权重系数,一般情况设定qi为1,这样可以减少一个参数变量,让控制器更易于调整。将式(8)写为向量的形式为:
其中,
y(t)=[yr(t+1),…,yr(t+P)],
Q=diag(q1,…,qp),R=diag(r1,…,rM)
Q,R分别为误差权重矩阵和控制权重矩阵,且都为对角矩阵。
将(6)代入(9)得:
其中A由模型向量(2)组成,表示为:
当前控制变量为(12)的第一个元素:
求得Δu(t)之后,对其进行积分,得到:
将Δu(t)作为控制量给到发动机,就可以推导出t+1时刻的控制量Δu(t+1),实现反复滚动优化控制。
步骤二,反馈校正。
在t时刻把u(t)施加到预测模型(4)上就可以得到下一时刻的预测输出此时可以将预测输出作为下一时间点t+1时刻的初始预测输出但是由于建模误差和外部扰动等原因,预测值会出现较大的偏差,若不及时进行反馈校正,误差会随着时间积累,甚至导致系统不可控,那我们需要对预测值进行实时校正。在下一时刻t+1,可以得到发动机的真实输出y(t+1)和预测输出可得预测误差:
利用反馈误差对预测输入进行校正就可以得到校正预测值:
为了进一步说明本发明中显式模型预测控制器的控制效果,通过几组仿真实验,来验证本发明中方法的有效性。
(1)数值仿真验证。
本发明选择H=0km,Ma=0和H=5km,Ma=0.3两个包线点,根据航空发动机控制理论,设计如说明书附图1所示的多变量控制结构,仿真结果如图说明书附图2-10所示,图中实线为目标值,虚线为发动机模型被控输出值。由图可知,航空发动机可实现从起动到慢车、中间等状态的稳态和过渡态控制,同时满足执行机构和输出约束。高压转子转速超调不超过±0.5%,稳态误差不超过±1%,落压比稳态控制精度不超过±1%。在图9-10中,转速和压比在起动阶段出现执行机构饱和现象,这是因为高度马赫数变化时,发动机模型非线性变强,但是整体上DMC表现了优异的鲁棒性,说明了设计的算法的可行性。
Claims (1)
1.一种航空发动机动态矩阵控制器的设计方法,其特征在于,步骤如下:
第一步,获得航空发动机辨识模型
首先使发动机开环启动到慢车稳态点,分别给主燃油和尾喷管面积一个小扰动,获取主燃油、尾喷管面积、转速、压比的增量输入输出数据,利用ident工具箱即可航空发动机的辨识模型;
第二步,设计动态矩阵控制器
步骤一,建立预测模型
动态矩阵控制基于被控对象的阶跃响应来建立非参数模型,其适用于渐进稳定的线性对象;对被控对象的单位采样值:
ai=a(iT),i=1,2,…,N (1)
式中,T为采样周期,N为建模时域,在时间t=NT之后,阶跃响应将趋于一个稳定值,即对象的动态信息可用一个有限区间:
a1,a2,a3,…,aN (2)
来表示,这个集合的参数构成DMC的模型参数,公式(2)为矩阵的模型向量;公式(2)对未来时刻的输出值进行预测;在t时刻假设对未来N个时刻的初始预测值为:
当k时刻有一控制量Δu(t)时,可利用线性叠加原理算出其在未来预测时刻内的输出值:
同理,从k时刻开始,在M个连续的控制增量
Δu(k),…,Δu(k+M-1) (5)
作用下得到未来各个时刻的输出值:
式(5)就是DMC所需要的预测模型;
步骤二,滚动优化
滚动优化是在每一时刻t,都要再一次确定该时刻之后M个控制增量(5),使得预测输出(4)都能够不断的接近给定的目标值:
yr(t+i),i=1,…,P (7)
其中,M,P分别代表控制步长和预测步长,M≤P≤N,保证可行性;
为了滚动计算每一时刻t的预测输出,必须确定最优的控制增量Δu(t),设计性能如下的性能指标函数:
其中,qi,rj为权重系数,设定qi为1;将式(8)写为向量的形式为:
其中,
y(t)=[yr(t+1),…,yr(t+P)],
Q=diag(q1,…,qp),R=diag(r1,…,rM)
其中,Q,R分别为误差权重矩阵和控制权重矩阵,且都为对角矩阵;
将(6)代入(9)得:
其中A由模型向量(2)组成,表示为:
当前控制变量为(12)的第一个元素:
求得Δu(t)之后,对其进行积分,得到:
将Δu(t)作为控制量给到发动机,推导出t+1时刻的控制量Δu(t+1),实现反复滚动优化控制;
步骤二,反馈校正
利用反馈误差对预测输入进行校正就得到校正预测值:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011282663.9A CN112363411A (zh) | 2020-11-14 | 2020-11-14 | 一种航空发动机动态矩阵控制器的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011282663.9A CN112363411A (zh) | 2020-11-14 | 2020-11-14 | 一种航空发动机动态矩阵控制器的设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112363411A true CN112363411A (zh) | 2021-02-12 |
Family
ID=74515228
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011282663.9A Withdrawn CN112363411A (zh) | 2020-11-14 | 2020-11-14 | 一种航空发动机动态矩阵控制器的设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112363411A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327669A (zh) * | 2020-11-14 | 2021-02-05 | 大连理工大学 | 一种航空发动机显式预测控制器的设计方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106444388A (zh) * | 2016-12-06 | 2017-02-22 | 杭州电子科技大学 | 一种焦炭炉炉膛压力的分布式pid型动态矩阵控制方法 |
CN106569476A (zh) * | 2016-10-20 | 2017-04-19 | 西安坤蓝电子技术有限公司 | 一种天线轴系固有频率及控制系统带宽的自我检测方法 |
CN106886151A (zh) * | 2017-04-17 | 2017-06-23 | 大连理工大学 | 一种航空发动机多工况下约束预测控制器的设计及调度方法 |
CN108694294A (zh) * | 2018-06-20 | 2018-10-23 | 西安交通大学 | 一种航空发动机高压转子装配偏心量预测方法 |
CN108762089A (zh) * | 2018-06-15 | 2018-11-06 | 大连理工大学 | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 |
CN108803336A (zh) * | 2018-06-28 | 2018-11-13 | 南京航空航天大学 | 一种航空发动机自适应lqg/ltr控制器设计方法 |
CN108958027A (zh) * | 2018-05-21 | 2018-12-07 | 东北大学 | 一种区间系统的动态矩阵控制方法 |
CN111177930A (zh) * | 2019-12-31 | 2020-05-19 | 厦门大学 | 基于混合序列的航空发动机系统辨识建模方法 |
CN111274692A (zh) * | 2020-01-16 | 2020-06-12 | 西安交通大学 | 航空发动机非线性控制系统建模方法 |
CN111856933A (zh) * | 2020-07-06 | 2020-10-30 | 大连理工大学 | 一种考虑不确定性的航空发动机自适应控制器设计方法 |
-
2020
- 2020-11-14 CN CN202011282663.9A patent/CN112363411A/zh not_active Withdrawn
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106569476A (zh) * | 2016-10-20 | 2017-04-19 | 西安坤蓝电子技术有限公司 | 一种天线轴系固有频率及控制系统带宽的自我检测方法 |
CN106444388A (zh) * | 2016-12-06 | 2017-02-22 | 杭州电子科技大学 | 一种焦炭炉炉膛压力的分布式pid型动态矩阵控制方法 |
CN106886151A (zh) * | 2017-04-17 | 2017-06-23 | 大连理工大学 | 一种航空发动机多工况下约束预测控制器的设计及调度方法 |
CN108958027A (zh) * | 2018-05-21 | 2018-12-07 | 东北大学 | 一种区间系统的动态矩阵控制方法 |
CN108762089A (zh) * | 2018-06-15 | 2018-11-06 | 大连理工大学 | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 |
CN108694294A (zh) * | 2018-06-20 | 2018-10-23 | 西安交通大学 | 一种航空发动机高压转子装配偏心量预测方法 |
CN108803336A (zh) * | 2018-06-28 | 2018-11-13 | 南京航空航天大学 | 一种航空发动机自适应lqg/ltr控制器设计方法 |
CN111177930A (zh) * | 2019-12-31 | 2020-05-19 | 厦门大学 | 基于混合序列的航空发动机系统辨识建模方法 |
CN111274692A (zh) * | 2020-01-16 | 2020-06-12 | 西安交通大学 | 航空发动机非线性控制系统建模方法 |
CN111856933A (zh) * | 2020-07-06 | 2020-10-30 | 大连理工大学 | 一种考虑不确定性的航空发动机自适应控制器设计方法 |
Non-Patent Citations (5)
Title |
---|
SI-XIN WEN等: "A Multivariable Robust Adaptive Control Scheme for Aero-engines", 《2020 39TH CHINESE CONTROL CONFERENCE (CCC)》 * |
XIN-YUE WANG等: "Controller Design of Aero-engines under the Distributed Architecture with Time Delays", 《2020 39TH CHINESE CONTROL CONFERENCE (CCC)》 * |
商国军: "基于非线性模型的航空发动机控制技术研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
汪大卫: "分数阶系统的动态矩阵控制算法研究", 《中国优秀硕士论文全文数据库 信息科技辑》 * |
郑斐华: "基于系统辨识的航空发动机建模研究", 《中国优秀硕士论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327669A (zh) * | 2020-11-14 | 2021-02-05 | 大连理工大学 | 一种航空发动机显式预测控制器的设计方法 |
CN112327669B (zh) * | 2020-11-14 | 2022-02-18 | 大连理工大学 | 一种航空发动机显式预测控制器的设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11454177B2 (en) | Method of aero-engine on-line optimization and multivariable control based on model prediction | |
CN108762089B (zh) | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 | |
US11392094B2 (en) | Steady state and transition state multi-variable control method of a turbofan engine based on an active disturbance rejection control theory (ADRC) | |
CN106647253B (zh) | 航空发动机分布式控制系统多性能鲁棒跟踪控制方法 | |
CN112286047B (zh) | 基于神经网络的narma-l2多变量控制方法 | |
CN110579962B (zh) | 基于神经网络的涡扇发动机推力预测方法及控制器 | |
CN112483261B (zh) | 一种航空发动机抗加力扰动方法 | |
DeCastro | Rate-based model predictive control of turbofan engine clearance | |
CN102411305A (zh) | 单旋翼直升机/涡轴发动机综合抗扰控制系统设计方法 | |
CN110219736A (zh) | 基于非线性模型预测控制的航空发动机直接推力控制方法 | |
CN109446605A (zh) | 涡轴发动机非线性动态逆控制方法及装置 | |
CN111880403A (zh) | 航空发动机最大推力状态容错二自由度μ控制器 | |
Liu et al. | Resilient H∞ control for uncertain turbofan linear switched systems with hybrid switching mechanism and disturbance observer | |
CN112363411A (zh) | 一种航空发动机动态矩阵控制器的设计方法 | |
CN111965981B (zh) | 一种航空发动机强化学习控制方法及系统 | |
Chen et al. | Model-Based Improved Advanced Adaptive Performance Recovery Control Method for a Commercial Turbofan Engine | |
CN115981156A (zh) | 一种时变输出约束下的航空发动机主动限制保护控制方法 | |
Wen et al. | A multivariable robust adaptive control scheme for aero-engines | |
CN110985216A (zh) | 一种含在线修正的航空发动机智能多变量控制方法 | |
CN114047692B (zh) | 一种涡扇发动机鲁棒容错抗干扰模型参考动态输出反馈控制方法 | |
CN114415506B (zh) | 航空发动机双模跟踪预测控制系统设计方法 | |
CN111856929B (zh) | 航空发动机容错增益调度二自由度h∞控制器 | |
CN111459028B (zh) | 航空发动机最大推力状态降保守性二自由度μ控制器 | |
CN114035429A (zh) | 一种基于干扰观测器的涡扇发动机切换系统的输出跟踪控制方法 | |
CN112327669B (zh) | 一种航空发动机显式预测控制器的设计方法 |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20210212 |
|
WW01 | Invention patent application withdrawn after publication |