[en] Numerical crack propagation schemes were augmented in an elegant manner by the X-FEM method. The use of special tip enrichment functions, as well as a discontinuous function along the sides of the crack allows one to do a complete crack analysis virtually without modifying the underlying mesh, which is of industrial interest, especially when a numerical model for crack propagation is desired. This paper improves the implementation of the X-FEM method for stress analysis around cracks in three ways. First, the enrichment strategy is revisited. The conventional approach uses a 'topological' enrichment (only the elements touching the front are enriched). We suggest a 'geometrical' enrichment in which a given domain size is enriched. The improvements obtained with this enrichment are discussed. Second, the conditioning of the X-FEM both for topological and geometrical enrichments is studied. A preconditioner is introduced so that 'off the shelf' iterative solver packages can be used and perform as well on X-FEM matrices as on standard FEM matrices. The preconditioner uses a local (nodal) Cholesky based decomposition. Third, the numerical integration scheme to build the X-FEM stiffness matrix is dramatically improved for tip enrichment functions by the use of an ad hoc integration scheme. A 2D benchmark problem is designed to show the improvements and the robustness. Copyright (c) 2005 John Wiley & Sons, Ltd.
Disciplines :
Engineering, computing & technology: Multidisciplinary, general & others Mathematics
Author, co-author :
Béchet, Eric ; Université de Liège - ULiège > Département d'aérospatiale et mécanique > Conception géométrique assistée par ordinateur
Minnebol, H.; Institut de Recherche en Génie Civil et Mécanique, UMR CNRS 6183, Ecole Centrale de Nantes, B.P. 92101, 44321 Nantes Cedex 3, France
Moes, N.; Institut de Recherche en Génie Civil et Mécanique, UMR CNRS 6183, Ecole Centrale de Nantes, B.P. 92101, 44321 Nantes Cedex 3, France
Burgardt, B.; SNECMA Moteurs, Rond Point René Ravaud—Réau, 77550 Moissy-Cramayel, France
Language :
English
Title :
Improved implementation and robustness study of the X-FEM for stress analysis around cracks
Publication date :
2005
Journal title :
International Journal for Numerical Methods in Engineering
ISSN :
0029-5981
eISSN :
1097-0207
Publisher :
John Wiley & Sons, Inc, Chichester, United Kingdom
Barsoum RS. Application of quadratic isoparametric elements in linear fracture mechanics. International Journal of Fracture 1974; 10:603-605.
Babuška I, Melenk JM. The partition of unity method. International Journal for Numerical Methods in Engineering 1997; 40:727-758.
Strouboulis T, Babuška I, Copps K. The design and analysis of the generalized finite element method. Computer Methods in Applied Mechanics and Engineering 2000; 181:43-69.
Moës N, Dolbow J, Belytschko T. A finite element method for crack growth without remeshing. International Journal for Numerical Methods in Engineering 1999; 46:131-150.
Belytschko T, Black T. Elastic crack growth in finite elements with minimal remeshing. International Journal for Numerical Methods in Engineering 1999; 44:601-620.
Strouboulis T, Copps K, Babuška I. The generalized finite element method: an example of its implementation and illustration of its performance. International Journal for Numerical Methods in Engineering 2000; 47:1401-1417.
Berntson J, Espelid TO, Genz A. An adaptive algorithm for the approximate calculation of multiple integrals. ACM Transactions on Mathematical Software 1991; 17:437-451.
Destuynder P, Djaoua M, Lescure S. Some remarks on elastic fracture mechanics. Journal de Mécanique Théorique et Appliquée 1983; 2:113-135.
Shih CF, Moran B, Nakamura T. Energy release rate along a three-dimensional crack front in a thermally stressed body. International Journal of Fracture 1986; 30:79-102.
Rice JR. A path independent integral and the approximate analysis of strain concentrations by notches and cracks. Journal of Applied Mechanics 1968; 35:379-386.
Moës N, Gravouil A, Belytschko T. Non-planar 3D crack growth by the extended finite element and level sets - Part I: mechanical model. International Journal for Numerical Methods in Engineering 2002; 53:2549-2568.
Gosz M, Moran B. An interaction energy integral method for computation of mixed-mode stress intensity factors along non-planar crack fronts in three dimensions. Engineering Fracture Mechanics 2002; 69:299-319.
Bui HD. Mécanique de la rupture fragile (in French). Masson: Paris, 1978.
Chessa J, Wang H, Belytschko T. On the construction of blending elements for local partition of unity enriched finite elements. International Journal for Numerical Methods in Engineering 2003; 57:1015-1038.
Burgardt B. Contribution à l'étude des méthodes des équations intégrales et à leur mise en oeuvre numérique en élastoplasticité (in French). Thèse de doctorat, Université de Nantes, 1999.
Moës N, Cloirec M, Cartraud P, Remacle J-F. A computational approach to handle complex microstructure geometries. Computer Methods in Applied Mechanics and Engineering 2003; 192:3163-3177.
Gravouil A, Moës N, Belytschko T. Non-planar 3D crack growth by the extended finite element and level sets - Part II: level set update. International Journal for Numerical Methods in Engineering 2002; 53:2569-2586.
Laborde P, Pommier J, Renard Y, Salaün M. High order extended finite element method for cracked domains. International Journal for Numerical Methods in Engineering 2004, submitted.
Béchet E, Minnebo H, Moës N, Burgardt B. Convergence and conditioning issues with X-FEM in fracture mechanics. WCCM6, 5-10 September 2004, Beijing.