考虑基体变形的混合润滑结合面接触特性

为研究混合润滑状态下粗糙表面基体变形对结合面接触特性的影响,建立了考虑基体变形的结合面接触刚度模型。首先,通过单微凸体-基体系统模型分别求解微凸体和基体的接触刚度,利用不动点迭代法获得微凸体真实变形量;其次,基于分形理论建立结合面固体接触刚度模型,通过固体接触刚度获得液体介质的接触刚度。根据仿真结果分析了基体变形、粗糙表面形貌以及润滑介质对结合面接触特性的影响规律。结果表明:当真实接触面积一定时,通过新模型计算的法向载荷小于忽略基体变形的模型;在接触前期,结合面的接触刚度主要由液体介质接触刚度主导,随着真实接触面积的增加,液体接触刚度占总刚度的比率越来越小,最后转变为固体的接触刚度主导结合面的接触刚度。该模型为研究混合润滑状态下结合面的接触特性提供了理论基础。     

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

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

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.     

Effect of surface roughness on contact between solid surfaces

The analysis of the mechanism of contact between two solids was carried out considering the distribution of the radii of curvature of asperity peaks. The analytical results show that the mean radius of curvature of asperity peaks has a considerable effect on the nature of the deformation of contact asperities, i.e. whether the contact is plastic or elastic, and more effect on the real area of contact than the variation of the distribution of the radii of curvature.The radii of curvature at the asperity peaks and the real area of contact between two smooth surfaces were measured for comparison with the theoretical results. The results for isotropic surfaces produced by buffing and sandpaper agree with the theory; the real area of contact increases with decreasing surface roughness.     

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

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

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

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

Development of a wear equation for polymer-metal sliding in terms of the fatigue and topography of the sliding surfaces

A wear equation has been derived using the concept of fatigue failure due to asperity interactions in the contact region between sliding bodies. One of the three principal stresses that arise in the contact zone under the effect of a normal as well as a tangential load is of tensile nature. It is this principal stress that has been considered to be responsible for the initiation and propagation of fatigue cracks. It is assumed that the deformation in the contact zone is of elastic nature and that both the contacting surfaces are covered with asperities that have spherical tips. The wear equation involves the asperity height distribution φ(z). The particular distribution for a sliding situation is determined from experimental studies of the topography of sliding surfaces. The wear equation indicates that the wear rate depends upon the fatigue properties of the weaker material, normal load, sliding speed, coefficient of friction, moduli of elasticity of the contacting materials, asperity density, asperity radius of curvature and the distribution and standard deviation of asperity heights. The variation of wear with these parameters as indicated by the wear equation is in agreement with the experimental studies already reported in the literature.     

Effects of Friction on the Contact and Deformation Behavior in Sliding Asperity Contacts

Finite-element analyses are carried out to study the effects of friction on the contact and deformation behavior of sliding asperity contacts. In the analysis, on elastic-perfectly-plastic asperity is brought in contact with a rigid flat at a given normal approach. Two critical values of the normal approach are used to describe the asperity deformation. One is the approach corresponding to the point of initial plastic yielding, and the other at the point of full plastic flow. Additional variables used to characterize the deformation behavior include the shape and size of the plastic zone and the asperity contact size, pressure, and load capacity. Results from the finite-element analysis show that the two values of critical normal approach decrease significantly as the friction in the contact increases, particularly the approach that causes plastic flow of the asperity. The size of the plastically deformed zone is reduced by the friction when the contact becomes fully plastic. The reduction is very considerable with a high friction coefficient, and the plastic deformation is largely confined to a small thin surface layer. For a low friction coefficient, the contact size, pressure and load capacity of the asperity are not very sensitive to the friction coefficient. For a moderate friction coefficient, the contact pressure is reduced and the junction size increased; the load capacity of the asperity is not significantly affected due to the compensating effects of the pressure reduction and the junction growth. For a high friction coefficient, the pressure-junction compensation is not longer sufficient and the asperity load capacity is reduced. The degree of the friction effects on these contact variables depends on the applied force or the normal approach. Although the analyses are conducted using a line-contact model, the authors believe that the effects of friction in sliding asperity contacts of three-dimensional geometry are essentially the same and the same conclusions would have been reached. These results may provide some guidance to the modeling of rough surfaces in boundary lubrication, in which the asperity friction coefficient can be high and vary significantly both in time and from one micro-contact to another.     

