| 查看: 3106 | 回复: 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  应力-应变曲线 |
回复此楼» 猜你喜欢
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 石墨烯与氧化石墨的XRD,石墨烯是否合成成功？
已经有48人回复
- 用石墨烯和蛋白纤维合成的纸本事大
已经有7人回复
- PRL: 计算机模拟显示石墨炔性能胜过石墨烯-ZZ
已经有14人回复
- 石墨烯结构文件的生成
已经有9人回复
- 用第一性原理能否模拟石墨烯拉曼谱的d-band？
已经有4人回复
- 碳纳米管轴向切割制备石墨烯疑问？
已经有11人回复
- 美找到量产石墨烯的简单方法
已经有20人回复
- 【求助】怎么实现单轴拉伸模拟？
已经有11人回复
- 【讨论】材料拉伸测试时单轴向是什么意思？
已经有10人回复
» 抢金币啦！回帖就可以得到:
-
祈福
+5/550
-
中科院招聘博士后及科研助理
+2/228
-
3,4-二(三甲基硅氧基)-1-丁烯（CAS：89543-84-0）
+1/87
-
乙二醇二甲基环硅醚（CAS：57915-66-9）
+1/87
-
1,3,2-Dioxaphospholane, 2-(2,2,3,3,3-pentafluoropropoxy)-
+1/86
-
丁烯二醇（CAS: 497-06-3）
+1/86
-
实验室富余
+1/76
-
新人投稿arxiv 需要cs.IT/CR大佬背书
+1/68
-
坐标重庆 诚征男友
+1/63
-
哈工大（深圳）何自开团队诚招2027级化学工程与技术专业博士生
+1/32
-
哈工大（深圳）何自开团队诚招化学工程与技术专业博士生
+1/28
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
+1/12
-
夏尽秋来，蓄力沉淀，静待硕果成文
+1/11
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
+1/10
-
希望大家都高中！！
+1/7
-
香港理工大学殷骏课题组诚聘1 名博士后以及2 名博士生
+1/5
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/5
-
氨氮空白偏高！求指教！
+1/4
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/3
-
香港中文大学（深圳）管君课题组 纳米光子学方向 招收硕士、博士、博士后
+1/2
2楼2013-08-12 13:51:27
3楼2017-05-10 09:19:30
