虚边界元法的理论分析.pdfVIP

第 17 卷第 1 期 计 算 力 学 学 报 V o l. 17 N o. 1 2000 年 2 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS February 2000 文章编号:(2000) 虚边界元法的理论分析 1 2 2 张耀明 , 　孙焕纯 , 　杨家新 ( 1 山东工程学院数理系, 淄博 255012; 2 大连理工大学力学系, 大连 116023) 摘　要: 通过平面位势问题, 对虚边界元方法的理论基础进行了探讨性研究。在此基础上, 给 出虚边界与真实边界间的距离选取法, 其实质是距离的选取依赖于边界的离散化。 关键词: 虚边界元; 平面位势问题 中图分类号: O 342　　　文献标识码:A 1　引　　言 虚边界元法[ 1～ 4 ] 的最大优点是能够避免求解边界量时遇到的奇异积分和求解边界附近区 域物理量时遇到的拟奇异积分的计算。平面和空间的数值算例均表明此方法具有较高的计算 精度和效率。虚边界元法的应用已相当广泛, 不仅应用于静力问题, 而且应用于动力问题的求 解。然而, 虚边界元法的理论依据是什么? 如何选取虚边界和真实边界间的距离等诸多问题还 有待进一步地研究。本文通过平面位势问题, 对虚边界元方法的理论基础进行了探讨性研究。 文[5 ]提出域外奇点法, 但它们是两种截然不同的方法。 在虚边界元法中, 真实边界和虚边界间的距离选取是一个十分重要的问题, 它直接影响计 算结果的精度和可靠性。本文通过对方法的分析, 给出距离选取公式, 其特点是距离的选取与 边界的划分密切相关。数值算例表明, 这些公式是完全正确的。 2　虚边界元方法概述[ 1～ 4 ] 设所研究的区域为 , 其边界为 , 设想 被嵌入在一个无限的区域中, 在 外的无限区 ( ) 域中有一延拓边界 ′这里称为虚边界, ′所围的区域为 ′, 且 ′ , 沿着边界 ′有一个 ( ) ( 待定的虚拟密度函数 x 在位势问题中为虚拟的分布源, 在弹性力学问题中为虚拟分布 ) 力 , 令此虚拟密度函数对真实边界所产生的位势或通量、位移或面力与边界条件一致。从而达 到求解待定密度的目的。以上称之为虚边界元法。 位势问题中的虚边界元方程为 ( ) ( ) ( ) ( ) u x = y u x , y d ′, 　x ∈ 1a ∫ y u (x ) u (x , y ) = (y ) d ′, 　x ∈ ( 1b) ∫ y n ′ n 　　计算内点位势和通量的边界积分表示式为 u (x ) = (y ) u (x , y ) d ′,