Fw:Re: ??? Addforce using the coordination of atom

Dear Shafat Mubin

Thanks for your remind!

I will see the fix_rigid.cpp later.But from the description,“The fix rigid command treats a collection of particles as a rigid body, computes its inertia tensor, sums the total force and torque on the rigid body each timestep due to forces on its constituent particles, and integrates the motion of the rigid body.”(6.14 from LAMMPS Users Manual).The time integration should consider the force correctly.If not,there may be a bug.

Using the dump command,I just only input the force about the single particle.So how to judge the force added to rigid bodies is wrong?Is there some way to input the force added to rigid bodies?

Best Wishes!

Chen Wang