# testing plate sinkage
processors	2 2 2

# variables
variable        d_max equal 0.0115
variable        r_max equal ${d_max}/2
variable        density_const equal 2480

# setting up simulation environment
atom_style	granular
atom_modify     map array
units           si
dimension	3
newton          off
boundary        f f f

communicate     single vel yes

# creating box
region          domain block -0.100 0.100 -0.100 0.100 0.000 0.500
create_box	2 domain

neighbor        0.002 bin
neigh_modify    delay 0

# material properties
fix             m1 all property/global youngsModulus peratomtype 70.e6 70.e6
fix             m2 all property/global poissonsRatio peratomtype 0.24 0.24
fix             m3 all property/global coefficientRestitution peratomtypepair 2 0.9 0.5 0.5 0.5
fix             m4 all property/global coefficientFriction peratomtypepair 2 0.16 0.45 0.45 0.45
fix             m5 all property/global coefficientRollingFriction peratomtypepair 2 0.09 0.09 0.09 0.09

# pair style
pair_style	gran model hertz tangential history rolling_friction cdt torsionTorque on
pair_coeff	* *

# timestep
timestep        0.0000001
fix             ts_check all check/timestep/gran 1000 0.2 0.2

# gravity
fix             gravi all gravity 9.81 vector 0.0 0.0 -1.0

# fix box walls
fix             boxwalls_x1 all wall/gran model hertz tangential history primitive type 2 xplane +0.100
fix             boxwalls_x2 all wall/gran model hertz tangential history primitive type 2 xplane -0.100
fix             boxwalls_y1 all wall/gran model hertz tangential history primitive type 2 yplane +0.100
fix             boxwalls_y2 all wall/gran model hertz tangential history primitive type 2 yplane -0.100
fix             boxwalls_z1 all wall/gran model hertz tangential history primitive type 2 zplane +0.000
fix             boxwalls_z2 all wall/gran model hertz tangential history primitive type 2 zplane +0.500

# particle distribution
fix             pts1 all particletemplate/sphere 91428863 atom_type 1 density constant ${density_const} radius constant ${r_max}
fix             pdd1 all particledistribution/discrete 91428989 1 pts1 1.0

# data dumping
dump            dmp all custom/vtk 20000 post3/platesinkage_test_*.vtk id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
restart         2000000 restart3/platesinkage_test_*.restart
thermo_style    one

# insertion
fix             insert all insert/pack seed 15485863 distributiontemplate pdd1 insert_every once vel constant 0 0 -1.0 overlapcheck yes all_in yes particles_in_region 1500 region domain ntry_mc 1000
fix             integr all nve/sphere

# run
run		2000000 upto

unfix		insert

run		4000000 upto
run		6000000 upto
run		8000000 upto
run		10000000 upto