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 | Sat, 04/13/2019 - 15:10
you are missing a fix
you are missing a fix
fix m5 all property/global characteristicVelocity scalar 2.
see LIGGGHTS-PUBLIC\examples\LIGGGHTS\Tutorials_public\contactModels\in.newModels
jonathancuba | Wed, 05/15/2019 - 23:12
characteristic Scalar velocity 2 ?
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?
richti83 | Thu, 05/16/2019 - 08:22
https://www.cfdem.com/comment
https://www.cfdem.com/comment/19783#comment-19783