基于快速多极边界元法声学敏感度分析在声屏障结构形状优化中应用障结构形状优化中应用.pdfVIP

基于快速多极边界元法的声学敏感度分析在声屏 基于快速多极边界元法的声学敏感度分析在声屏 基基于于快快速速多多极极边边界界元元法法的的声声学学敏敏感感度度分分析析在在声声屏屏 障结构形状优化中的应用 障结构形状优化中的应用 障障结结构构形形状状优优化化中中的的应应用用 郑昌军，陈海波，陈磊磊 （中国科学技术大学，合肥 230027） 摘要 摘要 摘摘要要：现今，交通噪声已成为城市噪声的主要来源，而声屏障作为纾缓交通噪声的重要工具已逐渐得到了广泛 的应用。本文基于快速多极边界元法的三维声学敏感度分析，开展声屏障的三维结构形状优化研究，通过对声 屏障顶端形状控制参数的优化从而得到一组降噪性能更优的声屏障。声学敏感度分析采用直接微分法，并设计 出相应的结构形状优化算法。为了克服边界元法在分析求解外声场问题时所遇到的解的非唯一性问题，本文采 用了Burton-Miller 方法，并通过常量单元离散来避免使用正则化等消除超奇异性的方法，从而提高了快速多极 边界元法的分析效率。 关键词 关键词 关关键键词词：边界元法；快速多极算法；声学敏感度；声屏障；结构形状优化 引 言 引 言 引引 言言 近年来，我国的城市交通噪声污染已日益严重，特别是由于一些市内主干公路、城际高速公路和轨 道交通难免会穿越人口密集区，所带来的噪声污染更是给周边居民带来了巨大的困扰。为了降低交通噪 声的影响，设置声屏障成为目前最为经济和有效的途径。虽然声屏障具有很好的降噪效果，但其性能取 决于其材料和尺寸形状等因素，不同的设计对其降噪效果有着很大的影响。根据国内外的相关研究经验， [1-5] 目前主要通过改进声屏障的材料吸声性能和几何形状来达到不增加高度而提高降噪性能的目的 。例 如，Okubo和Fujiwara研究了顶端为水车轮式的声屏障，探讨不同了的槽深、开口角度和直径对其降噪 [3] 效果的影响 ；Ishizuka 和Fujiwara基于前人的研究成果，对顶端具有不同形状和材料的声屏障的降噪效 [4] 果进行了系统的对比研究 。虽然对提高声屏障的降噪性能有着很好的指导意义，但是这些研究采用了 类似半经验的研式究方式，即根据经验或者前人的研究成果预先设计出数种不同形状或材料的声屏障， 然后通过数值计算得出的插入损失而决定孰优孰劣。若要进一步提高声屏障的降噪性能，半经验式的研 究方式往往是不够的，还需要对不同形状的声屏障的形状参数进行优化。而在目前常用的各种优化算法 [6] 中，基于梯度的优化方法较之其他方法通常具有更高的求解效率 。 在基于梯度的优化算法中，敏感度分析通常是重要而又关键的一步。声学形状敏感度分析给出了辐 射声场中的声学量因结构形状改变而引起的变化率，在形状优化过程中，这些变化率可以用来决定最优 搜索方向。目前国内外很多学者和研究人员利用有限元和边界元等数值方法提出了很多敏感度分析方法， 但由于边界元法更适于处理无限域或半无限域问题，加之具有降维和精度高等优点，令其在分析和评估 1 [1-5] 声屏障降噪性能中得到了广泛应用 。然而传统边界元法由于所形成的系统方程的系数矩阵是非对称满 2 阵，所以具有高内存占有量和高计算量的固有缺点。对于一个N自由度的问题，保存这个满阵需要O(N) 3 量级的内存空间。若采用高斯消元法之类的直接法求解这个问题需要 O(N)量级的计算量；而若使用