摩擦振动信号的EEMD和多重分形去趋势波动分析

摩擦振动信号的EEMD和多重分形去趋势波动分析李精明;魏海军;魏立队;孙迪;杨智远;梅立强【摘要】To investigate the variation rules of frictional vibration inrunning⁃in wear processes and to identify wear states through frictional vibration, we conducted experiments on a testing machine on the friction and wear of a pis⁃ton ring against a cylinder in a marine diesel engine. We decomposed the frictional vibration signals and acquired several intrinsic mode functions ( IMFs ) without mode mixing by using ensemble empirical mode decomposition ( EEMD) . Then, we analyzed the resynthesized characteristic signals of frictional vibration by utilizing the multi⁃fractal detrended fluctuation analysis ( MFDFA) algorithm to derive the MFDFA spectrum and its parameters. The results show that we can use EEMD to extract the weak characteristic signal of frictional vibration and then usethe MFDFA spectrum and its parameters to characterize the frictional vibration signals.%为了研究摩擦副磨合磨损过程中摩擦振动变化规律，实现通过摩擦振动识别摩擦副的磨合磨损状态，在摩擦磨损试验机上进行了船用柴油机缸套—活塞环摩擦副摩擦磨损试验。

应用总体经验模式分解对摩擦振动信号进行分解，获得若干个无模式混叠的本征模式分量。

利用多重分形去趋势波动分析（ Multifractal detrended fluctua⁃tion analysis，MFDFA）对重构获得的摩擦振动特征信号进行分析，得到摩擦振动信号的MFDFA谱图，并根据谱图求取摩擦振动信号的多重分形谱参数。

研究结果表明，总体经验模式分解能够实现微弱摩擦振动特征信号的提取， MFDFA谱图及其参数可以表征摩擦振动信号的特征。

【期刊名称】《哈尔滨工程大学学报》【年(卷),期】2016(037)009【总页数】6页(P1204-1208,1214)【关键词】总体经验模式分解;多重分形去趋势波动分析;谱参数;摩擦振动;Hurst指数;特征提取【作者】李精明;魏海军;魏立队;孙迪;杨智远;梅立强【作者单位】上海海事大学商船学院，上海201306; 大连海事大学轮机工程学院，辽宁大连116026;上海海事大学商船学院，上海201306;上海海事大学商船学院，上海201306;大连海事大学轮机工程学院，辽宁大连116026;上海海事大学商船学院，上海201306;上海海事大学商船学院，上海201306【正文语种】中文【中图分类】TH117.1摩擦副磨合过程中产生的摩擦振动现象，蕴含着反映磨合状态的信息[1]。

摩擦振动信号是微弱信号，往往埋没于背景噪声之中，未经处理的摩擦振动信号不能真实地反映摩擦振动特征[2]。

因此，如何对获得的摩擦振动信号进行消噪，以及根据提取的有用信号找出摩擦振动信号的特征，是摩擦振动研究的关键问题。

经验模式分解(empirical mode decomposition,EMD)是美籍华人Huang等[3] 1998年提出的一种自适应数据驱动的信号处理方法，能清晰地分辨出非平稳非线性复杂数据的内蕴模式，但存在模式混叠问题[4]。

Wu等[5]基于白噪声信号的统计特征，提出了总体经验模式分解(ensemble empirical mode decomposition, EEMD)，解决了EMD模式混叠问题，该方法成功地应用于信号处理[6]、故障诊断[7-8]等的研究。

多重分形去趋势波动分析(multifractal detrended fluctuation analysis, MFDFA)是Kantelhardt[9]于2002年提出的基于去趋势波动分析(detrended fluctuation analysis, DFA)的非稳态时间序列分析方法，可以有效地分析非线性非平稳信号的多重分形特征。

该方法已被应用到地质学[10]、脑电波信号处理[11]、交通流分析[12]、金融市场[13]等领域，取得很好的效果。

摩擦振动是摩擦副磨合磨损过程中产生的现象，其变化规律必然能反映系统摩擦学的状态和特征[14]。

李国宾等[2]应用谐波小波提取摩擦振动的特征，黄朝明等应用时频谱图[15]研究了特征参数与摩擦振动的关系，孙迪等应用奇异值分解[16]和混沌吸引子[17]探讨了摩擦振动在磨合磨损过程中的变化规律，通过摩擦振动识别磨合磨损状态和特征。

本文针对船用柴油机缸套—活塞环摩擦副摩擦磨损试验过程中获得的摩擦振动信号进行EEMD分解，获得若干个无模式混叠的IMF分量，根据摩擦振动信号的特点，选取包含摩擦振动信号特征的IMF分量重构摩擦振动特征信号。

应用多重分形去趋势波动分析算法对摩擦振动特征信号进行分析，得到摩擦振动信号的Hurst指数、多重分形标度指数以及多重分形谱。

并求取摩擦振动信号多重分形谱的宽度、极大值、维差以及驻点值，实现特征参数对摩擦副摩擦磨损状态的表征，为基于摩擦振动信号的机械摩擦副摩擦磨损行为的研究提供了新的途径。

