| 查看: 2593 | 回复: 1 | ||
zzhwise金虫 (小有名气)
|
[求助]
关于使用setforce施加压力相关问题 已有1人参与
|
|
我使用以下命令，拉动上平板沿着x方向运动，同时给其施加一个沿着z轴负方向的力，为什么它会穿过另一个平板 units real dimension 3 boundary p p p newton on atom_style atomic lattice fcc 4.08 region box block 0 11 0 11 -10 20 create_box 4 box region upboundary block 0 11 0 11 12 INF units lattice region lowboundary block 0 11 0 11 INF -3 units lattice region force block 0 11 0 11 9 12 units lattice region lower block 0 11 0 11 -3 0 units lattice region upper block 0 11 0 11 7 9 units lattice create_atoms 1 region lower create_atoms 1 region force mass 1 196.096655 group lower region lower group upper region upper group force region force group Rigid union force #拉动组 group FREEZ union lower neigh_modify every 1 delay 10 exclude group Rigid Rigid #刚体内相互作用不考虑 pair_style eam pair_coeff * * Au_u3_real.eam neighbor 0.2 bin # rcutoff+0.2 neigh_modify every 1 delay 10 timestep 2 thermo 200 #结构优化 min_style sd minimize 1.0e-4 1.0e-4 10000 1000 timestep 2 # fix 80 all nve velocity upper create 300 12345 mom yes rot no fix 10 upper nvt temp 300.0 300.0 100.0 #加温 fix 20 Rigid rigid group 1 Rigid force 1 off off on #刚体，xy方向无力，z方向不取消力 fix 21 Rigid move linear 0.1 0.0 NULL #右拉 fix 22 Rigid setforce 0.0 0.0 -1.0 #下压 fix 30 FREEZ setforce 0.0 0.0 0.0 #固定不动下版 dump 1 all xyz 500 model1.xyz run 10000  1.png  2.png |
回复此楼» 猜你喜欢
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有12人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
2楼2017-04-27 18:03:58
