A comprehensive method of thermoelastohydrodynamic (TEHD) lubrication analysis for dynamically loaded journal bearings is presented. An algorithm for mass conserving cavitation is included, and the effect of viscosity variation with the temperature is taken into account. The Reynolds equation in the film is solved using the finite element (FE) discretization. Thermal distortions as well as the elastic deformation of the bearing surfaces are computed using the FE method. The temperature of the lubrication film is treated as a time-dependent three-dimensional variable with a parabolic variation with respect to the film thickness. In order to compute the temperature of the film and its surrounding solid surfaces, a new heat flux conservation algorithm is proposed. An important element in this analysis is the consideration of thermal boundary layers for solids. It is known that the thermal transients on the film-solid interfaces and the dynamic loading have the same period (one cycle). However, beyond the thermal boundary layers, the time scale for thermal transient in the journal and bushing are several orders of magnitude greater than those for the oil film. The Fourier series approximates the instantaneous temperature fields in the solid boundary layers. In this way, the mean heat flux that passes into the solid can be computed and a steady-state heat conduction equation can be used to obtain thermal fields inside the solids. Finally, solving the complex problem of big-end connecting-rod bearing TEHD lubrication proves the efficiency of the algorithm. Oil film temperatures are found to vary considerably over the time and space.