1.1 试验材料与方法试验的设备采用CFT-I型摩擦磨损试验机(见图1)，采用销-盘摩擦副作为配副。

试样均用线切割机从船用柴油机缸套和活塞环截取，盘试样作为下式样，取自船用柴油机的缸套，φ30 mm，原始表面粗糙度Ra=1.72 μm，材质为合金铸铁，硬度HV300～400；销试样作为上试样，取自和缸套配对的活塞环，矩形截面尺寸3 mm×4 mm,原始表面粗糙度Ra=0.67 μm，材质为合金铸铁，硬度HV600～720。

上试样用专用夹具固定不动；下试样由专用夹具固定在台架上，台架经电机驱动的偏心机构和连杆机构将回转运动转变为往复运动，运动行程5 mm，电机转速600 r/min，选用船用Mobilgard-412润滑油进行滴油润滑。

载荷通过加载弹簧经销试样施加到盘试样上，施加的正压力为30 N，名义接触载荷2.5 MPa，试验时间600 min。

1.2 摩擦振动信号的采集应用NI公司生产的PXIe-1071信号采集系统采集摩擦振动信号，采样频率25 600 Hz，采样点数4 096，采样间隔2 min，每次采样时间0.16 s。

测量摩擦振动信号的传感器采用PCB PIEZOTRONICS公司生产的356A16型ICP三轴加速度传感器，灵敏度100 mV/g，量程50 g。

加速度传感器水平安装在盘试样下方，随盘试样一起做往复运动。

图2为缸套活塞环摩擦副摩擦磨损试验在初期、中期、末期采集的摩擦振动信号时域波形。

从图2可以看出，试验获得的摩擦振动信号是非线性非平稳信号，微弱的摩擦振动信号埋没于背景噪声中，信号波动复杂，时域波形体现不出变化，如果直接采用此信号来分析，则无法提取正确的摩擦振动特征。

EEMD方法的本质是将待分析信号与高斯白噪声叠加，再进行多次EMD分解，利用具有频率均匀分布统计特性的高斯白噪声使待分析信号在不同尺度上具有连续性，从而降低各IMF分量的模式混叠程度。

根据零均值高斯白噪声的特性，通过若干组IMF总体平均使加入的白噪声相互抵消，还原被分析信号。

EEMD算法归纳如下：1)初始化EMD总体平均次数M和加入的白噪声幅值系数k，令m=1。

2)执行第m次EMD分解：①对待分析信号x(t)加入一个给定幅度的高斯白噪声序列nm(t)，得到第m次加噪后的信号xm(t)：②用EMD分解xm(t)，得到一组IMFcj,m(j=1，2，…，I)，其中，cj,m为第m次分解得到的第j个IMF；③若m<M，则返回步骤2)，m=m+1；3)对M次分解得到的各IMF计算均值：4)输出作为EEMD分解得到的第j个IMF。

白噪声幅值系数k影响着信号的分解精度，通常k的取值范围是0.1～0.4。

总体平均次数M影响着信号的消噪能力和计算时间，当M增大到一定值后，EEMD对信号的消噪效果增强不明显，但计算时间显著增加。

应用总体经验模式分解对摩擦振动信号进行分解，白噪声幅值系数k取0.1，总体平均次数M取100。

图3是缸套活塞环摩擦副摩擦磨损试验初期的摩擦振动信号EEMD分解结果,分解得到8个IMF分量C1～C8和一个残差r8，限于篇幅，摩擦磨损试验中期、后期的摩擦振动EEMD分解结果图从略。

从图3可以看出，分解得到的IMF分量可使试验获得的原始摩擦振动信号在不同的分辨率下显现出来，摩擦振动信号具有频率高振幅小等特点[2]，因此选择前两个分量C1、C2重新合成摩擦振动特征信号。

图4为选择的IMF分量重构得到的摩擦振动特征信号，可以看出IMF分量重构信号的振幅随着磨合的进行呈现明显的趋势变化，摩擦振动的冲击信息显著出现。

3.1 多重分形去趋势波动分析理论及算法多重分形去趋势波动分析方法是Kantelhardt在去趋势波动分析的基础上提出的非稳态时间序列分析方法，MFDFA 算法步骤如下：1)计算时间序列{xi}(i=1，2，…，N)的离差序列Y(i):式中是均值。

2)将序列Y(i)划分成Ns=int(N/s)个区间，每个连续不重叠的区间均含有s个数据，当N不能整除s时，Y(i)会有数据剩余，为数据不丢失，再从序列反向开始重复这一分割过程，得到2Ns个等长小区间，包含序列的所有数据。

3)计算均方误差F2(s,v)，以区间(v=1，2，…，2Ns)为例，进行k阶多项式拟合：对于区间(v=1，2，…，Ns)：对于区间(v=Ns+1，…，2Ns)：4)计算q阶波动函数F(q,s):式中，q为非0实数。

5)阶数q依次取某数值，尺度s取不同值，重复步骤2)～4)，计算F(q,s)对q的双对数值，确定波动函数与尺度之间存在的幂率关系：式中，h(q)为Hurst指数。