基于表面细分技术边界谱方法.pdfVIP

2005年12月 西北工业大学学报 Dec． 2005 JournalofNorthwestern V01．23No．6 第23卷第6期 PolytechnicalUniversity 基于表面细分技术的边界谱方法 文立华，王卫祥 (西北工业大学航天学院，陕西西安 710072) 摘要：把表面细分技术与边界谱方法相结合，提出了一种基于表面细分技术的子波边界谱方法。 它能精确描述复杂图形边界，提高处理复杂边界问题的能力和计算精度，把表面细分技术建立的几 何模型直接用于边界谱方法的求解。数值算例表明，该方法是有效的。 关键词：表面细分，子波，边界谱方法 中图分类号：0422．2 文献标识码：A 边界积分方程的谱方法是一种全局近似技术， 尝试把表面细分技术用于求解边界积分方程的边界 比一般插值方法的边界元方法有更高的收敛性和精 谱方法中。 度，因而采用高阶多项式或级数逼近(或插值)的谱 方法在解边界积分方程中受到了重视。但是，无论是 1 表面细分 边界元方法还是目前的边界谱方法，它们的系统矩 阵都是非对称满阵。近几年，研究者发展了求解声辐 表面细分是从一个称为“控制网格或点域”的原 射和声散射的子波谱方法。由于子波具有紧支集和 始多面形或多边形开始，将其细分成多个所属的面 高阶消失矩特性，使得积分方程数值求解中的满矩 或多边形，按细分规则排列所有顶点，通过这种细分 阵变为稀疏矩阵，可获得O(NlogN)的算法，而且 算法，得到高精度的曲面与曲线。它的优点是能够快 用子波逼近函数有较快的收敛性，从而在工程分析 速而且精确地描述复杂的曲线和曲面。表面细分方 计算中大大提高了计算效率，降低计算成本。尽管求 法大体可以分为2类：①插值法，细分时点重新计 解未知量采用子波逼近，但边界描述仍然采用传统 算，原控制网点只是得到修正；②逼近法，这类方法 的边界元方法，因此有必要提高现有边界谱方法的 在控制网格细分的时候同时计算控制网点和新产生 描述精度和对复杂边界的描述能力。 的点。本文采用的表面细分方法是插值法。对于控 表面细分技术按照计算机图形学处理方法对复 制网点域z：，1个简单的二维表面细分算法为‘23 杂图形的边界用网格控制点进行描述，通过表面细 z扩1一z? (1) 分算法对边界表面网格进行细分，描述边界曲面和 1 曲线可以达到f”连续。如果把表面细分技术与边界 1U 2)(2) 谱方法相结合，能较大提高构建复杂表面或曲线模 式中，上标k表示细分的次数；式(1)表示原控制网 型的精度和能力，同时也可以把CAD技术与数值 分析模拟直接结合在一起，在边界积分方程求解问 中的控制网格点计算新产生的点。用这种细分方法 题中进行一体化建模。目前将表面细分技术应用于 得到的正弦曲线如图1。 有限元和边界积分方程，还处于初步的概念阶段，公 开的文献只有用于板壳问题的有限元分析[1~3]，用网格和细分网格都是由三角形组成，细分过程如 于求解边界积分方程的尚未见公开报道。因此，本文 图2。 收稿日期：20