Unloading an elastic-plastic contact of rough surfaces

A statistical model for the unloading of elastic-plastic contact of rough surfaces is presented for a single load-unload cycle. The hystereses of load-separation and load real contact area behavior are analyzed for a wide range of surface roughness and loading conditions. The residual topography of the unloaded rough surfaces is also analyzed and the new distribution functions of asperity heights and summit radii of curvature along with a corresponding GW residual plasticity index are presented. A new modified plasticity index (MPI) is suggested which considers the energy dissipation due to unrecovered plastic deformations. This MPI varies from zero for purely elastic contacts to unity for purely plastic contacts and hence, can better define the level of plasticity of contacting rough surfaces compared to the original GW plasticity index.     

粗糙颗粒的随机离散元模拟——随机法向接触定律

真实颗粒的力学性质会受到其随机粗糙表面的影响,然而在传统离散元模拟中通常假设颗粒具有光滑表面,因此有必要在定量考虑颗粒表面粗糙度的基础上改进离散元的接触模型。本文基于经典 Greenwood-Williamson(GW)模型通过理论分析和数值模拟提出了一种可以考虑颗粒表面粗糙度的法向接触定律;开发了基于 Newton-Raphson迭代的数值计算方法,通过输入颗粒重叠量和一系列表面粗糙系数计算总接触力;讨论了改进计算方法效率和准确性的相关问题。相对于 GW模型中接触关系的复杂积分表示,拟合得到新随机接触定律的表达式具有类似 Hertz定律的简单结构,只包含一个表征颗粒表面粗糙度标准偏差的新增参数,σ,可以方便的引入当前离散元模拟程序中进行计算。     

An adhesive wear model of fractal surfaces in normal contact

A generalized adhesive wear analysis that takes into account the effect of interfacial adhesion on the total load was developed for three-dimensional fractal surfaces in normal contact. A wear criterion based on the critical contact area for fully-plastic deformation of the asperity contacts was used to model the removal of material from the contact interface. The fraction of fully-plastic asperity contacts, wear rate, and wear coefficient are expressed in terms of the total normal load (global interference), fractal (topography) parameters, elastic–plastic material properties, surface energy, material compatibility, and interfacial adhesion characteristics controlled by the environment of the interacting surfaces. Numerical results are presented for representative ceramic–ceramic, ceramic–metallic, and metal–metal contact systems to illustrate the dependence of asperity plastic deformation, wear rate, and wear coefficient on global interference, surface roughness, material properties, and work of adhesion (affected by the material compatibility and the environment of the contacting surfaces). The analysis yields insight into the effects of surface material properties and interfacial adhesion on the adhesive wear of rough surfaces in normal contact.     

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.     

Adhesion in the contact of a spherical indenter with a layered elastic half-space

With the emergence of micro- and nano-technology, the contact mechanics of MEMS and NEMS devices and components is becoming more important. Thus it is important to gain a better understanding of the role of coatings and thin films on micro- and nano-scale contact phenomena, and to understand the interactions of measurement devices, such as an atomic force microscope (AFM), with layered media.More specifically, in this work the frictionless contact, with adhesion, between a spherical indenter and an elastic-layered medium is investigated. This configuration can be viewed as either a single contact model or as a building block of a multi-asperity rough surface contact model. As the scale decreases to the nano level, adhesion becomes an important issue. The presence of adhesion affects the relationships among the applied force, the penetration of the indenter, and the size of the contact area. This axisymmetric problem includes the effect of adhesion using a Maugis type of adhesion model. This model spans the range of the Tabor parameter between the JKR and DMT regions. The key parameters in this analysis are the elastic moduli ratio of the layer and the substrate, the dimensionless layer thickness, and the Maugis adhesion parameter. The results can be applied to a rigid or to an elastic indenter.     

粗糙表面统计接触模型的提出与发展

工程表面在微观尺度上是粗糙的，粗糙表面之间的接触状态对于多种物理现象都有重要影响，因此，粗糙表面间的接触建模方法一直是摩擦学领域研究的热点. 由Greenwood和Williamson提出的GW统计接触模型是最受认可的粗糙面接触模型，回顾了基于统计分析的粗糙面接触模型的发展，根据对GW模型主要缺点的改进，介绍了统计接触模型的研究现状，总结了统计接触模型未来可能的研究热点.      

Thermo-mechanical modelling of the contact between rough surfaces using homogenisation technique

In this study, thermo-mechanical behaviour of contacting rough surfaces has been modelled. Firstly, a numerical microscopic contact model that considers the properties of engineering surfaces has been developed. Geometrical characteristics of rough surfaces are deduced using the standard procedure for roughness and waviness parameter determination according to the so-called “motif” procedure. Secondly, an equivalent macroscopic contact model using a homogenisation technique has been presented. The interfacial behaviour of this model has been governed by the curves deduced from the microscopic model. The transition from microscopic to macroscopic scale was also validated.     

