近场声全息方法识别噪声源的实验研究

近场声全息方法识别噪声源的实验研究Ξ

于　飞　陈　剑　李卫兵　陈心昭

(合肥工业大学机械与汽车工程学院　合肥,230009)

摘　要　根据近场声全息(NA H)的原理,建立了全息实验所需要的采集、分析系统。针对影响重建精度较大的截止波数的选取问题,给出了较为详细的讨论,并提出一种不需先验知识的截止波数选取方法。最后通过对实测数据进行全息变换,重建结果表明:在采用提出的截止滤波选取方法后,NA H技术可以精确地对噪声源进行定位与识别,并且可以得到三维空间内的声压、质点振速和声强矢量等声学信息。

关键词:声源识别;近场声全息;实验研究;截止波数

中图分类号:TB532;TB533+.2

进行空间声场的可视化和噪声源的识别与定

位,对于噪声测量和控制工程具有非常重要的意义。上世纪80年代初提出的近场声全息技术(NA H),便是可视化空间声场和定位噪声源的一种强有力工具。近场声全息可以由一个测量面的声压标量数据,反演和预测另一面上的声压、质点振速、矢量声强等重要声场参量,受到了各国研究人员及一些相关公司的重视。近场声全息技术真正地将丰富的声学理论同噪声测量、控制工程紧密地结合起来[1～2]。20世纪80年代末,国内一些学者逐渐对此方法进行了研究:中科院武汉物理所对编磬表面振动模态做了研究[3～4];哈尔滨工程大学对基于边界元法的水下近场声全息也做了研究[5];清华大学汽车工程系对非近场声全息确定噪声源进行了研究[6～7];合肥工业大学机械工程学院对近场声全息方法识别噪声源作了一定的研究[8～9]。

近场声全息可以不受波长分辨率限制重建声场,但在此种全息过程中截止波数的选取对重建分辨率的影响非常大。文献[3]提出一种需要测量先验知识的优化滤波方法,而这种先验知识一般是不易获得的。本文根据截止波数的大小对重建结果的影响趋势,提出一种不需要先验和后验知识的截止波数选取方法。并根据近场声全息的原理,建立了全息实验所需要的采集、分析系统。采用提出的滤波参数选取方法后,对数据进行全息变换,得到了令人满意的重建结果。该优化截止波数选取方法的提出,有助于在实际工程中推进近场声全息技术在高分辨率识别噪声源、可视化声场等方面的应用。1　理论背景

由文献[1,8]可知,在稳态的三维空间声场中,一个平面(全息面)上声压的波数谱与另一个更靠近声源的平行面(声源面或重建面)上声压和质点法向振速的波数谱之间的关系为

P(k x,k y,z S)=P(k x,k y,z H)e-i k z(z H-z S)(1) V(k x,k y,z S)=k z P(k x,k y,z H)e-i k z(z H-z S) Θ0ck(2)式中　z H和z S分别为全息面和重建面的z坐标;k 为声波数;k x和k y分别为对应坐标x和y的波数;而k z与波数k x,k y之间的关系为

当k2x+k2y≤k2时

k z=k2-(k2x+k2y)(3)当k2x+k2y>k2时

k z=i(k2x+k2y)-k2(4) k z取值为式(3)时,对应的声波传播方式是以幅值不变、相位改变的传播波方式传播;当取值为式(4)时,对应的声波传播方式是以相位不变、幅值减小的倏逝波方式传播。倏逝波随全息面与重建面之间距离的增加,成指数倍地迅速衰减,对应的是高波数成分的声波。在非近场的声全息中,由于测量点位置与声源面之间距离过大造成倏逝波信息的丢失或被测量噪声所掩盖,全息重建的结果也就失去高频信息,这种高频信息类似于小波变换处理图像中的细节信息。

近场声全息技术除了能够由全息声压数据重建源面上的声压和法向振速之外,由Eu ler公式还能

第17卷第4期2004年12月

振　动　工　程　学　报

Jou rnal of V ib rati on Engineering

V o l.17N o.4

D ec.2004

Ξ国家自然科学基金资助项目(编号:50275044)及高等学校博士点科研基金资助项目(编号:20020359005)收稿日期:2004203203;修改稿收到日期:2004205231

够得到源面上的二维切向振速分布

U (k x ,k y ,z S )=k x P (k x ,k y ,z H )e

-i k z (z H -z S )

Θ0ck

(5)

W (k x ,k y ,z S )=k y P (k x ,k y ,z H )e

-i k z (z H -z S )

Θ0ck

(6)

式中　U 和W 分别代表声源面上振速x ,y 方向分

量u 和w 的波数谱。声强矢量是表征空间中声能量流向的重要物理量,经常用于噪声源的定位和识别。其中,有功声强表示声源各点对外声辐射强度的大小及声能量的流动方向,其定义为

I =1

2R e[p v 3]

(7)式中　“*”表示取共轭,R e 表示取实部。

2　实验测量

全息测量实验在半消声室中进行,实验室本底噪声小于28dB ,测量装置在消声室中的实物照片如图1所示。实验所用到的主要设备为:尺寸为60c m ×14c m ×1512c m 的音箱作为目标声源,16通道的信号采集器,16通道的信号调理器,传声器等等。全息面的大小为113m ×113m 的面,全息面距离音箱声源面的距离为12c m ;沿水平方向(x 轴)和竖直方向(y 轴)的全息测量间隔均为5c m ,同时以音箱两纸盆间的中间点为坐标系原点

。

图1　近场声全息测量系统在半消声室中的照片

通过信号发生设备,产生主频为400H z 的信号,经信号放大器处理后,传输到目标声源使之发出声音。首先在声场中布置一个自始至终不动的参考传声器,另一传声器线阵在全息网格点上做扫描式测量。最后由测得的参考传声器和扫描传声器信号,经过信号调理器的放大处理,由采集器采集到这些经过处理后的声信号。采集器的采样频率为5kH z ,每块采样长度为1024点的时域数据,对全息面上

每一点的声压信号连续采样32块。

假设参考传声器采集的时域声压信号为p r (t ),

扫描传声器采集的声压为p h (t )。为获得全息面上的复声压,可采用参考传声器互谱法获取全息面上各点的声压相位Η(Ξ),即

Η(Ξ)=arg [S r ,h (Ξ)

S r ,r (Ξ)]

(8)式中　S r ,h (Ξ)为参考声压p r (t )与扫描声压p h (t )的互谱;S r ,r (Ξ)为参考声压p r (t )的自谱;arg 表示取相位。全息面上的声压幅值可以通过扫描声压p h (t )的自谱得到

p (Ξ)=

S h ,h (Ξ)

(9)

由式(8)～(9),

最终可以得到全息测量面上的复声压的幅值和相位分布,结果如图2所示。

3　分析处理

由采集得到的声时域信号,经过自、互谱处理得到全息面上的复声压,再通过声全息变换可以得到

3

64第4期于　飞等:近场声全息方法识别噪声源的实验研究