接触热导率数值预测研究

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

Microscopic asperity contact and deformation of ultrahigh molecular weight polyethylene bearing surfaces

The effect of the roughness and topography of ultrahigh molecular weight polyethylene (UHMWPE) bearing surfaces on the microscopic contact mechanics with a metallic counterface was investigated in the present study. Both simple sinusoidal roughness forms, with a wide range of amplitudes and wavelengths, and real surface topographies, measured before and after wear testing in a simple pin-on-plate machine, were considered in the theoretical analysis. The finite difference method was used to solve the microscopic contact between the rough UHMWPE bearing surface and a smooth hard counterface. The fast Fourier transform (FFT) was used to cope with the large number of mesh points required to represent the surface topography of the UHMWPE bearing surface. It was found that only isolated asperity contacts occurred under physiological loading, and the real contact area was only a small fraction of the nominal contact area. Consequently, the average contact pressure experienced at the articulating surfaces was significantly higher than the nominal contact pressure. Furthermore, it was shown that the majority of asperities on the worn UHMWPE pin were deformed in the elastic region, and consideration of the plastic deformation only resulted in a negligible increase in the predicted asperity contact area. Microscopic asperity contact and deformation mechanisms may play an important role in the understanding of the wear mechanisms of UHMWPE bearing surfaces.     

Roughness and Lubricant Effect on 3D Atomic Asperity Contact

Large-scale molecular dynamics simulations were performed to study the sliding process of rough surfaces with and without lubricant. In the dry contact, a linear relationship has been observed between the load and the contact area for surfaces with large root mean square (RMS) roughness. However, it becomes nonlinear when the RMS is small. In the presence of adhesion, small roughness results in a large friction force when the surfaces are flattened and the contact area reaches 60 %. In order to confirm this observation, nonadhesive models have been established with an observation that under the combined influence from roughness and adhesion, the contact area plays a crucial role to determine whether the dry sliding is under the domination of roughness or adhesion. In the lubricated sliding, an increase in friction force has been found for the partially lubricated condition because the asperity contact still accounts for a great deal of resisting force. Besides, the lubricant exerts a comparable resisting force to the sliding.     

On the resolution dependence of micromechanical contact models

In this study we deal with the problem of normal contact between rough surfaces. A power spectral density of the asperity heights which obeys a power-law within a given range of wavelengths has been considered. Hence, the resolution dependence of the statistical parameters computed according to the random process theory, of the plasticity index and of the contact predictions provided by stochastic and fractal contact models is discussed and emphasized. It is demonstrated that the plasticity index diverges to infinity when the lower cut-off length of the system vanishes, whereas the slope of the real contact area versus normal load curve tends to zero. In this limit case stochastic and fractal approaches predict either a vanishing real contact area or an infinite normal pressure, regardless of the asperity deformation assumptions. On the other hand, when a non-zero lower cut-off length exists, finite contact predictions can be obtained by applying the contact models at that scale length. This possibility is numerically investigated and, comparing contact results of several models, it is shown that the plasticity index plays an important role for the characterization of the interface mechanical response.     

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.     

A deterministic model for line-contact partial elastohydrodynamic lubrication

A deterministic model for partial elastohydrodynamic lubrication (EHL) is presented in this paper. The modelling methodology adopts some of the concepts used in the stochastic modelling of partial EHL and some of the procedures for deterministic calculation of asperity pressures. The model is shown to be capable of simulating the basic process of asperity interaction and solid-to-solid contact within an EHL conjunction of rough surfaces. Deterministic results of transient partial EHL in line contacts are obtained when one pair or multiple pairs of asperities collide. The model may help to gain a fundamental understanding of the transient behaviour of asperity interactions in lubricated concentrated contacts of rough surfaces. Asperity pressures may be calculated more accurately than the conventional analyses under dry and static contact conditions. The work represents a first attempt in deterministic modelling of tribo-contacts operating in the mixed regime of micro-EHL and boundary lubrication. Future work will aim at developing more realistic models incorporating factors such as three-dimensional asperity contacts, asperity plastic deformation, thermal effects and the effect of tribo-chemistry.     

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.     

The Subsurface Stress Field Created by Three-Dimensionally Rough Bodies in Contact with Traction

