CN105758634A - 一种应用于直升机传动***振动信号的故障检测方法 - Google Patents
一种应用于直升机传动***振动信号的故障检测方法 Download PDFInfo
- Publication number
- CN105758634A CN105758634A CN201610105458.2A CN201610105458A CN105758634A CN 105758634 A CN105758634 A CN 105758634A CN 201610105458 A CN201610105458 A CN 201610105458A CN 105758634 A CN105758634 A CN 105758634A
- Authority
- CN
- China
- Prior art keywords
- signal
- approximate entropy
- fault detection
- mode function
- detection method
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种应用于直升机传动***振动信号的故障检测方法,包含1,对采集到的原始直升机传动***振动信号进行降噪处理,得到信号x(t);2,运用EMD算法对信号x(t)进行分解,得到内稟模态函数分量;3,求取内稟模态函数分量的近似熵,得到内稟模态函数分量的近似熵值;4,将内稟模态函数分量的近似熵值进行标量量化处理,最终得到标量量化的近似熵值;5,将标量量化的近似熵值带入模糊隶属度函数中,求取其隶属度函数值;6,通过求取隶属度函数值的最大值,以得到隶属于某一故障的隶属度程度,最终确定故障检测结果。通过本发明,能有效提高直升机传动***的安全性能。
Description
技术领域
本发明属于故障诊断方法应用领域,具体涉及一种应用于直升机传动***振动信号故障检测新方法。
背景技术
作为直升机上重要的动力传动装置——直升机传动***是直升机上的典型机械旋转***,是由大量相互联系但工作过程又彼此不同的部件组成,存在着强烈的非线性、非平稳性、不确定性等复杂***特征,导致了其易发故障,且一旦发生故障,将导致灾难性事故的发生。由于直升机传动***的工作环境恶劣,工作强度大等特点,导致其易发故障,且故障难以准确诊断,所以对直升机传动***故障进行准确诊断一直是提高航空武器装备安全运行的一个无法逾越的技术瓶颈,实现直升机传动***故障的准确、快速、可靠检测是保证直升机飞行安全和提高直升机出勤率重要手段之一。然而现有的多数直升机传动***故障检测方法都无法得出准确可靠的故障检测结果,且普遍存在虚警率高、诊断精度低的缺点。因此,亟需研究直升机传动***故障检测技术,提高直升机传动***的安全性能。
发明内容
本发明的发明目的在于提供一种应用于直升机传动***振动信号的故障检测方法,使直升机传动***故障检测能力更可靠、更准确性。
本发明的发明目的通过以下技术方案实现:
一种应用于直升机传动***振动信号的故障检测方法,包含以下步骤:
步骤1,对采集到的原始直升机传动***振动信号进行降噪处理,得到信号x(t);
步骤2,运用EMD算法对信号x(t)进行分解,得到内稟模态函数分量,记为u(i){i=1,2,…,k};
步骤3,求取内稟模态函数分量的近似熵,得到内稟模态函数分量的近似熵值;
步骤4,将内稟模态函数分量的近似熵值进行标量量化处理,最终得到标量量化的近似熵值;
步骤5,将标量量化的近似熵值带入模糊隶属度函数中,求取其隶属度函数值;
步骤6,通过求取隶属度函数值的最大值,以得到隶属于某一故障的隶属度程度,最终确定故障检测结果。
优选地,所述步骤1中采用缺省阈值降噪技术对原始直升机传动***振动信号进行降噪处理。
优选地,所述步骤2包含以下步骤:
步骤2.1,找出信号x(t)所有的极大值点,并将极大值点用三次样条函数插值成为上包络线v1(t);找出信号x(t)所有的极小值点,并将极小值点用三次样条函数插值成为下包络线v2(t),求出上包络线和下包络线的均值m1(t):
信号x(t)减去m1(t)得到去除低频的数据序列h1(t):
h1(t)=x(t)-m1(t)
步骤2.2,将数据序列h1(t)代替信号x(t),重复执行步骤2.1,得到数据序列h1n(t):
h1n(t)=h1(n-1)(t)-m1n(t),
其中,n代表步骤2.2重复执行的次数;
步骤2.3,当数据序列h1n(t)与数据序列h1(n-1)(t)的标准差SD低于0.3时,得到内稟模态函数c1(t)=h1n(t);
步骤2.4,将信号x(t)去除内稟模态函数c1(t)得到的剩余序列r1(t),再用剩余序列r1(t)代替信号x(t),重复步骤2.1-步骤2.3,直到得到的剩余序列rk(t)不再可分解,得到k个内稟模态函数,记为u(i){i=1,2,…,k}。
优选地,所述剩余序列rk(t)不再可分解的判断条件为以下任意一种:
①余项rk(t)足够小,接近于0;
②得到的rk(t)成单调递增或单调递减函数或极值点少于2个,无法形成包络线,无法分解得到IMF分量。
优选地,所述步骤3中,近似熵值通过以下计算过程得到:
步骤3.1,将序列u(i)按顺序组成m维矢量X(i),即
步骤3.2,对每一个i值计算矢量X(i)与其余矢量X(j)之间的距离:
式中:d[X(i),X(j)]为矢量X(i)与X(j)间的距离;
步骤3.3,按照设定的阈值r(r>0),对每一个i值统计d[X(i),X(j)]<r的数目与总的矢量个数k-m+1的比值,记做即
步骤3.4,先将取对数,再求其对所有i的平均值,记做
步骤3.5,重复过程(1)~(4),得到φm+1(r);
步骤3.5,得出的序列长度为K时的近似熵值的估计值:
ApEn(m.r,K)=φm(r)-φm+1(r)。
附图说明
下面结合附图和具体实施方式对本发明做进一步说明。
图1为本发明一种应用于直升机传动***振动信号的故障检测方法的流程示意图。
具体实施方式
如图1所示,为本发明一种应用于直升机传动***振动信号的故障检测方法,主要是依据采集到的直升机传动***非平稳振动信号,采用缺省阈值降噪技术,对采集到的非平稳振动信号进行数据降噪处理,得到本发明专利所需要的原始数据;然后,运用EMD方法(经验模态分解,EmpiricalModeDecomposition)对降噪过的非平稳振动信号进行分解,获取其内稟模态函数分量(IMF),从而使得将振动信号成分中的特征信息成分分解出来,利用获得的IMF分量,计算其近似熵值,最终可利用该近似熵值对直升机传动***进行故障检测;根据所获得的近似熵值,运用模糊理论算法可实现对直升机传动***故障发生情况进行准确有效地检测。具体包含以下步骤:
步骤1,根据缺省阈值降噪技术的原理,对采集到的直升机传动***振动信号进行降噪处理信号,得到信号x(t);
步骤2,运用EMD算法对信号x(t)进行分解,得到内稟模态函数分量,记为u(i){i=1,2,…,k}。主要分为以下几个实现过程:
过程(1):将步骤1所得到的降噪处理的过的信号x(t),找出信号x(t)所有的极大值点并将其用三次样条函数插值成为信号x(t)的上包络线v1(t):找出信号x(t)所有的极小值点并将其用三次样条函数插值成为信号x(t)的下包络线v2(t),求出上包络线和下包络线的均值m1(t):
信号x(t)减去m1(t)得到去除低频的新数据序列h1(t):
h1(t)=x(t)-m1(t)(2)
采用与计算信号x(t)的均值m1(t)相同的计算方法,计算h1(t)的包络均值m11(t),即找出数据序列h1(t)所有的极大值点并将其用三次样条函数插值成为数据序列h1(t)的上包络线v11(t):找出数据序列h1(t)所有的极小值点并将其用三次样条函数插值成为数据序列h1(t)的下包络线v21(t),
数据序列h1(t)减去新的包络均值m11(t)后得到更新的数据序列h11(t):
h11(t)=h1(t)-m11(t)(4)
为了得到振动数据序列频率最高的分量,需要不断地重复上述过程以得到新的数据序列h1n(t):
h1n(t)=h1(n-1)(t)-m1n(t)(5)
一般情况下,停止筛选的条件主要有3种:
①数据序列的过零点和极值点数量相等;
②符合基于幅值比的终止条件;
③h1n(t)与h1(n-1)(t)的标准差SD低于0.3。
本发明采用第三种终止条件,标准差的计算公式为:
最后得到第一阶,也就是最高频率成分的本征模函数c1(t)。
过程(2):信号x(t)去除内稟模态函数c1(t)得到的剩余信号r1(t),再对剩余信号r1(t)重复过程(1)的计算,直到最后得到的剩余序列rk(t)不再可分解,可以得到k个内稟模态函数。
其中:在上述EMD分解过程中,在过程(2)中需要确定剩余数据序列是否不可分解,一般认为当数据序列满足一下任意一个条件,就可以停止分解:
①余项rk(t)足够小,接近于0;
②得到的rk(t)成单调递增或单调递减函数或极值点少于2个,无法形成包络线,无法分解得到IMF分量。
步骤3,将步骤2得到的IMF分量记为u(i){i=1,2,…,k},预先给定模式维数m和相似容限r的值,则近似熵可以通过以下计算过程得到:
过程(1):将序列u(i)按顺序组成m维矢量X(i),即
过程(2):对每一个i值计算矢量X(i)与其余矢量X(j)之间的距离:
式中:d[X(i),X(j)]为矢量X(i)与X(j)间的距离。
过程(3):按照设定的阈值r(r>0),对每一个i值统计d[X(i),X(j)]<r的数目与总的矢量个数k-m+1的比值,记做即
过程(4):先将取对数,再求其对所有i的平均值,记做
过程(5):重复过程(1)~(4),得到φm+1(r)。
过程(6):理论上此序列的近似熵值为
一般而言,此极限值以概率1存在,但在实际工作中K不可能为∞,当K为有限值时,按上述过程得出的序列长度为K时的ApEn的估计值,记作:
ApEn(m.r,K)=φm(r)-φm+1(r)(11)
ApEn的值显然与m、r的取值有关。根据经验,通常取m=2,r=0.1~0.25SD(u)(SD表示序列{u(i){i=1,2,…,k}的标准差}),此时近似熵具有较为合理的统计特性。
步骤4,将步骤3得到的IMF分量的近似熵值进行标量量化处理,最终可得到标量量化的近似熵值;
步骤5,根据得到的量化近似熵值,将其带入模糊隶属度函数中,求取其隶属度函数值;
步骤6,根据步骤5得到的隶属度函数值,通过求取隶属度函数值的最大值,以得到其隶属于某一故障的隶属度程度,最终确定故障检测结果。
可以理解的是,对本领域普通技术人员来说,可以根据本发明的技术方案及其发明构思加以等同替换或改变,而所有这些改变或替换都应属于本发明所附的权利要求的保护范围。
Claims (5)
1.一种应用于直升机传动***振动信号的故障检测方法,包含以下步骤:
步骤1,对采集到的原始直升机传动***振动信号进行降噪处理,得到信号x(t);
步骤2,运用EMD算法对信号x(t)进行分解,得到内稟模态函数分量,记为u(i){i=1,2,…,k};
步骤3,求取内稟模态函数分量的近似熵,得到内稟模态函数分量的近似熵值;
步骤4,将内稟模态函数分量的近似熵值进行标量量化处理,最终得到标量量化的近似熵值;
步骤5,将标量量化的近似熵值带入模糊隶属度函数中,求取其隶属度函数值;
步骤6,通过求取隶属度函数值的最大值,以得到隶属于某一故障的隶属度程度,最终确定故障检测结果。
2.根据权利要求1所述的故障检测方法,其特征在于所述步骤1中采用缺省阈值降噪技术对原始直升机传动***振动信号进行降噪处理。
3.根据权利要求1所述的故障检测方法,其特征在于所述步骤2包含以下步骤:
步骤2.1,找出信号x(t)所有的极大值点,并将极大值点用三次样条函数插值成为上包络线v1(t);找出信号x(t)所有的极小值点,并将极小值点用三次样条函数插值成为下包络线v2(t),求出上包络线和下包络线的均值m1(t):
信号x(t)减去m1(t)得到去除低频的数据序列h1(t):
h1(t)=x(t)-m1(t)
步骤2.2,将数据序列h1(t)代替信号x(t),重复执行步骤2.1,得到数据序列h1n(t):
h1n(t)=h1(n-1)(t)-m1n(t),
其中,n代表步骤2.2重复执行的次数;
步骤2.3,当数据序列h1n(t)与数据序列h1(n-1)(t)的标准差SD低于0.3时,得到内稟模态函数c1(t)=h1n(t);
步骤2.4,将信号x(t)去除内稟模态函数c1(t)得到的剩余序列r1(t),再用剩余序列r1(t)代替信号x(t),重复步骤2.1-步骤2.3,直到得到的剩余序列rk(t)不再可分解,得到k个内稟模态函数,记为u(i){i=1,2,…,k}。
4.根据权利要求3所述的故障检测方法,其特征在于所述剩余序列rk(t)不再可分解的判断条件为以下任意一种:
①余项rk(t)足够小,接近于0;
②得到的rk(t)成单调递增或单调递减函数或极值点少于2个,无法形成包络线,无法分解得到IMF分量。
5.根据权利要求1所述的故障检测方法,其特征在于所述步骤3中,近似熵值通过以下计算过程得到:
步骤3.1,将序列u(i)按顺序组成m维矢量X(i),即
X(i)=[u(i),u(i+1),...,u(i+m-1)]
i=1~k-m+1;
步骤3.2,对每一个i值计算矢量X(i)与其余矢量X(j)之间的距离:
式中:d[X(i),X(j)]为矢量X(i)与X(j)间的距离;
步骤3.3,按照设定的阈值r(r>0),对每一个i值统计d[X(i),X(j)]<r的数目与总的矢量个数k-m+1的比值,记做即
步骤3.4,先将取对数,再求其对所有i的平均值,记做
步骤3.5,重复过程(1)~(4),得到φm+1(r);
步骤3.5,得出的序列长度为K时的近似熵值的估计值:
ApEn(m.r,K)=φm(r)-φm+1(r)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610105458.2A CN105758634B (zh) | 2016-02-26 | 2016-02-26 | 一种应用于直升机传动***振动信号的故障检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610105458.2A CN105758634B (zh) | 2016-02-26 | 2016-02-26 | 一种应用于直升机传动***振动信号的故障检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105758634A true CN105758634A (zh) | 2016-07-13 |
CN105758634B CN105758634B (zh) | 2019-06-21 |
Family
ID=56331286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610105458.2A Active CN105758634B (zh) | 2016-02-26 | 2016-02-26 | 一种应用于直升机传动***振动信号的故障检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105758634B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109000926A (zh) * | 2018-10-25 | 2018-12-14 | 中国计量大学 | 基于emd和近似熵的滚动轴承声发射信号特征提取方法 |
CN109340586A (zh) * | 2018-11-06 | 2019-02-15 | 三川智慧科技股份有限公司 | 一种供水管道泄露的检测方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644623A (zh) * | 2009-06-19 | 2010-02-10 | 湖南大学 | 基于多尺度形态学分析的齿轮故障诊断方法 |
US20120130682A1 (en) * | 2010-11-18 | 2012-05-24 | Jia Wang | Real-time detection system and the method thereof |
CN103234748A (zh) * | 2013-04-02 | 2013-08-07 | 北京工业大学 | 一种基于敏感imf的克林根贝尔格锥齿轮故障诊断方法 |
CN104316317A (zh) * | 2014-10-08 | 2015-01-28 | 西北工业大学 | 一种基于com组件的齿轮***多故障诊断方法 |
CN104748961A (zh) * | 2015-03-30 | 2015-07-01 | 中国矿业大学 | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 |
-
2016
- 2016-02-26 CN CN201610105458.2A patent/CN105758634B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644623A (zh) * | 2009-06-19 | 2010-02-10 | 湖南大学 | 基于多尺度形态学分析的齿轮故障诊断方法 |
US20120130682A1 (en) * | 2010-11-18 | 2012-05-24 | Jia Wang | Real-time detection system and the method thereof |
CN103234748A (zh) * | 2013-04-02 | 2013-08-07 | 北京工业大学 | 一种基于敏感imf的克林根贝尔格锥齿轮故障诊断方法 |
CN104316317A (zh) * | 2014-10-08 | 2015-01-28 | 西北工业大学 | 一种基于com组件的齿轮***多故障诊断方法 |
CN104748961A (zh) * | 2015-03-30 | 2015-07-01 | 中国矿业大学 | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109000926A (zh) * | 2018-10-25 | 2018-12-14 | 中国计量大学 | 基于emd和近似熵的滚动轴承声发射信号特征提取方法 |
CN109340586A (zh) * | 2018-11-06 | 2019-02-15 | 三川智慧科技股份有限公司 | 一种供水管道泄露的检测方法及*** |
CN109340586B (zh) * | 2018-11-06 | 2020-07-28 | 三川智慧科技股份有限公司 | 一种供水管道泄露的检测方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN105758634B (zh) | 2019-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109975013B (zh) | 基于ivmd-se的风电机组齿轮箱故障特征提取方法 | |
Zou et al. | Bearing fault diagnosis based on combined multi-scale weighted entropy morphological filtering and bi-LSTM | |
CN103968937B (zh) | 一种基于emd样本熵和fcm的配电开关机械状态诊断方法 | |
CN111397896B (zh) | 一种旋转机械故障诊断方法、***及存储介质 | |
CN105094111B (zh) | 基于联合降噪和经验模态分解的控制***健康状态分析方法 | |
CN102721545A (zh) | 一种基于多特征参量的滚动轴承故障诊断方法 | |
CN104897403A (zh) | 一种基于排列熵和流形改进动态时间规整的自适应故障诊断方法 | |
CN104316317A (zh) | 一种基于com组件的齿轮***多故障诊断方法 | |
CN107451004A (zh) | 一种基于定性趋势分析的道岔故障诊断方法 | |
Li et al. | A new nonlinear blind source separation method with chaos indicators for decoupling diagnosis of hybrid failures: A marine propulsion gearbox case with a large speed variation | |
CN110653661A (zh) | 基于信号融合和多重分形谱算法的刀具状态监测识别方法 | |
EP3702238B1 (en) | Diagnostic system and a method of diagnosing faults | |
CN104596780B (zh) | 一种动车组制动***传感器故障的诊断方法 | |
CN110779724B (zh) | 一种基于频域组稀疏降噪的轴承故障诊断方法 | |
CN109238455B (zh) | 一种基于图论的旋转机械振动信号监测方法及*** | |
CN111678699B (zh) | 一种面向滚动轴承早期故障监测与诊断方法及*** | |
Liu et al. | A dependence-based feature vector and its application on planetary gearbox fault classification | |
CN105043766A (zh) | 一种基于Hessian局部线性嵌入的轴承变工况故障诊断方法 | |
CN109000926A (zh) | 基于emd和近似熵的滚动轴承声发射信号特征提取方法 | |
CN105758634A (zh) | 一种应用于直升机传动***振动信号的故障检测方法 | |
CN107831012A (zh) | 一种基于Walsh变换与Teager能量算子的轴承故障诊断方法 | |
Han et al. | An intelligent fault diagnosis method of variable condition gearbox based on improved DBN combined with WPEE and MPE | |
Liu et al. | Application of EMD-WVD and particle filter for gearbox fault feature extraction and remaining useful life prediction | |
CN114964769A (zh) | 一种风电齿轮箱振动信号故障诊断方法 | |
CN105718733A (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 |