Shot peening (SP) process is a typical surface strengthening process for metal and metal matrix composites, which can significantly improve the fatigue life and strength. The traditional SP simulation model falls short as it only takes into account one or a few shots, proving insufficient for accurately simulating the entire impact process involving hundreds of shots. In this study, a random multiple shots simulation modeling methodology with hundreds of random shots is proposed to simulate the impact process of SP. In order to reduce the simulation error, the random function Rand of MATLAB is used to generate the shot distributions many times, and the shot distribution closest to the average number is selected and the three-dimension parametric explicit dynamics numerical simulation model is built using ABAQUS software. Orthogonal experiments are carried out to investigate the influences of shot diameter, incident impact velocity, and angle on the residual stress distribution, roughness, and specimen deformation. Results showed that the average relative errors of maximum residual compressive stress, roughness, and deformation of specimen between simulation model and experimental value are 30.99, 16.14, and 16.73%, respectively. The primary factors affecting residual stress and deformation is shot diameter, and the main factor affecting roughness is impact velocity.