| 查看: 1371 | 回复: 6 | ||
[求助]
vmd 已有1人参与
|
|
lammps 输出的lammpsstrj导进vmd中模型为什么会出现连续的移动？我是做拉伸模拟 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
帮忙看看fileCode
已经有7人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
alundilong
木虫 (正式写手)
- 应助: 129 (高中生)
- 金币: 2019.1
- 红花: 18
- 帖子: 309
- 在线: 56小时
- 虫号: 595468
- 注册: 2008-09-08
- 性别: GG
- 专业: 传热传质学
2楼2016-04-24 01:25:35
3楼2016-04-24 09:55:55
alundilong
木虫 (正式写手)
- 应助: 129 (高中生)
- 金币: 2019.1
- 红花: 18
- 帖子: 309
- 在线: 56小时
- 虫号: 595468
- 注册: 2008-09-08
- 性别: GG
- 专业: 传热传质学
4楼2016-04-25 04:51:00
|
# LAMMPS cu units metal dimension 3 boundary s s s atom_style atomic neighbor 2.0 bin neigh_modify every 1 read_data data.cu pair_style hybrid eam/alloy airebo 3.0 1 1 lj/cut 10.0 pair_coeff * * eam/alloy Cu_zhou.eam.alloy NULL Cu pair_coeff * * airebo CH.airebo C NULL pair_coeff 1 2 lj/cut 0.019996 3.225 timestep 0.001 thermo 1000 group lower id <= 133 group force id >= 2539 group boundary union lower force group mobile subtract all lower force velocity mobile create 300 825577 dist gaussian velocity boundary set 0.0 0.0 0.0 units box fix 2 lower setforce 0.0 0.0 0.0 fix 3 force setforce 0.0 0.0 0.0 minimize 1.0e-8 1.0e-8 10000 10000 reset_timestep 0 compute 2 mobile temp/com compute 3 all pe/atom compute 4 all ke/atom compute 5 all coord/atom 3.0 compute s1 all stress/atom NULL compute strec1 mobile reduce sum c_s1[2] variable strecy equal c_strec1/(ly*lz*lx*10000) compute strec2 mobile reduce sum c_s1[1] variable strecx equal c_strec2/(ly*lz*lx*10000) thermo_style custom step c_2 pe ke vol v_strecx v_strecy lx ly lz dump 6 all custom 1000 lashen.lammpstrj id type x y z c_3 fix 10 all nvt temp 300 300 0.01 fix_modify 10 temp 2 label loop variable i loop 300 displace_atoms force move 0.0 0.05 0.0 units box displace_atoms lower move 0.0 -0.05 0.0 units box run 100000 next i jump in.cu loop |
5楼2016-04-25 15:05:17
alundilong
木虫 (正式写手)
- 应助: 129 (高中生)
- 金币: 2019.1
- 红花: 18
- 帖子: 309
- 在线: 56小时
- 虫号: 595468
- 注册: 2008-09-08
- 性别: GG
- 专业: 传热传质学
【答案】应助回帖
|
abel loop variable i loop 300 displace_atoms force move 0.0 0.05 0.0 units box displace_atoms lower move 0.0 -0.05 0.0 units box 这几行放在了循环里，那么你的group force/lower的原子，在每个循环中都要平移，不是吗？ 楼主为什么不用fix deform拉升呢？ http://lammps.sandia.gov/doc/fix_deform.html https://icme.hpc.msstate.edu/med ... anowire_Deformation |
» 本帖已获得的红花（最新10朵）
6楼2016-04-26 04:08:12
送红花一朵 |
谢谢，那应该怎么修改一下in文件它才不会发生平移把displace atom这几行移出label循环中？我之前用的是deform拉伸，每次拉伸都不在中间断裂而且得到数值不准确所以我就换了一种拉伸方法 发自小木虫Android客户端 |
7楼2016-04-26 09:17:27
100