| 查看: 5393 | 回复: 3 | ||
[求助]
LAMMPS纳米压痕模拟
|
|
列位大神，我是分子动力学模拟小白。最近想用LAMMPS模拟纳米压痕过程，把LAMMPS自带的in.indent文件修改了一下到三维情况，以下是我的输入文件： # 2d indenter simulation dimension 3 boundary p s p atom_style atomic neighbor 5.6 bin neigh_modify delay 5 # create geometry lattice fcc 1 #定义一个region，名为box，类型为block,后面的边界坐标 region box block 0 5 0 5 -2.5 2.5 #在region之上创建一个模拟盒子，2为模拟中用到的分子种类， #box为上面定义的region的名字 create_box 2 box create_atoms 1 box mass 1 1.0 mass 2 1.0 # LJ potentials pair_style lj/cut 2.5 pair_coeff * * 1.0 1.0 2.5 # define groups region 1 block INF INF INF 0.5 INF INF group lower region 1 group mobile subtract all lower set group lower type 2 # initial velocities compute new mobile temp velocity mobile create 0.2 482748 temp new fix 1 all nve fix 2 lower setforce 0.0 0.0 0.0 fix 3 all temp/rescale 100 0.1 0.1 0.01 1.0 # run with indenter timestep 0.003 variable k equal 1000.0/xlat variable y equal "-13.0*ylat + step*dt*0.02*ylat" fix 4 all indent $k sphere 2.5 v_y 0 0.5 thermo 1000 thermo_modify temp new dump 1 all atom 250 dump.indent dump 2 all image 1000 image.*.jpg type type & zoom 1.6 adiam 1.5 dump_modify 2 pad 5 dump 3 all movie 1000 movie.mpg type type & zoom 1.6 adiam 1.5 dump_modify 3 pad 5 run 30000 # run without indenter unfix 4 run 30000 可以运行，但是结果不正确。计算一开始完全看不出有压下过程，然后原子的顺序就突然变乱了。。。如下图： 求教！！！  10000.jpg  30000.jpg  31000.jpg |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有76人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2017-04-29 16:33:16
3楼2017-05-11 20:44:25
4楼2018-01-23 09:37:41
50