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.