基于HOA的圆柱腔低频声场重构

王岩, 陈克安, 胥健

(西北工业大学 航海学院, 陕西 西安 710072)

封闭空间声场重构在噪声控制方案验证、产品声学设计等方面具有重要的应用，其中圆柱结构作为飞机、潜艇和空间站等航行器舱室结构的典型代表，研究其低频声场重构具有重要意义。例如对于飞机、潜艇等舱室结构的有源噪声控制[1]，降噪方案的有效性和降噪效果往往需要在实际状态下进行验证，如果能够利用扬声器在实验室舱段模型内重构真实运行状态下的舱内声场，这样就可以方便地评估各种降噪方案，从而极大地节省研发及测试成本，提高工作效率。目前已有的方法是在重构区域布置大量传声器采集声场信息，在最小二乘意义下对空间格林函数求伪逆得到次级声源的驱动函数，这称为声学逆问题方法[2]。近些年有学者重新整理并丰富了该理论，讨论了求解唯一性及不适定性等问题[3]，提出基于波束形成的正则化方法并用于飞机舱室声场重构[4-5]。该方法的缺点是只能有效重构声阵列附近的区域，要求布置大量的传感器并且会对原始声场产生干扰。

高阶Ambisonics[6-8]是主要应用于厅堂等空间声场的重构技术，其基本原理是将声场分解为一系列球谐函数(基函数)及其展开系数的乘积，根据模态匹配原理(mode matching)求解次级声源的驱动函数，进而重构三维声场。目前相关研究都以自由空间声传播特性为理论基础展开，对于三维自由空间格林函数，它可以理解为点声源的空间-时间传递函数，在重构过程中其波数域的谱分量呈现出空间低通特性，能够抑制次级声源离散化导致的混叠效应[8]。在实际房间内重构空间声场时，学者们提出了主动房间补偿滤波器[9]、房间传递函数参数化方法[10]以及利用高阶指向性次级声源[11]等方法减小壁面声反射的影响，并取得了较好的重构效果。而对于小尺度空间低频声场，目前有关基于HOA的声场重构研究还鲜有报道。对于圆柱声腔，前期研究表明其声模态可以由球阵有效地进行球谐函数分解，获取其展开系数[12]，并依据其稀疏特性，利用压缩感知技术由较少的空间测量点求解高阶展开系数[13]，这为本文利用HOA重构圆柱腔内声场奠定了良好的基础。

为了克服传统声学逆问题方法的缺陷，在整个三维空间内更好地重构低频声场的时间和空间特性，本文以圆柱结构为研究对象，提出利用HOA方法重构小尺度空间低频声场，并与声学逆问题方法做对比。由于圆柱腔内低频声场是有限个声模态的叠加，因此腔内次级声源的声传播特性将与以往自由场和房间扩散声场有很大不同，主要体现在波数域内格林函数的谱分量受声模态的空间分布及其阶数影响。本文以声模态理论为基础重新推导重构算法，重点讨论声腔内格林函数谱分量与次级声源数目、重构参数kr以及腔内声模态的关系，定义声场重构滤波器并讨论其特性。在此基础上，对比分析HOA与声学逆问题方法，从声场采集、影响因素及重构效果等方面进行讨论，最后通过实验验证本文提出方法的有效性及优势。

1 腔内HOA声场重构理论

首先介绍HOA声场重构的基本原理，根据圆柱声腔内声传播特性，推导腔内格林函数的谱分量表达式，进而推导次级声源驱动函数的最小二乘解。

1.1 基本原理

假设刚性壁圆柱声腔在点声源激励下形成稳态声场，腔内存在封闭曲面边界∂V,如果边界∂V包围的封闭空间V内不存在声源,那么V内的声场可以用单层势来描述[7]

P(kr)=∮∂VD(kr0)GN(kr|kr0)dS0

(1)

式中，r=(r,φ,z)表示V内任意一点,r0=(r0,φ0,z0)为边界∂V上的点,D(kr0)称为次级声源的驱动函数,GN(kr|kr0)为圆柱声腔的纽曼格林函数[14]。在球坐标系下,将重构区域V内部的声场P(kr),∂V上的次级声源驱动函数D(kr0)以及纽曼格林函数GN(kr|kr0)进行球谐函数展开[14]

(5)

根据模态匹配方法(mode-matching method)可以推导出驱动函数的谱分量为

(6)

将(6)式带入(3)式,即可求出次级声源的驱动函数D(kr0),进而利用(1)式重构三维声场。

(7)

式中，Hn(kr,ka)为传递因子[14]

Hn(kr,ka)=

(8)

圆柱腔纽曼格林函数GN(kr|kr0)的球傅里叶变换可以写为

(9)

它表示在波数域次级源与重构点之间谱分量的传递关系。根据圆柱腔GN(kr|kr0)的表达式[14],可以计算出格林函数的谱分量为

(10)

D=(G*G)-1G·H·P

(11)

2 腔内HOA声场重构实现

根据圆柱腔内低频声场HOA重构理论,讨论其实现过程。假设有一刚性壁圆柱壳体,长3 m,半径0.8 m,厚3 mm,圆柱左右两端面分为位于z=0 m和z=3 m处,声腔中轴线为z轴。在圆柱内部(r=0.3,φ=π,z=0.2)处有一个点声源在腔内形成稳态声场。为了便于讨论腔内声场重构特性,本文根据球面Fliege离散点分布形式布放次级声源。重构区域是中心为(x=0,y=0,z=1.5)m,半径r=0.6 m的球面包围的内部区域。

