Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2009 / Article

Research Article | Open Access

Volume 2009 |Article ID 132980 | 10 pages | https://doi.org/10.1155/2009/132980

Dual Boundary Element Method Applied to Antiplane Crack Problems

Academic Editor: Angelo Luongo
Received23 Feb 2009
Revised07 Jul 2009
Accepted03 Aug 2009
Published05 Nov 2009

Abstract

This paper is concerned with an efficient dual boundary element method for 2d crack problems under antiplane shear loading. The dual equations are the displacement and the traction boundary integral equations. When the displacement equation is applied on the outer boundary and the traction equation on one of the crack surfaces, general crack problems with anti-plane shear loading can be solved with a single region formulation. The outer boundary is discretised with continuous quadratic elements; however, only one of the crack surfaces needs to be discretised with discontinuous quadratic elements. Highly accurate results are obtained, when the stress intensity factor is evaluated with the discontinuous quarter point element method. Numerical examples are provided to demonstrate the accuracy and efficiency of the present formulation.

1. Introduction

The problem of a cracked body subjected to an antiplane shear loading had been studied extensively. Sih [1] provided analytical solutions for mode III cracks in infinite regions by using Westergaard stress functions and Muskhelishvili's method. Chiang [2] presented analytical solutions for slightly curved cracks in antiplane strain in infinite regions using perturbation procedures similar to those carried out for in-plane loading cases by Cotterell and Rice [3]. Zhang [4, 5] and Ma and Zhang [6] gave analytical solutions for a mode III stress intensity factor considering a finite region with an eccentric straight crack. Ma [7] provided analytical solutions for mode III straight cracks in finite regions using Fourier transforms and Fourier series. Smith [8] studied the elastic stress distribution in the immediate vicinity of a blunt notch. However, their solutions were concerned with specified geometries or boundary conditions. To deal with the complexities of general geometries and boundary conditions, an accurate and efficient numerical method is essential [9โ€“12].

Several numerical solutions had been devised for antiplane crack problems. Wallentin et al. [13] investigated the railway wheel crack problem numerically, based on Betti's reciprocity theorem. Guagliano and Vergani [14] described the experimental and numerical analysis of internal cracks in wheels under Hertzian loads. Paulino et al. [15] provided numerical solutions for a curved crack subjected to an antiplane shear loading in finite regions by using the boundary integral equation method. Ting et al. [16] provided numerical solutions for mode III crack problems by using the boundary element alternating method. Liu and Altiero [17] provided numerical solutions for mode III crack problems using the boundary integral equation with linear approximation on displacements and stresses. Barlow and Chandra [18] discussed the computational fatigue crack growth rate by using the crack opening displacement approach to calculate the stress intensity factors. Mews and Kuhn [19] provided numerical solutions for the traction free central crack problem by using Green's function, instead of the usual fundamental solution. Mir-Mohamad-Sadegh and Altiero [20] used the indirect boundary integral equation method to solve traction problems, using displacement-based formulations. Sun et al. [21] derived a new boundary integral equation to analyse cracked anisotropic bodies under antiplane shear. Also, for the further study, the crack front plastic deformation in a ductile material was introduced to apply the effective Dugdale strip yield model [22โ€“24]. In general, the boundary element method (BEM) is a well-established numerical technique for the analysis of linear fracture mechanics problems. However, the solution of general crack problems cannot be achieved with the direct application of the BEM, because the coincidence of the crack surfaces gives rise to a singular system of algebraic equations.

To overcome this shortcoming, we provide an efficient numerical procedure, based on the dual boundary element method (DBEM), for antiplane shear loading problems. The dual boundary element method seems to have certain apparent advantages for in-plane loading problems with a single region formulation. This method incorporates two independent boundary integral equations, the displacement and traction equations. Portela et al. [25] considered the effective numerical implementation of the two-dimensional DBEM for solving general in-plane fracture mechanics problems. W. H. Chen and T. C. Chen [26] proposed a different DBEM formulation for in-plane crack problems. Chen and Chen suggested the use of the displacement integral equation applied only on the outer boundary and the traction integral equation on one of the crack surfaces. In Chen and Chen's formulation, relative displacement of crack surfaces was used instead of the displacement. This reduces the degrees of freedom and hence the computational effort. This study uses an integral equation formulation that combines with the crack modelling strategy of quadratic boundary elements for antiplane crack problems. The stress intensity factor is calculated based on the near tip displacement method. More accurate results are obtained by placing discontinuous quarter point elements at crack tips [27], which correctly model the behaviour of the crack tip displacement. This is a similar technique to that used for continuous quarter point elements [28]. Numerical examples are provided to demonstrate the accuracy and efficiency of the present formulation.

