CN111623982A - 一种基于apewt和imomeda的行星齿轮箱早期故障诊断方法 - Google Patents
一种基于apewt和imomeda的行星齿轮箱早期故障诊断方法 Download PDFInfo
- Publication number
- CN111623982A CN111623982A CN202010541496.9A CN202010541496A CN111623982A CN 111623982 A CN111623982 A CN 111623982A CN 202010541496 A CN202010541496 A CN 202010541496A CN 111623982 A CN111623982 A CN 111623982A
- Authority
- CN
- China
- Prior art keywords
- signal
- fault
- frequency
- spectrum
- scale space
- 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
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000003745 diagnosis Methods 0.000 title claims abstract description 21
- 238000001228 spectrum Methods 0.000 claims abstract description 102
- 238000000354 decomposition reaction Methods 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000012216 screening Methods 0.000 claims description 7
- 238000005096 rolling process Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 12
- 230000003044 adaptive effect Effects 0.000 description 7
- 230000000737 periodic effect Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 5
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 4
- 238000005192 partition Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 201000004569 Blindness Diseases 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000006247 magnetic powder Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000005669 field effect Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000009347 mechanical transmission Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
Images
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/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/148—Wavelet transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
- G06F17/153—Multidimensional correlation or convolution
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Acoustics & Sound (AREA)
- Computing Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种基于APEWT和IMOMEDA的行星齿轮箱早期故障诊断方法。首先,计算行星齿轮箱故障振动信号的离散尺度空间谱,在尺度空间谱上自适应确定频段划分边界;然后,采用APEWT将故障信号自动分解为一系列模态分量,并从中选取敏感分量;其次,采用IMOMEDA对敏感分量进行解卷积处理,并对解卷积后的分量进行波形延伸得到提升信号;最后,从提升信号的包络谱中提取故障特征频率,并与齿轮箱理论故障特征频率值对比,从而识别故障类型。本发明能够清晰准确地提取出行星齿轮箱早期故障特征信息,不仅解决了EWT中参数无法自适应确定的问题,还解决了MOMEDA存在的严重边缘效应的问题,提高了故障诊断的精确度。
Description
技术领域
本发明属于旋转机械故障诊断与振动信号分析技术领域,涉及一种旋转机械早期故障诊断方法,具体涉及一种基于自适应无参经验小波变换(AdaptiveParameterlessEmpirical Wavelet Transform,APEWT)和提升的多点优化最小熵解卷积修正(Improved Multipoint Optimal Minimum Entropy Deconvolution Adjusted,IMOMEDA)的行星齿轮箱早期故障诊断方法,可用于行星齿轮箱中齿轮和滚动轴承的故障诊断。
背景技术
行星齿轮箱具有体积小,传动比大、承载力强、运行平稳和工作效率高等优点,已被广泛应用于直升机、风力发电、重型卡车和舰船等大型复杂机械设备中。然而,行星齿轮箱的工作环境通常比较恶劣,长期运行在高负荷、强冲击和高污染的工况下,极易导致齿轮箱中的关键性部件(太阳轮、行星轮、齿圈和滚动轴承等)发生故障。而行星齿轮作为大型复杂机电装备机械传动系统的重要组成部分,一旦发生故障,将加速设备性能恶化,轻则引起整个动力传输系统失效,重则导致机毁人亡。因此,若能在齿轮箱失效早期阶段实现故障溯源并及时消除隐患,对保证设备稳定安全运行,避免灾难性事故的发生意义十分重大。
当行星齿轮箱内的齿轮或滚动轴承出现局部损伤时,在运转过程中损伤部位会产生周期性、非平稳的冲击振动,振动向外传递过程中会受到复杂传递路径的影响逐渐减弱,同时齿轮箱中多对齿轮啮合振动相互耦合,导致从壳体上采集到的振动响应信号表现为强背景噪声、非平稳、非线性和耦合调制等特征。时频分析方法是分析非平稳、非线性信号的有力工具,例如:Wigner–Ville分布、短时傅里叶变换、小波变换等。然而,这些方法都具有各自的局限性。Wigner–Ville分布用于多分量信号时存在严重的交叉干扰项;短时傅里叶变换具有固定时频分辨率的缺陷;小波变换需要预先确定小波基函数和分解水平,使其缺乏自适应性。因此,非自适应信号处理方法对于实际齿轮箱故障信号的分析难以取得满意的效果。经验模式分解(Empirical Mode Decomposition,EMD)方法作为较早使用的自适应信号分解方法,可以根据信号自身的局部尺度特征自适应地分解为一系列固有模态函数之和,从而揭示信号的内部本质。虽然EMD方法具有很多优势,但存在模态混叠、端点效应以及理论基础缺乏等问题。为了弥补EMD方法的不足,出现了集合经验模式分解(EnsembleEmpirical Mode Decomposition,EEMD)和噪声辅助分析的总体局部平均分解(EnsembleLocal Mean Decomposition,ELMD)等改进方法,但EEMD和ELMD中添加白噪声的两个关键参数(噪声的振幅和集成次数)选取机理尚未明确,同时添加噪声后在信号重构过程中存在残余噪声污染和运算量大等问题。不同于上述的分析方法,经验小波变换(EmpiricalWavelet Transform,EWT)是近几年提出的一种新颖的自适应信号分解算法,它能有效地分解信号中的固有模态,具有较低的运算复杂度与完备的数学理论基础而受到学者们的青睐。但EWT算法在应用时,需要依靠使用者来预先设置分解所需的模态个数对被分析信号的频谱进行划分;此外,基于频谱极大值的频段边界划分策略容易使该算法陷入局部最优,从而导致频谱划分过于密集,分解结果出现严重的模态混叠现象。
此外,由于行星轮箱早期故障特征比较微弱且时常被机械系统多干扰源和噪声所淹没,直接从EWT分解得到的模态分量中难以提取满意的故障特征信息。为了有效、准确地提取齿轮箱早期故障特征信息,需要进一步消除模态分量中的噪声和其他无关干扰成分,增强微弱故障特征。齿轮箱中故障源信号的传递过程可以看作是故障源信号与信号传输通道的一个线性卷积混合过程,提取故障原始信号则可以看作是一个解卷积的过程。多点优化最小熵解卷积修正(Multipoint Optimal Minimum Entropy Deconvolution Adjusted,MOMEDA)是一种新的非迭代盲解卷积增强技术,该算法引入时间目标向量来确定待解卷积脉冲序列的位置与权重,并应用多点峭度值确定故障发生的周期,从而实现对连续多点故障冲击脉冲的增强提取。然而,在MOMEDA算法中解卷积信号存在严重的边缘效应,特别是当滤波器长度较大时,会造成信号中重要的故障信息丢失,从而限制MOMEDA的应用。
发明内容
鉴于上述问题,本发明的目的在于提供一种基于自适应无参经验小波变换和提升多点优化最小熵解卷积修正(APEWT-IMOMEDA)的行星齿轮箱早期故障诊断方法,针对强背景噪声下行星齿轮箱早期故障特征微弱难以提取、传统EWT算法中信号频谱不能自适应划分和MOMEDA算法中存在严重的边缘效应等问题。
本发明的技术方案:
一种基于APEWT和IMOMEDA的行星齿轮箱早期故障诊断方法,包括如下步骤:
步骤S1:将加速度传感器安装在行星齿轮箱壳体的正上方,用于采集故障振动信号x(n)(n=1,2,...,N),N为信号的数据点数;
步骤S2:对行星齿轮箱故障信号x(n)进行快速离散傅里叶变换得到信号频谱X(f),将高斯核函数与信号频谱X(f)进行卷积运算得到其离散尺度空间谱L(f,s);
其中,离散尺度空间谱的计算过程为:将步骤S1中采集到的故障振动信号进行快速离散傅里叶变换得到其频谱X(f),并与高斯核函数g(f,s)进行卷积运算得到尺度空间谱,则信号傅里叶频谱的离散尺度空间表示为:
式中,s为尺度参数;f代表频率;高斯核函数为P为高斯核函数的长度,取以使其逼近误差能够被忽略不计;如果尺度参数设定为固定值s0,那么尺度空间L(f,s0)的每个局部极大值对应一个单分量成分,该单分量的频段由相邻的两个局部极小值确定;选取尺度参数fch为故障特征频率;
步骤S3:在尺度空间谱L(f,s)上自适应确定频段的划分边界,具体过程为:首先搜索离散尺度空间谱中的局部极大值点,将极大值点中的局部极小值点作为尺度空间谱划分的初始边界点,对离散尺度空间谱进行初始频段划分;假设初始边界将离散尺度空间谱划分为M段,所对应的初始频率边界为ωinitial={ωi,i=0,1...M},其中ω0=0,ωM=fs/2,fs为信号的采样频率;然后,在初始划分的每个频段区间内,搜索幅值最大值点Lmax(n),其对应的频率值为ωmax(f),并且在每个频率区间内最大值点的左侧搜索极小值点Llmin(n),其对应的频率值为ωlmin(f);在最大值点的右侧搜索极小值点Lrmin(n),其对应的频率值ωrmin(f);其次,通过spline三次样条插值方法拟合离散尺度空间谱上所获取的极大值点ωmax和极小值点ωlmin及ωrmin,从而得到离散尺度空间谱的变化曲线Y(k)(k=1,2,...,N/2);最后,通过选取变化曲线Y(k)上的局部极小值点所对应尺度空间谱上的频率值(其中)作为尺度空间谱的最终频段划分的边界点;
步骤S4:由于尺度空间谱与信号傅里叶频谱的频率轴相互对应,基于步骤S3划分得到的频段边界,在信号傅里叶频谱中构建正交小波滤波器组,采用APEWT算法将行星齿轮箱故障振动信号分解为一系列模态分量z(n);
步骤S5:根据敏感分量筛选准则,从APEWT分解的模态分量中选取敏感分量zmax(n)作进一步的分析;
敏感分量筛选准则为:计算APEWT分解得到的各模态分量与原始故障信号之间的余弦相似度值CS,选取最大CS值对应的模态作为敏感分量进行下一步的处理;其中,余弦相似度值CS的计算公式如下:
式中,X1(n)和X2(n)分别代表不同的信号;
步骤S6:采用IMOMEDA算法对敏感分量zmax(n)进行解卷积处理,从而实现微弱故障特征的增强;
IMOMEDA算法输入参数具体为:滤波器长度设置为1500;选取长度为3的矩形窗作为窗函数;故障周期搜索范围设置为[0,fs/2],迭代步长为0.2;
对解卷积后信号进行波形延伸的具体策略为:找解卷积后信号左侧的局部极值点及其位置,确定对称中心,对信号进行波形延伸。
确定信号延伸的对称中心的具体方法为:设a、b和c分别为信号的前三个局部极值点,若a≥c>b,则选取极值点a作为对称中心;若c>a≥b,则选取极值点b作为对称中心;若b>c≥a,则选取极值点a作为对称中心;若b>a≥c,则选取极值点b作为对称中心。
对信号进行波形延伸具体方法为:在确定对称中心后,选取对称中心后面一定数量的样本点作为扩展对象,然后取对称中心的镜像,得到扩展的信号波形;
本发明与现有技术的优势在于具有以下特点:
(1)本发明不仅发挥了APEWT算法自适应分解信号得到单分量信号的优势,而且还利用IMOMEDA算法能够提取周期性脉冲序列的特性,解决了传统方法在行星齿轮箱早期微弱故障特征提取困难的问题。
(2)本发明的APEWT算法中频谱边界划分技术保证了频谱划分时各频段区间内包含着一个局部极大值,同时频段的边界由尺度空间谱变化曲线上的两个局部极小值所对应的频率确定,这满足了频谱区间边界划分的要求,能够自适应准确地划分出具有调制特性的频率区间,从而分离出含有故障冲击特征的模态分量;
(3)本发明在尺度空间谱上自适应划分APEWT所需的频谱边界,有效地克服了传统EWT算法划分频谱边界时出现频段划分过窄,以及分解的信号中存在模态混叠的问题。
(4)本发明采用一种新的敏感分量筛选准则,利用余弦相似度指标在APEWT算法分解得到的一系列模态分量信号中自动选取包含丰富故障冲击特征的分量作为敏感分量进行后续的解卷积增强处理,从而避免了模态分量选取的盲目性和主观性。
(5)本发明针对原始MOMEDA算法存在的边缘效应,提出了一种波形延伸策略能够对解卷积信号进行自适应补偿,使其恢复到与原信号相同的长度,有效地克服了边缘效应的影响,保留了信号中重要的微弱特征信息,从而显著提高了MOMEDA的解卷积增强性能。
附图说明
图1是本发明诊断方法流程图;
图2是本发明实施例中行星齿轮箱故障实验台示意图;
图中,1抗振底座、2驱动电动机、3弹性联轴器、4可拆卸轴承座、5行星齿轮箱、6刚性联轴器、7磁粉加载器。
图3是本发明实施列中太阳轮齿根裂纹故障振动信号的时域波形及其傅里叶频谱;
图中,(a)为故障振动信号的时域波形、(b)为故障振动信号的傅里叶频谱。
图4是本发明实施列中太阳轮故障信号的尺度空间谱表示;
图5是本发明实施列中太阳轮故障信号尺度空间谱中频段边界划分示意图;
图6是本发明实施列中太阳轮故障信号的傅里叶频谱中频段边界划分示意图;
图7是本发明实施列中采用APEWT方法对太阳轮故障信号进行自适应分解的结果;
图中,(a)分解的IMF1分量、(b)分解的IMF2分量、(c)分解的IMF3分量、(d)分解的IMF4分量、(e)分解的IMF5分量、(f)分解的IMF6分量。
图8是本发明实施列中各模态分量与原始太阳轮故障信号之间余弦相似度的直方图;
图9是本发明实施列中敏感分量的直接包络谱;
图10是本发明实施列中在不同情况下的信号波形延伸示策略意图;
图中,(a)当a≥c>b时的情况、(b)当c>a≥b时的情况、(c)当b>c≥a时的情况、(d)当b>a≥c时的情况。
图11是本发明实施列中采用IMOMEDA方法对敏感分量进行解卷积增强后的结果;
图中,(a)敏感分量解卷积后信号波形延伸结果、(b)提升信号的包络谱。
具体实施方式
以下结合附图和实施例对本发明的具体实施过程做进一步说明。
本发明提供一种基于APEWT-IMOMEDA的行星齿轮箱早期故障诊断方法,诊断方法的具体流程如图1所示,其主要步骤如下:
步骤S1:将加速度传感器安装在行星齿轮箱壳体的正上方采集故障振动信号x(n)(n=1,2,...,N)(N为信号的数据点数);
步骤S2:对行星齿轮箱故障信号x(n)进行快速离散傅里叶变换得到信号的频谱X(f),将高斯核函数与信号频谱X(f)进行卷积运算得到其离散尺度空间谱L(f,s);
步骤S3:在尺度空间谱L(f,s)上自适应确定频段的划分边界;
步骤S4:由于尺度空间谱与信号傅里叶频谱的频率轴相互对应,基于步骤S3划分得到的频段边界,在傅里叶频谱中构建正交小波滤波器组,采用APEWT算法将行星齿轮箱故障振动信号分解为一系列模态分量z(n);
步骤S5:根据敏感分量筛选准则,从APEWT分解的模态分量中选取敏感分量zmax(n)作进一步的分析;
步骤S6:采用IMOMEDA算法对敏感分量zmax(n)进行解卷积处理,从而实现微弱故障特征的增强;
实施例1:
搭建行星齿轮箱故障诊断实验台,如图2所示。该实验台主要是由1抗振底座、2驱动电动机、3弹性联轴器、4可拆卸轴承座、5行星齿轮箱、6刚性联轴器、7磁粉加载器等组成。其中行星齿轮箱的结构参数见表1。通过线切割技术在太阳轮的某个轮齿上沿着齿根方向加工宽为0.15mm,深为1mm的微小贯通裂纹作为故障。实验时将故障太阳轮安装在行星齿轮箱内,进行实验数据采集。加速度传感器安装在行星齿轮箱外壳正上方的测点上。电机输入转速为1200r/min(转频为20Hz),待电机运行平稳后,进行试验数据采集,数据采样频率设置为5120Hz,采样时长为2s。在该实验条件下,可以计算得到行星齿轮箱中的特征频率如表2所示。
表1行星齿轮箱中齿轮参数(单位/个)
本发明基于APEWT-IMOMEDA的行星齿轮箱早期故障诊断方法流程如图1所示,具体包括以下步骤:
步骤一:利用加速度传感器采集所述故障诊断实验台中行星齿轮箱太阳轮故障振动信号。图3为垂直正上方加速度传感器获取到的太阳轮齿根裂纹故障信号及其傅里叶频谱。从图3(a)故障信号的时域波形中难以观察到有规律的周期性冲击特征,同时,图3(b)故障信号的频谱中太阳轮故障特征频率及倍频处无明显的峰值,噪声与其他无关干扰频率的幅值较为突出。因此,通过传统的时域和频域分析方法无法从原始故障信号中提取表征太阳轮健康状态的特征信息。
步骤二:将步骤一中采集到的故障振动信号进行快速傅里叶变换得到信号频谱X(f),将高斯核函数与信号频谱X(f)进行卷积运算得到其离散尺度空间谱表示L(f,s),则故障信号傅里叶频谱的离散尺度空间表示为:
式中,s为尺度参数;f代表频率;高斯核函数为(其中e为自然底数);P为高斯核函数的长度,取以使其逼近误差能够被忽略不计;如果尺度参数设定为固定值s0,那么尺度空间L(f,s0)的每个局部极大值对应一个单分量成分,该单分量的频段由相邻的两个局部极小值确定;本发明选取尺度参数s=9fch 2,fch为故障特征频率;图4为在该尺度参数下得到的尺度空间谱。
步骤三:在尺度空间谱上自适应划分频率的边界,如图5所示。具体步骤如下:
(1)搜索尺度空间谱中的局部极大值点,在获取的极大值点中寻找局部极小值点作为尺度空间谱划分的初始边界点,对尺度空间谱进行初始频段划分。假设初始边界将尺度空间谱划分为M段,所对应的初始频率边界为ωinitial={ωi,i=0,1...M}(其中ω0=0,ωM=fs/2,fs为信号的采样频率),如图5中的虚线所示。
(2)在初始划分的每个频段区间内,搜索幅值最大值点Lmax(n),其对应的频率值为ωmax(f)(如图5中▲所示),并且在每个频率区间内最大值点的左侧搜索极小值点Llmin(n),其对应的频率值为ωlmin(f)(如图5中所示);在最大值点的右侧搜索极小值点Lrmin(n),其对应的频率值ωrmin(f)(如图5中所示)。
(3)通过spline三次样条插值方法拟合尺度空间谱上所获取的极大值点ωmax和极小值点ωlmin及ωrmin,从而得到尺度空间谱的变化曲线Y(k)(k=1,2,...,N/2),如图5中点线所示。最后,通过选取变化曲线Y(k)上的局部极小值点所对应尺度空间谱上的频率值(其中)作为尺度空间谱的最终频段划分的边界点,如图5中点划线所示;
(4)进一步地,如果在[ω0,ω1]之间有两个或者多个新频率边界ωnew,选择最接近ω1的新边界作为最终频谱边界划分的第一个频率边界,则最终频率划分边界可表示为否则,最终的频率划分边界ωfinal=ωnew。
步骤四:由于尺度空间谱与信号傅里叶频谱的频率轴相互对应,将步骤三尺度空间谱上划分得到的频段边界ωfinal平移到信号傅里叶频谱上。图6所示为原始故障信号傅里叶频谱上的频率划分情况。基于信号傅里叶频谱上的各频率区间构建经验小波变换的带通滤波器组,根据Littlewood-Paley和Meyer小波的构造思路,对于任何n>0,其经验尺度函数和经验小波函数分别如下列公式:
式中,β(x)是定义在[0,1]上的函数且满足β(x)+β(1-x)=1,这里β(x)=x4(35-84x+70x2-20x3);为使得为紧支撑结构,需要根据上述频率边界,同经典的小波变换方法,可以求得细节系数和近似系数如下式:
图7为采用APEWT算法对太阳轮故障振动信号进行自适应分解获得的结果;其中,图7(a)为IMF1分量、图7(b)为IMF2分量、图7(c)为IMF3分量、图7(d)为IMF4分量、图7(e)为IMF5分量、图7(a)为IMF6分量。
步骤五:计算各模态分量与原始故障信号之间的余弦相似度CS,选取最大CS值对应的模态作为敏感分量进行下一步的处理。其中,余弦相似度CS的计算公式如下:
式中,X1(n)和X2(n)分别代表不同的信号。
如图8所示为各模态分量信号与原始故障信号之间的余弦相似度CS的直方图,可以看到,第6个模态分量的CS值最大,因此将其作为敏感分量zmax(n)。敏感分量的直接包络谱如图9所示。由图9可知,在敏感分量的包络谱中3倍的行星架转频3fc占据主导地位,虽然在太阳轮故障特征率及其倍频处存在峰值,但受噪声和其他干扰频率的影响,故障特征频率并不明显且幅值较低,据此很难准确判断故障类型。
步骤六:采用IMOMEDA算法对敏感分量zmax(n)进行解卷积处理,从而实现微弱故障特征增强;其中,IMOMEDA算法输入参数具体为:滤波器长度设置为1500;选取长度为3的矩形窗作为窗函数;故障周期搜索范围设置为[0,fs/2],迭代步长为0.2。
MOMDEA算法的具体步骤如下:
(1)计算信号多点D-范数,使其取得最大值:
通过求解多点D-范数的最大值使得解卷积效果达到最佳,即
tn=Pn(T)=δround(T)+δround(2T)+...+δround(nT)
式中,δ为样本n的一个脉冲。在实际应用中,由于故障周期不一定是采样周期的倍数,因此非整数T应该进行四舍五入取整数值。
步骤七:对步骤六中解卷积后的敏感分量进行波形延伸的具体策略为:
(1)寻找信号的极值点;寻找解卷积后信号左侧的局部极值点及其位置。
(2)确定信号延伸的对称中心;根据信号在左边界的局部特征确定信号延伸时的对称中心,设a、b和c分别为信号的前三个局部极值点,若a≥c>b,如图10(a)所示,则选取极值点a作为对称中心;若c>a≥b,如图10(b)所示,则选取极值点b作为对称中心;若b>c≥a,如图10(c)所示,则选取极值点a作为对称中心;若b>a≥c,如图10(d)所示,则选取极值点b作为对称中心。
(3)信号波形延伸;选取对称中心后面一定数量的样本点作为扩展对象,然后取对称中心的镜像,得到扩展的信号波形。
根据上述信号波形延伸策略,敏感分量解卷积后的波形延伸结果如图11(a)中点线所示;信号波形延伸策略有效地解决了MOMEDA解卷积信号存在的严重边缘效应,使信号长度与原始信号相等,从而更好地保留了信号中重要的微弱特征信息。
步骤八:对步骤七中获得的提升信号进行包络解调分析;由图11(a)可知,通过IMOMEDA解卷积处理后,隐藏在动态信号中的周期性冲击特征被清晰地揭露出来,噪声及其他干扰成分被有效地抑制;同时,在图11(b)提升信号的包络谱中太阳轮故障特征频率fs及其倍频nfs(n=1,2,3,...,5)处存在较为明显的谱线,且在整个频谱中占主导地位。因此,上述分析结果表明太阳轮出现了局部故障,这与实验设置相一致。从而验证了本发明基于APEWT-IMOMEDA的早期故障诊断方法的有效性和可行性。
综上所述,本发明基于APEWT-IMOMEDA的早期故障诊断方法,能够成功地提取微弱的故障特征信息,实现行星齿轮箱早期故障的准确识别与诊断。其中,APEWT算法具有无参数化和自适应性,能实现信号频谱的自适应划分,无需使用者输入任何参数即可对复杂的行星齿轮箱故障信号进行自动分解,最终获取一系模态分量,分解结果有效地避免了模态混叠现象。依据敏感分量筛选准则,选取包含故障信息最多的模态分量作为研究对象,减少了分量选取的盲目性。此外,对敏感分量进行IMOMEDA解卷积处理,能够消除强背景噪声及其他无关成分的干扰,增强微弱的周期性故障冲击特征,同时也提高了行星齿轮箱早期故障的表征能力与诊断精度。根据本发明中IMOMEDA波形延伸策略,对解卷积后的敏感模态进行信号波形延伸可对缺失的信息进行自适应补偿,从而有效地克服了MOMEDA算法解卷积信号边缘效应的影响,提升了MOMEDA对微弱故障冲击成分的解卷积增强性能。
以上所述仅为本发明的较好实施例,凡依本发明申请专利范围所做的均等变化修饰,皆应属本发明的涵盖范围。
Claims (1)
1.一种基于APEWT和IMOMEDA的行星齿轮箱早期故障诊断方法,其特征在于,包括如下步骤:
步骤S1:将加速度传感器安装在行星齿轮箱壳体的正上方,用于采集故障振动信号x(n)(n=1,2,...,N),N为信号的数据点数;
步骤S2:对行星齿轮箱故障信号x(n)进行快速离散傅里叶变换得到信号频谱X(f),将高斯核函数与信号频谱X(f)进行卷积运算得到其离散尺度空间谱L(f,s);
其中,离散尺度空间谱的计算过程为:将步骤S1中采集到的故障振动信号进行快速离散傅里叶变换得到其频谱X(f),并与高斯核函数g(f,s)进行卷积运算得到尺度空间谱,则信号傅里叶频谱的离散尺度空间表示为:
式中,s为尺度参数;f代表频率;高斯核函数为P为高斯核函数的长度,取以使其逼近误差;如果尺度参数设定为固定值s0,那么尺度空间L(f,s0)的每个局部极大值对应一个单分量成分,该单分量的频段由相邻的两个局部极小值确定;选取尺度参数s=9fch 2,fch为故障特征频率;
步骤S3:在尺度空间谱L(f,s)上自适应确定频段的划分边界,具体过程为:首先搜索离散尺度空间谱中的局部极大值点,将极大值点中的局部极小值点作为尺度空间谱划分的初始边界点,对离散尺度空间谱进行初始频段划分;假设初始边界将离散尺度空间谱划分为M段,所对应的初始频率边界为ωinitial={ωi,i=0,1...M},其中ω0=0,ωM=fs/2,fs为信号的采样频率;然后,在初始划分的每个频段区间内,搜索幅值最大值点Lmax(n),其对应的频率值为ωmax(f),并且在每个频率区间内最大值点的左侧搜索极小值点Llmin(n),其对应的频率值为ωlmin(f);在最大值点的右侧搜索极小值点Lrmin(n),其对应的频率值ωrmin(f);其次,通过spline三次样条插值方法拟合离散尺度空间谱上所获取的极大值点ωmax和极小值点ωlmin及ωrmin,从而得到离散尺度空间谱的变化曲线Y(k)(k=1,2,...,N/2);最后,通过选取变化曲线Y(k)上的局部极小值点所对应尺度空间谱上的频率值作为尺度空间谱的最终频段划分的边界点,其中
步骤S4:由于尺度空间谱与信号傅里叶频谱的频率轴相互对应,基于步骤S3划分得到的频段边界,在信号傅里叶频谱中构建正交小波滤波器组,采用APEWT算法将行星齿轮箱故障振动信号分解为一系列模态分量z(n);
步骤S5:根据敏感分量筛选准则,从APEWT分解的模态分量中选取敏感分量zmax(n)作进一步的分析;
敏感分量筛选准则为:计算APEWT分解得到的各模态分量与原始故障信号之间的余弦相似度值CS,选取最大CS值对应的模态作为敏感分量进行下一步的处理;其中,余弦相似度值CS的计算公式如下:
式中,X1(n)和X2(n)分别代表不同的信号;
步骤S6:采用IMOMEDA算法对敏感分量zmax(n)进行解卷积处理,从而实现微弱故障特征的增强;
IMOMEDA算法输入参数具体为:滤波器长度设置为1500;选取长度为3的矩形窗作为窗函数;故障周期搜索范围设置为[0,fs/2],迭代步长为0.2;
对解卷积后信号进行波形延伸的具体策略为:找解卷积后信号左侧的局部极值点及其位置,确定对称中心,对信号进行波形延伸;
确定信号延伸的对称中心的具体方法为:设a、b和c分别为信号的前三个局部极值点,若a≥c>b,则选取极值点a作为对称中心;若c>a≥b,则选取极值点b作为对称中心;若b>c≥a,则选取极值点a作为对称中心;若b>a≥c,则选取极值点b作为对称中心;
对信号进行波形延伸具体方法为:在确定对称中心后,选取对称中心后面一定数量的样本点作为扩展对象,然后取对称中心的镜像,得到扩展的信号波形;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010541496.9A CN111623982B (zh) | 2020-06-15 | 2020-06-15 | 一种基于apewt和imomeda的行星齿轮箱早期故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010541496.9A CN111623982B (zh) | 2020-06-15 | 2020-06-15 | 一种基于apewt和imomeda的行星齿轮箱早期故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111623982A true CN111623982A (zh) | 2020-09-04 |
CN111623982B CN111623982B (zh) | 2021-03-26 |
Family
ID=72270404
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010541496.9A Active CN111623982B (zh) | 2020-06-15 | 2020-06-15 | 一种基于apewt和imomeda的行星齿轮箱早期故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111623982B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112098093A (zh) * | 2020-09-15 | 2020-12-18 | 丽水市特种设备检测院 | 一种轴承故障特征识别方法及识别系统 |
CN112284719A (zh) * | 2020-09-16 | 2021-01-29 | 北京信息科技大学 | 一种行星齿轮箱齿面磨损故障诊断方法及系统 |
CN112539933A (zh) * | 2020-12-16 | 2021-03-23 | 东南大学 | 一种基于共振基带宽傅立叶分解的齿轮箱故障诊断方法 |
CN113176460A (zh) * | 2021-04-28 | 2021-07-27 | 河海大学 | 一种基于改进经验小波变换的电能质量扰动信号检测方法 |
CN113310693A (zh) * | 2021-06-07 | 2021-08-27 | 华润电力技术研究院有限公司 | 一种机械故障检测方法、装置、设备及存储介质 |
CN113537102A (zh) * | 2021-07-22 | 2021-10-22 | 深圳智微电子科技有限公司 | 一种微震信号的特征提取方法 |
CN114112397A (zh) * | 2021-11-24 | 2022-03-01 | 郑州恩普特科技股份有限公司 | 一种滚动轴承故障诊断方法 |
CN114216676A (zh) * | 2021-11-30 | 2022-03-22 | 上海海事大学 | 一种时变工况下无转速计的行星齿轮箱复合故障诊断方法 |
CN114528525A (zh) * | 2022-01-11 | 2022-05-24 | 西南交通大学 | 基于最大重加权峭度盲解卷积的机械故障诊断方法 |
CN114757242A (zh) * | 2022-06-16 | 2022-07-15 | 中国空气动力研究与发展中心低速空气动力研究所 | 基于循环维纳滤波的直升机噪声增强方法以及检测方法 |
CN114894468A (zh) * | 2022-04-07 | 2022-08-12 | 大连理工大学 | 一种基于混沌检测的行星齿轮箱早期微弱故障诊断方法 |
CN115931358A (zh) * | 2023-02-24 | 2023-04-07 | 沈阳工业大学 | 一种低信噪比的轴承故障声发射信号诊断方法 |
CN116520095A (zh) * | 2023-07-03 | 2023-08-01 | 昆明理工大学 | 故障测距方法、系统以及计算机可读存储介质 |
CN117928951A (zh) * | 2024-03-21 | 2024-04-26 | 安徽大学 | 基于改进经验小波变换和包络谱能量比的故障诊断方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20180112258A (ko) * | 2017-04-03 | 2018-10-12 | 두산인프라코어 주식회사 | 선회장치의 고장진단 장치 및 방법 |
CN109187024A (zh) * | 2018-09-04 | 2019-01-11 | 温州大学激光与光电智能制造研究院 | 一种活塞式空压机曲轴箱滚动轴承故障诊断方法 |
CN109324604A (zh) * | 2018-11-29 | 2019-02-12 | 中南大学 | 一种基于多源信号的智能列车综合故障分析方法 |
CN110470475A (zh) * | 2019-09-04 | 2019-11-19 | 中国人民解放军空军工程大学航空机务士官学校 | 一种航空发动机中介轴承早期微弱故障诊断方法 |
CN111141520A (zh) * | 2020-02-24 | 2020-05-12 | 江南大学 | 一种基于改进经验小波变换的滚动轴承故障诊断方法 |
-
2020
- 2020-06-15 CN CN202010541496.9A patent/CN111623982B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20180112258A (ko) * | 2017-04-03 | 2018-10-12 | 두산인프라코어 주식회사 | 선회장치의 고장진단 장치 및 방법 |
CN109187024A (zh) * | 2018-09-04 | 2019-01-11 | 温州大学激光与光电智能制造研究院 | 一种活塞式空压机曲轴箱滚动轴承故障诊断方法 |
CN109324604A (zh) * | 2018-11-29 | 2019-02-12 | 中南大学 | 一种基于多源信号的智能列车综合故障分析方法 |
CN110470475A (zh) * | 2019-09-04 | 2019-11-19 | 中国人民解放军空军工程大学航空机务士官学校 | 一种航空发动机中介轴承早期微弱故障诊断方法 |
CN111141520A (zh) * | 2020-02-24 | 2020-05-12 | 江南大学 | 一种基于改进经验小波变换的滚动轴承故障诊断方法 |
Non-Patent Citations (4)
Title |
---|
ZHI-WU SHANG 等: "Rolling Bearing Fault Diagnosis Method Based on MOMEDA and IEWT", 《INTERNATIONAL JOURNAL OF INFORMATION AND MANAGEMENT SCIENCES》 * |
ZIYING ZHANG 等: "Gearbox Composite Fault Diagnosis Method Based on Minimum Entropy Deconvolution and Improved Dual-Tree Complex Wavelet Transform", 《ENTROPY》 * |
任学平 等: "基于VMD和MED的滚动轴承微弱故障特征提取", 《现代制造工程》 * |
郑近德 等: "自适应无参经验小波变换及其在转子故障诊断中的应用", 《中国机械工程》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112098093A (zh) * | 2020-09-15 | 2020-12-18 | 丽水市特种设备检测院 | 一种轴承故障特征识别方法及识别系统 |
CN112284719A (zh) * | 2020-09-16 | 2021-01-29 | 北京信息科技大学 | 一种行星齿轮箱齿面磨损故障诊断方法及系统 |
CN112539933A (zh) * | 2020-12-16 | 2021-03-23 | 东南大学 | 一种基于共振基带宽傅立叶分解的齿轮箱故障诊断方法 |
CN113176460A (zh) * | 2021-04-28 | 2021-07-27 | 河海大学 | 一种基于改进经验小波变换的电能质量扰动信号检测方法 |
CN113310693A (zh) * | 2021-06-07 | 2021-08-27 | 华润电力技术研究院有限公司 | 一种机械故障检测方法、装置、设备及存储介质 |
CN113537102A (zh) * | 2021-07-22 | 2021-10-22 | 深圳智微电子科技有限公司 | 一种微震信号的特征提取方法 |
CN114112397A (zh) * | 2021-11-24 | 2022-03-01 | 郑州恩普特科技股份有限公司 | 一种滚动轴承故障诊断方法 |
CN114216676A (zh) * | 2021-11-30 | 2022-03-22 | 上海海事大学 | 一种时变工况下无转速计的行星齿轮箱复合故障诊断方法 |
CN114528525A (zh) * | 2022-01-11 | 2022-05-24 | 西南交通大学 | 基于最大重加权峭度盲解卷积的机械故障诊断方法 |
CN114528525B (zh) * | 2022-01-11 | 2023-03-28 | 西南交通大学 | 基于最大重加权峭度盲解卷积的机械故障诊断方法 |
CN114894468A (zh) * | 2022-04-07 | 2022-08-12 | 大连理工大学 | 一种基于混沌检测的行星齿轮箱早期微弱故障诊断方法 |
CN114894468B (zh) * | 2022-04-07 | 2023-01-03 | 大连理工大学 | 一种基于混沌检测的行星齿轮箱早期微弱故障诊断方法 |
CN114757242A (zh) * | 2022-06-16 | 2022-07-15 | 中国空气动力研究与发展中心低速空气动力研究所 | 基于循环维纳滤波的直升机噪声增强方法以及检测方法 |
CN115931358A (zh) * | 2023-02-24 | 2023-04-07 | 沈阳工业大学 | 一种低信噪比的轴承故障声发射信号诊断方法 |
CN115931358B (zh) * | 2023-02-24 | 2023-09-12 | 沈阳工业大学 | 一种低信噪比的轴承故障声发射信号诊断方法 |
CN116520095A (zh) * | 2023-07-03 | 2023-08-01 | 昆明理工大学 | 故障测距方法、系统以及计算机可读存储介质 |
CN116520095B (zh) * | 2023-07-03 | 2023-09-12 | 昆明理工大学 | 故障测距方法、系统以及计算机可读存储介质 |
CN117928951A (zh) * | 2024-03-21 | 2024-04-26 | 安徽大学 | 基于改进经验小波变换和包络谱能量比的故障诊断方法 |
CN117928951B (zh) * | 2024-03-21 | 2024-06-11 | 安徽大学 | 基于改进经验小波变换和包络谱能量比的故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111623982B (zh) | 2021-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111623982B (zh) | 一种基于apewt和imomeda的行星齿轮箱早期故障诊断方法 | |
CN109029977B (zh) | 一种基于vmd-amckd的行星齿轮箱早期故障诊断方法 | |
CN107505135B (zh) | 一种滚动轴承复合故障提取方法及系统 | |
CN102998706B (zh) | 一种衰减地震数据随机噪声的方法及系统 | |
CN110598593B (zh) | 基于共振稀疏分解和FastICA算法的行星齿轮箱故障诊断方法 | |
CN110398364B (zh) | 基于共振稀疏分解和FastICA算法的行星齿轮箱故障诊断方法 | |
CN111507221B (zh) | 基于vmd和最大重叠离散小波包变换的齿轮信号去噪方法 | |
CN114216676B (zh) | 一种时变工况下无转速计的行星齿轮箱复合故障诊断方法 | |
CN112945546B (zh) | 一种齿轮箱复杂故障精确诊断方法 | |
KR102093929B1 (ko) | Cim 기반의 시스템 건전성 진단을 위한 장치 및 방법 | |
CN111638055B (zh) | 一种基于共振稀疏分解改进算法的齿轮箱故障诊断方法 | |
CN115493836B (zh) | 机械故障诊断方法及系统 | |
CN107966287B (zh) | 一种自适应机电装备微弱故障特征提取方法 | |
CN114881072B (zh) | 一种基于峰值包络谱的傅里叶分解信号降噪方法 | |
CN113074935A (zh) | 一种用于齿轮箱冲击故障特征的声学分离和诊断方法 | |
CN112098093A (zh) | 一种轴承故障特征识别方法及识别系统 | |
Wang et al. | Research on gearbox composite fault diagnosis based on improved local mean decomposition | |
CN107941511B (zh) | 一种基于信号时频分解的频率—峭度图的实现方法 | |
CN117370781A (zh) | 一种rv减速器复合缺陷特征提取方法、装置 | |
CN107610055B (zh) | 傅里叶变换光谱仪干涉图噪声检测及抑制方法 | |
CN112580451A (zh) | 一种基于改进emd和med的数据降噪方法 | |
CN116578840A (zh) | 一种基于eemd-isvd的滚动轴承振动信号降噪方法、装置、电子设备及介质 | |
CN111582248A (zh) | 一种基于svd的齿轮箱信号降噪方法 | |
CN113657268B (zh) | 一种应用于风电机组齿轮箱故障诊断的信号自动分解方法 | |
Zhang et al. | Adaptive singular value decomposition and its application to the feature extraction of planetary gearboxes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |