Hello,
How can I enable torsion torque for a primitive wall? Please see the error below. The ordering works for pair_style, but not for wall/gran. Thanks!
fix wall_z1 all wall/gran model hertz tangential history cohesion sjkr rolling_friction epsd torsionTorque on primitive type 1 zplane 0 # floor
ERROR: Fix wall/gran (id wall_z1): Unknown argument or wrong keyword order: 'primitive' (../fix_wall_gran_base.h:131)
% script, up to error
5 ## Preliminaries
6 units si
7 atom_style granular
8
9 # the product of the processors must equal P, the number of processors used when running the file: mpirun -np 1 ...
10 processors 1 1 1
11
12 # boundary conditions
13 boundary p f fm
14
15 #LIGGGHTS-Public shouuld produce the same answers for any newton flag setting, besides round off err
16 newton off
17
18 communicate single vel yes
19
20 soft_particles yes
21
22 ## Declare domain
23 region reg block -2.0 2.0 -1.0 10.0 0.0 20.0 units box
24
25 create_box 1 reg
26
27 neighbor 0.0625 bin
28
29 neigh_modify delay 0
30
31 ## Material Property variables
32 variable restiv equal 0.8
33 variable static_fric equal 0.3
34 variable rolling_fric equal 0.8
35 variable visc_damp equal 0.2
36 variable k_c equal 0
37 variable r_p equal 0.00125
38 ## Material properties
39 fix m1 all property/global youngsModulus peratomtype 5e6
40 fix m2 all property/global poissonsRatio peratomtype 0.3
41 fix m3 all property/global coefficientRestitution peratomtypepair 1 ${restiv}
42 fix m4 all property/global coefficientFriction peratomtypepair 1 ${static_fric}
43 fix m5 all property/global characteristicVelocity scalar .01
44 fix m6 all property/global cohesionEnergyDensity peratomtypepair 1 ${k_c}
45 fix m7 all property/global coefficientRollingFriction peratomtypepair 1 ${rolling_fric}
46 fix m8 all property/global coefficientRollingViscousDamping peratomtypepair 1 ${visc_damp}
47
48 ## Pair style
49 pair_style gran model hertz tangential history cohesion sjkr rolling_friction epsd torsionTorque on
50 pair_coeff * *
51
52 # timestep
53 timestep 0.00001
54 # integration
55 fix 1 all nve/sphere
56
57 # Earth gravity
58 fix 2 all gravity 9.81 vector 0.0 0.0 -1.0
59
60 ## Thermo settings
61 thermo_style custom step atoms ke vol
62 # thermo output
63 thermo 1000
64 thermo_modify lost ignore norm no
65
66 # floor
67 fix wall_z1 all wall/gran model hertz tangential history cohesion sjkr rolling_friction epsd torsionTorque on primitive type 1 zplane 0