Abstract: To explore the seismic response law and failure mode of the slope containing faults. In this study, we took the Woqian Landslide as an example, a two-dimensional numerical model of "the slope containing faults and the upper hard and the lower soft" was constructed based on the geological prototype survey using discrete element software, and the Wenchuan Earthquake wave was used as the dynamic load input from the bottom boundary. The results show that the seismic dynamic response inside the slope is strong at the fault zone, and the amplification coefficient of peak acceleration can reach 2.5~3.5 times, while the amplification coefficient of peak acceleration inside the slope far away from the fault zone is generally not more than 2.5 times, and the amplification coefficient of peak acceleration inside the slope has elevation amplification effect. Moreover, the spectrum analysis reveals that there are multi-step peak frequencies in the spectral response in the slope, but the high-frequency amplification effect greater than 6.7 Hz is dominant in the footwall of the fault, while the low-frequency amplification effect of 2-4 Hz is dominant in the hanging wall of the fault. The dynamic response analysis of the "upper hard and lower soft" lithology combination reveals that the amplification factor of the peak acceleration of the upper limestone is greater than that of the lower slate, and the dynamic response at the lithologic interface is stronger than that of the upper hard and lower soft rock. The simulation of the slope failure movement process shows that, with the fault and lithologic interface as the boundary, the failure mode of the landslide can be summarized as three stages, namely, 0~20 s for vibration tension fracture, 20~40 s for starting the migration, and 40~80 s for accumulation and compaction.