Abstract: The dynamic response of a space-axisymmetrical pile-fluid-saturated soil coupled system is studied. Firstly, based on the theory of elasticity and the Porous Media theory, the governing differential equations for the coupled system are presented. Then, the differential quadrature element method is developed. After that, the developed differential quadrature element method and the second-order backward difference scheme are applied to discretize the governing differential equations of the coupled system on spatial and temporal domains, respectively. Finally, the Newton-Raphson iteration method is used to derive the numerical solutions of the system with initial conditions, the dynamic response of the system and the effect of parameters are studied, and the validity of the analysis method is verified.