CN112016158B - 一种航空发动机高空小表速抖振抑制方法 - Google Patents
一种航空发动机高空小表速抖振抑制方法 Download PDFInfo
- Publication number
- CN112016158B CN112016158B CN202010697863.4A CN202010697863A CN112016158B CN 112016158 B CN112016158 B CN 112016158B CN 202010697863 A CN202010697863 A CN 202010697863A CN 112016158 B CN112016158 B CN 112016158B
- Authority
- CN
- China
- Prior art keywords
- engine
- aeroengine
- speed
- closed
- interference
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000005764 inhibitory process Effects 0.000 title description 2
- 238000012546 transfer Methods 0.000 claims abstract description 20
- 230000003321 amplification Effects 0.000 claims abstract description 8
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 8
- 239000000295 fuel oil Substances 0.000 claims abstract description 7
- 238000013461 design Methods 0.000 claims description 17
- 230000000694 effects Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000005316 response function Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000000452 restraining effect Effects 0.000 abstract description 5
- 230000002401 inhibitory effect Effects 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 238000013016 damping Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mechanical Engineering (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种航空发动机高空小表速抖振抑制方法,首先,获得航空发动机高空小表速燃油与转速的开环和闭环传递函数;其次,定义闭环控制对航空发动机干扰振动的放大函数;设计最优闭环控制增益以达到抗干扰性能指标;然后,检验最优闭环控制增益满足发动机动态和稳态性能要求;最后,确认最优闭环控制增益。本发明可消除或减轻航空发动机在高空低速飞行包线中的振动量,达到抑制高空小表速抖振的目的。
Description
技术领域
本发明属于航空发动机控制领域,具体涉及一种航空发动机高空小表速抖振抑制方法。
背景技术
作为高安全性和高可靠性的航空发动机系统,适航规章对其整机振动做出了严苛的规定,要求在整个飞行包线、整个转速和推力范围内,不应因过量振动致零部件应力过大或将过大振动力传递至飞行器。因此,必须对航空发动机的振动加以控制。然而,引起航空发动机振动的原因很多,转轴不对中、装配误差、轴承损伤、燃烧室激振、压气机喘振等均能引起发动机的过量振动。在工程实践中,通常是通过减小装配误差来减轻过量振动致零部件应力过大;而通过在传递路径上加装阻尼器来减小将过大振动力传递至飞行器。例如专利CN109489979A公开了一种发动机振动与机匣类零组件相关性的分析方法,其目的即是通过确定机匣类零组件的关键影响因子控制值,以达到降低发动机翻修成本和提升装配质量的目标;而专利CN110552995A则公开了一种航空发动机减振拉杆,通过消耗发动机辅助吊挂处的振动能量来隔离发动机和飞机机体。
在实际工程中,还经常遇到另一种振动问题:飞机在高空低速飞行中,发动机产生剧烈振动,甚至引发控制系统发出空中停车指令。这被称为“高空小表速抖振”的发动机振动问题极其严重但又非常特殊,主要体现在:(1)高空空气稀薄,点火器二次点火困难,使得空中停车成为引发飞行事故甚至坠毁的主要原因之一;(2)高空小表速抖振在实际工程中虽然经常出现,但是至今仍然没有完全明白其原因;(3)高空小表速表现在飞行高度高、马赫数低,相应地发动机进口的空气流量很小、空气温度和压力很低(压力5-10kPa、温度-40°、流量3-5kg/s),超出了高空台设计工作范围,因此难以通过地面高空台试验对其机理进行探究。注意到飞行过程中,大表速不抖振,但是小表速抖振,说明不是装配、破损等结构性因素,很可能是航空发动机作为一个动力系统本身的一个特性,因此需要从动力系统抗干扰特性的角度进行分析。本发明即从动力系统控制角度,提出一种航空发动机高空小表速抖振抑制方法。
发明内容
发明目的:针对航空发动机高空小表速抖振问题,本发明提出一种航空发动机高空小表速抖振抑制方法,可消除或减轻航空发动机在高空低速飞行包线中的振动量,达到抑制高空小表速抖振的目的。
技术方案:本发明所述的一种航空发动机高空小表速抖振抑制方法,包括以下步骤:
(1)获得航空发动机高空小表速燃油与转速的开环和闭环传递函数;
(2)定义闭环控制对航空发动机干扰振动的放大函数;
(3)设计最优闭环控制增益以达到抗干扰性能指标;
(4)检验最优闭环控制增益是否满足发动机动态和稳态性能要求;
(5)如果控制增益同时满足步骤(3)和步骤(4)提出的抗干扰性能指标和稳态性能指标,则确认该最优设计;否则返回步骤(3)和(4)。
进一步地,所述步骤(1)包括以下步骤:
(11)获得航空发动机高空小表速燃油与转速的发动机开环传递函数:
其中,GE(s)为发动机开环传递函数;s为拉普拉斯复变量;a、b、c和d为实系数;
(12)航空发动机在高空、低速飞行时,全权限电子控制器采用比例控制,闭环传递函数表示为:
其中,G(s)为发动机闭环传递函数,k为全权限电子控制器增益。
进一步地,所述步骤(2)的实现过程如下:
衡量闭环控制对航空发动机干扰振动的作用,定义干扰振动放大函数为Td(jω),表示为:
其中,GE(jω)和G(jω)分别为传递函数GE(s)和G(s)的频域响应函数;由于航空发动机为一惯性系统,因此可得:
对某一频率ω,Td(jω)为一复数,其幅值即表达了闭环控制对该三角函数干扰情况下的放大作用:
进一步地,所述步骤(3)的实现过程如下:
限制Td(jω)的幅值,即:
|Td(jω)|≤σ
其中,σ为一正实数,要求σ≤1,满足要求的全权限电子控制器增益k为:
进一步地,所述步骤(4)的实现过程如下:
将稳态性能定义为对某一阶跃转速要求的稳态误差ess,由拉普拉斯变换特性的终值定理可得:
稳态性能指标定义为:
|ess|≤τ
其中,τ为一正实数;即要求全权限电子控制器增益k的设计满足:
为同时满足稳态特性和步骤(3)要求的抗干扰特性,全权限电子控制器增益k的设计必须同时满足和
有益效果:与现有技术相比,本发明的有益效果:本发明针对航空发动机高空小表速抖振问题,通过所提出的上述方法,可消除或减轻航空发动机在高空低速飞行包线中的振动量,达到抑制高空小表速抖振的目的。
附图说明
图1为本发明的流程图;
图2为实施例中的稳态特性和抗干扰特性的折中设计示意图。
具体实施方式
下面结合附图和具体实施例,对本发明作进一步详细说明。应理解这些实例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图1所示,本发明提供一种航空发动机高空小表速抖振抑制方法,具体包括如下步骤:
步骤1:获得航空发动机高空小表速燃油与转速的开环和闭环传递函数。
首先采用系统辨识法或者小扰动法,获得航空发动机高空小表速燃油与转速的发动机开环传递函数:
其中:s为拉普拉斯复变量;GE(s)为发动机开环传递函数;a、b、c和d为实系数。
航空发动机在高空、低速飞行时,全权限电子控制器采用比例控制(积分控制的系数极小以避免积分饱和),因此闭环传递函数可表示为:
其中:G(s)为发动机闭环传递函数,k为全权限电子控制器增益。
本实施例中,采用典型的发动机传递函数其中ωn为发动机自然频率;ξ为发动机阻尼系数。
步骤2:定义闭环控制对航空发动机干扰振动的放大函数。
为衡量闭环控制对航空发动机干扰振动的作用,特定义干扰振动放大函数为Td(jω),且表示为:
其中:GE(jω)和G(jω)分别为传递函数GE(s)和G(s)的频域响应函数。由于航空发动机为一惯性系统,因此可得:
其中,j为纯虚数,对某一频率ω,Td(jω)为一复数,其幅值即表达了闭环控制对该三角函数干扰情况下的放大作用:
针对本实施例,则有:
步骤3:设计最优闭环控制增益以达到抗干扰性能指标
为了使上述干扰不引起过量振动,即产生抖振,必须限制Td(jω)的幅值,即:
|Td(jω)|≤σ;
其中σ为一正实数,通常要求σ≤1,即如果σ≤1,则干扰信号通过航空发动机闭环控制系统后,幅值会被衰减;而如果σ>1则会引起振动增强。显然,当σ超过某一临界值时,将会引起发动机强烈振动,形成高空小表速抖振。满足上述要求的全权限电子控制器增益k为:
该式奠定了航空发动机高空小表速抖振抑制的设计基础。针对本实施例,则需要:
若要求ω=ωn时σ=1,则有:
也就是说全权限电子控制器增益k应满足:
高空小表速时,发动机呈现过阻尼特性,取ξ=2则有全权限电子控制器增益k必须满足|k|≤1.03。
步骤4:检验最优闭环控制增益满足发动机动态和稳态性能要求
由于全权限电子控制器增益k的设计还要满足航空发动机对动态和稳态特性的要求,因此必须检验步骤(3)中的增益设计,满足动态和稳态性能指标。如稳态性能定义为对某一阶跃转速要求的稳态误差ess,由拉普拉斯变换特性的终值定理可得:
稳态性能指标定义为:
|ess|≤τ;
其中,τ为一正实数。也就是要求全权限电子控制器增益k的设计满足:
因此,为同时满足本步骤要求的稳态特性和步骤(3)要求的抗干扰特性,全权限电子控制器增益k的设计必须同时满足和
针对本实施例,稳态性能可计算得到:
根据适航要求,该稳态误差不应大于1%,即:
由此可到k≥99。
步骤5:确认最优闭环控制增益
如果全权限电子控制器增益k同时满足步骤(3)和步骤(4)提出的稳态性能指标和抗干扰性能指标,则确认该最优设计;如果不同时满足上述两个步骤提出的性能指标要求,则必须对相应指标要求修改,因此也就必须对上述性能进行折中考虑。
针对实施例中,步骤(3)要求的抗干扰特性的全权限电子控制器增益|k|≤1.03;而步骤(4)要求的稳态性能需要k≥99,显然两个性能指标不能同时满足。例如,取k=0.5,则满足|k|≤1.03,抗干扰特性得到改善;但是不满足稳态性能要求,甚至严重退化,见附图2。因此,必须返回步骤3和步骤4,根据本发明所述方法开展性能综合分析,才能确认性能指标达到期望值。本实施例中,由于稳态性能通常要求采用“高增益控制”方法,实际上本实施例中,不能使得两个性能指标同时满足,这就确认了航空发动机稳态和动态控制设计必须和发动机抗干扰特性做出平衡的新的“设计准则”,但是这在实际设计过程中经常被忽略。这也正是高空小表速抖振问题一直没有得到解决的重要原因。
总之,本发明针对航空发动机高空小表速抖振问题,提出一种可消除或减轻航空发动机在高空低速飞行包线中的振动量,达到抑制高空小表速抖振的目的;同时本发明也指出了抑制高空小表速抖振可能与发动机稳态性能的指标要求存在折中,这个折中的设计理念不仅解释了高空小表速抖振的机理,也对实际工程具有重要价值。
Claims (1)
1.一种航空发动机高空小表速抖振抑制方法,其特征在于,包括以下步骤:
(1)获得航空发动机高空小表速燃油与转速的开环和闭环传递函数;
(2)定义闭环控制对航空发动机干扰振动的放大函数;
(3)设计最优闭环控制增益以达到抗干扰性能指标;
(4)检验最优闭环控制增益是否满足发动机动态和稳态性能要求;
(5)如果控制增益同时满足步骤(3)和步骤(4)提出的抗干扰性能指标和稳态性能指标,则为最优设计;否则返回步骤(3)和(4);
所述步骤(1)包括以下步骤:
(11)获得航空发动机高空小表速燃油与转速的发动机开环传递函数:
其中,GE(s)为发动机开环传递函数;s为拉普拉斯复变量;a、b、c和d为实系数;
(12)航空发动机在高空、低速飞行时,全权限电子控制器采用比例控制,闭环传递函数表示为:
其中,G(s)为发动机闭环传递函数,k为全权限电子控制器增益;
所述步骤(2)的实现过程如下:
衡量闭环控制对航空发动机干扰振动的作用,定义干扰振动放大函数为Td(jω),表示为:
其中,ω为频率,GE(jω)和G(jω)分别为传递函数GE(s)和G(s)的频域响应函数;由于航空发动机为一惯性系统,因此可得:
其中,a、b、c和d为实系数,Td(jω)为一复数,其幅值即表达了闭环控制对干扰振动放大函数干扰情况下的放大作用:
所述步骤(3)的实现过程如下:
限制Td(jω)的幅值,即:
|Td(jω)|≤σ
其中,σ为一正实数,要求σ≤1,满足要求的全权限电子控制器增益k为:
所述步骤(4)的实现过程如下:
将稳态性能定义为对某一阶跃转速要求的稳态误差ess,由拉普拉斯变换特性的终值定理可得:
稳态性能指标定义为:
|ess|≤τ
其中,τ为一正实数;即要求全权限电子控制器增益k的设计满足:
为同时满足稳态特性和步骤(3)要求的抗干扰特性,全权限电子控制器增益k的设计必须同时满足和
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010697863.4A CN112016158B (zh) | 2020-07-20 | 2020-07-20 | 一种航空发动机高空小表速抖振抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010697863.4A CN112016158B (zh) | 2020-07-20 | 2020-07-20 | 一种航空发动机高空小表速抖振抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112016158A CN112016158A (zh) | 2020-12-01 |
CN112016158B true CN112016158B (zh) | 2024-06-25 |
Family
ID=73498646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010697863.4A Active CN112016158B (zh) | 2020-07-20 | 2020-07-20 | 一种航空发动机高空小表速抖振抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112016158B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108443022A (zh) * | 2018-03-01 | 2018-08-24 | 南京航空航天大学 | 变旋翼转速飞行器发动机扭振抑制方法及装置 |
CN108762089A (zh) * | 2018-06-15 | 2018-11-06 | 大连理工大学 | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110850740B (zh) * | 2019-11-21 | 2021-04-20 | 大连理工大学 | 基于改进型Smith预估器的航空发动机H∞控制方法 |
-
2020
- 2020-07-20 CN CN202010697863.4A patent/CN112016158B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108443022A (zh) * | 2018-03-01 | 2018-08-24 | 南京航空航天大学 | 变旋翼转速飞行器发动机扭振抑制方法及装置 |
CN108762089A (zh) * | 2018-06-15 | 2018-11-06 | 大连理工大学 | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112016158A (zh) | 2020-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9695752B2 (en) | Engine-induced aircraft cabin resonance reduction system and method | |
Lietzau et al. | Model based control concepts for jet engines | |
CN112016158B (zh) | 一种航空发动机高空小表速抖振抑制方法 | |
US11015531B2 (en) | Engine-induced aircraft cabin resonance reduction system and method | |
Chen et al. | Aero‐Engine Real‐Time Models and Their Applications | |
CN111881530B (zh) | 一种航空发动机减振优化设计方法 | |
CN105628393B (zh) | 一种航空发动机超温虚警的排故方法 | |
CN111597633B (zh) | 航空发动机与挂架耦合减振的刚度反馈设计方法 | |
Wang et al. | Study on adaptive torsional vibration suppression methods for helicopter/turboshaft engine system with variable rotor speed | |
CN112881025B (zh) | 一种航空发动机振动控制与能量收集的方法 | |
KR20210019661A (ko) | 가변식 이득 조절기를 이용한 구조연동시험 방법 | |
DeLaat et al. | The high stability engine control (HISTEC) program-flight demonstration phase | |
CN115826412A (zh) | 一种航空发动机的无扰控制方法、装置、系统和存储介质 | |
CN114355766A (zh) | 一种飞行器控制系统三回路控制参数的设计方法 | |
Stapelfeldt et al. | On the importance of engine-representative models for fan flutter predictions | |
CN113075882A (zh) | 基于非平衡态线性化的涡扇发动机多变量鲁棒变增益控制方法 | |
CN114995121B (zh) | 一种航空电动燃油泵自适应滑模预见流量控制器设计方法 | |
Rebeske Jr et al. | Acceleration of high-pressure-ratio single-spool turbojet engine as determined from component performance characteristics: effect of compressor interstage air bleed | |
CN115370512B (zh) | 一种内燃机实验用外部空气压强调节系统 | |
Linyuan et al. | Steady state control schedule optimization for a variable cycle engine | |
Wang | Aeroengine High‐Attitude/Low Mach Number Oscillations: Mechanism and Prevention Design | |
CN116677499A (zh) | 一种航空发动机慢车稳定工作控制方法 | |
CN118686701A (zh) | 发动机加减速防喘控制方法、装置、电子设备及存储介质 | |
CN108757188B (zh) | 船舶推进柴油机双重负荷的保护方法 | |
ROONEY et al. | Correlation full scale wind tunnel and flight measured aerodynamic drag |
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 |