基于逆边界元法机械噪声远场声压预测方法
【技术领域】
[0001] 本发明属于声压预测领域，尤其涉及一种适用于结构复杂振动体的远场声压预测 的，基于逆边界元法机械噪声远场声压预测方法。
【背景技术】
[0002] 目前，基于近场声全息技术的表面振动速度重建方法是利用靠近声源或振动表面 (d<<A)处测量声源的全息数据，其全息数据包含了随距离指数递减的倐逝波成分，从 而可以获得较高的分辨率，然而其要求测量距离d<<A过于苛刻，限制其应用。基于远 场声全息的测量方法是基于平面波假设，经过特殊传感器测量远离声源的声压，通过计算 传声器接收到的信号相位差重建表面振速。但由于远场声压信号记录较少的倐逝波，因此 分辨率受到波长的限制。
[0003] 直接测量振动体表面法向振动速度计算辐射声压技术是基于振动体辐射的声功 率与表面振动之间的关系利用振动信号去直接预估辐射噪声，表面法向振动速度和表面积 等一定，声功率估算的关键是如何确定出辐射系数，由于辐射系数不仅与部件的形状和边 界条件有关，还与振动频率相关，因此高精度计算十分困难。并且对于复杂的结构表面无法 实施。
【发明内容】
[0004] 本发明的目的提供一种适用于复杂结构表面的，高精度的，基于逆边界元法机械 噪声远场声压预测方法。
[0005] 基于逆边界元法机械噪声远场声压预测方法，包括以下几个步骤，
[0006] 步骤一：在目标声源最高分析频率的一个波长内布置传声器阵列，测试面大于目 标声源正投影面，一个波长内至少含有2个测量点，在目标声源附近放置参考传声器，测量 得到传声器阵列与参考传声器互谱后的场点复声压；
[0007] 步骤二：基于逆边界元法建立场点复声压与法向振动速度的传递关系，得到传递 矩阵；
[0008] 步骤三：对传递矩阵进行奇异值分解，得到法向振动速度；
[0009] 步骤四：根据法向振动速度预测远场声压py，根据边界积分方程建立远场声压与 法向振动速度的关系py=ATMyvn，ATMyS对应远场声压的传递矩阵。
[0010] 本发明基于逆边界元法机械噪声远场声压预测方法，还可以包括：
[0011] 1、场点复声压与法向振动速度的传递关系为：
[0012] pf=ATMvn
[0013] 其中，pf场点声压，vn为法向振动速度，ATM为传递矩阵。
[0014] 2、法向振动速度为：
[0015]
CN 105004416 A W ~P 2/4页
[0016] 其c
Up\为单位正交矩阵U、V中的列向量，人为正则化参 数，〇1> 〇 2>... > 〇nf>〇为传递矩阵的奇异值，对传递矩阵奇异值分解：
[0017] [ATM]t=U2Vt
[0018] 其中U、V为单位正交矩阵，即UUT=I，VVT=I，
[0019] 正则化参数A利用广义交叉检验法选取：
[0020]
[0021] 其中，[ATM]+为ATM的广义伪逆矩阵，xA为正则化后得到的正则解，traceO表示 矩阵的迹，1"单位矩阵，当g(A)取得最小值对应的A作为正则化参数。
[0022] 有益效果：
[0023] 一种基于逆边界元法的机械噪声远场声压预测方法，利用近场测得的声压重建声 源表面振动法向速度，可以识别声源的振动强度。
[0024] 利用重建的表面法向振动速度预测远场声压大小，避免了多次测量的不便，并且 可以利用重建的表面振动速度预测传声器不易安装点处的声压，给工程测量带来方便，并 且适用于复杂形状振动体远场声压预测。
[0025] 本发明目的在于克服现有的技术的不足，提供一种基于逆边界元法的机械噪声远 场声压预测方法，传递矩阵的广义伪逆矩阵的小奇异值引起的高误差敏感性采用截断奇 异值正则化进行抑制，其中正则参数利用广义交叉检验发选取，近场场点声压作为已知量， 重建振源表面振动速度，利用得到的表面振动速度预测远场声压，该方法适用于结构复杂 振动体的远场声压预测。
【附图说明】
[0026] 图1为本发明操作流程图。
【具体实施方式】
[0027] 下面将结合附图对本发明做进一步详细说明。
[0028] 本发明为解决其技术问题所采用的技术方案，如图1所示：
[0029] 步骤1 :建立声源、测试面模型的仿真模型。振动声源模型划分网格单元要求小于 最大分析频率对应的波长含有1/6。
[0030] 步骤2 :建立场点声压与机械表面振动速度传递关系矩阵。基于边界元法建立边 界积分方程
[0031]
[0032]其中：p⑵为场点P声压，p(Q)、v(Q)j别为声源表面声压和法向振动速度。为 建立表面声压与表面法向振动速度关系，将场点P移到声源表面，此时有 L1N丄USUU44丄ti fl p/J }J\
[0033]
⑵
[0034] 根据边界积分方程（1)、（2)得到场点声压与边界表面振动传递矩阵Pf=ATM*Vn， 其中为Pf测试场点声压矩阵，Vn为声源表面法向振动速度矩阵，ATM为传递矩阵。
[0035] 步骤3 :对传递矩阵采取奇异值分解技术
[0036] [ATM]t=U2Vt (3)
[0037] 其中U、V为单位正交矩阵，即UUT= I，VVT=I，奇异值矩阵2是非负、递减对角 阵，8卩〇1> 〇 2>... > 〇nf>〇。由于5：中的小奇异值导致传递矩阵对误差敏感性提 高，采用截断奇异值正则化将对解的贡献量较小的奇异值赋给零。
[0038]
[0039] 其^
4、\为单位正交矩阵U、V中的列向量，人为正则化参 n， 数，利用广义交叉检验法选取
[0040]
(5)
[0041] Vn=[ATM]+Pf，[ATM]+为ATM的广义伪逆矩阵，[ATM]+=U2+VT，
[0042]
[0043] 其中[ATM]+为ATM矩阵正则后的广义伪逆矩阵，xx为正则化后得到的正则解， trace()表示矩阵的迹，In单位矩阵。g(A)是关于A的函数。当g(A)取得最小值对应 的A作为正则化参数。利用声阵列传感器与参考传声器互谱后的场点复声压以重建表面 振动速度。
[0044] Vn= [ATM] +Pf (6)
[0045] 其中[ATM]+为矩阵正则后的广义伪逆矩阵，Vn为声源表面法向振动速度。
[0046] 步骤4 :根据表面振动速度预测远场声压py，根据边界积分方程建立远场声压与边 界表面振动速度的关系py=ATMyvn，利用重建的表面振动计算远场声压。
[0047] 本发明涉及一种基于逆边界元法的机械噪声远场声压预测方法，传统声压预测 是根据结构表面振动速度，利用边界元法预测远场声压，
[0048] 而基于逆边界元法声压预测是根据近场声压预测远场声压，根据逆边界元法建立 场点声压与表面振动速度的传递关系，通过传递矩阵的逆变换重建表面振动速度，在利用 边界元法计算场点声压。机械表面振动速度的反演过程由于不适定问题导致反演结果波 动，这里首先采用奇异值分解求广义逆矩阵，利用截断奇异值方正则化进行遏制结果波动， 其中截断参数选用广义交叉检验法选取，进而根据边界元法预测远场声压。基于逆边界元 法的机械噪声远场声压预测方法，也可以为：
[0049] 第一步：基于逆边界元法建立场点声压与声源表面法向振动速度关系，声传递矩 阵的采用截断奇异值处理。
[0050] 第二步：测量振动体场点声压，计算振动机械表面振动速度，根据边界元法计算远 场声压。根据边界元法利用表面振动速度重建远场声压P，根据边界积分方程（1)建立远场 声压与边界表面振动速度的关系py=ATMyvn，由振动速度作为已知量获得远场声压。
[0051] 结合实施例一对本发明做进一步详细说明。
[0052] 根据声源与测量场点位置建立声源、测试场点仿真模型。根据逆边界元法只需要 划分声源表面网格，要求振动声源模型面网格单元小于最大分析频率对应的波长含有1/6。
[0053] 编写场点声压与声源表面法向振动速度关系矩阵计算程序。程序的理论基础是基 于逆边界元法建立边界积分方程
[0054]
[0055] 为减少未知数，建立边界声压与边界法向振动速度将场点P移到边界表面，此时 有
[0056]
M " \ Uil J
[0057]根据边界积分方程得到场点声压与边界表面振动关系Pf=ATMn，ATM为声传递 矩阵。
[0058] 步骤4 :利用程序对传递矩阵采取奇异值分解技术求解广义逆矩阵，并采取正则 化处理，求解表面振动速度，其中传递矩阵
[0059] [ATM]t=U2Vt (9)
[0060] 小奇异值导致传递矩阵对误差敏感性提高，采用截断奇异值正则化将对解的贡献 量较小的奇异值赋零值
[0061] /=1
[0062] 其e
;中A利用广义交叉检验法选取
[0063]
[0064] 当g(A)取得最小值对应的A作为正则化参数。利用声阵列传感器与参考传声 器互谱后的场点复声压以重建表面振动速度。
[0065] Vn= [ATM] +Pf (12)
[0066] 步骤5 :根据表面振动速度预测远场声压py，根据边界积分方程（7)建立远场声压 与边界表面振动速度的关系Py=ATMy ?Vn，利用重建的表面法向振动速度作为已知量获得 远场声压。
【主权项】
1. 基于逆边界元法机械噪声远场声压预测方法，其特征在于：包括w下几个步骤， 步骤一：在目标声源最高分析频率的一个波长内布置传声器阵列，测试面大于目标声 源正投影面，一个波长内至少含有2个测量点，在目标声源附近放置参考传声器，测量得到 传声器阵列与参考传声器互谱后的场点复声压； 步骤二：基于逆边界元法建立场点复声压与法向振动速度的传递关系，得到传递矩 阵； 步骤=:对传递矩阵进行奇异值分解，得到法向振动速度； 步骤四：根据法向振动速度预测远场声压Py，根据边界积分方程建立远场声压与法向 振动速度的关系Py=ATMyV。，ATMy为对应远场声压的传递矩阵。2. 根据权利要求1所述的基于逆边界元法机械噪声远场声压预测方法，其特征在于： 所述的场点复声压与法向振动速度的传递关系为： Pf=ATMvn 其中，Pf场点声压，V。为法向振动速度，ATM为传递矩阵。3. 根据权利要求1所述的基于逆边界元法机械噪声远场声压预测方法，其特征在于： 所述的法向振动速度为：其中Ui、Vi为单位正交矩阵U、V中的列向量，A为正则化参数，0 1 > 〇2>w> 0M>0为传递矩阵的奇异值，对传递矩阵奇异值分解： [atm]t=usyT 其中u、V为单位正交矩阵，即UIJT=I，wT二I， 正则化参数A利用广义交叉检验法选取：其中，[ATM]+为ATM的广义伪逆矩阵，X,为正则化后得到的正则解，traceO表示矩阵 的迹，I。单位矩阵，当g(A)取得最小值对应的A作为正则化参数。
【专利摘要】本发明公开了基于逆边界元法机械噪声远场声压预测方法，属于声压预测领域。在目标声源最高分析频率的一个波长内布置传声器阵列，测试面大于目标声源正投影面，一个波长内至少含有2个测量点，在目标声源附近放置参考传声器，测量得到传声器阵列与参考传声器互谱后的场点复声压；基于逆边界元法建立场点复声压与法向振动速度的传递关系，得到传递矩阵；对传递矩阵进行奇异值分解，得到法向振动速度；根据法向振动速度预测远场声压py，根据边界积分方程建立远场声压与法向振动速度的关系py＝ATMyvn，ATMy为对应远场声压的传递矩阵。本发明能够适用于复杂结构表面的，具有高精度的优点。
【IPC分类】G06F19/00, G01H17/00
【公开号】CN105004416
【申请号】CN201510381902
【发明人】肖友洪, 吕国学, 季振林, 田新娜, 郑伟, 褚阵豪
【申请人】哈尔滨工程大学
【公开日】2015年10月28日
【申请日】2015年7月2日