Elastoplastic contact mechanics model of rough surface based on fractal theory

Because the result of the MB fractal model contradicts with the classical contact mechanics, a revised elastoplastic contact model of a single asperity is developed based on fractal theory. The critical areas of a single asperity are scale dependent, with an increase in the contact load and contact area, a transition from elastic, elastoplastic to full plastic deformation takes place in this order. In considering the size distribution function, analytic expression between the total contact load and the real contact area on the contact surface is obtained. The elastic, elastoplastic and full plastic contact load are obtained by the critical elastic contact area of the biggest asperity and maximun contact area of a single asperity. The results show that a rough surface is firstly in elastic deformation. As the load increases, elastoplastic or full plastic deformation takes place. For constant characteristic length scale G, the slope of load-area relation is proportional to fractal dimension D. For constant fractal dimension D, the slope of load-area relation is inversely proportional to G. For constant D and G, the slope of load-area relation is inversely proportional to property of the material ϕ, namely with the same load, the material of rough surface is softer, and the total contact area is larger. The contact mechanics model provides a foundation for study of the friction, wear and seal performance of rough surfaces.     

粗糙表面接触模型的研究进展

实际工程中的接触表面都不是绝对光滑的,当两物体相互接触时,真实表面的接触实际上是微凸体间的接触.粗糙表面间的接触行为对摩擦、磨损、润滑、密封和传热等有着重要的影响,是摩擦学研究的主要课题之一.介绍了Hertz弹性接触模型、统计学接触模型和分形接触模型,对粗糙表面接触模型的研究现状和进展作了分析和评述.由Hertz模型只能求解一些几何形状比较规则的物体,应用范围非常有限;各种统计学模型都是经过简化的理想模型,这些简化究竟能导致多大误差目前还没有精确的分析;分形模型利用了包含全部表面粗糙度信息的分形参数D和G,其对粗糙表面接触性质的预测不受仪器分辨率和取样长度的影响,使预测具有惟一性或确定性.     

新的粗糙表面弹塑性接触模型

提出一种新型的粗糙表面弹塑性微观接触模型.该模型的建立基于接触力学理论和接触微凸体由弹性变形向弹塑性变形及最终向完全塑性变形的转化皆是连续和光滑的假设.研究单个微凸体在载荷逐渐增加时的变形规律,并重点推出弹塑性变形区间的接触方程.在此基础上应用概率统计理论导出了粗糙表面的接触载荷、平均分离和实际接触面积之间的数学关系式.在不同的塑性指数和载荷条件下,该模型与GW弹性模型和CEB弹塑性模型就实际接触面积和法向距离的预测结果进行了对比.结果表明,在同样塑性指数和载荷条件下比GW模型预测的实际接触面积大但法向距离小,且两者的差距随塑性指数和载荷的增加而增大.因此该模型的预测结果更加符合人们的试验观察和直觉,能够更加科学和合理地描述两个粗糙表面的微观和宏观接触状态.     

Resolving the contradiction of asperities plastic to elastic mode transition in current contact models of fractal rough surfaces

The paper suggests a revision to the asperities plastic to elastic mode transition in the elastic–plastic contact model of fractal rough surfaces, offered by Majumdar and Bushan [A. Majumdar, B. Bushan (MB model) J. Tribol. 113 (1991) 1–11.]. According to the MB model, the contact mode of a single fractal asperity transfers from plastic to elastic when the load increases and the growing contact area exceeds a certain critical value, which is scale independent. This surprising result of the MB model is in contrast with classical contact mechanics where increasing contact area due to increased load is associated with a transition from elastic to plastic contact. The present study describes a revised elastic–plastic contact model of a single fractal asperity showing that, contrary to the MB model prediction, the critical contact area is scale dependent. The revised model also shows that a fractal asperity behaves as would be expected from classical contact mechanics namely, as the load and contact area increase a transition from elastic to plastic contact mode takes place in this order.     

依据各向异性分形几何理论的固定结合部法向接触力学模型

