## Example 5: constant volumn simple shear in a 2D granular material # New tricks: this is very similar to example 4 intruder. # Commands that appears in previous examples will not be repeatedly explained. # More command details in the documentation https://lammps.sandia.gov/doc/Commands.html # LAMMPS version: LAMMPS (7 Aug 2019) on linux dimension 2 atom_style sphere boundary p p p newton off comm_modify vel yes region reg block -5 5 -5 5 -0.5 0.5 units box create_box 3 reg read_data packingConfig add merge read_data boundaryConfig add append group movingBd type == 2 group fixedBd type == 3 group packing type == 1 neighbor 0.2 bin neigh_modify delay 0 exclude group movingBd fixedBd pair_style gran/hooke/history 400.0 NULL 10.0 NULL 0.75 1 pair_coeff * * timestep 0.001 fix 2 all enforce2d fix 1 packing nve/sphere compute ke all ke compute re all erotate/sphere thermo_style custom step atoms c_ke c_re thermo 100 thermo_modify lost warn norm no dump 2 all image 100 *.jpg type diameter zoom 1.6 view 0 0 dump_modify 2 pad 5 run 5000 fix hxFix movingBd move linear 0.7 0.0 NULL run 10000