Dear lammps users,
Now I’m trying to simulate the movement of a carbon nanotube in fluid. I want to treat the CNT as a rigid body using the “fix rigid” command. Also, I want the velocity of the carbon atoms vz (velocity along z direction) to be zero. This means the CNT will only translate and rotate on the x-y plane.
My code is like:
fix 1 cnt rigid single force * on on off torque * on on off
But vz is still not zero.
Also, I find “fix setforce” cannot be used along with “fix rigid”. Actually, if I delete the “fix rigid” command and use “velocity set” and “fix setforce”, vz of the carbon atoms can be set to zero.
I’m looking forward to having your help.