Abstract: The liquefiable foundation of the subway station may bring severe seismic damages to it. Accordingly, a nonlinear model for the liquefiable foundation-subway station interaction is constituted to calculate the large liquefaction deformation by using the constitutive model for sand. The space-time transformation of the dynamic pore pressure, the displacement and the acceleration of the liquefiable foundation are analyzed. The ascent of the subway station and the dynamic settlement of the side foundation are investigated. The distributions of liquefaction area and the displacement vectors of the foundation around the subway station are analyzed. The research findings can be helpful for analyzing the earthquake damages of liquefiable foundation around the underground structures or the ground structures.