In order to improve the efficiency and accuracy of numerical simulation of laser shock peening(LSP)

a continuous explicit-dynamic impact simulation strategy was proposed based on the traditional simulation strategy. The explicit-dynamic analysis was used directly to simulate the multiple shot peening. And then the implicit-static analysis was used to obtain the stable residual stress field after being balanced. The three-dimensional finite element model of a plate was established by ABAQUS. Based on this strategy

the distribution of residual stress field after repeated impacts was studied. After Python post-processing

the simulated values were in good agreement with the experimental measurements. It was shown that the maximum residual stress was -212.5 MPa; the mean value was -216.7 MPa

and the error was 1.9% when the laser power density was 1 GW/cm2. With the increase of laser power density from 1 GW/cm2 to 4 GW/cm2

the depth of residual stress layer increased from 0.7 mm to 1 mm. The simulation accuracy was effectively improved on the basis of greatly improving the simulation efficiency

providing a simulation idea for large area LSP numerical simulation of large-scale structures.