Hello,
I have been studying friction simulation and encounter some issues.
To aid understanding, I have attached a figure.
When applying a normal load to the top layer (as shown in the figure), I used the command "Fix nve". However, in this case, the structure of the top part did not stabilize.
I am wondering how to apply force in the normal direction in a friction simulation.
Below I briefly wrote down my LAMMPS input commands.
#############################
# "lofix" means bottom fix layer
# "hifix" means top support layer
# "lotherm" means bottom temperature control layer
# group
wall union lofix hifix
notwall subtract all wall
# minimization
min_style cg
minimization 1e-15 1e-15 5000 10000
velocity all create 300.0
# Equilibrium
velocity lofix set 0.0 0.0 0.0
fix 1 lofix setforce 0.0 0.0 0.0
fix 2 all npt temp 300.0 300.0 0.05 x 1.0 1.0 0.5 y 1.0 1.0 0.5
run
unfix 2
fix 2 all nvt temp 300.0 300.0 0.05
run
unfix 2
# Apply force
fix 1 all nve
fix 2 hifix spring/self 10
fix 3 hifix addforce 0.0 0.0 -10.0
run
unfix 1
# Move
fix 1 all nve
fix 4 hifix move linear 0.1 0.0 0.0 units box
fix 5 lotherm langevin 300.0 300.0 0.05 1234 tally yes
run
#############################
- Badges
- Science topic