首页 | 本学科首页   官方微博 | 高级检索  
     


Partition of unity for localization in implicit a posteriori finite element error control for linear elasticity
Authors:Carsten Carstensen  Jan Thiele
Affiliation:1. Department of Mathematics, Humboldt‐Universit?t zu Berlin, Unter den Linden 6, Berlin D‐10099, Germany;2. BMW Group, Hufelandstraβe 8a, München D‐80788, Germany
Abstract:The partition of unity for localization in adaptive finite element method (FEM) for elliptic partial differential equations has been proposed in Carstensen and Funken (SIAM J. Sci. Comput. 2000; 21 : 1465–1484) and is applied therein to the Laplace problem. A direct adaptation to linear elasticity in this paper yields a first estimator ηL based on patch‐oriented local‐weighted interface problems. The global Korn inequality with a constant CKorn yields reliability for any finite element approximation uh to the exact displacement u. In order to localize this inequality further and so to involve the global constant CKorn directly in the local computations, we deduce a new error estimator µL. The latter estimator is based on local‐weighted interface problems with rigid body motions (RBM) as a kernel and so leads to effective estimates only if RBM are included in the local FE test functions. Therefore, the excluded first‐order FEM has to be enlarged by RBM, which leads to a partition of unit method (PUM) with RBM, called P1+RBM or to second‐order FEMs, called P2 FEM. For P1+RBM and P2 FEM (or even higher‐order schemes) one obtains the sharper reliability estimate . Efficiency holds in the strict sense of . The local‐weighted interface problems behind the implicit error estimators ηL and µL are usually not exactly solvable and are rather approximated by some FEM on a refined mesh and/or with a higher‐order FEM. The computable approximations are shown to be reliable in the sense of . The oscillations are known functions of the given data and higher‐order terms if the data are smooth for first‐order FEM. The mathematical proofs are based on weighted Korn inequalities and inverse estimates combined with standard arguments. The numerical experiments for uniform and adapted FEM on benchmarks such as an L‐shape problem, Cook's membrane, or a slit problem validate the theoretical estimates and also concern numerical bounds for CKorn and the locking phenomena. Copyright © 2007 John Wiley & Sons, Ltd.
Keywords:partition of unity method  a posteriori error estimate  adaptive finite element method  Korn inequality  linear elasticity  incompressibility  locking
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号