24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wuzhihua1988

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 62
帖子: 33
在线: 42小时
虫号: 1717025

[交流] 石墨烯单轴拉伸模拟

这是用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

应力-应变曲线

回复此楼