CN114458611B - 一种基于压气机出口压力的失速与喘振机载识别方法 - Google Patents

一种基于压气机出口压力的失速与喘振机载识别方法 Download PDF

Info

Publication number
CN114458611B
CN114458611B CN202111432796.4A CN202111432796A CN114458611B CN 114458611 B CN114458611 B CN 114458611B CN 202111432796 A CN202111432796 A CN 202111432796A CN 114458611 B CN114458611 B CN 114458611B
Authority
CN
China
Prior art keywords
outlet pressure
compressor
surge
mark
stall
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
Application number
CN202111432796.4A
Other languages
English (en)
Other versions
CN114458611A (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.)
AECC Shenyang Engine Research Institute
Original Assignee
AECC Shenyang Engine Research Institute
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 AECC Shenyang Engine Research Institute filed Critical AECC Shenyang Engine Research Institute
Priority to CN202111432796.4A priority Critical patent/CN114458611B/zh
Publication of CN114458611A publication Critical patent/CN114458611A/zh
Application granted granted Critical
Publication of CN114458611B publication Critical patent/CN114458611B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04D—NON-POSITIVE-DISPLACEMENT PUMPS
    • F04D27/00—Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
    • F04D27/001—Testing thereof; Determination or simulation of flow characteristics; Stall or surge detection, e.g. condition monitoring
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • 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
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/10—Internal combustion engine [ICE] based vehicles
    • Y02T10/12—Improving ICE efficiencies

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Mechanical Engineering (AREA)
  • Control Of Positive-Displacement Air Blowers (AREA)

Abstract

本申请提供来一种基于压气机出口压力的失速与喘振机载识别方法，所述方法包括：根据所述压气机出口压力信号确定所述压气机出口压力的标准差及平均值，根据压气机出口压力的标准差及平均值确定压气机出口压力的紊流度；根据紊流度阈值判断是否生成失速标识；根据所述压气机出口压力信号确定所述压气机出口压力的变化率及平均值，根据压气机出口压力的变化率及平均值确定压气机出口压力的相对变化率；根据压气机出口压力的相对变化率判断是否生成喘振标识；若生成任意一个的失速标识和喘振标识，则生成失稳标识，即发动机处于失稳状态。该方法既能降低喘振识别时间、保证失速识别的可靠性，又能够排除停车、测压管路破裂等特殊情况的干扰而不误判。

Description

