The numerical solution of thermal elastohydrodynamics of grease lubrication in line contact is calculated based on combined solution of grease lubrication control equations from Ostwald′s model and energy equations by pressure-temperature calculation cycle.The effect of normal running parameters on pressure distribution and shape of oil film is analyzed,and the effect of thermal effect on the solutions is discussed under different running conditions.