基于各向异性分形几何理论,考虑微凸体变形特点、表面微凸体承受法向载荷的连续性和光滑性原理,以及区分微凸体分别处于弹性、塑性变形时的一个微凸体实际微接触面积,建立固定结合部法向接触力学模型。采用二变量Weierstrass-Mandelbrot函数模拟各向异性三维分形轮廓表面。推导出划分弹塑性区域的临界弹性变形微接触截面积、结合部量纲一法向载荷、结合部量纲一法向接触刚度的数学表达式。数值仿真结果表明：当表面形貌的分形维数、分形粗糙度一定时,真实接触面积随着结合部法向载荷的增大而增大;结合部法向接触刚度随着真实接触面积、结合部法向载荷、相关因子或材料特性参数的增大而变大;当分形维数由1变大时,结合部法向接触刚度随着分形维数的变大而增大;当分形维数增加到趋近于2时,结合部法向接触刚度有时却会随着分形维数的增加而降低。结合部法向接触力学模型的构建,有助于分析固定接触表面间的真实接触情况。     

基于分形理论的滑动摩擦表面接触力学模型

依据分形理论,考虑微凸体变形特征及摩擦作用的影响建立滑动摩擦表面接触力学模型。采用一个三次多项式来表达弹塑性变形微凸体的接触压力与接触面积的关系,从而满足在变形状态转变临界点处的微凸体接触面积与接触压力转化皆是连续和光滑的条件。推导出滑动摩擦表面临界弹性变形微接触面积、临界塑性变形微接触面积、量纲一真实接触面积的数学表达式。理论计算结果表明,表面形貌一定时,真实接触面积随着载荷的增大而增大;载荷一定时,真实接触面积随着特征尺度系数的增大而减小,随着分形维数的增大先增大后减小;当表面较粗糙时,摩擦因数对真实接触面积的影响很小;随着表面光滑程度的增大,摩擦因数对真实接触面积的影响增大,真实接触面积随着摩擦因数的增大而增大,特别是当摩擦因数较大时,真实接触面积增大的幅度也较大。接触力学模型的建立,为研究滑动摩擦表面间的摩擦磨损性能提供了依据。     

Transient Heat Conduction Between Rough Sliding Surfaces

When two rough bodies slide against each other, asperities on the opposing surfaces interact with each other, defining a transient contact and heat conduction problem. We represent each body by a Greenwood and Williamson asperity model with a Gaussian height distribution of identical spherical asperities. The heat transfer during a typical asperity interaction is analyzed, and the results are combined with the height distributions to determine the mean heat flux and the mean normal contact pressure as functions of the separation between reference planes in the two surfaces. We find that the effective thermal conductance is an approximately linear function of nominal contact pressure, but it also increases with the square root of the sliding speed and decreases with the 3/4 power of the combined RMS roughness. The results can be used to define an effective thermal contact resistance and division of frictional heat in macroscale (e.g., finite element) models of engineering components, requiring as input only the measured roughness and material properties.     

基于域扩展因子和微凸体相互作用的结合面接触刚度模型研究

结合面接触刚度直接影响了机械设备的整机动态特性,为了建立更为准确的接触刚度模型,以分形几何理论为基础,利用单一微凸体承受局部载荷时的弹性变形特性,并基于域扩展因子引入微接触截面积分布函数,推导了考虑表面微凸体相互作用影响的结合面接触刚度分形模型。为了验证所提出模型的准确性,通过三维非接触式测量,获得了试验试样的表面轮廓数据,并根据结构函数法,计算了各个试样的表面分形参数,进而将理论接触刚度与试验结果对比分析,结果表明：法向接触刚度的增长速率与粗糙面表面临界接触面积有关,临界接触面积决定了结合面内的弹性变形占比。考虑微凸体相互作用后,所提出模型的预测曲线更加符合试验中法向载荷与接触刚度的关系。     

大偏心下滑动轴承热弹流润滑状态及特性*

滑动轴承在大偏心条件下工作时，热效应及弹性变形使得油膜润滑状态发生变化，进而影响摩擦特性。为此建立耦合轴瓦弹性变形、轴颈轴瓦粗糙峰接触、油膜温度分布及黏温-黏压关系的滑动轴承混合润滑模型，采用有限差分法求解得到不同工况下油膜压力场、温度场分布，分析热效应及弹性变形对润滑状态转变及轴承各特性参数的影响；搭建实验台测量试件内表面温度分布，测试结果验证了计算模型的正确性。结果表明：大偏心时热效应和弹性变形使得油膜润滑状态出现转化；粗糙峰的接触使摩擦热增加，且在最小油膜处形成温度峰值；热效应和轴瓦弹性变形使得接触压力峰值集中在轴承两端，承载能力和摩擦力均有所下降。     

Effect of asperity interactions on rough surface elastic contact behavior: Hard film on soft substrate

