A Peierls perspective on mechanisms of atomicfriction

Stick-slip behavior observed from nanoscale asperity friction experiments are often modeled by the one-degree-of-freedom Tomlinson model, which is unable to explain the effects of lattice structure and interface defects, or by molecular simulations which suffer temporal limitations in modeling the velocity- and temperature-dependent frictional behavior. A Peierls-type model developed in this work views the atomic frictional process as the initiation and gliding passage of dislocations with diffused cores on the interface. As a consequence of loss-of-ellipticity instability, the occurrence of stick-slip behavior relies on the interaction among interface slip field, contact stress fields, and existing defects. The friction stress for commensurate interface under large contact area can be approximated from the Rice model of screw-dislocation nucleation from a planar crack tip. The spatially inhomogeneous nature of rate-limiting processes and the coupling effects between contact size and interface incommensurability are successfully determined, which cannot otherwise be tackled in the Tomlinsonmodel.     

Shakedown in elastic contact problems with Coulomb friction

Elastic systems with frictional interfaces subjected to periodic loading are sometimes predicted to ‘shake down’ in the sense that frictional slip ceases after the first few loading cycles. The similarities in behaviour between such systems and monolithic bodies with elastic–plastic constitutive behaviour have prompted various authors to speculate that Melan’s theorem might apply to them – i.e., that the existence of a state of residual stress sufficient to prevent further slip is a sufficient condition for the system to shake down.In this paper, we prove this result for ‘complete’ contact problems in the discrete formulation (i) for systems with no coupling between relative tangential displacements at the interface and the corresponding normal contact tractions and (ii) for certain two-dimensional problems in which the friction coefficient at each node is less than a certain critical value. We also present counter-examples for all systems that do not fall into these categories, thus giving a definitive statement of the conditions under which Melan’s theorem can be used to predict whether such a system will shake down.     

Homogenization in non-linear dynamics due to frictional contact

This work is devoted to a study of the classical homogenization process and its influence on the behavior of a composite under non-linear dynamic loading due to contact and friction. First, the general problem of convergence of numerical models subjected to dynamic contact with friction loading is addressed. The use of a regularized friction law allows obtaining good convergence of such models. This study shows that for a dynamic contact with friction loading, the classical homogenization process, coupled with an homogenization of the frictional contact, enables replacing the entire heterogeneous model by a homogenized one. The dynamic part of the frictional contact must be homogenized by modifying the dynamic parameter of the friction law. Modification of the dynamic parameter of the friction law is function of the type and regime of instability. A calculation of a homogenized friction coefficient is presented in view to homogenizing the static part of the frictional contact when the friction coefficient is not constant over the contact surface. Finally matrix and heterogeneities stresses in the heterogeneous models are identified by using the relocalization process and a frictional contact dynamic analysis of a homogeneous model.     

A hysteretic model of localized frictional contacts with instrumental stiffness

The steady-state dynamic response of a single-degree-of-freedom system comprising both a hysteretic element and a spring is considered. The Hertz–Cattaneo–Mindlin theoretical framework for modeling of local tangential contact with friction is applied in conjunction with the Masing model of hysteresis to describe the hysteretic behavior of the multiple localized frictional contact interface. The steady-state tangential displacement amplitude of a rigid body under harmonic tangential force excitation is approximately determined by means of the equivalent linearization technique, based on the harmonic balance principle. A special attention is paid to the evaluation of the frictional damping and the determination of the backbone curve of the Masing model from the dissipation-amplitude relation.

     

表面形貌变形对塑性成形滑动接触界面摩擦的影响

为了更好地理解塑性成形滑动接触界面的摩擦行为,构建了一种新型的摩擦试验装置,运用表面纹理化技术制备了两类表面形貌的1050铝材试件,在不同的接触压力和滑动速度条件下进行一系列拉伸摩擦试验.对试验前后试件三维表面形貌进行了测量;提取真实接触面积比、封闭空体面积比和开放空体面积比等三维表面参数,来描述试件表面形貌的变化.试验发现:摩擦系数随名义接触压力和滑动速度增加而逐渐减小;试件初始表面形貌对摩擦有明显的影响;试件表面形貌和参数随接触条件出现了规律性变化.基于机械流变模型的分析表明:随着试件表面形貌变形,不同的机理决定界面摩擦行为,摩擦系数对名义接触压力和滑动速度的依赖性可分别归因于微观塑性流体动压润滑效应和入口区流体动压牵引效应.     

