Galerkin 边界元法用于弹塑性分析的准高次元方法1).pdfVIP

第　29 　卷　第　6 　期 力 　学 　学 　报 　 Vol. 29 ,No. 6 1997 　年　11 　月 ACTA MECHANICA SINICA Nov. , 1997 Galerkin 边界元法用于弹塑性分析的 准高次元方法1) 刘清　 ( 北京自动测试技术研究所, 北京 100088) 岑章志 　　徐秉业 (清华大学工程力学系 , 北京 100084) 摘要 　提出了一种适用于 Galerkin 边界元法的高次单元插值方法和半解析半数值的积分方 法准高次元法 ,建立了有关数值模型和将该方法应用于结构弹塑性分析的算法模型. 有关算 例结果表明 ,本文建议的方法是切实可行的. 关键词 　边界元法 ,弹塑性分析 ,数值模型 引　言 近几十年来 ,B IEBEM 获得了很大进展 ,也得到了广泛的工程应用 ,成为有限元法的重要补 充. 边界元法的最大特点是它降低了求解问题的维数 , 同时由于它采用了基本解来描述真实状态 的位移、应力和应变场 ,一般它的求解量有较高的精度 , 因而可以说边界元法有精度较高和计算 量小的特点. 但是 ,必须考虑到边界元法最后形成的求解矩阵为非对称满阵 ,对其矩阵分析求解 的计算量是非常大的 ,特别是在采用高度协调的边界元时所涉及的积分运算量尤为如此. 在这种 背景下 ,一些学者进行了产生对称矩阵的边界元法的探索 ,最初见诸发表的是 1979 年 Sirtori 的 工作[1 ] ,加之后来 Costable 和 Wendland 的工作[2 ,3 ]和 Maier , Polizzoto 及 Novati 的工作[4～6 ] 构成 了用 Galerkin 边界元法处理弹塑性分析问题的基本思路. 中国学者余德浩在 80 年代进行了‘‘自 然边界元法’’方面的探索[10 ] ,该方法对于规则边界结构分析可以产生一个二重积分的、可以进 行解析计算的边界积分方程 ,因此降低了边界元法的求解难度 , 同时提高了数值结果的精度. Galerkin 边界元法可以得出一个对称的方程系数矩阵 ,这对方程求解和结构力学性态分析 都提供了很大方便. 但是 , Galerkin 边界元法需要一个二重积分才能求得其系数矩阵 ,给单元划 分和方程系数矩阵的建立带来了巨大困难 ,Sirtori 和 Maier 等人利用复变函数方法对位移的线 性元及面力和应力的常数元找到了相关的解析积分表达式 ,这解决了二重积分所带来的困难 ,为 Galerkin 边界元法理论确立了一条成功的算法. 理论分析与数值实例表明 ,这一算法无论从计算 时间上和计算精度上 ,在求解大部分二维问题时都不逊于传统边界元法. 但是 ,该算法在结构与 加载方式的适应性方面也存在着很多缺陷 ,这给 Galerkin 边界元法的推广发展与工程应用都带 来了很大困难 ,划分单元的工作也比较麻烦 ,无形中降低了这一算法的实用性. 本文建议了一种 二重插值的算法 ,又可称为‘‘准高次元法’’,即在初始插值时采用高次元格式对结构和位移、面 1) 国家自然科学基金资助项目. 　收到第一稿 收到修改稿. © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 746 　　　　力 　　学 　　学 　　报 1997 年 第 29 卷 力及应力进行描述 ,而在进行数值积分时则要进一步对单元进行细化 ,然后在细化的单元中进行 位移线性插值和面力、应力的常数插值 ,利用积分中值定理对高次插值函数取积分中值 ,而对核 函数继续采用解析积分. 本文介绍二重插值方法的一些基本公式和它在弹塑性分析问题上的应 用. 1 　Galerkin 边界元法中的二重插值与积分技术 对称型 Galerkin 边界元法用于求解弹塑性分析问题的基本公式是基于