单层传声器阵列原位测量半空间边界表面声阻抗的方法

：
1.本发明涉及利用单层传声器阵列原位采集半空间边界近场的全息声压，基于近场声全息和虚源原理重构计算边界表面声学量，进而获得边界表面声阻抗的方法，属于声学原位测量、全息测量、声场成像和近场声全息技术领域。


背景技术：

2.准确获取声场物理边界的表面声阻抗等声学特性参数是对其附近声场和声源进行正确控制和利用的前提。阻抗管法和混响室法是测量小尺寸和大尺寸材料样本声阻抗的标准方法，其实施过程需要具备特制的声学测量设备与标准的声学环境，测量结果不能准确反映被测材料在实际的应用环境中的声学特性和声学参数。
3.利用传声器阵列在半空间边界近场采集全息声压，基于近场声全息和虚源原理重构计算声场声学量，进而获得边界表面声阻抗，是在实际工程环境中原位测量边界表面声阻抗的有效方法。
4.而基于傅里叶声学法、统计最优近场声全息和等效源法的边界表面声阻抗测量方法，需要利用双层平面传声器阵列或者球面传声器阵列采集全息声压；基于逆边界元法的测量方法，则需要分别布置包围声源和包围边界面的两层传声器阵列。换言之，这些方法或者需要更多的传声器阵列，如双层传声器阵列；或者需要特殊几何形状的传声器阵列，如球面传声器阵列。


技术实现要素：