一种基于压气机出口压力的失速与喘振机载识别方法
技术领域
本申请属于航空发动机控制技术领域，特别涉及一种基于压气机出口压力的失速与喘振机载识别方法。
背景技术
喘振与失速是发动机压气机的两种主要失稳状态。喘振时，压气机内的气流在轴线方向进行低频率、高振幅的振荡，发动机整机表现为振动增大、发出低沉噪声、燃气温度升高、转速摆动或下降，轻则造成发动机工况急剧恶化，重则导致发动机停车或机械损伤，严重危及飞行安全。失速时，压气机内的气流出现较高频率的脉动，可能使叶片发生共振，是导致压气机叶片疲劳断裂的主要原因之一，而且失速极度发展之后往往导致喘振。
为了避免发动机出现严重损伤，发动机控制系统需要具备失速与喘振的识别功能。
目前的发动机失速与喘振识别方法中主要有：紊流度分析法、脉动压力方差法、压力变化率法、频谱分析法、小波分析法、李雅谱诺夫稳定性分析、相关积分法、自相关检测法等。其中，小波分析、李雅谱诺夫稳定性分析和相关积分法算法较为复杂，主要用于信号的事后分析；紊流度分析法、脉动压力方差法和压力变化率法、自相关检测法是基于信号时域波形幅值特征的识别方法；频谱分析是基于频域特征的识别方法。
不同方法所需的信号采样周期和计算量差异较大，针对发动机机载情况，受机载电子控制器计算能力的限制，计算量较小的紊流度分析法和压力变化率法是目前工程中常用的机载失稳识别方法。
而目前工程中常用的机载紊流度分析法和压力变化率法存在以下缺点：
1)单独使用紊流度分析法识别失速与喘振，识别喘振需要的时间较长，不利于喘振的及时处置；
2)单独使用压力变化率法识别失速与喘振，识别失速的可靠性不高；
3)常规的紊流度分析法和压力变化率法在停车、测压管路破裂、压气机/涡轮结构失效、意外熄火等特殊情况下容易误判。
发明内容
本申请的目的是提供了一种基于线段拟合在线识别间隙环节特征参数的方法及装置，以解决或减轻上述至少一个问题。
在一方面，本申请提供了一种基于压气机出口压力的失速与喘振机载识别方法，所述方法包括：
根据所述压气机出口压力信号确定所述压气机出口压力的标准差及平均值，根据压气机出口压力的标准差及平均值确定压气机出口压力的紊流度；
若当前周期的压气机出口压力的紊流度大于紊流度阈值则生成紊流度阈值标识，若自当前周期连续向后m1个周期均具有紊流度阈值标识，且自第m1周期连续向前m2个周期的压气机出口压力与压气机出口压力的平均值之差不全为负值且不全为正值，m1＞m2，则生成失速标识；
根据所述压气机出口压力信号确定所述压气机出口压力的变化率及平均值，根据压气机出口压力的变化率及平均值确定压气机出口压力的相对变化率；
若当前周期的压气机出口压力的相对变化率小于突降阈值，且向后m3个周期的压气机出口压力的相对变化率均小于突降阈值，则生成突降标识；若当前周期的压气机出口压力的相对变化率大于突降阈值，且向后m4个周期的压气机出口压力的相对变化率均大于突降阈值，则生成突升标识；若具有突升标识，且前m5个周期内具有突降标，则生成喘振标识；
若生成任意一个的失速标识和喘振标识，则生成失稳标识，即发动机处于失稳状态。
进一步的，所述压气机出口压力的标准差为：
Figure GDA0004221704240000031
式中，σ为标准差，n为样本容量，P3(i)(i＝1～n)为包括当前周期的最新n个周期P3的输入数值；P3为P3(i)(i＝1～n)的算术平均值。
进一步的，所述压气机出口压力的标准差能够简化为：
Figure GDA0004221704240000032
进一步的，所述压气机出口压力的紊流度为：
Figure GDA0004221704240000033
式中，k表示当前周期。
进一步的，所述压气机出口压力的变化率为：
Figure GDA0004221704240000034
其中，T为采样周期，Y(k)为压气机出口压力的变化率，X(k)为压气机出口压力的采集值，k为当前周期，k-1为前1个周期，k-2为前2个周期，k-a表示当前周期的前a个周期。
进一步的，所述压气机出口压力的相对变化率为：
Figure GDA0004221704240000035
另一方面，本申请提供了一种发动机控制系统，所述发动机控制系统至少包括一数据处理装置，所述数据处理装置执行如上任一所述的基于压气机出口压力的失速与喘振机载识别方法。
本申请的基于压气机出口压力的失速与喘振机载识别方法既能降低喘振识别时间，又能保证失速识别的可靠性，同时该方法能够排除停车、测压管路破裂、压气机/涡轮结构失效、意外熄火等特殊情况的干扰而不误判。
附图说明
为了更清楚地说明本申请提供的技术方案，下面将对附图作简单地介绍。显而易见地，下面描述的附图仅仅是本申请的一些实施例。
图1为本申请的识别方法整体流程示意图。
图2为本申请的紊流度分析逻辑示意图。
图3为本申请的压力变化率逻辑示意图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行更加详细的描述。
本申请中提出了一种基于压气机出口压力的失速与喘振机载识别方法，该方法通过综合运用紊流度分析法和压力变化率法，既能降低喘振识别时间，又能保证失速识别的可靠性，同时该方法能够排除停车、测压管路破裂、压气机/涡轮结构失效、意外熄火等特殊情况的干扰而不误判。
如图1所示，本申请的失速与喘振识别方法基于压气机出口压力(总压或静压不限，可用P3表示)，运用紊流度分析法和压力变化率法综合实现失速与喘振的识别，具体过程如下：
一、紊流度分析逻辑
紊流度能够反映信号的脉动程度，通过紊流度用作失速的判定参数。
1.1)标准差计算
在本申请一实施例中，压气机出口压力P3的标准差σ通过以下计算：
Figure GDA0004221704240000051
式中，n为样本容量，根据工程项目的具体情况确定；P3(i)(i＝1～n)为包括当前周期的最新n个周期P3的输入数值；
Figure GDA0004221704240000054
为P3(i)(i＝1～n)的算术平均值。
进一步的，上式计算量较大，为了降低计算量，本申请中将上式简化为：
Figure GDA0004221704240000052
简化的标准差公式的计算结果与原公式差异较小，并且数值变化趋势完全一致，不影响逻辑判断结果。
需要说明的是，样本容量n对标准差的影响较大，样本容量偏大会使紊流度对发动机加减速比较敏感，样本容量偏小会使紊流度对低频失速不敏感，故应根据工程项目的具体情况折衷确定。
1.2)平均值计算
压气机进口压力的平均值为包括当前周期的最新m个周期数值的算术平均值，m根据工程项目的具体情况确定，本处不在赘述。
1.3)紊流度计算
在本申请中，压气机进口压力的紊流度采用如下公式计算：
Figure GDA0004221704240000053
其中，k表示当前周期数值。
1.4)阈值及阈值比较
不同的信号处理、滤波方法对紊流度阈值影响较大，因此紊流度阈值应根据工程项目的具体情况确定。
阈值比较逻辑如下：
若压气机进口压力P3的紊流度(k)＞紊流度阈值，则紊流度阈值比较结果为“1”，即生成紊流度阈值标识，反之紊流度阈值比较结果为“0”，即不生产紊流度阈值标识。
1.5)确认
若(紊流度阈值比较结果连续m1个周期为“1”)并且(自第m1个周期起，向前统计m2(m1＞m2)个周期，这m2个周期的压气机进口压力与上述的压气机进口压力平均值(m个周期)之差不全为负值且不全为正值，则生成失速标识，失速标识置“1”，反之则不生成失速标识，失速标识置“0”。
需要说明的是，周期数m1与m2根据工程项目的具体情况确定。
本申请中，通过m2个周期的压气机进口压力与压气机进口压力平均值之差的检查逻辑可避免发动机熄火或压气机进口压力测压管破裂、压气机/涡轮结构失效使压气机进口压力信号下降引起紊流度超限而导致的误判。
二、压力变化率逻辑
压力相对变化率能够快速反映信号的大幅度变化，本申请中用作喘振的判定参数。
2.1)变化率计算
本申请中的压气机出口压力的变化率计算方法为：
Figure GDA0004221704240000061
式中，T为采样周期，根据工程项目的具体情况确定，Y(k)为压气机出口压力的变化率，X(k)为压气机出口压力的采集值，k为当前周期，k-1为前1个周期，k-2为前2个周期，k-a表示当前周期的前a个周期。此公式相当于对变化率进行了一定程度的滤波。
但在该实施例中，压气机出口压力的变化率仅取前4个周期进行计算，因此有：
Figure GDA0004221704240000071
2.2)平均值计算
该压气机出口压力平均值与紊流度分析中压气机进口压力的m个周期的算数平均值计算过程相同。
2.3)相对变化率计算
本申请中，压气机出口压力大相对变化率计算公式为：
Figure GDA0004221704240000072
2.4)阈值
压气机出口压力变化率阈值包括突降阈值和突升阈值，根据工程项目的具体情况确定。
2.5)突降比较
若压气机出口压力相对变化率(k)＜突降阈值，则比较结果为“1”，否则比较结果为“0”。
2.6)突降确认
若突降比较结果连续m3个周期为“1”，则生成突降标识，突降标识置“1”。否则不生成突降标识，突降标识置“0”。
本申请中，周期数m3根据工程项目的具体情况确定，通常取2、3个周期。
2.7)突升比较
若压气机出口压力相对变化率(k)＞突升阈值，则比较结果为“1”，否则比较结果为“0”。
2.8)突升确认
若突升比较结果连续m4个周期为“1”，则生成突升表示，突升标识置“1”，否则不生成突升标识，突升标识置“0”。
本申请中，周期数m4根据工程项目的具体情况确定，通常取2、3个周期。
2.9)降/升综合
若(突升标识为“1”)并且(前m5个周期内有突降标识置“1”事件)，则生成喘振标识，喘振标识置“1”，否则不生产喘振标识，喘振标识置“0”。
本申请中，周期数m5根据工程项目的具体情况确定，通常按最长喘振周期的一半左右确定。
需要说明的是，本申请中，通过降/升综合逻辑可避免发动机熄火或压气机出口压力测压管破裂、压气机/涡轮结构失效导致的压气机出口压力信号单调下降的误判。
上述失速标识和喘振标识任意一个置为“1”，则产生失稳标识，即发动机处于失稳状态。
本申请的基于压气机出口压力的失速与喘振机载识别方法既能降低喘振识别时间，又能保证失速识别的可靠性，同时该方法能够排除停车、测压管路破裂、压气机/涡轮结构失效、意外熄火等特殊情况的干扰而不误判。
此外，本申请提供了一种发动机控制系统，所述发动机控制系统至少包括一数据处理装置，所述数据处理装置执行如上任一所述的基于压气机出口压力的失速与喘振机载识别方法。
以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。