2. The Dual Boundary Integral Equation for Antiplane Problems

Consider a finite domain subjected to an arbitrary antiplane shear loading, where the only nonzero displacement component ๐‘ข๐‘ง in the ๐‘ง direction may be specified as follows [7]: โˆ‡2๐‘ข๐‘ง=0.(2.1) The Laplace equation (2.1) can be transformed into a boundary integral equation, as is typical with the BEM. The boundary integral formulation of the displacement component, ๐‘ข๐‘ง, at an internal point ๐ˆ, is given by [29] ๐‘ข๐‘ง(๎€œ๐ˆ)+ฮ“๐ป(๐ˆ,๐ฑ)๐‘ข๐‘ง(๎€œ๐ฑ)๐‘‘ฮ“(๐ฑ)=ฮ“๐บ(๐ˆ,๐ฑ)๐‘ก๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ),(2.2) where ๐‘ก๐‘ง represents the traction component at a boundary point ๐ฑ. ๐ป(๐ˆ,๐ฑ) and ๐บ(๐ˆ,๐ฑ) represent the fundamental traction and displacement solutions, respectively, which are given as 1๐ป(๐ˆ,๐ฑ)=โˆ’2๐œ‹๐‘Ÿ๐œ•๐‘Ÿ1๐œ•๐ง,๐บ(๐ˆ,๐ฑ)=๎‚€12๐œ‹๐œ‡ln๐‘Ÿ๎‚,(2.3) where ๐œ‡ is the shear modulus, ๐‘Ÿ is the distance between ๐ˆ and ๐ฑ, and ๐ง denotes the outward normal unit vector at the point ๐ฑ on the boundary ฮ“. If we consider a finite body with ๐ฟ cracks, (2.2) can be written as ๐‘ข๐‘ง๎€œ(๐ˆ)+ฮ“๐‘†๐ป(๐ˆ,๐ฑ)๐‘ข๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ)+๐ฟ๎“๐‘™=1๎€œฮ“+๐‘™๐ป๎€ท๐ˆ,๐ฑ+๎€ธ๐‘ข๐‘ง๎€ท๐ฑ+๎€ธ๐‘‘ฮ“(๐ฑ)+๐ฟ๎“๐‘™=1๎€œฮ“โˆ’๐‘™๐ป(๐ˆ,๐ฑโˆ’)๐‘ข๐‘ง(๐ฑโˆ’=๎€œ)๐‘‘ฮ“(๐ฑ)ฮ“๐‘†๐บ(๐ˆ,๐ฑ)๐‘ก๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ)+๐ฟ๎“๐‘™=1๎€œฮ“+๐‘™๐บ๎€ท๐ˆ,๐ฑ+๎€ธ๐‘ก๐‘ง๎€ท๐ฑ+๎€ธ๐‘‘ฮ“(๐ฑ)+๐ฟ๎“๐‘™=1๎€œฮ“โˆ’๐‘™๐บ(๐ˆ,๐ฑโˆ’)๐‘ก๐‘ง(๐ฑโˆ’)๐‘‘ฮ“(๐ฑ),(2.4) where ๐ฑ+ and ๐ฑโˆ’ are the field points located on upper and lower crack surfaces, respectively. Note that ฮ“๐‘† denotes the outer boundary of the body, ฮ“+๐‘™ the ๐‘™th upper crack boundary, ฮ“โˆ’๐‘™ the ๐‘™th lower crack boundary, and ฮ“=ฮ“๐‘†+โˆ‘๐ฟ๐‘™=1(ฮ“+๐‘™+ฮ“โˆ’๐‘™). Using the fact that ๐ป(๐ˆ,๐ฑ+)โˆฃฮ“+=๐ป(๐ˆ,๐ฑโˆ’)โˆฃโˆ’ฮ“โˆ’ and ๐บ(๐ˆ,๐ฑ+)โˆฃฮ“+=๐บ(๐ˆ,๐ฑโˆ’)โˆฃโˆ’ฮ“โˆ’, (2.4) can be simplified to ๐‘ข๐‘ง๎€œ(๐ˆ)+ฮ“๐‘†๐ป(๐ˆ,๐ฑ)๐‘ข๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ)+๐ฟ๎“๐‘™=1๎€œฮ“+๐‘™๐ป๎€ท๐ˆ,๐ฑ+๎€ธฮ”๐‘ข๐‘ง=๎€œ(๐ฑ)๐‘‘ฮ“(๐ฑ)ฮ“๐‘†๐บ(๐ˆ,๐ฑ)๐‘ก๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ)+๐ฟ๎“๐‘™=1๎€œฮ“+๐‘™๐บ๎€ท๐ˆ,๐ฑ+๎€ธฮ”๐‘ก๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ),(2.5) where ฮ”๐‘ข๐‘ง=๐‘ข๐‘ง(๐ฑ+)โˆ’๐‘ข๐‘ง(๐ฑโˆ’) and ฮ”๐‘ก๐‘ง=๐‘ก๐‘ง(๐ฑ+)โˆ’๐‘ก๐‘ง(๐ฑโˆ’), however ฮ”๐‘ก๐‘ง is always zero on the crack faces. As the internal point approaches the outer boundary, that is, as ๐ˆโ†’๐, the displacement equation becomes ๐‘(๐)๐‘ข๐‘ง๎€ง(๐)+ฮ“๐‘†๐ป(๐,๐ฑ)๐‘ข๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ)+๐ฟ๎“๐‘™=1๎€œฮ“+๐‘™๐ป๎€ท๐,๐ฑ+๎€ธฮ”๐‘ข๐‘ง๎€œ(๐ฑ)๐‘‘ฮ“(๐ฑ)=ฮ“๐‘†๐บ(๐,๐ฑ)๐‘ก๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ),(2.6) where โจ represents the Cauchy principle value integral and ๐‘(๐)=1/2, given a smooth boundary at the point ๐.

