纳米级随机粗糙表面微观滑动摩擦力的计算研究

表面形貌很大程度上决定了摩擦副的摩擦性能, 而所有的表面都不可能是绝对光滑的.由于摩擦表面形貌的随机性, 决定了实际的摩擦过程具有随机性的特点, 因此为了获得与随机形貌对应的摩擦特性, 建立合理的随机摩擦模型是必要的. 本文基于Lennard-Jones势能建立了纳米级随机粗糙表面和原子级光滑的刚性平面间的随机摩擦模型; 模型中, 界面势能由法向载荷和界面间平衡距离决定.通过数值计算的方法, 推导了微观滑动摩擦力的计算公式和摩擦力与法向载荷之间的关系. 研究结果表明摩擦力随着法向载荷的增加而增加, 但不是线性增长. 结果也说明界面间的表面势能可能是微观摩擦力的本质起源. 关键词： 随机粗糙表面 Lennard-Jones势能 微滑动摩擦力 微摩擦     

基于非连续能量耗散的滑动摩擦系数计算模型

分析了界面摩擦状态下能量非连续耗散过程,建立了简化条件下晶体材料界面摩擦滑动摩擦系数计算模型.结果表明:在弹性接触状态下,滑动摩擦系数与载荷及实际接触面积无关,当实际接触面积接近名义接触面积时,滑动摩擦系数随载荷增加而减小.在缓慢滑动时,滑动摩擦系数随滑动速度的增高而缓慢增大,相对滑动速度愈高,滑动摩擦系数增大趋势愈显著.滑动摩擦系数随晶格常数的增加而降低,而当晶格常数较大时,其变化对滑动摩擦系数影响较小.同时,滑动摩擦系数随原子的可能温升增加而增大.研究结论对工程应用及相关的理论研究具有一定的参考意义. 关键词： 滑动摩擦系数 非连续能量耗散 界面摩擦     

摩擦微观能量耗散机理的复合振子模型研究

提出无磨损界面摩擦微观能量耗散机理的复合振子模型，指出滑动摩擦过程同时存在整体做低频弹性振动的宏观振子和界面原子受激励产生热振动的微观振子，并在此基础上分析了宏观振子和微观振子对摩擦能量耗散的不同影响. 通过对界面原子的动力学分析，指出摩擦过程界面激励力的频率是能量转换的关键:在平衡力作用阶段，界面作用力的频率趋于零，因而可以直接作用到每个原子，力的作用效果是整体和均匀的；在失稳跳跃阶段,由于界面激励力的频率极高，造成摩擦界面原子获得的能量分布很不均匀，从而产生不可逆的能量耗散过程. 与目前通用的独立振子模型比较，复合振子模型能够更准确描述摩擦能量耗散过程，可为摩擦控制提供理论指导. 关键词： 摩擦 能量耗散机理 复合振子模型 独立振子模型     

界面摩擦过程非连续能量耗散机理研究

结合无磨损界面摩擦微观能量耗散机理的复合振子模型,运用量子理论建立了微观能量耗散的量子力学模型.分析表明：在滑动过程中,当界面原子从一种平衡态跳跃至另一种平衡态时,摩擦功以离散形式耗散为界面原子热振子,且界面吸收能量的能力是离散的;高能态界面较低能态界面吸收能量的能力强,表现为易于吸收界面势能.界面原子吸收和释放能量的离散性在宏观上表现为摩擦功耗散的非连续性,为从微观角度解释无磨损界面摩擦状态周期性变化提供了理论基础. 关键词： 摩擦 非连续能量耗散 复合振子模型     

声子摩擦能量耗散机理研究

以界面摩擦为研究对象,分析了黏滑过程中的能量积累和耗散问题.基于晶格热动力学理论,通过分析界面原子在周期性势场中跳跃前后的势能差,推导了界面原子温升公式.理论表明,界面温升与摩擦系统的接触状态和材料特性有关,界面交互势能是其中影响较大的因素之一.在滑动阶段初期,由于界面原子处于非热平衡状态,晶格的热振动将通过激发出新声子而耗散能量,从而使得非热平衡向平衡状态转变.通过引入量子力学和热力学理论,分析了界面摩擦能量的耗散规律.结果表明,当声子振动频率较大时,黏着阶段存储于界面振子上的弹性势能在滑动阶段就很快完全耗散,耗散时间远小于滑动阶段的时间. 关键词： 界面摩擦 黏滑 声子 温升     

温度对微界面摩擦影响的研究

