Hello:
Everyone, I want to simulate the tension behavior of a macroscopical material.So the periodic boundary conditions are used. The core commands used in the input file are like this:
fix 1 all npt 20 20 0.1 aniso 0.0 0.0 0.0 0.0 NULL NULL 0.1 drag 1.0
variable i loop 80
label loop
print step=$i
displace_box all z delta -0.1 0.1 units box
dump 1 all atom 100 gdc-$i.lammpstrj
run 1000
undump 1
next i
jump in.tension loop
Is is correct? Hoping for suggestions! Thanks!
|
You could do something like that. Or you
could use fix deform to apply the tension more
continuously.
Steve