Abstract: In order to obtain the residual stress field induced by laser shock processing(LSP), the numerical simulation method was used. A 2-D and axisymmetric finite element analysis(FEA)model of single laser shock processing on 35CD4 thick parts was established. History of the energies of material during dynamic analysis and surface dynamic stresses at different times were analyzed to validate the reasonability of the total time of dynamic analysis. The effect of mesh refinement and spatial distribution models of the loading on the simulation results were discussed. The results show that the element length should be around 0.03mm to get convergent results. The predicted results for single LSP with round laser spot are consistent with the available experimental data.