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

浅析有限元方法的发展 摘要：1965年“有限元”这个名词第一次在我国出现,到今天有限元在工程上得到广泛应用,经历了三十多年的发展历史,理论和算法都已经日趋完善。有限元法(Finite Element Method,简写为FEM)是求解微分方程的一种非常有效的数值计算方法,用这种方法进行波动数值模拟受到越来越多的重视。 关键字：有限元法，发展，应用 前言 有限元法是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。它是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。 1．有限元的发展历程 有限元法的发展历程可以分为提出(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年，他发表论文《基于变分原理的差分格式》，标志着有限元法在我国的诞生。冯康的这篇文章不但提出了有限元法，而且初步发展了有限元法。他得出了有限元法在特定条件下的表达式，独创了“冯氏大定理”并且初步证明了有限元法解的收敛性。虽然冯康创造的有限元法不成熟，但他能在当时的条件下独立提出有限元法已十分不易。对于他的这项成就，国内外专家学者和国家领导人都有很高的评价。 半个世纪来， 有限单元法蓬勃发展， 不仅已经成为结构分析中必不可少的工具，而且成为现象分析的一种手段。其应用已由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定问题、动力问题和波动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、渗流与固结理论、热传导与热应力问题、磁场问题以及建筑声学与噪音问题。不仅涉及稳态场问题，还涵盖材料非线性、几何非线性、时间维问题和断裂力学等。 已出现多种新单元（先后有等参元、高次元、不协调元、拟协调元、杂交元、 样条元、边界元、罚单元，还有半解析的有限条等不同单元）和求解方法（如半 带宽与变带宽消去法、超矩阵法、波前法、子结构法、子空间迭代法等） 。能解决各种复杂耦合问题的软件和软件系统不断涌现。 对网格自动剖分和网格自适应过程的研究，大大加强了有限元法的解题能力，使有限单元法逐渐趋于成熟。 有限元法作为一种离散化的数值解法，也已成为应用数学的一个新的分支。 2.有限元法解题思路 有限元法分析计算的思路和做法可归纳如下：1）物体离散化；2）单元特性分析；3）单元组集；4）求解未知结点位移。 1）物体离散化