A promising direction in the tailored control of quantum materials is the modification of their structure-function relationship via engineered laser-driven shock waves. In order to inform how the laser excitation should be tailored to optimize the effect of the shock waves, a model is desired for strongly focused nonlinear surface waves in crystals. A model was developed previously in the paraxial approximation for weakly diffracting nonlinear surface wave beams in isotropic solids [Shull et al., JASA 97, 2126 (1995)]. Harmonic generation and shock formation were investigated numerically in that work for sources with uniform amplitude distributions, with analytical solutions for the fundamental and second harmonic obtained for weakly focused Gaussian beams. Here, numerical results are reported for shock formation in focused nonlinear surface waves obtained using a model unrestricted by the paraxial approximation in order to accommodate strong focusing in isotropic solids. Removal of the restriction to isotropic solids is work in progress based on representation of the wave field in terms of its angular spectrum and a model for planar nonlinear surface waves in crystals [Hamilton et al., JASA 105, 639 (1999)]. [Work supported by the IR&D Program at Applied Research Laboratories.]