. An analytical model is developed for the prediction of residual stresses in burnishing. The model is simplified as a concentrated force pressing on elastic-plastic half-space using the solution to the Boussinesq-Flament problem. The treated material admits the elastic-plastic properties with hardening using a power law constitutive relation. Trial computation using Johnson-Cook model on AISI 1042 steel is presented and the results are verified with the experimental results given by Bouzid’s previous work. The residual stresses in the feed direction show the same trend with the experimental results while some differences still exist near the surface because of the concentrated normal force assumption and such stresses increase with the increase of burnishing force, decrease with the increase of depth and turn to zero beyond the plastic deformation boundary.