Soil liquefaction due to the earthquake causes serious damages and engineering problems, such as the reduction of the soil strength, large settlement of the ground surface, the flow of liquefied soil and the uplift behavior to the underground structures, and the large deformation induced by the uplift force threatens the stability and safety of the structures. In this paper, a FE-FD coupled method is used in the simulation, the cyclic elasto-plastic constitutive model and the updated lagrangian formulation are applied to deal with the material and geometrical nonlinearity of liquefied soil. The results show that after the earthquake, the exceed pore water pressure will still exist for some time and the structure has an obvious vertical uplift displacement related to the liquefied area and the flow of liquefied soil. The uplift displacement will decrease as the thickness of the upper liquefiable soil layer is reduced. The results can be regarded as a guidance and reference for the design of the large underground structures.