# pushing in mesh
processors	2 2 2

# reading in restart file
read_restart	restart3/platesinkage_test_12000000.restart

neighbor        0.002 bin
neigh_modify    delay 0
communicate     single vel yes
newton          off

# 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 = 0.1 microsecond
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

# fix sphere
fix             integr all nve/sphere

# insert cad
fix             plate all mesh/surface/stress file mesh/plate.stl type 1 rotate axis 1.0 0.0 0.0 angle 90 move -0.04 0.04 0.060 stress on
run             1

# 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
dump            dmpstress all mesh/gran/VTK 20000 post3/plate_*.vtk stress stresscomponents wear plate
restart         2000000 restart3/platesinkage_test_*.restart
thermo_style    one

variable        time equal step*dt

variable        fx equal f_plate[1]
variable        fy equal f_plate[2]
variable        fz equal f_plate[3]
variable        Mx equal f_plate[4]
variable        My equal f_plate[5]
variable        Mz equal f_plate[6]
variable        move_x equal f_plate[7]
variable        move_y equal f_plate[8]
variable        move_z equal f_plate[9]

fix             forceslog all print 20000 "${time},${fx},${fy},${fz},${Mx},${My},${Mz},${move_x},${move_y},${move_z}" append forces/forces.csv title "t,Fx,Fy,Fz,Mx,My,Mz,x,y,z" screen no

# move mesh
fix             move_plate all move/mesh mesh plate linear 0.0 0.0 -0.007

# run
run		14000000 upto
run		16000000 upto
run		18000000 upto
run		20000000 upto
run		22000000 upto
run		24000000 upto
run		26000000 upto
run		28000000 upto
run		30000000 upto