A steady turbulent flow over a flat plate is one of the basic problems of convective heat transfer processes, with the key theoretical significance and wide engineering applications. In this paper, the thermal boundary layer for steady turbulent flows over a flat plate is divided into the laminar sublayer and turbulent core zone. The temperature profile in respective zone is described by the cubic polynomial and 1/5 power function, respectively. Accordingly, the energy equations are established for the thermal boundary layer based on the integral method, and the analytical solutions of the integro-differential equation groups are obtained by employing the fourth-order Runge-Kutta method. Compared with the experimental results measured by Blackwell, Moffat and Kays, respectively, it is indicated that the analytical solutions in this paper are correct, and these theoretical solutions also agree well with those obtained from the classical Prandtl-Taylor’s turbulent two-layer theoretical model.