基于逆边界元法机械噪声远场声压预测方法技术

本发明专利技术公开了基于逆边界元法机械噪声远场声压预测方法，属于声压预测领域。在目标声源最高分析频率的一个波长内布置传声器阵列，测试面大于目标声源正投影面，一个波长内至少含有2个测量点，在目标声源附近放置参考传声器，测量得到传声器阵列与参考传声器互谱后的场点复声压；基于逆边界元法建立场点复声压与法向振动速度的传递关系，得到传递矩阵；对传递矩阵进行奇异值分解，得到法向振动速度；根据法向振动速度预测远场声压py，根据边界积分方程建立远场声压与法向振动速度的关系py＝ATMyvn，ATMy为对应远场声压的传递矩阵。本发明专利技术能够适用于复杂结构表面的，具有高精度的优点。

【技术实现步骤摘要】

本专利技术属于声压预测领域，尤其涉及一种适用于结构复杂振动体的远场声压预测 的，。
技术介绍
目前，基于近场声全息技术的表面振动速度重建方法是利用靠近声源或振动表面 (d<<A)处测量声源的全息数据，其全息数据包含了随距离指数递减的倐逝波成分，从 而可以获得较高的分辨率，然而其要求测量距离d<<A过于苛刻，限制其应用。基于远 场声全息的测量方法是基于平面波假设，经过特殊传感器测量远离声源的声压，通过计算 传声器接收到的信号相位差重建表面振速。但由于远场声压信号记录较少的倐逝波，因此 分辨率受到波长的限制。 直接测量振动体表面法向振动速度计算辐射声压技术是基于振动体辐射的声功 率与表面振动之间的关系利用振动信号去直接预估辐射噪声，表面法向振动速度和表面积 等一定，声功率估算的关键是如何确定出辐射系数，由于辐射系数不仅与部件的形状和边 界条件有关，还与振动频率相关，因此高精度计算十分困难。并且对于复杂的结构表面无法 实施。
技术实现思路
本专利技术的目的提供一种适用于复杂结构表面的，高精度的，基于逆边界元法机械 噪声远场声压预测方法。 ，包括以下几个步骤， 步骤一：在目标声源最高分析频率的一个波长内布置传声器阵列，测试面大于目 标声源正投影面，一个波长内至少含有2个测量点，在目标声源附近放置参考传声器，测量 得到传声器阵列与参考传声器互谱后的场点复声压； 步骤二：基于逆边界元法建立场点复声压与法向振动速度的传递关系，得到传递 矩阵； 步骤三：对传递矩阵进行奇异值分解，得到法向振动速度； 步骤四：根据法向振动速度预测远场声压py，根据边界积分方程建立远场声压与 法向振动速度的关系py=ATMyvn，ATMyS对应远场声压的传递矩阵。 本专利技术，还可以包括： 1、场点复声压与法向振动速度的传递关系为： pf=ATMvn 其中，pf场点声压，vn为法向振动速度，ATM为传递矩阵。 2、法向振动速度为：CN 105004416 A W ~P 2/4页 其cUp\为单位正交矩阵U、V中的列向量，人为正则化参 数，〇1> 〇 2>... > 〇nf>〇为传递矩阵的奇异值，对传递矩阵奇异值分解： t=U2Vt 其中U、V为单位正交矩阵，即UUT=I，VVT=I， 正则化参数A利用广义交叉检验法选取： 其中，+为ATM的广义伪逆矩阵，xA为正则化后得到的正则解，traceO表示 矩阵的迹，1"单位矩阵，当g(A)取得最小值对应的A作为正则化参数。 有益效果： 一种基于逆边界元法的机械噪声远场声压预测方法，利用近场测得的声压重建声 源表面振动法向速度，可以识别声源的振动强度。 利用重建的表面法向振动速度预测远场声压大小，避免了多次测量的不便，并且 可以利用重建的表面振动速度预测传声器不易安装点处的声压，给工程测量带来方便，并 且适用于复杂形状振动体远场声压预测。 本专利技术目的在于克服现有的技术的不足，提供一种基于逆边界元法的机械噪声远 场声压预测方法，传递矩阵的广义伪逆矩阵的小奇异值引起的高误差敏感性采用截断奇 异值正则化进行抑制，其中正则参数利用广义交叉检验发选取，近场场点声压作为已知量， 重建振源表面振动速度，利用得到的表面振动速度预测远场声压，该方法适用于结构复杂 振动体的远场声压预测。【附图说明】 图1为本专利技术操作流程图。【具体实施方式】 下面将结合附图对本专利技术做进一步详细说明。 本专利技术为解决其技术问题所采用的技术方案，如图1所示： 步骤1 :建立声源、测试面模型的仿真模型。振动声源模型划分网格单元要求小于 最大分析频率对应的波长含有1/6。 步骤2 :建立场点声压与机械表面振动速度传递关系矩阵。基于边界元法建立边 界积分方程其中：p⑵为场点P声压，p(Q)、v(Q)j别为声源表面声压和法向振动速度。为 建立表面声压与表面法向振动速度关系，将场点P移到声源表面，此时有 L1N丄USUU44丄ti fl p/J }J\⑵ 根据边界积分方程（1)、（2)得到场点声压与边界表面振动传递矩阵Pf=ATM*Vn， 其中为Pf测试场点声压矩阵，Vn为声源表面法向振动速度矩阵，ATM为传递矩阵。 步骤3 :对传递矩阵采取奇异值分解技术 t=U2Vt (3) 其中U、V为单位正交矩阵，即UUT= I，VVT=I，奇异值矩阵2是非负、递减对角 阵，8卩〇1> 〇 2>... > 〇nf>〇。由于5：中的小奇异值导致传递矩阵对误差敏感性提 高，采用截断奇异值正则化将对解的贡献量较小的奇异值赋给零。 其^4、\为单位正交矩阵U、V中的列向量，人为正则化参 n， 数，利用广义交叉检验法选取(5) Vn=+Pf，+为ATM的广义伪逆矩阵，+=U2+VT， 其中+为ATM矩阵正则后的广义伪逆矩阵，xx为正则化后得到的正则解， trace()表示矩阵的迹，In单位矩阵。g(A)是关于A的函数。当g(A)取得最小值对应 的A作为正则化参数。利用声阵列传感器与参考传声器互谱后的场点复声压以重建表面 振动速度。 Vn= +Pf (6) 其中+为矩阵正则后的广义伪逆矩阵，Vn为声源表面法向振动速度。 步骤4 :根据表面振动速度预测远场声压py，根据边界积分方程建立远场声压与边 界表面振动速度的关系py=ATMyvn，利用重建的表面振动计算远场声压。 本专利技术涉及一种基于逆边界元法的机械噪声远场声压预测方法，传统声压预测 是根据结构表面振动速度，利用边界元法预测远场声压， 而基于逆边界元法声压预测是根据近场声压预测远场声压，根据逆边界元法建立 场点声压与表面振动速度的传递关系，通过传递矩阵的逆变换重建表面振动速度，在利用 边界元法计算场点声压。机械表面振动速度的反演过程由于不适定问题导致反演结果波 动，这里首先采用奇异值分解求广义逆矩阵，利用截断奇异值方正则化进行遏制结果波动， 其中截断参数选用广义交叉检验法选取，进而根据边界元法预测远场声压。基于逆边界元 法的机械噪声远场声压预测方法，也可以为： 第一步：基于逆边界元法建立场点声压与声源表面法向振动速度关系，声传递矩 阵的采用截断奇异值处理。 第二步：测量振动体场点声压，计算振动机械表面振动速度，根据边界元法计算远 场声压。根据边界元法利用表面振动速度重建远场声压P，根据边界积分方程（1)建立远场 声压与边界表面振动速度的关系py=ATMyvn，由振动速度作为已知量获得远场声压。 结合实施例一对本专利技术做进一步详细说明。 根据声源与测量场点位置建立声源、测试场点仿真模型。根据逆边界元法只需要 划分声源表面网格，要求振动声源模型面网格单元小于最大分析频率对应的波长含有1/6。 编写场点声压与声源表面法向振动速度关系矩阵计算程序。程序的理论基础是基 于逆边界元法建立边界积分方程 为减少未知数，建立边界声压与边界法向振动速度将场点P移到边界表面，此时 有M " \ Uil J根据边界积分方程得到场点声压与边界表面振动关系Pf=ATMn，ATM为声传递 矩阵。 步骤4 :利用程序对传递矩阵采取奇异值分解技术求解广义逆矩阵，并采取正则 化处理，求解表面振动速度，其中传递矩阵 t=U2Vt (9) 小奇异值导致传递矩阵对误差敏感本文档来自技高网...

【技术保护点】
基于逆边界元法机械噪声远场声压预测方法，其特征在于：包括以下几个步骤，步骤一：在目标声源最高分析频率的一个波长内布置传声器阵列，测试面大于目标声源正投影面，一个波长内至少含有2个测量点，在目标声源附近放置参考传声器，测量得到传声器阵列与参考传声器互谱后的场点复声压；步骤二：基于逆边界元法建立场点复声压与法向振动速度的传递关系，得到传递矩阵；步骤三：对传递矩阵进行奇异值分解，得到法向振动速度；步骤四：根据法向振动速度预测远场声压py，根据边界积分方程建立远场声压与法向振动速度的关系py＝ATMyvn，ATMy为对应远场声压的传递矩阵。

【技术特征摘要】

【专利技术属性】
技术研发人员：肖友洪，吕国学，季振林，田新娜，郑伟，褚阵豪，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江;23