Claims (7)

1.一种基于压气机出口压力的失速与喘振机载识别方法，其特征在于，所述方法包括：
根据所述压气机出口压力信号确定所述压气机出口压力的标准差及平均值，根据压气机出口压力的标准差及平均值确定压气机出口压力的紊流度；
若当前周期的压气机出口压力的紊流度大于紊流度阈值则生成紊流度阈值标识，若自当前周期连续向后m1个周期均具有紊流度阈值标识，且自第m1周期连续向前m2个周期的压气机出口压力与压气机出口压力的平均值之差不全为负值且不全为正值，m1＞m2，则生成失速标识；
根据所述压气机出口压力信号确定所述压气机出口压力的变化率及平均值，根据压气机出口压力的变化率及平均值确定压气机出口压力的相对变化率；
若当前周期的压气机出口压力的相对变化率小于突降阈值，且向后m3个周期的压气机出口压力的相对变化率均小于突降阈值，则生成突降标识；若当前周期的压气机出口压力的相对变化率大于突降阈值，且向后m4个周期的压气机出口压力的相对变化率均大于突降阈值，则生成突升标识；若具有突升标识，且前m5个周期内具有突降标，则生成喘振标识；
若生成任意一个的失速标识和喘振标识，则生成失稳标识，即发动机处于失稳状态。
2.如权利要求1所述的基于压气机出口压力的失速与喘振机载识别方法，其特征在于，所述压气机出口压力的标准差为：
Figure FDA0004221704230000011
式中，σ为标准差，n为样本容量，P3(i)(i＝1～n)为包括当前周期的最新n个周期P3的输入数值；
Figure FDA0004221704230000021
为P3(i)(i＝1～n)的算术平均值。
3.如权利要求2所述的基于压气机出口压力的失速与喘振机载识别方法，其特征在于，所述压气机出口压力的标准差简化为：
Figure FDA0004221704230000022
4.如权利要求3所述的基于压气机出口压力的失速与喘振机载识别方法，其特征在于，所述压气机出口压力的紊流度为：
Figure FDA0004221704230000023
式中，k表示当前周期。
5.如权利要求1所述的基于压气机出口压力的失速与喘振机载识别方法，其特征在于，所述压气机出口压力的变化率为：
Figure FDA0004221704230000024
其中，T为采样周期，Y(k)为压气机出口压力的变化率，X(k)为压气机出口压力的采集值，k为当前周期，k-1为前1个周期，k-2为前2个周期，k-a表示当前周期的前a个周期。
6.如权利要求5所述的基于压气机出口压力的失速与喘振机载识别方法，其特征在于，所述压气机出口压力的相对变化率为：
Figure FDA0004221704230000025
7.一种发动机控制系统，其特征在于，所述发动机控制系统至少包括一数据处理装置，所述数据处理装置执行如权利要求1至6中任一所述的基于压气机出口压力的失速与喘振机载识别方法。
CN202111432796.4A 2021-11-29 2021-11-29 一种基于压气机出口压力的失速与喘振机载识别方法 Active CN114458611B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111432796.4A CN114458611B (zh) 2021-11-29 2021-11-29 一种基于压气机出口压力的失速与喘振机载识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111432796.4A CN114458611B (zh) 2021-11-29 2021-11-29 一种基于压气机出口压力的失速与喘振机载识别方法