2.1 声压谱分量

图1 重构半径r=0.4 m的球面声压谱分量

2.2 格林函数谱分量

格林函数谱分量的计算与次级声源的数目、重构频率与重构半径(即kr)以及声腔内声模态有关。在实际声场重构时,有限的次级声源数目会影响高阶谱分量的计算。选取球面分布49个Fliege离散点,计算n=0～7阶谱分量,如图2所示。

图2 圆柱声腔格林函数谱分量的幅值

2.3 重构滤波器与最小二乘解

图3 重构滤波器的滤波特性

3 腔内HOA与声学逆问题的对比

求解次级声源驱动函数的传统方法是声学逆问题方法,它利用传声器阵列在重构区域采集声场信息,根据Tikhonov正则化求解次级声源驱动函数的最小二乘解。本节将对这2种重构方法进行对比分析。

在声腔中心处半径为0.6 m的球面上按照Fliege离散点分布16个次级声源进行HOA声场重构,声场最高截断阶数N=4,并与声学逆问题方法做对比。对于HOA方法,利用声腔中心处半径为0.1 m的球阵采集声场信息,球阵表面分布49个传声器,计算次级声源的驱动函数。对于声学逆问题方法,在重构区域内X=0平面上均匀分布108个传声器采集声场信息,根据Tikhonov正则化求解次级声源驱动函数的最小二乘解。

根据这2种方法计算的驱动函数重构声场分布。图4给出了2种方法的重构效果对比,选择重构区域内X=0的横截面与Z=L/2的纵截面为重构面,即图中虚线圆圈内。可以看出,X=0的横截面上,2种方法都能较好地重构出真实声场,重构的相对误差如图5a)所示,只有在声模态节线位置重构误差较大,在整个重构面上HOA方法的重构相对误差比声学逆问题的小。对于Z=L/2的纵截面上,HOA方法依然能够较好地重构声场,只有在声模态节线位置重构误差较大,而声学逆问题方法的重构误差就更大,如图5b)所示。由此可见,HOA方法在圆柱腔低频声场的重构精度更高。

图4 重构声场与理论声场的声压幅值比较

图5 2种方法声压幅值重构的相对误差

4 实验验证

在半消声室内,利用HOA方法开展圆柱形飞机舱室低频声场重构的实验研究。舱室模型长度为6.4 m,直径为2.67 m,底板距离顶部最高点2 m,将一个十二面体声源放置于角落,由白噪声激励产生低频声场。将半径为0.1 m,球面均匀分布64个传声器的刚性球阵放置于声腔中心处进行声场测量,并将球面声压分解为一系列球谐函数及其展开系数,如图6所示。以声腔中心为圆心,半径为0.35 m的球面上分布16个次级声源,计算次级声源的驱动函数,并重构224 Hz时声腔中心水平面的声压分布(x=0 m,y=-0.35～0.35 m,z=-2.2～1.3 m),如图7所示。可以看出,HOA方法的重构结果与实际测量值有较好的吻合,而声学逆问题的重构声压幅值明显比实际测量值低,这表明HOA方法在圆柱低频声场重构中能够取得较好的效果,优于声学逆问题方法。

图6 实验模型及刚性球阵

图7 声腔水平面声压重构值与理论值的对比

5 结 论

本文提出利用HOA重构圆柱腔低频声场，并与传统的声学逆问题方法做对比。主要结论如下：

1) 利用球阵进行声场采集并将腔内声场分解为球谐函数的形式，得到的球面声压谱分量包含了腔内声场信息，因此可以利用该谱分量重构内部声场。这种声场采集方式避免了传统声阵列需要大量的传声器，减少了对原始声场的干扰。

2) 声腔内格林函数的谱分量受声模态分布影响，重构区域中心主要存在低阶谱分量并主要由低阶声模态贡献，而远离重构区域中心存在高阶谱分量并主要由高阶声模态贡献，在重构频率附近的声模态是谱分量最主要贡献者，因此在实际声场重构时，要根据重构参数kr合理选择声模态阶数，布置次级声源时要避开主要声模态的节线处。

3) 影响HOA重构精度的主要参数是截断阶数 。在声场采集时，截断阶数越高，获取的声场信息越丰富，但在声场重构时，重构滤波器会放大球面声压谱分量，且阶数越高，重构区域越大，放大作用越明显，这意味着由各种测量噪声引起的谱分量误差也会被严重放大。因此需要合理选择最优的截断阶数，即保证重构精度又满足求解稳定性。

4) 对点声源激励下的圆柱腔内声场进行仿真分析，结果表明HOA能够较好地重构整个三维声场，从声场采集、声场重构的精度及区域大小等方面都优于声学逆方法。

5) 在圆柱形舱室模型内对本文方法进行了实验验证，并与声学逆问题做对比。实验结果表明，HOA方法能够获得比逆问题更好的重构效果，具有良好的实际应用价值。

6) 最后需要指出的是，本文仅针对单个点声源激励下的圆柱腔内声场进行了仿真和实验研究，且并未考虑模态重叠等情况。对于实际中多声源激励，或结构辐射声等情况，还需进一步研究。