面心立方晶格原子模拟的人工界面条件及其应用.pdfVIP

北京力学会第18届学术年会论文集：计算力学 面心立方晶格原子模拟的人工界面条件及其应用宰 方明。 唐少强1 李志辉2 (1北京大学工学院，100871；2中国空气动力研究与发展中心，621000) 摘要：本文研究固体晶格原子模拟的界面处理。通过匹配色散关系，我们提出了面心立方 晶格原子模拟的人工界面条件，并将其应用于纳米压痕的数值模拟。数值结果与现有理论 一致。 关键词：面心立方晶格，原子模拟，人工界面条件，纳米压痕 一、简介 界面条件是同体品格原子模拟及多尺度计算‘11的核心问题之一。时间历史积分‘21是 一种准确的界面条件，但其涉及到长时间的卷积计算。变分边界条件‘31能碌著减少计算 量，但是其依赖时间步长且需要通过优化反射系数来待定系数。新近提出的匹配界面条 件‘41在计算量极小且在低维问题中表现良好。本文将该方法推广到三维面心晶格，并将 其应用于纳米压痕的数值模拟‘51。 人工界面条件的构造 3．1面心立方晶格的动力学方程 考虑最近邻相互作用，无量纲化的动力学方程可以写为 K．一 ‰，o一㈠] ，o．[亡柚‰，-一[{；寸 Ko．I。=[o * ‰。=㈦；] m‰。=融 通过匹配长波极限处的色散关系，我们得到如下形式的界面条件 ‘国家自然科学基金91016027)资助 Ⅵ．23 北京力学会第18届学术年会论文集：计算力学 u。，0+cl ，+u。“。) p川。-+D川。·+D。+-，·+【：，。-．．·]=鼠以。。+且(玑扎¨+以 其中 l 笪oo _—— 2 4 l Cl= ●一 o—4—i-o 4 2 l 0 0 l _—— ，曰。=(一2喜／F一2：√虿jI]，口．= 4 3．3数值结果 将上述界面条件应用到纳米眶痕的原子模拟中，所得结果与理论预测相一致。 I — o 一 图1．纳米压痕过程中的载荷．位移曲线。 三、结论 本文构造了面心立方晶格原子模拟的人工界面条件并将其用于纳米压痕的数值模 拟中。模拟所得结果与理论值一致。 参考 文献 KamLiu，Eduard S．Park．NanoMechanicsandMaterials．科学出版社，2006． 1 Wing G．Karpov，Hardold in and 2 Eduard Kam a1．Multiscaleconditions G．Karpov，WingLiu，et bou