Publications (2)

Publication Number Publication Date
CN114458611A CN114458611A (zh) 2022-05-10
CN114458611B true CN114458611B (zh) 2023-06-20

Family

ID=81405936

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111432796.4A Active CN114458611B (zh) 2021-11-29 2021-11-29 一种基于压气机出口压力的失速与喘振机载识别方法

Country Status (1)

Country Link
CN (1) CN114458611B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110005628A (zh) * 2019-03-27 2019-07-12 南京航空航天大学 基于异位方差分析的压气机气动失稳在线辨识方法和系统
CN110608187A (zh) * 2019-10-30 2019-12-24 江西理工大学 基于频率特征变化的轴流压气机失速喘振预测装置
CN112460061A (zh) * 2020-12-16 2021-03-09 珠海格力电器股份有限公司 一种离心式压缩机失速状态判定方法、装置及机组设备
CN112539192A (zh) * 2019-09-20 2021-03-23 中国航发商用航空发动机有限责任公司 燃气轮机、燃烧室、压气机失速的监测装置、监测方法以及计算机可读存储介质

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6823254B2 (en) * 2003-03-28 2004-11-23 Honeywell International, Inc. Method and system for turbomachinery surge detection

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110005628A (zh) * 2019-03-27 2019-07-12 南京航空航天大学 基于异位方差分析的压气机气动失稳在线辨识方法和系统
CN112539192A (zh) * 2019-09-20 2021-03-23 中国航发商用航空发动机有限责任公司 燃气轮机、燃烧室、压气机失速的监测装置、监测方法以及计算机可读存储介质
CN110608187A (zh) * 2019-10-30 2019-12-24 江西理工大学 基于频率特征变化的轴流压气机失速喘振预测装置
CN112460061A (zh) * 2020-12-16 2021-03-09 珠海格力电器股份有限公司 一种离心式压缩机失速状态判定方法、装置及机组设备

