致; 2、 pair_coeff 中原子类型要与 data 文件中的 atom types 数目一致,在 atom s 中对第二列原子类型数也要做相应的修改; 3、在 region 区域定义中,注意拉伸方向的限制范围的表示方式; 4、 velocity 中温度的数值是整数型不是浮点型; 5、 fix 中系统的选取一般是 nvt ,在弛豫时温度的设定可以适当低些以便于更快的达到能量最低状态。我在这里 Tstart Tstop 设为 0.01 K ,Tdamp 取为 0.1 ; 6、 compute 应力 strec 时(in 文件), 指针数组 x 方向即为 1, y 方向即为 2, z 向为 3; 7、想得到什么输出结果可以在 thermo_ style 中设置; 8、 dump 输出图片格式应为“ dump.*.cfg ”,后面的坐标输出要严格根据 lammps manual 设定,输出步一般可以比循环运行的 step 多 1或者与之相等。 9、对 graphene 在 pair_style 用的势函数是 tersoff ;对于 T ,在 pair_styl e 用的势函数则为 rebo 即 the second generation of Brenner 。相应的 in文件在附件中。三、计算结果 1、在 win7 下的计算: 打开文件所在路径: cmd/cd MS5.5/pare /lmp_win -in “ in 文件名”,然后点回车键,坐等计算结果就 OK ; 2、 Lammps 下的计算: cd 打开文件所在目录, mpd&( 打开 mpd 并行运算), mpirun -np 4 lmp_g++<xx(in 文件名)四、结果分析: 最后的 T 应力应变如图: 相应的最小二次拟合为: 石墨烯的应力应变曲线: 相应的最小拟合曲线: graphene-t 应力应变曲线:
全文预览