The stress components ๐œŽ๐‘–๐‘ง are obtained from differentiation of equation (2.5), followed by the application of Hooke's law. At an internal point ๐ˆ, these components are given by ๐œŽ๐‘–๐‘ง๎€œ(๐ˆ)+ฮ“๐‘†๐‘†๐‘–(๐ˆ,๐ฑ)๐‘ข๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ)+๐ฟ๎“๐‘™=1๎€œฮ“+๐‘™๐‘†๐‘–(๐ˆ,๐ฑ)ฮ”๐‘ข๐‘ง๎€œ(๐ฑ)๐‘‘ฮ“(๐ฑ)=ฮ“๐‘†๐ท๐‘–(๐ˆ,๐ฑ)๐‘ก๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ),(2.7) where ๐‘†๐‘–(๐ˆ,๐ฑ) and ๐ท๐‘–(๐ˆ,๐ฑ) contain derivatives of ๐ป(๐ˆ,๐ฑ) and ๐บ(๐ˆ,๐ฑ) in the ๐‘– direction, respectively, which are given as ๐‘†๐‘–๐œ‡(๐ˆ,๐ฑ)=2๐œ‹๐‘Ÿ2๎‚ธ๐œ•๐‘Ÿ๐œ•๐‘ฅ๐‘–๐œ•๐‘Ÿโˆ’๎‚ต๐›ฟ๐œ•๐ง๐‘–๐‘—โˆ’๐œ•๐‘Ÿ๐œ•๐‘ฅ๐‘—๐œ•๐‘Ÿ๐œ•๐‘ฅ๐‘–๎‚ถ๐‘›๐‘—๎‚น,๐ท๐‘–1(๐ˆ,๐ฑ)=โˆ’2๐œ‹๐‘Ÿ๐œ•๐‘Ÿ๐œ•๐‘ฅ๐‘–,(2.8) where ๐‘›๐‘– denotes the ๐‘–th component of the outward normal to the boundary at point ๐ฑ, and ๐›ฟ๐‘–๐‘— is the Kronecker delta. Again, by moving the source point ๐ˆ to the upper crack boundary ๐, and using ๐‘ก๐‘ง=๐œŽ๐‘–๐‘ง๐‘›๐‘–, we obtain the traction integral equation 12๐‘ก๐‘ง๎€œ(๐)+ฮ“๐‘†๐‘›๐‘–(๐)๐‘†๐‘–(๐,๐ฑ)๐‘ข๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ)+๐ฟ๎“ฮ“๐‘™=1+๐‘™๐‘›๐‘–(๐)๐‘†๐‘–(๐,๐ฑ)ฮ”๐‘ข๐‘ง=๎€œ(๐ฑ)๐‘‘ฮ“(๐ฑ)ฮ“๐‘†๐‘›๐‘–(๐)๐ท๐‘–(๐,๐ฑ)๐‘ก๐‘ง(๐ฑ)๐‘‘ฮ“(๐ฑ),(2.9) where represents the Hadamard principal value integral. Both Cauchy and Hadamard principal-value integrals in (2.6) and (2.9) are finite parts of improper integrals. To solve the finite part integrals, we can follow the method mentioned in Portela et al. [25].

