| 查看: 3127 | 回复: 2 | |||
[交流]
石墨烯单轴拉伸模拟
|
|
这是用airebo势函数，300K条件下的拉伸试验，下面附带有IN文件以及应力-应变曲线图。个人觉得问题有：1.驰豫过程中，温度是从300K开始升到两千多，然后快到驰豫结束时才又降为300K，不知道这个现象是否正常？；2 应力-应变曲线起伏波动比较大，没有一般文献的那么标准，不知道是否哪一快设置没有好？；3。拉伸断裂用VMD观察时，没有缓慢断裂而是比较突然，从应力应变曲线中可以相应看出。且断裂后，Y方向的长度是非常明显的变化，不知道是不是因为300K这个温度，使原子断裂后受热的扰动影响作用而导致。4.应变为零时，应力没法从零开始，大概是20GPA左右，不知道怎么可以能够改进（PS：平衡到最后，应力值就是保持在20几GPA) 希望各位能够不吝赐教。 dimension 3 > boundary p s s > > units metal > atom_style atomic > newton on > > neighbor 2.0 nsq > neigh_modify every 10 check yes > > read_data data.graphene > > # potentials > pair_style airebo 2.0 0 0 > pair_coeff * * CH.airebo C > > # groups > region Ledge block INF 3 INF INF INF INF units box > region Redge block 69 INF INF INF INF INF units box > > group left region Ledge > group right region Redge > group mobile subtract all left right > > thermo_style custom step lx ly lz press pxx pyy pzz pe temp > thermo 1000 > > # initial velocities > velocity left zero linear > velocity right zero linear > velocity mobile create 300.0 5481314 dist gaussian mom yes rot yes loop > local units box > > # nve equilibrate > dump 1 all atom 100000 grain.lammpstrj > > fix 0 right store/force > fix 1 left setforce 0.0 NULL 0.0 > fix 2 right setforce 0.0 NULL 0.0 > fix 3 mobile nvt temp 3000.0 300.0 0.1 > > timestep 0.001 > > min_style cg > minimize 1.0e-15 1.0e-5 1000 1000 > > run 1000 > unfix 3 > undump 1 > > # fracture > #fix 5 all temp/berendsen 300.0 300.0 0.1 > #fix 3 all nve > fix 3 all nvt temp 300.0 300.0 0.1 > run_style verlet > > dump 1 all atom 10000 grain.lammpstrj > dump_modify 1 append yes first yes sort id > > label loading > variable flag loop 150 > fix 6 all deform 1000 x erate 0.004 units box > run 1000 > unfix 6 > > next flag > jump in.graphene loading  应力-应变曲线 |
回复此楼» 猜你喜欢
面上意见出来了
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 石墨烯与氧化石墨的XRD,石墨烯是否合成成功？
已经有48人回复
- 用石墨烯和蛋白纤维合成的纸本事大
已经有7人回复
- PRL: 计算机模拟显示石墨炔性能胜过石墨烯-ZZ
已经有14人回复
- 石墨烯结构文件的生成
已经有9人回复
- 用第一性原理能否模拟石墨烯拉曼谱的d-band？
已经有4人回复
- 碳纳米管轴向切割制备石墨烯疑问？
已经有11人回复
- 美找到量产石墨烯的简单方法
已经有20人回复
- 【求助】怎么实现单轴拉伸模拟？
已经有11人回复
- 【讨论】材料拉伸测试时单轴向是什么意思？
已经有10人回复
2楼2013-08-12 13:51:27
3楼2017-05-10 09:19:30
投票: