SUMMARYThe problem of predicting the evolution of liquefied ground, modelled as a viscoplastic material, is addressed by combining a minimum principle for the velocity field, which characterizes such an evolution, and a time step integration procedure. Two different numerical schemes are then presented for the finite element implementation of this minimum principle, namely, the regularization technique and the decomposition-co-ordination method by augmented Lagrangian. The second method, which proves more accurate and efficient than the first, is finally applied to simulate the incipient flow failure and subsequent spreading of a liquefied soil embankment subject to gravity. The strong influence of liquefied soil residual shear strength on reducing the maximum amplitude of the ground displacement is particularly emphasized in such an analysis.