有限元法调研

有限元法概念

基本思想

有限元素方法是基于变分原理，既通过求解一个泛函取极小值的变分问题。

有限元素法是在变分原理的基础上吸收差分格式的思想发展起来的, 是变分问题中欧拉法的进一步发展。它是人们在尝试求解具有复杂区域, 复杂边界条件下的数学物理方程的过程中,找到的一种比较完美的离散化方法。

它比有限差分法的矩形网格划分方法在布局上更为合理，在处理复杂区域和复杂边界条件时更方便和适当。采用有限元素法还能使物理特性基本上被保持, 计算精度和收敛性进一步得到保证。正是由于有限元素法有这样一些优点, 尽管其计算格式比较复杂, 但仍然在很多场合代替了差分法。

不过需要指出的是：并不是所有有限差分法可以处理的问题都可以采用有限元素法。

泛函数

数学上，通常变量与变量间的关系称为函数，而泛函则是函数集合的函数。

例如，静电场的势函数ϕ(r)是定义在坐标空间的函数集，系统电场总能量 (ϕ(rU ))则是定义在该函数集中的一个泛函，可记为 (ϕ(rI ))。

离散化

一个连续的介质会被离散成数个简单的基本几何单元，即元素。这些元素通过节点互相联系。这个过程就是网格划分。

特点

（1） 能处理复杂区域和复杂边界条件的求解问题。

（2） 有限元素法是一套求解微分方程的系统化数值计算方法。

它比传统解法具有理论完整可靠，物理意义直观明确，解题效能强等优点。

（3） 对连续体的问题采用有限元素法, 是将连续问题离散化的数值求解方法。

应用范围

由于这种方法适应性强，形式单纯、规范，因而自50年代以来, 在计算机的配合下，有限元素法在物理和工程设计计算的许多领域得到了广泛的应用。该方法不仅适用于电磁场问题的求解, 也是对其它具有复杂边值问题的数学物理方程求解时的高效能的方法。

 

有限元素法的一般步骤

（1）    推导出与给定边界条件的偏微分方程等价的泛函表示。

（2）    把求解的区域用三角形元素划分为小的单元。然后对每个节点和三角形元素按照约定的规则分别进行编号。

（3）    利用公式

 

它的一般形式可以写为

 

由定义

 

公式

 

 

采用到如下三角形型函数的积分公式

 

可以得出

 

计算出各个三角形元素e的系数矩阵和 。

 

（4） 将各个三角形单元的系数矩阵 和

（5）装配成总矩阵(K)和(P)，形成有限元方程组，然后利用强加边界条件法对有限元方程组进行修正。

利用超松弛迭代法求解有限元方程组，则得到域内各个节点上的函数值。

 

 

有限元素法与有限差分法的比较

 数学方法的不同

有限元素法是基于数学上的变分原理，将所要求解的物理问题化为对泛函求极值的一个变分方程；再利用差分法中的区域划分的离散化方法，并通过元素划分所构造的插值函数，把求解连续的变分方程问题离散化为求解线性方程组。按照这样的有限元素方法来处理物理问题，就不再需要通过建立偏微分方程这一道步骤，并且其物理问题在离散化的整个过程中就始终具有明确的物理意义。

有限差分法来求解物理问题的数值解时，必须首先从物理模型出发，列出相应的偏微分方程及定解条件，然后通过网格划分将偏微分方程的求解问题离散化为对差分方程组的数值求解。

对区域的离散化方法不同

有限元素法采用的一般是三角形划分的方法。这样的划分对节点在区域内的配置方式比较任意，其配置方式可以根据边界条件的情况来选择。这样就可以在边界形状比较复杂时，仍然可以选 择边界节点完全处在区域的边界上，从而在边界上可以做到较好的逼近。特别是在由不同介电常数 的介质构成的静电场域内求解时，我们可以将节点取在不同介质区域的交界面上，并在电位梯度较大的区域，节点还可配置密一些，以实现较好的计算精度。一般在有限元素法中采用三角划分时，如果三角单元无限缩小，有限元近似解收敛于精确解，近似解的平均导数也收敛于精确解的导数。

在有限差分法中，通常采用的是矩形网络区域划分，因而很难实现网络节点在区域中的配置与边界（不同介质界面）的良好逼近。

计算精度比较

在有限元素法中，实践证明选择三角形不要太狭长，三角形元素边长比越接近 1，计算得到的数值质量越好。有限元素法的计算精度于三角形划分时最大三角形边长h有关，若精确解有二阶导数，则函数误差与 同阶，导数误差与h同阶（当h趋于无穷小时）。用有限元素法求解物理问题时，它是用统一的观点对区域内的节点和边界节点列出计算格式。这就使得各节点的计算精度总体上比较协调。此外，有限元素法的计算格式中的矩阵 (K)具有比较好的性质，即它是一个对称正定的大型稀疏矩阵。这就给求解有限元方程组带来方便。

有限差分法在采用直交网络时其计算精度与矩形最大边长h有关，此时列出的计算格式比有限元 素法简单方便。而有限差分法则是孤立地对微分方程及定解条件分别列差分方程，因而各节点精度总体上不够一致。

应用范围比较

在对边界形状规则的求解区域，自然采用有限差分法就比较合适。事实上，有限差分法的适用范围要比有限元素法广泛得多。有很多物理问题目前还不能用有限元素法求解，但是人们总是可以采用有限差分法。特别是在边界形状比较规则时，采用有限差分法是最合适的。

1.5      数学形式以及度计算机的要求

有限元素法的节点配置比较任意，计算格式要列出来就要复杂得多。但是这些计算格式都可以在电子计算机上自动形成，也容易将程序标准化，因而这并不会影响它的实际应用。但是有限元素法要求的计算机内存量比较大，需要准备输入的数据量也比较大，这是它的缺点之一。

 

---

写在最后：这是我之前做的一些技术初步调研，主要目的是了解和展示一些技术，并进行了一些个人总结，其中必然有疏漏。其中引用了不少论文，百科和他人博客，但因为时间太久大都忘了出处，如果侵权提示后删除。