CN104677486A - 基于转速脉冲重构的航空发动机振动信号相位测量方法 - Google Patents
基于转速脉冲重构的航空发动机振动信号相位测量方法 Download PDFInfo
- Publication number
- CN104677486A CN104677486A CN201310625488.2A CN201310625488A CN104677486A CN 104677486 A CN104677486 A CN 104677486A CN 201310625488 A CN201310625488 A CN 201310625488A CN 104677486 A CN104677486 A CN 104677486A
- Authority
- CN
- China
- Prior art keywords
- prime
- vibration signal
- signal
- tacho
- pulse
- 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
Abstract
基于转速脉冲重构的航空发动机振动信号相位测量方法,主要解决现有技术对振动信号相位测量误差大和不稳定的问题。所述测量方法,通过分析转速信号和振动信号,在本地重构两路正交的参考信号,并通过测角算法估计出振动信号的相位。本发明具有运算量小、稳定性好、且测速精度高的优点,可用于对航空发动机振动信号相位的精确测量。
Description
技术领域
本发明属于航空发动机振动监测技术领域,具体的说是一种航空发动机振动信号相位测量方法,用于在航空发动机高速运行时的复杂振动环境下,稳定并精确的实现对振动信号相位的测量。
背景技术
目前,振动信号分析是诊断航空发动机故障的一种有效方法,绝大多数结构强度方面的故障都与振动信号有密切的关系。发动机振动监测是状态监控与故障诊断的一项重要内容。通过相关的传感器,可以采集到发动机的各种振动信号,通过对振动信号的振幅、振动烈度、相位等参数的实时检测,结合振动信号的固有特征,可实时监测发动机的运转情况,避免大事故发生给企业造成损失。振动信号的相位信息对转子的动平衡、发动机的动态特性和故障特性有着重要的意义,因此研究振动信号相位测量方法具有非常重要的工程意义。
传统的振动信号相位测量方法一般先采用FFT运算将信号转换到频域,然后再实现对相位的测量。这种方法是等时间间隔采样,对转速稳定的稳态信号,其估计的相位较准确。但当转子运行不稳定时,比如启、停车过程,每周期内的采集点数是时变的,此时利用FFT计算非周期采样将会产生由频谱泄露引起的泄露误差,使测得的幅值、相角偏离实际值,尤其是相位测量误差更大。其改进方法是保证等角度采样,从而将非稳态信号通过恒定角增量采样转变为角域的稳定信号,有效地克服了“频谱泄露”现象。在现有方法中,一般通过实时改变采样率或通过对数字信号进行插值来实现等角度采样。实时改变采样率的方法通过转速信息控制实际的采样率实现,它是一个动态的平衡过程,其精度有限;数字信号插值的方法相对比较容易实现,在转速达到动平衡时,其精度较高,但当转速大动态变化时,需要插值的点数变大,会影响插值精度,从而影响对相位的估计精度。
发明内容
本发明的目的在于克服上述振动信号相位测量方法的不足,提出了一种基于转速脉冲重构的航空发动机转子振动信号相位测量方法,以避免由于频谱泄露、信噪比低等诸多因素的影响,提高测量精度。
本发明的技术解决方案是:一种基于转速脉冲重构的航空发动机振动信号相位测量方法,包括如下过程:
(1)监测转速脉冲,根据下式计算当前时刻转速脉冲的周期Tk,
Tk=tk-tk-1
1)
其中,tk为转速脉冲的第k个上升沿对应的时刻,k=1,2,3,...,表示上升沿的序号;
(2)计算振动信号在转速脉冲周期Tk内采集数据的点数,记为Nk,并将该周期内的振动信号记为Sk(t),tk-1<t≤tk;
(3)结合实际系统中的采样率,产生周期为Tk、点数为Nk的两路正交信号,作为参考信号;两路正交信号分别为相位为0°的正弦和余弦信号,分别记为Ik(t)和Qk(t),tk-1<t≤tk;
(4)将Ik(t)和Qk(t)作为两路本振,和振动信号Sk(t)进行混频,并将不同时刻的混频结果进行顺序拼接,得到新的数据I′k(t)和Q′k(t);
(5)通过低通滤波器滤去高频分量,滤波后的信号记为I″(t)和Q″(t);
(6)分别对I″(t)和Q″(t)进行积分,得到两路积分结果;
(7)通过下式估计振动信号的相位
(8)通过α,β滤波或者卡尔曼滤波对振动信号相位进行平滑,得到最终的相位估计结果ψ(t)。
上述步骤4)I′k(t)和Q′k(t)的计算方法是:
I′(t)=[...,Ik-2(t)*Sk-2(t),Ik-1(t)*Sk-1(t),Ik(t)*Sk(t),...], 2.1)
Q′(t)=[...,Qk-2(t)*Sk-2(t),Qk-1(t)*Sk-1(t),Qk(t)*Sk(t),...], 2.2)
其中,*表示混频。
上述步骤5)I″(t)和Q″(t)的计算方法是:
I″(t)=lpf[I′(t)], 3.1)
Q″(t)=lpf[Q′(t)], 3.2)
其中,lpf[x]表示对x进行低通滤波;
上述步骤6)具体积分方法是:
其中,T为积分时间,T=[Tk-L+1,Tk-L+2,...,Tk-1,Tk],L为整数,
当积分结果送给步骤(7)后,自动将积分结果清零,以便进行下一次积分;
上述步骤(3)中所述的参考信号,是通过步骤(1)中估计的周期Tk、步骤(2)中估计的点数Nk和实际系统中的采样率产生的。
上述步骤(4)中所述的混频,是将通过步骤(3)所产生的两路正交信号作为两路本振的。
上述步骤(1)至步骤(8)循环执行,实现对发动机振动信号相位的实时测量。
本发明与现有技术相比具有以下优点:
(1)本发明采用信号重构的方法,通过转速脉冲信号重构出与振动信号具有相同频率和相同采样点数的本地0相位参考信号,为测量振动信号的相位提供了稳定和标准的参考源,使估计的信号相位更加稳定和准确;
(2)相比通过频谱分析来估计相位的方法,本发明采用直接鉴相方法获得参考信号和振动信号之间的相位差,它的本质是每个转速脉冲周期为单位,若信噪比高,它能有效鉴别出每个转速脉冲周期Tk内的震动信号相位,灵敏高;若信噪比低,可通过增加积分时间T来等效降低滤波带宽,提高信噪比,获得高精度的相位估计;
(3)由于本发明重构的本地信号和振动信号频率相同,当振动信号存在较大的高次谐波分量时,对相位估计影响很小,方法的抗干扰能力强;
(4)本发明的运算量很小,信号重构方法简单有效,特别适合实时实现;
附图说明
图1是本发明方法的应用环境框图;
图2是本发明的振动信号相位测量流程框图;
图3是本发明方法中的同相参考信号重构示意图;
图4是本发明方法中的正交参考信号重构示意图;
图5是本发明方法等价后的震动信号相位示意图;
图6是本发明方法中振动信号经过低通滤波后的示意图;
图7是本发明方法中通过积分结果估计的相位示意图;
具体实施方式
本发明提出了一种基于转速脉冲重构的航空发动机振动信号相位测量方法,下面结合附图详细说明本发明的内容和效果。
参照图2,本振动信号相位方法包括如下步骤:
步骤1,监测转速脉冲,根据下式计算当前时刻转速脉冲的周期Tk,
Tk=tk-tk-1 1)
其中,tk为转速脉冲的第k个上升沿对应的时刻,k=1,2,3,...,表示上升沿的序号;
步骤2,计算振动信号在转速脉冲周期Tk内采集数据的点数,记为Nk,并将该周期内的振动信号记为Sk(t),tk-1<t≤tk;
步骤3,结合实际系统中的采样率,产生周期为Tk、点数为Nk的两路正交信号,作为本发明方法中的参考信号。它们分别为相位为0°的正弦和余弦信号,分别记为Ik(t)和Qk(t),tk-1<t≤tk;
图3和图4分别是本发明方法中的同相和正交参考信号重构示意图。如图所示,根据当前转速脉冲周期Tk,采样点数Nk和采样率即可重构出同相和正交参考信号。
图5是本发明方法中需测量的振动信号相位示意图。如图所示,当构造出参考信号后,振动信号相位的求解等价为求解振动信号和同相参考信号之间的相位差。该相位差可通过步骤4~步骤8实现。
步骤4,将Ik(t)和Qk(t)作为两路本振,和振动信号Sk(t)进行混频,并将不同时刻的混频结果进行顺序拼接,得到新的数据I′k(t)和Q′k(t),该过程可通过下式表示,
I′(t)=[...,Ik-2(t)*Sk-2(t),Ik-1(t)*Sk-1(t),Ik(t)*Sk(t),...], 2.1)
Q′(t)=[...,Qk-2(t)*Sk-2(t),Qk-1(t)*Sk-1(t),Qk(t)*Sk(t),...], 2.2)
其中,*表示混频;
步骤5,通过低通滤波器滤去高频分量,滤波后的信号记为I″(t)和Q″(t),该过程可以表示为,
I″(t)=lpf[I′(t)], 3.1)
Q″(t)=lpf[Q′(t)], 3.2)
其中,lpf[x]表示对x进行低通滤波;
图6是本发明方法中振动信号经过低通滤波后的示意图。如图所示,低通滤波后高频分量被滤除,只剩下决定振动信号相位的两路直流信号。
步骤6,按照下式分别对I″(t)和Q″(t)进行积分,得到两路积分结果,
其中,T为积分时间,T=[Tk-L+1,Tk-L+2,...,Tk-1,Tk],L为整数,可根据所需积分时间来设置。当积分结果送给步骤7后,自动将积分结果清零,以便进行下一次积分;
步骤7,通过下式估计振动信号的相位
图7是本发明方法中通过积分结果估计的相位示意图。如图所示,其中的相位通过式5)估计得到。
步骤8,通过α,β滤波或者卡尔曼滤波对振动信号相位进行平滑,得到最终的相位估计结果ψ(t);
步骤9,循环执行步骤1至步骤8,实现对发动机振动信号相位的实时测量。
本发明的效果可以通过以下仿真数据进一步的说明:
采样率设置为1MHz,积分时间为1s,振动信号相位真值设为0°,分别改变振动信号的频率变化率,信噪比和包含的谐波分量用以模拟转子动平衡、转速大动态变化、存在高次谐波分量以及它们混合存在的情况,不同情况下的相位测量结果和误差如表1、表2、表3和表4所示。
表1本发明在振动信号在转子动平衡情况下的处理结果
表2本发明在振动信号的转速大动态变化情况下的处理结果
由表1和表2可见,在转子动平衡和转速大动态变化的情况下,信噪比相同时,测量相位均值和误差标准差变化很小;信噪比不同时,低信噪比下的测量误差标准差较大,均值变化不大。这说明本发明方法的测量精度受振动信号动态的影响很小,受信噪比影响较大。
表3本发明转子动平衡时并存在高次谐波情况下的处理结果
表4本方法在转速大动态变化并存在高次谐波情况下的处理结果
相比表1和表2,表3和表4中添加了振动信号的高次谐波分量。对比相同信噪比下的测量结果发现,测量精度有所下降,这是由高次谐波分量引起的。由于本发明方法在步骤3中产生的参考信号仅仅是振动信号的基频,因此,高次谐波分量对测量结果的影响有限,该方法测量结果的稳定性好。
综上,本发明充分考虑到航空发动机振动信号相位测量的实际应用问题,根据转速信号和振动信号完全耦合的特征,通过转速信号的周期和振动信号的采样点来重构本地参考信号,再以本地参考信号为基准通过角度解算方法测量振动信号的相位,所获得的结果不但精确,而且很稳健。
Claims (7)
1.基于转速脉冲重构的航空发动机振动信号相位测量方法,包括如下过程:
(1)监测转速脉冲,根据下式计算当前时刻转速脉冲的周期Tk,
Tk=tk-tk-1 1)
其中,tk为转速脉冲的第k个上升沿对应的时刻,k=1,2,3,...,表示上升沿的序号;
(2)计算振动信号在转速脉冲周期Tk内采集数据的点数,记为Nk,并将该周期内的振动信号记为Sk(t),tk-1<t≤tk;
(3)结合实际系统中的采样率,产生周期为Tk、点数为Nk的两路正交信号,作为参考信号;两路正交信号分别为相位为0°的正弦和余弦信号,分别记为Ik(t)和Qk(t),tk-1<t≤tk;
(4)将Ik(t)和Qk(t)作为两路本振,和振动信号Sk(t)进行混频,并将不同时刻的混频结果进行顺序拼接,得到新的数据I′k(t)和Q′k(t);
(5)通过低通滤波器滤去高频分量,滤波后的信号记为I″(t)和Q″(t);
(6)分别对I″(t)和Q″(t)进行积分,得到两路积分结果;
(7)通过下式估计振动信号的相位
(8)通过α,β滤波或者卡尔曼滤波对振动信号相位进行平滑,得到最终的相位估计结果ψ(t)。
2.根据权利要求1所述的基于转速脉冲重构的航空发动机振动信号相位测量方法,其特征在于:所述步骤4)I′k(t)和Q′k(t)的计算方法是:
I′(t)=[...,Ik-2(t)*Sk-2(t),Ik-1(t)*Sk-1(t),Ik(t)*Sk(t),...], 2.1)
Q′(t)=[...,Qk-2(t)*Sk-2(t),Qk-1(t)*Sk-1(t),Qk(t)*Sk(t),...], 2.2)
其中,*表示混频。
3.根据权利要求1所述的基于转速脉冲重构的航空发动机振动信号相位测量方法,其特征在于:所述步骤5)I″(t)和Q″(t)的计算方法是:
I″(t)=lpf[I′(t)], 3.1)
Q″(t)=lpf[Q′(t)], 3.2)
其中,lpf[x]表示对x进行低通滤波。
4.根据权利要求1所述的基于转速脉冲重构的航空发动机振动信号相位测量方法,其特征在于:所述步骤6)具体积分方法是:
其中,T为积分时间,T=[Tk-L+1,Tk-L+2,...,Tk-1,Tk],L为整数,
当积分结果送给步骤(7)后,自动将积分结果清零,以便进行下一次积分。
5.根据权利要求1或2或3或4所述的基于转速脉冲重构的航空发动机振动信号相位测量方法,其特征在于:步骤(3)中所述的参考信号,是通过步骤(1)中估计的周期Tk、步骤(2)中估计的点数Nk和实际系统中的采样率产生的。
6.根据权利要求5所述的基于转速脉冲重构的航空发动机振动信号相位测量方法,其特征在于:步骤(4)中所述的混频,是将通过步骤(3)所产生的两路正交信号作为两路本振的。
7.根据权利要求6所述的基于转速脉冲重构的航空发动机振动信号相位测量方法,其特征在于:所述步骤(1)至步骤(8)循环执行,实现对发动机振动信号相位的实时测量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310625488.2A CN104677486B (zh) | 2013-11-27 | 2013-11-27 | 基于转速脉冲重构的航空发动机振动信号相位测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310625488.2A CN104677486B (zh) | 2013-11-27 | 2013-11-27 | 基于转速脉冲重构的航空发动机振动信号相位测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104677486A true CN104677486A (zh) | 2015-06-03 |
CN104677486B CN104677486B (zh) | 2017-10-24 |
Family
ID=53312840
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310625488.2A Active CN104677486B (zh) | 2013-11-27 | 2013-11-27 | 基于转速脉冲重构的航空发动机振动信号相位测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104677486B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105466550A (zh) * | 2015-12-04 | 2016-04-06 | 中国人民解放军国防科学技术大学 | 非均匀欠采样叶端定时振动信号重构方法及其装置 |
CN107449932A (zh) * | 2017-06-02 | 2017-12-08 | 中国航空规划设计研究总院有限公司 | 一种航空发动机转子转速的测量方法 |
CN109668735A (zh) * | 2018-12-11 | 2019-04-23 | 中国航空工业集团公司西安航空计算技术研究所 | 一种发动机转子相位基准确定方法、装置和电路 |
CN112746875A (zh) * | 2019-10-31 | 2021-05-04 | 中国航发商用航空发动机有限责任公司 | 航空发动机转子轴系复杂振动的主动控制系统、方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS57179625A (en) * | 1981-04-30 | 1982-11-05 | Hitachi Ltd | Method for diagnosing vibration in rotary machine |
US4607529A (en) * | 1984-03-21 | 1986-08-26 | John Morey | Vibration analysis |
CN101750198A (zh) * | 2008-11-28 | 2010-06-23 | 上海宝钢工业检测公司 | 旋转机械不同测点振动信号相对相位的测量方法 |
CN102565806A (zh) * | 2011-12-31 | 2012-07-11 | 北京握奇数据系统有限公司 | 一种激光测距方法及装置 |
CN103226200A (zh) * | 2013-03-22 | 2013-07-31 | 浙江理工大学 | 一种手持式激光测距装置 |
CN103234627A (zh) * | 2013-04-17 | 2013-08-07 | 国家电网公司 | 旋转机械振动信号整周期同步采样分析方法 |
-
2013
- 2013-11-27 CN CN201310625488.2A patent/CN104677486B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS57179625A (en) * | 1981-04-30 | 1982-11-05 | Hitachi Ltd | Method for diagnosing vibration in rotary machine |
US4607529A (en) * | 1984-03-21 | 1986-08-26 | John Morey | Vibration analysis |
CN101750198A (zh) * | 2008-11-28 | 2010-06-23 | 上海宝钢工业检测公司 | 旋转机械不同测点振动信号相对相位的测量方法 |
CN102565806A (zh) * | 2011-12-31 | 2012-07-11 | 北京握奇数据系统有限公司 | 一种激光测距方法及装置 |
CN103226200A (zh) * | 2013-03-22 | 2013-07-31 | 浙江理工大学 | 一种手持式激光测距装置 |
CN103234627A (zh) * | 2013-04-17 | 2013-08-07 | 国家电网公司 | 旋转机械振动信号整周期同步采样分析方法 |
Non-Patent Citations (4)
Title |
---|
周丽芹等: "一种转速信号测量和相位信息获取的方法", 《振动、测试与诊断》 * |
杨小牛等: "《软件无线电技术与应用》", 30 April 2010, 北京理工大学出版社 * |
王四季等: "发动机高速动平衡的振动相位分析", 《燃气涡轮试验与研究》 * |
闫磊: "正交混频及数字解调的激光测距系统研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105466550A (zh) * | 2015-12-04 | 2016-04-06 | 中国人民解放军国防科学技术大学 | 非均匀欠采样叶端定时振动信号重构方法及其装置 |
CN105466550B (zh) * | 2015-12-04 | 2018-08-28 | 中国人民解放军国防科学技术大学 | 非均匀欠采样叶端定时振动信号重构方法及其装置 |
CN107449932A (zh) * | 2017-06-02 | 2017-12-08 | 中国航空规划设计研究总院有限公司 | 一种航空发动机转子转速的测量方法 |
CN109668735A (zh) * | 2018-12-11 | 2019-04-23 | 中国航空工业集团公司西安航空计算技术研究所 | 一种发动机转子相位基准确定方法、装置和电路 |
CN112746875A (zh) * | 2019-10-31 | 2021-05-04 | 中国航发商用航空发动机有限责任公司 | 航空发动机转子轴系复杂振动的主动控制系统、方法 |
CN112746875B (zh) * | 2019-10-31 | 2022-08-19 | 中国航发商用航空发动机有限责任公司 | 航空发动机转子轴系复杂振动的主动控制系统、方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104677486B (zh) | 2017-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108871742B (zh) | 一种改进的无键相故障特征阶次提取方法 | |
Wang et al. | Rolling element bearing fault diagnosis via fault characteristic order (FCO) analysis | |
Yang et al. | Vibration signal analysis using parameterized time–frequency method for features extraction of varying-speed rotary machinery | |
CN105738696B (zh) | 全相位时移相位差频率估计方法及装置 | |
Brandt et al. | Main principles and limitations of current order tracking methods | |
Rodopoulos et al. | A parametric approach for the estimation of the instantaneous speed of rotating machinery | |
CN102353500B (zh) | 一种用于动平衡测量的不平衡信号提取方法 | |
CN101509945B (zh) | 正负序电量实时检测的方法 | |
CN104697436A (zh) | 一种基于傅里叶级数的圆感应同步器误差模型分析方法 | |
CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
CN104677486A (zh) | 基于转速脉冲重构的航空发动机振动信号相位测量方法 | |
CN102221639A (zh) | 正负序电流实时检测的方法 | |
CN108761117B (zh) | 便携式电流检测转速测试仪 | |
Tang et al. | Fractional frequency band entropy for bearing fault diagnosis under varying speed conditions | |
CN114061678B (zh) | 一种科氏流量计数字驱动方法 | |
CN102721462B (zh) | 旋转机械启停车过程波德图/奈奎斯特图的快速计算方法 | |
Chen et al. | Nonstationary signal denoising using an envelope-tracking filter | |
CN110736926A (zh) | 电机运行状态的特征参数提取方法和装置 | |
CN103018555A (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN104502707A (zh) | 一种基于三次样条插值的电力系统同步相量测量方法 | |
CN108398260B (zh) | 基于混合概率方法的齿轮箱瞬时角速度的快速评估方法 | |
Lv et al. | Generalized synchroextracting-based stepwise demodulation transform and its application to fault diagnosis of rotating machinery | |
CN105571874A (zh) | 一种发动机振动不平衡相位的实时测量方法 | |
CN107247155A (zh) | 一种汽轮发电机组转子瞬时转速测试装置及其测试方法 | |
CN114486252B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |