粗糙峰微接触及其对润滑的影响

根据椭圆弹塑性接触模型,提出了流量因子影响模型,给出了几个该模型的应用计算实例．计算结果表明,该模型与粗糙表面的实际润滑状况有较好的一致性．通过算例分析,发现粗糙峰弹塑性变形对流量因子的影响主要决定于表面综合属性、弹性和塑性变形量．     

单峰接触研究及其在分形表面接触中的应用

基于有限元方法,建立了弹塑性单峰的接触模型.粗糙峰为理想的弹塑性材料,为了考虑不同的材料特性对微凸体变形的影响,分别对9种不同的材料进行了分析.根据有限元计算结果,分析了接触面积,平均接触压力和接触力与变形干涉量之间的关系,并进行了经验公式的拟合.单峰接触所经历的4个不同的阶段,以及不同阶段之间的转化点均作了明确的表达.然后,根据分形理论,将单峰接触模型扩展到了三维的粗糙表面的接触,并提出了一个计算接触表面法向刚度的模型.通过与实验数据和以往模型的结果对比,证明本文中所提出的模型具有较高的精度.     

基于微凸体接触压力分布的粗糙表面法向接触建模

提出一种考虑微凸体弹塑性接触变形影响的粗糙表面法向接触力学模型。采用有限元模拟微凸体弹塑性接触过程,分析不同塑性屈服条件对微凸体接触载荷和实际接触面积的影响。再根据微凸体接触面上压力分布的变化规律,将微凸体的接触状态分为完全弹性接触阶段、弹塑性接触阶段、完全塑性接触阶段。分析接触面压力变化规律对微凸体法向接触载荷-变形的影响,再利用GW模型的数理统计分析的方法得到粗糙表面的法向接触载荷。将本文提出的模型与完全弹性模型、CEB模型、ZMC模型、KE模型、JG模型进行对比,并且研究了塑性指数对粗糙面接触载荷-平均高度距离的影响。结果表明,本文提出的模型能够更好地描述微凸体法向接触载荷与接触变形的变化趋势,模型预测粗糙表面法向载荷与ZMC、KE模型具有较好的一致性;粗糙面接触载荷随着平均接触距离增加而减少,随着塑性指数的增加,不同模型预测的法向接触载荷差异逐渐增大。     

分形粗糙表面弹塑性接触力学模型

为了建立更为精准的粗糙表面接触模型,基于分形理论提出一种修正的弹塑性接触模型,该模型同时考虑了微凸体的弹性、弹塑性、完全塑性三种变形状态;建立了临界第一弹塑性接触面积、临界第二弹塑性接触面积与临界弹性接触面积之间的关系,分别推导出三种变形阶段的接触载荷与接触面积之间的关系式;结合海洋岛屿面积分布密度函数,获得了接触表面上总接触载荷与真实接触面积之间的表达式。计算结果表明:单个微凸体的临界弹性接触面积不但与材料属性和分形参数有关,也与微凸体的尺寸有关;分形参数一定时,随着接触载荷的增大,真实接触面积也在增大,弹性接触面积与真实接触面积的比值逐渐减小。该模型的建立为进一步研究粗糙表面的摩擦、磨损、润滑提供了理论依据。     

表面粗糙度特征对齿轮接触区润滑特性的影响

大部分工程实际粗糙表面符合非高斯分布,并对齿轮接触副润滑特性有重要影响.将渐开线齿轮啮合过程中齿面接触等效为三维无限长线接触,建立了一个可分析直齿轮和斜齿轮的混合弹流润滑计算模型;采用基于快速傅里叶变换的数值仿真方法生成给定参数的非高斯粗糙表面;运用该模型对直齿轮和斜齿轮啮合过程进行分析,求得不同表面粗糙度特征齿轮在各个啮合点的油膜厚度、接触区载荷以及接触区比例的情况.结果表明:对于标准差相等的非高斯粗糙表面,偏度值对齿轮润滑状况的影响与工况紧密相关,在润滑良好的条件下,偏度值越小润滑状况越优;润滑恶劣的条件下,偏度值越大润滑状况越优;而在各种工况下,峰度值对齿轮润滑状况的影响都表现出峰度值越大润滑状况越优的特点.     

椭圆接触弹流润滑气穴现象的实验观察