Rapid sliding indentation with friction on a transversely isotropic thermoelastic half-space

A rigid insulated die slides at a constant sub-critical speed on a transversely isotropic half-space in the presence of friction. In a two-dimensional analysis of the dynamic steady-state, the coupled equations of thermoelasticity are invoked. All elements of the Coulomb friction model are strictly enforced, thus giving rise to auxiliary conditions, including two unilateral constraints.Robust asymptotic forms of an exact solution to a related problem with unmixed boundary conditions lead to analytical solutions for the sliding indentation problem. The solution expressions, abetted by calculations for zinc, show the role of frictional heating on the half-space surface. The effects of friction and sliding speed on contact zone size and location and average contact zone temperature are also studied.The analysis is aided by factoring procedures that simplify the complicated forms that arise in anisotropic elasticity. A scheme that renders expressions for roots of certain irrational functions analytic to within a single quadrature also plays a role.     

Double frictional receding contact problem for an orthotropic layer loaded by normal and tangential forces

With the increasing research in the field of contact mechanics, different types of contact models have been investigated by many researchers by employing various complex material models. To ascertain the orthotropy effect and modeling parameters on a receding contact model, the double frictional receding contact problem for an orthotropic bilayer loaded by a cylindrical punch is taken into account in this study. Assuming plane strain sliding conditions, the governing equations are found analytically using Fourier integral transformation technique. Then, the resulting singular integral equations are solved numerically using an iterative method. The weight function describing the asymptotic behavior of the stresses are investigated in detail and powers of the stress singularities are provided. To control the trustworthiness and correctness of the analytical formulation and to compare the resulting stress distributions and contact boundaries, a numerically efficient finite element method was employed using augmented Lagrange contact algorithm. The aim of this paper is to investigate the orthotropy effect, modeling parameters and coefficients of friction on the surface and interface stresses, surface and interface contact boundaries, powers of stress singularities, weight function and to provide highly parametric benchmark results for tribological community in designing wear resistant systems.

     

小腿残肢与接受腔界面的摩擦行为分析及测试

利用断层扫描数据,图像处理和反求技术建立了骨骼、软组织以及假肢接受腔的三维有限元模型,分别施加Heel Strike、Foot Flat、Mid-Stance、Heel Off和Toe Off五个典型步态时相的载荷,计算残肢皮肤在一个步态周期内的应力、滑移量变化,确定最大临界参数;借助UMT-II多功能摩擦磨损试验机模拟残肢皮肤和接受腔摩擦界面,研究皮肤在步态周期临界参数下的摩擦行为变化.有限元结果表明:一个步态周期内五个典型时相,软组织表面最大正应力和剪切力均发生在髌韧带处,在Heel off步态时相,最大临界正应力为384.3 k Pa,剪应力为102.1 k Pa;随着正压力的减小,摩擦力减小,摩擦系数增大;当正应力下降到39.5 k Pa临界值时,接触面发生相对滑动,摩擦系数达到最大值.摩擦学试验结果表明:法向载荷分别为17 N和7 N时,残肢皮肤均处于黏着状态;皮肤弹性变形越大,黏着程度越大.     

On indentation and initiation of fracture in glass

The influence of indenter elasticity on Hertzian fracture initiation at frictional dissimilar elastic contact has been examined experimentally and numerically. In flat float glass specimens initiation of cone cracks has been observed and fracture loads measured with steel and tungsten carbide indenters at monotonically increasing loading and during a load cycle. The observed effect of indenter elasticity on fracture loads was found to be qualitatively different from the one predicted by the Hertz contact theory. This discrepancy may be explained by the presence of interfacial friction. The friction coefficient between the indenters and the specimen was measured and a contact cycle at finite Coulomb friction has been analyzed numerically. The influence of the indenter elasticity and the friction coefficient on the surface maximum tensile stress has been investigated and the results concerning the influence of these parameters on the fracture loads as given based on a critical stress fracture criterion. The obtained computational results were found to be in better agreement with experimental findings as compared to the predictions based on the frictionless contact theory. A remaining quantitative discrepancy was attributed to the well-known fact that a Hertzian macro-crack initiates from pre-existing defects on the specimen’s surface. In order to account for the influence of the random distribution of these defects a Weibull statistics was introduced. The predicted critical loads corresponding to the 50% failure probability were found to be in close agreement with experimentally observed ones.     

A new computational algorithm for contact friction modeling of large plastic deformation in powder compaction processes

