《计算物理学》课件 第六章 有限元方法.pptx

第六章有限元方法;有限元方法是数值求解复杂微分方程的一种非常有效的方法，其基于变分原理。

将求解微分方程的问题转化成一个泛函求极值的变分问题。

与有限差分法的规则网格划分求解域不同，有限元方法采用剖分逼近的离散化方式，将求解域剖分为有限个基本块，称为“单元”，如三角形，然后统一编号并求解。

有限元方法在处理具有复杂区域或复杂边界条件下的数学物理问题更具优势，如飞行器的结构设计、电磁场问题、建筑物的结构受力分析等。

有限元方法是传统Rayliegh-Ritz变分方法的发展，并融合了差分法的优点，在问题的处理上统一、适应能力强，已广泛应用于科学研究和工程设计。;发展历史;1941年，加拿大工程师Hrennikoff在《应用力学学报》上发文介绍了膜和板模型，并将求解域离散化为一个格栅结构网络。这篇论文被认为是有限元法诞生的起点。

同年，Courant在美国数学学会上报告了关于用变分法数值求解圆柱体扭转问题的二阶偏微分方程的工作，首次提出了可在定义域内利用三角形分片地使用展开函数来表达其上的未知函数的Rayleigh-Ritz方法，给出了有限元方法的数学基础。论文发表于1943年。

20世纪50年代，J.H.Argyris、M.J.Turner、R.W.Clough、H.C.Martin等提出了矩阵刚度法，也就是最早的有限元方法，用来研究航空、机械、土木工程中的实际问题。

随后，M.J.Turner、R.W.Clough、H.C.Martin和L.J.Topp等人首次将飞机的机翼离散为三角形板块的集合来进行平面应力的分析，取得了成功。;1960年，Clough首次将这种方法命名为有限元方法，准确地概括了该方法。

1963年，J.F.Besseling、R.J.Melosh和R.E.Jones等证明变分原理是有限元方法的数学原理。

20世纪60年代，冯康提出了基于离散化的变分原理的差分格式，来数值求解椭圆型偏微分方程，独立提出了有限元分区域的思想。

钱伟长以及胡海昌提出了广义变分原理。

第一个有限元开源软件是由加州大学伯克利分校的E.L.Wilson于1958年开发的，是一个基于矩形平面应力有限元的自动化程序。;有限元方法的基础和原理;泛函和变分;?;?;?;变分原理;?;?;?;变分有限元方法;剖分逼近是有限元离散化的手段。

采用特???的结构单元将求解域剖分为有限个单元域，然后在每个单元内通过Rayliegh-Ritz方法计算相应的函数，通过单元上的插值逼近，得到一个结构简单的函数集，进而整合得到整个求解域的解。

这里的离散单元一般采用具有规则形状的几何结构，如三角形、四面体等。;有限元方法包括以下几个主要步骤：

（1）将微分方程的边值问题转化为求极值的变分问题，也就是推导出与微分方程等价的泛函形式。

（2）对求解区域进行离散化，也就是针对不同维度的问题，采用相应的单元将求解区域分为有限个单元，然后对所有单元和节点进行编号。

（3）利用公式计算出各个三角形元素的系数矩阵。;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;偏微分方程边值问题的有限元解法;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;?;有限元方法与有限差分法的异同

两种方法处理问题的思路不同。

有限差分法需要列出所要研究问题的微分方程及定解条件，然后通过通过规则网格划分将微分方程的求解离散化，得到差分方程组，再进行数值求解。

有限元方法是基于变分原理，将所要求解的微分方程转化为对泛函求极值的变分问题。再将求解区域划分为有限个单元的离散化方法，通过构造单元的插值函数，将变分方程问题离散化为线性方程组，再进行数值求解。;两种方法对区域的离散化方式也大不相同。

有限差分法通常采用矩形网络区域划分，如果所求解区域不规则，有限差分法就很难实现网络格点与边界的良好逼近，这时需要对边界采用一些近似方法来处理。

而有限元方法一般采用三角形划分，这样的划分方式能够实现对节点在区域内的任意配置，并且配置方式可以根据边界条件的情况来进行选择和优化，灵活方便。对于求解区域具有复杂边界形状的问题，有限元方法可以选择单元节点完全处在求解区域的边界上，从而可以使单元划分很好地逼近边界。;与有限差分法相比，有限元方法的适用范围小很多。

对大多数微分方程都可以采用有限差分法进行数值求解，尤其是对于有规则求解区域的问题。

对于有限元方法，应用比较广的是对椭圆型偏微分方程的数值求解，而很少被用于求解双曲型偏微分方程。;ANSYS