In this paper the fretting wear simulation technique with the press-fitted specimens have been developed, which can consider the running-in period of total wear process by adopting nonlinear wear coefficient. The amount of microslip and contact variable at press-fitted and at bending loaded condition of press-fitted shaft was analysed by applying finite element method. With the finite element analysis result, a numerical approach was applied to predict fretting wear based on modified Archard's equation and updating the change of contact pressure caused by local wear with influence function method. The predicted wear profiles of press-fitted specimens at the contact edge were compared with the experimental results obtained by rotating bending fatigue tests. It is shown that the predicted wear profiles considering the running-in period with nonlinear wear coefficent is consistent with experimental results than that with constant wear coefficient. Therefore, the fretting wear simulation technique proposed is feasible and efficient for numerical simulation of fretting wear on press fits at the initial stage of fatigue life.