Local Stress

Hello, I want to ask how to separate the stress per atom due to bonded and non bonded interaction using lammps. The non bonded interaction used is Lennard jones while the bonded are bond, angle and dihedral. Can you help me?

Did you have a look at this page? It seems from it that if you list bond or angle as keyword to the compute stress/atom command, it will return the corresponding stress. In your case, you just need to define as many compute stress/atom as you need.

1 Like