有限元方法_原创文档.pdfVIP

有限元方法

求解微分方程，特别是椭圆型边值问题的一种离散化方法，其基础是变分原理和剖

分逼近。有限元方法是传统的里茨－加廖金方法的发展，并融会了差分法的优点，处理

上统一，适应能力强，已广泛应用于科学与工程中庞大复杂的计算问题。

作为有限元方法出发点的变分原理，是表达物理基本定律的一种普遍形式。其表述

可概括如下：给出一个依赖物理状态v的变量J(v)（v是函数,J(v)在数学上称为泛函），

同时给出J(v)的容许函数集V,即一切可能的物理状态,则真实的状态是V中使J(v)达到

极小值的函数。剖分逼近是有限元离散化的手段，把问题的整体（即求解域）剖分为有

限个基本限个基本块块，称为，称为单单元元，然后通过单元上的插值逼近，得到一个结构简单的函数集，

称为称为有限元有限元空间空间它一般是容许函数集V的子集或有某种联系。有限元方法就是在这个

有限元空间中寻找J(v)的极小解作为近似解。

典型问题为具体说明有限元方法，讨论二维有界域Ω上的椭圆型方程

,(1)

变系数β表示介质不均匀。物理学中许多平衡态或定常态问题都可归结为这个典型

方程。与方程(1)相配的有如下三类边界条件：

第一类：；

第二类：；

第三类：。这里的φ、g及α均为定义在边界дΩ上的已

知函数，表示外法向导数，第二类边界条件是第三类当α=0时的特例。

为说明有限元方法能统一处理复杂的情况，假定讨论的问题是混合边值，并且介质

有间断,即дΩ分成Г和Г两部分，分别有边界条件

01

,(2)

,(3)

-+

β（x,y)有间断线，把Ω分为Ω,Ω两部分，在间断线上微分方程(1)无定义，而代之

以接触条件

,(4)

+-

及表示间断线上分别指向Ω及Ω的法向导数。

变分原理与微分方程(1)及附加条件(2)、(3)、(4)的边值问题相对应的是物理学

中的极小能量原理。构造中的极小能量原理。构造能量积分能量积分

并取J(v)的容许函数集V为一切满足边界条件(2)且一阶偏导数平方可积的函数,则

使J(v)达到极小值的u,即

,(6)

也必满足方程(1)及(2)、(3)、(4)。事实上，极小能量原理之类的变分原理是物理问

题的原始形式，微分方程是数学推导的结果。在变分问题中,只有边界条件(2)是强加到

容许函数集上的,边界条件(3)及间断介质的接触条件(4)都是极小解u自然满足的，这种

情况有利于离散化的统一处理。

剖分逼近几何剖分的基本单元可取为三角形、矩形、四边形、曲边形等等，其

中三角形最基本常用。

假定问题的求解区域为多边形，介质间断线

为折线，作三角剖分如图所示。在剖分中需注意

介质间断线与某些三角形的边重合，不同类边界

条件的交点与某些三角形的顶点重合。单元的顶

点称为网格结点，在дΩ上称