The displacement integral equation (2.6) and the traction integral equation (2.9) are the two main equations to solve for the displacement of the outer boundary and the relative displacement of the crack faces. Equation (2.6) is applied for collocation on the outer boundary where continuous quadratic elements are used, and (2.9) is applied on the upper crack faces which are modelled by discontinuous quadratic elements. By taking all the discretised nodes on the outer boundary ฮ“๐‘† and upper crack surfaces โˆ‘๐ฟ๐‘™=1ฮ“+๐‘™ at the source point ๐, the system of (2.6) and (2.9) for the multiple cracks problem can be written in a matrix form as ๎‚ธ๐‡1๐‡20๐’1๐’2๐ˆ๎‚นโŽกโŽขโŽขโŽฃ๐ฎ๐‘ง,๐‘†ฮ”๐ฎ๐‘ง,๐‘๐ญ๐‘ง,๐‘+โŽคโŽฅโŽฅโŽฆ=๎‚ธ๐†1๐ƒ1๎‚น๎€บ๐ญ๐‘ง,๐‘†๎€ป,(2.10) where ๐‡1, ๐‡2, ๐†1 and ๐’1, ๐’2, ๐ƒ1 are the corresponding assembled matrices from (2.6) and (2.9), respectively. The ๐ฎ๐‘ง,๐‘† and ๐ญ๐‘ง,๐‘† are the displacement and traction vectors on the outer boundary ฮ“๐‘†, respectively. ฮ”๐ฎ๐‘ง,๐‘ and ๐ญ๐‘ง,๐‘+ are the relative displacement vector and the traction vector on the upper crack faces.

3. Calculation of the Mode III Stress Intensity Factor

Near tip displacement extrapolation is used to evaluate the numerical values of the stress intensity factor. The relative displacements of the crack surfaces are calculated using the DBEM and are used in the near crack tip stress field equations to obtain the stress intensity factor. Due to the singular behaviour of the stress around the crack tip, it is reasonable to expect a better approximation by replacing the normal discontinuous quadratic element with a transition element possessing the same order of singularity at the crack tip. The discontinuous quarter point element method is used in the present formulation [27, 30]. The mode III stress intensity factor is evaluated as ๐พIII=๐œ‡4๎‚™2๐œ‹๐‘Ÿฮ”๐‘ข๐‘ง(๐‘Ÿ),(3.1) where ๐‘Ÿ is the distance from the crack tip to the nearest node on the upper crack face, and ฮ”๐‘ข๐‘ง(๐‘Ÿ) denotes the relative displacement in the antiplane direction.

4. Numerical Examples

In order to demonstrate the accuracy and efficiency of the technique previously described, and to illustrate possible applications, we now consider several examples. In all the numerical tests, the outer boundary is modelled by 24 continuous quadratic elements, and each crack discretization is carried out with three different meshes of 6, 8, and 10 discontinuous quadratic elements, respectively. The best accuracy is achieved with 6 elements, in which the crack discretization is graded, towards the tip, with ratios 0.25, 0.15, and 0.1. The plate is subjected to a uniform antiplane shear loading ๐œ, and the stress intensity factor is normalised with respect to ๐พ0โˆš=๐œ๐œ‹๐‘Ž,(4.1) where ๐‘Ž defines the half length of the crack. All computations are carried out under the condition of plane strain.

4.1. A Rectangular Plate Containing a Central Slant Crack

