CN105510032A - 基于谐躁比指导的解卷积方法 - Google Patents
基于谐躁比指导的解卷积方法 Download PDFInfo
- Publication number
- CN105510032A CN105510032A CN201510926945.0A CN201510926945A CN105510032A CN 105510032 A CN105510032 A CN 105510032A CN 201510926945 A CN201510926945 A CN 201510926945A CN 105510032 A CN105510032 A CN 105510032A
- Authority
- CN
- China
- Prior art keywords
- sigma
- filter
- iteration
- signal
- deconvolution
- 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
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/04—Bearings
- G01M13/045—Acoustic or vibration analysis
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
基于谐躁比指导的解卷积方法,首先对采集的振动信号进行截断和去均值处理;然后对截断后的信号直接进行解卷积处理,兼顾故障的冲击性和周期性特征,将峭度和谐躁比作为目标函数进行联立求迭代表达式,进而极大地减小了随机冲击的干扰;最后对解卷积后的信号进行包络分析,从包络谱中能提取故障特征频率,本发明在特征频率的提取过程中不需要人为参与,有利于实现故障特征提取和诊断监测的自动化,节约时间,效率更高。
Description
技术领域
本发明涉及机械设备故障诊断技术领域,特别涉及基于谐躁比指导的解卷积(Harmonics-to-NoiseRatioGuidedDeconvolution,HNRGD)方法。
背景技术
振动分析为现阶段机械设备故障诊断最为有效途径之一,机械设备的状态劣化往往表现为振动信息的变化或异常。峭度是描述信号冲击性的有效指标,当振动信号表现为明显的冲击性时其峭度值较大,被广泛应用于机械***监测诊断中。2007年最小熵解卷积(Minimumentropydeconvolution,MED)方法被Sawalhi和Randall等学者首次应用于滚动轴承的故障诊断领域,并取得了一定的效果。MED是一种不需要任何先验假设的信号时域盲解卷积技术,通过迭代选择一个(finiteimpulseresponse)FIR滤波器来最小化滤波信号的输出熵(即峭度最大化),旨在最小化噪声的同时提取出故障冲击,因此其在高噪声下也能得到理想的诊断结果。虽然MED对冲击性的增强与提取具有明显的效果,但是其目标函数仅仅是寻求滤波信号峭度值最大化,所以易受随机孤立冲击的干扰。即当故障信号中既存在周期性冲击序列和随机孤立冲击时,通过MED增强技术有可能只能增强孤立冲击,而对真正反映故障特征的周期性冲击序列没有效果。鉴于上述问题,2012年McDonald等学者提出了最大相关峭度解卷积(Maximumcorrelatedkurtosisdeconvolution,MCKD)算法,提出了相关峭度的概念,兼顾了故障的冲击性和周期性特征,减小了随机冲击的干扰,并成功地将其应用到齿根裂纹故障诊断上。然而MCKD方法需要提前准确估计故障特征频率,而在工程实际中,测速计往往存在误差,而且复杂的工况使得设备转速不可能保持恒定。并且当转速精准确定的前提下,MCKD也很难提取带有随机滑动的轴承故障冲击信号。这些不足给MCKD的应用推广带来了诸多不便。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供基于谐躁比指导的解卷积方法,实现机械设备不用预估故障特征频率等预处理,就能实现精确的故障诊断。
为实现上述目的,本发明方案采取的技术方案为:
基于谐躁比指导的解卷积方法,包括以下步骤:
步骤一:将振动加速度传感器吸附于对被测试滚动轴承的轴承座上,并对振动信号进行高频采样、截断和去均值处理,将振动信号记为x;接下来对振动信号x进行HNRGD处理;
步骤二:将峭度作为目标函数对滤波系数进行求偏导:
其中N为输出信号y的长度,使用其中f(l)为滤波器系数,l=1,2,……,L,L为滤波器长度,求导后的结果为:
再以谐躁比(Harmonics-to-NoiseRatio,HNR)为目标函数对滤波系数进行求偏导:
其中t为时间,T为周期,对目标函数进行离散处理,并使用 求导后的结果为:
联立(2)和(4),使用得到解卷积迭代表达式:
将上式写成矩阵的形式:
Af=b(6)
其中:
b--逆滤波器的输入振动信号x、输出y的互相关,b为L维列向量;
A--输入振动信号x的自相关,A为L×L维矩阵;
f--逆滤波器的滤波器系数,f为L维列向量;
首先计算自相关矩阵A;再假设逆滤波器的初始值f(0),设置滤波器长度L=100,并给定初始滤波器系数为[00……1-1……00],用y(0)和x(0)计算列向量b(1);然后求解新的滤波器系数f(1)=A-1b(1);通过对滤波器系数的计算进行收敛判断;如果通过比较后一次迭代与前一次迭代的谐躁比与峭度值,如果这两个参数都在增加那么迭代结束,如果并不是都在增加就更新滤波器系数进行下次迭代,直到满足收敛条件或者达到设定的最大迭代次数为30次,得到经过HNRGD处理后的信号yk;
步骤三:对HNRGD处理后的信号yk进行包络分析并得到包络谱,对包络谱进行分析,进而提取故障特征频率。
本发明相比于现有技术,具有以下有益效果:
a)本发明在传统MED方法的基础上,融入了谐躁比的概念,能兼顾故障的冲击性和周期性特征,极大地减小了随机冲击的干扰。
b)本发明不需要任何先验知识,也无需对***故障特征频率进行预估,方法具有鲁棒性。
c)本发明在特征频率的提取过程中不需要人为参与,有利于实现故障特征提取和诊断监测的自动化,节约时间,效率更高。
附图说明
图1为实施例美国西储大学轴承实验台示意图。
图2为本发明方法流程图。
图3为实例中对原信号进行截断处理后的振动信号x。
图4为实例中对截断后的信号进行HNRGD方法处理后的信号yk。
图5为实施例中经HNRGD方法处理后信号的包络谱。
图6为实施例中利用MCKD方法对截断信号进行改变采样率重采样之后的信号。
图7为实施例中从采样信号经MCKD滤波后的信号。
图8为实施例中经MCKD滤波后的信号的包络谱。
具体实施方式
下面结合附图和实施例对本发明做详细描述。
以下以美国西储大学轴承实验台为例进行说明,该试验台由电机、联轴器、转矩传感器、功率计以及轴承等组成,如图1所示,该实验台检测对象为驱动电机的驱动端轴承和风扇端轴承,轴承型号分别为:6205-2RSJEMSKF和6203-2RSJEMSKF,采样频率为12kHz。整个装置由电机驱动,将扭矩从电机沿转矩传感器传递到输出端的功率计。
本发明采用的数据是由美国西储大学提供的上述试验台测得的数据中编号为OR0146_3.mat的一组数据。根据W.A.SmithandR.B.Randall近期发表文章《RollingelementbearingdiagnosticsusingtheCaseWesternReserveUniversitydata:Abenchmarkstudy》描述,该组数据通过传统的方法无法进行故障信息的提取,因此该数据具有说服力。
利用基于谐躁比指导的解卷积方法对该试验台中驱动电机驱动端轴承进行诊断,应用本发明对实验数据进行分析并和传统MCKD方法进行对比。
如图2所示,基于谐躁比指导的解卷积方法,包括以下步骤:
步骤一:将振动加速度传感器吸附于对被测试滚动轴承的轴承座上,并对振动信号进行高频采样、截断和去均值处理,其中采样频率为12kHz,在使用数据时需要去掉起始噪声部分,截取整段信号中0.5-2.0s共1.5s的数据,如图所示3,将该1.5s振动信号记为x;接下来对振动信号x进行HNRGD处理;
步骤二:将峭度作为目标函数对滤波系数进行求偏导:
其中N为输出信号y的长度,使用其中f(l)为滤波器系数,l=1,2,……,L,L为滤波器长度,求导后的结果为:
再以谐躁比(Harmonics-to-NoiseRatio,HNR)为目标函数对滤波系数进行求偏导:
其中t为时间,T为周期,对目标函数进行离散处理,使用 求导后的结果为:
联立(2)和(4)使用得到解卷积迭代表达式:
将上式写成矩阵的形式:
Af=b(6)
其中:
b--逆滤波器的输入振动信号x、输出y的互相关,b为L维列向量;
A--输入振动信号x的自相关,A为L×L维矩阵;
f--逆滤波器的滤波器系数,f为L维列向量;
首先计算自相关矩阵A;再假设逆滤波器的初始值f(0),为与MCKD具有可比性,将滤波器长度设置与MCKD方法一样为100,并给定初始滤波器系数为[00……1-1……00],用y(0)和x(0)计算列向量b(1);然后求解新的滤波器系数f(1)=A-1b(1);通过对滤波器系数的计算进行收敛判断;如果通过比较后一次迭代与前一次迭代的谐躁比与峭度值,如果这两个参数都在增加那么迭代结束,如果并不是都在增加就更新滤波器系数进行下次迭代,直到满足收敛条件或者达到设定的最大迭代次数为30次;得到经过HNRGD处理后的信号yk,如图4所示;
步骤三:对HNRGD处理后的信号yk进行包络分析并得到包络谱,如图5所示,BPFO(ballpassfrequencyofouterrace)表示外圈故障特征频率,对包络谱进行分析,可以清楚的看到轴承外圈特征故障频率102.9435Hz以及其2倍频、4倍频、6倍频和7倍频,效果十分明显。
如图6所示为利用MCKD方法对振动信号x进行改变采样率重采样之后的效果,按照McDonald等提出的MCKD方法对重采样后的信号进行处理,对该方法的参数按照McDonald提出的推荐设置,其中参数设置为精确计算故障特征频率f=102.9435Hz,移位距离T=fs*1/f=12000/102.9435,移位次数M=3,滤波器长度设置为100,最大滤波次数设置为30,图7是经MCKD滤波后的信号,图8为其包络谱,从图8可以看出,在精确知道电机转速,轴承参数等一系列参数之后,该例中还必须知道轴承***外圈故障以及精确计算出其故障特征频率的众多前提下,才能艰难地从包络谱中看到故障特征频率的6倍频(617.6610Hz)附近的616Hz的频率成分。然而,在工程实际中,***经常处于变转速变载荷的复杂工况下,而且测速计往往存在误差,很难精确预估故障特征频率。因此相比于此,本发明提出的方法能避免以上提到的困难并能精确的提取故障特征频率,所以本发明优势明显。
Claims (1)
1.基于谐躁比指导的解卷积方法,其特征在于,包括以下步骤:
步骤一:将振动加速度传感器吸附于对被测试滚动轴承的轴承座上,并对振动信号进行高频采样、截断和去均值处理,将振动信号记为x;接下来对振动信号x进行HNRGD处理;
步骤二:将峭度作为目标函数对滤波系数进行求偏导:
其中N为输出信号y的长度,使用其中f(l)为滤波器系数,l=1,2,……,L,L为滤波器长度,求导后的结果为:
再以谐躁比(Harmonics-to-NoiseRatio,HNR)为目标函数对滤波系数进行求偏导:
其中t为时间,T为周期,对目标函数进行离散处理,使用 求导后的结果为:
联立(2)和(4)并使用得到解卷积迭代表达式:
将上式写成矩阵的形式:
Af=b(6)
其中:
b--逆滤波器的输入振动信号x、输出y的互相关,b为L维列向量;
A--输入振动信号x的自相关,A为L×L维矩阵;
f--逆滤波器的滤波器系数,f为L维列向量;
首先计算自相关矩阵A;再假设逆滤波器的初始值f(0),设置滤波器长度L=100,并给定初始滤波器系数为[00……1-1……00],用y(0)和x(0)计算列向量b(1);然后求解新的滤波器系数f(1)=A-1b(1);通过对滤波器系数的计算进行收敛判断;如果通过比较后一次迭代与前一次迭代的谐躁比与峭度值,如果这两个参数都在增加那么迭代结束,如果并不是都在增加就更新滤波器系数进行下次迭代,直到满足收敛条件或者达到设定的最大迭代次数为30次,得到经过HNRGD处理后的信号yk;
步骤三:对HNRGD处理后的信号yk进行包络分析并得到包络谱,对包络谱进行分析,进而提取故障特征频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510926945.0A CN105510032B (zh) | 2015-12-11 | 2015-12-11 | 基于谐噪比指导的解卷积方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510926945.0A CN105510032B (zh) | 2015-12-11 | 2015-12-11 | 基于谐噪比指导的解卷积方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105510032A true CN105510032A (zh) | 2016-04-20 |
CN105510032B CN105510032B (zh) | 2017-12-26 |
Family
ID=55718174
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510926945.0A Active CN105510032B (zh) | 2015-12-11 | 2015-12-11 | 基于谐噪比指导的解卷积方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105510032B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106525223A (zh) * | 2016-11-01 | 2017-03-22 | 苏州微著设备诊断技术有限公司 | 一种齿轮传动装置异响的下线检测方法 |
CN106564012A (zh) * | 2016-11-01 | 2017-04-19 | 苏州微著设备诊断技术有限公司 | 一种磨削加工颤振的检测方法 |
CN107560718A (zh) * | 2017-07-25 | 2018-01-09 | 苏州微著设备诊断技术有限公司 | 一种基于自相关峭度最大化的解卷积方法 |
CN109682601A (zh) * | 2019-03-04 | 2019-04-26 | 北京天泽智云科技有限公司 | 一种变转速工况下滚动轴承的早期故障识别方法 |
CN111024398A (zh) * | 2019-12-25 | 2020-04-17 | 西安交通大学 | 一种无需周期的最大相关峭度解卷积方法 |
CN112099226A (zh) * | 2020-03-06 | 2020-12-18 | 中国工程物理研究院激光聚变研究中心 | 一种用于丝靶瞄准的激光束引导方法 |
CN112326169A (zh) * | 2020-10-28 | 2021-02-05 | 未来振动(北京)测试技术有限公司 | 振动式异音异响检测*** |
CN113962250A (zh) * | 2021-09-14 | 2022-01-21 | 西南交通大学 | 自适应最小相关广义反卷积方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070027681A1 (en) * | 2005-08-01 | 2007-02-01 | Samsung Electronics Co., Ltd. | Method and apparatus for extracting voiced/unvoiced classification information using harmonic component of voice signal |
CN101667425A (zh) * | 2009-09-22 | 2010-03-10 | 山东大学 | 一种对卷积混叠语音信号进行盲源分离的方法 |
CN104251934A (zh) * | 2013-06-26 | 2014-12-31 | 华为技术有限公司 | 谐波分析方法和装置以及确定谐波间杂波的方法和装置 |
CN104330257A (zh) * | 2014-11-26 | 2015-02-04 | 天津工业大学 | 一种行星齿轮传动***故障诊断方法 |
CN104819766A (zh) * | 2015-05-13 | 2015-08-05 | 西安交通大学 | 基于谐噪比的包络解调频带确定方法 |
-
2015
- 2015-12-11 CN CN201510926945.0A patent/CN105510032B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070027681A1 (en) * | 2005-08-01 | 2007-02-01 | Samsung Electronics Co., Ltd. | Method and apparatus for extracting voiced/unvoiced classification information using harmonic component of voice signal |
CN101667425A (zh) * | 2009-09-22 | 2010-03-10 | 山东大学 | 一种对卷积混叠语音信号进行盲源分离的方法 |
CN104251934A (zh) * | 2013-06-26 | 2014-12-31 | 华为技术有限公司 | 谐波分析方法和装置以及确定谐波间杂波的方法和装置 |
CN104330257A (zh) * | 2014-11-26 | 2015-02-04 | 天津工业大学 | 一种行星齿轮传动***故障诊断方法 |
CN104819766A (zh) * | 2015-05-13 | 2015-08-05 | 西安交通大学 | 基于谐噪比的包络解调频带确定方法 |
Non-Patent Citations (2)
Title |
---|
YHONGHAO MIAO: "Sparse maximum harmonics-to-noise-ratio deconvolution for weak fault signature detection in bearings", 《MEAS. SCI. TECHNOL》 * |
雷亚国等: "行星齿轮箱故障诊断技术的研究进展", 《机械工程学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106525223A (zh) * | 2016-11-01 | 2017-03-22 | 苏州微著设备诊断技术有限公司 | 一种齿轮传动装置异响的下线检测方法 |
CN106564012A (zh) * | 2016-11-01 | 2017-04-19 | 苏州微著设备诊断技术有限公司 | 一种磨削加工颤振的检测方法 |
CN106564012B (zh) * | 2016-11-01 | 2018-08-21 | 苏州微著设备诊断技术有限公司 | 一种磨削加工颤振的检测方法 |
CN107560718A (zh) * | 2017-07-25 | 2018-01-09 | 苏州微著设备诊断技术有限公司 | 一种基于自相关峭度最大化的解卷积方法 |
CN109682601A (zh) * | 2019-03-04 | 2019-04-26 | 北京天泽智云科技有限公司 | 一种变转速工况下滚动轴承的早期故障识别方法 |
CN111024398A (zh) * | 2019-12-25 | 2020-04-17 | 西安交通大学 | 一种无需周期的最大相关峭度解卷积方法 |
CN111024398B (zh) * | 2019-12-25 | 2021-01-15 | 西安交通大学 | 一种无需周期的最大相关峭度解卷积方法 |
CN112099226A (zh) * | 2020-03-06 | 2020-12-18 | 中国工程物理研究院激光聚变研究中心 | 一种用于丝靶瞄准的激光束引导方法 |
CN112099226B (zh) * | 2020-03-06 | 2022-02-08 | 中国工程物理研究院激光聚变研究中心 | 一种用于丝靶瞄准的激光束引导方法 |
CN112326169A (zh) * | 2020-10-28 | 2021-02-05 | 未来振动(北京)测试技术有限公司 | 振动式异音异响检测*** |
CN113962250A (zh) * | 2021-09-14 | 2022-01-21 | 西南交通大学 | 自适应最小相关广义反卷积方法 |
CN113962250B (zh) * | 2021-09-14 | 2023-05-12 | 西南交通大学 | 自适应最小相关广义反卷积方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105510032B (zh) | 2017-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105510032A (zh) | 基于谐躁比指导的解卷积方法 | |
Wang et al. | Feature extraction of rolling bearing’s early weak fault based on EEMD and tunable Q-factor wavelet transform | |
CN104215456B (zh) | 一种基于平面聚类和频域压缩感知重构的机械故障诊断方法 | |
CN105424359A (zh) | 一种基于稀疏分解的齿轮和轴承混合故障特征提取方法 | |
CN105699082B (zh) | 一种稀疏化的最大谐噪比解卷积方法 | |
CN103424258A (zh) | 一种滚动轴承的故障诊断方法 | |
CN103674234B (zh) | 一种风力发电机组振动异常的状态预警方法及*** | |
CN110243605B (zh) | 多源时频脊线提取方法 | |
CN104006962A (zh) | 一种齿轮故障特征提取方法及*** | |
CN109765052B (zh) | 基于goa-asr的行星齿轮箱早期故障诊断方法 | |
Cong et al. | Research of singular value decomposition based on slip matrix for rolling bearing fault diagnosis | |
CN107506710A (zh) | 一种滚动轴承复合故障提取方法 | |
CN102998118B (zh) | 一种基于形态学滤波和复杂度测度的轴承定量诊断方法 | |
Sheng et al. | Applications in bearing fault diagnosis of an improved Kurtogram algorithm based on flexible frequency slice wavelet transform filter bank | |
CN105092241A (zh) | 一种齿轮局部故障诊断方法及*** | |
CN111896260B (zh) | NGAs同步优化小波滤波器与MCKD的轴承故障诊断方法 | |
Zhang et al. | Improved local cepstrum and its applications for gearbox and rolling bearing fault detection | |
CN105043767A (zh) | 一种滚动球轴承外圈剥落故障双冲击特征提取方法及*** | |
CN109063668A (zh) | 一种基于峰值保留降采样的冲击信号包络解调方法 | |
CN103267652B (zh) | 一种智能早期设备故障在线诊断方法 | |
CN107121283B (zh) | 一种基于离散随机分离算法的齿轮状态监测指标提取方法 | |
Zhang et al. | Fault diagnosis for gearbox based on EMD-MOMEDA | |
CN102305661A (zh) | 一种斜拉桥拉索振动信号的降噪处理方法 | |
Yan et al. | A bearing fault feature extraction method based on optimized singular spectrum decomposition and linear predictor | |
CN105987809A (zh) | 基于随机共振的离心式压缩机半开式叶轮裂纹检测方法 |
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 |