Modeling crack propagation path of anisotropic rocks using boundary element method

Chien Chung Ke, Chao Shi Chen, Cheng Yu Ku, Chih Hao Chen

Research output: Contribution to journalArticlepeer-review

15 Citations (Scopus)


In a cracked material, the stress intensity factors (SIFs) at the crack tips, which govern the crack propagatio and are associated with the strength of the material, are strongly affected by the crack inclination angle and the orientation with respect to the principal direction of anisotropy. In this paper, a formulation of the boundary element method (BEM), based on the relative displacements of the crack tip, is used to determine the mixed-mode SIFs of isotropic and anisotropic rocks. Numerical examples of the application of the formulation for different crack inclination angles, crack lengths, and degree of material anisotropy are presented. Furthermore, the BEM formulation combined with the maximum circumferential stress criterion is adopted to predict the crack initiation angles and simulate the crack propagation paths. The propagation path in cracked straight through Brazilian disc 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.

Original languageEnglish
Pages (from-to)1227-1253
Number of pages27
JournalInternational Journal for Numerical and Analytical Methods in Geomechanics
Issue number9
Publication statusPublished - 2009 Jun 25

All Science Journal Classification (ASJC) codes

  • Computational Mechanics
  • Materials Science(all)
  • Geotechnical Engineering and Engineering Geology
  • Mechanics of Materials


Dive into the research topics of 'Modeling crack propagation path of anisotropic rocks using boundary element method'. Together they form a unique fingerprint.

Cite this