In this paper, the large deformation frictional contact of powder forming process is modeled based on a new computational algorithm by imposing the contact constraints and modifying the contact properties of frictional slip. A simple and efficient numerical algorithm is presented for imposing the contact constraints and frictional contact properties based on the node-to-surface contact technique to simulate the large deformation contact problem in the compaction process of powder. The Coulomb friction law is used to simulate the friction between the rigid punch and the workpiece by the use of penalty approach. A double-surface cap plasticity model is employed together with the nonlinear contact friction algorithm within the framework of large FE deformation in order to predict the non-uniform relative density distribution during large deformation of powder die-pressing. Finally, the numerical schemes are examined for accuracy and efficiency in modeling of a set of powder components.     

双材料单伸臂叠合梁的界面滑移与自振特性分析

将叠合梁划分为接触区和分离区,接触区界面间的摩擦作用会对叠合梁的滑移、刚度和自振频率产生影响。本文给出了单伸臂叠合梁在均布荷载和集中力作用下,考虑叠合界面摩擦作用的滑移应变和滑移分布的表达式;推导了考虑叠合界面间摩擦力及摩擦力产生的抵抗弯矩共同作用下的截面刚度。假设梁按等波长和等刚度两种形式自由振动,运用传递矩阵法推导出了单伸臂叠合梁自振频率的计算公式。最后,基于一个钢-混凝土组合与叠合试验梁的算例分析了滑移沿梁长的分布及自振频率。结果表明:随界面间摩擦系数的增大,接触区界面的滑移值逐渐减小;摩擦效应对叠合梁自振频率的影响随阶次的升高越来越明显。     

Comparison and analysis of two Coulomb friction models on the dynamic behavior of slider-crank mechanism with a revolute clearance joint

The objective of this study is to investigate the effects of the Coulomb dry friction model versus the modified Coulomb friction model on the dynamic behavior of the slider-crank mechanism with a revolute clearance joint. The normal and tangential forces acting on the contact points between the journal and the bearing are described by using a Hertzian-based contact force model and the Coulomb friction models, respectively. The dynamic equations of the mechanism are derived based on the Lagrange equations of the first kind and the Baumgarte stabilization method. The frictional force is solved via the linear complementarity problem (LCP) algorithm and the trial-and-error algorithm. Finally, three numerical examples are given to show the influence of the two Coulomb friction models on the dynamic behavior of the mechanism. Numerical results show that due to the stick friction, the slider-crank mechanism may exhibit stick-slip motion and can balance at some special positions, while the mechanism with ideal joints cannot.     

Interfacial stress analysis and prediction of debonding for a thin plate bonded to a curved substrate

This paper focuses on the analytical and numerical modeling of the interface between a rigid substrate with simple constant curvature and a thin bonded plate. The interfacial behavior is modeled by independent cohesive laws in the normal and tangential directions, coupled with a mixed-mode fracture criterion. The newly developed analytical model determines the interfacial shear and normal stress distributions as functions of the substrate curvature, during the various behavioral stages of the interface prior to the initiation of debonding. The model is also able to predict the debonding load and the effective bond length. In the numerical model the interface is modeled by zero-thickness node-to-segment contact elements, in which both the geometrical relationships between the nodes of the discretized problem and the interface constitutive laws are suitably defined. Numerical results and comparisons between the predictions of the two models are presented.     

One-dimensional dynamically consistent gradient elasticity models derived from a discrete microstructure: Part 1: Generic formulation

This paper is the first in a series of two that focus on gradient elasticity models derived from a discrete microstructure. In this first paper, a new continualization method is proposed in which each higher-order stiffness term is accompanied by a higher-order inertia term. As such, the resulting models are dynamically consistent. A new parameter is introduced that accounts for the nonlocal interaction between variables of the discrete model and of the continuous model. When this parameter is set to proper values, physically realistic behavior is obtained in statics as well as in dynamics. In this sense, the proposed methodology is superior to earlier approaches to derive gradient elasticity models, in which anomalies in the dynamic behavior have been found. A generic formulation of field equations and boundary conditions is given based on Hamilton's principle. In the second paper, analytical and numerical results of static and dynamic response of the second-order model and the fourth-order model will be treated.     

原子尺度摩擦行为的分子动力学模拟

