FFT在平面机械结构声学灵敏度分析中的应用.pdfVIP

2008年12月 农 业 机 械 学 报 第39卷第12期 FFT在平面机械结构声学灵敏度分析中的应用* 张永斌 毕传兴 陈 剑 陈心昭 摘【要】 平面类或壳体类辐射体是组成复杂辐射结构的基础元件，对其进行声学灵敏度分析是进行复杂声学 结构优化设计的基础。将 F0urier变换引入平面辐射结构的声学灵敏度分析中，实现了声功率和空间任意点的声 学量关于声源的结构参数、材料特性等设计变量的灵敏度，又 由于FFT算法 比数值积分计算速度更快，从而改善 了基于Ryl。igh积分法的声学灵敏度分析，提高了其计算速度和工作效率。以简支铝板声源为例进行了验证，并 与Ray1eigh积分法计算结果进行比较，证 明了该方法的正确性和有效性。 关键词：声学灵敏度 平面辐射体 FFT 中图分类号 ：TB52 9 文献标识码 ：A v(z，)——声源面上的法向振速 引言 愚——波数 ——角频率 目前，声学灵敏度 的研究方法大都基于有限元 — — 媒质密度 c——声速 法和边界元法[卜 。有限元法 只适合用于 内声场 若用 “*”表示卷积，式(1)可简化为 问题的声学灵敏度计算，边界元法又要求声源必须 (z，， )=一iP (z，)*g(z，，) (2) 是闭合的，所以都不能用于平面声源的声学灵敏度 设设计变量为 ，则点(z，，)处的声压对设 分析。 目前 ，平面辐射体声学灵敏度分析 的主要方 计变量的灵敏度为 法是 Rayleigh积分法[ ]。但是该方法在应用时 需要进行数值积分，当声源较大或声源表面的离散 = 一 iP ————一 一fDL[——— *ggL(zz，，，，J)+十 点较多时，工作量大 ，计算时间长。 本文在Rayleigh积分 的基础上将 Fourier变换 1，(z…，))** 焉 ]I ((33)) 引入平面辐射体声学灵敏度计算过程，并以简支铝 定义空间内沿 z，方 向的Fourier变换为 板声源为例对此方法进行验证。 F(尼，忌)=j一。。j—o。_厂(z，)e一 ’dzd l 声学灵敏度分析 (4) 1．1 声压灵敏度计算 式中 是 尼——z、方 向的空间波数 在无限大障板中的平面声源，不考虑时间因素， 相应的逆变换为 其辐射声场可以用 Rayleh积分表示为 = 』』 ∥ 帆 (z，，)=一i II (z，)g(z—z，— ，)· 兰 (5) 对式(3)应用二维Fourier变换，得到 一 i (z ) (1) ———— 一 ：一 PF_1L[DL(走z，惫 ．’ 其中 R √(z—z)+(～ )+z 走= ／c 式中 g(z，，)——Green函数 G(忌， ，)+Vs(走，五)GD( ，愚，z)](6)