This study presents a numerical model for the thermal-elastohydrodynamic lubrication 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 according 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.

本研究提出了一种用于非均质材料在冲击运动中的热弹性流体动力润滑的数值模型，其中一个刚性球在具有非均质夹杂物的弹性半无限非均质固体的贫油非牛顿油覆盖平面表面上弹跳。研究了冲击-回弹过程和次表面不均匀夹杂物的微观响应。夹杂物根据 Eshelby 的等效夹杂物方法均质化。采用Elrod算法根据压力和油膜厚度的解来确定润滑不足，而非牛顿润滑剂的润滑剂速度和剪切率采用分离流法推导。受到恒定冲击质量、动量、讨论了能量以揭示初始跌落高度对冲击-回弹过程的影响。结果表明，当表面压力较高时，夹杂物会干扰次表面应力场并影响接触系统的动态响应。冲击功是应力峰值、最大水动力和恢复系数的决定性因素，而早期逼近过程的动力响应受落差控制。