一种应用于直升机传动系统振动信号的故障检测方法技术方案

本发明专利技术公开了一种应用于直升机传动系统振动信号的故障检测方法，包含1，对采集到的原始直升机传动系统振动信号进行降噪处理，得到信号x(t)；2，运用EMD算法对信号x(t)进行分解，得到内稟模态函数分量；3，求取内稟模态函数分量的近似熵，得到内稟模态函数分量的近似熵值；4，将内稟模态函数分量的近似熵值进行标量量化处理，最终得到标量量化的近似熵值；5，将标量量化的近似熵值带入模糊隶属度函数中，求取其隶属度函数值；6，通过求取隶属度函数值的最大值，以得到隶属于某一故障的隶属度程度，最终确定故障检测结果。通过本发明专利技术，能有效提高直升机传动系统的安全性能。

【技术实现步骤摘要】
一种应用于直升机传动系统振动信号的故障检测方法
本专利技术属于故障诊断方法应用领域，具体涉及一种应用于直升机传动系统振动信号故障检测新方法。
技术介绍
作为直升机上重要的动力传动装置——直升机传动系统是直升机上的典型机械旋转系统，是由大量相互联系但工作过程又彼此不同的部件组成，存在着强烈的非线性、非平稳性、不确定性等复杂系统特征，导致了其易发故障，且一旦发生故障，将导致灾难性事故的发生。由于直升机传动系统的工作环境恶劣，工作强度大等特点，导致其易发故障，且故障难以准确诊断，所以对直升机传动系统故障进行准确诊断一直是提高航空武器装备安全运行的一个无法逾越的技术瓶颈，实现直升机传动系统故障的准确、快速、可靠检测是保证直升机飞行安全和提高直升机出勤率重要手段之一。然而现有的多数直升机传动系统故障检测方法都无法得出准确可靠的故障检测结果，且普遍存在虚警率高、诊断精度低的缺点。因此，亟需研究直升机传动系统故障检测技术，提高直升机传动系统的安全性能。
技术实现思路
本专利技术的专利技术目的在于提供一种应用于直升机传动系统振动信号的故障检测方法，使直升机传动系统故障检测能力更可靠、更准确性。本专利技术的专利技术目的通过以下技术方案实现：一种应用于直升机传动系统振动信号的故障检测方法，包含以下步骤：步骤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，重复步骤3.1～步骤3.4，得到φm+1(r)；步骤3.6，得出的序列长度为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个，无法形成包络线本文档来自技高网...

【技术保护点】
一种应用于直升机传动系统振动信号的故障检测方法，包含以下步骤：步骤1，对采集到的原始直升机传动系统振动信号进行降噪处理，得到信号x(t)；步骤2，运用EMD算法对信号x(t)进行分解，得到内稟模态函数分量，记为u(i){i＝1,2，…，k}；步骤3，求取内稟模态函数分量的近似熵，得到内稟模态函数分量的近似熵值；步骤4，将内稟模态函数分量的近似熵值进行标量量化处理，最终得到标量量化的近似熵值；步骤5，将标量量化的近似熵值带入模糊隶属度函数中，求取其隶属度函数值；步骤6，通过求取隶属度函数值的最大值，以得到隶属于某一故障的隶属度程度，最终确定故障检测结果。

【技术特征摘要】
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)与数据序...

【专利技术属性】
技术研发人员：林泽力，何泳，沈勇，郑国，王景霖，单添敏，马振国，曹亮，何召华，
申请(专利权)人：中国航空工业集团公司上海航空测控技术研究所，
类型：发明
国别省市：上海;31