CN105954030A - 一种基于内秉时间尺度分解和谱峭度的包络分析方法 - Google Patents
一种基于内秉时间尺度分解和谱峭度的包络分析方法 Download PDFInfo
- Publication number
- CN105954030A CN105954030A CN201610492062.8A CN201610492062A CN105954030A CN 105954030 A CN105954030 A CN 105954030A CN 201610492062 A CN201610492062 A CN 201610492062A CN 105954030 A CN105954030 A CN 105954030A
- Authority
- CN
- China
- Prior art keywords
- signal
- decomposition
- data
- envelope
- frequency
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 63
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 51
- 230000003595 spectral effect Effects 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 45
- 238000001228 spectrum Methods 0.000 claims abstract description 20
- 238000001914 filtration Methods 0.000 claims abstract description 18
- 238000006467 substitution reaction Methods 0.000 claims abstract description 16
- 230000008707 rearrangement Effects 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 12
- 238000009499 grossing Methods 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000036039 immunity Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 9
- 238000005096 rolling process Methods 0.000 description 9
- 238000002474 experimental method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010892 electric spark Methods 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
-
- 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/04—Bearings
- G01M13/045—Acoustic or vibration analysis
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于内秉时间尺度分解和谱峭度的包络分析方法,该方法首先利用内秉时间尺度分解方法对原始信号进行分解,然后利用数据的重排和替代操作排除分解结果中的噪声分量和趋势项,接着再采用谱峭度方法对第一次滤波后的信号进行分析,得到最优滤波器的中心频率和带宽,然后利用该滤波器对第一次滤波后的信号再进行第二次滤波,然后采用三次样条迭代平滑包络分析方法对第二次滤波后的信号进行包络分析,最后根据包络谱确定旋转机械的故障类型。本发明适合于处理复杂的旋转机械故障信号,能够准确地判定出旋转机械的故障类型,具有良好的抗噪性和鲁棒性,便于工程应用。
Description
技术领域
本发明涉及旋转机械状态监测与故障诊断领域,具体涉及一种基于内秉时间尺度分解和谱峭度的包络分析方法。
背景技术
包络分析技术广泛应用于齿轮和滚动轴承的故障诊断中。现有的包络分析技术有下面三个缺陷:①现有的包络分析技术或者是直接对原始信号进行分析,或者是仅对原始信号进行简单的滤波后再进行分析,因此现有的方法容易受到噪声、趋势及其它成分的干扰,从而导致现有技术的分析精度较低;②现有的包络分析技术是以Hilbert变换为基础,而Hilbert变换要求被分析的信号必须是单分量的窄带信号,否则信号的频率调制部分将要污染信号的幅值包络分析结果,但是目前待分析的信号都不严格满足单分量且窄带的条件,这样就会导致现有技术因精度不高而容易出现误判问题;③由传统方法得到的包络谱存在着端点效应。
发明内容
本发明要解决的问题是针对以上不足,提出一种基于内秉时间尺度分解和谱峭度的包络分析方法,采用本发明的包络分析方法后,具有分析结果准确度和精确度高,并能准确地检测出旋转机械故障类型的优点。
为解决以上技术问题,本发明采取的技术方案如下:一种基于内秉时间尺度分解和谱峭度的包络分析方法,其特征在于,包括以下步骤:
步骤1:利用加速度传感器以采样频率fs测取旋转机械的振动信号x(k), (k=1, 2,…,N),N为采样信号的长度;
步骤2:采用内秉时间尺度分解算法将信号x(k)分解成n个分量和一个趋势项之和,即,其中,ci(k)代表由内秉时间尺度分解算法得到的第i个分量,rn(k)代表由内秉时间尺度分解算法得到的趋势项;
步骤3:对ci(k)执行重排操作和替代操作,经重排操作得到的数据用ci shuffle(k)表示,替代操作后得到数据用ci FTran(k)表示;
步骤4:对ci(k)、ci shuffle(k)和ci FTran(k)分别执行多重分形去趋势波动分析(Multifractal Detrended Fluctuation Analysis, MFDFA),得到广义Hurst指数曲线,ci(k)的广义Hurst指数曲线用Hi(q)表示;ci shuffle(k)的广义Hurst指数曲线用Hi shuffle(q)表示;ci FTran(k)的广义Hurst指数曲线用Hi FTran(q)表示;
步骤5:如果Hi(q) 与Hi shuffle(q)或Hi(q) 与Hi FTran(q)之间的相对误差小于5%,或者Hi(q) 、Hi shuffle(q) 和Hi FTran(q)三者都不随q而变化,则抛弃对应的ci(k)分量;
步骤6:对剩余的ci(k)分量求和,将该和记为信号经重排和替代滤波后的结果xf1(k);
步骤7:对xf1(k)执行谱峭度分析,求出信号峭度最大处所对应的中心频率f0和带宽B;
步骤8: 根据中心频率f0和带宽B对xf1(k)进行带通滤波,得到xf2(k);
步骤9:对信号xf2(k)执行三次样条迭代平滑包络分析,得到信号包络eov(k);
步骤10:对得到的信号包络eov(k)执行离散傅里叶变换得到包络谱,根据包络谱特征频率判断机器的故障类型。
一种优化方案,所述步骤2中内秉时间尺度分解算法包括以下步骤:
1) 对于任意信号xt,(t=1, 2, …,N),定义一个算子用于抽取低频基线信号,即:
其中是基线信号,是一个固有旋转分量, 假设是一个实值信号,代表xt的局部极值所对应的时刻,为方便起见定义;如果xt在某个区间上具有恒定值,考虑到邻近的信号存在着波动,我们仍然认为xt在这个区间上包含着极值,这时是该区间的右端点;为方便起见,定义,;假设和在上有定义,xk在上有定义,在区间上的连续极值点之间定义一个分段线性的基线信号抽取算子,即:
其中
这里参数是一个线性增益,,本例中;
2) 定义一个固有旋转分量抽取算子,即:
其中,为第1次迭代分解得到的分量;为第1次迭代分解得到的基线信号;在第1次迭代分解中,xt代表权利要求1所述步骤2中x(k);
3)再把作为新数据,重复上述步骤,可以分离出频率依次降低的固有旋转分量,直到基线信号变得单调为止;这样xt的整个分解过程可以写为:
;
其中代表第i次迭代分解得到的分量,代表第i次迭代分解得到的基线信号。
进一步地,所述步骤3中数据重排操作包括以下步骤:
随机打乱分量ci(k)的排列顺序。
进一步地,所述步骤3中数据替代操作包括以下步骤:
1) 对分量ci(k)执行离散傅里叶变换,获得分量ci(k)的相位;
2) 用一组位于(-π,π)区间内的伪独立同分布数来代替分量ci(k)的原始相位;
3) 对经过相位替代后的频域数据执行离散傅里叶逆变换得到数据ci IFFT(k),求取数据ci IFFT(k)的实部。
进一步地,所述步骤4中MFDFA方法包括以下步骤:
1)构造x(k)(k=1,2,…,N)的轮廓Y(i):
;
x(k)代表权利要求1所述步骤4中的ci(k)或ci shuffle(k)或ci FTran(k);
2)将信号轮廓Y(i)分成不重叠的NS段长度为s的数据,由于数据长度N通常不能整除s,所以会剩余一段数据不能利用;
为了充分利用数据的长度,再从数据的反方向以相同的长度分段,这样一共得到2NS段数据;
3)利用最小二乘法拟合每段数据的多项式趋势,然后计算每段数据的方差:
yv(i)为拟合的第v段数据的趋势,若拟合的多项式趋势为m阶,则记该去趋势过程为(MF-)DFAm;本例中,m=1;
4) 计算第q阶波动函数的平均值:
;
5)如果x(k)存在自相似特征,则第q阶波动函数的平均值Fq(s)和时间尺度s之间存在幂律关系:
当q=0时,步骤4)中的公式发散,这时H(0)通过下式所定义的对数平均过程来确定:
6)对步骤5)中的公式两边取对数可得ln[Fq(s)]=H(q)ln(s)+c(c为常数),由此可以获得直线的斜率H(q)。
进一步地,所述步骤7中的谱峭度方法包括以下步骤:
1)构造一个截止频率为fc=0.125+ε的低通滤波器h(n);ε>0,本例中fc=0.3;
2)基于h(n)构造通频带为[0, 0.25]的准低通滤波器h0(n)和通频带为[0.25, 0.5]的准高通滤波器h1(n),
;
3)信号ci k(n)经 h0(n)、 h1(n)滤波并降采样后分解成低频部分c2i k+1(n)和高频部分c2i+1 k+1(n),降采样的因子为2,再经多次迭代滤波后形成滤波器树,第k层有2k个频带,其中ci k(n)表示滤波器树中第k层上的第i个滤波器的输出信号,i=0,…, 2k-1,0≤k≤K-1,本例中K=8;c0 (n)代表权利要求1所述步骤7中xf1(k);
4)分解树中第k层上的第i个滤波器的中心频率fki和带宽Bk分别为
;
5)计算每一个滤波器结果ci k(n)( i=0,…, 2k-1) 的峭度;
6)将所有的谱峭度汇总,得到信号总的谱峭度。
进一步地,所述步骤9中的三次样条迭代平滑包络分析方法包括以下步骤:
1)计算信号z(k)的绝对值∣z(k)∣的局部极值;在第1次迭代中,z(k)代表权利要求1所述步骤9中xf2(k);
2)采用三次样条曲线拟合局部极值点得到包络线eov1(k);
3)对z(k)进行归一化处理得到;
4)第2次迭代:把z 1(k)重新作为新数据,重复执行上述步骤1)~3),得到;
5)第i次迭代:把z i-1(k) 重新作为新数据,重复执行上述步骤1)~3),得到;
6) 如果第n次迭代得到的z n (k)的幅值小于或等于1,则迭代过程停止,最后得到信号z(k)的包络为。
本发明采用以上技术方案,与现有技术相比,本发明具有以下优点:
1)利用内秉时间尺度分解对原始信号进行分解,然后利用数据的重排和替代操作排除其中的噪声和趋势分量,仅仅保留信号分量中的有用成分,从而避免了噪声和趋势分量对包络分析结果的影响,分析结果准确度和精确度高。
2)利用三次样条迭代平滑包络分析方法将信号包络与频率调制部分完全分离,能够避免频率调制部分对信号包络分析结果的影响,从而提高包络分析的精度。
3) 能够准确地检测出旋转机械的故障类型。
4) 由传统方法得到的包络谱存在端点效应,而由本发明得到的包络谱能够避免端点效应。
下面结合附图和实施例对本发明做进一步说明。
附图说明
附图1为本发明实施例中本发明方法的流程图;
附图2为本发明实施例中采用低通滤波器和高通滤波器对信号进行初步分解的示意图;
附图3为本发明实施例中采用树状滤波器结构快速计算谱峭度的示意图;
附图4为本发明实施例中具有内圈故障的滚动轴承振动信号;
附图5为本发明实施例中采用传统包络分析方法对内圈故障滚动轴承振动信号的分析结果;
附图6为本发明实施例中本发明对内圈故障滚动轴承振动信号的分析结果;
附图7为本发明实施例中具有外圈故障的滚动轴承振动信号;
附图8为本发明实施例中采用传统包络分析方法对外圈故障滚动轴承振动信号的分析结果;
附图9为本发明实施例中本发明对外圈故障滚动轴承振动信号的分析结果。
具体实施方式
实施例,如图1、图2、图3所示,一种基于内秉时间尺度分解和谱峭度的包络分析方法,包括以下步骤:
步骤1:利用加速度传感器以采样频率fs测取旋转机械的振动信号x(k), (k=1, 2,…,N),N为采样信号的长度;
步骤2:采用内秉时间尺度分解算法将信号x(k)分解成n个分量和一个趋势项之和,即,其中,ci(k)代表由内秉时间尺度分解算法得到的第i个分量,rn(k)代表由内秉时间尺度分解算法得到的趋势项;
步骤3:对ci(k)执行重排操作和替代操作,经重排操作得到的数据用ci shuffle(k)表示,替代操作后得到数据用ci FTran(k)表示;
步骤4:对ci(k)、ci shuffle(k)和ci FTran(k)分别执行多重分形去趋势波动分析(Multifractal Detrended Fluctuation Analysis, MFDFA),得到广义Hurst指数曲线,ci(k)的广义Hurst指数曲线用Hi(q)表示;ci shuffle(k)的广义Hurst指数曲线用Hi shuffle(q)表示;ci FTran(k)的广义Hurst指数曲线用Hi FTran(q)表示;
步骤5:如果Hi(q) 与Hi shuffle(q)或Hi(q) 与Hi FTran(q)之间的相对误差小于5%,或者Hi(q) 、Hi shuffle(q) 和Hi FTran(q)三者都不随q而变化,则抛弃对应的ci(k)分量;
步骤6:对剩余的ci(k)分量求和,将该和记为信号经重排和替代滤波后的结果xf1(k);
步骤7:对xf1(k)执行谱峭度分析,求出信号峭度最大处所对应的中心频率f0和带宽B;
步骤8: 根据中心频率f0和带宽B对xf1(k)进行带通滤波,得到xf2(k);
步骤9:对信号xf2(k)执行三次样条迭代平滑包络分析,得到信号包络eov(k);
步骤10:对得到的信号包络eov(k)执行离散傅里叶变换得到包络谱,根据包络谱特征频率判断机器的故障类型。
步骤2中内秉时间尺度分解算法包括以下步骤:
1) 对于任意信号xt,(t=1, 2, …,N),定义一个算子用于抽取低频基线信号,即:
其中是基线信号,是一个固有旋转分量, 假设是一个实值信号,代表xt的局部极值所对应的时刻,为方便起见定义;如果xt在某个区间上具有恒定值,考虑到邻近的信号存在着波动,我们仍然认为xt在这个区间上包含着极值,这时是该区间的右端点;为方便起见,定义,;假设和在上有定义,xk在上有定义,在区间上的连续极值点之间定义一个分段线性的基线信号抽取算子,即:
其中
这里参数是一个线性增益,,本例中;
2) 定义一个固有旋转分量抽取算子,即:
其中,为第1次迭代分解得到的分量;为第1次迭代分解得到的基线信号;在第1次迭代分解中,xt代表权利要求1所述步骤2中x(k);
3)再把作为新数据,重复上述步骤,可以分离出频率依次降低的固有旋转分量,直到基线信号变得单调为止;这样xt的整个分解过程可以写为:
;
其中代表第i次迭代分解得到的分量,代表第i次迭代分解得到的基线信号。
步骤3中数据重排操作包括以下步骤:
随机打乱分量ci(k)的排列顺序。
步骤3中数据替代操作包括以下步骤:
1) 对分量ci(k)执行离散傅里叶变换,获得分量ci(k)的相位;
2) 用一组位于(-π,π)区间内的伪独立同分布数来代替分量ci(k)的原始相位;
3) 对经过相位替代后的频域数据执行离散傅里叶逆变换得到数据ci IFFT(k),求取数据ci IFFT(k)的实部。
步骤4中MFDFA方法包括以下步骤:
1)构造x(k)(k=1,2,…,N)的轮廓Y(i):
;
x(k)代表权利要求1所述步骤4中的ci(k)或ci shuffle(k)或ci FTran(k);
2)将信号轮廓Y(i)分成不重叠的NS段长度为s的数据,由于数据长度N通常不能整除s,所以会剩余一段数据不能利用;
为了充分利用数据的长度,再从数据的反方向以相同的长度分段,这样一共得到2NS段数据;
3)利用最小二乘法拟合每段数据的多项式趋势,然后计算每段数据的方差:
yv(i)为拟合的第v段数据的趋势,若拟合的多项式趋势为m阶,则记该去趋势过程为(MF-)DFAm;本例中,m=1;
4) 计算第q阶波动函数的平均值:
;
5)如果x(k)存在自相似特征,则第q阶波动函数的平均值Fq(s)和时间尺度s之间存在幂律关系:
当q=0时,步骤4)中的公式发散,这时H(0)通过下式所定义的对数平均过程来确定:
6)对步骤5)中的公式两边取对数可得ln[Fq(s)]=H(q)ln(s)+c(c为常数),由此可以获得直线的斜率H(q)。
步骤7中的谱峭度方法包括以下步骤:
1)构造一个截止频率为fc=0.125+ε的低通滤波器h(n);ε>0,本例中fc=0.3;
2)基于h(n)构造通频带为[0, 0.25]的准低通滤波器h0(n)和通频带为[0.25, 0.5]的准高通滤波器h1(n),
;
3)信号ci k(n)经 h0(n)、 h1(n)滤波并降采样后分解成低频部分c2i k+1(n)和高频部分c2i+1 k+1(n),降采样的因子为2,再经多次迭代滤波后形成滤波器树,第k层有2k个频带,其中ci k(n)表示滤波器树中第k层上的第i个滤波器的输出信号,i=0,…, 2k-1,0≤k≤K-1,本例中K=8;c0 (n)代表权利要求1所述步骤7中xf1(k);
4)分解树中第k层上的第i个滤波器的中心频率fki和带宽Bk分别为
;
5)计算每一个滤波器结果ci k(n)( i=0,…, 2k-1) 的峭度;
6)将所有的谱峭度汇总,得到信号总的谱峭度。
步骤9中的三次样条迭代平滑包络分析方法包括以下步骤:
1)计算信号z(k)的绝对值∣z(k)∣的局部极值;在第1次迭代中,z(k)代表权利要求1所述步骤9中xf2(k);
2)采用三次样条曲线拟合局部极值点得到包络线eov1(k);
3)对z(k)进行归一化处理得到;
4)第2次迭代:把z 1(k)重新作为新数据,重复执行上述步骤1)~3),得到;
5)第i次迭代:把z i-1(k) 重新作为新数据,重复执行上述步骤1)~3),得到;
6) 如果第n次迭代得到的z n (k)的幅值小于或等于1,则迭代过程停止,最后得到信号z(k)的包络为。
试验1,利用具有内圈故障的滚动轴承振动数据对本发明所述算法的性能进行验证。
实验所用轴承为6205-2RS JEM SKF,利用电火花加工方法在轴承内圈上加工深度为0.2794mm、宽度为0.3556mm的凹槽来模拟轴承内圈故障,本实验负载约为0.7457kW,驱动电机转频约为29.5Hz,轴承内圈故障特征频率约为160Hz,采样频率为4.8KHz,信号采样时长为1s。
采集到的内圈故障信号如图4所示。
首先采用传统的包络分析方法对图4所示的信号进行分析,得到的分析结果如图5所示。从图5可以看出,轴承的故障特征完全被掩盖,因此传统的包络分析方法不能有效地提取轴承的故障特征;此外,从图5可以看出,包络谱的左端点存在着异常高值,这说明由传统方法得到的包络谱存在着端点效应。
采用本发明所提出的方法对图4所示的信号进行分析,得到的分析结果如图6所示。从图6可以看出,160Hz和320Hz所对应的谱线明显高于其它谱线,这两个频率分别对应轴承内圈故障特征频率的1倍频和2倍频,据此可以判断轴承具有内圈故障;从图6可以看出,由本发明得到的包络谱没有端点效应。
经多次实验表明,在负载和故障尺寸深度不变的情况下,本发明能够可靠识别的最小内圈故障尺寸宽度约为0.25 mm,而传统方法能够可靠识别的最小内圈故障尺寸宽度约为0.53mm,精度提高52.8%。
试验2,利用具有外圈故障的滚动轴承振动数据对本发明所述算法的性能进行验证。
实验所用轴承为6205-2RS JEM SKF,利用电火花加工方法在轴承外圈上加工深度为0.2794mm、宽度为0.5334mm的凹槽来模拟轴承外圈故障,本实验负载约为2.237 kW,驱动电机转频约为28.7Hz,轴承外圈故障特征频率约为103Hz,采样频率为4.8KHz,信号采样时长为1s。
采集到的外圈故障信号如图7所示。
首先采用传统的包络分析方法对图7所示的信号进行分析,得到的分析结果如图8所示。从图8可以看出,轴承的故障特征完全被掩盖,因此传统的包络分析方法不能有效地提取轴承的故障特征;此外,从图8可以看出,包络谱的左端点存在着异常高值,这说明由传统方法得到的包络谱存在着端点效应。
采用本发明所提出的方法对图7所示的信号进行分析,得到的分析结果如图9所示。从图9可以看出,103Hz和206Hz所对应的谱线明显高于其它谱线,这两个频率分别对应轴承外圈故障特征频率的1倍频和2倍频,据此可以判断轴承具有外圈故障;从图9可以看出,由本发明得到的包络谱没有端点效应。
经多次实验表明,在负载和故障尺寸深度不变的情况下,本发明能够可靠识别的最小外圈故障尺寸宽度约为0.35mm,而传统方法能够可靠识别的最小外圈故障尺寸宽度约为0.68mm,精度提高48.5%。
根据试验结果,分析后认为:
1) 传统的包络分析方法直接对原始信号进行包络分析,或者对仅经过简单处理后的原始信号进行包络分析,与传统的包络分析方法不同,本发明首先利用内秉时间尺度分解对原始信号进行分解,然后利用数据的重排和替代操作排除其中的噪声和趋势分量,仅仅保留信号分量中的有用成分,从而避免了噪声和趋势分量对包络分析结果的影响,提高了准确度和精确度。
2) 传统的包络分析方法以Hilbert变换为基础,而Hilbert变换要求被分析的信号必须是单分量的窄带信号,否则信号的频率调制部分将要污染信号的包络分析结果,但是目前待分析的信号都不严格满足单分量且窄带的条件,这样就会导致现有技术因精度不高而容易出现误判问题,与传统包络分析方法不同,本发明利用三次样条迭代平滑包络分析方法将信号包络与频率调制部分完全分离,能够避免频率调制部分对信号包络分析结果的影响,从而提高包络分析的精度。
3)能够准确地检测出旋转机械的故障类型。
4) 由传统方法得到的包络谱存在端点效应,而由本发明得到的包络谱能够避免端点效应。
5)各步骤作用:
第1)步:采集振动信号;
第2)步:将原始信号分解成不同分量和的形式,其中有些分量对应噪声和趋势项,有些分量对应有用信号;
第3)~5)步:对上述分解得到的信号执行重排操作和替代操作,剔除其中的噪声分量和趋势项,只保留有用信号;
第6)步:将剩余的有用信号求和,将该和作为信号经重排和替代滤波后的结果xf1(k);
第7)步:对滤波后的信号xf1(k)执行谱峭度分析,求出信号最大峭度处对应的中心频率f0和带宽B;
第8)步:根据中心频率f0和带宽B对xf1(k)进行带通滤波,得到信号xf2(k);
第9)步:计算信号xf2(k)的包络eov(k);
第10)步:对eov(k)执行离散傅里叶变换得到包络谱,根据包络谱判断轴承的故障类型。
本领域技术人员应该认识到,上述的具体实施方式只是示例性的,是为了使本领域技术人员能够更好的理解本发明内容,不应理解为是对本发明保护范围的限制,只要是根据本发明技术方案所作的改进,均落入本发明的保护范围。
Claims (7)
1.一种基于内秉时间尺度分解和谱峭度的包络分析方法,其特征在于,包括以下步骤:
步骤1:利用加速度传感器以采样频率fs测取旋转机械的振动信号x(k), (k=1, 2,…,N),N为采样信号的长度;
步骤2:采用内秉时间尺度分解算法将信号x(k)分解成n个分量和一个趋势项之和,即,其中,ci(k)代表由内秉时间尺度分解算法得到的第i个分量,rn(k)代表由内秉时间尺度分解算法得到的趋势项;
步骤3:对ci(k)执行重排操作和替代操作,经重排操作得到的数据用ci shuffle(k)表示,替代操作后得到数据用ci FTran(k)表示;
步骤4:对ci(k)、ci shuffle(k)和ci FTran(k)分别执行多重分形去趋势波动分析(Multifractal Detrended Fluctuation Analysis, MFDFA),得到广义Hurst指数曲线,ci(k)的广义Hurst指数曲线用Hi(q)表示;ci shuffle(k)的广义Hurst指数曲线用Hi shuffle(q)表示;ci FTran(k)的广义Hurst指数曲线用Hi FTran(q)表示;
步骤5:如果Hi(q) 与Hi shuffle(q)或Hi(q) 与Hi FTran(q)之间的相对误差小于5%,或者Hi(q) 、Hi shuffle(q) 和Hi FTran(q)三者都不随q而变化,则抛弃对应的ci(k)分量;
步骤6:对剩余的ci(k)分量求和,将该和记为信号经重排和替代滤波后的结果xf1(k);
步骤7:对xf1(k)执行谱峭度分析,求出信号峭度最大处所对应的中心频率f0和带宽B;
步骤8: 根据中心频率f0和带宽B对xf1(k)进行带通滤波,得到xf2(k);
步骤9:对信号xf2(k)执行三次样条迭代平滑包络分析,得到信号包络eov(k);
步骤10:对得到的信号包络eov(k)执行离散傅里叶变换得到包络谱,根据包络谱特征频率判断机器的故障类型。
2.根据权利要求1所述的一种基于内秉时间尺度分解和谱峭度的包络分析方法,其特征在于,所述步骤2中内秉时间尺度分解算法包括以下步骤:
1) 对于任意信号xt,(t=1, 2, …,N),定义一个算子用于抽取低频基线信号,即:
其中是基线信号,是一个固有旋转分量, 假设是一个实值信号,代表xt的局部极值所对应的时刻,为方便起见定义;如果xt在某个区间上具有恒定值,考虑到邻近的信号存在着波动,我们仍然认为xt在这个区间上包含着极值,这时是该区间的右端点;为方便起见,定义,; 假设和在上有定义,xk在上有定义,在区间上的连续极值点之间定义一个分段线性的基线信号抽取算子,即:
其中
这里参数是一个线性增益,,本例中;
2) 定义一个固有旋转分量抽取算子,即:
其中,为第1次迭代分解得到的分量;为第1次迭代分解得到的基线信号;在第1次迭代分解中,xt代表权利要求1所述步骤2中x(k);
3)再把作为新数据,重复上述步骤,可以分离出频率依次降低的固有旋转分量,直到基线信号变得单调为止;这样xt的整个分解过程可以写为:
;
其中代表第i次迭代分解得到的分量,代表第i次迭代分解得到的基线信号。
3.根据权利要求1所述的一种基于内秉时间尺度分解和谱峭度的包络分析方法,其特征在于,所述步骤3中数据重排操作包括以下步骤:
随机打乱分量ci(k)的排列顺序。
4.根据权利要求1所述的一种基于内秉时间尺度分解和谱峭度的包络分析方法,其特征在于:所述步骤3中数据替代操作包括以下步骤:
1) 对分量ci(k)执行离散傅里叶变换,获得分量ci(k)的相位;
2) 用一组位于(-π,π)区间内的伪独立同分布数来代替分量ci(k)的原始相位;
3) 对经过相位替代后的频域数据执行离散傅里叶逆变换得到数据ci IFFT(k),求取数据ci IFFT(k)的实部。
5.根据权利要求1所述的一种基于内秉时间尺度分解和谱峭度的包络分析方法,其特征在于:所述步骤4中MFDFA方法包括以下步骤:
1)构造x(k)(k=1,2,…,N)的轮廓Y(i):
;
x(k)代表权利要求1所述步骤4中的ci(k)或ci shuffle(k)或ci FTran(k);
2)将信号轮廓Y(i)分成不重叠的NS段长度为s的数据,从数据的反方向以相同的长度分段,得到2NS段数据;
3)利用最小二乘法拟合每段数据的多项式趋势,然后计算每段数据的方差:
yv(i)为拟合的第v段数据的趋势,若拟合的多项式趋势为m阶,则记该去趋势过程为(MF-)DFAm;
4) 计算第q阶波动函数的平均值:
;
5)如果x(k)存在自相似特征,则第q阶波动函数的平均值Fq(s)和时间尺度s之间存在幂律关系:
当q=0时,步骤4)中的公式发散,这时H(0)通过下式所定义的对数平均过程来确定:
6)对步骤5)中的公式两边取对数可得ln[Fq(s)]=H(q)ln(s)+c(c为常数),由此可以获得直线的斜率H(q)。
6.根据权利要求1所述的一种基于内秉时间尺度分解和谱峭度的包络分析方法,其特征在于:所述步骤7中的谱峭度方法包括以下步骤:
1)构造一个截止频率为fc=0.125+ε的低通滤波器h(n);
2)基于h(n)构造通频带为[0, 0.25]的准低通滤波器h0(n)和通频带为[0.25, 0.5]的准高通滤波器h1(n),
;
3)信号ci k(n)经 h0(n)、 h1(n)滤波并降采样后分解成低频部分c2i k+1(n)和高频部分c2i +1 k+1(n),降采样的因子为2,再经多次迭代滤波后形成滤波器树,第k层有2k个频带,其中ci k(n)表示滤波器树中第k层上的第i个滤波器的输出信号,i=0,…, 2k-1,0≤k≤K-1;
4)分解树中第k层上的第i个滤波器的中心频率fki和带宽Bk分别为
;
5)计算每一个滤波器结果ci k(n)( i=0,…, 2k-1) 的峭度;
6)将所有的谱峭度汇总,得到信号总的谱峭度。
7.根据权利要求1所述的一种基于内秉时间尺度分解和谱峭度的包络分析方法,其特征在于,所述步骤9中的三次样条迭代平滑包络分析方法包括以下步骤:
1)计算信号z(k)的绝对值∣z(k)∣的局部极值;在第1次迭代中,z(k)代表权利要求1所述步骤9中xf2(k);
2)采用三次样条曲线拟合局部极值点得到包络线eov1(k);
3)对z(k)进行归一化处理得到;
4)第2次迭代:把z 1(k)重新作为新数据,重复执行上述步骤1)~3),得到;
5)第i次迭代:把z i-1(k) 重新作为新数据,重复执行上述步骤1)~3),得到;
6) 如果第n次迭代得到的z n (k)的幅值小于或等于1,则迭代过程停止,最后得到信号z(k)的包络为。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610492062.8A CN105954030B (zh) | 2016-06-29 | 2016-06-29 | 一种基于内秉时间尺度分解和谱峭度的包络分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610492062.8A CN105954030B (zh) | 2016-06-29 | 2016-06-29 | 一种基于内秉时间尺度分解和谱峭度的包络分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105954030A true CN105954030A (zh) | 2016-09-21 |
CN105954030B CN105954030B (zh) | 2018-03-23 |
Family
ID=56902590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610492062.8A Expired - Fee Related CN105954030B (zh) | 2016-06-29 | 2016-06-29 | 一种基于内秉时间尺度分解和谱峭度的包络分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105954030B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291293A (zh) * | 2016-10-27 | 2017-01-04 | 西南石油大学 | 一种基于谱峭度与s变换的局放信号自适应去噪方法 |
CN108152363A (zh) * | 2017-12-21 | 2018-06-12 | 北京工业大学 | 一种抑端固有时间尺度分解的管道缺陷识别方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN202793793U (zh) * | 2012-08-30 | 2013-03-13 | 桂林电子科技大学 | 大型风力机组轴承故障诊断系统 |
CN103575523A (zh) * | 2013-11-14 | 2014-02-12 | 哈尔滨工程大学 | 基于FastICA-谱峭度-包络谱分析的旋转机械故障诊断方法 |
CN104198186A (zh) * | 2014-08-29 | 2014-12-10 | 南京理工大学 | 基于小波包与谱峭度结合的齿轮故障诊断方法及装置 |
CN104677632A (zh) * | 2015-01-21 | 2015-06-03 | 大连理工大学 | 利用粒子滤波与谱峭度的滚动轴承故障诊断方法 |
CN104792528A (zh) * | 2014-01-22 | 2015-07-22 | 中国人民解放军海军工程大学 | 一种自适应最优包络解调方法 |
KR101607047B1 (ko) * | 2015-01-12 | 2016-03-28 | 울산대학교 산학협력단 | 결함 검출을 위한 신호 분석 방법 및 장치 |
-
2016
- 2016-06-29 CN CN201610492062.8A patent/CN105954030B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN202793793U (zh) * | 2012-08-30 | 2013-03-13 | 桂林电子科技大学 | 大型风力机组轴承故障诊断系统 |
CN103575523A (zh) * | 2013-11-14 | 2014-02-12 | 哈尔滨工程大学 | 基于FastICA-谱峭度-包络谱分析的旋转机械故障诊断方法 |
CN104792528A (zh) * | 2014-01-22 | 2015-07-22 | 中国人民解放军海军工程大学 | 一种自适应最优包络解调方法 |
CN104198186A (zh) * | 2014-08-29 | 2014-12-10 | 南京理工大学 | 基于小波包与谱峭度结合的齿轮故障诊断方法及装置 |
KR101607047B1 (ko) * | 2015-01-12 | 2016-03-28 | 울산대학교 산학협력단 | 결함 검출을 위한 신호 분석 방법 및 장치 |
CN104677632A (zh) * | 2015-01-21 | 2015-06-03 | 大连理工大学 | 利用粒子滤波与谱峭度的滚动轴承故障诊断方法 |
Non-Patent Citations (1)
Title |
---|
林近山等: "多重分形去趋势波动分析在滚动轴承损伤程度识别中的应用", 《中国机械工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291293A (zh) * | 2016-10-27 | 2017-01-04 | 西南石油大学 | 一种基于谱峭度与s变换的局放信号自适应去噪方法 |
CN106291293B (zh) * | 2016-10-27 | 2018-09-18 | 西南石油大学 | 一种基于谱峭度与s变换的局放信号自适应去噪方法 |
CN108152363A (zh) * | 2017-12-21 | 2018-06-12 | 北京工业大学 | 一种抑端固有时间尺度分解的管道缺陷识别方法 |
CN108152363B (zh) * | 2017-12-21 | 2021-06-25 | 北京工业大学 | 一种抑端固有时间尺度分解的管道缺陷识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105954030B (zh) | 2018-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106168538B (zh) | 一种滚动轴承的itd、谱峭度和平滑迭代包络分析方法 | |
CN106198015B (zh) | 一种滚动轴承的vmd、谱峭度和平滑迭代包络分析方法 | |
CN106153339B (zh) | 一种基于变分模式分解滤波的包络分析方法 | |
CN106096313B (zh) | 一种基于奇异谱分解和谱峭度的包络分析方法 | |
CN106096200B (zh) | 一种基于小波分解和谱峭度的包络分析方法 | |
CN106096198B (zh) | 一种基于变分模式分解和谱峭度的包络分析方法 | |
CN106096199B (zh) | 一种滚动轴承的wt、谱峭度和平滑迭代包络分析方法 | |
CN105954030B (zh) | 一种基于内秉时间尺度分解和谱峭度的包络分析方法 | |
CN106053069B (zh) | 一种滚动轴承的ssd、谱峭度和平滑迭代包络分析方法 | |
CN106198009B (zh) | 一种滚动轴承的emd、谱峭度和平滑迭代包络分析方法 | |
CN105954031B (zh) | 一种基于奇异谱分解滤波的包络分析方法 | |
CN106053059B (zh) | 一种基于内秉时间尺度分解滤波的包络分析方法 | |
CN106153333B (zh) | 一种基于小波分解滤波的包络分析方法 | |
CN106053060B (zh) | 一种基于非线性模式分解滤波的包络分析方法 | |
CN106198010B (zh) | 一种基于局部均值分解滤波的包络分析方法 | |
CN106198012B (zh) | 一种基于局部均值分解和谱峭度的包络分析方法 | |
CN106053061B (zh) | 一种基于非线性模式分解和谱峭度的包络分析方法 | |
CN112683393A (zh) | 一种lmd设备故障诊断方法及系统 | |
CN106198017B (zh) | 一种滚动轴承的lmd、谱峭度和平滑迭代包络分析方法 | |
CN106096201B (zh) | 一种旋转机械的eemd和三次样条平滑包络分析方法 | |
CN106198016B (zh) | 一种滚动轴承的nmd、谱峭度和平滑迭代包络分析方法 | |
CN106124200B (zh) | 一种旋转机械的elmd和平滑迭代包络分析方法 | |
CN112697480A (zh) | 一种nmd设备故障诊断方法及系统 | |
CN112697478A (zh) | 一种wd设备故障诊断方法及系统 | |
CN112697469A (zh) | 一种lmd和gzc机器状态监测方法及装置 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180323 Termination date: 20210629 |
|
CF01 | Termination of patent right due to non-payment of annual fee |