Self-Adaptive Strategy for One-Dimensional Finite Element Method Based on EEP Method
-
摘要: 基于新近提出的一维有限元后处理超收敛算法——单元能量投影（EEP）法，将有限元自适应求解问题转化为对超收敛解答的自适应分段多项式插值问题；对于大多数问题，一步便可获得满意的有限元网格划分，在该网格上再次进行有限元计算，一般即可获得满足用户给定的误差限的有限元解答．即便未能完全满足精度要求，一般只需局部细分加密网格一至二步即可．该法简单实用、高效可靠，是一个颇具优势和潜力的自适应方法．以二阶椭圆型常微分方程模型问题为例，对该法的基本思想、实施策略及具体算法做一介绍，并给出有代表性的数值算例用以展示该法的优良性能和效果．Abstract: Based on the newly-developed element energy projection(EEP)method for computation of super-convergent results in one-dimensional finite element method(FEM),the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation.As a result, a satisfactory FEM mesh can be obtained,and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance.Even though the error tolerance was not completely satisfied,one or two steps of further local refinements would be sufficient.This strategy has been found to be very simple,rapid,cheap and efficient.Taking the elliptical ordinary differential equation of the second order as the model problem,the fundamental idea, implementation strategy and detailed algorithm were described.Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.
-
Key words:
- FEM /
- self-adaptive solution /
- super-convergence /
- element energy projection /
- ODE
-
[1] Babuska I,Rheinboldt W C.A posteriori error analysis of finite element method for one-dimensional problems[J].SIAM Journal on Numerical Analysis,1981,18(3):565—589. doi: 10.1137/0718036 [2] Zienkiewicz O C, Zhu J Z. The superconvergence patch recovery (SPR) and a posteriori error estimates,part 1: the recovery technique[J].Internat J Numer Methods Engrg,1992,33(7):1331—1364. doi: 10.1002/nme.1620330702 [3] Zienkiewicz O C, Zhu J Z.The superconvergence patch recovery (SPR) and a posteriori error estimates,part 2: error estimates and adaptivity[J].Internat J Numer Methods Engrg,1992,33(7):1365—1382. doi: 10.1002/nme.1620330703 [4] 林群，朱起定.有限元的预处理和后处理理论[M].上海：上海科学技术出版社，1994. [5] 陈传淼. 有限元超收敛构造理论[M].长沙：湖南科学技术出版社，2002. [6] Ascher U, Christiansen J,Russell R D.Algorithm 569, COLSYS: Collocation software for boundary value ODEs[D2]. ACM Trans Math Software,1981,7(2):223—229. [7] YUAN Si.The Finite Element Method of Lines[M].Beijing,New York: Science Press, 1993. [8] 袁驷，王枚，林永静，等. 有限元(线)法超收敛应力计算的新思路[J].工程力学(增刊). 长沙：全国结构工程大会，2002，112—118. [9] 袁驷，王枚.一维有限元后处理超收敛解答计算的EEP法[J]. 工程力学，2004，21(2):1—9. [10] Strang G, Fix G.An Analysis of the Finite Element Method[M].London: Prentice-Hall, 1973. [11] Douglas J, Dupont T. Galerkin approximations for the two point boundary problems using continuous piecewise polynomial spaces[J].Numer Math，1974,22(2):99—109. doi: 10.1007/BF01436724 -
点击查看大图
计量
- 文章访问数: 2523
- HTML全文浏览量: 64
- PDF下载量: 628
- 被引次数: 0
下载:
