A structural dynamic parameter identification method aided by a rPCK surrogate model comprises the following steps. Establish a finite element model that roughly reflects the structural system to be analyzed. Establish the dynamic parameter space sample set. The structural system response space sample set driven by the dynamic parameter space sample set is established by using the probabilistic finite element analysis. The robust polynomial Chaos Kriging surrogate model is obtained by mapping the dynamic parameter space sample set to the structural system response space sample set. The measured structural system response is used to drive the rPCK surrogate model, and then Bayesian inference is used to identify the structural dynamic parameters. The mean value of Bayesian posterior estimation is used as the estimated value of structural dynamic parameters. The proposed method creates conditions for establishing a high-fidelity finite element model of the actual engineering structural system.