How to fix a particle in a simulation and how to output damping force?

Submitted by jerry_1988 on Wed, 05/18/2011 - 19:46

Hi,

I may modify the LIGGGHTS code for a new contact model and do some basic benchmark tests. May anybody tell me how to fix the position and rotation of a particle in a simulation (its force and interaction with other particles remained ), and how to output the damping force of a particle?

Regards
Jerry

Forums:

JF's picture

JF | Wed, 05/18/2011 - 20:40

Hi Jerry,

To fix the position and rotation of a particle in a simulation, you can use: fix ID group-ID freeze (see the LIGGGHTS documentation).

Concerning the damping force of a particle, you can only output always with the fix command the force (Fx,Fy,Fz) of a particle.

Regards
JF

ckloss's picture

ckloss | Thu, 05/19/2011 - 12:59

Regarding damping force:
You can output per-particle total force straight-forward, and per-contact forces with a compute pair/gran/local.
But if you want to output spring/damper force separately per-contact, you would have to code this

Christoph