Railway axles are among critical components ensuring safe and efficient train operations. They are particularly susceptible to damage mechanisms such as fretting wear and fatigue. Fretting induced by high contact pressure and microslip between contact surface can significantly deteriorate fatigue strength at the contact edge of the press-fit section. Recent research has been conducted to enhance axle strength and reliability. However, fretting wear or microcrack formation at the wheel-press-fit zone of axles is still an active area of investigation. Accurately analyzing fretting wear is challenging due to its sensitivity to numerous factors such as changes in friction coefficient, influence of wear particles, and selection of an appropriate wear model. This paper aimed to establish a comprehensive analysis method for fretting wear in interference-fitted axles using finite element analysis (FEA) and numerical analysis techniques. Two wear models were applied in simulations: an Archard wear model and an energy-based wear model. Analysis results were compared with experimental data from rotating bending fatigue press-fit specimens. This comparison will help validate the proposed analysis method and assess the effectiveness and accuracy of different wear models in predicting fretting wear in press-fit axles.