CN103198193A - 基于一阶模态幅值斜率的压气机旋转失速预测方法及系统 - Google Patents
基于一阶模态幅值斜率的压气机旋转失速预测方法及系统 Download PDFInfo
- Publication number
- CN103198193A CN103198193A CN2013101273336A CN201310127333A CN103198193A CN 103198193 A CN103198193 A CN 103198193A CN 2013101273336 A CN2013101273336 A CN 2013101273336A CN 201310127333 A CN201310127333 A CN 201310127333A CN 103198193 A CN103198193 A CN 103198193A
- Authority
- CN
- China
- Prior art keywords
- rotating stall
- sensor
- order mode
- processing unit
- central processing
- 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
Images
Landscapes
- Control Of Positive-Displacement Air Blowers (AREA)
Abstract
本发明公开了一种基于一阶模态幅值斜率的压气机旋转失速预测方法及系统。本发明通过设置在压气机的壁面的传感器采集压力信号,将压力信号转变为数字信号,中央处理器根据压力分布得到流量分布,并分析得到一阶模态,进而得到预测指数,并与预定的阈值相比较,判断得到是否会发生旋转失速,从而能够提前预测旋转失速。本发明装置简单,容易实现,计算量小且速度快,适用于各种工况,减少了每种工况单独处理的复杂度,具有理论基础,并且在进入旋转失速之前发出预警信号,提前控制,解决了旋转失速控制困难的问题,从而避免损失。
Description
技术领域
本发明涉及航空涡扇发动机领域,具体涉及一种基于一阶模态幅值斜率的压气机旋转失速预测方法及系统。
背景技术
航空发动机被誉为飞机的“心脏”,没有高性能的航空发动机,就不会有先进的军用飞机和具市场竞争力的民用客机。在航空发动机领域,我国与发达航空大国之间尚存在代差,在可靠性、稳定性和效率等各种性能指标方面还有相当的差距。
旋转失速是发动机内流的系统性失稳之一,它限定了发动机的稳定工作区域。发动机包括进气口、压气机、涡流及出气口。旋转失速是沿压气机周向的非均匀轴向流动状态,而其低流速区沿着与转子转动相同方向转动。旋转失速一旦产生,发动机流量、压比和推力都会突然降低,从而导致涡轮过热、甚至烧毁,而发动机必须关闭再重新启动才能恢复到正常工况。
这里必须强调旋转失速一般都是亚临界Hopf分岔,并伴随有迟滞回线。这种亚临界特性使得稳定边界对于外界干扰(进气畸变)非常敏感,即使在线性稳定区域,较大干扰也可能触发这些流动不稳定性。这就决定了旋转失速的控制是很困难的。由于旋转失速发生的时间非常短,通常是几十毫秒量级的,所以对控制器的作动速率提出了很高的要求。同时,发动机进入旋转失速后,需要巨大的能量才能将其拉出,这就决定了控制器的带宽必须非常大。
基于旋转失速控制的难度,解决方法之一便是提前预测旋转失速,在其进入旋转失速之前发出预警信号,提前控制。
发明内容
针对旋转失速控制难度大,本发明提出了一种提前预测的方法,这种方法解决了在压气机进入旋转失速之后才进行控制的问题。
本发明的一个目的在于提出一种基于一阶模态幅值斜变化率的压气机旋转失速预测方法。
本发明的一种基于一阶模态幅值斜变化率的压气机旋转失速预测方法,包括以下步骤:
1)传感器采集压力信号:在压气机级前周向,沿壁面均匀地布置一周动压传感器,传感器的数目为N个,每相邻的两个传感器之间的角度为N≥3为自然数;
2)处理信号:数据处理系统将传感器采集到的信号进行信号处理,传送给中央处理器;
3)中央处理器得到当前的压力分布,第i个传感器所在位置的压力为Pi,通过热线矫正,利用伯努利方程,可以得到当前的流量分布φi,i为自然数,且1≤i≤N;
4)分析预测信号:利用空间离散傅里叶变换,计算得到干扰的一阶模态a1为
其中,W为傅里叶矩阵
5)提取了一阶模态之后,得到一个预测指数:
其中,J=|a1|2,是一阶模态幅值的平方,t为采样时间,根据实际情况,设定一个阈值S0,预测指数S与阈值S0相比较,当S>S0时,认为旋转失速将要发生,否则,则认为压气机工作在较安全的区域。
其中,在步骤1)中,传感器的数目为N个(N≥3为自然数),根据奈奎斯特采样定理,若要采集空间一阶模态,则需要三个以上的传感器。进一步,N为奇数,传感器的数目为奇数个,奇数个传感器可以比偶数个传感器少用一个传感器。
在步骤2)中,数据处理系统模数转换器和滤波器,对信号进行模数转换以及滤波。
由于传感器采集信号的时间为离散的时间点t1、t2……tm,因此预测指数S为离散值。在步骤5)中,预测指数S具体算法为:对于离散的时间点t1、t2……tm,得到相应的一阶模态幅值的平方的值分别为J1、J2……Jm,则其中,m为自然数,1≤j≤m-1的自然数。
本发明的另一个目的在于提供一种基于一阶模态幅值斜率的压气机旋转失速预测系统。
本发明的压气机旋转失速预测系统包括:N个传感器、N个数据处理系统、N条热线以及中央处理器;其中,N个传感器沿压气机的壁面均匀放置;N个传感器中的每一个分别连接至数据处理系统;N个数据处理系统中连接至中央处理器;N条热线连接至中央处理器,N≥3为自然数。
N个数据处理系统中的每一个包括模数转换器和滤波器。
进一步,传感器的数目为奇数个。
通常流体中的计算量都是相当大的,本发明采用动力学分析方法,只需要进行一次空间傅里叶变换和一次求斜率,计算速度非常快;对各种工况具有普适性,由于采用无量纲化,将各种复杂的工况归一化,减少了每种工况单独处理的复杂度;相对于以前的方法提前了许多,从实验结果可以看出,在旋转失速之前很多转,斜率已经开始增长;本方法可靠性高,具有理论基础,只要有旋转失速的现象存在,就会有一阶模态幅值斜率的增长。
本发明的优点:
本发明通过设置在压气机的壁面的传感器采集压力信号,分析得到一阶模态,进而得到预测指数,并与预定的阈值相比较,判断得到是否会发生旋转失速,从而能够提前预测旋转失速。本发明装置简单,容易实现,计算量小且速度快,适用于各种工况,减少了每种工况单独处理的复杂度,具有理论基础,并且在进入旋转失速之前发出预警信号,提前控制,解决了旋转失速控制困难的问题,从而避免损失。
附图说明
图1为本发明的基于一阶模态幅值斜率的压气机旋转失速预测方法的系统的示意图;
图2为本发明的一个实施例的传感器的设置的结构示意图;
图3为本发明的基于一阶模态幅值斜率的压气机旋转失速预测方法的流程图;
图4为本发明的一个实施例的进气口放置挡板的结构示意图;
图5为本发明的一个实施例的一阶模态幅值的平方的斜率随转数变化的曲线图。
具体实施方式
下面结合附图,通过实例对本发明做进一步说明。
如图1所示,本发明采用的系统包括:N个动压传感器T1~TN、N个数据处理系统A/D、N条热线HL1~HLN、以及中央处理器CPU;其中,N个动压传感器中的每一个分别连接至数据处理系统;N个数据处理系统中连接至中央处理器;N条热线连接至中央处理器。
在本实施例中,采用5个传感器,每相邻的两个传感器之间的角度为72°,传感器的设置如图2所示。
本实施例的一种基于一阶模态幅值斜率的压气机旋转失速预测方法,如图3所示,包括以下步骤:
1)传感器采集压力信号:在压气机级前周向,沿壁面均匀的布置一周动压传感器,传感器的数目为5个,每相邻的两个传感器之间的角度为72°;
2)处理信号:数据处理系统模数转换器和滤波器,对信号进行模数转换以及滤波,传送给中央处理器;
3)中央处理器得到当前的压力分布,第i个传感器所在位置的压力为Pi,通过热线矫正,利用伯努利方程,可以得到当前的流量分布φi,i为自然数,且1≤i≤5,对流量φi无量纲化为流量系数其中,Un表达式为Un=2πfR,f为转子转速,R为平均半径;
4)分析预测信号:利用空间离散傅里叶变换,可以得到干扰的一阶模态a1为
其中,W为傅里叶矩阵
Φ=[Φ1......Φ5]T;
5)提取了一阶模态之后,得到一个预测指数:
其中,J=|a1|2,是一阶模态幅值的平方,t为采样时间,对于离散的时间点t1、t2……tm,得到相应的一阶模态幅值的平方的值分别为J1、J2……Jm,其中,m为自然数,1≤j≤m的自然数,根据实际情况,设定一个阈值S0,当S>S0时,认为旋转失速将要发生,否则,则认为压气机工作在较安全的区域。
进气畸变是造成旋转失速的最主要原因,为了验证方法的可行性,我们在各种进气畸变的实验条件下进行了一系列的实验。
实验结果如图5所示,其中横轴是压气机动子的转数,每条线代表了不同的进气畸变大小。从实验结果可以看出,随着转数的增加,压气机原来有接近旋转失速,dJ/dt也越来越大,因此,本发明的预测方法能够预先判定是否会发生旋转失速。
最后需要注意的是,公布实施方式的目的在于帮助进一步理解本发明,但是本领域的技术人员可以理解:在不脱离本发明及所附的权利要求的精神和范围内,各种替换和修改都是可能的。因此,本发明不应局限于实施例所公开的内容,本发明要求保护的范围以权利要求书界定的范围为准。
Claims (7)
1.一种压气机旋转失速预测方法,其特征在于,所述预测方法包括以下步骤:
2)处理信号:数据处理系统将传感器采集到的信号进行信号处理,传送给中央处理器;
3)中央处理器得到当前的压力分布,第i个传感器所在位置的压力为Pi,通过热线矫正,利用伯努利方程,得到当前的流量分布φi,i为自然数,且1≤i≤N;
4)分析预测信号:利用空间离散傅里叶变换,计算得到干扰的一阶模态a1为
其中,W为傅里叶矩阵
5)提取了一阶模态之后,得到一个预测指数:
其中,J=|a1|2,是一阶模态幅值的平方,t为采样时间,根据实际情况,设定一个阈值S0,预测指数S与阈值S0相比较,当S>S0时,认为旋转失速将要发生,否则,则认为压气机工作在较安全的区域。
2.如权利要求1所述的预测方法,其特征在于,在步骤2)中,所述数据处理系统模数转换器和滤波器,对信号进行模数转换以及滤波。
5.一种压气机旋转失速预测系统,其特征在于,所述预测系统包括:N个传感器、N个数据处理系统、N条热线以及中央处理器;其中,N个传感器沿压气机的壁面均匀放置;N个传感器中的每一个分别连接至数据处理系统;N个数据处理系统中连接至中央处理器;N条热线连接至中央处理器,N≥3为自然数。
6.如权利要求5所述的预测系统,其特征在于,N个数据处理系统中的每一个包括模数转换器和滤波器。
7.如权利要求5所述的预测系统,其特征在于,所述传感器的个数N为奇数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310127333.6A CN103198193B (zh) | 2013-04-12 | 2013-04-12 | 基于一阶模态幅值斜率的压气机旋转失速预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310127333.6A CN103198193B (zh) | 2013-04-12 | 2013-04-12 | 基于一阶模态幅值斜率的压气机旋转失速预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103198193A true CN103198193A (zh) | 2013-07-10 |
CN103198193B CN103198193B (zh) | 2015-12-02 |
Family
ID=48720748
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310127333.6A Expired - Fee Related CN103198193B (zh) | 2013-04-12 | 2013-04-12 | 基于一阶模态幅值斜率的压气机旋转失速预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103198193B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425784A (zh) * | 2015-12-30 | 2016-03-23 | 中国航空综合技术研究所 | 一种发动机控制软件的检查方法 |
CN105629952A (zh) * | 2015-12-30 | 2016-06-01 | 中国航空综合技术研究所 | 一种发动机控制软件的测试方法 |
CN107576445A (zh) * | 2017-07-24 | 2018-01-12 | 西北工业大学 | 基于单路传感器信号的压气机失速先兆信号重构与识别方法 |
CN109214141A (zh) * | 2018-11-20 | 2019-01-15 | 西华大学 | 旋转失速预测方法及装置 |
CN109339954A (zh) * | 2018-10-08 | 2019-02-15 | 南京航空航天大学 | 一种基于航空发动机压气机部件气动失稳的主动控制方法 |
CN110382878A (zh) * | 2017-03-02 | 2019-10-25 | 柏林工业大学 | 确定用于预测压缩机中不稳定性的指标的方法和装置及其用途 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030061601A1 (en) * | 2001-09-26 | 2003-03-27 | Nec Corporation | Data processing apparatus and method, computer program, information storage medium, parallel operation apparatus, and data processing system |
US20070179763A1 (en) * | 2006-01-27 | 2007-08-02 | Ricardo, Inc. | Apparatus and method for compressor and turbine performance simulation |
US20080183441A1 (en) * | 2004-12-02 | 2008-07-31 | Mehdi Vahdati | Rotating stall |
US20080232950A1 (en) * | 2007-03-23 | 2008-09-25 | Johnson Controls Technology Company | Method for detecting rotating stall in a compressor |
CN101882177A (zh) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | 叶轮机械计入叶间相位角的气弹稳定性流固耦合预测方法 |
CN101887479A (zh) * | 2010-07-23 | 2010-11-17 | 华南理工大学 | 基于确定学习理论的轴流压气机旋转失速的快速诊断方法 |
CN101908088A (zh) * | 2010-07-22 | 2010-12-08 | 北京航空航天大学 | 一种基于时域双向迭代的叶轮机叶片颤振应力预测方法 |
-
2013
- 2013-04-12 CN CN201310127333.6A patent/CN103198193B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030061601A1 (en) * | 2001-09-26 | 2003-03-27 | Nec Corporation | Data processing apparatus and method, computer program, information storage medium, parallel operation apparatus, and data processing system |
US20080183441A1 (en) * | 2004-12-02 | 2008-07-31 | Mehdi Vahdati | Rotating stall |
US20070179763A1 (en) * | 2006-01-27 | 2007-08-02 | Ricardo, Inc. | Apparatus and method for compressor and turbine performance simulation |
US20080232950A1 (en) * | 2007-03-23 | 2008-09-25 | Johnson Controls Technology Company | Method for detecting rotating stall in a compressor |
CN101882177A (zh) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | 叶轮机械计入叶间相位角的气弹稳定性流固耦合预测方法 |
CN101908088A (zh) * | 2010-07-22 | 2010-12-08 | 北京航空航天大学 | 一种基于时域双向迭代的叶轮机叶片颤振应力预测方法 |
CN101887479A (zh) * | 2010-07-23 | 2010-11-17 | 华南理工大学 | 基于确定学习理论的轴流压气机旋转失速的快速诊断方法 |
Non-Patent Citations (5)
Title |
---|
WANG CONG等: "Identification and Rapid Detection of Rotating Stall via Deterministic Learing", 《PROCEEDINGS OF THE 29TH CHINESE CONTROL CONFERENCE》 * |
WANJING DONG等: "Construction of Lyapunov Function near Bifurcations with Application to Rotating stall", 《PROCEEDINGS OF THE 31TH CHINESE CONTROL CONFERENCE》 * |
周龙沙: "基于确定学习的航空发动机旋转失速快速检测的并行计算实现", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
童志庭: "轴流压气机叶尖泄漏锅、失速先兆、叶尖微喷气非定常关联性的实验研究", 《万方学位论文全文数据库》 * |
谢锋等: "基于小波分析的航空发动机旋转失速检测方法", 《航空动力学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425784A (zh) * | 2015-12-30 | 2016-03-23 | 中国航空综合技术研究所 | 一种发动机控制软件的检查方法 |
CN105629952A (zh) * | 2015-12-30 | 2016-06-01 | 中国航空综合技术研究所 | 一种发动机控制软件的测试方法 |
CN110382878A (zh) * | 2017-03-02 | 2019-10-25 | 柏林工业大学 | 确定用于预测压缩机中不稳定性的指标的方法和装置及其用途 |
CN110382878B (zh) * | 2017-03-02 | 2020-12-08 | 柏林工业大学 | 确定用于预测压缩机中不稳定性的指标的方法和装置及其用途 |
CN107576445A (zh) * | 2017-07-24 | 2018-01-12 | 西北工业大学 | 基于单路传感器信号的压气机失速先兆信号重构与识别方法 |
CN107576445B (zh) * | 2017-07-24 | 2019-09-24 | 西北工业大学 | 基于单路传感器信号的压气机失速先兆信号重构与识别方法 |
CN109339954A (zh) * | 2018-10-08 | 2019-02-15 | 南京航空航天大学 | 一种基于航空发动机压气机部件气动失稳的主动控制方法 |
CN109214141A (zh) * | 2018-11-20 | 2019-01-15 | 西华大学 | 旋转失速预测方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103198193B (zh) | 2015-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103198193B (zh) | 基于一阶模态幅值斜率的压气机旋转失速预测方法及系统 | |
DK179081B1 (da) | Overvågning af en vindmølles vingefrekvenser | |
CN104345118B (zh) | 固体推进剂多靶线动态燃烧性能测试系统及方法 | |
US10496775B2 (en) | Method and system for use in dynamically configuring data acquisition systems | |
EP3401600B1 (en) | Systems and methods related to detecting flame blowout occurrences in gas turbines | |
CN110848166A (zh) | 轴流压气机喘振频率预测方法 | |
CN104888396A (zh) | 一种检测消防水泵在快速启动过程中瞬态故障的方法 | |
CN105351153A (zh) | 风电机组变桨控制器的性能评估方法 | |
CN105829983B (zh) | 用于检测机器的当前的损伤状态的设备以及方法 | |
US9447735B2 (en) | Method of controlling a turbomachine | |
CN104712378A (zh) | 火电机组主蒸汽压力闭环节能控制方法和系统 | |
CN114156910A (zh) | 一种基于信号同源的agc负荷解耦控制方法和自动寻优装置 | |
GB2522847A (en) | Method and system for detecting a flow blockage in a pipe | |
US10954812B2 (en) | Gas turbine blade flutter monitoring and control system | |
CN104200107A (zh) | 不同飞行迎角对航空发动机失速/喘振特性影响的适航审定方法 | |
CN102900564B (zh) | 用于轴流压缩系统叶顶喷气稳定性控制的控制器 | |
JP2003314305A (ja) | ガスタービン制御装置、ガスタービンシステム及びガスタービン制御方法 | |
Zhao et al. | Experimental investigation of transient characteristics of mild surge and diffuser rotating stall in a centrifugal compressor with vaned diffuser | |
US8342010B2 (en) | Surge precursor protection systems and methods | |
CN107560864B (zh) | 用于燃烧器中结垢监测和预测的方法与装置 | |
Wang et al. | Triaxial gas turbine performance analysis for variable power turbine inlet guide vane control law optimization | |
CN110761946B (zh) | 一种风机塔架振动控制方法及装置 | |
CN109869227B (zh) | 基于大数据运算中心的柴油发电机组 | |
Li et al. | Automatic efficiency optimization of an axial compressor with adjustable inlet guide vanes | |
CN112555050B (zh) | 一种分布式超音速空气流生成装置运行状态监测及急停方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151202 Termination date: 20180412 |