ERROR: gran model hook/stiffness tangential history

Submitted by teenut567 on Wed, 06/18/2014 - 17:18

Hi everybody,
i have a problem with implementing the "gran model hook/stiffness tangential history" model:
If i replace "hook" with "hertz" curiously the program works fine.

Here's the first part of my input-script:
atom_style granular
boundary f f f
newton off
atom_modify sort 0 0

communicate single vel yes
units si
processors 5 2 1

region reg block -0.2 0.2 -0.1 0.3 -0.006 0.0 units box
create_box 2 reg
region box block -0.2 0.2 -0.0 0.3 -0.006 0.0 units box
region box2 block -0.02 0.02 -0.08 0.0 -0.006 0.0 units box

neighbor 0.006 bin
neigh_modify delay 0

#Material properties required for new pair styles

fix m1 all property/global kn peratomtypepair 2 400e3 600.0e3 600.0e3 600.0e3
fix m2 all property/global kt peratomtypepair 2 400e3 600.0e3 600.0e3 600.0e3
fix m3 all property/global gamman_abs peratomtypepair 2 1.4 2. 4 2.4 2.4
fix m4 all property/global gammat_abs peratomtypepair 2 1.4 2. 4 2.4 2.4
fix m5 all property/global coefficientFriction peratomtypepair 2 0.6 0.4 0.4 0.4

#New pair style
pair_style gran model hook/stiffness tangential history
pair_coeff * *

timestep 5.0e-6

fix 1 all nve/sphere
fix 2 all gravity 9.81 vector 0.0 -1.0 0.0

#import mesh
fix cad1 all mesh/surface file meshes/geometry.stl type 2 #scale 0.01
fix cad2 all mesh/surface file meshes/geometry2.stl type 2 scale 0.01

#use the imported mesh as granular wall
fix granwalls all wall/gran model hook/stiffness tangential history mesh n_meshes 2 meshes geometry geometry2

# Create lid
#fix lid all wall/gran/hook/history primitive type 1 zplane 0.0
fix lid all wall/gran model hook/stiffness tangential history primitive type 2 yplane 0.0
...

And this is the error message of my output-script:

LIGGGHTS (Version LIGGGHTS-PUBLIC 3.0.1, compiled 2014-06-16-17:23:14)
Created orthogonal box = (-0.2 -0.1 -0.006) to (0.2 0.3 0)
2 by 5 by 1 MPI processor grid
ERROR: Invalid pair style (../force.cpp:287)

I hope anybody of you can hely me with my problem.