`Mathematical Problems in EngineeringVolume 2010, Article ID 947818, 22 pageshttp://dx.doi.org/10.1155/2010/947818`
Research Article

## Analysis of Crack Propagation Path on the Anisotropic Bi-Material Rock

Department of Resources Engineering, National Cheng Kung University, Tainan, Taiwan

Received 29 July 2010; Accepted 5 September 2010

Copyright © 2010 Chao-Shi Chen 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.

#### Abstract

This paper presents a single-domain boundary element method (SDBEM) for linear elastic fracture mechanics analysis in the 2D anisotropic bimaterial. In this formulation, the displacement integral equation is collocated on the uncracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. The complete fundamental solution (Green's function) for anisotropic bi-materials was also derived and implemented into the boundary integral formulation so the discretization along the interface can be avoided except for the interfacial crack part. A special crack-tip element was introduced to capture exactly the crack-tip behavior. A computer program with the FORTRAN code has been developed to effectively calculate the stress intensity factors, crack initiation angle, and propagation path of an anisotropic bi-material. This SDBEM program has been verified having a good accuracy with the previous researches. In addition, a rock of type (1)/(2) disk specimen with a central crack was made to conduct the Brazilian test under diametrical loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation.

#### 1. Introduction

In rock masses, the interbed construction suffers from cracking which is caused by various factors. Of greater concern are those cracks that develop as a result of initiation and propagation path, leading to significant change in the failure resistance of the structure. Many relative researches and discussions are started like wildfire and never stopped.

Because of the discontinuities of weak interbed on laccoliths (rock mass), there are many damages that occurred which are due to earth stress efforts or under geotechnical engineering, such as slope slip, tunnel collapse, deep excavation collapse, and crack openings which are due to deep well drilling. Crack is one of the fracture models to cause those damages because the crack openings and propagation on the field affect the rock mass structure stability. This paper will discuss the behavior of crack propagation on the basis of the theory of fracture mechanics. We defined the rock sample on interbed as “Bimaterial rock”.

Theoretically, interfacial crack problems in isotropic Bimaterials were studied [13] where the authors showed that the stresses possess the singularity of . Rice [4] re-examined the elastic fracture mechanics concepts for the isotropic interfacial crack and introduced an intrinsic material length scale so that the definition of the stress intensity factor (SIFs) possesses the same physical significance as those for the homogeneous cracks. Clements [5] and Willis [6] studied interfacial crack problems in dissimilar anisotropic materials. They showed that the oscillatory behavior of the stresses and the phenomenon of interpenetrating of the crack faces were also present near the crack tip for anisotropic interface cracks. Recent studies on interfacial cracks in anisotropic materials have been conducted by many authors [721], and different definitions for the stress intensity factor exist. By introducing a characteristic length, however, the definition given by Gao et al. [15], and Wu [12, 13] is consistent with Rice’s general definition [4] and appears to be more explicit than other definitions.

The study of fracture initiation and propagation in anisotropic rocks is subjected to Brazilian loads. A numerical procedure based on the SDBEM and maximum tensile stress criterion has been developed to predict the angle of crack initiation and the path of crack propagation in anisotropic rocks. Crack propagation in an anisotropic homogeneous rock disc under mixed mode I-II loading is simulated by an incremental crack extension with a piecewise linear discretization. A computer program, which can automatically generate a new mesh, has been developed to simulate the crack propagation process. Some experimental observations of crack initiation angles and crack propagation were obtained by conducting diametrical loading of initially cracked discs of a gypsum/cement. It was found that the numerical analysis could predict relatively well the direction of crack initiation and the crack propagation path.

#### 2. Theoretical Background

##### 2.1. Green’s Function in Anisotropic Bimaterials

With three complex analytical functions , one can, in general, express displacements, stresses, and tractions as follows [11, 22, 23]: where , denotes the real part of a complex variable or function, a prime denotes the derivative, the three complex numbers and the elements of the complex matrices B and A are functions of the elastic properties [11, 22, 23].

Assume that the medium is composed of two joined dissimilar anisotropic and elastic half-planes. We let the interface be along the x-axis and the upper and lower half-planes occupied by materials () and (), respectively (Figure 1).

Figure 1: Definition of the coordinate systems within an anisotropic Bimaterial.

For concentrated force acting at the point () in material () (), we express the complex vector function as [11] where the vector function with the argument having the generic form .

In (2.2), is a singular solution corresponding to a point force acting at the point () in an anisotropic infinite plane with the elastic properties of material (). This singular solution can be expressed as [11, 23] where , is the point force vector, and H is given by

There are two unknown vector functions to be solved in (2.2), that is, and . While the former is analytic in the upper (material ()) half-plane, the latter is analytic in the lower (material ()) half-plane. These expressions can be found by requiring continuity of the resultant traction and displacement across the interface, along with the standard analytic continuation arguments. Following this approach and after some complex algebraic manipulation, the complex vector functions in materials () and () are obtained as

In (2.6), the special subscripts () and () are used exclusively to denote that the corresponding matrix or vector is in material () () and material () (), respectively.

Similarly, for a point force in material (2.1) (), these complex functions can be found as where the vector functions are the infinite plane solution given in (2.4) but with the elastic properties of material ().

With the complex function given in (2.6) and (2.7), Green’s functions of the displacement and traction can be obtained by substituting these complex functions into (2.1). These Green’s functions have four different forms depending on the relative location of the field and source points.

It is noteworthy that these Green’s function can be used to solve both plane stress and plane strain problems in anisotropic Bimaterials. Although the isotropic solution cannot be analytically reduced from these Green’s functions one can numerically approximate it by selecting a very weak anisotropic (or nearly isotropic) medium [24, 25].

##### 2.2. BEM Formulation for 2D Cracked Anisotropic Bimaterials

In this section, we present an SDBEM formulation in which neither the artificial boundary nor the discretization along the un-cracked interface is necessary. This SDBEM formulation was widely used recently by Chen et al. [26], for homogeneous materials and is now extended to anisotropic Bimaterials.

The displacement integral equation applied to the outer boundary results in the following form ( only, Figure 2): where = 1,2, and are Green’s tractions and displacements, and are the boundary displacement and tractions, respectively, are quantities that depend on the geometry of the boundary and are equal to for a smooth boundary, and and are the field points and the source points on the boundary of the domain. has the same outward normal as . Here, the subscripts and denote the outer boundary and the crack surface, respectively.

Figure 2: Geometry of a two-dimensional cracked domain.

The traction integral equation (for being a smooth point on the crack) applied to one side of the crack surfaces is ( only) where is the outward normal at the crack surface and is the fourth-order stiffness tensor.

Equations (2.8) and (2.9) form a pair of boundary integral equations [24, 27, 28] and can be used for the calculation of SIFs in anisotropic Bimaterials. The main feature of the BEM formulation is that it is a single-domain formulation with displacement integral (2.8) being collocated on the un-cracked boundary only and traction integral (2.9) on one side of the crack surface only. For problems without cracks, one needs (2.8) only, with the integral on the crack surface being discarded. Equation (2.8) then reduces to the well-known displacement integral on the un-cracked boundary being discarded. We emphasize here that since bimaterial Green’s functions are included in (2.8), discretization along the interface can be avoided, with the exception of the interfacial crack part which will be treated by the traction integral equation presented by (2.9).

It is well known that a cracked domain poses certain difficulties for BEM modeling (Cruse, 1988 [29]). Previously, fracture mechanics problems in isotropic or anisotropic bimaterials were mostly handled by the multidomain method in which each side of the crack surface is put into different domains and artificial boundaries are introduced to connect the crack surface to the un-cracked boundary. For the bimaterial case, discretization along the interface is also required if one uses the Kelvin-type (infinite domain) Green’s functions.

##### 2.3. Crack-Tip Modeling

In fracture mechanics analysis, especially in the calculation of the SIFs, one needs to know the asymptotic behavior of the displacements and stresses near the crack-tip. In our BEM analysis of the SIFs, we propose to use the extrapolation method of the crack tip displacements. We therefore need to know the exact asymptotic behavior of the relative crack displacement behind the cracktip. This asymptotic expression has different forms depending on the location of the cracktip. In this paper, two cases will be discussed, that is, a crack-tip within the homogeneous material and an interfacial crack-tip. Inclined cracks terminating at the interface will be discussed in a future paper.

###### 2.3.1. A Crack Tip within a Homogeneous Material

The mixed mode stress intensity factors (SIFs) for anisotropic media can be determined by using the extrapolation method of the relative crack displacement (RCD), combined with a set of the shape functions. The RCD is defined as [24]

For this case, the relation of the RCDs at a distance behind the crack tip and the SIFs can be found as [24, 25, 30] where

Substituting the RCDs into (2.10) and (2.11), we obtain a set of algebraic equations in which the SIFs and can be solved.

###### 2.3.2. An Interfacial Crack Tip

For this case, the relative crack displacements at a distance behind the interfacial crack tip can be expressed, in terms of the three SIFs, as [15] where , and are the relative parameters in material () and material (). Utilizing (2.5), we defined the matrix of material as where and are two real matrices, and then utilizing these two matrices, we define matrix as and the characteristic relative to material Then, we used the characteristic obtain to define oscillation index as where is a identity matrix.

Utilizing the relationship between the characteristic and oscillation index , we define constant as where is the characteristic distance along material interface to crack tip.

Comparing (2.10) with (2.13), we noticed that while the relative crack displacement behaves as a square root for a crack tip within a homogeneous medium, for an interfacial crack tip, its behavior is , a square-root feature multiplied by weak oscillatory behaviors.

Equation (2.13) can be recast into the following form, which is more convenient for the current numerical applications: where is the characteristic length and is a matrix function with its expression given by Again, in order to capture the square-root and the weak oscillatory behavior, we construct a crack-tip element with tip at in terms of which the relative crack displacement is expressed as

In meshing the 2D anisotropic bimaterial problem (as shown in Figure 3), we assume that the interface is along the X-axis and the upper () and lower () half-planes are occupied by materials () and (), respectively. The corner of outer boundary is processed by the discontinuous elements Type I and Type III; the continuous element Type V is to deal with all outer smooth boundary; internal crack surface is processed by crack surface elements Type II; and crack tip element Type VI is to process crack tip problem. In order to avoid the oscillatory behavior of the interface, we mesh an anisotropic problem (as shown in Figure 3()) of bimaterial in which neither the interfacial elements nor the discretization along the un-cracked interface is necessary, with the exception of the interfacial crack part which will be treated by traction integral Equation (2.9) and (2.10) of RCD.

Figure 3: Mesh 2D bimaterial problem with the five types of quadratic elements.
##### 2.4. Crack Initiation and Fracture Propagation

In fracture mechanics, there are three criteria commonly used to predict the crack initiation angle: the maximum tensile stress criterion, or -criterion [31], the maximum energy release rate criterion, or G-criterion [32], and the minimum strain energy density criterion, or S-criterion [33]. Among them, the -criterion has been found to predict well the directions of crack initiation compared to the experimental results for polymethylmethacrylate [34, 35] and brittle clay [36]. Because of its simplicity, the -criterion seems to be the most popular criterion in mixed mode I-II fracture studies [37]. Therefore, the -criterion was used in this paper to determine the crack initiation angle for anisotropic plates.

For anisotropic materials, the general form of the elastic stress field near the crack tip in the local Cartesian coordinates of Figure 4 can be expressed in terms of the two stress intensity factors and a follows [30]:

Figure 4: Crack tip coordinate system and stress components.

Using coordinate transformation, the stress fields near the crack tip in the polar coordinates () of Figure 4 are

If the maximum -criterion is used, the angle of crack initiation, , must satisfy

A numerical procedure was applied to find the angle when is a maximum for known values of the material elastic constants, the anisotropic orientation angle , and the crack geometry.

In this paper, the process of crack propagation in an anisotropic homogeneous plate under mixed mode I-II loading is simulated by incremental crack extension with a piecewise linear discretization. For each incremental analysis, crack extension is conveniently modeled by a new boundary element. A computer program has been developed to automatically generate new data required for analyzing sequentially the changing boundary configuration. Based on the calculation of the SIFs and crack initiation angle for each increment, the procedure of crack propagation can be simulated. The steps in the crack propagation process are summarized as follows (Figure 5):(1)compute the SIFs using the proposed BEM;(2)determine the angle of crack initiation based on the maximum tensile stress criterion;(3)extend the crack by a linear element (of length selected by the user) along the direction determined in step ;(4)automatically generate the new BEM mesh;(5)repeat all the above steps until the crack is near the outer boundary.

Figure 5: Process of crack propagation by increasing the number of linear elements.

#### 3. Numerical Examples

In this section, a lot of numerical examples are presented to verify the formulation and to show the efficiency and versatility of the present SDBEM for problems related to fracture in anisotropic Bimaterial.

##### 3.1. Stress Intensity Factors (SIFs)
###### 3.1.1. Horizontal Crack in Material ()

A horizontal crack under a uniform pressure is shown in Figure 6. The crack has a length , and is located at a distance to the interface. The Poisson ratios for both materials () and () the same, that is, = = 0.3, while the ratio of the shear module ratio varies. A plane stress condition is assumed. In order to calculate the SIFs at crack tip () or (), 20 quadratic elements were used to discretize the crack surface. The results are given in Table 1 for various values of the shear module ratio. They are compared to the results given by Isida and Noguchi [38], using a body force integral equation method and those by Yuuki and Cho [39], using a multidomain BEM formulation. As can be observed from this table, the results compare quite well.

Table 1: Comparison of the SIFs (horizontal crack).
Figure 6: A horizontal crack under uniform pressure within material () of an infinite Bimaterial.
###### 3.1.2. Vertical Crack in Material ()

Consider a vertical crack in material () subjected to far-field horizontal stresses as shown in Figure 7. The crack has a length and is located at a distance to the interface. The Poisson ratios for both materials () and () are the same, that is, = 0.35, = 0.3, and the ratio of the shear module is the same, that is, = 23.077, while the ratio of the crack length located at a distance () varies. A plane stress condition is assumed. In order to calculate the SIFs at crack tip () or (), 20 quadratic elements were used to discretize the crack surface. The results are given in Table 2 for various values of the crack length located at a distance. They are compared to the results given by Isida and Noguchi [38] using a body force integral equation method and those by Cook and Erdogan [41] using a Wiener-Hopf technique and an asymptotic analysis. As can be observed from this table, the results compare quite well.

Table 2: Comparison of the SIFs (vertical crack).
Figure 7: A vertical crack in material (): an interface under far-field stresses.
###### 3.1.3. Vertical Crack Intersecting an Interface

Consider a vertical crack intersecting an interface and subjected to far-field horizontal stresses as shown in Figure 8. The horizontal far field stresses applied in materials () and () are, respectively, and (=). The Poisson ratio and are assumed to be equal to 0.3 and the shear modulus ratio is assumed to vary. The distance of the crack tip () and () to the interface are the same, that is, == , the half-length of the crack. Again, a plane stress condition is assumed and 20 quadratic elements were used to discretize the crack surface. The SIFs at the crack tips () and () are listed in Table 3 for several values of the shear modulus ratio, and are compare to those given by Isida and Noguchi [38]. Again, the results between the two numerical analyses compare quite well.

Table 3: Comparison of the SIFs (vertical crack intersecting an interface).
Figure 8: A vertical crack intersecting an interface under far-field stresses.
###### 3.1.4. An Interface Kinked Crack in Infinite Bimaterials

Consider an interface kinked () crack in infinite Bimaterials subjected to far-field tensile stresses as shown in Figure 9. The far-field tensile stresses applied in materials () and () are, respectively, and (=). The kink crack length is ; the main crack length is . Kink crack’s tip is (). The Poisson ratios and are assumed to be equal to 0.3. The Young module in material () is 1 GPa, and according to the shear module ratio with material () and material () ( = 0.25), we can get the Young module of material (). Then according to equation , we can get the relation between kink length and stress intensity factors. After numerical analysis, we compared with the results of Isida and Noguchi [38], using a body integral equation force method, as shown in Table 4, the results the two numerical analyses and found that they compare quite well.

Table 4: Comparison of the SIFs (interfacial kinked crack).
Figure 9: Interfacial kinked crack within infinite Bimaterials.
###### 3.1.5. Interfacial Horizontal Crack in Infinite Bimaterials

An interfacial crack along the x-axis of length is shown in Figure 10. The crack surface is under a uniform pressure and the materials can be either isotropic or anisotropic. Twenty quadratic elements were used to discretize the crack, and the characteristic length is assumed as .

Figure 10: Interfacial crack within infinite Bimaterials.

The SIFs at the crack tip of an interfacial crack were also calculated for the anisotropic Bimaterial case. The anisotropic elastic properties in material () were assumed to be those of glass/epoxy with = 48.26 GPa, = 17.24 GPa, = 6.89 GPa, and = 0.29. For material (), a graphite/epoxy with = 114.8 GPa, = 11.7 GPa, = 9.66 GPa, and = 0.21 was selected [42]. The material axis in material () and material () makes angles and , respectively, with respect to the horizontal direction (Figure 10). While the material axis in material () was assumed to be along the horizontal direction (i.e., = 0), the -axis in material () makes an angle with respect to the horizontal direction. The interfacial SIFs at crack tip () obtained by the present method are listed in Table 5 and compared to the exact solutions proposed by Wu [12]. A very good agreement is found between the numerical analysis and the exact solution.

Table 5: Comparison of the SIFs for infinite anisotropic problem (interfacial crack).
###### 3.1.6. Interfacial Horizontal Crack in Finite Bimaterials

The example is included as a comparison with the literature in order to demonstrate the accuracy of SDBEM approach for an interfacial crack in anisotropic bimaterial plate. The geometry is that of rectangular plate and is shown in Figure 11. For the comparison, crack length is taken as = 2, = , and = 0.4, and static tensile loading is applied on the upper and the lower boundary of the plate. Plane stress condition is assumed. The anisotropic elastic properties for materials () and () are given in Table 6. The normalized complex stress intensity factors at crack tip () or () are listed in Table 7 together with those from the work of Cho et al. [43], who used a multidomain BEM formulation and the results from Wünsche et al. [21] for a finite body. The outer boundary and interfacial crack surface were discretized with 40 continuous and 20 discontinuous quadratic elements, respectively. It is obvious from Table 7 that these are very close to those obtained by the other researchers.

Table 6: Elastic properties for materials () and ().
Table 7: Comparison of the normalized complex SIFs for finite anisotropic problem (interfacial crack).
Figure 11: An interfacial crack within finite rectangular plate of Bimaterial.
##### 3.2. Crack Initiation Angles
###### 3.2.1. Comparison of Numerical Predictions of Crack Initiation Angles with Experimental Results

In this section, we compared the numerical predictions of crack initiation angles with experimental results to verify the formulation. Al-Shayea [44] conducted uniaxial pressure test on limestone rock discs 98 mm and 84 mm in diameter and 22 mm in thickness with 30 mm notch. The crack orientation angle between the crack plane and the tensile stress was varied. Figure 12 shows the variation of the crack initiation angle with the crack angle determined numerically and experimentally. A good agreement is found between the experimental results of Al-Shayea [44] and our numerical predictions.

Figure 12: A variation of crack initiation angle with the crack angle .

#### 4. Comparison of Numerical Predictions of Crack Propagation Paths with Experimental Results

To demonstrate the proposed SDBEM procedure when predicting crack propagation in the anisotropic Bimaterials under mixed mode I–II loading, the propagation path in a CSTBD specimen is numerically predicted and compared with the actual laboratory observations. In these experiments, a crack initially inclined with respect to the applied stress is allowed to grow under concentrated diametrical loading (as shown in Figure 13). The Brazilian tests on CSTBD specimens with a diameter of 7.4 cm, a thickness of 1 cm, and a crack length of 2.2 cm are conducted to observe the actual propagation paths and are compared with the numerical predictions. The anisotropic elastic properties for rocks of types () and () are given in Table 8. The ratios of and are 1.635 and 4.301, respectively. Since the value of = 1.635, this rock of type () can be classified as a slightly anisotropic rock. Two specimens with the material inclination angle = and , defined as the SSD-1 and SSD-2, have crack angles = and , respectively. After Brazilian tests with cracked bimaterial specimens, the paths of crack propagation for SSD-1 and SSD-2 are shown in Figures 14(a) and 15(a), respectively. It can be observed that the crack propagates nearly perpendicular to the crack surface in the initial stage and then rapidly approaches toward the loading point. The proposed SDBEM procedure is also used to simulate crack propagation in the CSTBD specimens. The outer boundary and crack surface are discretized with 28 continuous and 20 discontinuous quadratic elements, respectively. Figures 14(b) and 15(b) are the comparisons examples; it can be concluded that the proposed SDBEM is capable of predicting the crack propagation in anisotropic bimaterial rocks.

Table 8: Elastic properties for the rocks of type () and ().
Figure 13: Geometry of a cracked straight through brazilian disc (CSTBD) specimen of anisotropic Bimaterial rock under diametrical loading.
Figure 14: Comparison between experimental observations and numerical predictions for specimen SSD-1: (a) photograph of specimen SSD-1 after failure ( = and = ) and (b) propagation of a crack at the center of a CSTBD specimen with = and = .
Figure 15: Comparison between experimental observations and numerical predictions for specimen SSD-2: (a) photograph of specimen SSD-2 after failure ( = and = ) and (b) propagation of an interfacial crack at the center of a CSTBD specimen with = and = .

#### 5. Conclusions

This paper shows that the mixed mode stress intensity factors of anisotropic Bimaterial rock under diametrical loading can be successfully determined by the SDBEM.

A new SDBEM procedure based on the maximum tensile stress failure criterion was developed to predict the crack initiation direction and the crack propagation path in anisotropic Bimaterial rock discs under mixed mode loading. A good agreement was found between crack initiation angles and propagation paths predicted with the SDBEM and experimental observations reported by previous researchers on anisotropic materials. Numerical simulations of crack initiation and propagation in CSTBD specimens of type ()/() were also found to compare well with the experimental results. Since the present method is simple and can be used for curved cracks, it will be straightforward to extend the current SDBEM formulation to analyze fracture propagation in 2D anisotropic Bimaterials, which is currently under investigation by the authors.

#### References

1. M. L. Williams, “The stresses around a faule or crack in dissimilar media,” Bulletin of Seismological Society of America, vol. 49, no. 2, pp. 199–204, 1959.
2. A. H. England, “A crack between dissimilar media,” Journal of Applied Mechanics, vol. 32, pp. 400–402, 1965.
3. J. R. Rice and G. Sih, “Plane problems of cracks in dissimilar media,” Journal of Applied Mechanics, Transactions ASME, vol. 32, pp. 418–423, 1965.
4. J. R. Rice, “Elastic fracture mechanics concepts for interfacial crack,” Journal of Applied Mechanics, Transactions ASME, vol. 55, no. 1, pp. 98–103, 1988.
5. D. L. Clements, “A crack between dissimilar anisotropic media,” International Journal of Engineering Science, vol. 9, no. 2, pp. 257–265, 1971.
6. J. R. Willis, “Fracture mechanics of interfacial cracks,” Journal of the Mechanics and Physics of Solids, vol. 19, no. 6, pp. 353–368, 1971.
7. T. C. T. Ting, “Explicit solution and invariance of the singularities at an interface crack in anisotropic composites,” International Journal of Solids and Structures, vol. 22, no. 9, pp. 965–983, 1986.
8. T. C. T. Ting, “Interface cracks in anisotropic bimaterials,” Journal of the Mechanics and Physics of Solids, vol. 38, no. 4, pp. 505–513, 1990.
9. J. L. Bassani and J. Qu, “Finite crack on bimaterial and bicrystal interfaces,” Journal of the Mechanics and Physics of Solids, vol. 37, no. 4, pp. 435–453, 1989.
10. V. K. Tewary, R. H. Wagoner, and J. P. Hirth, “Elastic Green's function for a composite solid with a planar interface,” Journal of Materials Research, vol. 4, no. 1, pp. 113–123, 1989.
11. Z. Suo, “Singularities, interfaces and cracks in dissimilar anisotropic media,” Proceedings of the Royal Society of London. Series A, vol. 427, no. 1873, pp. 331–358, 1990.
12. K. C. Wu, “Stress intensity factor and energy release rate for interfacial cracks between dissimilar anisotropic materials,” Journal of Applied Mechanics, Transactions ASME, vol. 57, pp. 882–886, 1990.
13. W. Kuang-Chong, “Explicit crack-tip fields of an extending interface crack in an anisotropic bimaterial,” International Journal of Solids and Structures, vol. 27, no. 4, pp. 455–466, 1991.
14. J. Qu and Q. Li, “Interfacial dislocation and its applications to interface cracks in anisotropic bimaterials,” Journal of Elasticity, vol. 26, no. 2, pp. 169–195, 1991.
15. H. Gao, M. Abbudi, and D. M. Barnett, “Interfacial crack-tip field in anisotropic elastic solids,” Journal of the Mechanics and Physics of Solids, vol. 40, no. 2, pp. 393–416, 1992.
16. M. S. Wu and H. Zhou, “Interactions of kinked interfacial cracks,” International Journal of Solids and Structures, vol. 36, no. 2, pp. 241–268, 1999.
17. E. Pan and B. Amadei, “Boundary element analysis of fracture mechanics in anisotropic bimaterials,” Engineering Analysis with Boundary Elements, vol. 23, no. 8, pp. 683–691, 1999.
18. S. J. Ping, L. X. Hui, and C. Y. Heng, “A complex variable boundary element method for solving interface crack problems,” International Journal of Fracture, vol. 96, no. 2, pp. 167–178, 1999.
19. W.-Y. Tian, K.-T. Chau, and Y.-H. Chen, “J-integral analysis of the interaction between an interface crack and parallel subinterface cracks in dissimilar anisotropic materials,” International Journal of Fracture, vol. 111, no. 4, pp. 305–325, 2001.
20. Z. C. Ou and Y. H. Chen, “Near-tip stress fields and intensity factors for an interface crack in metal/piezoelectric bimaterials,” International Journal of Engineering Science, vol. 42, no. 13-14, pp. 1407–1438, 2004.
21. M. Wünsche, Ch. Zhang, J. Sladek, V. Sladek, S. Hirose, and M. Kuna, “Transient dynamic analysis of interface cracks in layered anisotropic solids under impact loading,” International Journal of Fracture, vol. 157, no. 1-2, pp. 131–147, 2009.
22. S. G. Lekhnitskii, Theory of Elasticity of An Anisotropic Elastic Body, Holden-day, San Francisco, Calif, USA, 1963.
23. T. C. T. Ting, Anisotropic Elasticity:Theory and Application, Oxford University Press, New York, NY, USA, 1996.
24. E. Pan, “A general boundary element analysis of 2-D linear elastic fracture mechanics,” International Journal of Fracture, vol. 88, no. 1, pp. 41–59, 1997.
25. P. Sollero and M. H. Aliabadi, “Fracture mechanics analysis of anisotropic plates by the boundary element method,” International Journal of Fracture, vol. 64, no. 4, pp. 269–284, 1993.
26. C. S. Chen, Characterization of deformability, strength, and fracturing of anisotropic rocks using Brazilian tests, Ph.D. thesis, Department of Civil Engineering, University of Colorado, 1996.
27. J. T. Chen and H.-K. Hong, “Review of dual boundary element methods with emphasis on hypersingular integrals and divergent series,” Applied Mechanics Reviews, vol. 52, no. 1, pp. 17–32, 1999.
28. S. L. Crouch and A. M. Starfield, Boundary Element Method in Solid Mechanics, George Allen and Unwin, London, UK, 1983.
29. T. A. Cruse , Boundary element analysis in computational fracture mechanics, Kluwer Academic, Dordrecht, 1988.
30. G. C. Sih, P. C. Paris, and G. R. Irwin, “On cracks in rectilinearly anisotropic bodies,” International Journal of Fracture Mechanics, vol. 1, no. 3, pp. 189–203, 1965.
31. F. Erdogan and G. C. Sih, “On the crack extension in plates under loading and transverse shear,” Journal of Basic Engineering, vol. 58, pp. 519–527, 1963.
32. K. Palaniswamy and W. G. Knauss, “Propagation of a crack under general, in-plane tension,” International Journal of Fracture Mechanics, vol. 8, no. 1, pp. 114–117, 1972.
33. G. C. Sih, “Strain-energy-density factor applied to mixed mode crack problems,” International Journal of Fracture, vol. 10, no. 3, pp. 305–321, 1974.
34. C. W. Woo and L. H. Ling, “On angle crack initiation under biaxial loading,” Journal of Strain Analysis for Engineering Design, vol. 19, no. 1, pp. 51–59, 1984.
35. H. A. Richard, “Examination of brittle fractured criteria for overloapping mode I and mode II loading applied to cracks,” in Application of Fracture Mechanics to Materials and Structures, G. C. Sih et al., Ed., pp. 309–316, Martinus Nijhoff, Leiden, The Netherlands, 1984.
36. L. E. Vallejo, “The brittle and ductile behavior of a material containing a crack under mixed-mode loading,” in Proceedings of the 28th US Symposium on Rock Mechanics (USRMS '87), pp. 383–390, University of Arizona, Tucson, Ariz, USA, 1987.
37. B. N. Whittaker, R. N. Singh, and G. Sun, Rock Fracture Mechanics: Principles, Design and Applications, Elsevier, New York, NY, USA, 1992.
38. M. Isida and H. Noguchi, “Plane problems of arbitrarily oriented cracks in bonded dissimilar materials,” Transactions of the Japan Society of Mechanical Engineers, vol. 49, pp. 36–45, 1983.
39. Y. Ryoji and C. Sang-Bong, “Efficient boundary element analysis of stress intensity factors for interface cracks in dissimilar materials,” Engineering Fracture Mechanics, vol. 34, no. 1, pp. 179–188, 1989.
40. M. Isida and H. Noguchi, “Plane elastostatic problems of bonded dissimilar materials with an interface crack and arbitrarily located cracks,” Transactions of the Japan Society of Mechanical Engineers, vol. 49, pp. 137–146, 1983.
41. T. S. Cook and F. Erdogan, “Stresses in bonded materials with a crack perpendicular to the interface,” International Journal of Engineering Science, vol. 10, no. 8, pp. 677–697, 1972.
42. J. F. Dwyer and E. Pan, “Edge function analysis of stress intensity factors in cracked anisotropic plates,” International Journal of Fracture, vol. 72, no. 4, pp. 327–342, 1995.
43. S. B. Cho, K. R. Lee, Y. S. Choy, and R. Yuuki, “Determination of stress intensity factors and boundary element analysis for interface cracks in dissimilar anisotropic materials,” Engineering Fracture Mechanics, vol. 43, no. 4, pp. 603–614, 1992.
44. N. A. Al-Shayea, “Crack propagation trajectories for rocks under mixed mode I-II fracture,” Engineering Geology, vol. 81, no. 1, pp. 84–97, 2005.