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 | Fri, 06/06/2014 - 16:31
I hopefully have the chance
I hopefully have the chance to check this before the 3.0.2 release and let you know in that case!
Christoph
rberger | Sat, 06/07/2014 - 20:00
typo
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.