The purpose of this article is to study the axial stress distribution in the interference-fit bolted joints. Firstly, by integrating the complex potential function theory with the Coulomb friction criterion, an analytical model for interfacial friction and a through-thickness compressive stress distribution equation within the clamping zone were established. Secondly, a finite element model of the interference-fit bolted joint plate was developed. Finally, the validity of the finite element model was verified against the analytical solutions, followed by a systematic investigation into the effects of interference level, preload, friction coefficient, and Poisson's ratio on stress distribution. Results indicate that within the linear elastic regime, both interference level and preload significantly govern the compressive stress distribution, the peak of axial stress diminishes with the increasing of interference or reduced preload, the friction coefficients paly important role on the stress distribution on the upper surface of the plate, the uniformity of stress around hole varying depths is enhanced by decreasing Poisson's ratio.