粗糙表面的分形特征与分形表达研究

得用触针轮廓仪和数据采集系统对磨削和车削表面的粗糙轮廓曲线进行了测量，并就粗糙表面的分形特征作了分析与讨论，同时还提出了粗糙表面的特征粗糙度概念及其定义，并将其用表面粗糙度水平的表述。     

基于共轭梯度法和快速傅立叶变换的粗糙表面微动接触数值研究

实际工程表面多为粗糙表面,研究粗糙表面的表面形貌对微动接触中压力和应力的影响具有重要意义。本文研究接触过程中法向载荷保持不变,切向载荷为周期性的交变载荷。首先,建立接触算法和模型,算法核心是利用共轭梯度法CGM(Conjugate gradient method)计算微动接触中的表面压力及切向应力并利用快速傅里叶变换FFT(Fast Fourier transform)加快计算速度。然后,对单峰表面、正弦表面和随机粗糙表面的接触进行数值研究。结果表明,表面幅值对切向载荷-位移曲线以及接触过程中的能量耗散有影响,表面幅值越大,相同切向载荷作用下产生的切向位移越大,能量耗散也越大。     

Lattice–Boltzmann simulations for analysing the detachment of micron-sized spherical particles from surfaces with large-scale roughness structures

Fully resolved numerical simulations of a micron-sized spherical particle residing on a surface with large-scale roughness are performed by using the Lattice–Boltzmann method. The aim is to investigate the influence of surface roughness on the detachment of fine drug particles from larger carrier particles for transporting fine drug particles in a DPI (dry powder inhaler). Often the carrier surface is modified by mechanical treatments for modifying the surface roughness in order to reduce the adhesion force of drug particles. Therefore, drug particle removal from the carrier surface is equivalent to the detachment of a sphere from a rough plane surface. Here a sphere with a diameter of 5 μm at a particle Reynolds number of 1.0, 3.5 and 10 are considered. The surface roughness is described as regularly spaced semi-cylindrical asperities (with the axes oriented normal to the flow direction) on a smooth surface. The influence of asperity distance and size ratio (i.e. the radius of the semi-cylinder to the particle radius, R_c/R_d) on particle adhesion and detachment are studied. The asperity distance is varied in the range 1.2 < L/R_d < 2 and the semi-cylinder radius between 0.5 < R_c/R_d < 0.75. The required particle resolution and domain size are appropriately selected based on numerical studies, and a parametric analysis is performed to investigate the relationship between the contact distance (i.e. half the distance between the particle contact points on two neighbouring semi-cylinders), the asperity distance, the size ratio, and the height of the particle centroid from the plane wall. The drag, lift and torque acting on the spherical particle are measured for different particle Reynolds numbers, asperity distances and sizes or diameters. The detachment of particles from rough surfaces can occur through lift-off, sliding and rolling, and the corresponding detachment models are constructed for the case of rough surfaces. These studies will be the basis for developing Lagrangian detachment models that eventually should allow the optimisation of dry powder inhaler performance through computational fluid dynamics.     

粗糙表面法向接触刚度的分形模型

提出了以往有关粗糙表面法向接触刚度理论研究工作的缺陷与不足,并在一定的前提假设下,基于球体与平面的接触理论和粗糙表面的分形接触理论,从理论上给出了具有尺度独立性的粗糙表面法向接触刚度分形模型,并进行了数字仿真研究。     

Effects of oxidation and surface roughness on contact angle

Contact angle is known to be a parameter that effects boiling. This study was undertaken to measure contact angle of high and low surface tension fluids on copper and aluminum surfaces.Data were taken for polished, oxidized, and rough surfaces. A simple, yet fairly accurate method of measuring the static equilibrium contact angle of a solid/liquid interface is presented. The principles of a line light source and tilting plate were modified and then combined in the design of this apparatus. The angles obtained and their variation with the solid surface properties were in good agreement with previously published data. The contact angle of distilled water o of the organic fluids and refrigerants tested were in the range of 2–5°. Roughness and oxidation reduce the contact angle. If the depth of the roughness is less than 0.5 μm contact angle. The apparatus is fairly simple in construction, is inexpensive, and has good reproductibity. The measured angles were then compared to those measured with the sessile drop method.     

Statistical model of nearly complete elastic rough surface contact

