摘要:为研究织构表面动压润滑性能,基于Reynolds方程和织构深度方程建立单一织构不可压缩牛顿体润滑条件下动压润滑模型,分析不同织构参数对织构动压润滑性能的影响;并采用正交试验法研究织构类型、织构横截面形状、织构深度和织构分布角度对织构动压润滑性能影响的主次顺序。数值计算结果表明:圆形、椭圆形和条形3种织构中圆形织构性能最优,且织构类型的影响与织构长宽比有关;织构截面形状和织构分布角度的影响与织构深度有密切关系;4种织构参数对承载力影响的主次顺序为:织构深度>织构截面形状>织构类型>织构分布角度,对摩擦系数影响的主次顺序为:织构深度>织构截面形状>织构分布角度>织构类型,因而织构深度对动压润滑性能的影响最大,而织构深度在最小油膜厚度左右时织构动压润滑性能最好。
关键词:表面织构;动压润滑性能;正交试验
《大连工业大学学报》(双月刊)创刊于1981年,是由大连工业大学主办的以反映大连轻工业大学科研成果为主的自然科学综合性学术期刊。
人類从大自然的生物表面并非绝对光滑,而是存在一定的凹坑或凸起,但却具有良好的抗磨减阻效果,提出仿生表面织构技术,以提高摩擦副表面的润滑及摩擦学性能。不同工况下的实验[1-3]及不同仿真模型下的理论研究[4-6]分析了织构对润滑及摩擦学性能的影响,研究结果表明,基于不同润滑状态下的润滑机理,表面织构的存在对摩擦副的润滑减磨性能有显著的改善。但研究过程中也发现,并非所有参数织构都有润滑减磨效果,织构参数设计不当可能对摩擦性能带来一定的负面影响 [10-11]。因而,在证实表面织构能有效改善摩擦副表面润滑及摩擦学性能的基础上,织构参数优化成为仿生织构研究的重点。
织构参数中,对其动压润滑性能存在影响的主要包括织构类型、织构横截面形状、织构深度、织构分布角度,不同参数对动压润滑性能影响规律并不相同。Tala-Ighil等[12]基于雷诺方程建立织构润滑理论模型,分析了球形、圆柱形和平行六面体织构对润滑性能的影响差异,结果表明织构几何形状对承载能力有重要的影响,且平行六面体织构表面摩擦学性能最优。Nanbu等[13]基于虚拟织构模型和数值仿真方法建立不同截面形状微织构模型,分析织构截面形状对润滑性能的影响,结果指出,相对于平底形状织构,微楔形和微台阶型截面形状织构有利于提升润滑油膜厚度,改善织构表面润滑性能。Han等[14]基于N-S方程建立单一微球形织构三维动压润滑模型,分析织构深度对织构动压润滑性能的影响,数值仿真结果表明,存在最优轴深度使得承载力最大和摩擦系数最小,优选无量纲凹坑深度变化范围为0.80~2.00。Yu等[15]开展了织构分布方向对法向接触平面流体动压形成效应影响的理论研究分析圆形、椭圆形和三角形凹坑织构在不同分布角度下无量纲平均油膜压力,分析结果表明,织构分布方向对接触表面的承载力具有显著影响,侧边垂直于运动方向的三角形织构和长轴垂直于滑动方向的椭圆形织构表面承载力和平均油膜压力均优于圆形凹坑织构。从上述研究可看出,不同织构参数的改变均对织构润滑性能存在一定程度的影响,且任意织构参数均存在最优值以获得最优织构润滑性能。但并没有任何研究分析不同织构参数对织构润滑性能影响的主次顺序,继而分析影响织构润滑性能的主要因素。
因此,本文基于雷诺方程和织构深度方程建立单一织构牛顿体动压润滑条件下的润滑理论模型,通过数值计算分析織构类型、织构横截面形状和织构分布角度对织构动压润滑性能的影响及其规律;并采用正交实验法开展织构类型、织构横截面形状、织构分布角度和织构深度对动压润滑性能影响的正交数值计算,研究4种织构参数对织构动压润滑性能影响的主次顺序,为合理选择优选织构参数及织构参数优选方法提供更多理论指导。
1 模型建立
1.1 几何模型
织构化表面示意图,如图1所示。对于周期性分布的表面织构,为简化计算,本文选择单一织构单元为研究对象进行分析。织构单元长为L,宽为W,织构直径大小为D。织构表面横截面示意图如图1(b)所示, h0为摩擦副间隙,hp为织构深度,U为摩擦副两表面相对滑动速度,本文所选择的织构单元参数如表1所示。
1.2 数学模型
自1886年Reynolds首次阐述流体动压产生机理以来,Reynolds方程被广泛应用于流体动压润滑研究中[16]。基于基本假设:1)忽略流体介质体积力的影响;2)不考虑流体边界滑移的影响;3)油膜厚度方向,油膜压力不变;4)流体为牛顿体,符合牛顿内摩擦定律;5)流体流动状态为层流,没有涡流和紊流;6)忽略流体惯性力的影响,可得到简化后的Reynolds方程用于求解织构表面流体动压,简化后的Reynolds方程可表述为:
式中:x、y为方向坐标;p为油膜动压;h为油膜厚度;u为摩擦副表面相对滑动速度;η为润滑介质黏度。
不同参数表面织构其主要通过影响摩擦副表面油膜厚度而对油膜压力分布造成影响。本文中,不同织构类型、织构横截面形状下润滑油膜厚度方程如表2和表3所示,且表3仅仅列出圆形织构不同截面形状下油膜厚度方程。
此外,为了研究织构不同分布角度下的影响,在编程过程中对织构旋转一定角度后的新坐标,可采用如式(2)所示的坐标变换公式,图2为椭圆织构由粗实线旋转到细虚线后的示意图。
对于单一织构研究对象,由于在周期分布方向油膜压力呈现出周期性分布,因此油膜压力为周期性分布边界条件,另一方向为有限宽边界,则油膜压力假设等于大气压,继而织构单元油膜压力边界条件为:
1.3 求解过程
基于计算时间和计算精度的考虑,对织构单元采用256×256网格节点进行网格化,在结合油膜压力边界方程的基础上采用有限差分法和高斯赛德尔迭代对雷诺方程进行离散求解,忽略织构空化效应的影响,计算得到织构单元表面油膜压力分布。油膜压力收敛判据为:
式中:Pi,j为网格节点(i,j)处油膜压力;k为迭代次数。
由于织构化表面处于流体动压润滑状态,因而表面承载力由流体动压提供,摩擦力则由流体间的剪切应力产生,因此织构表面承载力和摩擦系数计算式可表示为:
1.4 试验方案设计
1.4.1 单因数分析
织构参数对织构润滑性能影响单因素分析主要研究在不同织构深度下,织构类型、织构横截面形状和织构分布角度对织构动压润滑性能的影响,并指出织构类型、织构横截面形状和织构分布角度对动压性能的影响规律,为织构参数优选提供基本优选思路。织构参数如表4所示。
1.4.2 正交试验方案设计
在单因素分析的基础上,通过设计正交单元试验,分析织构类型、织构横截面形状、织构深度和织构分布角度4种织构参数对织构润滑性能的主次影响,找出影响织构润滑性能的主要因素,继而找到通过改变织构参数获取更优织构润滑性能的方法;此外,找出织构参数最优组合,已获得最优织构润滑性能。正交仿真试验各因数及各因素水平如表5所示,正交试验方案则如表6。
2 结果与讨论
2.1 单因素影响分析
织构面积比为20%,织构深度在16 μm范围内,在相同织构深度下圆形、椭圆形和条形3种矩形截面织构表面承载力大小对比,如图3所示。由图3可看出,4种织构深度下,表面承载力的大小关系均表现为:圆形织构表面>条形织构表面>椭圆形织构表面,因此,3种不同织构中,圆形织构性能最优,其次为条形织构,椭圆形织构性能最差。此外,从图3中也可看出,圆形织构表面承载力相比于条形和椭圆形织构表面承载力要大很多,条形及椭圆形织构表面承载力的差距则很小。而圆形织构与条形及椭圆形织构间的较大差距为织构在长宽比上的差距,进而织构类型对承载力的影响与织构长宽比有较大的关系,织构类型优选可考虑不同织构类型长宽比进行。
圆形织构、织构面积为30%、织构深度在16 μm范围内时,三角形、抛物线和矩形织构截面对承载力的影响,如图4所示。从图4可知,织构深度为
1 μm、2 μm和4 μm条件下,不同截面织构表面承载力的大小关系为:矩形截面织构>抛物线截面织构>三角形截面织构,而织构深度为 6 μm时,则三角形截面织构>抛物线截面织构>矩形截面织构,因此,织构截面形状对织构动压润滑性能的影响与织构深度有较大的关系。相同情况下,织构最大深度与织构平均深度的大小关系,如图5所示。由图5可知,3种不同截面织构在4种相同织构深度下,其平均织构深度的大小关系为:矩形截面织构>抛物线截面织构>三角形截面织构。图4为摩擦副表面最小油膜厚度为2 μm条件下表面承载力,图5中4种不同织构深度下,织构平均深度最接近2 μm的分别为矩形截面织构(1 μm)、矩形截面织构(2 μm)、抛物线截面织构 (4 μm),三角形截面织构(6 μm)。因此,对比图4和图5可知,不同截面形状织构对动压润滑性能的影响规律基本呈现为织构平均深度越接近摩擦副表面油膜厚度,织构的动压润滑性能越好。继而证实了织构截面形状对动压润滑性能的影响与织构深度有关的推论。