Fretting fatigue occurs when two contact bodies undergo small oscillatory relative motion due to cyclic loading. It leads to fretting damage, which reduces the service life of the bodies in contact. This can be explained by the high stresses generated at the contact surface between the two bodies. Therefore, numerical analysis, such as Finite Element Analysis (FEA), would be useful to understand the fretting fatigue phenomenon and to investigate techniques to reduce its effect on fatigue lifetime of mechanical components. In this paper, FEA of fretting fatigue Aluminum specimen is carried out in order to study the stress distribution and crack initiation location. Two commercial FEA packages, namely ANSYS and ABAQUS are used to analyze the specimen. The stress distribution along the contact surface between the two bodies obtained using both codes is compared and analyzed.