二阶非自伴两点边值问题Galerkin有限元后处理超收敛解答计算的EEP法.pdfVIP

第24 卷第2期 计 算 力 学 学 报 Vol . 2 4 ，No . 2 2007 年4 月 Chinese Journal of Computational Mechanics April 2007 文章编号：1007-4708（2007 ）02-0 142-06 二阶非自伴两点边值问题Galerkin 有限元 后处理超收敛解答计算的EEP 法 * 袁 驷 ， 林永静 （清华大学 土木工程系，北京 100084 ） 摘 要：将一维Ritz 有限元法超收敛计算的EEP（单元能量投影）法推广到二阶非自伴常微分方程两点边值问 题Galerkin 有限元法的超收敛计算。在对精确单元的研究中，发现与Ritz 有限元法不同，只要检验函数采用伴 随算子方程的解，无论试函数取何形式，在结点处都可得到精确的解函数值。对近似单元的研究表明，EEP 法 同样适用于Galerkin 有限元法，不仅保留了简便易行、行之有效、效果显著的特点，同时也保留了EEP 法的特有 优点，如：任一点的导数和解函数的误差与结点值的误差具有相同的收敛阶。 关键词： Galerkin 有限元；非自伴问题；一维问题；超收敛；单元能量投影 中图分类号：0242 . 21 文献标识码：A 问题： 1 引 言 有限元法是一种非常有效且通用的数值分析 {Lu E- （p u ） + r u + g u =f ，0 x 1 （1a ） u（0 ）= 0 ， u （1 ）= 0 （1b ） 计算方法。在常规的位移型有限元法中，应力（导 式中p 0 ，r 、g 0 ，（） = c（）/ cx 。其中L 为线性 数）精度比位移（函数）精度呈数量级的下降，引发 非自伴随微分算子，其伴随算子L * 为 了超收敛计算的研究，成为近年来有限元领域研究 L *1 E- （p 1 ） - （r 1 ） + g 1 （2 ） ［1-3 ］ ［7 ］ 中的热点与难点之一 。 定义双线性形及线性形如下 ： 1 文献［4 ］对于二阶自伴随两点边值问题的Ritz a（u ，1 ）= J（p u 1 + r u 1 + g u1 ）cx 有限元法提出了一种简单高效的单元能量投影 0 1 EEP（Element Energy Projection ）法，可以给出任一 （f ，1 ）= Jf 1cx （3 ） 0 单元上任一点处的超收敛的位移和应力，并都能获