以微观界面摩擦为研究对象,分析了温度变化对材料摩擦性能的影响.基于Towle剪切强度-温度经验公式和晶格热动力学理论,推导出摩擦力与温度之间的理论计算公式.理论分析表明:当界面温度低于材料的德拜温度时,摩擦力随着温度的增加而降低.理论计算结果与原子力显微镜实验结果对比,发现二者趋势一致,表明本文提出的理论和方法可行. 关键词： 界面摩擦 真实接触面积 温度 摩擦力     

基底支撑刚度梯度变化对石墨烯层间摩擦力的影响

基于纳米摩擦能耗理论,利用分子动力学方法建立了公度接触下支撑刚度梯度变化的石墨烯层间摩擦力模型,分析了基底质心刚度和支撑刚度梯度变化对基底和薄片各接触区摩擦能耗的贡献.结果表明:软边界区始终贡献驱动力;硬边界区贡献的摩擦力最大,且随着支撑刚度的增大,硬边界区对总摩擦的贡献比也越高.各接触区的摩擦力是薄片和基底之间的褶皱势和接触区产生的法向变形差两部分的共同作用.前者是公度接触下阻碍滑移的界面势垒和刚度梯度方向上不同刚度支撑原子热振动引起的势梯度;后者是接触边界过渡区两侧原子的非对称变形和自由度约束突变引起的非平衡边界势垒相耦合的结果.本文对研究公度接触下刚度梯度支撑的纳米器件的相对运动规律有指导意义.     

高通量计算二维材料界面摩擦

建立了基于第一性原理方法研究二维材料界面摩擦的高通量计算程序,该程序实现了自动化批量建模、批量提交任务、多任务并发计算,以及计算结果自动收集、处理和图像绘制,使用该程序可以节省时间.采用此程序计算了不同层间距离下双层氮化硼和双层石墨烯的滑移势能面,及层间界面摩擦力和摩擦系数.研究发现,随着层间距离减小,双层氮化硼界面的平均摩擦力近似线性增大,摩擦系数为0.11—0.17,双层石墨烯界面摩擦力先增大后减小再增大,其摩擦系数在12 nN载荷下达到最小值(0.014),这些结果与已有研究结果一致,验证了该计算程序的可靠性.此外还研究了表面氢化和氟化对双层氮化硼界面摩擦的影响,发现氟化氮化硼/氮化硼界面的摩擦系数更低.     

最大静摩擦力初探

为什么最大静摩擦力总是大于滑动摩擦力呢?摩擦力是分子力的表现,当两个面相互接触时,两个面的部分表层分子产生了相互作用,就会出现接触力(一般来说,接触力是倾斜的)。把这个接触力正交分解,法向作用为弹力,切向为摩擦力.摩擦力的大小f=μ×N,N为弹力,μ为摩擦因数,它与两个因     

一类厄尔尼诺-南方涛动耦合振子动力学模型的震荡近似解

研究了一类厄尔尼诺和南方涛动(ENSO)耦合振子动力学模型. 利用奇摄动理论的参数变值法和平均法, 得到了对应ENSO耦合振子模型方程的震荡近似解. 关键词： ENSO振子 奇摄动 近似解     

Regular and reverse nanoscale stick-slip behavior: Modeling and experiments

We recently proposed a new nanoscale friction model based on the bristle interpretation of single asperity contacts. The model is mathematically continuous and dynamic which makes it suitable for implementation in nanomanipulation and nanorobotic modeling. In the present paper, friction force microscope (FFM) scans of muscovite mica samples and vertically aligned multi-wall carbon nanotubes (VAMWCNTs) arrays are conducted. The choice of these materials is motivated by the fact that they exibit different stick-slip behaviors. The corresponding experimental and simulation results are compared. Our nanoscale friction model is shown to represent both the regular and reverse frictional sawtooth characteristics of the muscovite mica and the VAMWCNTs, respectively.     

Friction and normal forces of model friction modifier additives in simulations of boundary lubrication

ABSTRACT

Interaction forces between solid surfaces are often mitigated by adsorbed molecules that control normal and friction forces at nanoscale separations. Molecular dynamics simulations were conducted of opposing semi-ordered monolayers of united-atom chains on sliding surfaces to relate friction and normal forces to imposed sliding velocity and inter-surface separation. Practical examples include adsorbed friction-modifier molecules in automatic transmission fluids. Friction scenarios in the simulations had zero, one, or two fluid layers trapped between adsorbed monolayers. Sliding friction forces increased with sliding velocity at each stable separation. Lower normal forces were obtained than in most previous nanotribology molecular simulations and were relatively independent of sliding speed. Distinguishing average frictional force from its fluctuations showed the importance of system size. Uniform velocities were obtained in the sliding direction across each adsorbed film, with a gradient across the gap containing trapped fluid. The calculated friction stress was consistent with measurements reported using a surface forces apparatus, indicating that drag between an adsorbed layer and trapped fluid can account sufficiently for sliding friction in friction modifier systems. An example is shown in which changes in molecular organisation parallel to the surface led to a large change in normal force but no change in friction force.     

