In order to fully understanding the distribution of residual stress after riveting and the relationship between residual stress and riveting process parameters during riveting, Finite Element Method was used to establish a riveting model. Quasi-static method to solve the convergence difficulties was adopted in riveting process. The riveting process was divided into six stages according to the stress versus time curves. The relationship of residual stress with rivet length and rivet hole clearance were established. The results show numerical simulation is effective for riveting process and can make a construction for the practical riveting.