本发明属于滚磨光整加工领域，提供了一种用于模拟滚磨光整加工工件表面残余应力的方法，解决了现有方法仅适用于喷丸强化这类颗粒运动特征相对简单的情况，不能模拟滚磨光整加工这种颗粒运动复杂的工况的问题。本方法将离散元和有限元相结合，利用Python搭建离散元软件EDEM和有限元软件ABAQUS的桥梁，将EDEM仿真得到的滚抛磨块对工件的作用特征提取出来并导入ABAQUS有限元软件中进行残余应力模拟。该方法利用Python工具库将离散元和有限元相结合，解决了ABAQUS中粒子生成器*Particle难以模拟滚磨光整加工复杂运动工况的问题，能够有效反映滚磨光整加工后工件表面残余应力的分布状态，并揭示滚抛磨块对工件表面的作用机制和残余应力的演变规律。