快速多极边界元法在薄板结构中的应用.pdfVIP

维普资讯 第 28卷 第 2期 燕 山大学学报 Vb1．28 NO．2 2004年3月 JoumalofYanshanUniversity M ar． 2004 快速多极边界元法在薄板结构中的应用 赵丽滨 姚振汉 摘 ’要 基于 Taylor级数多极展开研究了边界元快速多极算法 (FM—BEM)，并将它应用于薄板结构 。算例分析 表 明FM—BEM 的计算时间和存储空间明显少于常规边界元迭代解法。随着问题规模的增大，这种优势将更加突 出。 关键词 边界元法，Taylor级数，数值方法，多极展开。 ApplicationofFastMultipoleBoundaryElement MethodforThinPlateStructures Abstract BasedonthemultipoleexpansionusingTaylorseries，fastmultipoleboundaryelementmethod (FM—BEM) isstudiedandappliedtothinplatestructures．Numericalexamplesshow thatFM —BEM haslesscomputationcostand memoryrequirementthanthestandardBEM withitemtivesolve~ Furthermore，thedifferenceincreaseswiththesizeof problem． Keywords boun daryelementmethod，Taylorseries，numericalmethod，multipoleexpansion． 0 引言 计算表面积与体积 比较大的薄型结构 问题成为可 能。本文将基于 Taylor级数多极展开的边界元快速 边界元方法与有限元方法相比，具有两方面突 多极算法应用于薄板结构，算例分析表明：在保证 出的特点：一是降低了求解问题的维数、具有较高 较高计算精度的同时，它的计算时间和存储空间明 的计算精度；二是边界元法求解代数方程组的系数 显少于常规边界元迭代解法，问题的规模越大这种 矩阵是满阵，一般情况下还是非对称矩阵，这给边 差距也越大。 界元法的解题规模带来了很大的限制。虽然近年来 发展了一些迭代算法，如 GMRES，Bi—CGSTAB 1 三维弹性静力学问题边界积分方程 等方法，可以将运算量由直接求解的D()降低到 不计体积力情况的三维弹性静力学边界积分 了DC?v2)的量级，但计算量仍然很大，这依然阻碍 方程为 了边界元的进一步发展。 近年来，国际上发展起来的快速多极边界元法 c ) )+ ，y)uj(y)d／ly)= 为边界元法研究者展示了广阔的应用前景。这种算 法用节点集群的多极展开式来近似表示核函数与 ， y)tj(y)dl-~) (1) 远场边界变量乘积的边界积分，从而将边界元法的 其中，厂为边界， ，Y分别为边界上核函数的源点 满阵方程组化为具有稀疏矩阵的方程组。对存储空 和边界场点，c )是与 处边界几何特征有关的系 间和求解时间两方面的需求都从原来的DC?v2)减少