On the basis of generating the stationary Gaussian bridge deck irregularity by the inverse fourier transform method, the relationship between the phase angle and the high-order central moment of the random excitation process was established, and then the non-Gaussian characteristics of the bridge deck irregularity were changed by continuously adjusting the phase angle, the random excitation of non-Gaussian bridge deck irregularity with given power spectral density, skewness and kurtosis was generated by iteration.At the same time, based on the measured data of deck roughness, the stochastic process of deck roughness was reconstructed by phase modulation method.Numerical examples and reconstruction results of measured bridge deck irregularities show that the non-Gaussian bridge deck irregularities generated by phase modulation can meet the requirements of given target power spectral density, kurtosis and skewness, and the simulation accuracy is high.The area of amplitude distribution and the maximum value of amplitude increase with the increase of kurtosis, and the amplitude distribution area of the positive or negative deviation of the deck irregularity increases with the increase of the skewness.The engineering application results show that the vehicle-bridge vibration response excited by super-Gaussian bridge deck irregularity is the largest, followed by Gaussian bridge deck irregularity, and the vehicle-bridge vibration response excited by sub-Gaussian bridge deck irregularity is relatively small.The sensitivity of vehicle vibration to non-Gaussian bridge deck irregularity excitation is greater than that of bridge vibration.