Abstract: With the aid of integral transforms, a method is presented for solving the problem of scattering of plane harmonic compression and shear waves by a line of discontinuity or crack of finite width embedded in an elastic medium of infinite extent. When the incoming waves are applied in an arbitrary direction, the scattered-wave field may be determined by separating the crack-surface boundary conditions into functions even and odd with respect to the variable along the line crack. The problem is reduced to the evaluation of a system of coupled Fredholm integral equations with special emphasis placed 011 finding the near-field solution which consists of a knowledge of the detailed structure of the displacements and stresses in a small region around the crack vertex. Dynamic stress-intensity factors, the critical values of which govern the condition of crack propagation, are defined and found to be dependent on the incident wave length and Poisson's ratio of the medium. At certain wave lengths, they are larger than those encountered under static loading. Such information is of particular importance in perdicting the fracture strength of structures subjected to oscillating loads. Introduction. Although the scattering of waves by obstacles of different shapes has been the subject of many past investigations in various branches of physics [1]—[3], to the authors' knowledge none of these investigations analyzed, in detail, the singular behavior of the stresses near a scatterer in the form of a line of discontinuity or finite crack. The main reason for this omission is the lack of an effective mathematical method for obtaining the near-field solution, which is of considerable theoretical interest and has innumerable applications in the field of fracture mechanics as well as in electromagnetic and acoustic theory. A popular approach to the diffraction of waves from obstacles has been that of separation of variables, where the formal solution of the wave equation is given by an infinite series of orthogonal functions. Such an approach, however, is effective only for obstacle shapes adapted to those coordinate systems in which the wave equation is separable. For this reason, the dynamic stress concentrations around circular and parabolic obstacles have received considerable attention in the past. A comprehensive survey of the literature in a field as wide and diversified as the propagation of elastic waves is clearly beyond the scope of this paper. In recent years, the Mow-Pao-Thau school [4]—[6] has * Received January 8, 1968; revised version received March 29, 1968. The research described in this paper was sponsored by the U. S. Navy under Contract Nonr-610(06) with the Office of Naval Research, Washington, D. C. 194 G. C. SIH AND J. F. LOEBER [Vol. XXVII, No. 2 published a number of papers on this subject. References to other work can be found in [4]-[6], It is well known that problems involving diffraction of plane harmonic, horizontally polarized shear waves (SH-waves) by a semi-infinite crack can be formulated in terms of integral equations, and solved by the Wiener-IIopf technique [7]. As pointed out by Sih [8], however, since the static limit of the semi-infinite crack solution is zero, it is not possible to estimate the precise magnification of the stresses due to dynamic effects. To overcome this shortcoming, Loeber and Sih [9] proposed to add another characteristic dimension into the problem, namely the crack width, and managed to obtain the exact behavior of the crack-front displacement and stress fields for the case of SH-waves diffracted by a finite or internal crack. Ang and Knopoff [10] have attempted to solve the internal crack problem earlier but their method yields results which are restricted to low frequencies and to distances far away from the crack. In elastodynamics, the farfield crack solution is not useful in the sense that it offers no information to the development of the theories of crack propagation. Generally speaking, the far-field solution can always be determined by the standard method of Wiener-Hopf [7] in a straightforward manner. On the other hand, considerable difficulty is encountered when the WienerHopf method is applied to find the near-field solution. One of the difficulties arises from the factorization of certain functions into functions analytic in the upper and lower half planes. The problem of the diffraction of electromagnetic waves' incident upon a slit has also been treated by Schmeltzer and Lewin [11] using the function-theoretic approach. Their results are left in terms of several complicated integrals the evaluation of which becomes a problem in itself, particularly in seeking the analytical form of the solution in the vicinity of the slit. Having discussed the previous work related to crack problems of SH-waves, it is natural to follow the discussion with a few remarks concerning the diffractions of plane harmonic compression waves (P-waves) and vertically polarized shear waves (SY-waves) by a line crack. Although both Miles [12] and Papadopoulos [13] have investigated crack problems of this type, their work discusses only the qualitative character of the displacement potentials without any explicit information given as to the nature of the local stress distribution. The mathematical description of these problems is somewhat complex because the scattered waves, caused by the line crack, are composed of both compression and shear waves even though the input wave may be of one type, either the Por SVwaves. For this and other reasons, the near-field solution of waves scattered by a crack with finite width is yet to be found. The purpose of this paper, aside from obtaining the stress solution close to the crack point, is to offer a method of solution for solving diffraction problems involving Pand SY-waves incident upon a line of discontinuity. The method can handle different types of boundary conditions2 on the line of discontinuity. For illustration, only the case of a traction-free crack will be considered. An important conclusion is that within certain ranges of wave lengths the dynamic stress distribution around the crack is quite sensitive to changes in the wave number. This is displayed graphically for different values of the irThe scattering of plane-polarized electromagnetic waves by a screen in a fluid medium is mathematically analogous to the SH-wave crack problem in elastodynamics. 2By following the steps outlined in this paper, it is clear that the problems of a rigid and rigidsmooth strip can be solved in the same way. 1969] WAVE PROPAGATION IN AN ELASTIC SOLID 195 Poisson's ratio. The knowledge gained in this investigation is believed to add further impetus to the understanding of the propagation of cracks under fluctuating loads. Field equations and input waves. Consider the propagation of elastic waves, produced by the action of oscillating compressional and shear forces, which vary harmonically in time and are applied in the .xy-plane containing a through crack. In the plane, there arise both compressional and shear waves, and the resulting displacements can be expressed in terms of two scalar functions
and *p each of which depend upon x, y, and t. The rectangular components of the displacement vector are ux = d/dx + Sip/dy, uy = d/dy — d\\p/dx, (1) Substituting Eq. (1) into the equations of motion under the conditions of plane strain, the following wave equations on (/> and d2 2(d~\\p d2\\p\\ d2 \\p .. C\\^? + W2) = ^• c*\\d? + w) = d?\" () In Eq. (2), cl and c2 stand, respectively, for the velocities of compression (irrotational) and shear (equivoluminal) waves in an infinitely extended elastic medium; they are given by c, = [(X + 2m)/p]1/2, c2 = (m/p)1/2 (3) with p being the mass density. As usual, in the case of generalized plane stress the Lame constant X in Eq. (3) is to be replaced by 2\\/i/(X + 2/x), while the shear modulus of elasticity u remains unchanged. From the stress and displacement relations, it is found that (d^ ay \\ \\5x'2 dx dy/ ' vxx = XV > + 2fx +2\" @ £i) ■ « = (o _ ?Jk i