Abstract: The dynamic interaction of pile-soil-superstructure system is a complex process, especially in inclined liquefaction-induced lateral spreading site. Due to the permanent large deformation of the ground generated during the earthquake, nonlinear reactions such as sliding and cracking may occur between piles and soil. Therefore, the simulation of pile-soil interaction is crucial. In order to explore the influence of nonlinear contact between pile and soil on the dynamic response of the inclined liquefied site-group pile-supertructure system, the finite element numerical models considering the pile-soil interaction springs and the direct binding between the pile nodes and soil nodes were established based on OpenSees. The results indicate that when considering the Pyliq springs for pile-soil interaction, the amplitude of soil acceleration is slightly lower, the constraint of the pile foundation on the soil becomes significantly weaker, and the residual displacement of the soil increases. Meanwhile, models with Pyliq springs can better simulate the curvature response of piles, while models that directly bind pile nodes and soil nodes overestimate the curvature of the pile top, making it difficult to accurately estimate the most unfavorable position of pile bending resistance. The influence of pile-soil interaction springs on the dynamic response of the upper structure is relatively small.