基于逆边界元法识别机械表面振动强度的方法
【技术领域】
[0001] 本发明涉及的是一种振动强度预测方法。
【背景技术】
[0002] 目前，基于空间Fourier变换和逆变换法的近场声全息的声源识别有实现过程简 单，应用方便特点，但由于空间Fourier变换和逆变换法只适用于较简单结构，对任意形状 的振动结构无法实施。
[0003] 统计能量分析法将系统分解为多个子系统，运用能量流关系式对复合的、谐振的 组装结构进行动力特性、振动响应级及声辐射的理论评估，然而它的分析频率范围限制在 高频，对低频效果较差。
【发明内容】
[0004] 本发明的目的在于提供适用于任意形状、中低频振动机械结构表面振动强度预测 的基于逆边界元法识别机械表面振动强度的方法。
[0005] 本发明的目的是这样实现的：
[0006] 本发明基于逆边界元法识别机械表面振动强度的方法，其特征是：
[0007] (1)在振动声源福射声场布置参考传声器与传声器阵列来测量福射声场的场点声 压，二者至振动声源的距离小于分析最高频率对应的波长，参考传声器个数1-2个，传声器 阵列的测量场点个数满足一个声波波长内含有3-5个测点，测量面积大于振动声源正投影 面；
[0008] (2)参考传声器与传声器阵列测量得到的声压利用互谱法得到传声器阵列对应的 测量场点处频域上的复声压，该复声压作为求解机械表面振动强度输入的已知参数；
[0009] (3)基于逆边界元法计算振动表面法向振动速度：
[0010] 建立动表面与测量点积分方程为
[0011]
【主权项】
1.基于逆边界元法识别机械表面振动强度的方法，其特征是： (1) 在振动声源福射声场布置参考传声器与传声器阵列来测量福射声场的场点声压， 二者至振动声源的距离小于分析最高频率对应的波长，参考传声器个数1-2个，传声器阵 列的测量场点个数满足一个声波波长内含有3-5个测点，测量面积大于振动声源正投影 面； (2) 参考传声器与传声器阵列测量得到的声压利用互谱法得到传声器阵列对应的测量 场点处频域上的复声压，该复声压作为求解机械表面振动强度输入的已知参数； (3) 基于逆边界元法计算振动表面法向振动速度： 建立动表面与测量点积分方程为
其中：c (P)为场点复声压系数，c (P)括号内P表不福射声场场点，p (P)为场点复声压， P (Q)、v (Q)"分别为振动声源边界复声压和边界法向振动速度，Q为振动声源边界点，基本解
为基本解的法向导数，k为场点声压系数，r为辐射声场场点P与振动声源边 界点Q的距离，
，P为空气密度，ω = 2 π f，振动声源边界为S，C(P)可表示为
其中：V振动声源辐射声场，将振动声源边界s离散
其中η为振动声源边界S离散后单元数，|2为高斯积分点，Nk是ξ i、|2的函数， 角标k为单元节点数，1为高斯积分取的点数，J为雅克比； 由此得到场点复声压矩阵Pf与振动声源边界复声压矩阵P s及边界法向振动速度矩阵 Vn关系： CPf= HP s+GVn, 其中矩阵C、H、G为系数矩阵； 为得到振动声源边界复声压P (Q)及边界法向振动速度V (Q) n关系，将场P点移到振动 Ρψ 声源边界S上，由于ψ与ι出现奇异，采用广义极坐标变换来消除奇异；设S离散为三角 an 形单元，将移到边界单元上的P点作为三角形的顶点，利用极坐标变换
? 其中上标1、2、3表示以P为顶点的三角形按逆时针各点编号，P所在的点上标为1，下 标1、2三角形各点的两个坐标，此时有
其中 此时1、η2为高斯积分点；由以上关系得到声场复 , ， 声压与边界法向振动速度关系 Ps=ATMs-Vn, 其中ATMs为振动声源边界复声压与边界法向振动速度声学的传递矩阵，由式CP f= HPs+GVn以及P s= ATM s · Vn得到场点复声压与振动边界法向振动速度关系 Pf= ATMf · Vn, 其中为Pf场点复声压矩阵，ATMf为场点复声压与边界法向振动速度声学传递矩阵，由 此得边界法向振动速度 Vn = ATM丨· Pi 其中/Γ/'Μ,为传递矩阵的逆矩阵，根据式已=嫌3乂，利用场点复声压计算得到表面法向振 动速度;根据等式
得到的振动声源法向振动速度与边界声压关系计算得到边界声压Ps; (4)根据表面法向振动速度计算振动体辐射声功率 W = / s0. 5re (p (Q) sv (Q) n) ds, 其中re表示边界声压与边界法相振动速度实部，根据振动机械辐射声功率评判振动 强度大小。
【专利摘要】本发明的目的在于提供基于逆边界元法识别机械表面振动强度的方法，首先建立振动机械边界声压和边界法向振动速度与场点声压的传递矩阵，在求解边界声压与边界法向振动速度关系矩阵时出现的奇异积分采用广义极坐标变换来消除。利用传声器阵列测得的机械振动辐射声场复声压作为已知，根据机械振动边界声压、边界法向振动速度与场点声压的传递矩阵，重建边界法相振动速度，进而获得边界声压，利用边界声压与边界法向振动速度计算机械振动辐射声功率，实现机械表面振动强度的预测。本发明利用机械辐射声场复声压作为计算的已知量，避免了直接接触测量表面振动速度。振动机械形状没有限制，可以对任意形状机械振动强度计算。
【IPC分类】G01H1-00
【公开号】CN104764518
【申请号】CN201510146310
【发明人】肖友洪, 吕国学, 田新娜, 郑伟, 褚阵豪
【申请人】哈尔滨工程大学
【公开日】2015年7月8日
【申请日】2015年3月31日