Abstract: The numerical simulation method was developed based on the Archard method and the Arbitrary Lagrangian-Eulerian adaptive meshing (ALE adaptive meshing) within the Abaqus. The evolution of the surface profile with fretting cycles was predicted based on the numerical method, and the influences of fretting wear on the contact pressure, the frictional shear stress and the slip amplitude were investigated in detail. The results show that the evolution of the surface profile can be predicted accurately by the numerical method. It is found that the peak value of contact pressure at the contact rim increases and the location of it moves inward towards the contact zone gradually due to the fretting wear. The location of the maximum frictional shear stress shifts from the stick-slip boundary to the worn-unworn boundary. Meanwhile, the width of open regime and the slip amplitude on the slip regime increase with increasing cycles.