Dear Sir,
I am new LIGGGHT user and i am trying to do a simulation in which we will pour the particles of different sizes but I am facing the following problems:
1. I can only see one type of particle in simulation whereas in input script i have used two different particles of different radius.
2.During running the input file the following warning is coming
"Less insertions than requested, the particle size distribution you wish may not be pictured".
3.The particles shown in simulation appears to float, but i want to apply some kind of force on them to compact them, so how can I modify my input script to include this compactive force.
The mass of sample is 4.581 kg containing two different particles of different sizes is to be poured into a cylinder of radius 7cm and 30 cm height.
My input script is:
#Particle inseration
atom_style granular
atom_modify map array
boundary m m m
newton off
echo both
communicate single vel yes
units si
region reg cylinder z 0.0 0.0 0.07 0.0 0.30 units box
create_box 1 reg
neighbor 0.0002 bin
neigh_modify delay 0
#Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 7.e10
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
#New pair style
pair_style gran/hertz/history 1 0 #0 für keine Kohesion
pair_coeff * *
timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
fix zwalls all wall/gran/hertz/history 1 0 zplane 0.0 0.3 1
fix cylwalls all wall/gran/hertz/history 1 0 zcylinder 0.07 1
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.00575
fix pts2 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.009
fix pdd1 all particledistribution/discrete 1. 2 pts1 0.29 pts2 0.71
#region and insertion
group nve_group region reg
region bc cylinder z 0.0 0.0 .07 0.20 .30 units box
fix ins nve_group pour/dev mass 4.581 20 distributiontemplate pdd1 vol 0.5 200 massflowrate 0 vel uniform 0. 0. 0. 0. 0.0 region bc
#apply nve integration to all particles that are inserted as single particles
fix integr nve_group nve/sphere
#insert the first particles
run 1
dump dmp all custom 350 post/dump.test id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
run 100000 upto
unfix ins