采用多光束干涉测量技术对椭圆滚子-玻璃盘形成的椭圆接触气穴现象进行了实验观察.实验在传统光弹流实验机上进行,只研究了纯滑动的情形.结果表明:气穴形状及其发展过程依赖于椭圆接触区短轴与卷吸方向之间的夹角.在一定的载荷-速度条件下,出现片状气穴区;载荷一定时,椭圆接触区短轴与卷吸方向之间的夹角越大,形成片状气穴区需要的速度越高.卷吸速度增大,气穴区向接触中心靠近;载荷增大,气穴区远离接触中心.椭圆接触区短轴与卷吸方向之间的夹角越大,气穴区越是远离接触中心,且气穴区位置受速度影响越显著.     

直流磁场下销盘摩擦接触区的电磁感应现象

为了研究直流磁场对45钢销盘摩擦副的摩擦磨损特性的影响,根据磁场销盘摩擦试验机的结构和销盘摩擦副在摩擦过程中的实际接触情况,建立了二维微凸峰接触静态磁场和瞬态磁场有限元模型,分析了销盘摩擦接触区的电磁感应现象,得出以下结论:磁感应强度B在摩擦接触区分布不均,在微凸峰接触点区域的磁感应强度B值较大;摩擦试验中,在销盘磨痕和磨屑的微凸峰接触区将产生较高频率的动态磁化,同时在微凸峰上产生较大的感应电流,这些现象促进了销盘磨痕表面和磨屑的氧化.     

粗糙表面的弹塑性接触研究

建立了综合载荷作用下粗糙表面弹塑性接触的确定性模型,考虑了微凸峰接触的弹塑性变形阶段,数值求解得到实际接触面积、压力分布和微凸峰塑性形变.分析了实际接触面积与法向载荷的关系,并研究了点接触的椭圆参数对上述关系的影响.建立了结点增长模型,分析了结点增长与滑动摩擦系数的关系以及滑动摩擦系数随椭圆参数的变化.结果表明:随着法向载荷增大,实际接触面积与法向载荷之间的非线性关系愈加显著;椭圆参数越大,实际接触面积越小,选择较小的椭圆参数可降低平均接触压力;结点增长的速率随滑动摩擦系数增大而增大;表面剪切作用力使最大Mises应力值升高,疲劳裂纹的发生源向表面靠近;重载时选择较小的滚动轴承沟曲率半径系数有利于减小摩擦功耗.     

点接触弹塑性流体动力润滑研究

建立了点接触混合润滑模型,根据下表面应力分布迭代求解出下表层的塑性应变,将下表面塑性应变等效转化为本征应变,结合半无限体内本征应变对弹性场的应力扰动解法求解残余应力,表面塑性变形根据本征应变采用半解析方法求解.计算结果表明:本混合润滑模型在塑性计算模块、弹塑性流体动力润滑计算均表现出了很好的准确性以及高效性;本模型能够模拟真实机加工粗糙表面下弹塑性混合润滑问题;能够模拟由全膜润滑、混合润滑、边界润滑以及干接触全工况下的润滑情况,当滚动速度逐渐减小时,平均油膜厚度逐渐减小,接触区由全膜润滑转变为混合润滑,最终演变干接触.     

考虑摩擦系数和微凸体相互作用的粗糙表面接触热导分形模型

基于三维分形理论，建立了考虑摩擦系数和微凸体相互作用的粗糙表面接触热导分形模型，并且考虑了微凸体的弹性变形、弹塑性变形和完全塑性变形. 通过该模型，分析了摩擦系数、分形维数、分形粗糙度和接触载荷对热接触热导的影响. 研究结果表明:接触热导随着摩擦系数和分形粗糙度的增大而减小，随着分形维数和接触载荷的增大而增大. 该研究为开展接合面的热传递提供了一定的理论基础.      

微凸体碰撞对接触应力应变的影响

基于有限元法建立了单对半球状微凸体的接触、碰撞模型,研究了碰撞对微凸体在完全弹性、弹塑性和完全塑性三种接触状态下的Mises应力分布和等效塑性应变(PEEQ)的影响.结果表明:接触变形状态不同时,微凸体的Mises应力分布在碰撞过程具有不同的变化过程,但碰撞结束后,Mises应力值都降低,仅在完全塑性接触时,个别节点的Mises应力峰值大于屈服极限.碰撞使微凸体的等效塑性应变迅速增加到峰值,并保持到碰撞结束,与接触变形状态无关.完全弹性、弹塑性接触状态和完全塑性接触状态下,微凸体等效塑性应变的峰值分别位于接触中心区域和应变峰值环,这意味着裂纹可能产生的部位也不同.     

Hydrodynamic lubrication in fully plastic asperity contacts

