有限元
在数学中,有限元法(FEM,Finite Element Method)是一种为求解偏微分方程边值问题近似解的数值技术。求解时对整个问题区域进行分解,每个子区域都成为简单的部分,这种简单部分就称作有限元。
它通过变分方法,使得误差函数达到最小值并产生稳定解。类比于连接多段微小直线逼近圆的思想,有限元法包含了一切可能的方法,这些方法将许多被称为有限元的小区域上的简单方程联系起来,并用其去估计更大区域上的复杂方程。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的(较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。
中文名:有限元
外文名:Finite Element
所属领域:数学
应用学科:数学
方 法:有限元法
应用范围:热传导 电磁学 土力学
国内发展:冯康先生在50年代提出
目录
计算简介
英文:Finite Element
有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域—飞机结构静、动态特性分析中应用的一种有效的数值分析方法,随后很快广泛地应用于求解热传导、电磁场、流体力学等连续性问题。
方法
有限元法分析计算的思路和做法可归纳如下:
物体离散化
将某个工程结构离散为由各种单元组成的计算模型,这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来;单元节点的设置、性质、数目等应视问题的性质,描述变形形态的需要和计算精度而定(一般情况单元划分越细则描述变形情况越精确,即越接近实际变形,但计算量越大)。所以有限元中分析的结构已不是原有的物体或结构物,而是同新材料的由众多单元以一定方式连接成的离散物体。这样,用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理,则所获得的结果就与实际情况相符合。
选择位移模式
在有限单元法中,选择节点位移作为基本未知量时称为位移法;选择节点力作为基本未知量时称为力法;取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化,所以,在有限单元法中位移法应用范围最广。
当采用位移法时,物体或结构物离散化之后,就可把单元总的一些物理量如位移,应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常,有限元法我们就将位移表示为坐标变量的简单函数,这种函数称为位移模式或位移函数。
分析力学性质
根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等,找出单元节点力和节点位移的关系式,这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式,从而导出单元刚度矩阵,这是有限元法的基本步骤之一。
等效节点力
物体离散化后,假定力是通过节点从一个单元传递到另一个单元。但是,对于实际的连续体,力是从单元的公共边传递到另一个单元中去的。因而,这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去,也就是用等效的节点力来代替所有作用在单元上的力。
单元组集
求解
解有限元方程式得出位移。这里,可以根据方程组的具体特点来选择合适的计算方法。
通过上述分析,可以看出,有限单元法的基本思想是"一分一合",分是为了就进行单元分析,合则为了对整体结构进行综合分析。
发展概况
943年 courant在论文中取定义在三角形域上分片连续函数,利用最小势能原理研究St.Venant的扭转问题。
1960年 clough的平面弹性论文中用“有限元法”这个名称。
1965年冯康发表了论文“基于变分原理的差分格式”,这篇论文是国际学术界承认我国独立发展有限元方法的主要依据 [1] 1970年 随着计算机和软件的发展,有限元发展起来。
1975年 谢干权发表论文“三维弹性问题的有限单元法”,标志着我国学者在与世隔绝的情况下,独立发展出真正能应用于三维实践的有限元方法和有限元工程软件 。文中还在全世界率先得到了三维有限元的超收敛结果。[2]
涉及的内容
有限元所依据的理论,单元的划分原则,形状函数的选取及协调性。
限元法涉及
数值计算方法及其误差、收敛性和稳定性。
应用范围
求解的情况
杆、梁、板、壳、块体等各类单元构成的弹性(线性和非线性)、弹塑性或塑性问题(包括静力和动力问题)。能求解各类场分布问题(流体场、温度场、电磁场等的稳态和瞬态问题),水流管路、电路、润滑、噪声以及固体、流体、温度相互作用的问题。
视频
有限元分析数值解的精度与性质
参考文献
- ↑ [冯康, 基于变分原理的差分格式, 应用数学与计算数学,1965, 2 (4): 238-262]
- ↑ 中国知网.1975,引用日期2017-02-27