Some numerical solutions for Lamb's problem

Abstract

We consider a version of Lamb's Problem in which a vertical time-dependent point force acts on the surface of a uniform half-space. The resulting surface disturbance is computed as vertical and horizontal components of displacement, particle velocity, acceleration, and strain. The goal is to provide numerical solutions appropriate to a comparison with observed wave forms produced by impacts onto granite and onto soil. Solutions for step- and delta-function sources are not physically realistic but represent limiting cases. They show a clear P arrival (larger on horizontal than vertical components) and an obscure S arrival. The Rayleigh pulse includes a singularity at the theoretical arrival time. All of the energy buildup appears on the vertical components and all of the energy decay, on the horizontal components. The effects of Poisson's ratio upon vertical displacements for a step-function source are shown. For fixed shear velocity, an increase of Poisson's ratio produces a P pulse which is larger, faster, and more gradually emergent, an S pulse with more clear-cut beginning, and a much narrower Rayleigh pulse. For a source-time function given by cos²(πt/T), −T/2 ≦ T/2, a × 10 reduction in pulse width at fixed pulse height yields an increase in P and Rayleigh-wave amplitudes by factors of 1, 10, and 100 for displacement, velocity and strain, and acceleration, respectively. The observed wave forms appear somewhat oscillatory, with widths proportional to the source pulse width. The Rayleigh pulse appears as emergent positive on vertical components and as sharp negative on horizontal components. We show a theoretical seismic profile for granite, with source pulse width of 10 µsec and detectors at 10, 20, 30, 40, and 50 cm. Pulse amplitude decays as r⁻¹ for P wave and r^−$12$ for Rayleigh wave. Pulse width broadens slightly with distance but the wave form character remains essentially unchanged.