Also Published As

Publication number Publication date
CN114458611A (zh) 2022-05-10

Similar Documents

Publication Publication Date Title
CN110131193B (zh) 航空发动机喘振故障监测方法和系统
US8152496B2 (en) Continuing compressor operation through redundant algorithms
JP4646528B2 (ja) 切迫センサ故障を検出するための方法及び装置
CN110735669B (zh) 一种航空燃气涡轮发动机旋转失速判断方法及装置
US6231306B1 (en) Control system for preventing compressor stall
US6438484B1 (en) Method and apparatus for detecting and compensating for compressor surge in a gas turbine using remote monitoring and diagnostics
Sun et al. Cyclic spectral analysis of vibration signals for centrifugal pump fault characterization
JP5337853B2 (ja) 発電プラントの診断装置、及び発電プラントの診断方法
US7409854B2 (en) Method and apparatus for determining an operating status of a turbine engine
US9250153B2 (en) Methods and systems for monitoring health of blades
CN108699966B (zh) 增压器的喘振检测方法以及喘振检测装置
JPH08503757A (ja) 圧縮機を監視し制御するための方法及び装置
JP4174031B2 (ja) ターボ機械のサージング限界またはブレード損傷の警告
US8869603B2 (en) Debris detection in turbomachinery and gas turbine engines
US6010303A (en) Apparatus and method of predicting aerodynamic and aeromechanical instabilities in turbofan engines
CN112580267A (zh) 基于多分支特征融合网络的航空发动机喘振预测方法
CN106678069A (zh) 离心式压缩机防喘振发生的检测方法
US20150044020A1 (en) Method for avoiding pump surges in a compressor
WO2019202917A1 (ja) 観測装置、観測方法及びプログラム
CN114458611B (zh) 一种基于压气机出口压力的失速与喘振机载识别方法
US8601861B1 (en) Systems and methods for detecting the flame state of a combustor of a turbine engine
US20210285457A1 (en) Surging precursor detecting device, method of detecting surging precursor, and program
CN111720218B (zh) 涡轮发动机的信号响应监测
CN114962305B (zh) 压气机失稳在线检测方法、装置、系统、设备及介质
JP4523826B2 (ja) ガスタービン監視装置及びガスタービン監視システム

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