The aim of this paper is to propose a two-stage theory-based analytical method for the dynamic performance of pile groups in layered poroelastic saturated cross-anisotropic soils induced by moving loadings. Among them, the free-field vibrational analysis of saturated soils is performed by the analytical element-layer approach (ALEA) and Fourier transformation. Based on the free-field response, the boundary element (BE) solution for the soil resistance at the soil-pile interface is derived utilizing the two-stage theory. Simultaneously, the finite element (FE) solutions for the pile shaft resistance and deformation of pile groups are derived based on the Timoshenko beam theory. Finally, the FE-BE coupled dynamic equation for deformations and internal loadings of the soil-pile system is obtained. Thereafter, the reliability of the proposed method is validated by comparing with existing solutions and FE data from ABAQUS. Based on the derived solutions, a comprehensive parametric study is performed to examine the effects of loading amplitude, force speed, soft soil-layer stiffness, soil anisotropy, and pile length on the dynamic responses of pile groups.