Firstly, consider a rectangular plate containing a central slant crack as shown in Figure 1. The crack has length 2๐‘Ž and makes an angle ๐œƒ with the horizontal direction. For a horizontal crack (๐œƒ=0โˆ˜), the normalised mode III stress intensity factor is calculated for various ratios of ๐‘Ž/โ„Ž and ๐‘Ž/๐‘ค and compared to those given in [17, 31] (see Table 1). The largest difference between these does not exceed 1.65 per cent. Further, the normalised mode III stress intensity factor is calculated for โ„Ž/๐‘ค=2, while the crack slanted an angle ๐œƒ with the various ratios of ๐‘Ž/๐‘ค. Three cases are considered, where ๐œƒ=30โˆ˜, 45ยฐ, and 60ยฐ, respectively. The results obtained are presented in Figure 2. As it can be seen, when the ratio of ๐‘Ž/๐‘ค increases, the stress intensity factor increases due to edge effect.


๐‘Ž โˆถ โ„Ž 1 โˆถ 0 . 2 5 1 โˆถ 0 . 5 1 โˆถ 1 1 โˆถ 2 1 โˆถ 4

๐‘Ž โˆถ ๐‘ค Present 1.909 1.724 1.689 1.688 1.661
1 โˆถ 1 . 2 Reference [17] 1.897 1.723 1.689 1.686 1.686
Reference [31] 1.900 1.725 1.691 1.689 1.689

๐‘Ž โˆถ ๐‘ค Present 1.796 1.467 1.371 1.361 1.361
1 โˆถ 1 . 4 Reference [17] 1.780 1.460 1.369 1.359 1.358
Reference [31] 1.782 1.463 1.370 1.361 1.360

๐‘Ž โˆถ ๐‘ค Present 1.784 1.405 1.257 1.236 1.236
1 โˆถ 1 . 6 Reference [17] 1.771 1.399 1.254 1.233 1.233
Reference [31] 1.773 1.401 1.256 1.235 1.235

๐‘Ž โˆถ ๐‘ค Present 1.792 1.384 1.179 1.131 1.129
1 โˆถ 2 . 0 Reference [17] 1.770 1.377 1.176 1.127 1.126
Reference [31] 1.772 1.379 1.178 1.130 1.128

For the case where ๐‘Ž/๐‘ค=1/50, which could be considered as the case of infinite geometry since ๐‘Žโ‰ช๐‘ค, we compare the results with the analytical results for the latter as given in [32]. The results are plotted in Figure 3. Excellent agreement is observed; the maximum error is around 0.02 per cent.

4.2. A Rectangular Plate Containing Two Identical Collinear Cracks

As shown in Figure 4, the second example is a rectangular plate containing two identical collinear cracks. 2๐‘Ž is the length of the inclined crack and 2๐‘‘ is the distance between the centre of the cracks. The geometric parameters are โ„Ž/๐‘ค=2 and ๐‘Ž/๐‘ค=1/50. Figure 4 displays the variations of normalised mode III stress intensity factors at tip ๐ด and tip ๐ต versus different ratios of ๐‘Ž/๐‘‘. Due to the interaction between the two cracks, the computed normalised mode III stress intensity factor at tip ๐ด is always larger than that at tip ๐ต. Hence, as the crack centre distance ๐‘‘ decreases, the difference of stress intensity factor increases. There is excellent correlation between the computed results using the present method and those from analytical solutions; the difference between these results does not exceed 0.03 per cent at tip ๐ด or 0.09 per cent at tip ๐ต.

4.3. An Infinite Plate Containing Two Parallel Cracks

The third example is an infinite plate (โ„Ž/๐‘ค=2, ๐‘Ž/๐‘ค=1/50) containing two parallel cracks, as shown in Figure 5. 2๐‘Ž is the length of the two identical cracks and 2๐‘‘ is the distance between the cracks. The computed results are compared with the published results in [32]. The results of normalised mode III stress intensity factor for different ๐‘  are plotted in Figure 5, where ๐‘ =๐‘Ž/(๐‘Ž+๐‘‘). The effect of the interaction of cracks on the mode III stress intensity factor is observed. The largest difference between the present and the published results does not exceed 0.65 per cent.

5. Conclusions

An efficient and accurate dual boundary element technique has been successfully developed for the analysis of two dimensional cracks subjected to an antiplane shear loading. The dual boundary equations are the usual displacement boundary integral equation and the traction boundary integral equation. When the displacement equation is applied on the outer boundary and the traction equation is applied on one of the crack surfaces, a general crack problem can be solved in a single region formulation. The discontinuous quarter point elements are used for evaluating the mode III stress intensity factor, which correctly describes the ๐‘Ÿ1/2 behaviour of the near tip displacements. This, therefore, allows accurate results for mode III stress intensity factors to be calculated.

