计算物理 ›› 2012, Vol. 29 ›› Issue (5): 753-758.

• 论文 • 上一篇    下一篇

纳米接触行为的三维多尺度数值模拟

段芳莉, 郭其超   

  1. 重庆大学机械传动国家重点实验室, 重庆 400030
  • 收稿日期:2012-03-08 修回日期:2012-05-22 出版日期:2012-09-25 发布日期:2012-09-25
  • 作者简介:段芳莉(1971-),女,博士,副教授,主要从事分子模拟和多尺度研究,E-mail:flduan@cqu.edu.cn
  • 基金资助:
    国家自然科学基金(50875271);机械传动国家重点实验室科研业务费(SKLMT-ZZKT-2012MS18)资助项目

Multiseale Modeling of Three-dimensional Contacts

DUAN Fangli, GUO Qichao   

  1. State Key Laboratory of Mechanical Transmission, Chongqing University, Chongqing 400030, China
  • Received:2012-03-08 Revised:2012-05-22 Online:2012-09-25 Published:2012-09-25

摘要: 应用分子动力学与有限元耦合的桥域多尺度算法,模拟三维刚性球形压头与光滑基体表面的纳米尺度接触行为,并与全原子分子模拟结果比较.考察在一定载荷下的系统弛豫行为、两种模型桥接区位移和应力的连续性、法向力和接触面积随压头位移变化等,结果表明:一定外载荷下,桥域多尺度算法能较快达到平衡状态,且压头的振荡幅度更小,系统初始温度为0 K时该算法的相对误差最小.在准静态加载过程中,该算法能够将原子区的位移、应力等连续的过渡到连续介质区,具有较好的耦合效果;法向力-压头位移和接触半径-压头位移曲线几乎与分子模拟结果重合,表明算法具有较高的计算精度.

关键词: 桥域多尺度方法, 纳米接触行为, 多尺度模拟, 分子动力学模拟

Abstract: A bridging domain method,in which molecular models and continuum models are coupled,is performed to study threedimensional contact of a rigid spherical tip and a smooth fiat substrate.It is compared with full-atom molecular simulations.We focus on relaxation behavior under a given load,continuity of displacement and stress at bridging domain,normal force and contact radius.It shows that the bridging domain method quickly reaches equilibrium and has weaker oscillations under a given external load.It has small relative errors as initial temperature of system is at 0 K.The bridging domain method can transmit continuously displacement and stress from molecular domain to continuum domain,and has goad coupling effect.Normal force-displacement and contact radius-displacement are nearly coincide with molecular simulations.Its computational accuracy is high.

Key words: bridging domain method, nanoseale contact, muhiscale simulation, molecular dynamics simulation

中图分类号: