Mathematical Problems in Engineering
Volume 2012 (2012), Article ID 581493, 26 pages
Research Article

Two-Dimensional Fracture Mechanics Analysis Using a Single-Domain Boundary Element Method

1Geotechnical Engineering Research Center, Sinotech Engineering Consultants, Inc., Taipei 110, Taiwan
2Department of Business Administration, Fu Jen Catholic University, New Taipei City 242, Taiwan
3Department of Resources Engineering, National Cheng Kung University, Tainan 701, Taiwan

Received 13 March 2012; Accepted 21 March 2012

Academic Editor: Kue-Hong Chen

Copyright © 2012 Chien-Chung Ke et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


This work calculates the stress intensity factors (SIFs) at the crack tips, predicts the crack initiation angles, and simulates the crack propagation path in the two-dimensional cracked anisotropic materials using the single-domain boundary element method (BEM) combined with maximum circumferential stress criterion. The BEM formulation, based on the relative displacements of the crack tip, is used to determine the mixed-mode SIFs and simulate the crack propagation behavior. Numerical examples of the application of the formulation for different crack inclination angles, crack lengths, degree of material anisotropy, and crack types are presented. Furthermore, the propagation path in Cracked Straight Through Brazilian Disc (CSTBD) specimen is numerically predicted and the results of numerical and experimental data compared with the actual laboratory observations. Good agreement is found between the two approaches. The proposed BEM formulation is therefore suitable to simulate the process of crack propagation. Additionally, the anisotropic rock slope failure initiated by the tensile crack can also be analyzed by the proposed crack propagation simulation technique.