Abstract: Method for simulating chatter-free milled surface topography includes: performing dynamics modeling on milling system and building delay differential equation with multiple delays to represent dynamic model; constructing state transition matrix between two adjacent cutting tool rotation periods by extending GRK method; predicting stable milling parameter zones based on Floquet theory; calculating vibration displacements on discretized time nodes during one cutting tool rotation period by means of fixed mapping point theorem; constructing motion trajectories of cutting tool edges in normal and feed directions with regard to milled surface of workpiece; predicting surface topography of workpiece by using spline interpolation to densify motion trajectories of cutting tool edges that participate in generating final surface of workpiece; calculating surface location error of milled surface based on motion trajectories of cutting tool edges in normal direction with regard to milled surface of workpiece; calculating