Based on the linear contact thermal elastohydrodynamic grease lubrication numerical calculation model,a linear contact thermal elastohydrodynamic grease lubrication numerical calculation model considering solid particles was established by modifying the relevant parameters of a single spherical solid particle.The multigrid method was used to solve the lubrication indexes such as pressure,film thickness and average temperature rise of lubricating oil film.The pressure,film thickness and temperature rise distribution of lubricating oil film at different particle velocity,size radius and central position were obtained and compared.The results show that the solid particles in grease are prone to cause sudden changes in oil film pressure and temperature rise.With the movement of solid particles to the center of oil film and the increase of central velocity and radius of particles,the overall distribution of pressure,film thickness and average temperature rise move to the entrance area,especially the influence of particle radius on oil film pressure,film thickness and average temperature rise is particularly significant.Therefore,solid particles,especially those with relatively large radius,should be avoided as far as possible in practical work.