无单元法与有限元法及工程中对比分析.pdfVIP

无单元法与有限元法及工程中的对比分析 马连军 南昌大学建筑工程学院，南昌（330029 ） E-mail：hxh501@126.com 摘 要：本文分别介绍了三维无单元伽辽金法和非线性有限单元法及其相关理论。无单元伽 辽金法是在移动最二乘法的基础上，建立全计算域的高阶连续可导的插值函数，只利用节点 信息来建立离散模型，基于变分原理来推导无单元伽辽金法的平衡方程。非线性弹性有限单 元法是在划分单元基础上，构造连续可导的插值函数，根据变分原理，由最小位能原理方程 推导平衡方程。论文阐述了移动最小二乘法的基本原理和权函数、节点影响半径的选取问题， 同时推导全量形式的非线性弹性本构关系，并且编制了三维无单元伽辽金法和非线性弹性全 量有限单元法的程序，对悬臂梁算例进行分析，通过计算值与理论值比较表明，新兴的无单 元伽辽金法在某种领域上，计算精度高于有限单元法，已经成为有限单元法的有力补充。 关键词：移动最小二乘法，三维无单元伽辽金法，非线性本构关系，非线性弹性有限单元法， 程序设计 中图分类号：TU375 文献标识码：A 1.前言 [1] 有限单元法 是工程数值分析领域的最广泛、重要、成熟的数值方法,伴随计算机科学的 飞速发展，越来越发挥它的强大功能，具有对复杂几何构形的适应性、各种物理问题的可应 用性、建立严格理论基础上的可靠性及适合计算机实现的高效性。它建立于20 世纪,为解决 重大的科学工程问题作出卓越的贡献.随着科学的发展,有限单元法遇到了许多的困难问题. 例如动态裂纹扩展、材料破坏及失效、应变局部化、高速撞击，以及大变形问题，利用有限 单元法分析时出现网格畸变等问题。 无单元法[2]是一种新兴的数值方法,它只需要节点信息。与传统的有限单元法相比，无单 元法的近似函数不依赖于网格，在分析大变形问题时具有优势，同时它克服了有限元需要划 分大量的单元及单元重构的缺点，减少了工作量，具有精度高、收敛快的优点。比较早的无 单元法是光滑粒子动力学法[3][4]，又称SPH 法。为在分析中得到更佳的结果，将移动最小二 [5] [6] 乘法引入 Galerkin 法中 ，后来经改进与保留提出无单元 Galerkin 法 ,即无单元伽辽金法 （EFGM ）。它是通过几个不相关的数值来模拟函数。无单元伽辽金法是与有限单元法相类 似的无单元数值方法，它采用移动最小二乘法拟合形函数，由拟合函数的自由度远远少于 插值函数的基点数，具有较高的精度，又能保证应力场的连续性，同时该方法的前后处理比 较方便。 2. 三维无单元伽辽金法(EFGM) 2.1 移动最小二乘法的基本原理 移动最小二乘法是用加权最小二乘法来近似模拟场函数。设场函数为u(x) ,其中 x (x , y , z )T 为场点，即变量。给定场函数u(x) 在节点上的值为： u(x ) u i 1,2,3,……, n (1) i i 用移动最小二乘法，场函数可近似的表示为： - 1 - m h T u (x) = ∑p (x )a (x ) p (x )a (x ) (2)