Abstract: The immersed tunnel is buried in the shallow surface of seabed, which is quite different from the land tunnel. It is of great scientific significance and engineering application value to study the seismic stability of immersed tunnel. Within the framework of Biot's dynamic consolidation equations, a viscoelastic-plastic stress-strain hysteresis curve for soil was constructed, and a real-time coupled shear-volumetric strain increment model was introduced. Based on the FLAC3D computational platform, an effective stress method for the sand liquefaction process was implemented, and the validity of this method was verified. Taking a typical borehole in Bohai sea as site condition, the seismic stability of immersed tunnel in sandy seabed was numerically simulated, and the uplift mechanism of immersed tunnel in liquefied seabed was studied. The results show that the decrease of the friction resistance of the side wall, the increase of the buoyancy of the bottom and the large deformation of seabed together lead to the uplift of the immersed tunnel. The asymmetric acceleration of the input ground motion leads to the asymmetric response phenomenon on both sides of the tunnel, and then leads to the separation of vertical displacement time history on both sides.