In the area of homogeneous, isotropic, linear elastic rough surface normal contact, many classic statistical models have been developed which are only valid in the early contact when real area of contact is infinitesimally small, e.g., the Greenwood–Williamson (GW) model. In this article, newly developed statistical models, built under the framework of the (i) GW, (ii) Nayak–Bush and (iii) Greenwood’s simplified elliptic models, extend the range of application of the classic statistical models to the case of nearly complete contact. Nearly complete contact is the stage when the ratio of the real area of contact to the nominal contact area approaches unity. At nearly complete contact, the non-contact area consists of a finite number of the non-contact regions (over a finite nominal contact area). Each non-contact region is treated as a mode-I “crack”. The area of each non-contact region and the corresponding trapped volume within each non-contact region are determined by the analytical solutions in the linear elastic fracture mechanics, respectively. For a certain average contact pressure, not only can the real area of contact be determined by the newly developed statistical models, but also the average interfacial gap. Rough surface is restricted to the geometrically-isotropic surface, i.e., the corresponding statistical parameters are independent of the direction of measurement. Relations between the average contact pressure, non-contact area and average interfacial gap for different combinations of statistical parameters are compared between newly developed statistical models. The relations between non-contact area and average contact pressure predicted by the current models are also compared with that by Persson’s theory of contact. The analogies between the classic statistical models and the newly developed models are also explored.     

A slip-line plasticity analysis of abrasive wear of a smooth and soft surface sliding against a rough (fractal) and hard surface

Abrasive wear of a soft and smooth surface sliding against a rough (fractal) and hard surface was analyzed by the slip-line theory of plasticity. The analysis is based on a slip-line model of a rigid spherical asperity (wear particle) plowing through a soft surface and removing material by microcutting. Integration of this single-contact model into a three-dimensional contact mechanics analysis of an abrasive surface exhibiting multi-scale roughness described by fractal geometry yielded relationships of the abrasive wear rate and wear coefficient in terms of the interfacial shear strength (adhesion), topography (fractal) parameters of the hard/rough surface, elastic–plastic material properties of the soft/smooth surface, and total normal load. Analytical results from the single-contact analysis provide insight into the deformation of a perfectly plastic material caused by the abrasive action of a rigid asperity/wear particle under different normal load and interfacial friction (adhesion) conditions. The dependence of the abrasive wear rate and wear coefficient on normal load (global interference), roughness of the abrasive surface, elastic–plastic material properties of the abraded surface, and interfacial shear strength (lubrication effect) is interpreted in the context of numerical results obtained for representative ceramic/ceramic, ceramic/metallic, and metal/metal sliding systems.     

Tabor数、粘着数与微尺度粘着弹性接触理论

微电子机械系统（ＭＥＭＳ）等领域的飞速发展,促使我们迈进了一个表面效应在许多现象 中占主导地位的研究领域．本文重点介绍在ＭＥＭＳ中经常遇到的微尺度粘着弹性接触的相关理论． 通过对两个无量纲数——Ｔａｂｏｒ数μ（以及其相应形式）和粘着数θ的分析,以及考虑它们对于粘 着力的影响,指出了粘着弹性接触理论中所隐含的尺度效应,随着特征尺度的减小,粘着弹性接触中 的表面效应愈加明显．     

A friction contact stiffness model of fractal geometry in forced response analysis of a shrouded blade

To investigate the nonlinear vibration behavior of a shrouded blade with friction dynamic contact interface, a friction contact stiffness model is proposed to describe the friction force at different rough interfaces and different normal loads. In the proposed model, the friction contact interface is discretized to a series of friction contact pairs and each of them can experience stick, slip, or separate states. Fractal geometry is used to simulate the topography of contact surfaces. The contact stiffness is calculated using the Hertz contact theory and fractal geometry, which is related to contact interfaces parameters including normal load, roughness, Young??s modulus, and Poisson??s ratio. The trajectory tracking method is used to predict the friction force and it is not necessary to judge the transition condition among stick, slip, and separate states. It is suitable for complicated periodic motion of the contact interfaces. The forced response of a real shrouded blade is predicted using the proposed model and the multi-harmonic balance method. The effect of surface roughness, initial normal load, and contact area on the forced response of a shrouded blade is studied. It is shown that contact stiffness increases with normal load and fractal dimension. The resonant amplitude is sensitive to the initial normal load and contact surface roughness. The response can be influenced by the contact area, which is an important parameter for blade designers.     

A slip-line plasticity analysis of sliding friction of rough surfaces exhibiting self-affine (fractal) behavior

A plasticity analysis of sliding friction of rough (fractal) surfaces sliding against smooth surfaces was developed based on a slip-line model of a rigid spherical asperity (wear particle) plowing and cutting through a soft semi-infinite medium. Solutions of the fraction of fully plastic asperity microcontacts responsible for the evolution of friction and energy dissipation were obtained in terms of the total normal load (global interference), interfacial adhesion characteristics, topography (fractal) parameters of the hard surface, and elastic–plastic material properties of the soft surface. This was accomplished by incorporating the slip-line model of a single microcontact into a friction analysis of sliding surfaces demonstrating multi-scale roughness. Numerical results provide insight into the effects of global interference (normal load), fractal parameters (surface roughness) of the hard surface, interfacial shear strength (adhesion), and material properties of the soft surface on plastic deformation at the microcontact level, global coefficient of friction, and frictional energy dissipated during sliding.