Computer Calculations Of Explosion-Produced Craters

A fundamental goal of the Plowshare Program is to predict crater dimensions when an explosive of known yield is buried at a given depth in a given medium. In this paper a technique is presented that seems to calculate adequately, from first principles, the mound and cavity growth that occur during nuclear and high-explosive cratering events. The technique features a standard, numerical approach to high- intensity stress-wave propagation coupled with a unique model of material behavior in brittle failure. A schematic of stress-wave propagation is presented in Fig. 1. The equation of motion provides a functional relation between the applied stress field and the resulting acceleration of each point in the medium. Accelerations, when allowed to act over a small time increment [A]t, produce new velocities, velocities produce displacements, displacements produce strains, and strains produce a new stress field. Time is incremented by [A]t, and the cycle is repeated. The time increment [A]t is determined by an independent stability condition which requires the increment to be smaller than the time necessary far a compressional wave to travel across the smallest zone. We assume that failure occurs in a brittle material when either the tensile strength or the maximum allowable strain energy of distortion is exceeded. When either of these two conditions is violated, the code uses the cracked equation of state to determine the stress. A preshot testing program is used to determine the consolidated and cracked equations of state. In situ properties to be determined by field logging are density and elastic velocity (compressional and shear ve-