Isogeometric boundary element methods for three dimensional static fracture and fatigue crack growth

被引:194
作者
Peng, X. [1 ]
Atroshchenko, E. [2 ]
Kerfriden, P. [1 ]
Bordas, S. P. A. [1 ,3 ,4 ,5 ]
机构
[1] Cardiff Univ, Inst Mech Mat & Adv Mfg, Cardiff CF24 3AA, S Glam, Wales
[2] Univ Chile, Dept Mech Engn, Santiago 8370448, Chile
[3] Univ Luxembourg, Fac Sci Technol & Commun, 6 Rue Richard Coudenhove Kalergi, L-1359 Luxembourg, Luxembourg
[4] Res Unit Engn Sci, Campus Kirchberg G 007, Luxembourg, Luxembourg
[5] Univ Western Australia, 35 Stirling Hwy, Crawley, WA 6009, Australia
基金
英国工程与自然科学研究理事会; 欧洲研究理事会;
关键词
Isogeometric analysis; NURBS; Linear elastic fracture; Boundary element method; Crack growth; Fracture simulations directly from CAD; EXTENDED FINITE-ELEMENT; STRESS INTENSITY FACTORS; SHAPE OPTIMIZATION; INTEGRAL METHOD; LEVEL SETS; COMPUTATIONAL DOMAIN; METHOD XIBEM; PART I; 3D; PROPAGATION;
D O I
10.1016/j.cma.2016.05.038
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
We present a novel numerical method to simulate crack growth in 3D, directly from the Computer-Aided Design (CAD) geometry of the component, without any mesh generation. The method is an isogeometric boundary element method (IGABEM) based on non-uniform rational B-splines (NURBS). NURBS basis functions are used for the domain and crack representation as well as to approximate the physical quantities involved in the simulations. A stable quadrature scheme for singular integration is proposed to enhance the robustness of the method in dealing with highly distorted elements. Convergence studies in the crack opening displacement is performed for a "penny-shaped crack and an elliptical crack. Two approaches to extract stress intensity factors (SIFs): the contour M integral and the virtual crack closure integral are compared using dual integral equations. The results show remarkable accuracy in the computed SIFs, leading to smooth crack paths and reliable fatigue lives, without requiring the generation of any mesh from the CAD model of the component under consideration. (C) 2016 Elsevier B.V. All rights reserved.
引用
收藏
页码:151 / 185
页数:35
相关论文
共 128 条
[1]   A well-conditioned and optimally convergent XFEM for 3D linear elastic fracture [J].
Agathos, Konstantinos ;
Chatzi, Eleni ;
Bordas, Stephane P. A. ;
Talaslidis, Demosthenes .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2016, 105 (09) :643-677
[2]  
Aimi A., 2015, APPL MATH COMPUT
[3]  
[Anonymous], 2008, ANAL FINITE ELEMENT
[4]   THE POST-PROCESSING APPROACH IN THE FINITE-ELEMENT METHOD .2. THE CALCULATION OF STRESS INTENSITY FACTORS [J].
BABUSKA, I ;
MILLER, A .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 1984, 20 (06) :1111-1129
[5]   FURTHER APPLICATION OF QUADRATIC ISOPARAMETRIC FINITE-ELEMENTS TO LINEAR FRACTURE MECHANICS OF PLATE BENDING AND GENERAL SHELLS [J].
BARSOUM, RS .
INTERNATIONAL JOURNAL OF FRACTURE, 1975, 11 (01) :167-169
[6]   Numerical schemes for the Hamilton-Jacobi and level set equations on triangulated domains [J].
Barth, TJ ;
Sethian, JA .
JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 145 (01) :1-40
[7]   A simple approach to the numerical simulation with trimmed CAD surfaces [J].
Beer, G. ;
Marussig, B. ;
Zechner, J. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2015, 285 :776-790
[8]   Enriched finite elements and level sets for damage tolerance assessment of complex structures [J].
Bordas, S ;
Moran, B .
ENGINEERING FRACTURE MECHANICS, 2006, 73 (09) :1176-1201
[9]   Three-dimensional crack initiation, propagation, branching and junction in non-linear materials by an extended meshfree method without asymptotic enrichment [J].
Bordas, Stephane ;
Rabczuk, Timon ;
Zi, Goangseup .
ENGINEERING FRACTURE MECHANICS, 2008, 75 (05) :943-960
[10]   Derivative recovery and a posteriori error estimate for extended finite elements [J].
Bordas, Stephane ;
Duflot, Marc .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2007, 196 (35-36) :3381-3399