An analytical model is proposed in this paper to predict the local bond strength τf by incorporating the heterogeneity at interface regions for deformed reinforcing bars centrally anchored in concrete. The rib width on bar surface is introduced as an interfacial characteristic parameter G in the proposed model accounting for the heterogeneity. Both the τf and local interfacial fracture energy GIIf from each specimen are linked to the G and determined analytically if the maximum pull-out loads Fmax are given from the test. It is found that the predicted τf is larger than the maximum average bond stress τavg-max. The discrepancy between the two values is reduced with the increasing of L/G. Moreover, as the L/G increases, the predicted τf shows a certain decrease and the reduction becomes smaller with stronger interfacial homogeneity. The predicted GIIf is significantly increased because of the weaker boundary effect. The validity of the proposed model is verified upon the comparisons between the predicted Fmax using the determined τf and GIIf and the experimental ones with the only failure mode of bar pull-out. Moreover, the bars here can be steel or FRP (fiber-reinforced polymer) bars, and the concrete refers to all types of cementitious materials.