How to add particles in cylinder disc (closely packing) from Particle size distribution

Submitted by chow2003 on Sun, 06/02/2019 - 20:43

Hello Guys,

I'm new to Liggghts. I created code from Particle size distribution, but I'm not understanding how to make it in group and insert particles randomly in cylinder disc.

Please help me.. thank you

chow2003 | Tue, 06/04/2019 - 16:43

help me to insert particles inside the cylinder plate from data of Particle size distribution and I got Error of 'volume expectancy is too small'

### Initialization

# Preliminaries
units si
atom_style sphere
boundary f f f
newton off
communicate single vel yes

# Domain
region reg block -0.320 0.320 -1.960 1.460 -0.290 0.335 units box
create_box 3 reg
Created orthogonal box = (-0.32 -1.96 -0.29) to (0.32 1.46 0.335)
1 by 1 by 1 MPI processor grid

### Setup

# Material and interaction properties required
fix m1 all property/global youngsModulus peratomtype 2.5e7 2.5e7 2.5e7
fix m2 all property/global poissonsRatio peratomtype 0.25 0.25 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 3 1 0.5 1 0.5 1 1 1 1 1
fix m4 all property/global coefficientFriction peratomtypepair 3 1 0.5 1 0.5 1 1 1 1 1

###Particle Insertion

region bc1 cylinder z 2 3 5 -5.0 1 units box
fix pts1 all particletemplate/sphere 12347 atom_type 1 density constant 1800.0 radius constant 4.7e-5
fix pts2 all particletemplate/sphere 12373 atom_type 1 density constant 1800.0 radius constant 3.45e-5
fix pts3 all particletemplate/sphere 12377 atom_type 1 density constant 1800.0 radius constant 2.55e-5
fix pts4 all particletemplate/sphere 12379 atom_type 1 density constant 1800.0 radius constant 1.9e-5
fix pts5 all particletemplate/sphere 12391 atom_type 1 density constant 1800.0 radius constant 1.4e-5
fix pts6 all particletemplate/sphere 12401 atom_type 1 density constant 1800.0 radius constant 1e-5
fix pts7 all particletemplate/sphere 12409 atom_type 1 density constant 1800.0 radius constant 7.5e-6
fix pts8 all particletemplate/sphere 12413 atom_type 1 density constant 1800.0 radius constant 5.5e-6
fix pts9 all particletemplate/sphere 12421 atom_type 1 density constant 1800.0 radius constant 4e-6
fix pts10 all particletemplate/sphere 12433 atom_type 1 density constant 1800.0 radius constant 3e-6
fix pts11 all particletemplate/sphere 12437 atom_type 1 density constant 1800.0 radius constant 2e-6
fix pts12 all particletemplate/sphere 12451 atom_type 1 density constant 1800.0 radius constant 1.5e-6
fix pts13 all particletemplate/sphere 12457 atom_type 1 density constant 1800.0 radius constant 1e-6
fix pts14 all particletemplate/sphere 12569 atom_type 1 density constant 1800.0 radius constant 5e-7
fix pdd1 all particledistribution/discrete/numberbased 15485959 14 pts1 0.007 pts2 0.045 pts3 0.1004 pts4 0.1322 pts5 0.1244 pts6 0.1007 pts7 0.0794 pts8 0.0655 pts9 0.0608 pts10 0.0601 pts11 0.0556 pts12 0.0449 pts13 0.0397 pts14 0.0843
ERROR: Fix particletemplate/sphere (id pts1): Volume expectancy too small. Change 'volume_limit' if you are sure you know what you're doing (../fix_template_sphere.cpp:436)
numberbased 15485959 14 pts1 0.007 pts2 0.045 pts3 0.1004 pts4 0.1322 pts5 0.1244 pts6 0.1007 pts7 0.0794 pts8 0.0655 pts9 0.0608 pts10 0.0601 pts11 0.0556 pts12 0.0449 pts13 0.0397 pts14 0.0843
ERROR: Fix particletemplate/sphere (id pts1): Volume expectancy too small. Change 'volume_limit' if you are sure you know what you're doing (../fix_template_sphere.cpp:436)