The importance of accurate time-integration in the numerical modelling of P-wave propagation

The numerical dissipation characteristics of the Newmark and generalised-α time-integration schemes are investigated for P-wave propagation in a fully saturated level-ground sand deposit, where higher frequencies than those for S-waves are of concern. The study focuses on resonance, which has been shown to be of utmost importance for triggering liquefaction due to P-waves alone. The generalised-α scheme performs well, provided that the time-step has been carefully selected. Conversely, the dissipative Newmark method can excessively damp the response, changing radically the computed results. This implies that a computationally prohibiting small time-step would be required for Newmark to provide an accurate solution.