浅析有限元方法的发展与应用

浅析有限元方法的发展与应用

1965年“有限元”这个名词第一次在我国出现,到今天有限元在工程上得到广泛应用,经历了三十多年的发展历史,理论和算法都已经日趋完善。有限元法(Finite Element Method,简写为FEM)是求解微分方程的一种非常有效的数值计算方法,用这种方法进行波动数值模拟受到越来越多的重视。

绪论

有限元法是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。它是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。

一、有限元的发展历程

有限元法的发展历程可以分为提出(1943)、发展(1944-1960)和后期(1961-二十世纪九十年代)三个阶段。有限元法是受内外动力的综合作用而产生的。

1943年，柯朗在《美国数学学会公报》(Bulletin of The American Mathematical Society)上发表了《平衡和振动问题的变分解法》(Variational Methods for The Solution of Problems of Equilibrium And Vibration)一文，这篇文章实际上是他1941年在美国数学学会演讲的书面稿，在其中柯朗提出了有限元法的核心思想。大约与柯朗同时，工程师阿格瑞斯在另一个领域独立地提出了有限元法。柯朗和阿格瑞斯各自在数学和工程学领域独立提出了有限元法，他们分别开创了有限元法的数学传统和工程学传统。

有限元法被提出来以后，经过一段时间的沉寂期，在二十世纪五十年代和六十年代初有了很大的发展。主要表现为在代数表达形式、单元划分、单元类型选择和解的收敛性研究上取得的突破。1960年，克劳夫在《平面应力分析中的有限元》(The Finite Element in Plan Stress Analysis)的论文中，第一次从数学上说明了将定义域划分成有限的单元能够成功的原因:他表明对一些特定类型的单元来说，随着单元尺寸的减小，近似解将收敛到精确解，这就在某些情况下证明了有限元法的收敛性。并第一次提出了“有限元法”这个名称，这个名称一直沿用至今，标志着有限元法早期发展阶段的结束。

有限元法后期阶段的发展有国外和国内两条线索。在国外的发展表现为: 第一，建立了严格的数学和工程学基础；第二，应用范围扩展到了结构力学以外的领域；第三，收敛性得到了进一步研究，形成了系统的误差估计理论；第四，发展起了相应的商业软件包。

在国内，我国数学家冯康在特定的环境中独立于西方提出了有限元法。1965年，他发表论文《基于变分原理的差分格式》，标志着有限元法在我国的诞生。冯康的这篇文章不但提出了有限元法，而且初步发展了有限元法。他得出了有限元法在特定条件下的表达式，独创了“冯氏大定理”并且初步证明了有限元法解的收敛性。虽然冯康创造的有限元法不成熟，但他能在当时的条件下独立提出有限元法已十分不易。对于他的这项成就，国内外专家学者和国家领导人都有很高的评价。

半个世纪来，有限单元法蓬勃发展，不仅已经成为结构分析中必不可少的工具，而且成为现象分析的一种手段。其应用已由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定问题、动力问题和波动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、渗流与固结理论、热传导与热应力问题、磁场问题以及建筑声学与噪音问题。不仅涉及稳态场问题，还涵盖材料非线性、几何非线性、时间维问题和断裂力学等。

已出现多种新单元（先后有等参元、高次元、不协调元、拟协调元、杂交元、样条元、边界元、罚单元，还有半解析的有限条等不同单元）和求解方法（如半带宽与变带宽消去法、超矩阵法、波前法、子结构法、子空间迭代法等）。能解决各种复杂耦合问题的软件和软件系统不断涌现。对网格自动剖分和网格自适应过程的研究，大大加强了有限元法的解题能力，使有限单元法逐渐趋于成熟。有限元法作为一种离散化的数值解法，也已成为应用数学的一个新的分支。

二、有限元法解题思路

有限元法分析计算的思路和做法可归纳如下：1.物体离散化；2.单元特性分析；3.单元组集；4.求解未知结点位移。

1.物体离散化

将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算精度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。

2.单元特性分析

(1)选择位移模式，通常，有限元法我们就将位移表示为坐标变量的简单函数。这种函数称为位移模式或位移函数。

(2)分析单元的力学性质，据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵。

(3)计算等效节点力，物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代替所有作用在单元上的力。

3.单元组集

利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程。

4.求解未知结点位移

最后利用已求出结点位移计算各个单元的应力，并经后处理软件整理、显示计算结果。有限元法是应用局部的近似解来建立整个定义域的解的一种方法。先把注意力集中在单个单元上，进行上述所谓的单元分析。基本前提是每一单元要尽可能小，以致其边界值在整个边界上的变化也是小的。这样，边界条件就能取某一在结点间插值的光滑函数来近似，在单元内也容易建立简单的近似解。因此，比起经典的近似法，有限元法具有明显的优越性。

三、有限元方法的实施过程

有限元法在计算机软件中应用尤为突出，可以说有限元法与生俱来就是为了使用计算来进行复杂的计算。有限元方法的实施过程可以分为三个步骤：

1.前处理。将整体结构或其一部分简化为理想的数学模型，用离散化的网格代替连续的实体结构。

2.计算分析。分析计算结构的受力、变形及特性。

3.将计算结果进行整理和归纳。

对于有限元程序使用者而言，第一步和第三步的工作量最大，一个有限元程序的好坏，在很大程度上取决于第一步的前处理和第三步的后处理功能是否强大。

前处理：对于第一步的前处理而言，要根据计算的目的和所关心的区域，将结构模型化、离散化。需要给出下列信息：

(1)节点的空间位置。

(2)单元与节点的连接信息。

(3)结构的物质特性和材料参数。