微观随机粗糙表面接触有限元模型的构建与接触分析

基于相关文献提出粗糙表面模拟方法,通过软件工具在ANSYS中建立微观粗糙表面接触有限元模型,利用该模型分析载荷对弹塑性变形和接触面积的影响。结果表明:随着正压力的增大,粗糙表面上不断地有微凸峰与平面发生弹性接触变形,接触斑点(或接触斑点群)的数目逐渐增加,斑点中心区域的弹性变形很快达到最大,微凸峰负荷变形的同时也使斑点四周区域受到挤压;初始接触时,轮廓高度较大的微凸峰率先发生弹性变形,随着压力的增大,金属材料所受应力达到屈服极限同时粗糙表面的弹性变形和塑性变形的集中区域不断增加,真实接触面积不断增大;接触区数目的增多和接触区面积的增加都可以导致接触面上真实接触面积增加;随着压力的增大,真实接触面积的增大并不是由于接触区数目的增多,而是微观接触区面积的增大。     

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.     

Deformation due to contact between a rough surface and a smooth ball

Theoretical and experimental results are presented to evaluate the deformation behavior of the contact between a real rough flat surface and a smooth ball. There are three deformation responses: plastic deformation of the asperities only, plastic deformation of the bulk only and combined plastic deformation of both the asperities and the bulk. The effects of the surface roughness and the Hertzian contact parameters on the effective contact pressure are presented. The experimental results confirmed the theoretical prediction very well. For a given Hertzian contact situation the surface roughness plays an important role in controlling the deformation behavior of the contacting surfaces. A criterion is presented to predict the deformation behavior of contacting engineering surfaces.     

Thermomechanical analysis of elastoplastic medium in sliding contact with fractal surface

The effects of mechanical and thermal surface loadings on deformation of elastic–plastic semi-infinite medium were analyzed simultaneously by using the finite element method. Rigid rough surface of a magnetic head and smooth surface of an elastic–plastic hard disk were chosen to perform a comprehensive thermo-elastic–plastic contact analysis at the head–disk interface (HDI). A two-dimensional finite element model of a rigid rough surface characterized by fractal geometry sliding over an elastic–plastic medium was then developed. The evolution of deformation in the semi-infinite medium due to thermomechanical surface loading is interpreted in terms of temperature, von Mises equivalent stress, and equivalent plastic strain. In addition to this, the effects of friction coefficient, sliding, and interference distance on deformation behavior were also analyzed. It is shown that frictional heating increases not only the contact area but also the contact pressure and stresses.     

Effect of the cut-off frequency on rough-point and flat-surface contacts?

In the past years, contact between two bodies has been studied from various ways that do not consider the cut-off frequency effect on the contact mechanism. This paper reports the correlation between rough point contact and flat surface contact at different cut-off frequencies of filter. The similarity and difference between the two types of contact mechanisms are presented for materials with linear or elasticperfectly plastic deformation. The conjugate gradient method (CGM) is used for analysing the rough point contact, while the rough flat surface contact is studied with an improved CGM in which the influence coefficient for the elastic deformation of the rough flat surface is obtained with finite element method. Numerical results show that for the above two types of contacts, their von Mises stress and maximum shear stress are greatly affected by the cut-off frequency of a high-pass filter. Moreover, a decrease in the cut-off frequency leads to an increase in the contact area and a decrease in the approach for the rough flat surface contact, while the opposite variations is for the point contact between rough bodies with the small radii.     

Criterion for surface contact deformation of metals

In most engineering applications, bulk plastic deformation of the surface is avoided. There is, however, no criterion for determining whether or not bulk plastic deformation occurs during the contact between rough surfaces. This paper presents a criterion for predicting the deformation behaviour of the contact between a real rough surface against a curved body. In order to verify the proposed criterion, experimental investigations were carried out for three different materials. It was found that the experimental results confirmed the proposed criterion very well. For a given operational contact condition, a contact system where bulk deformation does not occur can be designed using this criterion. If the criterion value C_m is less than 0.52 plastic deformation occurs on asperity level only and there is no bulk deformation, but if C_m is greater than 0.52 there are two possible deformation modes, i.e. plastic deformation of the bulk and plastic deformation in both the asperity and the bulk. Copyright © 2006 John Wiley & Sons, Ltd.     

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

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

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.     

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

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

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

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

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

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

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

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

Contact analysis of three-dimensional rough surfaces under frictionless and frictional contact

A methodology for surface and sub-surface stress calculation of nominally flat on flat rough surface contact has been developed. This methodology is applicable for both large area contact (Hertzian contact) and small area of asperity contact (point load contact) with and without surface friction. A total of nine rough surfaces are generated by the computer with specified standard deviation of surface heights, σ, of 0.3, 1.0 and 3.0 nm, and correlation length, β^*, of 0.1, 0.5 and 0.9 μm. Under the typical applied load at the magnetic head slider-disk interface, small numbers of contact points are obtained and the deformation is purely elastic. Since these contact points are scattered and isolated, asperity contact behaves like point load contact. As β^* becomes larger, more adjacent points will be in contact at a certain contact spot and this is especially true at small σ. All the cases of flat on flat rough surface contact yield maximum von Mises stress on and near the surface at both frictionless and frictional contacts; no local maximum occurs in the sub-surface. In general, the friction effect in the vicinity of contact point is to increase the stress magnitude, while outside this region it also alters the stress distribution. For a surface of small β^* and large σ at high load of 1000 times of the nominal pressure at the head-disk interface, the contact pressure reaches the hardness at a few contact points and plastic deformation takes place in the near surface.     

粗糙表面塑性变形对弹流润滑性能的影响

