lammps案例：空位缺陷对石墨烯力学性能影响模拟

大家好，我是小马老师。
本文介绍使用lammps软件模拟带有空位缺陷的石墨烯拉伸过程。
石墨烯缺陷包含单空位（SV）、双空位（DV）、SW等多种缺陷，这些缺陷对石墨烯的力学性能产生一定的影响，相关的模拟论文也比较多，模拟方法基本类似。
在之前的文章中，介绍过使用velocity和deform两种方式对石墨烯进行拉伸，这也是两种比较常用的拉伸方式。
在velocity拉伸一文并没有给出应力计算公式，鉴于粉丝在后台对velocity拉伸应力计算咨询较多，因此，本文以单空位缺陷石墨烯拉伸为例，详细介绍velocity拉伸方式以及应力应变的计算方法，并给出全部in文件代码。
模拟的主要流程：
（1）建模
在lammps中使用velocity方式拉伸，在拉伸方向的边界条件不能周期性边界（p），本例拉伸方向为x方向，因此，x方向边界条件为s。
以lattice custom方式自定义石墨烯晶格，以此建立石墨烯结构。
也可以用其他软件建模，使用read_data读入模型。
通过以下命令删除中间部分原子，生成一个空位缺陷。
group del id 382
delete_atoms group del
（2）边界层的设置