基于分子动力学-格林函数法的微凸体接触数值分析

表面接触是摩擦的先决条件,其真实接触面积、压应力大小、空间分布等一直是接触力学关注的核心问题.采用分子动力学-格林函数法(GFMD)模拟粗糙面的接触过程,验证了其在大规模接触分析中的高效及准确性,同时探讨了由微球体组成的粗糙面的接触力学特性,并分析了分子尺度下的结果和传统力学模型计算结果的差异.结果表明,单个微凸体接触结果和分子动力学-格林函数法模拟所得非常接近,误差在5%以内.数值模拟发现,在微凸体高度符合高斯分布的情况下,接触面积和接触力成线性关系;在相同接触面积下,微凸体模型得出的接触力偏高,是上限值.微凸体模型没有考虑微凸体间的相互影响,实际是高估了弹性体的刚度;实际接触过程中微凸体相互影响,微凸体对临域形变影响尤其大,使接触区域更加离散.GFMD模型可以准确计算数十亿量级别分子、原子接触过程中真实接触面积及分布,为后续摩擦、滑移等分析提供可靠的参考.     

基于分子动力学-格林函数法的分形粗糙表面摩擦行为研究

任何物体间的表面摩擦均可看成是粗糙面间的摩擦,且大部分粗糙表面具有分形特征.为探究分形粗糙表面的摩擦行为,利用分子动力学-格林函数法(GFMD)建立微观分形粗糙表面模型,采用位移加载控制分形粗糙表面的接触和摩擦过程,并根据广度优先搜索算法识别接触团簇分布.之后分别计算原子尺度、接触团簇尺度和界面尺度下的最大摩擦系数与摩擦力,并利用影响矩阵法研究摩擦过程中接触团簇之间的相互作用,分析接触团簇之间的距离和面积对相互作用的影响.结果表明:在摩擦过程中,摩擦系数从小尺度到大尺度逐渐减小;摩擦力随位移呈现周期性波动,接触团簇并非同时达到最大摩擦力,而是发生局部滑移,整体滑移模型预测的摩擦力是分子模拟结果的上限值;所提出的影响矩阵法可以较好地模拟接触团簇之间的相互作用,利用影响矩阵计算得到的最大摩擦力与GFMD模型结果基本一致,而不考虑局部滑移影响计算得到的最大摩擦力比GFMD模型结果大20%,并且接触团簇之间的相互作用与距离成反比,与面积成正比.结果可为粗糙表面的界面分析和优化提供理论依据.