An improved elastic micro-contact model of rough surfaces accounting for asperity interactions is proposed. The contact behavior of a single asperity system is composed of a stiffer hemi-spherical asperity deformation and bellowing softer substrate deformation, which is then extended to rough surface contact including asperity interactions. Using the solution of substrate deformation, normal positions of individual asperities are adjusted during quasi-static contact, from which surface interactive forces are obtained. Analytical simulations are performed using the proposed rough surface contact model, whose results are compared to Greenwood-Williamson-based models and with experimental measurements.     

热力耦合下微接触点塑性应变沿深度变化分析

在考虑粗糙实体弹塑性变形、热力耦合、微凸体间相互作用和摩擦热流耦合等影响下,运用有限元法数值模拟具有三维分形特性的粗糙面与刚性平面间滑动摩擦过程,分析了粗糙实体接触凸点塑性变形随深度变化情况。发现：在速度的突变和闪点温度形成时,摩擦接触表层等效塑性应变增大明显;在这一摩擦表层,过不同接触点的纵向剖面塑性应变沿深度分布不同：有的是接触表面塑性变形最大,有的是在接触微凸体表面下某一深度塑性变形最严重,而接触凸点表面的塑性应变稍小些。这与相关文献用SEM研究干摩擦后金属摩擦表层变形照片后发现的结果一致。滑动摩擦过程中,金属粗糙摩擦接触表层塑性变形的不断累积,将会导致材料表层中的夹杂或微观缺陷周围萌生微孔和裂纹源。     

粗糙表面微凸体对液黏传动的影响

为研究液黏传动过程中粗糙表面的承载特性,将分形理论引入到两粗糙表面摩擦过程之中,分析传动过程中混合摩擦和边界摩擦两阶段的微凸体承载过程,考虑微凸体弹塑性变形,对M-B模型进行修正,建立修正的微凸体承载模型。建立基于修正M-B模型的微凸体承载模型。通过数值仿真得到有效面积系数、分形参数对液黏调速离合器传动过程的影响规律;对修正的微凸体承载模型的计算结果与M-B模型的计算结果进行对比分析。结果表明:微凸体接触载荷和传递转矩随着面积比的增大而增大,当有效面积系数与尺度系数增大时,接触载荷与传递转矩均有所增大;分形维数为1.5时,微凸体接触载荷与传递转矩最小且随面积比的变化最为缓慢;在整个接触区域内,弹性变形区域的面积、接触载荷以及传递转矩最大,其次是弹塑性变形区域,塑性变形区域最小;考虑弹塑性变形时,微凸体接触载荷与传递转矩均有所下降;修正M-B模型和M-B模型间的修正系数范围在25%以内,修正系数随着有效面积系数、尺度系数的增大而增大,随着分形维数的增大而减小。     

尺度相关的分形结合面法向接触刚度模型

基于分形理论,利用双变量Weierstrass-Mandelbrot函数模拟三维分形结合面,建立尺度相关的三维分形结合面法向接触刚度模型。推导出各等级微凸体发生弹性、弹塑性以及完全塑性变形的存在条件。确定结合面上各等级微凸体的面积分布密度函数,推导出法向接触刚度和法向接触载荷的解析表达式。计算结果表明：当结合面上的微凸体只能发生弹性变形,即自身等级小于弹性临界等级的微凸体,该部分微凸体引起的法向接触刚度和对应法向载荷关系呈非线性。当微凸体的等级大于弹性临界等级,在结合面接触过程中,微凸体弹性变形引起的法向接触刚度与对应的法向载荷关系为线性,非弹性变形引起的法向接触刚度与法向载荷关系为非线性。微凸体的等级范围对结合面的刚度影响较大,在相同的法向载荷作用下,高等级微凸体的结合面产生较高的法向接触刚度,即结合面越平整,结合面的法向刚度越高。     

On the Relation of Load to Average Gap in the Contact Between Surfaces with Longitudinal Roughness

A computer simulation model for the contact between longitudinally-oriented rough surfaces has been formulated. This model closely duplicates the actual surf ace contact deformation behavior by taking into account the elastic interactions between the asperities. There were no assumptions made about the shapes, or any deformation behavior of the asperities, except for their obeying the laws of elasticity. The plastic deformations on the high asperity peaks were taken into account by setting a ceiling on their contact pressures at the material hardness value. The simulations used real surface profiles which were digitized from unworn circumferentially ground steel surfaces. Each pair of these profiles was mathematically combined to form an equivalent rough profile pressing against an infinitely rigid flat and having the appropriately adjusted elastic modulus. A total of 28 different pairs of profiles were used in the simulations. Each contacting pair was subjected to 30 different load levels and the local contact pressures and deformations were calculated. The contact simulations yielded some important mathematical relationships between parameters, such as the real area of contact, average gap, and average asperity load through statistical curve fitting. Two analytical functions were generated to relate the average load to average gap and the real area of contact to load.     