应用大规模分子动力学方法,模拟了具有原子级光滑和原子级粗糙形貌的刚性球形探头与弹性平面基体的干摩擦行为,研究了无/有粘附条件下的载荷与摩擦力、载荷与真实接触面积,以及摩擦力与真实接触面积之间的关系,对纳米尺度下的摩擦行为规律进行了分析。几种系统的真实接触面积－载荷关系都与相应的连续力学接触模型定性的一致,它们分别是Hertz光滑表面接触模型、Greenwood-Williamson粗糙表面接触模型和Maugis-Dugdale粘着接触模型。无论是由光滑表面还是粗糙表面构成的摩擦系统,在无粘附条件下摩擦力与载荷成正比,而摩擦力与真实接触面积之间没有一个简单的关系;在粘附条件下摩擦力与真实接触面积成正比,而摩擦力与载荷之间表现为Maugis-Dugdale模型预测的亚线性关系。我们的研究表明,当表面作用从无粘附到粘附时,控制摩擦力的决定因素从载荷转变为接触面积,摩擦行为从载荷控制摩擦转变为粘着控制摩擦。     

Influence of Tribological Processes on a Chaotic Motion of a Bush in a Cylinder-Bush System

In this work periodic and chaotic dynamics of a bush being in a contact with a rotating shaft is analysed using the classical friction and abrasive wear models with the inclusion of the frictional heat generation. First, an analytical chaos prediction using the Melnikov method (without tribologic processes) is given, and then the analytical predictions are verified numerically. Then numerical analysis of the system including tribologic processes is carried out.     

含摩擦与碰撞平面多刚体系统动力学线性互补算法

基于接触力学理论和线性互补问题的算法, 给出了一种含接触、碰撞以及库伦干摩擦, 同时具有理想定常约束(铰链约束) 和非定常约束(驱动约束) 的平面多刚体系统动力学的建模与数值计算方法. 将系统中的每个物体视为刚体, 但考虑物体接触点的局部变形, 将物体间的法向接触力表示成嵌入量与嵌入速度的非线性函数,其切向摩擦力采用库伦干摩擦模型. 利用摩擦余量和接触点的切向加速度等概念, 给出了摩擦定律的互补关系式; 并利用事件驱动法, 将接触点的黏滞-滑移状态切换的判断及黏滞状态下摩擦力的计算问题转化成线性互补问题的求解. 利用第一类拉格朗日方程和鲍姆加藤约束稳定化方法建立了系统的动力学方程, 由此可降低约束的漂移, 并可求解该系统的运动、法向接触力和切向摩擦力, 还可以求解理想铰链约束力和驱动约束力. 最后以一个类似夯机的平面多刚体系统为例, 分析了其动力学特性, 并说明了相关算法的有效性.      

Non-linear micropolar and classical continua for anisotropic discontinuous materials

Non-linear Cosserat and Cauchy anisotropic continua equivalent to masonry-like materials, like brick/block masonry, jointed rocks, granular materials or matrix/particle composites, are presented.An integral procedure of equivalence in terms of mechanical power has been adopted to identify the effective elastic moduli of the two continuous models starting from a Lagrangian system of interacting rigid elements. Non-linear constitutive functions for the interactions in the Lagrangian system are defined in order to take into account both the low capability to carry tension and the friction at the interfaces between elements. The non-linear problem is solved through a finite element procedure based on the iterative adjustment of the continuum constitutive tensor due to the occurrence of some limit situation involving the contact actions of the discrete model.Differences between the classical and the micropolar model are investigated with the aid of numerical analyses carried out on masonry walls made of blocks of different size. The capability of the micropolar continuum to discern, unlike the classical continuum, the behaviour of systems made of elements of different size is pointed out. It is also shown that for anisotropic materials, even in the elastic case, the micropolar solution in general does not tend to the classical solution when the size of the elements vanishes.     

Correlation mechanism of friction behavior and topological properties of the contact network during powder compaction

The effect of friction behavior on the compacted density is significant, but the relationship between the topological properties of the contact network and friction behavior during powder compaction remains unclear. Based on the discrete element method (DEM), a DEM model for die compaction was established, and the Hertz contact model was modified into an elastoplastic contact model that was more suitable for metal-powder compaction. The evolution of the topological properties of the contact network and its mechanism during powder compaction was explored using the elastoplastic contact model. The results demonstrate that the friction behavior between the particles is closely related to the topological properties of the contact network. Side wall friction results in smaller clustering coefficient (CC) and excess contact (EC) in the lower region near the side wall. Corresponding to this phenomenon, the upper region near the side wall has more high-stress particles when the major principal stress threshold was considered, and the CC and EC are significantly higher than those in the other regions. This study provides a theoretical basis for improving powder compaction behavior.