Error :Could not locate a fix/property storing value(s) for characteristicVelocity as requested by model hooke

Submitted by Rohith tvn on Tue, 04/09/2019 - 13:40

Im doing tutorial and got an error and not able to solve. here is my code can anyone help me with it.
units si
atom_style sphere
atom_modify map array
boundary f f f
newton off
communicate single vel yes
#processors 2 2 3

# Declare domain
region domain block -0.138 0.138 -0.138 0.138 -0.0045 0.43 units box
create_box 2 domain

# Neighbor listing
neighbor 0.003 bin
neigh_modify every 1 check no

# Define the physics
pair_style gran model hooke tangential history
pair_coeff * *

### Setup

# Material and interaction properties required
fix m1 all property/global youngsModulus peratomtype 2.5e7 2.5e7
fix m2 all property/global poissonsRatio peratomtype 0.25 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.5 0.5 0.5 0.5
fix m4 all property/global coefficientFriction peratomtypepair 2 0.2 0.175 0.175 0.5

# Particle insertion
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 1000 &
radius constant 0.0015
fix pdd1 all particledistribution/discrete 15485867 1 pts1 1.0
region factory cylinder z 0 0 0.10 0.4 0.41 units box
fix ins all insert/rate/region seed 86028157 distributiontemplate pdd1 &
nparticles 150000 particlerate 150000 insert_every 1000&
overlapcheck yes vel constant 0. 0. 10 region factory ntry_mc 10000

# Import mesh from cad:
fix cad1 all mesh/surface file cad_files_liggghts_tutorial/binflow/hopper.stl type 2 scale 0.001

# Use the imported mesh as granular wall
fix geometry all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad1

# Create stopper for funnel
fix stopper all wall/gran model hertz tangential history primitive type 2 zplane 0.0

# Detailed settings

# Integrator
fix integrate all nve/sphere

# Gravity
fix grav all gravity 9.81 vector 0.0 0.0 -1.0

# Time step
timestep 0.00001

# Thermodynamic output settings
thermo_style custom step atoms ke cpu
thermo 10000
thermo_modify norm no lost ignore

# Check time set and intialize dump file
fix ctg all check/timestep/gran 1000 0.1 0.1
run 10000
unfix ctg

# Dump output
dump dumpstl all stl 10000 posth/dump*.stl
dump dmp all custom 10000 posth/dump.1 id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius

### Execution and further settings

# Fill the hopper
run 150000 upto

# Pull the stopper
unfix stopper

run 300000

richti83's picture

richti83 | Sat, 04/13/2019 - 15:10

you are missing a fix

fix m5 all property/global characteristicVelocity scalar 2.

see LIGGGHTS-PUBLIC\examples\LIGGGHTS\Tutorials_public\contactModels\in.newModels

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

jonathancuba | Wed, 05/15/2019 - 23:12

hello I have a question, which means this command:
Fix m5 all properties / global characteristic Scalar velocity 2.
It allows to increase the speed twice?