Molecular origin of friction

The wearless friction originating from molecular interactions has been discussed in this paper. We find that the frictional properties are closely related to the structural match of two surfaces in relative motion. For the surfaces with incommensurate structure and week inter-surface interaction, zero static and kinetic friction can be achieved. In a sliding considered as in a quasi-static state, the energy dissipation initiates when interfacial particles move in a discontinuous fashion, which gives rise to a finite kinetic friction. The state of superlubricity is a result of computer simulations, but the prediction will encourage people to look for a technical approach to realizing the state of super low friction.     

Static friction of sinusoidal surfaces: a discrete dislocation plasticity analysis

Discrete dislocation plasticity simulations are carried out to investigate the static frictional response of sinusoidal asperities with (sub)-microscale wavelength. The surfaces are first flattened and then sheared by a perfectly adhesive platen. Both bodies are explicitly modelled, and the external loading is applied on the top surface of the platen. Plastic deformation by dislocation glide is the only dissipation mechanism active. The tangential force obtained at the contact when displacing the platen horizontally first increases with applied displacement, then reaches a constant value. This constant is here taken to be the friction force. In agreement with several experiments and continuum simulation studies, the friction coefficient is found to decrease with the applied normal load. However, at odds with continuum simulations, the friction force is also found to decrease with the normal load. The decrease is caused by an increased availability of dislocations to initiate and sustain plastic flow during shearing. Again in contrast to continuum studies, the friction coefficient is found to vary stochastically across the contact surface, and to reach locally values up to several times the average friction coefficient. Moreover, the friction force and the friction coefficient are found to be size-dependent.     

Frictional effect of ultrasonic-vibration on upsetting

The ultrasonic-vibration ring compression test and finite element analysis were performed on aluminum alloy specimens to explore the frictional effect of superimposing ultrasonic-vibration during upsetting. The extrapolated compression test was first adopted to obtain the frictionless material properties for finite element analysis. Experimental results of extrapolated compression test also indicate that ultrasonic-vibration can reduce the compressive force when friction is eliminated and can increase the temperatures of a material at the same time.The following results of the hot extrapolated compression test and the hot ring compression test reveal that increasing temperature by ultrasonic-vibration may reduce the flow stress and increase the interfacial friction. Finally, finite element analysis was conducted to derive the friction calibration curves and to evaluate the friction factor.     

基于滑动势能面的二维材料原子尺度摩擦行为的量化计算

近年来,二维材料优异的摩擦特性成为人们关注的焦点,然而目前缺乏理论上对其摩擦力进行快速、有效、精确的计算预测方法.本文提出采用密度泛函理论计算真实体系的滑动势能面,利用得到的"数值型势能面"替代传统的解析势函数,并结合Prandtl-Tomlinson模型,量化求解具有复杂形状势能面的真实二维材料体系的摩擦行为.基于该方法,揭示了原子力显微镜实验中观察到的石墨烯Moire纹超晶格结构的双周期"黏-滑"摩擦现象;理论预测了二维材料异质结构的层间超低摩擦现象,相对于同质材料,其静摩擦力和滑动摩擦力均成数量级降低,发现势能面起伏和驱动弹簧刚度均会影响层间相对滑动路径,进而对层间的摩擦行为产生影响.该方法同样可拓展到其他van der Waals作用主导的界面摩擦体系.     

Effects of non-flat contact and interference on self-assembled monolayers under sliding friction

The nanotribology mechanism of alkanethiol self-assembled monolayers (SAM) chemisorbed on a gold surface under a non-flat contact by a tilt plane was studied using molecular dynamics (MD) simulations. The molecular trajectories, tilt angles, normal forces, shear forces, and frictional coefficient of the SAM were evaluated during the friction and relaxation processes for various parameters, including the tilt angle of the slider, interference magnitude, and SAM length. At the nanoscale, the magnitude of interface interactional forces is strongly dependent on the magnitude of the contact area, not on the surface geometry. The contact area and the exerted normal force of the SAM increase with decreasing the tilt angle of the slider at the same contact interference. In contrast, the periods in both normal force and shear force are gradually delayed as the tilt angle of the slider increases. Once the contact interference increases, the normal force and shear force increase together. During the sliding friction process with a smaller tilt slider angle, SAM molecules can maintain a better collective ordered structure. Short SAM molecules are more sensitive to a compressive loading and react to a larger normal force under the same contact interference due to the deformation of a larger tilt angle and decrease in chain length. The friction coefficient of SAM is significantly more dependent on the tilt angle of the slider than the contact interference.