A numerical model for non-linear coupled analysis of the seismic response of liquefiable soils