基于错位叠加法的同步液压马达冲击故障成分的降噪提取.pdf

第48 卷第2 期2024年2 月doi:10.11832/j.issn.1000-4858.2024.02.018基于错位叠加法的同步液压马达冲击故障成分的降噪提取液压与气动Chinese Hydraulics&PneumaticsVol.48 No.2February.2024鲁晓，郭鹏伟，宁大勇，侯交义（大连海事大学船舶与海洋工程学院，辽宁大连1 1 6 0 0 0）摘要：同步液压马达经常工作于高负载条件下，对其进行在线故障监测可有效避免严重安全事故和经济损失。为提高故障监测的准确性,提出了一种基于改进的错位叠加法的同步液压马达故障成分自适应降噪以及提取方法，并以同步液压马达的齿轮磨损状态的声信号为对象进行了研究。实验表明，该方法对同步液压马达齿轮磨损状态的冲击故障成分具有良好的降噪和提取效果。关键词：同步液压马达；错位叠加法；冲击故障；降噪提取中图分类号：TH137文献标志码：B文章编号：1 0 0 0-48 58(2 0 2 4)0 2-0 1 3 7-0 7Denoising Extraction of Impact Fault Components of SynchronousHydraulic Motor Based on Displacement Superposition MethodLU Xiao,GUO Peng-wei,NING Da-yong,HOU Jiao-yi(School of Naval Architecture and Ocean Engineering,Dalian Maritime University,Dalian,Liaoning 116000)Abstract:Synchronous hydraulic motor often works under high load condition and on-line fault monitoring caneffectively avoid serious safety accidents and economic losses.In order to improve the accuracy of fault monitoring,an adaptive noise reduction and extraction method of fault components of synchronous hydraulic motor based onimproved dislocation superposition method is proposed in this paper,and the acoustic signals of gear wear state ofsynchronous hydraulic motor are studied.The experimental results show that this method has good noise reductionand extraction effect on the impact fault sound signal of synchronous hydraulic motor gear wear state.Key words:synchronous hydraulic motor,dislocation superposition method,impact fault,the noise reductionextraction引言同步液压马达广泛用于大负载多缸液压同步系统中，目前尚缺乏对其进行在线诊断的方法。对同步液压马达的故障成分进行有效的降噪提取将会为诊断其故障打下良好的基础。本研究从旋转机械故障诊断方向人手,研究了液压马达的故障诊断。目前,对于旋转机械故障诊断的研究趋势是通过处理信号进行故障诊断，而声信号作为一种以空气为介质的振动信号，具有非接触式采集、安装简便、信号易于获取等优点，可用于提取故障特征并进行诊断 1-3 国内外许多学者都对于机械故障诊断领域中声信号的特征提取进行了研究。王广逸等 4 利用非接触式检测技术，通过分析汽车发动机运行时的异常响动来检测故障;杜正昱等 5提出一种加权平均时变滤波经验模态分解及频谱幅度调制的滚动轴承声信号特征提取方法，实验表明该方法相较于常用的快速谱峭度方法能够更加有效地提取声信号故障特征频率；汪欣等 6 设计了一种基于声信号和一维卷积神经网络结构的电机故障诊断方法，实现了对电机故障声信号的收稿日期：2 0 2 2-0 2-1 6作者简介：鲁晓（1 9 9 8 一）,男,内蒙古呼和浩特人,硕士,主要从事旋转机械故障诊断领域的研究工作。修回日期：2 0 2 3-0 6-2 6138有效分类;YAO Jiachi 等 7 提出了一种基于声信号的行星齿轮箱的故障诊断方法，利用傅立叶分解方法分解声信号，再采用能量综合特征参数和时间和包络谱峰度来克服噪声和微弱的声信号；最后采用随机森林分类算法建立故障诊断方法，实验显示该方法的故障诊断准确率极高；STEVEKM等 8 本研究提出了一种利用发动机缸盖声信号进行发动机燃烧故障分类的自动模型，故障诊断的准确度和识别速度都更为良好。在常用的降噪方法中,小波变换作为一种广泛应用的信号处理工具，在时域和频域上拥有强大的多分辨率分析能力。于红梅 9 将完整齿轮、裂纹齿轮和缺齿齿轮3 种工况下的齿轮声学信号进行小波分析并构建了深度自编码网络和模糊推理相结合的诊断系统，实现了对齿轮故障的诊断和识别。武城同 1 0 模拟并采集了风电机齿轮箱的断齿2 0%、断齿40%齿面点蚀的3 种故障并采用小波包变换与Hilbert变换相结合的方法对声信号进行处理，可用于齿轮故障的早期检测。本研究针对同步液压马达的冲击故障成分的降噪提取方法，利用错位叠加法（DSM)对同步液压马达的齿轮磨损状态的冲击故障声信号进行叠加降噪并提取出来,其研究方法对于之后判别同步液压马达的故障种类打下了基础。1 DSM的基本原理DSM的数学公式如下：&(n)=K+1AK1Zx(n+hL)K=0式中，（n)为原始信号；（n）为经过DSM处理后的信号，也称为叠加信号；K为叠加次数；L为叠加长度（一般为待处理信号的周期）。图1 9 1 是DSM对信号进行处理的图示说明,其中A是要获取的目标信号;B是干扰信号；C是待处理的信号，即A和B的混合信号；对信号C使用式（1）进行叠加，设叠加次数K=10，叠加长度L为信号A的周期得到的叠加信号为信号D。显然，相比于原始的待处理信号C,在经过DSM算法提取后生成的信号D中，需要的目标信号A的成分比例有了明显的增加，而干扰信号B的成分比例被有效地削减了。用皮尔森相关系数检测目标信号在经过DSM处理后与原始目标信号的相似性，以此评估DSM的处理效果。根据Cauchy-Schwarz不等式，皮尔液压与气动森系数的范围在0 1 之间,数值越接近0,意味着两者的线性相关性越低,反之,数值越接近1,两者线性相关性越强。2-2020-20202020020-202同步液压马达故障成分的降噪及提取2.1同步液压马达的故障声信号在齿轮式液压马达中，存在着各种各样的误差和振动，这就导致了其声信号并非是标准的周期信号，而是属于一种准周期信号，准周期信号的示意图如图2所示。其中1 是平均周期,和l,分别表示每个周期长度的偏差值。(1)1+A2值采样点数图2 准周期信号的示意图Fig.2 Schematic diagram of quasi-periodic signal由于故障信号是准周期信号，,这就会导致在DSM的过程中每一段信号的叠加步长的不同。这可能会使冲击故障分量的相位产生偏差,进而会导致故障分量信息在叠加过程中被破坏。根据冲击故障特征，人为地将故障信号划分为“激变区”和“平稳区”，如图3 所示。平稳区包含的故障能量和信息较少，而激变区包含较多能量和信息。因为激变区具有短暂性、周期性、能量大等特点,而且不易被背景噪声淹没,所以选择将第48 卷第2 期200400200400400200400图1 DSM方法示意图Fig.1 Schematic diagram of DSM1+Al600600600600采样点数800800800800100010001000100012001200120012002024年第2 期激变区图3 故故障信号的激变区与平稳区示意图Fig.3Schematic diagram of shock region andstable region of fault signal激变区的声信号作为冲击故障成分提取出来 1 2-1 4。2.2基于错位叠加法的降噪和提取方法选择与冲击故障成分波形相似的母小波进行连续小波变换得到小波系数，并对小波系数进行硬阈值化处理,最后利用阈值小波系数进行加窗函数提取冲击故障成分的特征点位置，叠加降噪并提取冲击故障成分。具体步骤如下：（1）导人原始的故障信号S（n）,选择与故障信号激变区相似的母小波。由于本研究把冲击故障信号的“激变区”作为冲击故障成分，因此应该选择与冲击故障信号“激变区”的波形相似程度较高的小波函数作为小波基函数。通过对小波函数的筛选,本研究选择了dbN小波系中的db5小波作为冲击故障信号连续小波变换的小波基函数,如图4所示。对信号 S(n)进行连续小波变换,得到小波系数WS(a,b)。1.00.50-0.5-1.00图4db5波形图Fig.4db5 waveform（2）设置一个硬阈值函数来消除含有干扰成分的较小的小波系数，进而保留含有冲击故障成分的较大的小波系数。硬阈值函数的数学模型,如式(2)所示：=Ws(a,b),1 Wus(a,b)/nn(W,S(a,b),入)=lo,式中，m（WS(a,b),)为硬阈值函数处理后的小波系数,记为m;WS(,b)为小波系数;入为设置的硬液压与气动平稳区阅值大小。（3）为了提取每个周期的故障成分的特征点位置,选取两个斜对的时频点（ti,fi），（t 2,f)形成矩形窗R,从阈值化后的小波时频图中框选出一个时频区采样点块,该时频区块记作m（o,bo）,且选取时需要满足tit2,fi0且 nHCi+2=0,则取Ci+1（t i+1，f i+1)为冲击故障成分的特征点。条件二,若nHi=0,InHGi+1|0,nHci+2|0,ImHci+n|0,HCi+1=0时,取 max（|mH Ci+1 l,1H i+2 l,，|n H c i+n l 1 对应的C点为冲击故障成分的特征点。最后把冲击故障成分的特征点全部赋值于Q（t i，Ji）,故冲击故障成分每个周期特征点的点集Q1,Q2，,Q,其中j=1,2,3,。利用每个周期的冲击故障信号的特征点位置估计故障信号的周期,则故障信号246XI W,S(a,b)/(2)139nt8的周期大约为 T,=t2-tu,T,=t3-t22，,T-1=t,-t-1)C-1)。可以通过同步液压马达实时流速与同步液压马达排量V计算理论故障周期lo来判断特征点选取是否有效，公式如下：lo=fsV1000q（4）由于连续小波变换时是根据激变区的波形选取的母小波,且冲击故障成分的特征点是依据冲击故障信号的激变区确定的,即,所有的特征点均处于激变区的某个位置。根据已获取得特征点Q（t i，Ji),可以按照式(5)得到其对应的采样点数：ni=ti:fs式中，n;为故障特征点对应的采样点数;t,为冲击故障60(4)(5)140成分的特征点的时间。根据采样点就可以确定每个周期的冲击故障成分的范围：从第n；个采样点向前、向后分别扩展(-1)个采样点,该信号段记作W（n），则W(n,)段信号包括故障成分。n,=n j-+1,n;+-1,为激变区的长度。（5）由于背景噪声的干扰，每个周期的故障成分特征点在激变区的位置都存在差异，根据故障特征点的位置和故障成分的长度（激变区的长度）进一步计算得到,偏移量的范围为0 T-1,如果不进行激变区位置的调整，直接进行叠加计算，偏移量T会影响错位叠加效果，导致叠加结果不准确，甚至会破坏冲击故障成分。为了解决激变区偏移的问题,取任意一段含有故障成分的W（n j）为基准信号,记作Ws（n s）。根据偏移量的范围为0-1,把W（n;)再向前、向后分别扩展-1个采样点，得到扩展信号段w（n i），其中n,=n j-2+2,ni+2g-2。从w(n,)中第1个采样点连续截取长度为2 1 个采样点记作信号段（n）,将,设置为W（n）在W（n,）上的偏移量,将w(n,)偏移后的信号记录为W(n;+T,)。将Tj的子集表示为L,将L设置为1 0:2-2。计算出W（n;+T)中偏移量,的最佳偏移量,进一步就可以得到最佳的起始叠加点。其中叠长度为2 1 个采样点。遍历L中T所有的值,分别计算基准信号段Ws（n s）与Wj（n j+,）的皮尔森相关系数p（Ws（n s）,W(n+,)）,选取p（Ws（n s）,W（n;+,)最大值时所对应的T,作为最佳偏移量,如式(6)所示：(ri,r,.,i)=.armax.p(W(ns),W(n;+T,)式中，argmax，为最大值自变量点集。选取连续的2 K个W（n+T）信号，根据序号j的奇偶性将信号段分成奇偶两组，并对每组按照序号由低到高进行不同次数的叠加运算。方程式如下：Dx=+w(n+)+W(ns+t)+.+Wk-(n2k-1+T2k-I)K+IW(na+Ti)+W(n+T.)+Dw2K.+Wk(n2k+T2k)p(K)=p(DwIK,Dw2k)液压与气动其中,Dwlk为奇数组叠加后的信号;Dw2k为偶数组叠加后的信号;K为叠加次数；W2k-1（n 2 K-1+T 2 k-1）为含有故障成分的奇数组信号段；Wk（n 2 k+T k）为含有故障成分的偶数组信号段;p(K)为奇数组叠加后的信号与偶数组叠加后的信号的相关性。设置阈值M,当p(K)M时,停止迭代,并输出叠加结果。2.3参数说明1）硬阈值入的设置设置硬阈值入是为了去除干扰噪声的小波系数，当入的值选取过小会导致背景噪声产生小波系数可能去除不完全，进而使故障周期性不明显；当入的值过大不仅会去除较小的小波系数，而且会去除部分含有故障成分的较大的小波系数，导致故障信号的周期呈现跳跃性,但是最终还是能够找到故障成分的位置，不影响叠加效果。因此硬阈值入的值尽量取较大。2）矩形窗R的设置矩形窗R是由2 个斜对的时频点（ti,fi），（t,f)构成的,则矩形窗 R的长l=t-ti,宽w=f-fi。考虑到冲击故障信号的频率分布和冲击故障成分的长度Q,激变区的时间为：t。=fs矩形窗从t=0时刻开始,因此t=0。由于冲击故障声信号的特征频率一般分布在中低频，故f可以设置为0 50 0 0 Hz之间，fz5000Hz。为了防止矩形窗过长或过短引发的矩形窗框选多个故障信号周期和计算量增加的问题,其长度的一般取值范围为2 t。14t。TjeL,j=1,2,3,.3）激变区的长度(6)激变区的长度一般是人为选取故障成分的高幅值区域的长度。一般的经验公式如式（9），f。为声信号的采样频率：(9)441004）阈值M的设置随着叠加次数K不断越大，相关系数p不会一直增加，当p一般增大到0.9 附近就会趋于平稳，不再有明显的增大的趋势 1 5。为了使故障成分达到更好的降噪效果,在这里取M=0.90。3试验台的搭建试验台由同步液压马达液压系统、声音传感器、数(7)据采集卡和计算机组成，根据实际情况与试验条件，制第48 卷第2 期(8)Js3002024年第2 期图5同步液压马达液压系统实物图Fig.5 Physical diagram of hydraulic system ofsynchronous hydraulic motor造了齿轮磨损故障状态。同步液压马达液压系统实物,如图5所示。其中声音传感器放置于同步液压系统的同步液压马达旁，用于采集从同步液压马达中产生的故障声信号。声信号经过声音传感器传人NI-USB-6341数据采集卡，将模拟信号转换成计算机可识别的数字信号并通过计算机接口将数据实时导人计算机中。同步液压马达的型号为FMA-4R2.1S,参数如表1 所示。表1 同步液压马达参数Tab.1Parameters of synchronous hydraulic motor参数数值排量/mL r-12.17最小输人流量（7 50 r/min)/Lmin=l6.51标准输人流量（1 50 0 r/min)/Lmin=l13.23最大输人流量（3 0 0 0 r/min)/Lmin=126.04最大出口压力/MPa214实验结果图6 是同步液压马达发生齿轮磨损故障时的声信号。此时同步液压马达出口压力6 MPa，流速q为7.5L/min。已知采样频率为2 0 0 0 0 Hz,根据经验公式(9）,可以得到冲击故障成分的长度大约为1 50 个采样点。偏移量T，的子集L设置为1 0:2 9 8）。0-20图6 齿轮磨损故障时的实际声信号Fig.6 Actual acoustic signals of gear wear failure液压与气动矩形窗从t=0时刻开始选取,即t=0。由于冲击故障声信号的特征频率一般分布在中低频,所以f设置区间为0 50 0 Hz,f5000Hz。为了防止矩形窗过长或过短引发的矩形窗框选多个故障信号周期和计算量增加的问题,其长度的一般取值范围为2 t。l4t。,这里设置矩形窗R的2 个斜对的时频点分别为（0,0），（0.0 7,50 0 0），设置阈值M=0.9使故障成分达到很好的降噪效果。图7 为齿轮磨损故障信号的小波时频图。根据小波时频图的最大小波系数，设置硬阈值函数的阅值为8。10000r9000F8000F7000F500040003000200010000图7齿轮磨损故障信号小波变换后的小波时频图Fig.7Wavelet time-frequency diagram of gear wearfault signal after wavelet transform通过小波系数硬阈值处理与冲击故障特征点提取算法,可以提取断齿每个周期的冲击故障成分的特征点,并计算2 个特征点之间的差值用来估算故障信号的周期,如表2 所示。表2 特征点的位置Tab.2Positions of feature points序号1采样点1630序号5采样点8546序号9采样点14141通过式(5)计算得到理论故障周期l。=1 3 8 8。通过采样点估算的周期如表3 所示（第1 3 个特征点位置20004000采样点数1411614121086420510244066995910182636000800015t/s5773711343111967610000为第2 2 450 个采样点）。可以发现，估算的故障周期大多围绕1 3 8 8 这一理论故障周期上下浮动，同时存在一些更大的周期，但都接近1 3 8 8 的2 倍或3 倍，这是因20325471578127301221063142为硬阈值设置为8 后，一些特征点的小波系数因为噪声的影响没能达到硬阈值,因此被忽略。所以,提取出的故障特征点是有效的。表3 根据特征点的位置估算的故障周期Tab.3Estimates failure period based onlocation of characteristic points序号1估算周期2776序号5估算周期1413序号9估算周期4122在故障信号中，利用故障特征点分别向前、向后扩展1 49 个采样点得到j个截取信号段W（n,），任意选取一段为基准信号,如图8 所示,选取s=1为基准信号，即W（n）为基准信号。-1-20通过对信号段W（n,）分别向前、向后扩展1 49 个采样点得到信号段W（n,），并从信号W（n）中连续截取2 9 8 个采样点的长度得到信号段W（n），将T,设置为W（n）在W（n,）上的偏移量,W（n）偏移后的信号记录为W（n,+,）。T,的子集表示为L,L为(O：298)。通过计算W（n）与（n,+,）的相关性得到（n;+T）中偏移量,的最佳偏移量,如表3 所示。表4最佳偏移量Tab.4Optimum offsetj11j51779131液压与气动根据W（n,+T,）最佳偏移量,可以确定最佳的起始叠加点为（n,+,），叠加长度为基准信号的长度，即2 9 8 个采样点。利用起始叠加点（n;+）与叠加长度2 9 8 个采样点得到W（n;+T），选取连续的2 K个w（n,+)信号段，根据序号i的奇偶性将w（n,+)分成奇数组和偶数组，并对每组按照序号由低到2313671384671384138710111413138750100图8 选取的基准信号Fig.8 Selected reference signal2148615010171第48 卷第2 期4高进行不同次数的叠加运算,并计算奇数组与偶数组叠加后生成的叠加信号的相关系数p（K），其随叠加次1389数增加的变化曲线，如图9 所示。80.9514110.90120.8513870.8020.750.700.650.600图9 随叠加次数的变化曲线Fig.9Curve with number of superposition150200采样点数3131715211149一偶数组叠加信号一奇数组叠加信号12K25030041748152121513由图可知，当叠加次数K=3时，p（3）=0.9 0 9 5，大于设置的阈值M,但p（4）=0.8 7 6 8,故叠加5次的叠加信号为齿轮磨损故障成分，p（5）=0.9 2 7 0,如图10所示，其中实线为奇数组叠加信号，虚线为偶数组叠加信号。一偶数组叠加信号2F一奇数组叠加信号，/0-2E05结论本研究围绕旋转机械故障诊断的课题,以同步液压马达中的齿轮这一典型的旋转零部件为研究对象，提出了一种基于连续小波变换与错位叠加法的旋转机械冲击故障声信号的降噪提取方法，与改进的DSM相比,该方法摆脱了编码器的限制,提高了DSM的适用性,并利用该方法实现了对同步液压马达齿轮磨损状态的冲击故障声信号的降噪提取。450100图1 0 获取的故障成分Fig.10Obtained fault components5150200采样点数2503002024年第2 期对于已完成的工作，仍然存在一定的改进空间,可以从以下两个方面进行进一步的研究：（1）同步液压马达的型号繁多,如果条件允许，可以选择更为大型的同步液压马达进行研究。（2）不同工况和不同冲击故障下,该方法中硬阈值的设置也不尽相同,所提取出的故障成分也有差别，如何设置一个合适的硬阈值是该方法需要改进的地方。参考文献：1 NING D,CONG Y.Shocking Fault Component of AbnormalSound Signal in the Fault Engine Extract Method Based onLinear Superposition Method and Cross-Correlation AnalysisJ.Advances in Mechanical Engineering,2015,7(8):1-9.2GLOWACZ A.Fault Diagnosis of Single-phase Induction MotorBased on Acoustic Signals J.Mechanical Systems andSignal Processing,2019,117:65-80.3CARLETTI E,MICCOLI G,PEDRIELLI F,et al.Vibroac-oustic Measurements and Simulations Applied to External GearPumps.An Integrated Simplified Approach J.Archives ofAcoustics,2016,41(2):285-296.4王广逸，龚清阳.汽车发动机机械故障的非接触式检测技术研究 J.内燃机与配件,2 0 2 2，（2）：1 3 7-1 3 9.WANG Guangyi,CONG Qingyang.Research on Non-contactDetection Technology of Automobile Engine MechanicalFault J.Internal Combustion Engine&Parts,2022,(2):137 139.5杜正昱，马洁.基于WATVFEMD-SAM的轴承声信号故障诊断 J.组合机床与自动化加工技术，2 0 2 1，（1 0）：2 3-27.DU Zhengyu,MA Jie.Fault Diagnosis of Bearing AcousticSignal Based on WATVFEMD-SAM J.Modular MachineTool&Automatic Manufacturing Technique,2021,(10):23-27.6汪欣，毛东兴，李晓东.基于声信号和一维卷积神经网络的电机故障诊断研究J.噪声与振动控制,2 0 2 1,41(2):125-129.WANG Xin,MAO Dong xing,LI Xiaodong.Motor FaultDiagnosis Using Microphones and One-dimensional Convol-引用本文：鲁晓，郭鹏伟，宁大勇，等.基于错位叠加法的同步液压马达冲击故障成分的降噪提取 J.液压与气动，2 0 2 4,48（2）：1 3 7 1 43.LU Xiao,GUO Pengwei,NING Dayong,et al.Denoising Extraction of Impact Fault Components of Synchronous Hydraulic Motor Basedon Displacement Superposition Method J.Chinese Hydraulics&Pneumatics,2024,48(2):137-143.液压与气动utional Neural Network J.Noise and Vibration Control,2021,41(2):125 129.7YAO Jiachi,LIU Chao,SONG Keyu,et al.Fault Diagnosisof Planetary Gearbox Based on Acoustic Signals J.Applied Acoustics,2021,(181):108151.8STEVE K,YU ZHANG.Acoustic-based Engine Fault Dia-gnosis Using WPT,PCA and Bayesian Optimization J.Applied Sciences,2020,(10):6890.9于红梅.基于深度自编码网络与模糊推理相结合的矿用齿轮箱故障诊断方法 J.机床与液压,2 0 2 0,48（9）：1 8 1-186.YU Hongmei.Fault Diagnosis Method of Mine Gear BoxBased on Deep Self-coding Network Combined with FuzzyReasoning J.Machine Tools and Hydraulics,2020,48(9):181 186.【1 0 武成同.基于声信号的风电机齿轮箱故障诊断研究 D.大连：大连海事大学，2 0 1 8.WU Chengtong.Research on Fault Diagnosis of WindTurbine Gear Box Based on Sound Signal D.Dalia:Dalian Maritime University,2018.11 孙鸿宇.基于声信号同步液压马达故障诊断方法及试验研究 D.大连：大连海事大学,2 0 2 0.SUN Hongyu.Fault Diagnosis Method and ExperimentalResearch of Synchronous Hydraulic Motor Based on AcousticSignal D.Dalian:Dalian Maritime University,2020.12 1NING Dayong,JIANG Yuhua,SUN Hongyu,et al.Separ-ation Method of Impulsive Fault Component for GasolineEngine Based on Acoustic Signal Analysis J.Shock andVibration,2019,(2019):1-15.13 NING D a y o n g,SU N Ch a n g l e,CO NG Yo n g j u n,e t a l.Extraction of Fault Component from Abnormal Sound inDiesel Engines Using Acoustic Signals J.MechanicalSystems and Signal Processing,2016,(75):544-555.14 GAO Xuejin,WEN Huanran,WANG Pu.A New Methodof Bearing Fault Diagnosis Based on LMD and WaveletDenoising C/2017 29th Chinese Control and DecisionConference Chongqing,2017:4155-4162.15HAI Qiu,JAY L,JING Lin,et al.Wavelet Filter-basedWeak Signature Detection Method and Its Application onRolling Element Bearing Prognostics J.Journal of Soundand Vibration,2006,289(4-5):1066-1090.143