有限元之我见.pdfVIP

有限元之我见 雒海涛 北京索为高科系统技术有限公司 2013/1/29 会工作，更要会思考 一、 概述 写此文档的目的是理顺一下自己对有限元的一个理解，长期以来对有限元的理解有些概 念还是很模糊的，对后面的学习也造成了一些影响，希望在这个梳理过程中能对自己以后继 续深入学习有限元相关的知识起到帮助；另一方面这些理解都是基于本人的理解，难免有偏 颇之处，也希望对有限元有深入理解的朋友能对文中出现的错误之处指正，文中很多内容也 是来源于网络。 二、 为什么要用有限元 在力学中，目前发展而言，很多情况下我们只能对很少的构件得出精确的应力解答，如 轴对称问题，无限域问题等。有限元就是对构件划分单元，利用泛函的思想求得近似解，对 实际问题的应用更加广泛。 这儿我的一个理解是在求解力学问题的时候，尤其是求解复杂系统的时候很难对一些大 型计算求得其解析解，进而采用数值解。关于解析解和数值解这儿带两句，因为有些人可能 没有学过数值分析的课程，就拿一个方程组而言，一般来说两个方程两个未知数的方程组来 说是很容易求得其两个未知数的值的，但是对于大规模的系统需要罗列很多方程，有些是有 重复的要去掉，很多未知数那么求解量就很大了，这个时候才用数值求解的方式，给定初始 值然后给每个参数一个步长增量来求其近似解就会更容易一些。 在说到有限元上，有限元就是将连续的求解域离散为一组单元的组合体，然后对一 组单元的组合体进行求解，根据各种方程建立起代数方程组就需要用到了数值求解问 题。 三、 什么是有限元 有限元的原理：将连续的求解域离散为一组单元的组合体，用在每个单元内假设的近 似函数来分片的表示求解域上待求的未知场函数，近似函数通常由未知场函数及其导数在单 元各节点的数值插值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度 问题。 有限元的概念：有限元法（finite element method ）是一种高效能、常用的计算方法。 有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊 松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。自从1969 年以 来，某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得 了有限元方程，因而有限元法可应用于以任何微分方程所描述的各类物理场中，而不再要求 这类物理场和泛函的极值问题有所联系。基本思想：由解给定的泊松方程化为求解泛函的极 值问题。 既然概念里面提到了拉普拉斯方程和泊松方程，这儿先解释一下这两个方程是干嘛的， 我的理解是这两个方程是以势函数的方式描述了电场、引力场和流场，也就是说这两个方程 是用函数的形势描述真实存在的物理场，或者说真实物理场的数学表达。拉普拉斯和泊松方 程的定义如下： 四、 有限元法结算物理问题的一般步骤 步骤1：剖分: 将待解区域进行分割，离散成有限个元素的集合．元素（单元）的形状原则上是任意的．二 维问题一般采用三角形单元或矩形单元，三维空间可采用四面体或多面体等．每个单元的顶 点称为节点（或结点）． 步骤2：单元分析: 进行分片插值，即将分割单元中任意点的未知函数用该分割单元中形状函数及离散网格 点上的函数值展开，即建立一个线性插值函数 步骤3 ：求解近似变分方程 用有限个单元将连续体离散化，通过对有限个单元作分片插值求解各种力学、物理问题 的一种数值方法。有限元法把连续体离散成有限个单元：杆系结构的单元是每一个杆件；连 续体的单元是各种形状（如三角形、四边形、六面体等）的单元体。每个单元的场函数是只 包含有限个待定节点参量的简单场函数，这些单元场函数的集合就能近似代表整个连续体的 场函数。根据能量方程或加权残量方程可建立有限个待定参量的代数方程组，求解此离散方 程组就得到有限元法的数值解。有限元法已被用于求解线性和非线性问题，并建立了各种有 限元模型，如协调、不协调、混合、杂交、拟协调元等。 五、 方法、算法和求解器 提到有限元就会提到某某算法，现在很多人把方法和算法混淆在一起，根据个人理解这 儿对方法、算法进行整理一下。 六、 方法 我觉得方法是求解实际物理问题的办法，例如有限元方法，这儿除了有限元方法还有很 多其他的方法，列举如下： 显示/ 隐式有限元法：所谓显式和隐式，是