A numerical simulation technique for calculating the complete subsurface stress field for three-dimensionally rough bodies in sliding contact is described. The stresses are calculated using real digitized three-dimensional surface profiles. The effects of the surface roughness and the sliding friction are presented. Using an existing contact simulation code, the digitized surfaces are mathematically pressed together and the real areas of contact and the asperity pressures are calculated. The surfaces are assumed to remain elastic throughout the contact simulation process. The shear forces at the asperity contact interfaces are assumed to be proportional to their calculated normal pressures. The subsurface stresses are then determined with these known normal and tangential forces at the surface.     

The residual stress distribution in a plastically deformed model asperity

When rough metallic surfaces come into contact, plastic deformation may occur locally, even at the lightest loads. This plastic deformation is thought to be an important element in a wide range of contact failure mechanisms, including fatigue and nearly all forms of wear.In this paper a simple model of asperity plastic deformation is presented. The model is based on slip line field theory and is used to calculate residual and full-load stress distributions at fully plastic asperity contacts for normal and moderate tangential loads.Measurements of surface residual stress were carried out using two different techniques on a range of plastic contacts of various materials and geometries. The results show agreement with the main predictions of the theory.     

Modelling of static friction in rubber–metal contact

A static friction model for contact between rough rubber and metal surfaces is developed. This model is based on the contact of a viscoelastic–rigid asperity couple. Single asperity contact is modelled in such a way that the asperities stick together in a central region and slip over an annulus at the edge of the contact. The slip area increases with increasing tangential load. Consequently, the static friction force is the force when the slip area is equal to the contact area. Using the model, the traction distributions, contact area, tangential and normal displacement of two contacting asperities are calculated. The single asperity model is then extended to multi-asperity contact, suitable for rough surfaces. This model allows calculation of the above-mentioned parameters for two rough surfaces (a rubber and a metal one) subjected to normal and tangential loads. A parametric study will be presented. The results are qualitatively in good agreement with those found in literature.     

Random process model of rough surfaces in plastic contact

The plastic contact of a rough surface and a hard, smooth flat is analyzed by modeling the rough surface as an isotropic, Gaussian, random process. The applicability of this model to the contact of two rough surfaces is discussed, and it is shown that the model is appropriate.It is not necessary to analyze interactions of asperity pairs, with the attendant questions of their misalignment, the shape of their caps, etc. Instead, a model involving the interaction of the continuous surfaces is developed, which implicitly takes into account these geometrical factors, as well as allowing for the possibility of the coalescence of microcontacts as the normal pressure is increased.Approximate relations between the density of finite contact patches, their mean area and mean circumference, and the normal pressure/hardness ratio are derived. These relations depend not only on the density and height distribution of maxima, but also on the shape of the Power Spectral Density of the surface. Many surfaces of interest are likely to give rise to multiply-connected contact patches at all except very high separations. The density of holes appearing within the contact patches as well as their area is estimated.Results are derived for surfaces that may be partitioned into two components, one with a large r.m.s. value and a narrow roughness spectrum, and the other with a small r.m.s. value and an arbitrary spectrum. For these surfaces, the density of holes at small separations becomes equal to the density of finite contact patches; the area of the holes remains small, however. It is conjectured that for surfaces that may not be partitioned in this manner, conventional models of contact are inapplicable. Specifically, the contact patches are likely to be perforated by holes at all separations, the hole area being a significant fraction of the contact area. Unit events such as the contact or collision of asperities also appear to become meaningless     

考虑蠕变特性的橡胶微凸体与金属表面接触特性分析

金属-橡胶接触广泛存在于密封结构中,密封接触表面上微凸体间的相互作用会直接影响整个密封界面的接触特性,进而影响其密封性能。基于粗糙密封界面的单个微凸体,考虑橡胶的蠕变特性,采用理论分析和仿真研究相结合的方式研究橡胶微凸体与金属表面的接触特性。通过橡胶蠕变特性的实验结果,构建橡胶蠕变计算模型;构建半球微凸体与金属平板间的有限元模型,进行考虑蠕变特性的仿真,分析其接触特性,并与Hertz接触理论的计算值进行对比。结果表明：在蠕变阶段,接触半径、法向变形量和最大等效蠕变应变均随蠕变时间的增加而增大,最大接触压力随蠕变时间增大而减小,这均可能导致密封性能的下降;随压力载荷的增大,接触半径、法向变形量、最大接触压力和最大等效蠕变应变均增大,但增大的趋势逐渐减小;橡胶微凸体与金属表面间的等效模量随蠕变时间的增加而减小,随压力载荷增大而增大。