基于自适应交叉近似边界元法的快速声学灵敏度分析.pdfVIP

第31卷第4期 计 算 力 学 学 报 Vol_31，No．4 2O14年 8月 ChineseJournalofComputationalMechanics August2014 文章编号 ：1007—4708(2014)04—0413—06 基于 自适应交叉近似边界元法的 快速声学灵敏度分析 郑 昌军 ， 陈磊磊 ， 陈海波 (中国科学技术大学 近代力学系，合肥 230027) 摘 要：基于 自适应交叉近似边界元法构造一组快速声学灵敏度分析方法，其中，声学灵敏度分析分别采用直接 微分法和伴随变量法；而 自适应交叉近似算法被用以克服常规边界元法的高计算量和高存储量的固有缺点。自 适应交叉近似算法在迭代求解之前对边界元 系数矩阵进行压缩存储，可 以在 降低存储量 的同时提高求解效率 。 在声学灵敏度分析 中，通过直接使用求解未知边界状态值时保存 的压缩系数矩阵，可 以进一步提高求解效率。数 值算例验证 了所构造的方法 的计算精度和求解效率，以及在大规模声场 问题 的最优化分析 中的应用潜力。 关键词 ：边界元法；自适应交叉近似 ；声学灵敏度 ；直接微分法；伴 随变量法 中图分类号 ：0242．2 文献标志码 ：A doi：10．7511／js1x2O14O4OOl 量内积分解。将 ACA算法用于边界元子矩阵的 1 引 言 快速近似与存储 ，可以大幅提高求解速度并降低存 边界元法 BEM (BoundaryElementMethod) 储量，并且该算法是一种纯数值技术的矩阵压缩方 较之有限元法具有一些独特优点，如可对模型进行 法，与物理背景无关 ，因此与 FMA相 比更 易于应 降维处理，求解精度更高，可 以方便地处理无限域 用。另外，由于可以将压缩分解形成的矩阵存储于 或半无 限域问题等。但 由于边界元法所形成 的系 内存中，便于再次迭代求解 ，因而更适于需要多次 统方程的系数矩阵是非对称稠密矩阵，故而导致 了 迭代步骤的求解过程。 常规边界元法 的高计算量和高存储量的缺点。对 在基于梯度的最优化过程 中，灵敏度的计算最 于一个 N 自由度 的问题，存储量级为 0(N )，若 复杂且最耗时 ，也是决定最优化成败的关键部分。 采用直接方法求解需O(N。)量级 的计算量；而若 因此 ，设计高效且精确的灵敏度分析算法成为结构 采用迭代求解算法，则仍需要 0(N )量级 的计算 优化研究中的核心内容之一 。由于边界元法在求 量。因此 ，常规边界元法不适于处理大规模工程问 解开域问题 中的优势，因此已广泛应用于声学灵敏 题 。 度分析中。灵敏度分析方法包括直接微分法DDM 为了提高求解效率，自2O世纪8O年代末至今， (DirectDifferentiationMethod)和伴 随 变量法 将快速多极算法 FMA(FastMultipoleAlgorithm)[妇 AVM(AdjointVariableMethod)Elo3。其 中，直 接 与边界元法相结合形成 的快速多极边界元法有了 微分法利用直接微分得到 的灵敏度边界积分方程 很大的进展，并取得 了广泛应用口_8]。快速多极边 求解 出未知的边界灵敏度值和域内关注点处的灵 界元法的计算量与边界未知量呈线性关系，但是该 敏度值 ，然后代人 目标 函数的灵敏度表达式中求出 算法与特定物理 问题相关 ，需要针对不 同物理场 的