Abstract: To explore the influences of pile-soil system on site seismic response, the secondary development is carried out based on the finite element software ABAQUS, and an equivalent linearization model is embedded to characterize the dynamic nonlinear constitutive relationship of soils. A three-dimensional finite element seismic response analysis model with pile foundation is established. By analyzing the seismic response of the pile-soil interaction system, the numerical results are compared with the centrifugal model test ones to verify the feasibility of the numerical model. The results show that: (1) Under earthquake action, the bending moment value of the pile body is the highest at the middle and lower positions, and a "waist cinching" phenomenon appears near the surface. (2) As the burial depth of the soil layer decreases, the peak displacement gradually increases, and the maximum value appears at the surface. (3) The peak acceleration amplification factor of soils exhibits a trend of first decreasing and then increasing with the decrease of the burial depth. By calculating the cumulative absolute velocity amplification factor of different soil layers, it is shown that the cumulative sustained effects of seismic motion are the greatest at the surface.