The finite element method (FEM) was used to study the flat cylindrical indentation fatigue behavior using a kinematic hardening model (A-F model). This study was motivated by the experimental work of the preceding paper [B.X. Xu and Z.F. Yue, J. Mater. Res.21, 1793 (2006)], in which there were obvious similarities in the behavior of conventional fatigue specimens and indentation fatigue specimens. It is proposed that the A-F model can predict the indentation fatigue behavior. Generally, the experimental behavior of the indentation fatigue testing can be explained by the FEM analysis. In addition, the effect of residual stress on the indentation depth per cycle was studied. The effect of friction between the indenter and the specimen and evolution of von Mises stress beneath the indenter was also investigated. Numerical results showed that the effect of friction on the indentation depth propagation can be neglected. Further analysis showed that the steady-state indentation depth per cycle increases with increasing compressive residual stress and decreasing tensile residual stress.