I am using the second bond to stretch linear chain(s) in a part of the
simulation domain. It would be enough to apply force to the head and
to the tail of the chain but I did not find how to "synchronize" the
forces to conserve linear momentum. They should be simultaneously
applied and released ones certain condition is met.
It should be more clear from the patch I am attaching. It adds
directory USER-STRETCHING, a draft of documentation in
doc/bond_stretching.[txt|html] and simple test cases in
examples/polymer-stretching-test
I also put documentation and animation online
http://slitvinov.github.com/lammps/bond_stretching/
http://slitvinov.github.com/lammps/bond_stretching/GIF/st.gif
0001-Add-bond_style-stretching-in-USER-STRETCHING.-Add-ex.patch (26.8 KB)