Spectral evolution equations are used to perform numerical studies of nonlinear surface acoustic waves in the (111) plane of several nonpiezoelectric cubic crystals. Nonlinearity matrix elements which describe the coupling of harmonic interactions are used to characterize velocity waveform distortion. In contrast to isotropic solids and the (001) plane of cubic crystals, the nonlinearity matrix elements usually cannot be written in a real-valued form. As a result, the harmonic components are not necessarily in phase, and dramatic variations in waveforms and propagation curves can be observed. Simulations are performed for initially monofrequency surface waves. In some directions the waveforms distort in a manner similar to nonlinear Rayleigh waves, while in other directions the velocity waveforms distort asymmetrically and the formation of shocks and cusped peaks is less distinct. In some cases, oscillations occur near the shocks and peaks because of phase differences between harmonics. A mathematical transformation based on the phase of the matrix elements is shown to provide a reasonable approximation of asymmetric waveform distortion in cases where the matrix elements have similar phase.