I’m attempting to simulate the bending of a nanobeam using the indent feature in LAMMPS. I followed the documentation for the fix indent command, but for some reason, it doesn’t seem to have any effect on the beam. Even after running the simulation for a considerable number of timesteps, there’s no observable change. Additionally, upon inspection, I noticed that all components of the indent vector are zero. Below is the code I’m using. I’m hoping to receive some advice on how to address this issue. Thank you.
indent.ln (3.1 KB)
Many commands in LAMMPS have a “units” keyword. This is not a problem when not using a lattice command, but when setting a lattice command, this makes a difference and beginner users often overlook this.
Thanks for your answer. I am not sure if I understand your tip correctly, but I have added “units box” to the “fix 3 middle indent 1000.0 sphere 29.8705 v_y 581.117 13.5775 units box”. Again, nothing changed and the problem still exists. Or I may misunderstood your tip.
My tip is correct, but that is not the only problem with your input.
Another is just right in front of you. Unlike the more subtle one that I reminded out of, this is so obvious, that you’ll have to be able to find yourself. Just look at the thermo output.