Dear Professor,
My name is Hank, and I am a user of LAMMPS.
Would you mind if I ask you one question about LAMMPS.
My question is that I hope to increase the maximum of touching neighbors for sphere particles but I don’t know how to do that. I have looked up the manual and searched on the Internet. They all say that I can modify the definition of “MAXTOUCH” in the file “fix_shear_history.cpp”. But I can’t find it in my fix_shear_history.cpp. Could you please teach me how to do that?
The attachment is the file fix_shear_history.cpp in the version of my LAMMPS.
And the script below is my code:
units cgs
atom_style sphere
boundary f f f
newton off
communicate single vel yes
lattice fcc 1.0
region boxin1 cylinder z 0 0 5 0 250 units lattice
create_box 1 boxin1
create_atoms 1 random 10000 5 boxin1 units box
group bal region boxin1
set region boxin1 diameter 1
neighbor 0.2 bin
neigh_modify delay 0
pair_style gran/hooke/history 200000.0 NULL 100.0 NULL 0.5 1
timestep 0.0001
fix wigz bal wall/gran 2000000.0 NULL 10000.0 NULL 0.5 1 zcylinder 5
fix wigx bal wall/gran 2000000.0 NULL 10000.0 NULL 0.5 1 zplane 0 250
fix 1 bal nve/sphere
fix g bal gravity 1000.0 spherical 0.0 -180.0
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol
thermo 100
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
dump id all custom 500 dump.10k_1 x y z
run 100000
With best wishes,
Hank
fix_shear_history.cpp (13.3 KB)