5.本发明克服上述边界表面声阻抗原位测量方法的局限性，提出只需要单层传声器阵列的原位测量边界表面声阻抗的方法。
6.本发明在声源中心及其关于边界表面的镜像点处，各使用一组球面波基函数的叠加来分别表述声源辐射与边界反射声场，建立边界近场的半空间声场的数学模型，构造半空间场点声压与边界表面声压、边界流体介质质点法向振速之间的传递函数，利用单层传声器阵列采集的声压分布，重构边界表面声压和法向振速分布。
7.本发明根据重构的边界表面声压和法向振速分布，计算边界表面声阻抗，也可以更进一步计算反射系数和吸声系数等其它表征边界声学特性的参数。
8.本发明所考察的边界表面为平面，通过调整声源方位，改变声波入射角，进而计算不同入射角条件下的边界表面声阻抗。
9.本发明的单层传声器阵列原位测量半空间边界表面声阻抗的方法，包含以下内容：
10.s1.在边界的半空间一侧布置声源和单层传声器阵列，进行声压全息测量；
11.在半空间中布置声源，将声源中心记为o1；在边界表面近场，平行于边界表面布置单层的平面传声器阵列，阵列面与边界之间的距离为ha；将阵列面中心在边界的投影记为点o，将点o和声源中心o1的连线与点o处边界法线之间的夹角定为声波入射角θ
inc
；布置完
成之后，启动声源发出声波，单层传声器阵列采集全息声压；采集完成之后，关闭声源，调整声源位置以改变声波入射角θ
inc
；重复进行全息声压采集，获取不同入射角θ
inc
条件下的全息声压；
12.当流体介质为水介质时，声源和传声器阵列更换为水介质中的声源和水听器阵列；
13.s2.建立边界一侧半空间声场的数学模型；
14.以点o为坐标原点，边界所在平面为x-y坐标平面，建立全局坐标系；将o1关于边界的镜像点记为o2，分别以o1和o2为坐标原点建立局部坐标系，将半空间场点x在两局部坐标系中的坐标分别记为x1≡(r1,θ1,φ1)和x2≡(r2,θ2,φ2)；对于稳态声场，在两局部坐标系中，各以一组球面波函数的线性叠加来表述声源辐射和边界反射在该场点处的声压响应，其表达式分别为：
[0015][0016]
和
[0017][0018]
其中，ψ
1,j
(x1；ω)和ψ
2,j
(x2；ω)分别为以o1和o2为坐标原点的球面波基函数；c
1,j
(ω)和c
2,j
(ω)分别为两组基函数展开项系数；ω为声波角频率；j为基函数展开项序数；j为基函数展开项项数；在球面坐标系中，ψj的表达式为：
[0019][0020]
其中，为第一类球汉克尔函数，k＝ω/c为声波波数，c为声速；为球谐函数；在式(1)～式(3)中，整数n，l和j满足关系式j＝n2+n+l+1，其中，-n≤l≤n，0≤n≤n，n为n的截断值；
[0021]
场点x处的半空间声压p
half
(x；ω)为声源辐射与边界反射的线性叠加，可表示为：
[0022][0023]
s3.以部分测点的全息测量值作为输入，重构其余测点的声压值，以声压重构误差为最小为原则确定最优基函数展开项数；
[0024]
将阵列的声压测点坐标记为m为声压测点数目；根据式(4)，阵列采集的一组声压值可以表示为如下矩阵形式：
[0025][0026]
其中，球面波函数展开项矩阵为：
[0027][0028]
系数向量{c(ω)}
2j
×1为：
[0029][0030]
其中，上标t表示向量的转置；
[0031]
按照隔点取点的方式，将声压测点分为两组，第一组测点坐标记为记为记为第二组测点坐标记为第二组测点坐标记为其中，和分别表示向上取整和向下取整；根据式(5)，建立第二组测点的声压重构值与第一组测点的声压测量值之间的数学关系：
[0032][0033]
其中，为第二组测点的声压重构值组成的向量，为第一组测点的声压测量值到第二组测点的声压重构值的传递矩阵：
[0034][0035]
其中，上标表示对矩阵求伪逆：
[0036][0037]
其中，上标h表示矩阵的共轭转置；
[0038]
设定基函数展开项数j的取值上限为j
max
，即1≤j≤j
max
；对于该范围内的任意j，利用式(8)～式(10)，重构第二组测点的声压值，并计算第二组测点的声压重构值与测量值之间的相对误差：
[0039][0040]
其中，||
·
||2为向量2-范数；从1～j
max
遍历所有的j，将ε最小值对应的展开项数确定为最优展开项数j
opt
；
[0041]
s4.建立全息测量声压与边界表面声压的数学关系，重构边界表面声压；
[0042]
设定基函数展开项数为j
opt
，则式(5)化为如下形式：
[0043][0044]
根据式(12)，建立全息测量声压与边界表面声压的数学关系，则边界表面的声压可由下式重构计算：
[0045][0046]
其中，为边界表面的重构点坐标，s＝1,2,
…
,s，s为重构点数目，为全息测量声压到边界表面声压的传递矩阵：
[0047][0048]
其中，为两组球面波函数在第s个重构点的展开项组成的矩阵，矩阵元素的布置形式参见式(6)；
[0049]
s5.建立全息测量声压与边界表面流体介质质点法向振速的数学关系，重构边界质点法向振速；
[0050]
场点x处的声压p
half
(x；ω)与质点振速v
half
(x；ω)满足欧拉方程：
[0051][0052]
其中，ρ0为流体介质密度，为梯度算子，根据式(12)和式(15)，建立全息测量声压与边界表面流体介质质点法向振速的数学关系，边界表面的质点法向振速可由下式重构计算：
[0053][0054]
其中，为全息测量声压到边界表面质点法向振速的传递矩阵：
[0055][0056]
其中，n为边界表面在重构点的单位法向矢量；
[0057]
s6.计算边界表面声阻抗率；
[0058]
边界表面的声阻抗率定义为边界表面声压p
half
(x；ω)和质点法向振速v
nhalf
(x；ω)的比值；利用重构的一组声压和法向振速计算各重构点处的表面声阻抗率，并取平均值，得到边界表面声阻抗率为：
[0059][0060]
本发明的有益效果：
[0061]
1.本发明基于近场声全息结合传声器阵列原位采集物理边界近场声压信息，实现了边界表面声阻抗的原位测量与计算，从而不再需要依赖标准的阻抗测量设备和测量环境来实施边界材料样本声阻抗的测量，同时，本发明获取的边界表面声阻抗更加符合边界材料在实际的工程应用环境中的声学特性。
[0062]
2.本发明利用单层传声器阵列采集全息声压信息，相比于其它利用双层传声器阵列或球面传声器阵列的测量方法，降低了传声器阵列的设计、加工和应用于原位全息声压
测量时在现场布置和安装的难度。
附图说明：
[0063]
图1为半空间中声源和传声器阵列布置示意图；
[0064]
图2为确定最优球面波基函数展开项数时测点分组示意图；
[0065]
图3为边界表面声阻抗率计算值与理论值的对比。
具体实施方式
[0066]
下面结合附图具体说明本发明的实施方式。
[0067]
本发明的单层传声器阵列原位测量半空间边界表面声阻抗的方法的实施，按以下步骤进行：
[0068]
步骤1，在边界的半空间一侧布置声源和单层传声器阵列，进行声压全息测量。
[0069]
如图1所示，在半空间中布置声源，将声源中心记为o1。在边界表面近场，平行于边界表面布置单层的平面传声器阵列，阵列面与边界之间的距离为ha。将阵列面中心在边界的投影记为点o，将点o和声源中心o1的连线与点o处边界法线之间的夹角定为声波入射角θ
inc
。布置完成之后，启动声源发出声波，单层传声器阵列采集全息声压；采集完成之后，关闭声源，调整声源位置以改变声波入射角θ
inc
；重复进行全息声压采集，获取不同入射角θ
inc
条件下的全息声压。
[0070]
当流体介质为水介质时，声源和传声器阵列更换为水介质中的声源和水听器阵列。
[0071]
步骤2，建立边界一侧半空间声场的数学模型。
[0072]
如图1所示，以点o为坐标原点，边界所在平面为x-y坐标平面，建立全局坐标系。将o1关于边界的镜像点记为o2，分别以o1和o2为坐标原点建立局部坐标系，将半空间场点x在两局部坐标系中的坐标分别记为x1≡(r1,θ1,φ1)和x2≡(r2,θ2,φ2)。对于稳态声场，在两局部坐标系中，各以一组球面波函数的线性叠加来表述声源辐射和边界反射在该场点处的声压响应，其表达式分别为：
[0073][0074]
和
[0075][0076]
其中，ψ
1,j
(x1；ω)和ψ
2,j
(x2；ω)分别为以o1和o2为坐标原点的球面波基函数；c
1,j
(ω)和c
2,j
(ω)分别为两组基函数展开项系数；ω为声波角频率；j为基函数展开项序数；j为基函数展开项项数。在球面坐标系中，ψj的表达式为：
[0077][0078]
其中，为第一类球汉克尔函数，k＝ω/c为声波波数，c为声速；为球谐函数。在式(1)～式(3)中，整数n，l和j满足关系式j＝n2+n+l+1，其中，-n≤l≤n，0≤n≤n，n为n的截断值。
[0079]
场点x处的半空间声压p
half
(x；ω)为声源辐射与边界反射的线性叠加，可表示为：
[0080][0081]
步骤3，以部分测点的全息测量值作为输入，重构其余测点的声压值，以声压重构误差为最小为原则确定最优基函数展开项数。
[0082]
将阵列的声压测点坐标记为m为声压测点数目。根据式(4)，阵列采集的一组声压值可以表示为如下矩阵形式：
[0083][0084]
其中，球面波函数展开项矩阵为：
[0085][0086]
系数向量{c(ω)}
2j
×1为：
[0087][0088]
其中，上标t表示向量的转置。
[0089]
如图2所示，按照隔点取点的方式，将声压测点分为两组，第一组测点坐标记为如图2所示，按照隔点取点的方式，将声压测点分为两组，第一组测点坐标记为第二组测点坐标记为第二组测点坐标记为第二组测点坐标记为其中，和分别表示向上取整和向下取整。根据式(5)，建立第二组测点的声压重构值与第一组测点的声压测量值之间的数学关系：
[0090][0091]
其中，为第二组测点的声压重构值组成的向量，为第一组测点的声压测量值到第二组测点的声压重构值的传递矩阵：
[0092][0093]
其中，上标表示对矩阵求伪逆：
[0094][0095]
其中，上标h表示矩阵的共轭转置；
[0096]
设定基函数展开项数j的取值上限为j
max
，即1≤j≤j
max
。对于该范围内的任意j，利
用式(8)～式(10)，重构第二组测点的声压值，并计算第二组测点的声压重构值与测量值之间的相对误差：
[0097][0098]
其中，||
·
||2为向量2-范数。从1～j
max
遍历所有的j，将ε最小值对应的展开项数确定为最优展开项数j
opt
。
[0099]
步骤4，建立全息测量声压与边界表面声压的数学关系，重构边界表面声压。
[0100]
设定基函数展开项数为j
opt
，则式(5)化为如下形式：
[0101][0102]
根据式(12)，建立全息测量声压与边界表面声压的数学关系，则边界表面的声压可由下式重构计算：
[0103][0104]
其中，为边界表面的重构点坐标，s＝1,2,
…
,s，s为重构点数目，为全息测量声压到边界表面声压的传递矩阵：
[0105][0106]
其中，为两组球面波函数在第s个重构点的展开项组成的矩阵，矩阵元素的布置形式参见式(6)。
[0107]
步骤5，建立全息测量声压与边界表面流体介质质点法向振速的数学关系，重构边界质点法向振速。
[0108]
场点x处的声压p
half
(x；ω)与质点振速v
half
(x；ω)满足欧拉方程：
[0109][0110]
其中，ρ0为流体介质密度，为梯度算子，根据式(12)和式(15)，建立全息测量声压与边界表面流体介质质点法向振速的数学关系，边界表面的质点法向振速可由下式重构计算：
[0111][0112]
其中，为全息测量声压到边界表面质点法向振速的传递矩阵：
[0113][0114]
其中，n为边界表面在重构点的单位法向矢量。
[0115]
步骤6，计算边界表面声阻抗率。
[0116]
边界表面的声阻抗率定义为边界表面声压p
half
(x；ω)和质点法向振速v
nhalf
(x；ω)的比值。利用重构的一组声压和法向振速计算各重构点处的表面声阻抗率，并取平均值，得到边界表面声阻抗率为：
[0117][0118]
实施例：以脉动球作为声源，对边界一侧的半空间声场进行数值仿真计算。在图1所示的全局坐标系中，边界坐标为z＝0，其表面声阻抗率满足声阻抗模型：
[0119]
z0＝0.436(1+i)(σe/f)
0.5
+19.48iαe/f
ꢀꢀꢀ
(19)
[0120]
其中，σe为边界有效流阻率，取38kpasm-2
；αe为边界孔隙度随边界深度的降低率，取15m-1
；f为声波频率。
[0121]
单层传声器阵列面与边界平行，其几何中心位于z轴，到边界距离为ha＝0.3m；阵列孔径为0.15m
×
0.15m，其上均布6
×
6个声压测点，相邻测点间距为0.03m，测点分布示意图参见图2。脉动球心坐标设为此时，声波入射角为θ
inc
＝0
°
。脉动球半径为a＝0.05m，其表面质点径向振动速度v0＝0.01m/s。空气介质密度为ρ0＝1.29kg/m3，空气中的声速为c＝343m/s。模拟传声器测量误差的影响，对测量值加入信噪比为30db的高斯白噪声。
[0122]
使用本发明方法，首先利用传声器阵列的声压测量值重构边界表面声压和质点法向振速分布，然后计算边界表面的声阻抗率。考察频率范围为f＝1000hz～5000hz。不同频率条件下，边界表面声阻抗率计算值与理论值的对比如图3所示。
[0123]
观察图3发现，当声波入射角为θ
inc
＝0
°
时，本发明方法计算得到的不同频率条件下的边界表面声阻抗率与其理论值相符。当调整声源位置，声波入射角取θ
inc
＝0
°
～90
°
范围内的其它值时，本发明方法也能够计算得到与理论值相符的声阻抗率。结果表明，本发明的单层传声器阵列原位测量半空间边界表面声阻抗的方法，通过在边界近场原位采集全息声压，能够精确获取边界表面声阻抗。
[0124]
本说明书实施例所述的内容仅仅是对发明构思的实现形式的案例之一。本发明的保护范围包含但不限于实施例所陈述的具体形式和参数，也包含本领域技术人员根据本发明构思所能够想到的同等技术手段。