A line contact inlet zone analysis is carried out for the hydrodynamic lubrication in a fully plastic asperity contact. A governing equation of the central film thickness i.e. the film thickness in the fully plastic contact area is derived. An equation predicting this film thickness is also derived. It is found that for the fully plastic contact, under relatively light loads the prediction accuracy for the central film thickness is good, while at the load heavy enough the prediction equation greatly overestimates the central film thickness and the central film thickness solved from the analytical governing equation is significantly low showing the asperity in boundary layer lubrication. For the fully plastic contact, the central film thickness is nearly half of that obtained based on the elastic contact assumption for relatively light loads or even lower for heavier loads. The hydrodynamic lubrication is found difficult to form in the fully plastic asperity contact for the carried load heavy enough or the significantly low sliding speed between the asperities. To achieve a high hydrodynamic lubrication film thickness in the fully plastic asperity contact it is recommended to employ a high sliding speed or a high fluid viscosity. However, in the fully plastic asperity contact, the potential hydrodynamic load-carrying capacity is limited and much smaller than that based on the elastic contact assumption or predicted by conventional line contact elasto-hydrodynamic lubrication theory.     

Elastic-plastic contact model for rough surfaces based on plastic asperity concept

A mathematical formulation for the contact of rough surfaces is presented. The derivation of the contact model is facilitated through the definition of plastic asperities that are assumed to be embedded at a critical depth within the actual surface asperities. The surface asperities are assumed to deform elastically whereas the plastic asperities experience only plastic deformation. The deformation of plastic asperities is made to obey the law of conservation of volume. It is believed that the proposed model is advantageous since (a) it provides a more accurate account of elastic-plastic behavior of surfaces in contact and (b) it is applicable to model formulations that involve asperity shoulder-to-shoulder contact. Comparison of numerical results for estimating true contact area and contact force using the proposed model and the earlier methods suggest that the proposed approach provides a more realistic prediction of elastic-plastic contact behavior.     

A Molecular Dynamics Model for Single Adhesive Contact

The normal adhesive contact between a pair of asperities is performed using molecular dynamics. To simplify the problem, the equivalent contact problem of sphere–plane interaction is solved. Displacements in the contact zone are very small compared to the asperity size, therefore, the computational model is focused on the neighborhood of the contact area. The adhesion between the asperity and the plane is calculated as a sum of interactions between atoms of the asperity and the plane. A computational experiment of pull-on and pull-off is carried out to study the influence of the adhesion on the formation of the contact forces and deformations. The numerical results are compared with theoretical predictions.     

考虑粗糙结合面弹塑性接触的黏滑摩擦建模

提出一种同时考虑粗糙面上微凸体弹性变形和塑性接触的切向黏滑摩擦建模方法。采用Hertz弹性理论和Mindlin解描述弹性接触微凸体的切向载荷和相对变形的关系;采用AF(Abbott-Firstone)塑性理论和Fujimoto模型描述塑性接触微凸体切向载荷和相对变形的关系。再利用GW(Greenwood-Williamson)模型统计分析方法建立粗糙表面切向载荷和相对变形之间的关系。将模型与仅考虑微凸体弹性接触情况的模型进行对比,并研究了不同塑性指数对切向载荷和相对变形关系的影响。结果表明：与完全弹性接触模型相比,本文模型引入了塑性接触理论,能够更好地描述粗糙表面切向载荷和相对变形关系,并且考虑不同接触条件下弹性变形微凸体和塑性变形微凸体对切向接触载荷的贡献,在微滑移阶段,主要由弹性接触变形影响,而在进入宏观滑移阶段之后,切向行为主要由塑性变形影响。界面切向载荷由黏着和滑移接触作用共同决定,随着切向变形的增加,滑移接触力逐渐增加,而黏着接触力先增加后减少,反映了界面由微滑移逐渐向宏滑移演化的过程。随着塑性指数的增加,粗糙面上发生塑性接触的微凸体数目逐渐增加,切向黏滑行为主要受到塑性接触特征的控制。     

Elasto-plastic contact of rough surfaces: a mixed-lubrication model for the textured surface analysis

An improved asperity contact model for two rough surfaces with misalignment is presented in this study. The contact model is statistical and can account for the elasto-plastic behavior of interacting asperities. By combining the improved asperity contact model and the average flow Reynolds equation together, a mixed-lubrication model is developed to understand the effect of surface texturing. By comparing with the results of the purely elastic asperity contact model, it is found that the improved asperity contact model can predict the contact force and actual contact area more accurately, particularly under high load conditions. Moreover, comparing with the elasto-plastic model with an equivalent rough surface against a plane, the improved contact model can consider the effect of permitting misalignment of two rough surfaces. This is beneficial for analyzing the performance of the textured piston ring/liner system, especially when asperities contact and wear happen.     

