# Preliminaries

units si
atom_style granular
atom_modify map array
boundary m m m
communicate single vel yes
newton off #default, used for multicore systems

# Declare domain
variable xmin equal -5
variable xmax equal 5

variable ymin equal -5
variable ymax equal 5

variable zmin equal -5
variable zmax equal 5
region reg block ${xmin} ${xmax} ${ymin} ${ymax} ${zmin} ${zmax} units box
create_box 1 reg
neighbor 0.002 bin #default
neigh_modify delay 0 #default

 # Import Mesh
 fix deck all mesh/surface file Mesh/Deck.stl type 1 scale 0.001
 fix walls all wall/gran model hertz tangential history mesh n_meshes 1 meshes deck

# Material properties

fix m1 all property/global youngsModulus peratomtype 2.e7
fix m2 all property/global poissonsRatio peratomtype 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.2
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5

# Contact Physics

pair_style  gran model hertz tangential history #Hertzian without cohesion
pair_coeff  * *

# Integrator

fix integrate all nve/sphere

# Gravity
 fix grav all gravity 9.81 vector 0. 1. 0.
 

 # Particle insertion
 
fix pts all particletemplate/sphere 15485863 atom_type 1 density constant 1000 radius constant 0.02
fix pdd all particledistribution/discrete 15485863 1 pts 1.0
fix factory all mesh/surface file Mesh/Inseration.stl type 1 scale 0.01


fix ins all insert/stream seed 86028121 distributiontemplate pdd mass 1000 massrate 1000 overlapcheck yes all_in yes vel constant 0 0.75 0 insertion_face factory extrude_length 0.1
  # Timestep
 timestep 0.001
 
 run 10000