一种计算周期结构板声学散射系数的方法

摘要

本发明提出了一种计算周期结构板声学散射系数的方法，首先推导适合于周期结构板节点声压计算的系统方程；另一方面，将周期结构板及相同尺寸的参考板用节点进行建模，然后利用移动最小二乘法构建形函数，结合两者获得节点处声压差后，即可求得接收点声压，进而求得方向散射系数和平均散射系数。本发明将无网格法引入到周期结构的声学散射系数数值计算之中，避免了传统数值方法中因为网格的存在而导致的一系列问题。本发明有良好的自适应性，在需要提高计算频率上限时，可局部地增加节点密度，而不需对模型进行重新划分。通过与测量实验对比，验证了本发明具有较高的精度，因而在周期结构散射系数数值计算中具有广阔应用前景。

主权项

1.一种计算周期结构板声学散射系数的方法，其特征在于：包括以下步骤：步骤1：将待计算的周期结构板划分为包含n个均匀分布节点的周期结构板节点模型；建立一个纯平的参考板，参考板与周期结构板的平面投影形状相同，将所述参考板也划分为包含n个均匀分布节点的参考板节点模型；步骤2：在以周期结构板节点模型几何中心为球心的一个半球面上均匀设置各不少于100个声源点及接收点，所述半球面变径不小于周期结构板几何中心到周期结构板边缘最大距离的两倍；也在以参考板节点模型几何中心为球心的一个半球面上设置相同位置和数量的声源点及接收点；步骤3：利用移动最小二乘法，建立周期结构板节点模型和参考板节点模型的形函数；步骤4：分别计算周期结构板节点模型和参考板节点模型对应的系统方程，得到周期结构板节点模型及参考板节点模型上节点处的声压差p_d；所述系统方程通过离散Helmholtz微分方程得到，系统方程的形式为：(C+D)·p_d＝F式中C、D为n×n阶的系数矩阵，p_d为为n×1阶向量，指周期结构板节点模型或参考板节点模型上所有节点的声压差、F为n×1阶向量的载荷矩阵；其中c(P_i)＝α/4π，i＝1,...,nα为节点P_i处表面的立体角；m＝1,…,ni＝1,…,nG(P_i,Q)和G(P_i,r₀)均为格林函数，具体表示为表示P_i与声场中 任意一点Q之间的距离，且Q点与声源点和接受点不重合；为周期结构板节点模型或参考板节点模型的形函数中对应P_i及Q的元素，n_Q表示周期结构板节点模型或参考板节点模型上经过点Q的法向向量，表示周期结构板节点模型或参考板节点模型上经过点P_i的法向向量；Γ表示周期结构板节点模型或参考板节点模型的表面边界；ρ为空气密度，ω为圆频率，k为波数，j为虚数单位；r₀表示声源点，q_ω(r₀)表示声源点r₀的声源强度；步骤5：利用下式和步骤4得到的周期结构板节点模型及参考板节点模型上节点处的声压差p_d，分别计算周期结构板节点模型和参考板节点模型对应的接收点的声压：p_ω(R)＝-B·p_d其中p_ω(R)表示接受点R处的声压，N_RQ为周期结构板节点模型或参考板节点模型的形函数中对应R及Q的元素，G(R,Q)为格林函数；步骤6：根据下式计算每个声源点的方向散射系数：式中：θ和分别代表声源点相对于周期结构板节点模型几何中心的俯仰角与方位角；s为周期结构板节点模型中接收点的数量；θ_l和分别代表第l个接收点的俯仰角与方位角；是周期结构板节点模型所对应的第l个接收点的声压；是参考板节点模型所对应的第l个接收点的声压；*代表复共轭；步骤7：得到所有声源点的方向散射系数之后，按照下式计算平均散射系数δ_r:。