# used for offical lammps, developed by deyun
#========================================================================================
variable skindistance equal 0.0001 #say 10% of max particles
variable interval_fix equal 10000
#========================================================================================
variable clumpratio equal 0.199327526 # be careful, this is estimated from matlab, should be case by case
#========================================================================================
units		si # (This means m, kg, s for base units.)
boundary	p p p 
atom_style	bpm/sphere
dimension	3
newton		on off
echo            both
comm_modify     vel yes 
processors      * * * 
#========================================================================================
#read_data       rubber_sand100small12.data
read_restart   restart_filecheck.137800000
#========================================================================================
neighbor	${skindistance} bin 
group rubber type 1
group sand   type 2
neigh_modify exclude molecule/intra all  delay 0 every 1 one 200 page 2000
bond_style      bpm/rotational 
bond_coeff      1 1e4 1e4 5e5 5e5 1.0e12 1.0e12 1.0e12 1.0e12 1 1 0.0005 0.0005 
special_bonds   lj 0 1 1 coul 1 1 1
#========================================================================================
pair_style granular 
pair_coeff * * hertz/material 7e6 0.5 0.5 tangential mindlin_rescale/force NULL 1.0 0.5
timestep	1e-8
run_style	verlet
#========================================================================================  
compute     7 all pressure NULL pair 
#========================================================================================
variable 	step equal step
variable	xlength equal lx
variable	ylength equal ly
variable	zlength equal lz
variable    xxstress equal c_7[1]
variable    yystress equal c_7[2]
variable    zzstress equal c_7[3]
variable    xystress equal c_7[4]
variable    xzstress equal c_7[5]
variable    yzstress equal c_7[6]
#========================================================================================
fix         9 all press/berendsen x 30000 30000 10 y 30000 30000 10 z 30000 30000 10
fix	        3 all nve/bpm/sphere
#========================================================================================
fix 		2 all print ${interval_fix} "${step} ${xxstress} ${yystress} ${zzstress} ${xystress} ${xzstress} ${yzstress} ${xlength} ${ylength} ${zlength}" file check_mean_stresses_iso.txt screen no
#========================================================================================
restart     100000 restart_filecheck
#=========================================================================================
run		    999999999
#=========================================================================================