Dear All,

The Doc. page says lammps calculates the force as F= -K(r-R)^2. I am looking forward to know why there is power 2 for r-R term. I would appreciate it if someone could help me.

Best Regards,

Carmelo

Dear All,

The Doc. page says lammps calculates the force as F= -K(r-R)^2. I am looking forward to know why there is power 2 for r-R term. I would appreciate it if someone could help me.

Best Regards,

Carmelo

Dear All,

The Doc. page says lammps calculates the force as F= -K(r-R)^2. I am

looking forward to know why there is power 2 for r-R term. I would

appreciate it if someone could help me.

it is there because that is how the fix indent force is computed.

please see the source code, if you have doubts.

axel.

If you wanted a different force between

the indenter and particles, you would need

to derive a new fix from fix indent and implement

whatever formula you wish.

Steve