Abstract: Based on the existing experiment results, the fretting fatigue contact geometry of a riveted aluminum specimen was studied using the finite element method. In the numerical model, the direct restrict method and the Stick-Slip friction model were introduced. A two-step tension load was suggested to realize the pre-tighten force. The contact states of the specimen’s up-surface were analyzed, and the results were in well agreement with the observations of the existing experiments. Also, the distributions of stress amplitude and the mean stress of the inner and outer contact edges on the specimen’s up-surface under fatigue loads were analyzed. By comparing the analyzed results, the influences of the contact friction coefficient, fastening forces and remote stress on crack initiation were discussed, which would be useful for the next experimental researches.