CONDENSED MATTER THIRD SERIES, VOLUME 51, NUMBER 22 1 JUNE 1995-II Computationally efficient representation for elastostatic and elastodynamic Green's functions for anisotropic solids V. K. Tewary Materials Reliability Diuision, National Institute of Standards and Technology, 325 Broadway, Boulder, Colorado 80303 (Received 5 July 1994) A computationally efficient representation of the three-dimensional elastostatic and elastodynamic Green's functions for anisotropic solids is derived by solving the Christoffel equation in terms of delta functions. The representation is also applicable to other wave equations. The method is applied to cal- culate the transient and the static displacement field due to a point source in infinite and semi-infinite an- isotropic cubic solids. For elastodynamic calculations in anisotropic solids, our representation saves the computational time by a factor of about 1000 over the conventional Fourier-Laplace representation. In the elastostatic case, the computational efficiency of our method is much more than the conventional Fourier representation but comparable to the methods of Barnett and Barnett and Lothe in specific cases. I. INTRQDUCTIGN We derive an integral representation for the elasto- dynarnic Green's function in terms of a5 function. In this representation, even in the general anisotropic 3D (three-dimensional} case, only a 1D integration needs to be done numerically. The integration does not involve oscillatory functions or singularities. This reduces the problem of numerical convergence. The CPU time re- quired for calculating the elastodynamic Green's function is only about, ' of that required for the conventional Fourier-Laplace representation. ' Qur technique is also applicable to other wave equations such as the elec- tromagnetic or acoustic wave equation. We also calculate the elastostatic Green's function by taking the static limit of the elastodynamic response of a solid. Qur method is much more efticient than the con- ventional Fourier representation. However, the compu- tational eKciency of our method is comparable to that of Barnett" for an infinite solid and to that of Barnett and Lothe for calculations of surface displacements in a semi-infinite solid. The elastodynamic Green's function is useful for calcu- lating physical properties of solids involving long- wavelength phonons. Presently, there is a strong in- terest in wave-form-based ultrasonics for nondestructive characterization of anisotropic materials. These tech- niques measure the response of a material to an elastic pulse. which is very well modeled in terms of the elasto- dynamic Green's function. The elastostatic Green's func- tion is useful for calculating stress distribution in solids containing defects and discontinuities. ' In 1attice stat- ics calculations and the computer simulation of lattice de- fects, the elastostatic Green's function is needed to fix the asymptotic limit of the lattice distortions. ' The Green's function for an isotropic solid is usually' calculated by using the Fourier-I. aplace representation in the wave-vector — frequency space. This is adequate for traditional materials which could be approximated as iso- tropic. Modern composite materials are highly aniso- tropic. In this case the conventional Fourier representa- tion requires a 4D integration over oscillatory functions for the elastodynamic Green's function and 3D integra- tion for the elastostatic Green's function. In many appli- cations, we need to calculate the Green's function for bounded solids (containing free surfaces of interfaces). The Fourier-Laplace representation of the Green's func- tion in such cases is poorly convergent and its evaluation is CPU intensive. A powerful technique for calculating the 3D elastostat- ic Green's function was suggested by Barnett. This tech- nique is very well suited for calculation of the Green's function for infinite solids. " For this case, the computa- tional e%ciency of our method is about the same as that of Barnett. An extension of Stroh's formalism has been developed by Barnett and Lothe for bounded solids, but it is applicable only to 2D problems. Qur representation retains its simplicity even for 3D bounded solids. For ex- ample, our expression for the Green's function for a semi-infinite solid satisfies the free-surface boundary con- 51 15 695