求助铜薄膜拉伸的问题

[sweetivy]

我将薄膜的左端固定，右端施加一个x方向的位移，进行拉伸。可是最后结果右端没有移动，而且薄膜从左端固定的地方断开。 请帮忙看一下in文件哪有错
误，谢谢。


units metal
boundary p p s

atom_style atomic
lattice fcc 3.61
region box block 0 20 0 20 0 2
create_box 1 box
create_atoms 1 box

pair_style eam
pair_coeff * * Cu_u3.eam

neighbor 1.0 bin
neigh_modify delay 5

region left block INF 1 INF INF INF INF
region right block 19 INF INF INF INF INF
group left region left
group right region right
group boundary union right left


fix 1 all nvt 0.01 0.01 100 drag 0.2

run 10000

velocity boundary set 0.0 0.0 0.0 units box
fix 2 boundary setforce 0.0 0.0 0.0


thermo_style custom step atoms pxx
thermo 500


dump 1 all custom 100 dumpcu.lammpstrj tag type x y z

# for loop
label forloop
variable i loop 200


displace_atoms right move 0.1444 0 0 units box

fix 5 all nve


timestep 0.005


run 1000

next i
jump in.cu forloop

[357921449]

不好意思，我是初学者，不知道问题出在哪！

---

网易邮箱，中国第一大电子邮件服务商