ERROR MODEL GRAN/HERTZ/HISTORY/STIFFNESS - LIGGGHTS 3.0.0

Submitted by franciscohgt on Thu, 06/05/2014 - 20:48

Hi everyone

I am implementing the great model / hertz / history / stiffness in LIGGGHTS 3.0.0, the code is as follows:

units si
atom_style sphere
atom_modify map array
boundary m f m
newton off
communicate single vel yes
processors 2 2 2

region reg block -2.7 2.7 -0.1 1.1 -2.7 2.7 units box
create_box 1 reg

neighbor 0.001 bin
neigh_modify delay 0

#Material properties required for new pair styles

fix m1 all property/global kn peratomtype 1 1.0e5
fix m2 all property/global kt peratomtype 1 1.0e5
fix m3 all property/global gamman peratomtype 1 0.3
fix m4 all property/global gammat peratomtype 1 0.3
fix m5 all property/global coefficientRestitution peratomtypepair 1 0.38
fix m6 all property/global coefficientFriction peratomtypepair 1 0.3

#New pair style
pair_style gran model hertz stiffness tangential history

pair_coeff * *

timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0

#geometry
fix cad1 all mesh/surface file cilindro.stl type 1
fix cad3 all mesh/surface file plano0.stl type 1
fix cad4 all mesh/surface file plano1.stl type 1

fix granwalls all wall/gran model hertz stiffness tangential history mesh n_meshes 3 meshes cad1 cad3 cad4

fix pts all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.05
fix pdd all particledistribution/discrete 1. 1 pts 1.0

region bc cylinder y 0. 0. 2. 0. 1.

fix ins all insert/pack seed 1001 distributiontemplate pdd insert_every once overlapcheck yes all_in yes volumefraction_region 0.1530694 region bc

fix integrate all nve/sphere

compute rke all erotate/sphere
thermo_style custom step atoms ke c_rke vol
thermo 200
thermo_modify norm no lost ignore
compute_modify thermo_temp dynamic yes

fix ctg all check/timestep/gran 1 0.1 0.1
run 1
unfix ctg

dump dumpstl all mesh/stl 200 post/dump*.stl cad1
dump dmp all custom 200 post/dump*.cilindro id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius

run 100000

----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

And I have the following error:

LIGGGHTS (Version LIGGGHTS-PUBLIC 3.0.0, compiled 2014-05-30-02:38:19 by francisco based on LAMMPS 23 Nov 2013)

units si
atom_style sphere
atom_modify map array
boundary m f m
newton off
communicate single vel yes
processors 2 2 2

region reg block -2.7 2.7 -0.1 1.1 -2.7 2.7 units box
create_box 1 reg
Created orthogonal box = (-2.7 -0.1 -2.7) to (2.7 1.1 2.7)
2 by 2 by 2 MPI processor grid

neighbor 0.001 bin
neigh_modify delay 0

#Material properties required for new pair styles

fix m1 all property/global kn peratomtype 1 1.0e5
fix m2 all property/global kt peratomtype 1 1.0e5
fix m3 all property/global gamman peratomtype 1 0.3
fix m4 all property/global gammat peratomtype 1 0.3
fix m5 all property/global coefficientRestitution peratomtypepair 1 0.38
fix m6 all property/global coefficientFriction peratomtypepair 1 0.3

#New pair style
pair_style gran model hertz stiffness tangential history

pair_coeff * *

timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0

#geometria
fix cad1 all mesh/surface file cilindro.stl type 1
#fix cad2 all mesh/surface file lifter.stl type 1
fix cad3 all mesh/surface file plano0.stl type 1
fix cad4 all mesh/surface file plano1.stl type 1

fix granwalls all wall/gran model hertz stiffness tangential history mesh n_meshes 3 meshes cad1 cad3 cad4
ERROR: Fix wall/gran (id granwalls): Need to use define style 'mesh' or 'primitive' (../fix_wall_gran.cpp:242)

can someone help please?

thanks

ckloss's picture

ckloss | Fri, 06/06/2014 - 16:31

I hopefully have the chance to check this before the 3.0.2 release and let you know in that case!

Christoph

rberger's picture

rberger | Sat, 06/07/2014 - 20:00

First of all, there are a few typos:

pair_style gran model hertz stiffness tangential history

should be

pair_style gran model hertz/stiffness tangential history

and

fix granwalls all wall/gran model hertz stiffness tangential history mesh n_meshes 3 meshes cad1 cad3 cad4

should be

fix granwalls all wall/gran model hertz/stiffness tangential history mesh n_meshes 3 meshes cad1 cad3 cad4

But there should be proper error messages... as Christoph said, we'll put it on our list for 3.0.2.