边界元法边界层处理的一种新方法——拟多连域法.pdfVIP

第 14 卷　第 1 期 应 用 力 学 学 报 . 14　 . 1 V o l N o 1997年3月 CH INESE JOURNAL OF APPL IED M ECHAN ICS M ar. 1997 边界元法边界层处理的一种新方法 ——拟多连域法 王元丰　　　　　邹永超 ( ) ( ) 北方交通大学　北京　100044 　 哈尔滨船舶工程学院　哈尔滨　150001 摘　　要 提出了处理边界元法边界层问题的新方法。这种方法在求得边界未知量后, 在距 边界较远的域中构造一工作边界, 将实际边界与工作边界看作一多连域问题求得工 作边界的未知量。再将工作边界与边界层边界视为一多连域问题求解, 可得到满足精 度要求的边界附近点的位移与应力。这种方法理论简洁、计算方便、有效、精度高, 对 于需求边界上多点的未知量问题很具优越性。 关键词: 边界元法; 边界层; 拟多连域法 1　引　言 应用边界元法求解域内点的位移、应力时, 由于积分的核函数具有奇异性, 当观察点移到 边界近傍, 即当观察点与边界单元的距离小于单元的长度时, 如仍采用一般的内点未知量求解 方法, 解的精度将下降, 特别是当观察点接近边界时, 解呈发散状。这就是边界元法的边界层效 应问题。 尽管边界层效应的存在对边界元法求域内点的位移和应力有很大限制, 但有关这方面的 文献并不很多[ 17 ] , 处理边界层效应问题的方法大体有三种: 第一种是把靠近内点的边界单元 局部细分; 第二种是改善基本解的积分精度; 第三种方法是基于刚体运动准则, 利用一些特解 间接确定奇异积分。基于一般多连域问题的边界元求解法, 提出了求边界层上点未知量的拟多 连域法。 2　拟多连域法 对于图1所示的一般单连域问题为了求边界层点的位移, 可先在距边界较远的域中构造一 工作边界 N 。将实际工作边界 与工作边界 N 看作一多连域问题的内、外边界。在由边界积 分方程求得边界的所有未知量后, 按多连域问题求工作边界上点位移与面力增量。 在 和 N 上分别取A 、B 两点, 这样对于由 ∪N 组成的多连域问题的边界积分方程可 来稿日期: 修回日期: © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 第 1 期 边界元法边界层处理的一种新方法——拟多连域法 125 ( ) 写为 不计体力 B m m m C V m = (V P - P V ) d + (V P lm ∫l m m l ∫l m A m m - P m lV ) d + (V l P m ∫ N A