Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2013 (2013), Article ID 721656, 14 pages
Research Article

Two-Dimensional Stress Intensity Factor Analysis of Cracks in Anisotropic Bimaterial

1Department of Resources Engineering, National Cheng Kung University, No. 1 University Road, Tainan City 70101, Taiwan
2Graduate Institute of Applied Geology, National Central University, No. 300 Jhongda Road, Jhongli City 32001, Taiwan
3Graduate Institute of Geophysics, National Central University, No. 300 Jhongda Road, Jhongli City 32001, Taiwan
4Geotechnical Engineering Research Center, Sinotech Engineering Consultants, Inc., Taipei 110, Taiwan

Received 11 June 2012; Revised 3 November 2012; Accepted 9 November 2012

Academic Editor: Giuseppe Carbone

Copyright © 2013 Chia-Huei Tu 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 paper presents a 2D numerical technique based on the boundary element method (BEM) for the analysis of linear elastic fracture mechanics (LEFM) problems on stress intensity factors (SIFs) involving anisotropic bimaterials. The most outstanding feature of this analysis is that it is a singledomain method, yet it is very accurate, efficient, and versatile (i.e., the material properties of the medium can be anisotropic as well as isotropic). A computer program using the BEM formula translation (FORTRAN 90) code was developed to effectively calculate the stress intensity factors (SIFs) in an anisotropic bi-material. This BEM program has been verified and showed good accuracy compared with the previous studies. Numerical examples of stress intensity factor calculation for a straight crack with various locations in both finite and infinite bimaterials are presented. It was found that very accurate results can be obtained using the proposed method, even with relatively simple discretization. The results of the numerical analysis also show that material anisotropy can greatly affect the stress intensity factor.