大家好,我是小马老师。
本文介绍lammps模拟计算界面结合力的一个案例。
本案例参考论文《The interfacial mechanical properties of functionalized graphene- polymer nanocomposites》,对论文作者表示感谢。
模型共分2层,上层为聚乙烯模型,下次为单层石墨烯,模拟过程共分3部分:建模、弛豫、拉伸。
(1)建模
使用MS建立聚乙烯和石墨烯模型,具体方法可参考之前的推文,使用build layer方法建立聚乙烯和石墨烯组合模型,为方便压缩,本案例中共建2层石墨烯。
msi2lmp命令将layer.car转换为layer.data文件。
(2)弛豫
参考论文弛豫过程,先进行300K、1atm的npt系综下驰豫,然后进行300Knvt系综驰豫,使聚乙烯达到合适的密度。
(3)拉伸
固定聚乙烯上层,以0.01A/ps的速度向下移动石墨烯,使石墨烯与聚乙烯分离,并记录石墨烯受力。
模拟得到的最大拉力为13.38nN,界面最大应力为702MPa,与论文中的727MPa几乎一致。
以上案例代码下载见公众号。
公众号案例代码以及学习交流,请加QQ群:754749935
扫描关注微信公众号:lammps加油站。