Highlights• We implement an isogeometric BEM routine for fracture by use of dual boundary integral equations.• We propose a singular integration scheme to improve the quadrature accuracy for elements with high aspect ratios.• We investigate the approaches to compute stress intensity factors based on a NURBS representation of the crack surfaces.• We outline a geometric algorithm to propagate the crack based on the fatigue Paris law.
AbstractWe 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.