摘要
基于纳米摩擦能耗理论,利用分子动力学方法建立了公度接触下支撑刚度梯度变化的石墨烯层间摩擦力模型,分析了基底质心刚度和支撑刚度梯度变化对基底和薄片各接触区摩擦能耗的贡献.结果表明:软边界区始终贡献驱动力;硬边界区贡献的摩擦力最大,且随着支撑刚度的增大,硬边界区对总摩擦的贡献比也越高.各接触区的摩擦力是薄片和基底之间的褶皱势和接触区产生的法向变形差两部分的共同作用.前者是公度接触下阻碍滑移的界面势垒和刚度梯度方向上不同刚度支撑原子热振动引起的势梯度;后者是接触边界过渡区两侧原子的非对称变形和自由度约束突变引起的非平衡边界势垒相耦合的结果.本文对研究公度接触下刚度梯度支撑的纳米器件的相对运动规律有指导意义.
According to the molecular dynamics simulations and the mechanism of energy dissipation of nanofriction, we construct a model system with a flake sliding in commensurate configuration on a monolayer suspended graphene anchored on a bed of springs. The system is to analyze the contributions of different regions(T1–T7) of the graphene flake to friction force, with the substrate characterized by different stiffness gradients and midpoint stiffness.The results indicate that the soft region of contact(T1) always contributes to the driving force, whereas the hard region(T7) leads to the biggest friction force on all column atoms of the flake. Moreover, as the support stiffness increases, when the stiffness gradient and the midpoint stiffness are equal to 1.34 nN/nm^2 and 12 nN/nm, respectively,the contribution ratio of T7 to the total friction increases from 33% to 47%, which is approximately 4–15 times greater than those of each column atoms in T3–T6. The results also indicate that the energy barrier decreases with the increase of support stiffness along the stiffness gradient direction of the substrate, which induces the resistance forces on the relative motion to decrease. Meanwhile, the amplitude of the thermal atomic fluctuation is higher in the softer region while lower in the harder one. This difference in amplitude leads to the considerable potential gradient that ultimately causes the driving force. Finally, for a given point at the end of the flake(T1 or T7), the intensity of the van der Waals potential field is mainly determined by the nearest substrate atoms at that point. Part of these nearest atoms lie inside the contact region while the others do not. Consequently, the thermal vibration of the atoms inside the contact region is different from that of the atoms outside the confinement. The different thermal vibrations induce the greater edge barriers. In addition, T1 lies in the soft edge region and T7 in the hard one. As a result, the normal deformations of these two regions are always different, and therefore they also generate the driving force.At these points, the results reported here suggest that the friction force in each contact region is caused by the coupling of the energy barrier and the elastic deformation between the graphene surfaces. The former contribution,i.e.the energy barrier, includes the interfacial potential barrier in commensurate state which is against the sliding of the surfaces with respect to each other, and the potential gradient caused by the different vibration magnitudes of the substrate atoms against the different spring stiffness in the direction of stiffness gradient. The latter contribution, i.e.the elastic deformation, is the unbalanced edge energy barrier resulting from the asymmetrical deformation and the different degrees of freedom between the edge atoms of the slider and atoms inside and outside the contact area of the substrate. Results of this paper are expected to be able to provide theoretical guidance in considering the influence of stiffness gradient on friction between commensurate surfaces and in designing the nanodevices.
作者
董赟
段早琦
陶毅
Gueye Birahima
张艳
陈云飞
Dong Yun;Duan Zao-Qi;Tao Yi;Gueye Birahima;Zhang Yan;Chen Yun-Fei(School of Mechanical Engineering,Jiangsu Key Laboratory for Design and Manufacture of Micro-Nano Biomedical Instruments,Southeast University,Nanjing 211189,China;School of Mechanical and Electronical Engineering,Lanzhou University of Technology,Lanzhou 730050,China)
出处
《物理学报》
SCIE
EI
CAS
CSCD
北大核心
2019年第1期210-219,共10页
Acta Physica Sinica
基金
国家自然科学基金(批准号:51665030
51435003
51575104)
甘肃省自然科学基金(批准号:17JR5RA114)
国家留学基金(批准号:201708625074)
东南大学优秀博士学位论文基金(批准号:YBJJ1749)资助的课题~~
关键词
分子动力学
纳米摩擦
公度接触
刚度梯度
molecular dynamics simulations
nanofriction
commensurate configuration
stiffness gradient