Elastohydrodynamic lubrication (EHL) is a common mode of fluid-film lubrication in which many machine elements operate. Its thermal behavior is an important concern especially for components working under extreme conditions such as high speeds, heavy loads, and surfaces with significant roughness. Previous thermal EHL (TEHL) studies focused only on the cases with smooth surfaces under the full-film lubrication condition. The present study intends to develop a more realistic unified TEHL model for point contact problems that is capable of simulating the entire transition of lubrication status from the full-film and mixed lubrication all the way down to boundary lubrication with real machined roughness. The model consists of the generalized Reynolds equation, elasticity equation, film thickness equation, and those for lubricant rheology in combination with the energy equation for the lubricant film and the surface temperature equations. The solution algorithms based on the improved semi-system approach have demonstrated a good ability to achieve stable solutions with fast convergence under severe operating conditions. Lubricant film thickness variation and temperature rises in the lubricant film and on the surfaces during the entire transition have been investigated. It appears that this model can be used to predict mixed TEHL characteristics in a wide range of operating conditions with or without three-dimensional (3D) surface roughness involved. Therefore, it can be employed as a useful tool in engineering analyses.