固-液结合面法向动态接触刚度及阻尼的理论模型与实验分析

固-液接触状态广泛存在于机床核心单元关键零部件的接触运动副中,精确获得固-液结合面法向接触刚度及阻尼参数是高档数控机床产品在研发阶段就存在的一个关键理论与技术问题,并且仍然尚未根本解决。固-液结合面在介观层面上表现为两个粗糙表面的接触,在微观层面上表现为微凸体之间的接触,并在中/重载荷作用下微凸体可能会发生弹性/弹塑性/塑性变形。为了揭示静动态外载荷对固-液结合面接触刚度及阻尼的影响,分别基于GW模型、KKE模型和AF模型对接触微凸体弹性/弹塑性/塑性变形展开研究,并结合流体动力润滑REYNOLDS方程,建立了考虑接触微凸体弹性/弹塑性/塑性变形的固-液结合面接触刚度及阻尼模型。并对其进行实验验证,结果表明：随着预载荷的增大固-液结合面法向动态接触刚度表现出先减小后增大的规律,当接触载荷小于某阈值时动态接触刚度较大,反之静态接触刚度较大;法向动态接触刚度随着法向相对位移幅值的增大而增大,在低载荷时呈线性规律,而高载荷呈非线性规律;法向动态接触刚度随激振频率增大呈线性增大,且载荷越大线性斜率越小。对于法向接触阻尼,随着法向相对位移幅值和接触载荷增大呈非线性增大,随着激振频率增大几乎不变。精确获得固-液结合面法向接触刚度和阻尼及其关键因素的影响规律,对机械系统的分析、设计、优化以及静、动态性能控制都具有重要的理论意义。     

Adhesion-induced instabilities in elastic and elastic-plastic contacts during single and repetitive normal loading

Adhesive interaction in spherical contacts was modeled with the Lennard-Jones (L-J) potential. Elastic adhesive contact was analyzed by the equivalent system of a rigid sphere with reduced radius of curvature and a half-space of effective elastic modulus. The critical gap at the instant of abrupt surface contact (jump-in) and separation (jump-out) was determined from the deformed surface profile of the elastic half-space and geometrical relationships. A finite element model of a rigid sphere and an elastic-plastic half-space was used to examine elastic-plastic adhesive contact. Surface adhesion was modeled by nonlinear springs with a force-displacement relationship governed by the L-J potential. The evolution of the interfacial force and the central gap distance as well as the occurrence of jump-in and jump-out instabilities were investigated in terms of the Tabor parameter, plasticity parameter, and dimensionless maximum normal displacement. The force-displacement response due to several approach-retraction cycles was interpreted in the context of elastic and plastic shakedown behaviors using dimensionless parameters.     

Strain gradient plasticity analysis of elasto-plastic contact between rough surfaces

From a microscopic point of view, the real contact area between two rough surfaces is the sum of the areas of contact between facing asperities. Since the real contact area is a fraction of the nominal contact area, the real contact pressure is much higher than the nominal contact pressure, which results in plastic deformation of asperities. As plasticity is size dependent at size scales below tens of micrometers, with the general trend of smaller being harder, macroscopic plasticity is not suitable to describe plastic deformation of small asperities and thus fails to capture the real contact area and pressure accurately. Here we adopt conventional mechanism-based strain gradient plasticity (CMSGP) to analyze the contact between a rigid platen and an elasto-plastic solid with a rough surface. Flattening of a single sinusoidal asperity is analyzed first to highlight the difference between CMSGP and J₂ isotropic plasticity. For the rough surface contact, besides CMSGP, pure elastic and J₂ isotropic plasticity analysis is also carried out for comparison. In all cases, the contact area A rises linearly with the applied load, but with a different slope which implies that the mean contact pressures are different. CMSGP produces qualitative changes in the distributions of local contact pressures compared with pure elastic and J₂ isotropic plasticity analysis, furthermore, bounded by the two.     

Finite element analysis of large contact deformation of an elastic–plastic sinusoidal asperity and a rigid flat

Normal contact deformation of an asperity and a rigid flat is studied within an axisymmetric finite element model. The asperity features a sinusoidal profile and is elastic–plastic with linear strain hardening. Influences of geometrical (asperity height and width) and loading (the maximum interference) parameters on frictionless contact responses are explored for both loading and unloading. Dimensionless expressions for contact size and pressures covering a large range of interference and asperity ratio values are obtained in power-law forms. Results show the mean contact pressure after fully-plastic contact reaches a plateau only for small asperity ratios, while it continues increasing for large asperity ratios. The residual depth is found to be associated with plastically dissipated energy.