基于真应力-应变关系的粗糙表面法向接触模型

提出一种粗糙表面的法向弹塑性接触分析的建模方法。基于微凸体的弹塑性有限元接触模型,分别研究了40Cr、45和Q235三种钢材料的微凸体与刚性平面的法向接触特性。有限元模型中采用三种材料的真应力-应变关系,考察了不同强化特性对微凸体接触性质的影响。建立了微凸体在弹性、弹塑性、塑性变形阶段统一的接触变量变化规律的表达式。在此基础上应用概率统计理论建立粗糙表面法向弹塑性接触模型。所建立的接触模型中微凸体接触变量的变化规律完全基于弹塑性有限元模型的计算结果,无需将微凸体的变形过程区分为不同的变形阶段,避免了接触变量在各阶段采用不同函数表达式带来的连续性和光滑性问题,以及在弹塑性阶段采用插值函数的随意性问题。通过与其他接触模型的计算结果相比较,证明了所提出接触模型的合理性。     

The Effect of Scale-Dependent Hardness on Elasto-Plastic Asperity Contact between Rough Surfaces

Statistical methods are used to model elasto-plastic contact between two rough surfaces using a recent finite element model of elasto-plastic hemispherical contact and also recent advances in strain gradient modeling. The elasto-plastic hemispherical contact model used to model individual asperities accounts for a varying hardness effect due to deformation of the contact geometry that has been documented by other works. The strain gradient model accounts for changes in hardness due to scaling effects. The contact between surfaces with hypothetical material and surface properties, such as the elastic modulus, yield strength, and roughness are modeled. A model is also constructed to consider a variable asperity contact radius to evaluate if the strain gradient model will affect it differently. The models produce predictions for contact area, contact force, and surface separation. The strain gradient effects decrease the real area of contact and increase the average contact load in comparison to the model without these effects. The strain gradient model seems to have a larger influence on the predictions of contact load and area than does considering a variable asperity contact radius for the cases considered in this work.     

Elastic–plastic adhesive contact of rough surfaces based on plastic asperity concept

The paper describes an analysis of adhesive contact between rough surfaces with small-scale surface asperities using an elastic–plastic model of contact deformation based on fictitious plastic asperity concept developed by Abdo and Farhang [Int. J. Non-Linear Mech. 40 (2005) 495]. The model considers simultaneous occurrence of elastic and plastic behaviours for an asperity. The well-established elastic adhesion index and plasticity index are used to consider the different contact conditions that arise as a result of varying load and material parameters. The load-separation behaviour for different combinations of these parameters is obtained. Comparison with previous elastic–plastic model that was based on elastic-then-plastic assumption is made showing significant differences.     

考虑硬度变化的结合面法向接触刚度分形模型*

基于分形几何理论,考虑微凸体因应变硬化而造成弹塑性变形阶段硬度随变形量变化而变化,建立结合面第一、第二弹塑性变形阶段单次加载刚度分形模型。推导出在计入硬度变化的情况下,单个微凸体在弹塑性变形阶段法向接触刚度与接触面积之间的关系式,进而得出结合面在弹塑性变形阶段法向接触刚度与接触面积、接触载荷之间量纲为一的关系式,并通过仿真分析得出相关参数对结合面法向接触刚度的影响。仿真结果显示:考虑硬度变化时,结合面量纲一法向接触刚度的值与法向实际接触载荷、实际接触面积之间存在关系;结合面法向接触刚度随着分形维数D的增大而增大;分形维数一定时,结合面法向接触刚度随表面长度尺度参数G值增大而增大。     

接触热导率数值预测研究

假设粗糙表面轮廓高度服从高斯分布,并将两粗糙表面间的接触等效为粗糙表面与光滑刚性表面的接触。由均方根粗糙度和平均凸峰坡度得出接触点的平均半径和单位面积内的接触点个数,分别建立了用圆柱体和圆锥台来表现接触粗糙峰的两种有限元模型,对接触热导率进行了数值预测。通过对比发现,用圆柱体来表现粗糙峰的模型计算出的结果与试验吻合较好,且其误差随界面载荷的增大而减小。