4有限元外推方法.doc

4 有限元外推方法 提高有限元方法的计算精度有多种方法，主要有以下几种： (1) 加密剖分方法； (2) 采用高阶有限方法； (3) 样条插值有限元方法； (4) 棱边有限元法； (5) 互补变分方法； (6) 有限元外推方法。 从简单方便，易学易用和实际应用效果来看，有限元外推法是一种值得推广使用的方法。 有限元外推方法基于有限元的解答的误差分析。从有限元误差分析的理论和分析方法比较复杂，是中科院数学研究所的林群老师和他的博士生们的研究成果，对我们来说主要是学习研究的结果、讨论它的应用，是如何理解定义和应用数学定理的问题。从原理上讲，它是一种优化方法，从线性有限元解答的搭配比例来实现优化。 下面简单阐述它的基本定义和定理，也是我们要用到的部分。 4.1 外推法的概念 我们用某种算法求解某一问题，当计算了步或计算了次之后，得到了个点或个逼近的数据，，，…，，如果根据这个已得到的序列，经过某中组合运算，可得到对原问题更精确的解答，我们就称之为外推法，这是一种数值优化算法。现在以圆周率的计算为例来说明。 如果以表示内接于一个半径为1的圆的正边形一边的边长，那么圆周率的近似值可按下式计算： 下面表中列出了以正六边形开始，边数不断加倍后得到的圆周率的近似值以及精确值： n (边数) 6 1 3 12 03.1058285 24 03.1322856 48 03.1393502 96 0.065438166 3.1410319 192 0.032723463 3.1414524 …… …… …… 精确值 3.141592654 不断的成倍增加正边形边数，可以得到不断增加精确度的圆周率。如果利用相邻两个近似值和，作以下的线性组合 来作为计算圆周率的近似值公式，例如：取正6边形和正12边形的相应圆周率近似值计算，有： 这相当于正96边形计算所达到的精度。如果取正96边形和正192边形对应之和值来作组合计算，可得： 小数点后第七位与圆周率的精确值相吻合，它比的近似值要精确得多。 由以上的简单例子可以看出，经过适当的组合加工，可将较粗造的近似值通过上面的线性组合式“外推”出精确度高得多的近似值，这种计算方法称为外推法。 4.2 有限元外推法 有限元外推法跟以上简单例子相似，也是从较简单的、剖分单元数较少的线性有限元出发，通过组合的外推计算，得到计算精度相当于剖分单元数多得多的或用高阶单元的复杂的有限元计算所得结果。 有限元外推法是我们国家首先开始研究，并取得了一定的成果，首先从二维场的研究开始。 图1 定义1：对于有限元三角形网络剖分，若其中任意两个相邻三角单元构成——近似平行四边形，则称该有限元网络剖分是强正规的 其中h为单元的线度，为与无关的正常数。 这一强正规要求的条件是比较高，它在剖分中有较强的局限性。 为扩展它的应用需要有如下拓宽的条件如下定义： 定义2：对于给定区域，可用光滑曲线将其划分为个无重叠的曲边三角形子区域，使。在每一个分区域上的有限元三角网格剖分，使满足强正规条件，则称在区域上网格部分是分片强正规的。 图2 基于上面的定义，提出关于线性有限元的外推定理。 定理1：设二阶线性椭圆型边值问题的解满足，在凸区域上的有限元三角网格剖分是强正规或分片强正规的。将域进行粗剖分，对应粗剖分的各单元边的中点连接起来得到加密部分，如图2所示。记和分别为和的节点集，和分别表示对应剖分和的线性有限元数值解，并设，。则第一类边值问题三角形线性有限元的解适合于如下的展式： (1) (2) 式中，系数的选择对电磁场问题来说，可取为2。 定理2：设二阶线性椭圆边值问题的解，对凸求解区域的有限元剖分是强正规或分片强正规的，则存在一个与无关的函数，使线性有线元解具有渐近误差展式 式中与分别表示线性有限元解和的插值函数，为中任一单元中远离角点的节点坐标。 由以上定理可知：对于二阶线性椭圆边值问题的解，其线性有限元解的误差展式反映出它的解具有与其单元尺寸相当的1阶收敛形式，而按定理1所示出的展式是比加密剖分还高1阶的解的收敛形式。即是说可用定理1的(1)、(2)两式，去掉误差项就成为以粗分和细分的两种有限元解直接推求更高精度解的计算式，可称之为递推公式。 对所计算的区域进行三角形粗部分和加密剖分，若所给边值问题的线性有限元解分别为 (3a) (3b