[lammps-users] What's wrong with the result of the tension behavior?


Hello Everyone:
A tension behavior is simulated with lammps, the input file is next:
relax******
fix 1 all npt 20 20 0.1 aniso 0.0 0.0 0.0 0.0 0.0 0.0 0.1 drag 1.0
timestep 0.001
run 10000
undump 1
unfix 1
tension*
fix 1 all npt 20 20 0.1 aniso 0.0 0.0 0.0 0.0 0.0 0.0 0.1 drag 1.0
variable i loop 15
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 500
undump 1
next i
jump in.tension loop

The stress-strain curve is in the annex. Generally ,when the strain is zero, the stress should be zero too. But in the curve, the case is not so, the stress shows a vibration at the beginning. Who can tell me what’s wrong with the curve? Thanks !
|

1.bmp (307 KB)

Hi Chen,

I’m not quite sure why the curve behaves so, but
when you do tension along the z-axis with fix npt,
you probably need to turn off stress control along z-axis
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

Hope this help.

L. Wan

2009/9/2 云俊 陈 <yunchen1986@…215…>