LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] The force in fix indent
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] The force in fix indent


From: Steve Plimpton <sjplimp@...24...>
Date: Tue, 4 Jul 2017 18:29:41 -0600

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

On Tue, Jul 4, 2017 at 9:28 AM, Axel Kohlmeyer <akohlmey@...24...> wrote:
On Tue, Jul 4, 2017 at 7:02 AM, 万 京 <luckwan16@...4215...> wrote:
> 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.


>
>
> Best Regards,
>
>
> Carmelo
>
>
>
>
> ------------------------------------------------------------------------------
> Check out the vibrant tech community on one of the world's most
> engaging tech sites, Slashdot.org! http://sdm.link/slashdot
> _______________________________________________
> lammps-users mailing list
> lammps-users@...12...396...sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/lammps-users
>



--
Dr. Axel Kohlmeyer  akohlmey@...43...4...  http://goo.gl/1wk0
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.

------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
lammps-users mailing list
lammps-users@...6297....sourceforge.net
https://lists.sourceforge.net/lists/listinfo/lammps-users