Shrink fit joining a solid shaft and a cylinder requires an accurate estimate of the residual contact pressure to transmit high powers, as in the case of gas turbines. Such torques require deformations of materials beyond their elastic limits. This paper presents an analytical model that analyses the stresses in a shrink fit assembly consisting of a solid shaft and a cylindrical hub operating in the elastic-plastic range. Assuming the hub to have a nonlinear work-hardening material behavior, the distribution of stresses as functions of the interference, and the effect of geometry on contact pressure and on interference are exhibited. To validate the analytic method, Finite Element Method was used.Nomenclature: