第五讲有限元素方法.PDFVIP

第五章 有限元素方法 5.1 有限元素方法的基本思想 特点： （1） 能处理复杂区域和复杂边界条件的求解问题。 （2 ） 有限元素法是一套求解微分方程的系统化数值计算方法。 它比传统解法具有理论完整可靠，物理意义直观明确，解题效能强等优点。 （3 ） 对连续体的问题采用有限元素法, 是将连续问题离散化的数值求解方法。 应用范围： 由于这种方法适应性强，形式单纯、规范，因而自50 年代以来, 在计算机的配合下， 有限元素法在物理和工程设计计算的许多领域得到了广泛的应用。该方法不仅适用于 电磁场问题的求解, 也是对其它具有复杂边值问题的数学物理方程求解时的高效能的 方法。 基本思想： 有限元素方法是基于变分原理，既通过求解一个泛函取极小值的变分问题。 有限元素法是在变分原理的基础上吸收差分格式的思想发展起来的, 是变分问题中欧拉法的 进一步发展。它是人们在尝试求解具有复杂区域, 复杂边界条件下的数学物理方程的过程中, 找到的一种比较完美的离散化方法。 它比有限差分法的矩形网格划分方法在布局上更为合理，在处理复杂区域和复杂边界条件时 更方便和适当。采用有限元素法还能使物理特性基本上被保持, 计算精度和收敛性进一步得 到保证。正是由于有限元素法有这样一些优点, 尽管其计算格式比较复杂, 但仍然在很多场合 代替了差分法而受到计算物理工作者的偏爱。 不过需要指出的是：并不是所有有限差分法可以处理的问题都可以采用有限元素法。 泛函： 数学上，通常变量与变量间的关系称为函数，而泛函则是函数集合的函数，也就是函数的函数。 ( ( )) U r ϕ 例如，静电场的势函数ϕ(r)是定义在坐标空间的函数集，系统电场总能量 则是定义在该函 ( ( )) I r ϕ 。 数集中的一个泛函，可记为 例子： 为了进一步说明有限元素方法的基本思想，我们考虑一个确定静电势的问题，该场域的 r 介质中放置了一个球形金属导体，球形金属导体的半径为 , 球外距离球中心 r 处的电位为ϕ(r) 。 0 当这个系统处在电荷平衡的状态下时，金属导体上的电荷分布应当是均匀的, 导体表面是等电 位的。我们按照通常的作法，把从导体表面到无穷远处的球面之间的空间, 作为导体外的全空间。 假定在这个导体外的空间中的体电荷密度到处为零。 则在此空间中的能量为 2 2 +∞ +∞ +∞ ε ε ∂ϕ ∂ϕ 2 ⎛ ⎞ 2 ⎛ ⎞ 2 (U) E dV 4 2r dr r dr