CN101251445A - 旋转机械碰摩声发射信号的分形特征分析方法 - Google Patents
旋转机械碰摩声发射信号的分形特征分析方法 Download PDFInfo
- Publication number
- CN101251445A CN101251445A CNA2008100238086A CN200810023808A CN101251445A CN 101251445 A CN101251445 A CN 101251445A CN A2008100238086 A CNA2008100238086 A CN A2008100238086A CN 200810023808 A CN200810023808 A CN 200810023808A CN 101251445 A CN101251445 A CN 101251445A
- Authority
- CN
- China
- Prior art keywords
- dimension
- delta
- signal
- frame
- fractal
- 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
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明涉及旋转机械碰摩声发射信号的分形特征分析方法。本发明技术方案是从摩擦点的AE传感器提取AE信号,AE信号上叠加高斯白噪声和非平稳噪声,采样、量化、分帧,将每一帧含有噪声的数据等分成长为m的短时帧数据段;用基于波形的分形维数算法计算各短时帧数据段的分维值,取正中间序号的数据作为中值滤波的输出。本发明解决了现有方法存在的计算量非常大、Katz维存在着精度不高、短时帧分维曲线波动较大的缺点。本发明在强噪声环境下,对数波长维比Katz维和盒维具有更强区分噪声和AE信号的能力,并且AE信号帧的对数波长维比Katz维和盒维更接近于纯净AE信号的分维,从而能够有效地检测并反映出碰摩AE事件的发生,且精确度高,计算量减少。
Description
技术领域
本发明涉及旋转机械碰摩声发射信号的分析方法,特别涉及旋转机械碰摩声发射信号的分形特征分析方法。
背景技术
声发射技术AE(Acoustic Emission)以其灵敏度高、频响范围宽、信息量大、动态检测等特点,在碰摩故障检测中显示了其优越性。碰摩发生时伴有较强的AE信号产生。由AE特征参数及其变化可以判断碰摩的产生及强度,确定碰摩点的位置。但是由于AE源的多样性和噪声干扰的复杂性,AE的特征参数常常无法反映设备的真实状态,尤其在强噪声环境下,AE的特征识别与提取一直是个难题。利用分形理论对碰摩故障进行了理论计算与实验研究,将其作为定量描述AE信号的特征参数。可以有效的识别设备的状态和预警,取得良好的效果。但在实际应用中,分形维算法的计算量是限制其应用的瓶颈。
在本发明之前,目前常用的分形维方法有盒维,关联维和Katz维方法等。对于盒维来说,由于随着***维数的增加,需要的盒子数量指数增加,需要用于计算的数据长度也相应增加,使得高维***的收敛非常困难。而关联维严重依赖于两个参数:最小延迟时间和最小嵌入维。如果最小延迟时间太大,达到吸引子的大小,则所有矢量都是相关的,太小则计算出的相关维就是噪声的相关维。对于AE信号而言,由于每帧的最小延迟时间各不相同,因此每一帧都必须计算其最小延迟时间才能计算出准确的关联维,计算量非常大。Katz维不需要计算极限,对参数依赖性也不强,而且在抗噪声能力上优于其他算法,但Katz维存在着精度不高的缺点。
发明内容
本发明的目的就在于解决上述现有技术的缺陷,从分形维的定义出发,以降低计算量和提高精度为目的,设计出了一种基于波形长度的分形维计算方法——波长对数法。
本发明的技术方案是:
旋转机械碰摩声发射信号的分形特征分析方法,其步骤为:
首先对采样量化后的含有AE信号的噪声进行分帧处理,再将每一帧信号分成等长度的短时帧数据段,采用基于波形的分形维特征分析算法计算每一段的分形维,最后对每一帧进行中值滤波。
(1)AE信号的提取与预处理:
(1-1)通过碰摩声发射实验台接近摩擦点的第一AE传感器或位于联轴器另一端的第二AE传感器提取AE信号;
(1-2)在连续碰摩AE信号上叠加高斯白噪声和非平稳噪声;
(1-3)采样,量化,分帧;
(1-4)将每一帧含有噪声的数据等分成长为m的短时帧数据段;
(2)计算各帧的分维值:
(2-1)采用基于波形的分形维数算法即对数波长维算法,计算各短时帧数据段的分维值:
其中 δ为盒覆盖曲线的尺度,δ0为最小采样间隔,δ为最小采样间隔δ0的整数倍,即δ=kδ0,L(δ)为对应的曲线长度。
(2-2)将每段即Di-v,...,Di-1,Di,Di+1,...,Di+v按数值大小排列,取正中间序号的数据作为中值滤波的输出:
在所述的旋转机械碰摩声发射信号的分形特征分析方法中,所述的采样频率为5MHz,量化时的A/D精度为12位,分帧时每帧的帧长为512个采样点,50%的帧重叠,m必须为奇数。
在所述的旋转机械碰摩声发射信号的分形特征分析方法中,通常对更接近于真实声源信号的第一AE传感器提取的AE信号进行分析。
本发明的优点和效果在于:
1.在强噪声环境下,对数波长维比Katz维和盒维具有更强区分噪声和AE信号的能力,并且AE信号帧的对数波长维比Katz维和盒维更接近于纯净AE信号的分维,从而能够有效地检测并反映出碰摩AE事件的发生。
2.在计算复杂度、精确度的稳定性上,对数波长维也优于Katz维和盒维,从而有效地解决了采用分形理论对碰摩故障计算量大的问题。
3.利用中值滤波技术,减轻了由噪声造成的短时帧分维曲线波动较大的问题。
本发明的其他优点和效果将在下面继续描述。
附图说明
图1——碰摩声发射实验台示意图。
图2——连续碰摩时的声发射信号波形图。
图3——当δ较小时,盒内AE波形曲线图。
图4——对数波长维、Katz维和盒维的分维均值比较图。
图5——对数波长维、Katz维和盒维的计算复杂度比较图。
图6——白噪声环境(0dB)下连续碰摩AE信号的三种帧分维变化曲线图。
图7——有色噪声环境(0dB)下连续碰摩AE信号的三种帧分维变化曲线图。
具体实施方式
下面结合附图和实施例,对本发明所述的技术方案作进一步的阐述。
一.AE信号的提取与预处理
1.AE信号的提取
本试验采用3支承2跨转子***,如图1所示。第一轴承2、第二轴承6和第三轴承8都是流体动力润滑滑动轴承,摩擦点4靠近电动机1,第一AE传感器3接近摩擦点,第二AE传感器7位于联轴器5的另一端。通过第一AE传感器或第二AE传感器提取出的AE信号波形,本例中从第一AE传感器提取AE信号波形,可用于分析声信号传播的衰减及经过非连续介质耦合后信号畸变情况。在这里,我们对更接近真实声源信号的第一AE传感器3提取的信号进行分析。实验中设定采样频率为5MHz,连续碰摩AE信号波形如图2所示。
2.AE信号的预处理
在连续碰摩AE信号上叠加高斯白噪声和非平稳噪声(噪声源由英国TNO感知学会所属的荷兰RSRE研究中心提供)。对含噪信号进行采样,量化,分帧。分帧时每帧的帧长为512个采样点,50%的帧重叠。
考虑到噪声的影响可能导致各帧的分维曲线波动较大。在分帧后,将每一帧含有噪声的数据再等分成长为m的短时帧数据段(m为奇数),这样便于在计算出各短时数据帧分维值以后对其进行中值滤波。
二.基于波形的分形维方法
由分形维定义可知:
其中Nδ(F)为尺度为δ的盒覆盖曲线的最小个数,li(δ)为边长为δ的第i盒内曲线长度。
令: (式2)
分析AE信号波形,当δ趋向于0时,盒内曲线主要可以分为图3所示的几种情况。其中(a)、(b)为单调曲线,(c)、(d)为包含一个极值点的曲线,(e)、(f)包含多个极值点的曲线。为了分析的方便,当δ趋向于0时,将盒内曲线长度等效为:
k1,k2分别为至多包含一个极值点曲线和包含多个极值曲线的等效长度与盒边长之比。
则:
M为包含多个极点的盒数,则
则(式1)可改写为:
用泰勒级数展开(δ0为接近0的非常小的值):
对(式4)第2项化简可得:
(式5)第3项中,当δ0→0,L(δ0)的变化趋于缓慢,令L′(δ0)=r,当δ0→0,包含多个极点的盒数趋向于0。则:
根据(式4)、(式5)和(式6)可得:
当δ0→0时,上式中 和 均→0,进一步化简得:
为了进一步降低计算量,考虑到当δ0≈δ→0,包含多个极点的盒数M也趋向于0,可得:
其中:
(式8)即为基于波形的分形维数算法,称之为对数波长维。式中δ0可为最小采样间隔δ的整数倍,令δ0=kδ,L(δ0)为对应的曲线长度。a,b,k参数采用分形布朗曲线来确定。
由(式8)计算出各短时帧数据段的分维值,即Di-v,...,Di-1,Di,Di+1,...,Di+v,并按数值大小排列,取正中间序号的数据作为中值滤波的输出:
由于中值滤波是一种非线性滤波,不需要信号的统计特性,因此可以很好地减弱随机干扰和脉冲干扰的影响,而对频谱基本无影响。
三.性能评价
为了比较对数波长维的性能,将它与盒维和Katz维进行比较,图4为比较结果,该比较过程中的分形布朗曲线是计算机通过15次迭代产生的,共215=32786个点,程序将该曲线分帧,每帧512个点,计算每帧的分形维,并将各帧的分维值取平均得到最终结果。从图4中可以看出,对数波长维的精确度最高,其次是盒维,最后是Katz算法。
图5为对数波长维、Katz维和盒维的计算复杂度比较结果,m为中值滤波器长度,假设曲线有N个点,计算对数波长维和Katz维的步长为1,计算盒维的最小盒长也为1,共有k个点参加盒维的曲线拟合,其盒长分别为1,l1,l2,...,lk。从图5中可以看出,盒维的加法次数、乘法次数、非线性运算次数均要明显大于Katz维和对数波长维;而对数波长维在加法、乘法、比较、非线性运算方面在三者中均处于最小。
图6和图7分别给出了白噪声环境(0dB)和有色噪声环境(0dB)下连续碰摩AE信号的三种帧分维变化的曲线图,从图中可以看出:
(1)含噪AE信号中,对数波长法帧分维曲线相对比较平滑,在图6和图7中有多处分维值要明显小于噪声的分维值,可判断出这几处有AE信号产生;
(2)分析三种帧分维曲线的噪声段可看出对数波长法波动最小,而Katz维、盒维的波动比较大,可以得出对数波长维要比Katz维、盒维更稳定;
(3)含噪AE信号的信号段对数波长维明显要小于噪声段,而盒维则两者差别不大,在有色噪声环境下,盒维和Katz为几乎不能区分噪声和AE信号。
图6,7说明,在强噪声情况下,对数波长法更具有区分出噪声和碰摩AE信号的能力,并且AE信号帧的对数波长维比Katz维和盒维更接近于纯净AE信号的分维。
以上结果表明,在强噪声环境下,对数波长维能够反映碰摩AE事件的发生;对数波长维比Katz维和盒维具有更强区分噪声和AE信号的能力;在计算复杂度、精确度的稳定性上也优于Katz维和盒维。该算法的良好的性能使其为碰摩AE信号的特征识别与分析提供了一条新的途径,是一种处理碰摩AE信号比较好的分维计算和分析方法。
本发明请求保护的范围并不仅仅局限于本具体实施方式的描述。
Claims (3)
1.旋转机械碰摩声发射信号的分形特征分析方法,其步骤为:
首先对采样量化后的含有AE信号的噪声进行分帧处理,再将每一帧信号分成等长度的短时帧数据段,采用基于波形的分形维特征分析算法计算每一段的分形维,最后对每一帧进行中值滤波;
(1)AE信号的提取与预处理:
(1-1)通过碰摩声发射实验台接近摩擦点的第一AE传感器或位于联轴器另一端的第二AE传感器提取AE信号;
(1-2)在连续碰摩AE信号上叠加高斯白噪声和非平稳噪声;
(1-3)采样,量化,分帧;
(1-4)将每一帧含有噪声的数据等分成长为m的短时帧数据段;
(2)计算各帧的分维值:
(2-1)采用基于波形的分形维数算法即对数波长维算法,计算各短时帧数据段的分维值,
其中 δ为盒覆盖曲线的尺度,δ0为最小采样间隔,δ为最小采样间隔δ0的整数倍,即δ=kδ0,L(δ)为对应的曲线长度;(2-2)将每段即Di-v,...,Di-1,Di,Di+1,...,Di+v按数值大小排列,取正中间序号的数据作为中值滤波的输出:
2.根据权利要求1所述的旋转机械碰摩声发射信号的分形特征分析方法,其特征为,所述的采样频率为5MHz,量化时的A/D精度为12位,分帧时每帧的帧长为512个采样点,50%的帧重叠,m必须为奇数。
3.根据权利要求1所述的旋转机械碰摩声发射信号的分形特征分析方法,其特征为,步骤(1-1)中从第一AE传感器提取AE信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100238086A CN101251445B (zh) | 2008-04-16 | 2008-04-16 | 旋转机械碰摩声发射信号的分形特征分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100238086A CN101251445B (zh) | 2008-04-16 | 2008-04-16 | 旋转机械碰摩声发射信号的分形特征分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101251445A true CN101251445A (zh) | 2008-08-27 |
CN101251445B CN101251445B (zh) | 2010-06-30 |
Family
ID=39954946
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100238086A Expired - Fee Related CN101251445B (zh) | 2008-04-16 | 2008-04-16 | 旋转机械碰摩声发射信号的分形特征分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101251445B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101566519B (zh) * | 2009-05-22 | 2011-02-09 | 东南大学 | 基于模态波和高斯混合模型的转子碰摩声发射识别方法 |
CN101477194B (zh) * | 2009-02-17 | 2011-07-06 | 东南大学 | 一种转子碰摩声发射源定位方法 |
CN101782475B (zh) * | 2010-02-08 | 2011-07-20 | 天津工业大学 | 基于风力发电机组振动的叶片故障诊断方法 |
CN103854658A (zh) * | 2012-11-29 | 2014-06-11 | 沈阳工业大学 | 基于短时分形维数增强法的钢板腐蚀声发射信号去噪方法 |
CN106596088A (zh) * | 2016-12-13 | 2017-04-26 | 东南大学 | 基于近场声源聚焦定位的碰摩声发射故障位置识别方法 |
CN106706282A (zh) * | 2016-11-04 | 2017-05-24 | 东南大学 | 一种基于傅里叶分解的旋转机械故障诊断方法 |
CN110717274A (zh) * | 2019-10-11 | 2020-01-21 | 哈尔滨工业大学 | 一种基于分形布朗运动的噪声描述方法 |
CN113740010A (zh) * | 2021-09-07 | 2021-12-03 | 湘潭大学 | 基于挠度影响线Katz1分形维数的梁结构损伤识别方法 |
CN115436058A (zh) * | 2022-08-30 | 2022-12-06 | 苏州大学 | 一种轴承故障特征提取的方法、装置、设备以及存储介质 |
-
2008
- 2008-04-16 CN CN2008100238086A patent/CN101251445B/zh not_active Expired - Fee Related
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477194B (zh) * | 2009-02-17 | 2011-07-06 | 东南大学 | 一种转子碰摩声发射源定位方法 |
CN101566519B (zh) * | 2009-05-22 | 2011-02-09 | 东南大学 | 基于模态波和高斯混合模型的转子碰摩声发射识别方法 |
CN101782475B (zh) * | 2010-02-08 | 2011-07-20 | 天津工业大学 | 基于风力发电机组振动的叶片故障诊断方法 |
CN103854658A (zh) * | 2012-11-29 | 2014-06-11 | 沈阳工业大学 | 基于短时分形维数增强法的钢板腐蚀声发射信号去噪方法 |
CN106706282A (zh) * | 2016-11-04 | 2017-05-24 | 东南大学 | 一种基于傅里叶分解的旋转机械故障诊断方法 |
CN106596088A (zh) * | 2016-12-13 | 2017-04-26 | 东南大学 | 基于近场声源聚焦定位的碰摩声发射故障位置识别方法 |
CN110717274A (zh) * | 2019-10-11 | 2020-01-21 | 哈尔滨工业大学 | 一种基于分形布朗运动的噪声描述方法 |
CN110717274B (zh) * | 2019-10-11 | 2024-05-10 | 哈尔滨工业大学 | 一种基于分形布朗运动的噪声描述方法 |
CN113740010A (zh) * | 2021-09-07 | 2021-12-03 | 湘潭大学 | 基于挠度影响线Katz1分形维数的梁结构损伤识别方法 |
CN115436058A (zh) * | 2022-08-30 | 2022-12-06 | 苏州大学 | 一种轴承故障特征提取的方法、装置、设备以及存储介质 |
CN115436058B (zh) * | 2022-08-30 | 2023-10-03 | 苏州大学 | 一种轴承故障特征提取的方法、装置、设备以及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN101251445B (zh) | 2010-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101251445B (zh) | 旋转机械碰摩声发射信号的分形特征分析方法 | |
Liu et al. | Time-frequency representation based on robust local mean decomposition for multicomponent AM-FM signal analysis | |
Patel et al. | Defect detection in deep groove ball bearing in presence of external vibration using envelope analysis and Duffing oscillator | |
CN103575523B (zh) | 基于FastICA-谱峭度-包络谱分析的旋转机械故障诊断方法 | |
Wang et al. | Fault diagnosis of diesel engine based on adaptive wavelet packets and EEMD-fractal dimension | |
CN106096242B (zh) | 一种基于改进emd分解的尾水管压力脉动综合评价方法 | |
CN103091096A (zh) | 基于eemd和小波包变换的早期故障敏感特征提取方法 | |
CN105628381A (zh) | 一种基于改进的局部均值分解的往复压缩机轴承故障诊断方法 | |
CN105116442A (zh) | 岩性油气藏弱反射地震信号的重构方法 | |
Sheng et al. | Applications in bearing fault diagnosis of an improved Kurtogram algorithm based on flexible frequency slice wavelet transform filter bank | |
Han et al. | Roller bearing fault diagnosis based on LMD and multi-scale symbolic dynamic information entropy | |
CN113343928B (zh) | 变速路段高速铁路钢轨波磨检测方法及装置、计算机设备 | |
Shi et al. | The VMD-scale space based hoyergram and its application in rolling bearing fault diagnosis | |
CN108645620A (zh) | 一种基于信息熵和多尺度形态学的滚动轴承早期故障诊断方法 | |
Wang et al. | Fault feature extraction of hydraulic pump based on CNC de-noising and HHT | |
Zhang et al. | Fault diagnosis for gearbox based on EMD-MOMEDA | |
Wang et al. | A novel hybrid approach for damage identification of wind turbine bearing under variable speed condition | |
CN114021616A (zh) | 用于液压马达故障诊断的方法、装置及存储介质 | |
CN114325740A (zh) | 一种应用于激光多普勒测速仪的速度测量方法与*** | |
Wang et al. | Cyclic correlation density decomposition based on a sparse and low-rank model for weak fault feature extraction of rolling bearings | |
CN103674550A (zh) | 一种滚动轴承静电监测信号实时混合去噪方法 | |
CN107941511B (zh) | 一种基于信号时频分解的频率—峭度图的实现方法 | |
Jing et al. | Stability analysis in micro milling based on p‑leader multifractal method | |
US8995230B2 (en) | Method of extracting zero crossing data from full spectrum signals | |
Xu et al. | Time and frequency domain scanning fault diagnosis method based on spectral negentropy and its application |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100630 Termination date: 20140416 |