CN102879195B - 齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法 - Google Patents
齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法 Download PDFInfo
- Publication number
- CN102879195B CN102879195B CN201210361690.4A CN201210361690A CN102879195B CN 102879195 B CN102879195 B CN 102879195B CN 201210361690 A CN201210361690 A CN 201210361690A CN 102879195 B CN102879195 B CN 102879195B
- Authority
- CN
- China
- Prior art keywords
- noise reduction
- damage
- frequency coefficient
- row
- time domain
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法。该方法根据多小波系数之间的相关性,以Stein无偏风险估计误差最小作为约束条件,自适应地选取最优的邻域分块长度和阈值,能够有效消除噪声干扰,准确提取信号特征;利用所构造的自适应邻域分块阈值降噪方法分析连轧机传动齿轮箱振动信号,获得降噪后的重构信号,通过分析时域脉冲波形的周期性,有效提取出连轧机传动齿轮箱损伤特征。该检测方法结果可靠,实时性好,简单易行,通用性强,适用于连轧机等重载设备传动齿轮箱的损伤诊断。
Description
技术领域
本发明涉及机械设备故障诊断技术,特别涉及一种齿轮箱损伤的故障诊断方法。
背景技术
齿轮箱在运行中萌生的故障(即早期故障),具有症状不明显,特征信息微弱,往往被机械设备运行过程中的强背景噪声所淹没,从而使得传动齿轮箱早期故障的动态监测、诊断与故障预示的难度不断增大。因此,如何对采集到的振动信号进行降噪处理,突出或提取有用特征信息是故障诊断和故障预示中的一个关键课题。
基于相邻小波系数之间具有相关性的特点,现有的降噪处理方法一般为分块阈值降噪,它是将相邻的小波系数作为一个整体进行阈值处理。其不足之处是:(1)分块阈值降噪方法根据经验选择邻域分块长度,往往产生很大的偏差;(2)在小波分解的每一层均使用全局阈值,而全局阈值并不是最优的。
发明内容
本发明的目的在于提供一种可解决背景技术所存在的问题的齿轮箱损伤的时域诊断方法,该方法首先利用获得的最优邻域分块长度和阈值,对多小波分解后的高频系数进行阈值处理;其次把低频系数和阈值处理后的高频系数进行多小波反变换,得到重构信号;然后分析信号中时域脉冲波形的周期性和规律性,并结合齿轮箱机理分析,有效提取出齿轮箱损伤特征。
为达到以上目的,本发明是采取如下技术方案予以实现的:
一种齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法,其特征在于:
(1)对含噪信号进行预处理,得到矢量输入信号;
(2)对矢量输入信号进行多小波分解,得到高频系数和低频系数;
(3)根据判断条件对高频系数进行阈值处理:符合条件时采用自适应邻域分块阈值算法,否则采用逐点比较算法,得到降噪后的高频系数;
(4)对低频系数和经过阈值处理后的高频系数进行多小波反变换,重构降噪后的矢量输出信号;
(5)对矢量输出信号进行后处理,得到一维降噪结果;
(6)通过分析降噪结果中的时域脉冲波形的周期性,提取出传动齿轮箱的损伤特征进行故障诊断。
其中,步骤3)所述的自适应邻域分块阈值算法包括下述步骤:
1)设含噪信号S是所关心的故障特征F与噪声E的叠加,可以表示为:
S=F+E
以2×2维多小波变换为例,S=[s(1),s(2)]T表示经过预处理后的2维振动信号,噪声E服从多元正态分布N(0,Vj),在第k个邻域分块中,定义变量作为阈值降噪的基本变量,将其中向量化的系数转化为标量,矩阵Vj为多小波分解第j层的高频系数的协方差矩阵;
定义:
mad(y)=1.4826·median(|y-median(y)|)
式中,函数median()用于寻找数列y中的中位数,矩阵Vj可表达为:
式中,a1=1/mad(row1),a2=1/mad(row2),b1=mad(a1·row1+a2·row2),b2=mad(a1·row1-a2·row2),其中row1与row2代表多小波高频系数的2个分量;
式中,μ=2logn1,n1为多小波第j层分解的数据长度;
该估计算子称为SureBlock算子,当Td≤γd时,该算子退化为分块长度为1的James-Stein估计算子。
上述方案中,步骤(6)的具体实施方式为:
计算降噪结果中时域冲击波形的平均时间间隔作为冲击波形的周期,根据输入轴转速和齿轮箱传动比计算出各转轴的旋转周期,通过比较冲击波形的周期和齿轮箱各轴的旋转周期,判断齿轮箱是否存在损伤,并准确诊断出损伤齿轮所在的转轴。
由于本发明在齿轮箱故障诊断中采用多小波自适应分块阈值降噪时域诊断技术,有效提取出连轧机传动齿轮箱损伤特征。本发明具有以下区别于传统方法的显著优势:
1)根据多小波系数之间的相关性,以Stein(Stein C.Estimation of themean of a multivariate normal distribution.The Annals of Statistics,1981,9(6):1135-1151)无偏风险估计误差最小作为约束条件,在多小波分解的每一层自适应地选取最优的邻域分块长度和阈值,能够有效消除噪声干扰,准确提取信号特征;
2)利用所构造的自适应邻域分块阈值降噪方法分析连轧机传动齿轮箱振动信号,获得降噪后的重构信号,通过分析时域脉冲波形的周期性,有效提取出连轧机传动齿轮箱损伤特征。
3)本发明结果可靠,实时性好,简单易行,通用性强,适用于连轧机等重载设备传动齿轮箱的损伤诊断。
附图说明
下面结合附图及具体实施方式对本发明作进一步的详细说明。
图1为连轧机传动齿轮箱结构及测点布置示意图。
图2为连轧机传动齿轮箱振动信号时域波形(a)及频谱(b)。
图3为本发明多小波自适应分块阈值降噪算法流程图。
图4为连轧机传动齿轮箱信号多小波自适应分块阈值降噪结果。
图5为连轧机传动齿轮箱损伤照片。
具体实施方式
一种齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法,按以下步骤实施:
(1)对含噪信号进行预处理,得到矢量输入信号;
(2)对矢量输入信号进行多小波分解,得到高频系数和低频系数;
(3)根据判断条件对高频系数进行阈值处理:符合条件时采用自适应邻域分块阈值算法,否则采用逐点比较算法,得到降噪后的高频系数;
(4)对低频系数和经过阈值处理后的高频系数进行多小波反变换,重构降噪后的矢量输出信号;
(5)对矢量输出信号进行后处理,得到一维降噪结果;
(6)通过分析降噪结果中的时域脉冲波形的周期性,提取出齿轮箱的损伤特征进行故障诊断。
以下予以分述:
带钢热轧机组中的精轧机由双电机驱动,其传动齿轮箱为单级减速器,传动比是22/65,连轧机传动齿轮箱结构及测点布置如图1所示。监测系统采用速度传感器对轧机1~7#测点的轴承进行实时监测。高速小齿轮和低速大齿轮的转速分别是4.55Hz和1.54Hz。采集6#测点轴承的振动信号,采样频率为5120Hz,数据长度为4096。
传动齿轮箱振动信号时域波形及频谱如图2(a)和2(b)所示,图2(a)中横坐标表示时间,单位为s;纵坐标表示振动幅值,单位为mm/s。从振动信号时域波形中可以隐约看到周期性冲击信号,其时间间隔为0.22s,和高速轴转频一致,局部冲击信号被设备运行时产生的噪声和背景噪声淹没。图2(b)中横坐标表示频率,单位为Hz;纵坐标表示振动幅值,单位为mm/s。频谱中可以看到100Hz、497.5Hz及601.3Hz的谱峰,对应于减速箱的齿轮啮合频率100Hz及倍频成分。
根据判断条件对高频系数进行阈值处理:符合条件时采用自适应邻域分块阈值算法,否则采用逐点比较算法,如图3所示。设含噪信号S是所关心的故障特征F与噪声E的叠加,可以表示为:
S=F+E
以2×2维多小波变换为例,S=[s(1),s(2)]T表示经过预处理后的2维振动信号,噪声E服从多元正态分布N(0,Vj),在第k个邻域分块中,定义变量作为阈值降噪的基本变量,将其中向量化的系数转化为标量,矩阵Vj为多小波分解第j层的高频系数的协方差矩阵;
定义:
mad(y)=1.4826·median(|y-median(y)|)
式中,函数median()用于寻找数列y中的中位数,矩阵Vj可表达为:
式中,a1=1/mad(row1),a2=1/mad(row2),b1=mad(a1·row1+a2·row2),b2=mad(a1·row1-a2·row2),其中row1与row2代表多小波高频系数的2个分量;
式中,μ=2logn1,n1为多小波第j层分解的数据长度;
该估计算子称为SureBlock算子,当Td≤γd时,该算子退化为分块长度为1的James-Stein估计算子(Cai T,Zhou H.A data-driven block thresholdingapproach to wavelet estimation.The Annals of Statistics,2009,37:569-595)。
使用所构造的James-Stein估计算子处理高频系数,得到降噪后的高频系数;其次对低频系数和经过阈值处理后的高频系数进行多小波反变换,重构降噪后的矢量输出信号;对矢量输出信号进行后处理,得到一维降噪结果。连轧机传动齿轮箱信号降噪结果如图4所示。
当齿轮箱中齿轮出现损伤时,如图5所示,齿轮每转一周,损伤轮齿的啮合在振动信号中产生一组冲击波形,冲击波形的周期和齿轮的旋转周期一致。计算降噪结果中时域冲击波形的平均时间间隔作为冲击波形的周期,根据输入轴转速和齿轮箱传动比计算出各转轴的旋转周期,通过比较冲击波形的周期和齿轮箱各轴的旋转周期,判断齿轮箱是否存在损伤,并准确诊断出损伤齿轮所在的转轴。
为了验证本发明所述方法的正确性,在带钢热轧机组的精轧机上进行实验。如图1所示,精轧机由双电机驱动,其传动齿轮箱为单级减速器,传动比是22/65。实验中,传动齿轮箱的高速轴小齿轮存在两处由于高温熔焊导致的局部胶合破坏故障,且两处故障大概相距1/3个圆周,即相隔约8~10个齿。图5(a)中为局部胶合故障,图5(b)中为在整个齿宽范围内的严重胶合故障。监测系统采用速度传感器对轧机测点1~7#的轴承进行实时监测。高速小齿轮和低速大齿轮的转速分别是4.55Hz和1.54Hz。采集测点6轴承的振动信号,采样频率为5120Hz,数据长度为4096。如图2(a)所示,振动信号时域波形中可以隐约看到周期性冲击信号,其时间间隔为0.22s,和高速轴转频一致,局部冲击信号被设备运行时产生的噪声和背景噪声淹没。
多小波自适应分块阈值降噪算法的流程如图3所示,可以归纳为以下步骤:1)对原始信号进行预处理,得到矢量输入信号;2)对矢量输入信号进行多小波分解;3)根据判断条件对高频系数采用自适应邻域分块阈值方法进行降噪处理,得到降噪后的样本;4)对低频系数和阈值处理后的高频系数进行多小波反变换,重构降噪后的矢量输出信号;5)后处理得到一维输出信号。选择GHM多小波分解预处理后的多维振动信号,分解层数选为5层。根据Stein无偏风险估计误差最小准则选择最优的邻域分块长度和阈值,对每一层细节系数采取邻域分块阈值降噪方法进行降噪处理,获得降噪后的重构信号,通过分析时域脉冲波形的周期性,有效提取出连轧机传动齿轮箱的损伤特征。如图4所示,降噪结果中交替出现了两种强弱不等的周期冲击波形I1和I2,其周期均对应于高速轴旋转周期0.22s,反映齿轮箱高速轴小齿轮存在两处损伤程度不同的局部故障。其中,大冲击I2反映小齿轮存在一处严重损伤;小冲击I1反映小齿轮存在一处较严重损伤。交替出现的周期冲击波形I1和I2的时间间隔约为小齿轮旋转周期的1/3,说明两处局部故障相距约为小齿轮圆周的1/3。
Claims (2)
1.一种齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法,其特征在于:
(1)对含噪信号进行预处理,得到矢量输入信号;
(2)对矢量输入信号进行多小波分解,得到高频系数和低频系数;
(3)根据判断条件对高频系数进行阈值处理:符合条件时采用自适应邻域分块阈值算法,否则采用逐点比较算法,得到降噪后的高频系数;
(4)对低频系数和经过阈值处理后的高频系数进行多小波反变换,重构降噪后的矢量输出信号;
(5)对矢量输出信号进行后处理,得到一维降噪结果;
(6)通过分析降噪结果中的时域脉冲波形的周期性,提取出齿轮箱的损伤特征进行故障诊断;
其中,步骤(3)所述的自适应邻域分块阈值算法包括下述步骤:
1)设含噪信号S是所关心的故障特征F与噪声E的叠加,可以表示为:S=F+E
以2×2维多小波变换为例,S=[s(1),s(2)]T表示经过预处理后的2维振动信号,噪声E服从多元正态分布N(0,Vj),在第k个邻域分块中,定义变量作为阈值降噪的基本变量,将其中向量化的系数转化为标量,矩阵Vj为多小波分解第j层的高频系数的协方差矩阵;
定义:
mad(y)=1.4826·median(|y-median(y)|)
式中,函数median()用于寻找数列y中的中位数,矩阵Vj可表达为:
式中,a1=1/mad(row1),a2=1/mad(row2),b1=mad(a1·row1+a2·row2),b2=mad(a1·row1-a2·row2),其中row1与row2代表多小波高频系数的2个分量;
在自适应邻域分块阈值降噪算法中,第k个邻域分块中多小波系数组合为其中L为邻域分块长度;阈值表示为:
式中,μ=2logn1,n1为多小波第j层分解的数据长度;
该估计算子称为SureBlock算子,当Td≤γd时,该算子退化为分块长度为1的James-Stein估计算子。
2.如权利要求1所述的一种齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法,其特征在于,步骤(6)的具体实施方式为:计算降噪结果中时域冲击波形的平均时间间隔来作为冲击波形的周期,根据输入轴转速和齿轮箱传动比计算出各转轴的旋转周期,通过比较冲击波形的周期和齿轮箱各转轴的旋转周期,判断齿轮箱是否存在损伤,确定出损伤齿轮所在的转轴。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210361690.4A CN102879195B (zh) | 2012-09-25 | 2012-09-25 | 齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210361690.4A CN102879195B (zh) | 2012-09-25 | 2012-09-25 | 齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102879195A CN102879195A (zh) | 2013-01-16 |
CN102879195B true CN102879195B (zh) | 2014-03-12 |
Family
ID=47480606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210361690.4A Expired - Fee Related CN102879195B (zh) | 2012-09-25 | 2012-09-25 | 齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102879195B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103745085B (zh) * | 2013-12-16 | 2017-01-11 | 西安交通大学 | 一种旋转机械振动信号的数据驱动阈值降噪方法 |
CN103837345B (zh) * | 2014-03-25 | 2016-06-22 | 上海电机学院 | 齿轮箱故障诊断方法及装置 |
CN105258940B (zh) * | 2015-11-30 | 2018-02-09 | 上海无线电设备研究所 | 机械故障定量提取的标准化多小波与多小波包变换方法 |
CN106596111A (zh) * | 2016-12-10 | 2017-04-26 | 国网福建省电力有限公司 | 基于多小波的水电机组故障特征自适应提取方法 |
CN109708872A (zh) * | 2017-10-20 | 2019-05-03 | 株洲中车时代电气股份有限公司 | 一种列车齿轮箱联轴节故障诊断方法、装置及系统 |
CN108195932B (zh) * | 2017-11-30 | 2020-10-30 | 北京宇航系统工程研究所 | 一种飞行器管路损伤超声导波定量评估方法 |
CN112098526B (zh) * | 2020-07-23 | 2021-11-19 | 西安交通大学 | 基于激光超声技术的增材制件近表面缺陷特征提取方法 |
CN113804388B (zh) * | 2021-09-15 | 2024-04-02 | 西安因联信息科技有限公司 | 基于时域分析的机械设备旋转冲击故障检测方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008292288A (ja) * | 2007-05-24 | 2008-12-04 | Mitsubishi Electric Engineering Co Ltd | 減速機の軸受診断装置 |
CN102507186A (zh) * | 2011-11-01 | 2012-06-20 | 西安交通大学 | 一种基于特色参数的行星齿轮箱状态监测与故障辨识方法 |
CN102620928A (zh) * | 2012-03-02 | 2012-08-01 | 燕山大学 | 基于小波半软阈值和emd的风电齿轮箱故障诊断方法 |
CN102636347A (zh) * | 2012-04-24 | 2012-08-15 | 西安交通大学 | 一种用于变转速齿轮箱的振动信号时域同步平均方法 |
-
2012
- 2012-09-25 CN CN201210361690.4A patent/CN102879195B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008292288A (ja) * | 2007-05-24 | 2008-12-04 | Mitsubishi Electric Engineering Co Ltd | 減速機の軸受診断装置 |
CN102507186A (zh) * | 2011-11-01 | 2012-06-20 | 西安交通大学 | 一种基于特色参数的行星齿轮箱状态监测与故障辨识方法 |
CN102620928A (zh) * | 2012-03-02 | 2012-08-01 | 燕山大学 | 基于小波半软阈值和emd的风电齿轮箱故障诊断方法 |
CN102636347A (zh) * | 2012-04-24 | 2012-08-15 | 西安交通大学 | 一种用于变转速齿轮箱的振动信号时域同步平均方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102879195A (zh) | 2013-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102879195B (zh) | 齿轮箱损伤的多小波自适应分块阈值降噪时域诊断方法 | |
Yin et al. | Weak fault feature extraction of rolling bearings based on improved ensemble noise-reconstructed EMD and adaptive threshold denoising | |
Ding et al. | Multi-objective iterative optimization algorithm based optimal wavelet filter selection for multi-fault diagnosis of rolling element bearings | |
CN103645052B (zh) | 一种风电机组齿轮箱远程在线状态监测及寿命评估方法 | |
CN105241666B (zh) | 一种基于信号稀疏表示理论的滚动轴承故障特征提取方法 | |
Liang et al. | Maximum average kurtosis deconvolution and its application for the impulsive fault feature enhancement of rotating machinery | |
Cheng et al. | Gear fault identification based on Hilbert–Huang transform and SOM neural network | |
CN111089726B (zh) | 一种基于最优维数奇异谱分解的滚动轴承故障诊断方法 | |
CN103499445B (zh) | 一种基于时频切片分析的滚动轴承故障诊断方法 | |
CN103439110B (zh) | 滚动轴承早期微弱故障诊断方法 | |
Zhang et al. | Multi-fault diagnosis of gearbox based on resonance-based signal sparse decomposition and comb filter | |
Strömbergsson et al. | Bearing monitoring in the wind turbine drivetrain: A comparative study of the FFT and wavelet transforms | |
CN102937522B (zh) | 一种齿轮箱复合故障诊断方法及系统 | |
Hao et al. | Morphological undecimated wavelet decomposition for fault diagnostics of rolling element bearings | |
CN108388908B (zh) | 基于k-svd和滑窗降噪的滚动轴承冲击故障诊断方法 | |
Chen et al. | An improved singular value decomposition-based method for gear tooth crack detection and severity assessment | |
CN103983452B (zh) | 利用混合域特征向量和灰色关联分析对行星齿轮箱进行故障模式识别的方法 | |
CN107817106A (zh) | 基于贝叶斯残余变换‑奇异值分解‑高斯混合隐马尔科夫模型框架的滚动轴承故障诊断方法 | |
CN105698922A (zh) | 基于改进emd和谱峭度法变压器振动故障特征提取方法 | |
Zhang et al. | Adaptive filtering enhanced windowed correlated kurtosis for multiple faults diagnosis of locomotive bearings | |
CN107886078A (zh) | 一种基于分层自适应阈值函数的小波阈值降噪方法 | |
CN104614769A (zh) | 一种压制地震面波的聚束滤波方法 | |
Zhang et al. | Fault diagnosis for gearbox based on EMD-MOMEDA | |
Yan et al. | A bearing fault feature extraction method based on optimized singular spectrum decomposition and linear predictor | |
Yao et al. | Fault detection of complex planetary gearbox using acoustic signals |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140312 Termination date: 20160925 |