题解无网格法及和有限元法的比较的研究初步.pdfVIP

摘要 近几十年来，有限元法(FEM)由于其通用性和灵活性已经成为工程数值领域的主要方法， 但是有限元法在分析大变形、不连续性等问题时存在缺陷。无网格法(MLM)正是基于这些缺 陷提出的。本文在分析了无网格法的研究历史和现状、分类和优缺点的基础上，基于算例分 别采用无网格伽辽金法(EFGM)和有限元法进行数值模拟，对两者在结构体离散、刚度矩阵 建立、等效节点荷载施加、边界条件的引入等方面进行了比较分析，并指出了各自的优缺点。 基于一维算例，在相同条件下，分别对无网格伽辽金法中影响其精度的下述因素：布点 方案、基函数、权函数、形函数以及支撑域半径等方面进行分析，得到了各因素的经验参数， 具有一定的指导意义。相比有限元法单元的划分，无网格法布点更为简单，后处理也比较方 便。权函数的选择非常关键，在几种常用的权函数中，选取四次样条型权函数得到的结果最 为精确。运用点插值法建立形函数相比移动最小二乘无网格伽辽金法具有点插值性，更容易 施加边界条件。支撑域半径的选取决定着支撑域的大小，也是影响计算效率和精度的一个重 要因素。 文章简述了无网格法中另外的比较成熟的方法——配点法，并基于一维算例在相同条件 下分别与无网格伽辽金法计算结果进行了对比。相比于无网格伽辽金法，配点法计算效率高、 位移边界条件容易实现、但系数矩阵是非对称的、稳定性较差。 本文简单介绍了克立格法，并将其应用到求解无网格法形函数的过程中，建立了克立格 无网格伽辽金(KEFG)法(简称克立格无网格法)。用克立格无网格伽辽金法分析了一维杆件 算例，并与相同材料参数条件下运用无网格伽辽金法得到的结果进行对比，结果表明该方法 具有很好的稳定性和点插值性，方便施加位移边界条件，在固定端的误差相对于无网格伽辽 金法要小，且计算量远小于无网格伽辽金法。 最后运用无网格伽辽金法求解了土体的二维固结问题，并与相同条件下运用有限元法得 到的结果进行比较。结果表明两种方法计算结果比较吻合，说明该方法可以有效地模拟工程 问题。 关键字： 无网格伽辽金法；有限元法；配点法；克立格法。 I Abstract In recent decades, the finite element has become a major field of engineering numerical methods, as it is universal and flexible. But defects of the finite element method exist in analyzing large deformation, non-continuity and so on. The meshless method is proposed basing on these defects. The element-free Galerkin method (EFGM) and finite element method (FEM) are adopted to numerically simulate a numerical example basing on the development history and present situation, classification, advantages and disadvantages analysis of the meshless in this text. And comparisons are also made between the results by this two methods basing on the numerical examples in structure discrete, stiffness matrix establishment, application of equivalent node loads and boundary conditions application, etc. And theirs advantages and disadvantages are pointed out. Under the same conditions, the influencing facts of EFGM such as nod