三维无网格法与有限元法工程应用对比分析研究(排版文例).doc

三维无网格法与有限元法工程应用对比分析研究 徐镇凯 马连军 （南昌大学 建筑工程学院 江西 南昌 330029） 摘 要：本文阐述了三维无网格伽辽金法和非线性有限单元法的基本理论。无网格伽辽金法是在移动最小二乘法的基础上，建立全计算域的高阶连续可导的插值函数，只利用节点信息来建立离散模型的平衡方程。非线性有限单元法是在离散单元的基础上，构造连续可导的插值函数，进而建立平衡方程。论文推导了其相关的理论，编制了计算程序，通过对算例的分析表明：无网格伽辽金法的计算精度高于有限单元法，已经成为有限单元法的有力补充。 关键词：移动最小二乘法，变分原理，三维无网格伽辽金法，非线性本构关系，非线性弹性有限单元法，程序设计 3D Element-free Galerkin method and Nonlinear Finite Element Method of application and contrastive research in engineering XU Zhen-kai，MA Lian-jun (Civil Engineering and architecture school, Nanchang University, Nanchang 330029,china) Abstract The relevant theories of three-dimensions Element-free Galerkin method(EFGM) and Nonlinear Finite Element method(FEM) are introduced in this paper. EFGM is based on moving least square method(MLSM). EFGM formulate the discrete model with using only a set of nodes, so the pre-processing of the method is simplified and the approximating function is continuous with high-order differentials over the entire domain. Nonlinear elastic finite element method is based on deriving elements, formulate continuous with high-order differentials functions and equilibrium equations etc. Furthermore a example is analyzed in details by using promgram. The result show that new EFGM is accurate than FEM for problems in engineering. Key words: moving least square method; variational principle; three-dimensions Element-free Galerkin method; increment-type stress-strain relation; Nonlinear elastic finite element method; program design 1 前言 有限元法[1]是工程数值分析领域的较广泛、重要、成熟的数值方法,伴随计算机科学的飞速发展，越来越发挥它的强大功效，它具有对复杂几何构形的适应性、各种物理问题的可应用性、建立严格理论基础上的可靠性及适合计算机实现的高效性。有限单元法为解决重大的工程问题作出卓越的贡献，但是随着科技的发展,有限单元法也暴露出自身的一些缺陷和实际应用中的困难问题：（1）对于三维复杂结构的建模耗时，前处理困难；（2）其计算结果通常是位移连续，应力、应变不连续，需要后处理修匀；（3）当结构材料不可压缩时出现“体积自锁”现象；（4）在处理网格方面，出现网格畸变问题。在具体应用中，如动态裂纹扩展和结构大变形等方面存在困难问题。 无网格法[2]是一种新兴的数值方法,它只需要节点信息。比较早的无网格法是光滑粒子动力学法[3][4]，又称SPH法，为得到更佳的计算结果，将移动最小二乘法引入Galerkin法中[5]，提出无网格Galerkin法[6],即无网格伽辽金法（EFGM）。EFGM将结构体离散成一系列的节点，在整个计算域上建立高阶连续可导的插值函数，进而在节点的影响域上通过影响半径建立节点与计算点的关系。与传统的有限单元法相比，无