[lammps-users] Problem with the force sum

LAMMPS LAMMPS Mailing List Mirror
1

Dear LAMMPS users,

I met a very strange problem in my result. My single rigid body is constituted with three groups which are adjacent. When I output the ‘fcm’ of the rigid body, the fcm is zero. However, when I sum up three respective ‘fcm’ of each group. The result is not zero.

I think these two should be obviously same. Why not?

Your kind reply would be most appreciated.

Best
Yajie

2

I found the reason, because the force are double counted for the atoms at the adjacent line when I sum up them together. Please discard my question.

Best
Yajie