目的 提高高温合金GH3039激光冲击强化仿真建模的效率。方法 利用Python脚本语言对Abaqus进行二次开发，利用插件对高温合金GH3039激光冲击强化过程进行仿真分析。采用侧倾固定Ψ法，通过实验测量激光冲击强化后的残余应力，并对仿真结果进行验证，分析不同激光工艺参数作用下高温合金GH3039表面和深度方向残余应力的分布规律。结果 仿真插件界面简洁，操作性强，结果准确。在其他参数不变的情况下，残余压应力受到光斑尺寸的影响较大。相较于光斑直径为4、2 mm，在光斑直径为6 mm时，其中心位置残余压应力分别提高了4.3%、53%。随着光斑尺寸的增大，表面残余压应力增大，且变化梯度减小，深度方向的残余压应力增大。随着激光能量的增加，表面残余压应力增大，且变化梯度增大，残余压应力峰值位于中心区域附近，在激光能量为6、7、8 J时，残余压应力层的平均厚度分别为0.55、0.67、0.82 mm，深度方向残余压应力层增厚。随着冲击次数的增加，冲击区域表面残余压应力平均值高于单次冲击，且波动梯度增大，冲击1、2、3次后残余压应力层的平均厚度分别为0.55、0.71、0.85 mm，深度方向残余压应力层深度增大。结论 利用Python脚本语言对ABAQUS进行二次开发，提高了仿真建模的效率，可为快速预测不同激光工艺参数下高温合金GH3039残余应力的分布规律提供参考。

In the simulation process of laser shock peening of superalloy GH3039 with different laser process parameters, one or several laser parameters need to be frequently changed for repeated modeling, which reduces the efficiency of simulation. The method of redevelopment of Abaqus with Python script can realize the parametric modeling of simulation process and significantly improve the efficiency of simulation analysis. With superalloy GH3039 as the research object, a simulation plug-in for laser shock peening of superalloy GH3039 was established by Python language, and the correctness of the plug-in results was verified by experiments. Through the comparison between the experimental results and the simulation data, it was known that the plug-in simulation results had high accuracy. Based on this plug-in, the distribution of residual stress in superalloy GH3039 with different laser process parameters was studied, and the distribution of residual stress in surface direction and depth direction of superalloy GH3039 with different laser spot sizes, laser energy and laser impact times was analyzed. The spot size had a greater effect on the distribution law of residual stress in the surface direction. Compared with the maximum residual stress under the spot diameter of 4 mm and 2 mm, the maximum residual stress under the spot diameter of 6 mm increased by 4.3% and 53% respectively. With the increase of the spot size, the residual compressive stress in the surface direction increased and the variation gradient decreased, and the residual compressive stress in the depth direction increased. This was because the large spot transmitted the shock wave in the form of plane wave and small spot transmitted shock wave in the form of spherical wave, while plane wave attenuated slowly when transmitting laser energy. With the increase of energy, the peak value of laser shock wave pressure changed. On the premise that other process parameters remained unchanged, the superalloy GH3039 was subject to laser shock peening with 6 J, 7 J and 8 J laser energy respectively, and the peak pressure was 2.3 GPa, 2.5 GPa and 2.68 GPa respectively. The average value of residual compressive stress in the impact area increased with the increase of laser energy, and the fluctuation gradient of residual compressive stress increased. The peak value of residual compressive stress was near the central region. This was because the surface wave of the circular spot converged to the center of the laser impact area, causing local reverse plastic strain and reducing the residual compressive stress level. The average depth of the plastic layer subject to 6 J, 7 J and 8 J laser energy was 0.55 mm, 0.67 mm and 0.82 mm, and the depth of the plastic layer increased. The samples were impacted once, twice and three times at the same position respectively. The average value of residual compressive stress in the surface direction of multiple impact peening was higher than that of single impact peening, but the difference was small between the average value of residual compressive stress in the second impact and the third impact. This was because the material was hardened after laser impact peening. With the increase of impact times, the degree of hardening and the fluctuation gradient of residual compressive stress both became larger. With the increase of impact times, the residual compressive stress increased. The average depth of plastic layer impacted once, twice and three times was 0.55 mm, 0.71 mm and 0.85 mm. Multiple impacts can increase the depth of plastic layer.