Abstract: The vertical coupling vibration of beam-pile-soil system is investigated. First, the piles are assumed to be viscoelastic rods, and the plane strain model is employed to simulate the dynamic interactions between piles and surrounding soil layers. The dynamic beam-pile interactions are simplified as vertical point loads. On this basis, the matrix equation for solving the governing equations of motions for beam and piles is constructed in the frequency domain. The quasi-analytical solution in the time domain is then obtained by the discrete inverse Fourier transform. A comparison with numerical simulation is conducted to verify the rationality of the present solution. Finally, the effects of the geometric parameters of beam and pile defect on the dynamic response of beam-pile-soil system are discussed.