乏油与表面形貌：解锁纯滑动热弹流润滑的关键密码

乏油与表面形貌：解锁纯滑动热弹流润滑的关键密码一、引言1.1研究背景与意义在现代工业中，机械设备的高效稳定运行对生产效率和产品质量起着决定性作用。润滑作为减少机械部件摩擦与磨损、降低能耗、延长设备使用寿命的关键手段，一直是机械工程领域的研究重点。热弹流润滑（ThermalElastohydrodynamicLubrication，简称TEHL）作为一种特殊的润滑状态，在高负荷、高速等严苛工况下，能够在摩擦副表面形成具有承载能力的润滑油膜，有效减少金属表面之间的直接接触，从而降低磨损和摩擦系数，对于保障机械设备的可靠运行具有至关重要的意义，被广泛应用于航空航天、汽车制造、精密机械等众多领域。例如在航空发动机中，热弹流润滑确保了高速旋转的轴承和齿轮等部件的正常工作，减少了因摩擦和磨损导致的故障，提高了发动机的性能和可靠性；在汽车的变速箱和发动机系统中，热弹流润滑同样发挥着关键作用，有助于提升汽车的动力传输效率和燃油经济性。然而，在实际工况中，润滑油的供应情况和摩擦副表面的形貌特征往往会对热弹流润滑性能产生显著影响。乏油现象在机械设备的润滑系统中并不罕见，由于润滑油的泄漏、挥发、消耗过快或补充不及时等原因，都可能导致润滑区域内润滑油供应不足，形成乏油状态。乏油会使润滑膜厚度变薄，甚至出现局部干涸区域，进而加剧摩擦副表面的磨损，增加设备的能耗和故障率。表面形貌则是指物体表面的微观几何形状，任何零部件的接触表面都存在一定程度的粗糙度，且在工作过程中，由于磨损、塑性变形等因素，表面形貌会发生动态变化。这些微观几何特征直接影响着润滑油膜的形成、分布和承载能力，进而对热弹流润滑状态产生重要影响。深入研究乏油和表面形貌对纯滑动热弹流润滑的影响，对于提升工业设备的性能和可靠性具有重要的现实意义。从设备性能提升的角度来看，通过掌握乏油和表面形貌的影响规律，可以优化润滑系统的设计和维护策略。在润滑系统设计阶段，根据不同工况下乏油的可能性，合理选择润滑油的类型、供应量以及润滑方式，确保在各种情况下都能提供足够的润滑。在设备维护方面，能够更加准确地预测设备的磨损情况，及时采取措施进行调整和修复，避免因润滑不良导致的设备故障和停机，从而提高生产效率，降低维护成本。对乏油和表面形貌影响的研究也有助于开发新型的润滑技术和表面处理工艺。通过改进润滑油配方，提高其在乏油条件下的润滑性能；采用先进的表面处理技术，优化表面形貌，增强表面的耐磨性和润滑性，进一步提升设备的整体性能。从润滑理论发展的角度而言，乏油和表面形貌因素的引入，为热弹流润滑理论的完善提供了新的研究方向。传统的热弹流润滑理论通常假设润滑区域内润滑油充足且表面光滑，然而实际工况与这一假设存在较大差异。研究乏油和表面形貌对热弹流润滑的影响，可以使理论模型更加贴近实际情况，填补理论与实际之间的差距。这不仅有助于深入理解热弹流润滑的物理机制，还能够为润滑理论的发展提供更加坚实的实验和理论基础，推动润滑学科的不断进步。对乏油和表面形貌影响的研究还能够促进多学科的交叉融合。涉及到流体力学、材料科学、摩擦学、传热学等多个学科领域的知识，通过跨学科的研究方法，可以拓展研究思路，发现新的现象和规律，为解决复杂的工程润滑问题提供创新的方法和手段。1.2国内外研究现状1.2.1乏油对热弹流润滑影响的研究进展乏油对热弹流润滑的影响是润滑领域的重要研究课题，国内外学者对此开展了大量研究。国外方面，早期研究主要聚焦于乏油条件下润滑膜的形成与维持机制。学者们通过理论分析和数值模拟，建立了一系列模型来描述乏油状态下润滑油的流动和分布情况。这些研究初步揭示了乏油导致润滑膜厚度减小的原因，为后续研究奠定了理论基础。在实验研究上，利用先进的测量技术，如光干涉技术、超声测量技术等，对乏油润滑膜的厚度和压力分布进行了精确测量，获得了许多重要的实验数据，进一步验证和完善了理论模型。国内在乏油对热弹流润滑影响的研究起步相对较晚，但发展迅速。众多高校和科研机构积极参与其中，在理论、数值模拟和实验研究方面都取得了显著成果。理论研究上，一些学者针对国内工业设备的实际工况，对传统的热弹流润滑理论进行了改进和拓展，考虑了更多实际因素，如润滑油的流变特性、接触表面的微观几何形状等对乏油热弹流润滑的影响，使理论模型更加符合实际情况。数值模拟方面，运用计算流体力学（CFD）和有限元方法（FEM）等先进的数值计算手段，对乏油条件下的热弹流润滑过程进行了深入模拟，能够准确预测润滑膜的厚度、压力和温度分布等关键参数，为润滑系统的设计和优化提供了有力的技术支持。在实验研究中，自主研发了多种实验装置，开展了大量不同工况下的乏油热弹流润滑实验，获得了丰富的实验数据，为理论和数值模拟结果的验证提供了可靠依据。尽管国内外在乏油对热弹流润滑影响的研究上取得了诸多成果，但仍存在一些不足之处。在理论模型方面，虽然考虑了多种因素，但对于一些复杂的实际工况，如润滑油在高温、高压下的化学反应以及多相流等情况，模型的准确性和适用性还有待进一步提高。数值模拟中，计算精度和计算效率之间的矛盾依然突出，对于大规模、复杂几何形状的润滑系统，模拟计算的时间成本较高，难以满足实际工程快速设计和优化的需求。实验研究上，由于实验条件的限制，部分实验结果难以完全反映实际工况下的乏油热弹流润滑现象，实验数据的普适性和可靠性需要进一步验证和完善。1.2.2表面形貌对热弹流润滑影响的研究进展表面形貌对热弹流润滑的影响同样受到国内外学者的广泛关注。国外研究起步较早，在基础理论和实验技术方面处于领先地位。在基础理论研究中，从微观角度深入分析了表面粗糙度、纹理方向、波纹度等表面形貌参数对热弹流润滑性能的影响机制。通过建立微观接触力学模型，研究了表面微凸体之间的相互作用以及对润滑油膜形成和承载能力的影响，揭示了表面形貌与润滑性能之间的内在联系。实验技术上，采用原子力显微镜（AFM）、扫描电子显微镜（SEM）等先进的表面测量技术，能够精确测量表面形貌的微观特征，并结合光干涉热弹流润滑实验技术，实时观测表面形貌对润滑膜厚度、压力和温度分布的影响，为理论研究提供了直观的实验证据。国内在表面形貌对热弹流润滑影响的研究方面也取得了长足的进步。在理论研究中，结合国内的实际工程需求，对国外的理论模型进行了改进和创新，提出了一些具有自主知识产权的理论模型和计算方法。考虑了表面形貌的随机性和多尺度特性，以及表面形貌与润滑介质之间的耦合作用，使理论模型更加全面地反映实际情况。数值模拟方面，开发了一系列针对表面形貌的数值计算软件和算法，能够高效准确地模拟不同表面形貌下的热弹流润滑过程。通过数值模拟，深入研究了表面形貌参数对润滑性能的影响规律，为表面形貌的优化设计提供了理论依据。实验研究上，搭建了多种先进的实验平台，开展了大量表面形貌对热弹流润滑影响的实验研究。不仅对传统的表面形貌参数进行了研究，还对一些新型表面形貌，如仿生表面形貌、微纳结构表面形貌等进行了探索性研究，取得了许多有价值的研究成果。然而，目前表面形貌对热弹流润滑影响的研究仍存在一些需要完善的地方。在理论研究中，对于复杂表面形貌的数学描述和建模还不够完善，难以准确反映实际表面的多样性和复杂性。数值模拟中，对于表面形貌与润滑介质之间的强耦合作用的模拟还存在一定的误差，需要进一步改进模拟算法和模型。实验研究上，虽然能够测量表面形貌的微观特征，但对于表面形貌在实际工况下的动态变化及其对热弹流润滑性能的实时影响，还缺乏有效的监测手段和研究方法。1.3研究内容与方法1.3.1研究内容本研究聚焦于乏油和表面形貌对纯滑动热弹流润滑的影响，具体涵盖以下几个方面：乏油条件下热弹流润滑特性研究：深入探究不同乏油程度，即润滑油供应不足的不同状态，对润滑膜厚度的影响规律。通过理论分析和数值模拟，构建考虑乏油因素的润滑膜厚度计算模型，分析在不同工况下，如不同载荷、速度条件时，乏油导致润滑膜厚度的变化趋势，明确乏油程度与润滑膜厚度之间的定量关系。研究乏油对润滑膜压力分布的影响，揭示在乏油状态下，压力在接触区域内的重新分布机制，找出压力集中区域和压力变化的关键位置，为理解润滑膜的承载能力变化提供依据。探讨乏油对润滑膜温度分布的影响，考虑润滑油的热传递特性以及乏油时摩擦生热的变化，分析温度升高对润滑油粘度和润滑性能的影响，评估因乏油导致的温度异常升高可能引发的润滑失效风险。表面形貌对热弹流润滑特性的影响研究：针对表面粗糙度这一重要的表面形貌参数，研究不同粗糙度幅值和频率对热弹流润滑性能的影响。利用数值模拟方法，模拟不同粗糙度表面的润滑过程，分析粗糙度如何影响润滑油膜的形成和破裂，以及对摩擦系数、磨损率等性能指标的影响规律。探讨表面纹理方向对热弹流润滑的影响，研究当表面纹理方向与滑动方向一致、垂直或呈一定夹角时，润滑膜的压力、厚度和温度分布的变化情况，分析纹理方向对润滑油的流动特性和承载能力的影响机制，为表面纹理的优化设计提供理论指导。对表面波纹度进行研究，分析不同波纹度幅值和波长下热弹流润滑特性的变化，探讨波纹度如何与润滑油膜相互作用，影响润滑膜的稳定性和润滑效果，为表面形貌的综合评价和优化提供参考。乏油与表面形貌耦合作用下的热弹流润滑特性研究：考虑乏油和表面形貌同时存在时，二者之间的相互作用对热弹流润滑性能的影响。研究在乏油条件下，表面形貌的变化如何加剧或缓解乏油对润滑性能的不利影响；分析表面形貌对乏油区域内润滑油的流动和分布的调控作用，以及乏油如何改变表面形貌与润滑油膜之间的相互作用机制。建立考虑乏油和表面形貌耦合作用的热弹流润滑模型，通过数值模拟和实验验证，深入研究耦合作用下润滑膜的厚度、压力、温度分布以及摩擦系数、磨损率等性能指标的变化规律，为实际工程中复杂工况下的润滑设计提供理论依据和技术支持。1.3.2研究方法为实现上述研究内容，本研究将综合运用实验研究、数值模拟和理论分析三种方法：实验研究：搭建专门的热弹流润滑实验平台，该平台能够精确模拟纯滑动工况，具备加载系统以施加不同的载荷，运动控制系统以调节滑动速度，以及润滑系统用于控制润滑油的供应，从而实现对乏油状态的模拟。采用高精度的测量仪器，如光干涉膜厚测量仪，用于测量润滑膜的厚度，该仪器利用光的干涉原理，能够精确测量出纳米级别的膜厚变化；采用压力传感器测量润滑膜的压力分布，确保测量的准确性和可靠性。利用原子力显微镜（AFM）和扫描电子显微镜（SEM）等先进设备，对摩擦副表面形貌进行精确测量和表征，获取表面粗糙度、纹理方向、波纹度等参数，为后续的数值模拟和理论分析提供真实可靠的表面形貌数据。在不同的乏油程度和表面形貌条件下，进行多组热弹流润滑实验，获取润滑膜厚度、压力、温度以及摩擦系数、磨损率等实验数据。通过对实验数据的分析，直观地了解乏油和表面形貌对热弹流润滑性能的影响规律，验证数值模拟和理论分析的结果。数值模拟：基于计算流体力学（CFD）和有限元方法（FEM），建立考虑热效应、弹性变形、乏油和表面形貌的热弹流润滑数值模型。在模型中，充分考虑润滑油的流变特性，如粘度随温度和压力的变化关系，以及摩擦副材料的弹性特性，确保模型能够准确反映实际的热弹流润滑过程。利用商业软件或自主开发的程序，对建立的数值模型进行求解，模拟不同工况下热弹流润滑的过程，得到润滑膜的厚度、压力、温度分布等详细信息。通过改变模型中的乏油程度和表面形貌参数，系统地研究它们对热弹流润滑性能的影响，分析数值模拟结果，揭示热弹流润滑的内在物理机制，为实验研究和理论分析提供有力的支持。理论分析：基于经典的热弹流润滑理论，如Reynolds方程、弹性力学理论和传热学理论，建立考虑乏油和表面形貌影响的热弹流润滑理论模型。在理论模型中，引入适当的修正系数和假设，以考虑乏油对润滑油供应和流动的影响，以及表面形貌对润滑膜形成和承载能力的影响。运用数学分析方法，对理论模型进行求解，得到润滑膜厚度、压力、温度等参数的解析表达式或近似解。通过对理论解的分析，探讨乏油和表面形貌对热弹流润滑性能的影响规律，揭示热弹流润滑的基本物理原理，为实验研究和数值模拟提供理论指导，同时也为工程应用提供简化的计算方法和设计准则。二、相关理论基础2.1热弹流润滑理论2.1.1热弹流润滑的基本方程热弹流润滑理论是在考虑表面弹性变形和热效应的基础上，对流体动力润滑理论的进一步拓展，旨在揭示在高负荷、高速等复杂工况下，润滑油膜的形成、分布以及承载能力等特性。其基本方程涵盖了雷诺方程、能量方程、弹性变形方程等，这些方程相互关联，共同描述了热弹流润滑的物理过程。雷诺方程是热弹流润滑理论的核心方程之一，它基于流体力学的连续性方程和动量方程推导得出，用于描述润滑油在两接触表面间的流动和压力分布。在稳态、等温且不考虑惯性力和体积力的条件下，对于二维无限宽平行平板间的层流流动，雷诺方程的经典形式为：\frac{\partial}{\partialx}(\frac{\rhoh^3}{12\eta}\frac{\partialp}{\partialx})+\frac{\partial}{\partialy}(\frac{\rhoh^3}{12\eta}\frac{\partialp}{\partialy})=U\frac{\partial(\rhoh)}{\partialx}其中，x和y为坐标方向，\rho为润滑油密度，h为油膜厚度，\eta为润滑油动力粘度，p为油膜压力，U为两表面的相对速度。在热弹流润滑中，由于温度和压力的变化会显著影响润滑油的粘度和密度，因此需要对上述方程进行修正，以考虑这些因素的影响。引入粘压效应和温粘效应后，润滑油粘度\eta通常采用Roelands粘压温方程来描述：\eta=\eta_0\exp\left\{(\ln\eta_0+9.67)\left[-1+\left(1+\frac{p}{p_0}\right)^z\left(\frac{T_0}{T}\right)^s\right]\right\}其中，\eta_0为常压和参考温度T_0下的润滑油粘度，p_0为参考压力，z为粘压指数，s为温粘指数。润滑油密度\rho可通过Dowson-Higginson密压方程进行修正：\rho=\rho_0\left(1+\frac{0.6p}{1+1.7p}\right)其中，\rho_0为常压下的润滑油密度。考虑粘压效应和温粘效应后的雷诺方程，能够更准确地描述热弹流润滑中油膜压力的分布情况，为后续分析提供了重要的理论基础。能量方程用于描述热弹流润滑过程中的能量守恒关系，主要考虑了润滑油的对流换热、传导换热以及摩擦生热等因素。其一般形式为：\rhoc_p\left(u\frac{\partialT}{\partialx}+v\frac{\partialT}{\partialy}+w\frac{\partialT}{\partialz}\right)=k\left(\frac{\partial^2T}{\partialx^2}+\frac{\partial^2T}{\partialy^2}+\frac{\partial^2T}{\partialz^2}\right)+\eta\left(\left(\frac{\partialu}{\partialz}\right)^2+\left(\frac{\partialv}{\partialz}\right)^2\right)其中，c_p为润滑油的定压比热容，T为润滑油温度，u、v、w分别为润滑油在x、y、z方向的速度分量，k为润滑油的热导率。在热弹流润滑中，能量方程的求解对于准确分析油膜温度分布至关重要。由于油膜厚度通常远小于其他方向的尺寸，可对能量方程进行适当简化，如忽略y方向的导热和对流项，以及z方向的对流项等，从而得到适用于热弹流润滑分析的简化能量方程。通过求解能量方程，可以得到油膜温度分布，进而根据温粘关系确定润滑油粘度的变化，为雷诺方程的求解提供准确的粘度参数，实现热弹流润滑问题的耦合求解。弹性变形方程用于描述接触表面在油膜压力作用下的弹性变形情况。根据弹性力学理论，对于二维线接触问题，接触表面的弹性变形\delta(x)可通过Boussinesq公式计算：\delta(x)=\frac{2}{\piE'}\int_{-\infty}^{\infty}p(x')\ln\left|\frac{1}{x-x'}\right|dx'其中，E'为综合弹性模量，p(x')为油膜压力分布。在热弹流润滑中，接触表面的弹性变形会显著影响油膜厚度的分布，进而影响油膜压力和温度分布。因此，弹性变形方程与雷诺方程和能量方程相互耦合，需要进行迭代求解。通过不断更新油膜压力和弹性变形，直至满足收敛条件，才能得到准确的热弹流润滑解。在实际计算中，通常采用数值方法，如有限元法、边界元法等，对弹性变形方程进行离散求解，以提高计算精度和效率。除了上述主要方程外，热弹流润滑理论还涉及载荷平衡方程，用于确保油膜压力能够承受外部载荷，维持接触表面的平衡；以及膜厚方程，用于描述油膜厚度与接触表面形状、弹性变形以及相对运动等因素之间的关系。这些方程共同构成了热弹流润滑理论的基本框架，为深入研究热弹流润滑现象提供了坚实的理论基础。2.1.2热弹流润滑的求解方法由于热弹流润滑基本方程的非线性和强耦合特性，其求解过程较为复杂，通常需要借助数值方法来实现。目前，常用的求解方法包括多重网格法、有限元法、边界元法以及有限差分法等，每种方法都有其独特的优势和适用范围。多重网格法是一种高效的数值求解方法，它通过在不同尺度的网格上进行迭代计算，能够显著提高计算效率和收敛速度。在多重网格法中，首先在较粗的网格上求解方程，得到一个近似解，然后将该解作为初始值，在更细的网格上进行迭代细化，逐步逼近精确解。通过这种方式，能够有效地避免在单一细网格上求解时可能出现的计算量过大和收敛困难的问题。多重网格法在热弹流润滑求解中具有广泛的应用，尤其是对于大规模、复杂几何形状的润滑问题，其优势更为明显。在求解点接触或线接触热弹流润滑问题时，多重网格法能够快速准确地得到油膜压力、膜厚和温度分布等关键参数，为工程设计和分析提供了有力的支持。有限元法是一种基于变分原理的数值计算方法，它将求解区域离散为有限个单元，通过对每个单元进行分析，将连续的物理问题转化为离散的代数方程组进行求解。在热弹流润滑中，有限元法能够灵活处理各种复杂的边界条件和几何形状，对油膜的流动、传热以及接触表面的弹性变形等进行精确模拟。利用有限元软件，如ANSYS、ABAQUS等，可以方便地建立热弹流润滑模型，通过选择合适的单元类型和材料参数，对热弹流润滑过程进行数值模拟。有限元法的优点在于其计算精度高、适应性强，但计算量较大，对计算机硬件性能要求较高。对于一些复杂的热弹流润滑问题，如考虑表面形貌、多相流等因素时，有限元法能够提供较为准确的分析结果，但需要合理优化计算参数和网格划分，以提高计算效率。边界元法是一种基于边界积分方程的数值方法，它只需对求解区域的边界进行离散，从而降低了问题的维数，减少了计算量。在热弹流润滑中，边界元法主要用于求解弹性变形和油膜压力分布，通过将边界积分方程离散化，得到关于边界节点未知量的代数方程组，进而求解出边界上的物理量。边界元法的优势在于其对边界条件的处理较为简便，能够有效地处理无限域和半无限域问题。在求解热弹流润滑问题时，边界元法可以与其他方法，如有限元法相结合，充分发挥各自的优势，提高计算效率和精度。例如，在处理接触表面的弹性变形时，采用边界元法可以减少计算量，而在处理油膜的流动和传热时，采用有限元法可以提高计算精度。有限差分法是一种将求解区域划分为网格，通过对微分方程进行差分近似，将其转化为代数方程进行求解的方法。在热弹流润滑中，有限差分法常用于离散雷诺方程、能量方程和弹性变形方程等。该方法的优点是算法简单、易于实现，但对于复杂的几何形状和边界条件，其网格划分和差分格式的选择较为困难，可能会影响计算精度和稳定性。在早期的热弹流润滑研究中，有限差分法被广泛应用，随着计算机技术和数值算法的发展，虽然其在一些复杂问题的求解中逐渐被其他方法所取代，但在一些简单模型的分析和初步计算中，有限差分法仍然具有一定的应用价值。在实际应用中，为了获得更准确、高效的求解结果，往往会综合运用多种方法。结合多重网格法和有限元法，利用多重网格法加速有限元计算的收敛速度；或者将边界元法与有限差分法相结合，充分发挥边界元法在处理边界问题上的优势和有限差分法在计算效率上的长处。随着计算机技术的不断进步和数值算法的持续创新，热弹流润滑的求解方法也在不断发展和完善，为热弹流润滑理论的深入研究和工程应用提供了更强大的技术支持。2.2乏油润滑理论2.2.1乏油的定义与判定准则乏油是指在润滑过程中，由于润滑油供应不足，导致润滑区域内润滑油量低于维持良好润滑状态所需的最低量，从而使润滑性能下降的一种润滑状态。乏油现象在机械设备的运行中较为常见，如在航空发动机的高速旋转部件、汽车变速箱的齿轮啮合处等，都可能因润滑油的泄漏、挥发或补充不及时而出现乏油情况。当航空发动机长时间高负荷运行时，润滑油的消耗速度加快，如果润滑系统的供油能力不足，就容易导致轴承等部件处于乏油状态。判定乏油状态通常需要综合考虑多个参数，其中油膜厚度和供油率是两个关键指标。从油膜厚度角度来看，当实际测量得到的油膜厚度低于某一临界值时，可认为进入乏油状态。这一临界值与具体的工况条件、摩擦副材料以及润滑油的性能等因素密切相关。在高速轻载的工况下，临界油膜厚度可能相对较小；而在低速重载的工况下，为了保证良好的润滑效果，临界油膜厚度则需要更大。对于某一特定的机械部件，在正常润滑状态下，油膜厚度通常保持在一定的范围内，如在某实验中，当油膜厚度低于0.5μm时，摩擦系数明显增大，磨损加剧，此时可判定该部件处于乏油状态。供油率也是判定乏油状态的重要依据。供油率是指单位时间内供应到润滑区域的润滑油量与维持完全润滑所需润滑油量的比值。当供油率小于1时，说明润滑油供应不足，存在乏油的可能性。在实际应用中，通常将供油率小于0.8作为判定乏油状态的一个参考标准。在齿轮传动系统中，通过计算单位时间内进入啮合区域的润滑油量与理论上完全覆盖齿面所需润滑油量的比值，来确定供油率。当供油率低于0.8时，齿面间的油膜容易出现破裂，导致金属直接接触的概率增加，从而影响齿轮的正常运行和寿命。除了油膜厚度和供油率外，还可以通过观察润滑区域内润滑油的分布情况来辅助判定乏油状态。在乏油条件下，润滑油往往无法均匀地分布在摩擦副表面，会出现局部干涸或油膜不连续的现象。通过光学显微镜或其他可视化技术，可以直接观察到这些现象，从而进一步确认乏油状态的存在。利用荧光标记的润滑油，在显微镜下观察其在摩擦副表面的分布情况，当发现明显的荧光缺失区域，即表示该区域可能处于乏油状态。2.2.2乏油对润滑性能的影响机制乏油对润滑性能的影响是多方面的，主要体现在油膜厚度、摩擦系数和磨损等方面。从油膜厚度方面来看，乏油会导致润滑膜厚度显著减小。在正常润滑状态下，润滑油能够在摩擦副表面形成连续且具有一定厚度的油膜，有效地分隔两个表面，减少金属直接接触的概率。当出现乏油时，由于润滑油供应不足，无法及时补充油膜在运动过程中的消耗，导致油膜厚度逐渐变薄。在某乏油润滑实验中，随着供油率从1.0逐渐降低到0.5，油膜厚度从初始的1.2μm迅速下降到0.3μm，降幅达到75%。这是因为在乏油条件下，润滑油的流动和分布受到限制，难以在接触区域内形成稳定的压力分布，从而无法支撑起足够厚度的油膜。油膜厚度的减小使得摩擦副表面的微凸体更容易相互接触，增加了表面粗糙度对润滑性能的影响，进一步降低了润滑效果。乏油对摩擦系数的影响也十分显著。随着油膜厚度的减小，摩擦副表面的微凸体之间的接触面积和接触力增大，导致摩擦系数急剧上升。在边界润滑状态下，摩擦主要由表面微凸体之间的直接接触引起，此时摩擦系数较大。乏油使得润滑状态向边界润滑转变，从而增加了摩擦阻力。研究表明，在乏油状态下，摩擦系数可能会比正常润滑状态下增大数倍。在某滑动轴承的乏油实验中，正常润滑时摩擦系数约为0.05，当进入乏油状态后，摩擦系数迅速升高到0.2以上，这不仅增加了机械设备的能耗，还可能导致零部件的温度升高，影响设备的正常运行。摩擦系数的增大还会使机械系统的效率降低，如在齿轮传动系统中，由于乏油导致的摩擦系数增大，会使传递相同扭矩所需的驱动力增加，从而降低了传动效率。磨损是乏油对润滑性能影响的另一个重要方面。乏油状态下，油膜厚度的减小和摩擦系数的增大使得摩擦副表面的磨损加剧。表面微凸体之间的直接接触会产生较高的局部应力，导致材料表面发生塑性变形、疲劳剥落等磨损现象。磨损产生的颗粒还可能进入润滑油中，形成磨粒磨损，进一步加速表面的磨损过程。在乏油条件下，磨损率通常会比正常润滑状态下高出数倍甚至数十倍。在某机械零件的乏油磨损实验中，正常润滑时磨损率约为0.1mg/h，而在乏油状态下，磨损率迅速上升到1mg/h以上，这使得零件的使用寿命大大缩短，增加了设备的维护成本和停机时间。严重的磨损还可能导致零部件的尺寸精度下降，影响设备的性能和可靠性。乏油还会对润滑膜的压力分布和温度分布产生影响。在乏油条件下，润滑膜的压力分布会变得不均匀，出现压力集中的现象。这是因为润滑油供应不足，无法在整个接触区域内形成均匀的压力支撑，导致压力集中在局部区域，进一步加剧了表面的磨损和损坏。乏油还会使摩擦生热增加，由于润滑油的冷却作用减弱，无法及时带走摩擦产生的热量，导致润滑膜温度升高。温度的升高会使润滑油的粘度降低，进一步削弱油膜的承载能力，形成恶性循环，最终导致润滑失效。2.3表面形貌表征与分析2.3.1表面形貌的参数表征表面形貌是指物体表面具有的微观几何形状特征，其参数表征对于理解表面的物理性质和功能特性至关重要。粗糙度是描述表面微观不平度的重要参数，它反映了表面微观几何形状的不规则程度。常用的粗糙度参数包括轮廓算术平均偏差（Ra）、均方根偏差（Rq）等。Ra是指在一个取样长度内，轮廓偏距绝对值的算术平均值，其计算公式为：R_a=\frac{1}{l}\int_{0}^{l}|z(x)|dx其中，l为取样长度，z(x)为轮廓偏距，即表面轮廓上各点到基准线的距离。Rq则是指在一个取样长度内，轮廓偏距的均方根值，计算公式为：R_q=\sqrt{\frac{1}{l}\int_{0}^{l}z^{2}(x)dx}在实际测量中，通常使用轮廓仪、原子力显微镜（AFM）等设备来获取表面轮廓数据，进而计算粗糙度参数。轮廓仪通过触针扫描表面，将表面的微观起伏转化为电信号，经过数据处理得到表面轮廓曲线，再根据上述公式计算粗糙度参数。AFM则是利用微悬臂探针与表面原子间的相互作用力来扫描表面，能够获得更高分辨率的表面形貌信息，对于纳米级粗糙度的测量具有独特优势。波度是指表面上比粗糙度波长更长的几何形状误差，它对零件的运动精度、密封性等性能有重要影响。波度参数主要包括波高（Wz）、波距（Wl）等。波高是指在一个取样长度内，波峰与波谷之间的最大高度差；波距则是指相邻两个波峰或波谷之间的距离。在测量波度时，通常需要先对表面轮廓数据进行滤波处理，分离出波度成分，再进行参数计算。可以采用高斯滤波器、样条滤波器等对原始轮廓数据进行处理，根据ISO标准，通常将截止波长设置为1mm左右，将波长大于1mm的成分归为波度，小于1mm的成分归为粗糙度。经过滤波后的轮廓数据，即可用于计算波度参数，评估表面的波度特征。纹理方向是表面形貌的另一个重要特征，它描述了表面微观几何形状的方向性。纹理方向对润滑油的流动、接触表面的摩擦和磨损等都有显著影响。在一些机械加工表面，如磨削表面、铣削表面等，往往存在明显的纹理方向。纹理方向的表示方法有多种，常见的有基于傅里叶变换的方法和基于灰度共生矩阵的方法。基于傅里叶变换的方法通过对表面轮廓数据进行傅里叶变换，分析频谱图中能量分布的方向性来确定纹理方向；基于灰度共生矩阵的方法则是通过计算表面图像中不同方向上灰度值的共生概率，来提取纹理方向信息。在实际应用中，可根据具体情况选择合适的方法来表征纹理方向，为研究表面形貌对热弹流润滑的影响提供重要依据。除了上述参数外，表面形貌还可以通过其他一些参数进行表征，如表面支承率曲线、表面功能参数等。表面支承率曲线反映了表面在不同高度下的支承面积比例，对于分析表面的承载能力和磨损特性具有重要意义；表面功能参数则是根据表面的具体功能需求，定义的一些能够反映表面性能的参数，如用于表征表面润滑性能的储油坑体积、用于表征表面密封性能的泄漏率等。这些参数从不同角度全面地描述了表面形貌的特征，为深入研究表面形貌对热弹流润滑的影响提供了丰富的信息。2.3.2表面形貌的分析方法表面形貌的分析方法多种多样，每种方法都有其独特的优势和适用范围，轮廓算术平均偏差法是一种常用的表面形貌分析方法，它基于轮廓算术平均偏差（Ra）这一参数来评估表面粗糙度。通过测量表面轮廓上各点到基准线的距离，计算其绝对值的算术平均值，得到Ra值。该方法简单直观，能够快速地对表面粗糙度进行量化评价。在机械加工领域，常使用Ra值来衡量零件表面的加工质量，判断是否符合设计要求。在汽车发动机零部件的加工中，对活塞、气缸等表面的粗糙度有严格的要求，通过测量Ra值可以确保这些零部件的表面质量，保证发动机的正常运行。但轮廓算术平均偏差法也存在一定的局限性，它对表面轮廓的细节特征反映不够全面，对于一些具有复杂表面形貌的零件，可能无法准确地评估其表面质量。功率谱密度法（PowerSpectralDensity，简称PSD）是一种基于频域分析的表面形貌分析方法。它通过对表面轮廓数据进行傅里叶变换，将时域信号转换为频域信号，得到功率谱密度函数。功率谱密度函数描述了表面轮廓在不同频率成分上的能量分布情况，能够反映表面形貌的波长特性。通过分析功率谱密度函数，可以获取表面粗糙度、波度等信息，以及表面形貌的自相似性和分形特征。在研究材料表面的磨损过程中，功率谱密度法可以用于分析磨损前后表面形貌的变化，了解磨损机制。当材料表面受到磨损时，其表面形貌的频率成分会发生改变，通过对比磨损前后的功率谱密度函数，可以揭示磨损对表面形貌的影响规律。功率谱密度法还可以用于表面形貌的模拟和合成，为表面设计和优化提供支持。除了上述两种方法外，还有许多其他的表面形貌分析方法，如基于图像处理的方法、基于分形理论的方法等。基于图像处理的方法通过对表面图像进行灰度变换、边缘检测、特征提取等操作，来分析表面形貌的特征。利用图像的灰度值分布来计算表面粗糙度，或者通过边缘检测算法来识别表面的纹理方向和形状特征。这种方法具有直观、快速的特点，适用于对表面形貌进行定性分析。基于分形理论的方法则是将表面形貌看作是一种分形结构，通过计算分形维数等参数来描述表面的复杂程度和自相似性。分形维数越大，表面形貌越复杂，自相似性越强。在研究岩石表面的粗糙特性时，基于分形理论的方法可以有效地描述岩石表面的不规则性，为岩石力学性能的研究提供重要依据。在实际应用中，通常会根据具体的研究目的和表面形貌的特点，选择合适的分析方法。对于简单的表面形貌，轮廓算术平均偏差法可能就能够满足分析需求；而对于复杂的表面形貌，可能需要综合运用多种分析方法，如结合功率谱密度法和基于图像处理的方法，从不同角度全面地分析表面形貌的特征，为热弹流润滑研究提供更准确、详细的表面形貌信息。三、乏油对纯滑动热弹流的影响研究3.1实验研究3.1.1实验装置与材料实验采用的点接触试验机是整个研究的核心设备，其结构设计精密，能够准确模拟纯滑动工况。该试验机主要由加载系统、运动系统和润滑系统组成。加载系统配备高精度的载荷传感器，可实现0-1000N的连续加载，精度达到±0.1N，确保了施加的载荷能够精确控制并实时监测。运动系统由伺服电机驱动，通过高精度的滚珠丝杠和直线导轨，能够实现0.01-10m/s的稳定滑动速度，速度控制精度可达±0.001m/s，为研究不同速度下的热弹流润滑提供了可靠的条件。润滑系统采用了先进的微量供油装置，可精确控制润滑油的供给量，实现不同乏油程度的模拟，最小供油流量可低至0.01mL/min，保证了实验中对乏油工况的精确调控。实验选用的润滑油为某型号的矿物油，其在常温下的动力粘度为100mPa・s，密度为0.85g/cm³，具有良好的抗氧化性能和热稳定性，能够在较宽的温度范围内保持稳定的润滑性能。该润滑油的粘压系数为2.2×10⁻⁸Pa⁻¹，温粘系数为0.03/℃，这些参数对于研究热弹流润滑过程中润滑油的流变特性至关重要。为了进一步探究润滑油添加剂对乏油热弹流润滑的影响，还准备了添加了不同比例抗磨添加剂和极压添加剂的润滑油样本。抗磨添加剂能够在金属表面形成一层保护膜，减少磨损；极压添加剂则在高负荷条件下分解，与金属表面发生化学反应，生成具有高抗压强度的反应膜，提高润滑油的承载能力。试件的选择对于实验结果的准确性和可靠性同样关键。选用的试件材料为GCr15轴承钢，其具有高强度、高硬度和良好的耐磨性，洛氏硬度达到HRC60-62，弹性模量为210GPa，泊松比为0.3。试件的表面经过精密加工，粗糙度Ra控制在0.05-0.1μm之间，以确保实验过程中表面形貌的一致性和稳定性。在实验前，对试件进行了严格的清洗和脱脂处理，去除表面的油污和杂质，保证润滑效果不受外界因素干扰。为了研究表面形貌对乏油热弹流润滑的影响，还通过特殊的加工工艺制备了具有不同粗糙度和纹理方向的试件，包括通过磨削加工得到的不同粗糙度等级的试件，以及通过激光刻蚀制备的具有特定纹理方向的试件，为后续的实验研究提供了多样化的表面形貌条件。3.1.2实验方案设计本实验旨在全面探究乏油对纯滑动热弹流的影响，通过系统地改变油膜厚度、载荷、滑动速度等关键参数，精心设置了一系列不同乏油程度的工况，以获取全面且准确的实验数据。在油膜厚度方面，通过调节润滑系统的供油流量来实现不同油膜厚度的控制。设置了5个不同的供油流量水平，分别为0.05mL/min、0.1mL/min、0.2mL/min、0.4mL/min和0.8mL/min，对应的油膜厚度范围预计在0.05-0.5μm之间。根据前期的理论分析和预实验结果，这个油膜厚度范围能够涵盖从严重乏油到接近充分润滑的不同状态，为研究乏油对热弹流润滑的影响提供了丰富的工况条件。在每个供油流量下，进行多次重复实验，每次实验持续时间为30分钟，以确保实验数据的稳定性和可靠性。实验过程中，每隔5分钟记录一次油膜厚度、摩擦系数、温度等关键参数，通过对这些数据的分析，研究油膜厚度变化对热弹流润滑性能的影响规律。载荷参数设置了3个水平，分别为200N、500N和800N，以模拟不同的工作负荷条件。不同的载荷会导致接触表面的压力分布和变形情况发生变化，进而影响热弹流润滑的性能。在每个载荷水平下，结合上述不同的供油流量，进行全面的实验测试。在200N载荷下，研究不同乏油程度时，润滑膜的承载能力和稳定性；在500N和800N载荷下，进一步探究高负荷工况下乏油对热弹流润滑的影响，分析随着载荷增加，乏油状态对摩擦系数、磨损量等参数的影响趋势是否发生变化。对于每个载荷和供油流量的组合，同样进行多次重复实验，每次实验持续时间为30分钟，确保实验结果能够准确反映不同工况下的实际情况。滑动速度设置为0.1m/s、1m/s和5m/s三个水平，以研究不同速度下乏油对热弹流润滑的影响。滑动速度的变化会影响润滑油的流动特性和摩擦生热情况，从而对热弹流润滑性能产生重要影响。在0.1m/s的低速条件下，润滑油的流动相对缓慢，更容易出现乏油现象，研究此时乏油对润滑性能的影响，有助于了解低速工况下的润滑特性；在1m/s的中速条件下，分析乏油与速度的耦合作用对热弹流润滑的影响规律；在5m/s的高速条件下，由于摩擦生热增加，润滑油的粘度和性能会发生变化，研究此时乏油对热弹流润滑的影响，对于高速运转设备的润滑设计具有重要指导意义。在每个滑动速度下，与不同的载荷和供油流量进行组合实验，每次实验持续时间为30分钟，通过对实验数据的分析，揭示滑动速度与乏油程度对热弹流润滑性能的交互影响机制。为了确保实验结果的准确性和可靠性，每个工况下的实验均重复进行5次，取平均值作为最终实验结果。在每次实验前后，对试件表面进行仔细的清洗和检查，确保表面状态不受上一次实验的影响。实验过程中，实时监测并记录实验环境的温度和湿度，保证实验环境的稳定性。同时，对实验设备进行定期校准和维护，确保设备的各项性能指标符合实验要求，从而为实验数据的准确性提供有力保障。3.1.3实验结果与分析对不同乏油程度下的摩擦系数实验数据进行分析，结果表明，随着乏油程度的加剧，即供油流量的减小，摩擦系数呈现明显的上升趋势。在供油流量为0.8mL/min时，摩擦系数相对稳定，维持在0.05左右，此时润滑状态良好，润滑油能够在摩擦副表面形成连续且稳定的油膜，有效降低了表面间的摩擦。当供油流量减小到0.05mL/min时，摩擦系数急剧上升至0.2以上，这是因为乏油导致油膜厚度显著减小，表面微凸体之间的直接接触增多，从而使摩擦阻力大幅增加。在不同载荷和滑动速度条件下，这种趋势依然存在，但上升幅度有所差异。在高载荷和高速条件下，摩擦系数的上升更为显著。在800N载荷和5m/s滑动速度下，当供油流量从0.8mL/min减小到0.05mL/min时，摩擦系数从0.06迅速上升至0.35，这是由于高载荷和高速会加剧表面间的摩擦和磨损，而乏油进一步削弱了油膜的承载和润滑能力，导致摩擦系数急剧增大。磨损量的实验数据显示，乏油对磨损量的影响十分显著。随着乏油程度的增加，磨损量呈现指数增长。在充分供油（供油流量为0.8mL/min）时，磨损量较小，约为0.1mg/h，此时油膜能够有效地保护表面，减少磨损。当供油流量降低到0.1mL/min时，磨损量增加到0.5mg/h以上；当供油流量进一步减小到0.05mL/min时，磨损量急剧上升至2mg/h以上。这是因为乏油使得油膜无法有效隔离表面，微凸体之间的相互作用增强，导致表面材料的磨损加剧。不同载荷和滑动速度下，磨损量也有所不同。高载荷和高速会使磨损量进一步增大，且乏油对磨损量的影响在高载荷和高速条件下更为突出。在500N载荷和1m/s滑动速度下，当供油流量为0.8mL/min时，磨损量为0.2mg/h；当供油流量减小到0.05mL/min时，磨损量增加到1.5mg/h；而在800N载荷和5m/s滑动速度下，相同供油流量变化时，磨损量从0.3mg/h增加到3mg/h以上，表明高载荷和高速会加速乏油状态下的磨损过程。油膜厚度的实验结果表明，乏油会导致油膜厚度明显减小。随着供油流量的降低，油膜厚度逐渐变薄。在供油流量为0.8mL/min时，油膜厚度约为0.4μm；当供油流量减小到0.05mL/min时，油膜厚度降低到0.08μm以下。在不同载荷和滑动速度下，油膜厚度也会发生变化。高载荷会使接触表面变形增大，从而减小油膜厚度；高速则会使润滑油的剪切率增加，导致油膜厚度有所降低。在200N载荷和0.1m/s滑动速度下，供油流量为0.8mL/min时油膜厚度为0.45μm；在800N载荷和5m/s滑动速度下，相同供油流量时油膜厚度减小到0.35μm。乏油对油膜厚度的影响在高载荷和高速条件下更为显著，因为高载荷和高速会加剧润滑油的消耗和流失，使得乏油状态下油膜厚度的减小更加明显。温度的实验数据显示，随着乏油程度的加剧，润滑区域的温度显著升高。在充分供油时，温度相对较低，约为40℃。当供油流量减小到0.1mL/min时，温度升高到50℃以上；当供油流量进一步减小到0.05mL/min时，温度急剧上升至70℃以上。这是因为乏油导致摩擦系数增大，摩擦生热增加，同时润滑油的冷却作用减弱，无法及时带走热量，从而使温度升高。不同载荷和滑动速度下，温度变化也有所不同。高载荷和高速会使温度升高更快，且乏油对温度的影响在高载荷和高速条件下更为明显。在500N载荷和1m/s滑动速度下，当供油流量为0.8mL/min时，温度为45℃；当供油流量减小到0.05mL/min时，温度升高到75℃；而在800N载荷和5m/s滑动速度下，相同供油流量变化时，温度从50℃升高到90℃以上，表明高载荷和高速会加剧乏油状态下的温度上升，对润滑性能产生更为不利的影响。3.2数值模拟研究3.2.1数学模型建立基于热弹流润滑理论，建立考虑乏油影响的数学模型。雷诺方程作为描述润滑油膜压力分布的核心方程，在考虑乏油影响时，需对其进行修正。传统雷诺方程假设润滑油在整个接触区域内连续分布，而在乏油条件下，润滑油供应不足，导致油膜在部分区域不连续。引入油膜破裂函数F来考虑乏油对雷诺方程的影响，修正后的雷诺方程为：\frac{\partial}{\partialx}(\frac{\rhoh^3}{12\eta}\frac{\partialp}{\partialx})+\frac{\partial}{\partialy}(\frac{\rhoh^3}{12\eta}\frac{\partialp}{\partialy})=U\frac{\partial(\rhoh)}{\partialx}+\frac{\partial}{\partialx}(\frac{\rhoh^3}{12\eta}F\frac{\partialp}{\partialx})+\frac{\partial}{\partialy}(\frac{\rhoh^3}{12\eta}F\frac{\partialp}{\partialy})其中，F为油膜破裂函数，当处于充分润滑状态时，F=1；当出现乏油时，F的值根据乏油程度在0到1之间变化。通过这种方式，修正后的雷诺方程能够更准确地描述乏油条件下润滑油膜的压力分布。能量方程用于描述热弹流润滑过程中的能量守恒，在考虑乏油影响时，需考虑乏油导致的摩擦生热增加以及润滑油冷却作用减弱等因素。由于乏油使得油膜厚度减小，摩擦系数增大，摩擦生热相应增加。同时，润滑油供应量的减少导致其冷却能力下降，无法及时带走摩擦产生的热量。因此，在能量方程中增加摩擦生热项和冷却项的修正。摩擦生热项可表示为\mupU，其中\mu为摩擦系数，p为油膜压力，U为相对速度；冷却项的修正则根据润滑油的供应量和热传递特性进行调整。修正后的能量方程为：\rhoc_p\left(u\frac{\partialT}{\partialx}+v\frac{\partialT}{\partialy}+w\frac{\partialT}{\partialz}\right)=k\left(\frac{\partial^2T}{\partialx^2}+\frac{\partial^2T}{\partialy^2}+\frac{\partial^2T}{\partialz^2}\right)+\eta\left(\left(\frac{\partialu}{\partialz}\right)^2+\left(\frac{\partialv}{\partialz}\right)^2\right)+\mupU-\frac{q}{h}其中，q为单位面积的热流量，与润滑油的供应量和热传递系数有关，h为油膜厚度。通过这样的修正，能量方程能够更准确地反映乏油条件下的热传递过程。弹性变形方程用于描述接触表面在油膜压力作用下的弹性变形，在考虑乏油影响时，由于乏油导致油膜压力分布的变化，进而影响接触表面的弹性变形。根据弹性力学理论，接触表面的弹性变形\delta(x,y)与油膜压力p(x,y)之间的关系可通过Boussinesq公式计算：\delta(x,y)=\frac{2}{\piE'}\iint_{-\infty}^{\infty}p(x',y')\ln\left|\frac{1}{\sqrt{(x-x')^2+(y-y')^2}}\right|dx'dy'在乏油条件下，油膜压力p(x,y)由修正后的雷诺方程求解得到，将其代入弹性变形方程，即可计算出考虑乏油影响的接触表面弹性变形。通过迭代求解修正后的雷诺方程、能量方程和弹性变形方程，能够得到考虑乏油影响的热弹流润滑问题的数值解。边界条件的设定对于数学模型的求解至关重要。在入口边界，假设润滑油的压力为环境压力p_{in}，温度为环境温度T_{in}，油膜厚度为初始油膜厚度h_{in}。在出口边界，采用雷诺边界条件，即压力的二阶导数为零\frac{\partial^2p}{\partialx^2}=0，温度的二阶导数也为零\frac{\partial^2T}{\partialx^2}=0。在接触表面的两侧边界，假设压力和温度的法向梯度为零\frac{\partialp}{\partialn}=0，\frac{\partialT}{\partialn}=0，其中n为边界的法向方向。初始条件设定为润滑油膜在初始时刻均匀分布，压力和温度分布均匀，为后续的迭代求解提供初始值。3.2.2模拟结果与讨论通过数值模拟，深入研究乏油对压力分布、膜厚分布、温度分布的影响。从压力分布模拟结果来看，在充分润滑条件下，油膜压力在接触区域内呈现出典型的弹流润滑压力分布特征，即压力在接触中心区域较高，向两侧逐渐减小，形成一个近似椭圆形的压力分布区域。当出现乏油时，压力分布发生明显变化。随着乏油程度的增加，油膜压力在入口区域迅速下降，这是因为润滑油供应不足，无法在入口区域形成足够的压力支撑。压力分布的峰值向接触中心移动，且峰值压力增大。这是由于乏油导致油膜厚度减小，接触表面的微凸体之间的接触更加紧密，使得压力集中在更小的区域内，从而导致峰值压力升高。在严重乏油情况下，压力分布趋近于干接触状态下的赫兹压力分布，这表明乏油对油膜压力分布的影响非常显著，严重影响了润滑膜的承载能力。膜厚分布的模拟结果显示，乏油对膜厚的影响十分明显。在充分润滑时，油膜厚度在接触区域内相对均匀，且能够保持一定的厚度，有效地分隔接触表面，减少金属直接接触的概率。随着乏油程度的加剧，油膜厚度逐渐减小。在入口区域，由于润滑油供应不足，油膜厚度迅速变薄，甚至出现局部干涸区域。在接触中心区域，油膜厚度也明显降低，这是因为乏油导致油膜压力分布的变化，无法支撑起足够厚度的油膜。当乏油程度达到一定程度时，油膜厚度的最小值可能会降低到纳米级别，使得接触表面的微凸体之间直接接触的概率大大增加，从而加剧磨损和摩擦。通过模拟不同乏油程度下的膜厚分布，还发现膜厚分布的不均匀性随着乏油程度的增加而增大，这进一步说明了乏油对润滑膜稳定性的不利影响。温度分布的模拟结果表明，乏油会导致润滑区域的温度显著升高。在充分润滑条件下，润滑油能够有效地带走摩擦产生的热量，使得润滑区域的温度保持在较低水平。当出现乏油时，由于润滑油供应量减少，其冷却作用减弱，无法及时带走摩擦生热，导致温度迅速升高。在入口区域，由于油膜厚度减小，摩擦系数增大，摩擦生热增加，温度升高更为明显。随着乏油程度的增加，温度分布的最大值逐渐向接触中心移动，且最大值显著增大。在严重乏油情况下，润滑区域的最高温度可能会超过润滑油的闪点，导致润滑油氧化变质，进一步降低润滑性能，甚至引发润滑失效。温度升高还会导致润滑油粘度降低，从而削弱油膜的承载能力，形成恶性循环，加剧乏油对润滑性能的影响。通过对不同工况下的数值模拟结果进行对比分析，还发现乏油对压力分布、膜厚分布和温度分布的影响程度与载荷、滑动速度等因素密切相关。在高载荷条件下，乏油对压力分布和膜厚分布的影响更为显著，因为高载荷会使接触表面的变形增大，进一步加剧乏油对油膜的破坏作用。在高速条件下，乏油对温度分布的影响更为突出，由于高速运转会产生更多的摩擦热，而乏油使得散热困难，导致温度急剧升高。这些模拟结果为深入理解乏油对纯滑动热弹流的影响机制提供了重要依据，也为实际工程中的润滑系统设计和优化提供了有价值的参考。3.3乏油影响的理论分析从润滑机理角度深入剖析，乏油时润滑油供应不足会对热弹流润滑性能产生多方面的关键影响。在润滑油的流动与分布方面，当出现乏油状况时，由于单位时间内进入润滑区域的润滑油量减少，润滑油无法像充分润滑时那样顺畅地在摩擦副表面形成连续且均匀的油膜。在接触区域的入口处，由于润滑油供应不及时，油膜的连续性容易被破坏，出现局部的油膜变薄甚至干涸的现象。这会导致润滑油在接触区域内的流动阻力增大，难以有效地填充到接触表面的微凸体之间，从而影响润滑油膜的承载能力。在一些高速运转的机械部件中，如航空发动机的轴承，乏油会使润滑油在进入接触区域时无法迅速形成足够厚度的油膜，导致局部的金属直接接触，增加了磨损和摩擦的风险。从油膜承载能力来看，润滑油供应不足会使油膜的承载能力显著下降。在热弹流润滑中，油膜的承载能力主要依赖于润滑油在接触区域内形成的压力分布。当乏油时，由于润滑油无法充分填充接触区域，压力分布会发生改变。在入口区域，由于润滑油量不足，无法形成足够的压力梯度，导致压力迅速下降，无法有效地支撑外部载荷。随着乏油程度的加剧，压力分布的峰值向接触中心移动，且峰值压力增大，但整体的承载面积减小，使得油膜的平均承载能力降低。这就如同在一座桥梁上，原本均匀分布的支撑力量由于部分支撑缺失而变得集中在少数区域，虽然这些区域的支撑力可能会增大，但整体桥梁的承载能力却下降了。在实际的机械系统中，如汽车的变速箱齿轮，乏油会使齿面间的油膜承载能力下降，导致齿轮在传递动力时更容易出现疲劳磨损和胶合等失效形式。乏油还会对润滑油的热传递性能产生重要影响。在热弹流润滑过程中，润滑油不仅起到润滑作用，还承担着散热的功能。当润滑油供应不足时，其带走摩擦产生热量的能力减弱。摩擦生热会使润滑区域的温度升高，而高温又会导致润滑油的粘度降低，进一步削弱油膜的承载能力，形成恶性循环。由于乏油导致油膜厚度减小，润滑油与接触表面之间的热交换面积也减小，使得热量更难以散发出去。在一些高负荷、高速运转的机械设备中，如大型发电机组的轴承，乏油会导致温度迅速升高，如果不能及时解决，可能会引发轴承的烧蚀和损坏，严重影响设备的正常运行。乏油对润滑油的剪切特性也有影响。在乏油条件下，由于油膜厚度减小，润滑油所承受的剪切应力增大。当剪切应力超过润滑油的承受极限时，润滑油的分子结构会发生破坏，导致其润滑性能下降。润滑油的粘度会发生变化，不再符合理想的流变特性，进一步影响油膜的形成和稳定性。在一些精密机械中，如光学仪器的传动部件，对润滑油的剪切稳定性要求较高，乏油会使润滑油的剪切特性发生改变，从而影响仪器的精度和可靠性。四、表面形貌对纯滑动热弹流的影响研究4.1实验研究4.1.1实验装置与材料本实验采用的点接触试验机是专门为研究表面形貌对纯滑动热弹流影响而设计的，其结构紧凑，性能稳定，能够精确模拟实际工况中的纯滑动过程。试验机主要由加载系统、运动系统、润滑系统和测量系统组成。加载系统采用高精度的液压加载装置，可实现0-2000N的连续加载，加载精度达到±0.5N，能够满足不同载荷条件下的实验需求。运动系统由直流伺服电机驱动，通过精密的滚珠丝杠和直线导轨，实现了0.05-15m/s的稳定滑动速度，速度控制精度可达±0.002m/s，确保了滑动过程的稳定性和准确性。润滑系统配备了先进的微量喷油装置，可精确控制润滑油的供给量和供给方式，保证润滑油能够均匀地分布在接触表面，为实验提供稳定的润滑条件。测量系统是本实验装置的关键部分，它采用了多种先进的测量技术，以获取准确的实验数据。采用光干涉膜厚测量仪来测量润滑膜的厚度，该仪器利用光的干涉原理，能够精确测量出纳米级别的膜厚变化，测量精度可达±0.01μm。使用高精度的压力传感器来测量润滑膜的压力分布，传感器的测量范围为0-100MPa，精度达到±0.1MPa，能够实时监测压力的变化情况。为了测量表面形貌参数，配备了原子力显微镜（AFM）和扫描电子显微镜（SEM）。AFM能够对表面进行纳米级的三维成像，测量表面粗糙度等微观参数；SEM则可以提供高分辨率的表面图像，用于观察表面的微观结构和纹理特征。实验选用的润滑油为某高性能合成润滑油，其在常温下的动力粘度为80mPa・s，密度为0.88g/cm³，具有良好的抗氧化性能、抗磨损性能和热稳定性。该润滑油的粘压系数为2.5×10⁻⁸Pa⁻¹，温粘系数为0.025/℃，在不同的工况条件下能够保持较为稳定的润滑性能。为了研究润滑油添加剂对表面形貌热弹流润滑的影响，还准备了添加了不同比例抗磨添加剂和极压添加剂的润滑油样本。抗磨添加剂能够在金属表面形成一层保护膜，减少磨损；极压添加剂则在高负荷条件下分解，与金属表面发生化学反应，生成具有高抗压强度的反应膜，提高润滑油的承载能力。试件材料选用了45#钢，其具有良好的综合机械性能，硬度为HB220-250，弹性模量为200GPa，泊松比为0.3。试件的表面经过精密加工，粗糙度Ra控制在0.05-0.1μm之间，以确保实验过程中表面形貌的一致性和稳定性。在实验前，对试件进行了严格的清洗和脱脂处理，去除表面的油污和杂质，保证润滑效果不受外界因素干扰。为了研究不同表面形貌对热弹流润滑的影响，通过特殊的加工工艺制备了具有不同粗糙度和纹理方向的试件。利用磨削工艺制备了不同粗糙度等级的试件，粗糙度Ra分别为0.1μm、0.2μm、0.4μm和0.8μm；采用激光刻蚀技术制备了具有特定纹理方向的试件，纹理方向与滑动方向的夹角分别为0°、45°和90°，为后续的实验研究提供了多样化的表面形貌条件。4.1.2实验方案设计本实验旨在全面研究表面形貌对纯滑动热弹流的影响，通过系统地改变表面粗糙度、纹理方向等关键参数，精心设置了一系列不同表面形貌的工况，以获取全面且准确的实验数据。在表面粗糙度方面，设置了4个不同的粗糙度等级，分别为Ra=0.1μm、Ra=0.2μm、Ra=0.4μm和Ra=0.8μm，以研究粗糙度对热弹流润滑性能的影响。在每个粗糙度等级下，进行多次重复实验，每次实验持续时间为30分钟，以确保实验数据的稳定性和可靠性。实验过程中，每隔5分钟记录一次油膜厚度、摩擦系数、温度等关键参数，通过对这些数据的分析，研究粗糙度变化对热弹流润滑性能的影响规律。在Ra=0.1μm的粗糙度下，观察润滑膜的形成和稳定性；在Ra=0.8μm的粗糙度下，分析粗糙度对摩擦系数和磨损量的影响，以及对油膜厚度和温度分布的影响趋势。纹理方向参数设置了3个水平，分别为纹理方向与滑动方向夹角为0°、45°和90°，以模拟不同的纹理方向工况。不同的纹理方向会影响润滑油的流动特性和接触表面的摩擦状态，进而对热弹流润滑性能产生重要影响。在每个纹理方向下，结合上述不同的粗糙度等级，进行全面的实验测试。在纹理方向与滑动方向夹角为0°时，研究润滑油在纹理沟槽中的流动情况，以及对润滑膜承载能力的影响；在夹角为45°和90°时，分析纹理方向对摩擦系数、磨损量等参数的影响，以及对油膜厚度和温度分布的影响规律。对于每个纹理方向和粗糙度等级的组合，同样进行多次重复实验，每次实验持续时间为30分钟，确保实验结果能够准确反映不同工况下的实际情况。为了研究表面形貌与载荷、滑动速度的交互作用对热弹流润滑性能的影响，还设置了不同的载荷和滑动速度水平。载荷设置为300N、600N和900N三个水平，以模拟不同的工作负荷条件；滑动速度设置为0.2m/s、1m/s和5m/s三个水平，以研究不同速度下表面形貌对热弹流润滑的影响。在每个载荷和滑动速度水平下，与不同的表面粗糙度和纹理方向进行组合实验，每次实验持续时间为30分钟。在300N载荷和0.2m/s滑动速度下，研究表面形貌对热弹流润滑性能的影响；在900N载荷和5m/s滑动速度下，分析表面形貌与高载荷、高速工况的交互作用对热弹流润滑性能的影响，揭示表面形貌在不同工况下对热弹流润滑的影响机制。为了确保实验结果的准确性和可靠性，每个工况下的实验均重复进行5次，取平均值作为最终实验结果。在每次实验前后，对试件表面进行仔细的清洗和检查，确保表面状态不受上一次实验的影响。实验过程中，实时监测并记录实验环境的温度和湿度，保证实验环境的稳定性。同时，对实验设备进行定期校准和维护，确保设备的各项性能指标符合实验要求，从而为实验数据的准确性提供有力保障。4.1.3实验结果与分析不同表面粗糙度下的摩擦系数实验数据表明，随着表面粗糙度的增大，摩擦系数呈现明显的上升趋势。在表面粗糙度Ra=0.1μm时，摩擦系数相对较低，稳定在0.06左右，此时表面较为光滑，润滑油能够在表面形成连续且稳定的油膜，有效地降低了表面间的摩擦。当表面粗糙度增大到Ra=0.8μm时，摩擦系数急剧上升至0.2以上，这是因为粗糙表面的微凸体增多，表面间的实际接触面积增大，润滑油膜的连续性受到破坏，导致表面微凸体之间的直接接触增多，从而使摩擦阻力大幅增加。在不同载荷和滑动速度条件下，这种趋势依然存在，但上升幅度有所差异。在高载荷和高速条件下，摩擦系数的上升更为显著。在900N载荷和5m/s滑动速度下，当表面粗糙度从Ra=0.1μm增大到Ra=0.8μm时，摩擦系数从0.08迅速上升至0.35，这是由于高载荷和高速会加剧表面间的摩擦和磨损，而粗糙表面进一步削弱了油膜的承载和润滑能力，导致摩擦系数急剧增大。磨损量的实验数据显示，表面粗糙度对磨损量的影响十分显著。随着表面粗糙度的增加，磨损量呈现指数增长。在表面粗糙度Ra=0.1μm时，磨损量较小，约为0.15mg/h，此时油膜能够有效地保护表面，减少磨损。当表面粗糙度增大到Ra=0.4μm时，磨损量增加到0.6mg/h以上；当表面粗糙度进一步增大到Ra=0.8μm时，磨损量急剧上升至2.5mg/h以上。这是因为粗糙表面的微凸体在相对运动过程中更容易相互碰撞和挤压，导致表面材料的磨损加剧。不同载荷和滑动速度下，磨损量也有所不同。高载荷和高速会使磨损量进一步增大，且表面粗糙度对磨损量的影响在高载荷和高速条件下更为突出。在600N载荷和1m/s滑动速度下，当表面粗糙度为Ra=0.1μm时，磨损量为0.2mg/h；当表面粗糙度增大到Ra=0.8μm时，磨损量增加到1.8mg/h；而在900N载荷和5m/s滑动速度下，相同表面粗糙度变化时，磨损量从0.3mg/h增加到3.5mg/h以上，表明高载荷和高速会加速粗糙表面的磨损过程。油膜厚度的实验结果表明，表面粗糙度对油膜厚度有显著影响。随着表面粗糙度的增大，油膜厚度逐渐减小。在表面粗糙度Ra=0.1μm时，油膜厚度约为0.45μm；当表面粗糙度增大到Ra=0.8μm时，油膜厚度降低到0.1μm以下。在不同载荷和滑动速度下，油膜厚度也会发生变化。高载荷会使接触表面变形增大，从而减小油膜厚度；高速则会使润滑油的剪切率增加，导致油膜厚度有所降低。在300N载荷和0.2m/s滑动速度下，表面粗糙度为Ra=0.1μm时油膜厚度为0.5μm；在900N载荷和5m/s滑动速度下，相同表面粗糙度时油膜厚度减小到0.3μm。表面粗糙度对油膜厚度的影响在高载荷和高速条件下更为显著，因为高载荷和高速会加剧润滑油的消耗和流失，使得粗糙表面下油膜厚度的减小更加明显。温度的实验数据显示，随着表面粗糙度的增大，润滑区域的温度显著升高。在表面粗糙度Ra=0.1μm时，温度相对较低，约为42℃。当表面粗糙度增大到Ra=0.4μm时，温度升高到55℃以上；当表面粗糙度进一步增大到Ra=0.8μm时，温度急剧上升至75℃以上。这是因为粗糙表面导致摩擦系数增大，摩擦生热增加，同时润滑油膜厚度减小，其冷却作用减弱，无法及时带走热量，从而使温度升高。不同载荷和滑动速度下，温度变化也有所不同。高载荷和高速会使温度升高更快，且表面粗糙度对温度的影响在高载荷和高速条件下更为明显。在600N载荷和1m/s滑动速度下，当表面粗糙度为Ra=0.1μm时，温度为48℃；当表面粗糙度增大到Ra=0.8μm时，温度升高到80℃；而在900N载荷和5m/s滑动速度下，相同表面粗糙度变化时，温度从52℃升高到95℃以上，表明高载荷和高速会加剧粗糙表面下的温度上升，对润滑性能产生更为不利的影响。纹理方向对摩擦系数的影响实验结果表明，当纹理方向与滑动方向夹角为0°时，摩擦系数相对较低；随着夹角的增大，摩擦系数逐渐增大。在纹理方向与滑动方向夹角为0°时，润滑油能够在纹理沟槽中顺利流动，形成较为稳定的油膜，从而降低摩擦系数；当夹角为90°时，纹理对润滑油的流动产生较大阻碍，油膜的稳定性受到影响，导致摩擦系数增大。在不同粗糙度和载荷、滑动速度条件下，这种趋势依然存在，但变化幅度有所不同。在高粗糙度和高载荷、高速条件下，纹理方向对摩擦系数的影响更为显著。纹理方向对磨损量的影响实验数据显示，随着纹理方向与滑动方向夹角的增大，磨损量呈现先减小后增大的趋势。在夹角为45°左右时，磨损量相对较小，这是因为此时纹理对润滑油的引导作用和对表面微凸体的支撑作用达到较好的平衡，能够有效减少磨损；当夹角为0°或90°时，磨损量相对较大，因为在0°时，纹理对表面微凸体的支撑作用不足，而在90°时，纹理对润滑油的阻碍作用较大，都不利于减少磨损。不同粗糙度和载荷、滑动速度下，磨损量的变化趋势相似，但具体数值有所差异。纹理方向对油膜厚度的影响实验结果表明，纹理方向与滑动方向夹角为0°时，油膜厚度相对较大；随着夹角的增大，油膜厚度逐渐减小。这是因为在夹角为0°时，润滑油能够充分利用纹理沟槽形成较厚的油膜；当夹角增大时，润滑油在纹理沟槽中的流动受到阻碍，难以形成稳定且较厚的油膜。在不同粗糙度和载荷、滑动速度条件下，油膜厚度的变化趋势一致，但受粗糙度和工况条件的影响，具体数值有所不同。纹理方向对温度的影响实验数据显示，随着纹理方向与滑动方向夹角的增大，温度呈现先降低后升高的趋势。在夹角为45°左右时，温度相对较低，因为此时纹理对润滑油的流动和散热起到较好的促进作用；当夹角为0°或90°时，温度相对较高，因为在0°时，润滑油的散热效果不佳，而在90°时，纹理对润滑油的阻碍导致摩擦生热增加，都使得温度升高。不同粗糙度和载荷、滑动速度下，温度的变化趋势相似，但受工况条件的影响，温度升高的幅度有所不同。4.2数值模拟研究4.2.1数学模型建立基于热弹流润滑理论，构建考虑表面形貌影响的数学模型。在雷诺方程中，表面形貌主要通过影响油膜厚度和润滑油的流动特性来对热弹流润滑产生作用。对于具有微观粗糙度的表面，油膜厚度可表示为h(x,y,t)=h_0(x,y,t)+S(x,y,t)，其中h_0(x,y,t)为光滑表面的油膜厚度，S(x,y,t)为表面形貌函数，它描述了表面微观几何形状的起伏情况。将该油膜厚度表达式代入雷诺方程，得到考虑表面形貌影响的雷诺方程：\frac{\partial}{\partialx}(\frac{\rho(h_0+S)^3}{12\eta}\frac{\partialp}{\partialx})+\frac{\partial}{\partialy}(\frac{\rho(h_0+S)^3}{12\eta}\frac{\partialp}{\partialy})=U\frac{\partial(\rho(h_0+S))}{\partialx}在这个方程中，表面形貌函数S(x,y,t)的引入使得方程更加复杂，因为它增加了方程的非线性程度。表面形貌的不规则性会导致润滑油在流动过程中受到额外的阻力和扰动，从而改变油膜压力的分布。当表面存在粗糙度时，润滑油在流经微凸体和微凹谷时，流速和压力会发生变化，这种变化通过雷诺方程中的各项参数体现出来。能量方程在考虑表面形貌影响时，主要关注表面形貌对摩擦生热和热传递的影响。由于表面形貌的存在，接触表面的实际接触面积和接触压力分布发生改变，进而影响摩擦生热的大小和分布。粗