应用沙龙 | 电动零部件异响分析参数方案(二)
本帖最后由 BK声学与振动 于 2023-7-28 10:48 编辑作者:金鹏
HBK中国区应用服务经理
调制(Modulation)、包络分析(Envelope)
如果噪声中出现了明显的调幅(AM)现象,除了上篇提及的抖动度和粗糙度以外,BK Connect软件还提供包络分析(Envelope)。包络分析能够将调制信号从载波信号中分离出来。如下图所示,红色为载波信号,黑色为调制信号,经过包络分析后,即可得到右下角的调制频率fm。包络分析除了用于调制信号的解调分析以外,还适用于电机轴承的故障检测,能从轴承振动噪声信号中识别出微弱的故障信号,在出现故障的早期阶段即能发现故障。
https://picx.zhimg.com/80/v2-669c498c7f2958b6d4be54070dd65b6f_720w.jpg?source=d16d100b
如下图右图的例子中,对载波所在频段(125±100Hz)进行带通滤波后,对包络线时间历程曲线进行FFT分析。Y轴幅值为包络线的频谱幅值。在频谱中发现了明显的峰值,说明信号中存在调幅现象。在噪声信号的FFT分析结果中,如下图左图所示,我们发现在126Hz和120Hz有相邻的峰值,其中幅值更高的126Hz为载波频率,幅值较低的120Hz为调制信号频率。根据调制信号幅值和载波信号幅值的比值,可计算出调制度(Degree of Modulation),为0.082/0.395≈21%。
https://pic1.zhimg.com/80/v2-63400474a9280e0f9bc539f6de42afcc_720w.png?source=d16d100b
如果在分析之前,无法确定载波频率的频率范围,则可以先进行不同频段的包络分析,调查这些频段的数据中,是否存在调制。如下图例子所示,X轴横坐标为调制频率,左侧Y轴纵坐标为包络分析的载波频带中心频率,Z轴颜色坐标为调制信号幅值。从图中可以发现2500Hz载波频带和80Hz调制频率处有明显峰值(图中黄色圆圈位置),同时在其他载波频带,如500Hz处,也存在多个调制频率的峰值(图中橙色圆圈位置),说明这些位置存在明显的调制现象。
https://pic1.zhimg.com/80/v2-db7745554dda3cb6a856c667dd1d0963_720w.png?source=d16d100b
纯音比(Tone-to-noise Ratio)、突出比 (ProminenceRatio)、音调(Tonality)
由于电机和传动机构的结构特点,旋转机构经常发出一些纯音。这些纯音可能让电机噪声听起来尖锐刺耳、不和谐,令人烦躁。在分析纯音问题时,需要考虑宽频信号的掩蔽效应。当一个宽频信号中存在一个纯音时,纯音的声压级幅值超过纯音所在临界频带声压级幅值6个dB以上,人耳才能听到这个纯音。如果小于6dB,则人耳听不到这个纯音,它被周围频带的宽频声音所掩蔽。下图为纯音比和突出比的图例和定义。针对人耳的这种特性,BK Connect软件提供了多种纯音分析参数,包括纯音比(Tone-to-noise Ratio) 、突出比(Prominence Ratio)和音调(Tonality),这些参数能自动判断噪声中是否存在明显的纯音。
https://picx.zhimg.com/80/v2-082e091b57922de61742855357e23eea_720w.jpg?source=d16d100b
FFT、1/3倍频程(1/3 Octave)、临界频带(Critical Band)
FFT和1/3倍频程多用于表述声音能量在频域上的分布。它们的纵坐标为声音幅值,比如声压级dB、A计权声压级dB(A)等,横坐标为频率。FFT和1/3倍频程最明显区别之一是频率坐标的分辨率。FFT的频率分辨率在各个频率都是一致的(如下图左图),而1/3倍频程的频率分辨率是中心频率的固定百分比(约23%),频率越低,每个柱状图代表的频段越窄,也就是频率分辨率越小,反之频率越高,频带越宽,频率分辨率越大(如下图中图)。FFT和1/3倍频程的本质区别是噪声幅值的计算方法,FFT利用傅里叶变换原理,进行时域频域转换,得到各频率的噪声幅值。1/3倍频程采用时域滤波器的方法,对时域信号进行带通滤波,得到各频段噪声幅值。
FFT和1/3倍频程主要表述声音的能量分布,没有考虑心理声学中的频域掩蔽效应,因此与人对声音的主观感受存在一定的差距。为了减少这种差距,在进行心理声学分析时,通常使用临界频带(单位bark)表述不同频率的噪声,临界频带将人耳可听频段分为0-24bark(如下图右图)。
https://pic1.zhimg.com/80/v2-4f0c810a069e9b7f53f14a76bcd5da87_720w.png?source=d16d100b
更多参数
如果噪声信号具有非稳态特性,甚至是瞬态特性,除了上述常见的参数以外,还经常使用其他一些参数。他们由常规参数衍生出来的,这些参数也常用于评价电动零部件的异响。
百分位数
在研究噪声问题时,噪声中除了与自身形成机理有关的特征成分以外,还经常包含随机成分,这使得噪声数据经常分布在一定数值范围内。为了方便量化噪声数据的分布特性,引入统计学的一些方法,比如百分位数,常见的百分位数有第1、第50(也称为中位数)、第99百分位数等。
以连续重复测试的三组数据为例,在下图的响度时间历程中,由于被测物运行状态等因素的影响,无法保证三组数据的最大值是最稳定的、一致的。表格中为三组数据的第0至第15百分位数汇总,其中第0百分位数(第三行,即响度最大值)标准差为0.33Sone,而第3、4百分位数(第四、五行)标准差均明显小于第0百分位数的标准差,因此第3、4百分位数具有更高的一致性。
https://pic1.zhimg.com/80/v2-310d1f1637d1909520e2ce8c60bec3d7_720w.png?source=d16d100b
https://picx.zhimg.com/80/v2-ab6b42566ec87672adc0f37de1a02696_720w.png?source=d16d100b
(特征响度积分面积)百分位频率
使用临界频带表述不同频率的响度时,我们将结果称为特征响度,如下图结果,横坐标为临界频带(bark),纵坐标为响度。在研究噪声的频域分布特性时,可以利用特征响度曲线计算出积分面积。再利用积分面积百分位数所对应的频率,来表述频谱分布特性。百分位数的结果为特征频率,单位Bark。以下图为例,左右两图为两个不同声音的特征响度数据,两个声音的总响度分别为13.00 Sone和12.74 Sone,比较接近,但是它们的频谱分布是不同的,例如,在9 Bark附近,左图数值小于右图,而在16.5 Bark附近,左图数值大于右图。以第50百分位频率为例,低于某个特征频率的积分面积,占总积分面积的50%,也可以描述为,某个临界频带以上和以下特征频率积分面积相等,即特征响度“重心”,此临界频带对应的特征频率即为第50百分位频率。在上述两组数据中,左右两图的第50百分位频率分别为12.5 Bark和11 Bark(图中红色竖线),这意味着左图的特征响度“重心”比右图更偏右说明特征响度主要集中在较高频段。除了使用第50百分位频率以外,也可以使用其他百分位数(比如第70百分位数)表示特征响度的分布特性。
https://pic1.zhimg.com/80/v2-f88365d535270344a1909dfc7a9dd74f_720w.png?source=d16d100b
对于瞬态噪声信号,在瞬态噪声的时间历程中,频率分布是随时间变化的。以电机停止运行时的特征响度时域历程为例,如下图左侧图形所示,在0.09-0.25秒(蓝框内的时间范围)出现瞬态冲击噪声,在不同时刻下,特征响度是不同的。利用各个时刻的特征响度计算出百分位频率,从而反映出各个时刻的频率分布变化趋势。下图右侧图形为第70百分位频率的时间历程,其中的0秒为冲击噪声的峰值时刻,作为起始时刻。第70百分位频率曲线随时间逐步降低,意味着特征响度中的频率成分,随着时间逐渐向更低频率转移。
https://picx.zhimg.com/80/v2-5de734e5f0c8c1c09e703dc5925adc9e_720w.png?source=d16d100b
颤声(Warble)
当电机处于匀速运行时,由于负载或阻力的影响,容易引起噪声的低频颤动,如下图所示,在响度时间历程数据中,从2.4秒至5.1秒(蓝框内),响度值出现了明显的低频颤动。如果此类颤动的频率在2-8Hz(也有文献规定为1-10Hz)时,通常将之称为颤音(Warble)。下图的颤音现象中,包含了多个周期的极大值和极小值,如图所示的6个周期的极值。
https://picx.zhimg.com/80/v2-867ac8e30ffdc353945ad869595868e4_720w.png?source=d16d100b
未完待续......
下期文章中,我们将继续聊聊其他更多参数,敬请期待。
全新《HBK电声测量解决方案》已新鲜出炉点击这里,查看 / 下载完整解决方案
您还可以通过如下方式联系我们,了解更多产品与应用详情:
邮箱:cn.info@hbkworld.com
网址:www.bksv.com/zh
电话:400-900-3165(周一至周五9:00-18:00)
点击这里,咨询B&K产品信息:https://www.bksv.com/zh/request-a-quote
页:
[1]