This study presents a numerical model for the thermal-elastohydrodynamic lubri-cation of heterogeneous materials in impact motion,in which a rigid ball bounces on a starved non-Newtonian oil-covered plane surface of an elastic semi-infinite heterogeneous solid with inhomogeneous inclusions.The impact-rebound process and the microscopic response of the subsurface inhomogeneous inclusions are investigated.The inclusions are homogenized accord-ing to Eshelby\'s equivalent inclusion method.The Elrod algorithm is adopted to determine the lubrication starvation based on the solutions of pressure and film thickness,while the lubricant velocity and shear rate of the non-Newtonian lubricant are derived by using the separation flow method.The dynamic response of the cases subjected to constant impact mass,momentum,and energy is discussed to reveal the influence of the initial drop height on the impact-rebound process.The results imply that the inclusion disturbs the subsurface stress field and affects the dynamic response of the contact system when the surface pressure is high.The impact energy is the decisive factor for the stress peak,maximum hydrodynamic force,and restitution coefficient,while the dynamic response during the early approaching process is controlled by the drop height.