References

  1. G. C. Sih, โ€œStress distribution near internal crack tips for longitudinal shear problems,โ€ Journal of Applied Mechanics, vol. 32, pp. 51โ€“58, 1965. View at: Google Scholar | Zentralblatt MATH | MathSciNet
  2. C. R. Chiang, โ€œSlightly curved cracks in antiplane strain,โ€ International Journal of Fracture, vol. 32, no. 4, pp. R63โ€“R66, 1986. View at: Publisher Site | Google Scholar
  3. B. Cotterell and J. R. Rice, โ€œSlightly curved or kinked cracks,โ€ International Journal of Fracture, vol. 16, pp. 155โ€“169, 1980. View at: Publisher Site | Google Scholar
  4. X. S. Zhang, โ€œThe general solution of a central crack off the center line of a rectangular sheet for mode III,โ€ Engineering Fracture Mechanics, vol. 28, no. 2, pp. 147โ€“155, 1987. View at: Publisher Site | Google Scholar
  5. X. S. Zhang, โ€œA tearing mode crack located anywhere in a finite rectangular sheet,โ€ Engineering Fracture Mechanics, vol. 33, no. 4, pp. 509โ€“516, 1989. View at: Publisher Site | Google Scholar
  6. S. W. Ma and L. X. Zhang, โ€œA new solution of an eccentric crack off the center line of a rectangular sheet for mode-III,โ€ Engineering Fracture Mechanics, vol. 40, no. 1, pp. 1โ€“7, 1991. View at: Publisher Site | Google Scholar
  7. S. W. Ma, โ€œA central crack of mode III in a rectangular sheet with fixed edges,โ€ International Journal of Fracture, vol. 39, no. 4, pp. 323โ€“329, 1989. View at: Publisher Site | Google Scholar
  8. E. Smith, โ€œA comparison of Mode I and Mode III results for the elastic stress distribution in the immediate vicinity of a blunt notch,โ€ International Journal of Engineering Science, vol. 42, no. 5-6, pp. 473โ€“481, 2004. View at: Publisher Site | Google Scholar
  9. P. H. Geubelle and M. S. Breitenfeld, โ€œNumerical analysis of dynamic debonding under anti-plane shear loading,โ€ International Journal of Fracture, vol. 85, no. 3, pp. 265โ€“282, 1997. View at: Publisher Site | Google Scholar
  10. E. Heymsfield, โ€œInfinite domain correction for anti-plane shear waves in a two-dimensional boundary element analysis,โ€ International Journal for Numerical Methods in Engineering, vol. 40, no. 5, pp. 953โ€“964, 1997. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  11. S. M. Mkhitaryan, N. Melkoumian, and B. B. Lin, โ€œStress-strain state of a cracked elastic wedge under anti-plane deformation with mixed boundary conditions on its faces,โ€ International Journal of Fracture, vol. 108, no. 4, pp. 291โ€“315, 2001. View at: Publisher Site | Google Scholar
  12. Z.-G. Zhou and L. Ma, โ€œTwo collinear Griffith cracks subjected to anti-plane shear in infinitely long strip,โ€ Mechanics Research Communications, vol. 26, no. 4, pp. 437โ€“444, 1999. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  13. M. Wallentin, H. L. Bjarnehed, and R. Lundén, โ€œCracks around railway wheel flats exposed to rolling contact loads and residual stresses,โ€ Wear, vol. 258, no. 7-8, pp. 1319โ€“1329, 2005. View at: Publisher Site | Google Scholar
  14. M. Guagliano and L. Vergani, โ€œExperimental and numerical analysis of sub-surface cracks in railway wheels,โ€ Engineering Fracture Mechanics, vol. 72, no. 2, pp. 255โ€“269, 2005. View at: Publisher Site | Google Scholar
  15. G. H. Paulino, M. T. A. Saif, and S. Mukherjee, โ€œA finite elastic body with a curved crack loaded in anti-plane shear,โ€ International Journal of Solids and Structures, vol. 30, no. 8, pp. 1015โ€“1037, 1993. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  16. K. Ting, K.-K. Chang, and M.-F. Yang, โ€œAnalysis of mode III crack problems by boundary element alternating method,โ€ in Fatigue, Flaw Evaluation and Leak-Before-Break Assessments, vol. 280, pp. 21โ€“26, American Society of Mechanical Engineers, 1994. View at: Google Scholar
  17. N. Liu and N. J. Altiero, โ€œAn integral equation bethod applied to mode III crack problems,โ€ Engineering Fracture Mechanics, vol. 41, pp. 578โ€“596, 1992. View at: Publisher Site | Google Scholar
  18. K. W. Barlow and R. Chandra, โ€œFatigue crack propagation simulation in an aircraft engine fan blade attachment,โ€ International Journal of Fatigue, vol. 27, no. 10–12, pp. 1661โ€“1668, 2005. View at: Publisher Site | Google Scholar
  19. H. Mews and G. Kuhn, โ€œAn effective numerical stress intensity factor calculation with no crack discretization,โ€ International Journal of Fracture, vol. 38, no. 1, pp. 61โ€“76, 1988. View at: Publisher Site | Google Scholar
  20. A. Mir-Mohamad-Sadegh and N. J. Altiero, โ€œSolution of the problem of a crack in a finite plane region using an indirect boundary-integral method,โ€ Engineering Fracture Mechanics, vol. 11, no. 4, pp. 831โ€“837, 1979. View at: Publisher Site | Google Scholar
  21. Y.-Z. Sun, S.-S. Yang, and Y.-B. Wang, โ€œA new formulation of boundary element method for cracked anisotropic bodies under anti-plane shear,โ€ Computer Methods in Applied Mechanics and Engineering, vol. 192, no. 22-23, pp. 2633โ€“2648, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  22. J. Vrbik, B. M. Singh, J. Rokne, and R. S. Dhaliwal, โ€œPlastica deformation at the tip of an edge crack,โ€ Zeitschrift für Angewandte Mathematick und Mechanik, vol. 81, pp. 642โ€“647, 2001. View at: Publisher Site | Google Scholar
  23. H. T. Danyluk, B. M. Singh, and J. Vrbik, โ€œPlastic zones in an orthotropic plate of finite width containing a Griffith crack,โ€ International Journal of Fracture, vol. 75, no. 4, pp. 307โ€“322, 1996. View at: Publisher Site | Google Scholar
  24. X. Jin, S. Chaiyat, L. M. Keer, and K. Kiattikomol, โ€œRefined Dugdale plastic zones of an external circular crack,โ€ Journal of the Mechanics and Physics of Solids, vol. 56, no. 4, pp. 1127โ€“1146, 2008. View at: Publisher Site | Google Scholar
  25. A. Portela, M. H. Aliabadi, and D. P. Rooke, โ€œThe dual boundary element method: effective implementation for crack problems,โ€ International Journal for Numerical Methods in Engineering, vol. 33, no. 6, pp. 1269โ€“1287, 1992. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  26. W. H. Chen and T. C. Chen, โ€œAn efficient dual boundary element technique for a two-dimensional fracture problem with multiple cracks,โ€ International Journal for Numerical Methods in Engineering, vol. 38, pp. 1739โ€“1756, 1995. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  27. A. Saez, R. Gallego, and J. Dominguez, โ€œHypersingular quarter-point boundary elements for crack problems,โ€ International Journal for Numerical Methods in Engineering, vol. 38, pp. 1681โ€“1701, 1995. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  28. J. Martinez and J. Dominguez, โ€œOn the use of quarter-point boundary elements for stress intensity factor computations,โ€ International Journal for Numerical Methods in Engineering, vol. 20, no. 10, pp. 1941โ€“1950, 1984. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  29. C. A. Brebbia and J. Dominguez, Boundary Elements: An Introductory Course, Computational Mechanics, Southampton, UK, 1989. View at: MathSciNet
  30. P. Fedelinski, M. H. Aliabadi, and D. P. Rooke, โ€œA single-region time domain BEM for dynamic crack problems,โ€ International Journal of Solids and Structures, vol. 32, no. 24, pp. 3555โ€“3571, 1995. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  31. S. W. Ma, โ€œA central crack in a rectangular sheet where its boundary is subjected to an arbitrary anti-plane load,โ€ Engineering Fracture Mechanics, vol. 30, no. 4, pp. 435โ€“443, 1988. View at: Publisher Site | Google Scholar
  32. H. Tada, P. Paris, and G. Irwin, The Stress Analysis of Cracks Handbook, Paris Productions Incorporated, St. Louis, Mo, USA, 1985.

Copyright © 2009 Wei-Liang Wu. 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.


More related articles

722ย Views | 524ย Downloads | 7ย Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.