一种阻抗边界条件下封闭腔体内声压级的快速计算方法

摘要

本发明提供了一种阻抗边界下腔体内部声压级的计算方法。该计算方法将阻抗边界看作由点声源集合形成的面声源，基于刚性边界条件的腔体的声压分布和阻抗边界条件实现了腔体内声压级的快速计算。

说明书

技术领域

本发明涉及声学技术领域，具体涉及是一种阻抗边界下的腔体内部声压级的计算方法。

背景技术

随着汽车行业的发展，制造厂商开始关注车内的声学环境设计。对腔体内声压级分布的研究可用来指导汽车内的吸声材料布放和声源位置选择。

传统上，腔体声压级的计算通常使用傅里叶级数分解或腔板耦合的计算方法，前者计算精度较高，但运算量大且数学形式复杂，很难用于指导车内声学环境的优化和声学品质的提高；后者的形式相对简单，但仅适用于板振动的低频段，而无法适用于阻抗边界主导的中频段。

发明内容

本发明的目的在于提一种适用于阻抗边界且计算方便、易于实用的计算方法具有重要意义，本发明将阻抗边界看作由点声源集合形成的面声源，基于刚性边界条件的腔体的声压分布和阻抗边界条件实现了腔体内声压级的快速计算。

为实现上述目的，本发明提供如下技术方案：

本发明提供的一种阻抗边界下的腔体内部声压级的计算方法，具体包括如下步骤：

步骤1：直角坐标系下，刚性壁边界条件下的矩形腔体内任意一点对应的声压p(x)可以写成模态函数与其系数乘积的和：

其中，x＝(x

步骤2：刚性壁边界条件下n阶的模态函数如下式：

其中，L

n

对于某声源分布q

其中，ρ

A

ξ

步骤3：阻抗边界的表达如下：

其中，Z(y)表示法向声阻抗率，将阻抗边界的声压p(y)按照式(1) 展开：

其中，q

步骤4：

在所述腔体内声场分布，将阻抗边界对声场的影响可以看成声源，所述声源的振动速度和阻抗成反比；并且结合步骤3中的式(7) 令

用矩阵表达得出：

a＝[I+Z

其中，a＝[a

步骤5：由步骤1-4，得出，腔体内任一点x的声压级为：

p(x)＝a

其中，ψ＝[ψ

本发明还提供了一种计算机存储介质，其上存储有计算机程序，该程序被处理器执行时实现任意一步骤所述方法的步骤。

与现有技术相比，本发明具有以下几个方面的有益效果：

(1)简化阻抗边界下腔体计算速度快，通过程序一键完成。

(2)数学形式简单，物理意义明确，便于腔体内声品质的优化。

附图说明

图1为本发明的腔体几何示意图；

图2为本发明的计算结果与仿真结果比较图；

具体实施方式

下面结合具体实施方式对本发明的技术方案作进一步详细地说明。

实施例1

本实施例提供了一种阻抗边界下的腔体内部声压级的计算方法，包括如下步骤：

步骤1：直角坐标系下，刚性壁边界条件下的矩形腔体内任意一点对应的声压p(x)可以写成模态函数与其系数乘积的和：

其中，x＝(x

步骤2：刚性壁边界条件下n阶的模态函数如下式：

其中，L

n

对于某声源分布q

其中，ρ

A

ξ

步骤3：阻抗边界的表达如下：

其中，Z(y)表示法向声阻抗率，将阻抗边界的声压p(y)按照式(1) 展开：

其中，q

步骤4：

在所述腔体内声场分布，将阻抗边界对声场的影响可以看成声源，所述声源的振动速度和阻抗成反比；并且结合步骤3中的式(7) 令

用矩阵表达得出：

a＝[I+Z

其中，a＝[a

步骤5：由步骤1-4，得出，腔体内任一点x的声压级为：

p(x)＝a

其中，ψ＝[ψ

实施例2

下面以计算腔体内一点处的1～500Hz频段声压级为例，结合实施例1的具体运算步骤，对本发明的实施做详细说明。

腔体几何尺寸为3.4m×1.6m×1.2m，接收点坐标为(1m,1m,1m)。如图1所示，腔体内共振角频率为：

计算频率上限为500Hz，系统需满足

声阻尼ξ

其中δ为狄拉克函数，x

对比解析结果与仿真结果，可以看出本方法具有较好的准确度。

如图2所示，本发明计算结果与仿真结果比较下面以该方法计算腔体内一点处的1～500Hz频段声压级为例，对本发明的实施做详细说明。