基于波叠加方法的半自由声场全息理论

《基于波叠加方法的半自由声场全息理论》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、http://www.paper.edu.cn基于波叠加方法的半自由声场全息理论李卫兵陈剑于飞陈心昭合肥工业大学噪声振动工程研究所（230009）hf_lwb@sina.com摘要：在半自由声场中，实际全息测量声压为全息面上的直达声压和反射声压叠加；而常规声全息技术要求全息面声压只包含直达声压，这样就不能直接用常规全息方法来重建与预测半自由声场。以波叠加方法为全息变换算法，在充分考虑反射声压的情况下，建立了基于波叠加方法的半自由声场全息重建与预测理论，解决了半自由声场全息重建与预测问题，拓宽了全息技术的应用范围。

2、数值仿真的结果充分证明了基于波叠加方法的半自由声场全息理论的正确性和可行性，以及常规全息技术在半自由声场重建与预测过程中的局限性。关键词：近场声全息半自由声场波叠加方法1.引言上世纪80年代初，美国宾夕法尼亚大学学者E.Ｇ.Ｗilliams等提出了基于空间声场变[1-2]换的近场声全息。近场声全息是在紧靠被测声源物理表面的测量面上记录全息数据，然后通过变换技术重建三维空间声压场、振速场、声强矢量场，并能预报远场指向性。由于是近场测量，所以除了记录了传播波成分外，还能记录随传播距离按指数规律衰减的倏逝波成分，由于倏

3、逝波含有振动体细节信息，所以理论上可获得不受波长限制的高分辨率图像，测[1-8]量覆盖了从源出来的一个大的方位角，有指向性的源也能够被不失信息地检测出来。声全息是一种有效而快捷的噪声源辨识技术，只需要测量面上的复声压数据，就可以在很宽的频带范围内对声源特性进行研究。它对大型复杂结构的振动和噪声辐射特性研究、噪声源的识别与定位以及结构强度评价都是一种极为有效的方法，有助于对结构振动、噪声进行有效控制，在工程上具有很高的应用价值和应用前景。由于常规声全息技术只适用于自由声场，所以全息测量面上存在反射声的问题严重限制了

4、全息技术的应用。在文献[9]所提到的三个亟待解决的问题中就包括了全息面测量声压中包含反射声压的问题。针对这个问题，国内外许多研究者在实验中都采取一些措施来削弱反[1,3,4][5]射声的影响，比如在全消声室中进行测量，或者通过挡板将地面反射声与直达声隔开，[6]或将声源放置在离地面很高的地方进行测量。虽然这些方法对声源定位有一定作用，但是并不能准确地预测整个声场的辐射特性，给声源特性判别带来不便，不利于进行噪声源的控制。在声辐射问题中,为了寻求边界元方法的有效替代方法,Koopmann等提出了更容易理解[10-1

5、1]和实施的波叠加方法来计算声辐射问题.波叠加方法的基本思想是：任何物体辐射的声场都可以由置于该辐射体内部的，若干个不同大小源强的简单源产生的声波场叠加得到。本文以波叠加方法作为全息变换算法，在充分考虑反射声压的基础上，建立了反射面为刚性和-1-http://www.paper.edu.cn非刚性条件下的半自由声场全息重建与预测模型，解决了半自由声场环境下的全息重建与预测问题，拓宽了全息技术的应用范围。数值仿真的结果验证了该半自由声场全息理论的正确性、有效性和可行性，以及常规全息技术在半自由声场全息重建与预测过程

6、中的局限性。1.波叠加积分方程−iωt在理想流体媒质微小扰动形成的时谐声场中,去除时间的相关性e后,声场中任意一点r上的复声压p(r)必满足Helmholtz方程：22∇p(r)+kp(r)=0（1）式中p(r)为点r上的复声压；k=ω/c=2π/λ为声波数，c为声速，λ为声波长，ω为角频率。在如图1所示的声辐射外问题中，S是声辐射体的闭合表面，其外部区域记为E，内部区域记为D。域E中点r上的声压，可以通过解方程（1）得到∂gr,(r)SC(r)p(r)=[p(r)−iρcku(r)gr,(r)]dS（2）E∫S

7、SSSrS∂nrS式中，ikrg(r,r)=4/1(πr)e,r=

8、r−r

9、（3）SS当点r分别在D,S,E上时,系数C(r)分别为0,0.5,1；u(r)为边界表面点r上ESS的法向振速，n为点r处的外法线矢量。rSSrsrroΩn图.1辐射体与声场各个域之间的位置关系图在图1所示的内声辐射问题中,边界面仍为S,假设在域D内有一个连续分布的声源体Ω,此时应用质量守恒定律,可以得到一个修正的Helmholtz方程,22∇p(r)+kp(r)=iρckqo(r)（4）式中,⎧q(r),r∈Ωqo(r)=⎨（5）⎩,

10、0r∉ΩIr∈D于是可以得到方程（4）的解,即域D中点r上的复声压为∂gr,(r)SC(r)p(r)=[−p(r)+iρcku(r)gr,(r)]dS+iρckq(r)g(r,rd)Ω（6）I∫SS∂nSSrS∫ΩoorS-2-http://www.paper.edu.cn式中,当点r分别在D、S、E上时,系数C(r)分别取值为1、0.5、0。I由于上述的两种辐射问题对应