有限元法无单元法及自然单元法之比较研究.PDF

第24 卷 第5 期 岩石力学与工程学报 Vol.24 No.5 2005 年3 月 Chinese Journal of Rock Mechanics and Engineering March ，2005 有限元法、无单元法及自然单元法之比较研究 1 1 2 卢 波 ，葛修润 ，孔祥礼 (1. 中国科学院 武汉岩土力学研究所岩土力学重点实验室，湖北 武汉 430071 ； 2. 吉林省公路勘测设计院，吉林 长春 130021) 摘要：无单元法与有限元法的根本区别在于形函数的构造方法不同，因而得到的形函数形式也不同，但可统一于 单位分解。另外，采用 Galerkin 法对偏微分方程求解时需要对弱形式在求解域内积分，研究表明有限元法和无单 元法中对弱形式的积分可统一于单位分解积分方法。无单元法解决了有限元法中前处理复杂以及难以有效处理的 诸如裂纹扩展和大变形等问题；然而自身也存在着计算量大、边界条件和不连续的场变量及其导数处理上的困难。 自然单元法是新近出现的一种求解PDE 的数值方法，它兼有无单元的特性和有限元的优点，可以认为是介于两者 之间的一种极具发展前途的数值方法。 关键词：数值分析；有限元法；无单元法；自然单元法；单位分解；弱形式；单位分解积分 中图分类号：O 241.8 文献标识码：A 文章编号：1000 –6915(2005)05 –0780 –07 COMPARATIVE STUDY ON FINITE ELEMENT METHOD，MESHLESS METHOD AND NATURAL ELEMENT METHOD 1 1 2 LU Bo ，GE Xiu-run ，KONG Xiang-li (1. Institute of Rock and Soil Mechanics ，The Chinese Academy of Sciences，Wuhan 430071 ，China； 2. Jilin Provincial Highway Survey and Design Institute ，Changchun 130021，China) Abstract ：The fundamental difference between the finite element method(FEM) and the meshless method(MM) lies in the difference of the construction methods of the shape function. But the two kinds of shape functions both serve as a partition of unity. Integration of the weak form over the problem domain is needed when solving the partial differential equation using the Galerkin method. The newly presented quadrature method based on partition of unit provides a unified approach for the integration of the weak form in FEM and MM. MM is advantageous over the traditional FEM in the problems involved with th