###################################Header for General commands#######################################################

atom_style	granular		#Simulation of particles

boundary	f f f			#fixed boundaries -> particles will be deleted if leaving the simulation box

units		si

communicate	single vel yes	#default

newton		off				#default








































###################################System variables##################################################################

#Definition of boundaries
variable xmin equal -0.00005
variable xmax equal 0.001

variable ymin equal -0.00005
variable ymax equal 0.001

variable zmin equal -0.00005
variable zmax equal 0.002

#Definition of the timestep

variable dt equal 1e-5		#timestep = 0.0001 second; Each iteration step represents 0.0001 seconds.


































###################################Specific variables for current simulation#########################################

variable natoms equal 2 	#1 -> particle #2->Die 

####variable for material properties####

####Young Modulus####
variable	youngmodulus1	equal 	1e8		#N/mm²
variable	youngmodulus2	equal 	1e8		#N/mm²

####Poission ratio####
variable 	poission1	equal	0.3
variable 	poission2	equal	0.3

####variable for contact properties####

####coefficient of restitution####
variable	CoR11	equal	0.6
variable	CoR12	equal	0.6	
variable	CoR21	equal	0.6
variable	CoR22	equal	0.0

####sliding friction coefficient####
variable	sf11	equal	0.3
variable	sf12	equal	0.36
variable	sf21	equal	0.36
variable	sf22	equal	0.0

####rolling friction coefficient####
variable	rf11	equal	0.8
variable	rf12	equal	0.8
variable	rf21	equal	0.8
variable	rf22	equal	0.0





####variable for particle####

#Number of particle radius
variable	nradii	equal	2

variable	radius1	equal	0.000035	#m
variable	radius2	equal	0.00007	#m

variable	frac1	equal	0.4		#40%
variable	frac2	equal	0.6		#60%

variable	density	equal	4000	#kg/m³









































####filling parameters####

variable	filltime	equal	5		#seconds #filling time# 

variable	fillmass	equal	36e-7	#kg #total mass#

variable	fillmassrate	equal	${fillmass}/${filltime}	#kg/s

variable	fillsteps	equal	${filltime}/${dt}		#Transform time to iteration steps

####settle time####

variable	settletime	equal	1 #second

variable	settlesteps	equal	${settletime}/${dt}		#Transform time to iteration steps













###################################Definition of simulationbox#######################################################

region		reg block ${xmin} ${xmax} ${ymin} ${ymax} ${zmin} ${zmax} units box

create_box	2 reg

neighbor	0.004 bin		#default

neigh_modify	delay 0		#default

#################################Definition of the contact models####################################################

pair_style gran model hertz tangential history rolling_friction epsd2	#contact model

pair_coeff	* *				#default

timestep ${dt}

fix		integrator all nve/sphere			#default

fix		gravi all gravity 9.81 vector 0.0 0.0 -1.0		#gravity of 9.81 m/s² in negative z direction














###################################Definition of Material properties#################################################


fix 		m1 all property/global youngsModulus peratomtype ${youngmodulus1} ${youngmodulus2}

fix 		m2 all property/global poissonsRatio peratomtype ${poission1} ${poission2}

fix 		m3 all property/global coefficientRestitution peratomtypepair ${natoms} 0.6 0.6 0.6 0.6

fix 		m4 all property/global coefficientFriction peratomtypepair ${natoms} ${sf11} ${sf12} ${sf21} ${sf22}

fix 		m5 all property/global coefficientRollingFriction peratomtypepair ${natoms} ${rf11} ${rf12} ${rf21} ${rf22}



















###################################Generation and Loading of the Geometry .stl#######################################

fix			die 	all mesh/surface file Die.stl 	type 2  scale 0.001		#load mesh from STL file. Type 2 for geometry. Scale down to transform mm to meters


fix 		walls 	all wall/gran model hertz tangential history rolling_friction epsd2 mesh n_meshes 1 meshes die








































###################################Generation and Insertion of the particles#########################################

fix		pts1 all particletemplate/sphere 10487 atom_type 1 density constant ${density} radius constant ${radius1} volume_limit 1e-20	

fix		pts2 all particletemplate/sphere 11887 atom_type 1 density constant ${density} radius constant ${radius2} volume_limit 1e-20	

fix 	pdd1 all particledistribution/discrete 32452867 ${nradii} pts1 ${frac1} pts2 ${frac2} 

fix		ins_mesh all mesh/surface/planar file Insertionsface.stl type 1 scale 0.001 

fix 	ins all insert/stream seed 86028157 distributiontemplate pdd1 &
        mass ${fillmass} massrate ${fillmassrate} overlapcheck yes vel constant 0 0 -0.75 &
        insertion_face ins_mesh extrude_length 0.0009


		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		
		







###################################Dumping of the data for post-processing to visualize############################

shell 	mkdir post
shell test -w post || echo "Warning: Cannot write to post directory"

#Definition of the dumptime

variable dumptime equal 0.005	# Every 0.05 s 1 image 

variable dumpstep equal ${dumptime}/${dt}	#Transform to iteration steps


dump	dmpparticle all custom ${dumpstep} post/particles_*.liggghts id type x y z vx vy vz fx fy fz radius mass
dump	dmpdie all mesh/stl ${dumpstep} post/Die*.stl die














































####################################RUN the simulation filling###########################################################

run	 ${fillsteps}

unfix	ins

##################################RUN the simulation settling#############################################################

run	 ${settlesteps}