利用弹塑性动压润滑理论来分析线接触条件下粗糙表面的塑性变形对润滑性能的影响。首先计算总变形量,然后根据弹塑性理论对总变形量进行修正获得弹塑性变形量。计算中采用直接迭代法对雷诺方程进行数值求解,用多重网格法得到弹塑性变形量。为定量分析塑性变形的影响,采用余弦粗糙度来代替实际的粗糙表面,通过改变余弦粗糙度的峰高、波长和相位来表征粗糙度的变化,以此来分析粗糙度对弹塑性流体润滑的影响。计算结果表明:考虑塑性因素后,在接触的高压区膜厚变小,油膜分布曲线变得平坦,同时接触区域的油膜压力也明显变小;此外塑性变形因素的存在使得接触区的压力和油膜分布宽度增加,第二压力峰和最小膜厚点后移。     

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

利用弹塑性有限元和单纯形法求解弹塑性接触模型，分别模拟了屈服强度呈梯度变化的渗氮钢、未经处理的匀质材料和硬涂层材料粗糙表面的弹塑性接触行为。与未经处理的匀质材料相比，渗氮钢可承受更大接触载荷。在相同载荷作用下，渗氮钢表面粗糙峰接触面积较小，平均间距较大，接触体内材料不易发生屈服，从而显著提高接触性能。和硬涂层材料相比，渗氮钢接触体内等效von Mises应力分布平缓，没有应力突变。最后讨论了渗氮层和硬涂层的厚度对粗糙表面接触特性的影响。     

二维多粗糙峰涂层表面的弹塑性接触力学分析

应用有限元方法对二维多粗糙峰涂层表面的弹塑性接触力学行为进行了分析。对不同涂层材料弹性模量、不同屈服极限、不同涂层厚度及不同表面形貌的粗糙表面与刚性平面的弹塑性接触问题进行了模拟,分析了这几种因素对接触压力、接触面积、表面轮廓变形及应力场分布的影响。     

Asperity micro-contact models as applied to the deformation of rough line contact

Different statistical micro-contact models including Greenwood–Williamson (GW), Chang–Etsion–Bogy (CEB), Zhou–Maietta–Chang (ZMC), Kogut–Etsion (KE) and Jackson–Green (JG) are employed together with the bulk deformation of the bounding solids to predict dry rough line-contact characteristics such as the apparent pressure profile, contact width and real area of contact. The approach involves solving the micro-contact models and separation formulas simultaneously. Comparison of different contact models reveals that the use of elastic–plastic micro-contact models predicts a lower maximum normal pressure and a greater contact width and real contact area compared to the GW model. Further, based on the results of numerical simulations, useful relationships are provided for the prediction of the maximum contact pressure, contact width, real area of contact and pressure distribution.     

管路连接件粗糙表面接触力学与密封性能分析

为研究真实管道连接件在拧紧力矩作用下的接触状态和密封机制,从而有效控制管路密封性能,以真实管路连接件为研究对象,对粗糙接触表面真实接触形貌数据进行实测,通过逆向建模工程建立粗糙接触表面有限元模型,利用有限元仿真分析方法得出真实接触面积比与平均接触应力的关系。结果表明,在开始施加位移载荷时同时发生弹性与塑性变形,随着载荷的增加,发生塑性变形区域越来越大,形成真实接触区域部分也逐渐增加。基于逾渗理论对接触界面进行栅格化模型模拟,通过寻径算法得出接触界面达到密封时所需要的最小真实接触面积比,结合仿真结果得到管路连接件接触表面达到密封时所需要的最小接触应力值。研究结果为管路连接件密封性能仿真分析提供重要密封评判指标,也进一步通过仿真分析得到管接头密封条件下的最小拧紧力矩提供参考。     

A modified fractal microcontact model developed for asperity heights with variable morphology parameters

The cross-sections formed by the contact asperities of two rough surfaces at an interference are actually island-shaped, rather than having the commonly assumed circular or elliptic contour. These island-shaped contact area contours show fractal behavior with fractal dimension D_s of the two-dimensional profile. The three-dimensional surface fractal dimension for the topography of asperity heights is defined as D and the topothesy is defined as G. In Mandelbrot's study, the relationship between D and D_s was given as D = D_s + 1 if these two fractal dimensions are obtained before contact deformation. In the present study, D, G, and D_s are considered to be varying with the mean separation between two contact surfaces. The D–D_s relationships for the contacts at the elastic, elastoplastic, and fully plastic deformation regimes are derived and the inceptions of the elastoplastic and fully plastic deformation regimes are redefined using the equality of two expressions established in two different ways for the number of contact spots (N). A revised elastic–plastic contact model of a single fractal asperity is also proposed. The revised model shows that a fractal asperity behaves according to classical contact mechanics, but not those predicted by the MB model. The contact parameters, including the total force and the real contact area, were evaluated when the size distribution functions (n) for the three deformation regimes were available. The results indicate that both the D and D_s parameters in these deformation regimes increased with increasing mean separation.     

Modeling of fretting wear evolution in rough circular contacts in partial slip

This paper presents a numerical model that maps the evolution of contact pressure and surface profile of Hertzian rough contacting bodies in fretting wear under partial slip conditions. The model was used to determine the sliding distance of the contacting surface asperities for one cycle of tangential load. The contact pressure and sliding distance were used with Archard's wear law to determine local wear at each surface asperity. Subsequently, the contact surface profile was updated due to wear. The approach developed in this study allows for implementation of simulated and/or measured real rough surfaces and study the effects of various statistical surface properties on fretting wear. The results from this investigation indicate that an elastic–perfectly plastic material model is superior to a completely elastic material model. Surface roughness of even small magnitudes is a major factor in wear calculations and cannot be neglected.