Extended finite element method with edge-based strain smoothing (ESm-XFEM) for linear elastic crack growth

被引:200
作者
Chen, L. [2 ,3 ]
Rabczuk, T. [4 ]
Bordas, S. P. A. [1 ]
Liu, G. R. [5 ]
Zeng, K. Y. [3 ]
Kerfriden, P. [1 ]
机构
[1] Cardiff Univ, iMAM, Cardiff Sch Engn, Cardiff CF24 3AA, S Glam, Wales
[2] Inst High Performance Comp, Singapore 138632, Singapore
[3] Natl Univ Singapore, Dept Mech Engn, Singapore 117576, Singapore
[4] Bauhaus Univ Weimar, Inst Struct Mech, D-1599423 Weimar, Germany
[5] Univ Cincinnati, Cincinnati, OH 45221 USA
基金
英国工程与自然科学研究理事会;
关键词
Fracture analysis; Numerical method; Edge-based smoothed finite element method; Extended finite element method; Stress intensity factor; Convergence rate; CONFORMING NODAL INTEGRATION; PARTICLE METHODS; LEVEL SETS; FEM; PARTITION;
D O I
10.1016/j.cma.2011.08.013
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
This paper presents a strain smoothing procedure for the extended finite element method (XFEM). The resulting "edge-based" smoothed extended finite element method (ESm-XFEM) is tailored to linear elastic fracture mechanics and, in this context, to outperform the standard XFEM. In the XFEM, the displacement-based approximation is enriched by the Heaviside and asymptotic crack tip functions using the framework of partition of unity. This eliminates the need for the mesh alignment with the crack and re-meshing, as the crack evolves. Edge-based smoothing (ES) relies on a generalized smoothing operation over smoothing domains associated with edges of simplex meshes, and produces a softening effect leading to a close-to-exact stiffness, "super-convergence" and "ultra-accurate" solutions. The present method takes advantage of both the ES-FEM and the XFEM. Thanks to the use of strain smoothing, the subdivision of elements intersected by discontinuities and of integrating the (singular) derivatives of the approximation functions is suppressed via transforming interior integration into boundary integration. Numerical examples show that the proposed method improves significantly the accuracy of stress intensity factors and achieves a near optimal convergence rate in the energy norm even without geometrical enrichment or blending correction. (C) 2011 Elsevier B.V. All rights reserved.
引用
收藏
页码:250 / 265
页数:16
相关论文
共 66 条