about applying the radial load

Dear all,
I am running an MD simulation with lammps. The cross section of my model is a ring. I want to apply radial load on it as shown in the figure. At first I tried to use the setforce command. But the number of group is limited. I wonder if there are any other commands can be used to achieve the same effect.
360反馈意见截图16171106243652.jpg
May I have your suggestions ? Many thanks for your reply .

Fix setforce allows you to use atom-style variables

for the x,y,z components of force. So you can write 3 variables

which use the coords of an atom to generate a radial

force.

Steve

360反馈意见截图16171106243652.jpg

Dear all,
I am running an MD simulation with lammps. The cross section of my model is a ring. I want to apply radial load on it as shown in the figure. At first I tried to use the setforce command. But the number of group is limited. I wonder if there are any other commands can be used to achieve the same effect.
May I have your suggestions ? Many thanks for your reply .

You should be able to use fix setforce with an atom style variable. No groups needed.

Axel