Fix indent force constant

Dear Lammps users,

I’m trying to simulate nanoindentation on Li-Si alloys of various compositions, using fix indent command. In this simulation I need to assume the indenter as an infinitely rigid body. My question is, how to determine the force constant (K) for this simulation? Should it vary for different materials?



That’s not really a LAMMPS Q. The K is an input

to fix indent. It’s up to you to choose a value

good for your model. The fact that the indenter

is rigid is a separate issue from the value of K.

Fix indent is modeling a rigid sphere, etc.