Vol.21 No.4 程远方等：球形颗粒随机排列过程的计算机模拟 ·391· 2 Finney JL.Random Packings and the Structure of Simple 9 Cundall PA.Formulation of a Three-dimensional Distinct Liquids.Proc R Soc,1970,319A:479 Element Model?Part I.Int J Rock Mech,Min Sci Ge- 3 Bordia R K.A Theoretical Analysis of Random Packing omech Abstr,1988,25(3):107 Densities of Mono-sized Spheres in Two and Three Di- 10 Hart R,Cundall P A,Lemos J.Formulation of a Three-di- mensions.Script Metallurgica,1984,18:725 mensional Distinct Element Model?Part II.Int J Rock 4 Sutherland D N.Random Packing of Circles in Plane.J Mech,Min Sci Geomech Abstr,1988,25(3):117 Coll Interface Sci,1977,60:96 11 Lange F F.Effect of Interparticle Potentials on Particle 5 German R M.Particle Packing Characteristics.New Jersey: Packing for Ceramic Processing,Powder &Grains 93. Princeton,1989 Thorm T,ed.Balkema:Rotterdam,1998.187 6 Chen Jianhong.Computer Simulation of 2D Particle Pack- 12 Cundall P A.Distinct Element Models of Rock and Soil ing,Sintering Processes and Microstructure-Properties Re- Structure.In:Analytical and Computational Methods in lations:[PHD Thesis].Alfred:Alfred University,1990 Engineering Rock Mechanics,Brown E T,ed.London: 7 Zheng Jinmin.Computer Simulation of Particle Packing Allen Unwin,1987.129 and Sintering:[PhD Thesis].Alfred:Alfred University, 13 Bouvard D,Lange F F.Relation between Percolation and 1991 Particle Coordination in Binary Powder Mixtures.Acta 8 Cundall P A,Strack O D C.A Distinct Numerical Model metall Mater,1991,39(12):3083 for Granular Assemblies.Geotechnique,1979,29:47 Computer Simulation of Random Packing of Spherical Particles Cheng Yunafang,Guo Shijiun,Lai Heyi Material Science and Engineering School,UST Beijing,Beijing 100083,China ABSTRACT A computer simulation for evolution of random structures of spherical particles in two dimen- sions has been performed by the distinct element method.The forces among the particles considered were gravity,friction force and Van Der Waals Interaction(VDWI).The motion of particles involved translation and rotation.For 100 um particles,if friction and VDWI were not considered,the packing assembly was a mixture of random packing and order packing with packing density of 0.8696,if friction and VDWI were taken into account,the packing density was 0.821 3.To study the procedure of packing,packing densities with different time,particle trajectory and coordination number were examined.The packing structures of binary system with different conditions were also given. KEYWORDS particle random packing;distinct element method;coordination number;packing density从, 1 . 2 1 N o . 4 程 远方 等 : 球 形颗粒 随机 排 列过 程的计 算机 模拟 · 3 9 1 · 2 Fi n n e y J L . R a n d o m P a e k i n g s an d ht e St ur e utr e o f S im p l e L i qu id s . P ocr R S o e , 19 7 0 , 3 19A : 4 7 9 3 B o dr i a R K . A T h e o er ti c a l nA a ly s i s o f R a n d o m P a e k i n g D e n s i t i e s o f M o n o 一 s i z e d SPh e r e s i n l切 0 an d T hr e e D i - m e n s i o n s . S e ir P t M e at ll u rg l e a , 19 84 , 1 8 : 7 2 5 4 S u ht e r lan d D N . R an d o m P a e ki n g o f C i r e l e s i n P lan e . J C o ll I n t e far e e S e i , 19 7 7 , 6 0 : 9 6 5 G e mr an R M . P art i e l e P a e k i n g C har a e t e r i s t i e s . N e w J e r s e y : P r i n e 以 o n , 1 9 89 6 C h e n J 1 a n h o n g . C o m P ut e r Si m u l at i o n o f 2D P art i e l e P a e k - i n g , Si n t e ir n g P r o e e s s e s an d M i e r o s trU e ut r e 一 P r o P e rt i e s R e - lat i o ns : [ P l刃D T h e s i s ] . A l fer d : A l fer d U n i v e r s iyt , 19 90 7 Z h e n g l i钊m i n . C om P ut e r S而u lat i o n o f p art i e l e p a e k i n g 叨d Si n t e ir gn : [ P hD hT e s i s ] . A l fre d : A l fre d U n i v e r s i ty, 19 9 1 8 C un d a ll P A , S tr a c k 0 D C . A D i st i n e t N u m e r i e a l M o d e l for G anr u lar A s s e m b li e s . G e Oet e hn i q u e , 19 7 9 , 2 9 : 4 7 9 C u n d a l l P A . F o rl l l u l a t i o n o f a T hr e e 一 di m e n s i o n a l D i s t i n e t E l e m e n t M o d e l? P 叭 1 . I n t J R o e k M e e h , M i n S e i & G e - o m e e h A b s tr, 1 9 8 8 , 2 5 ( 3) : 10 7 10 H叭 R , C un d a l l P A , L e m o s J . F o mr u l a t i o n o f a T hr e e 一 d i - m e n s i o n a l D i s ti n e t E l e m e in M o d e l? P a rt 11 . I n t J R o e k M e e h , M i n S e i & G e o m e e h A b s t几 1 9 8 8 , 2 5 (3) : 1 1 7 1 1 L an g e F F . E fe e t o f I n t e rp a rt i e l e P o t e n t i a l s o n P art i e l e P a e k i n g of r C e r am i e P r o e e s s i n g , P o w d e r & G r a i n s 9 3 . T h o mr T, e d . B a l k e m a : R o t e r d am , 1 9 9 8 . 1 8 7 1 2 C un d a l l P A . D i s t i n e t E l e m e n t M o d e l s o f R o e k an d s o il S trU e ut r e . I n : A n a ly ti e a l an d C o m P u at t i o n a l M e ht o d s i n E n g i n e e r i n g oR e k M e c h an i e s , B or wn E 工 e d . L o n d o n : A l l e n & U n w i n , 19 8 7 . 12 9 13 B o u v ar d D , L an g e F .F eR l a t i o n b e wt e e n P e r e o lat i o n a n d P art i e l e C o o r d i n at i o n i n B i n a勺 P o w d e r M ixt u r e s . A c at m e at ll M a t e r, 19 9 1 , 3 9( 1 2 ) : 30 8 3 C o m P u t e r S im ul at i o n o f R an d o m P a e k i n g o f S P h e r i e a l P art i e l e s hC e ng uY n af ng, G u o hS ij i u ,n La i H办 , i M a t e ir ia S e i e cn e an d nE g in e e r in g S e h o l , U S T B e ij in g , B e ij ing 10 0 0 8 3 , C h in a AB S T R A C T A e o mP ut e r s im u l at i o n of r e v o lut i o n o f r an d o m s trU e trU e s o f s Ph e ir c a l Part i c l e s i n tw o d im e n - s ion s h as b e e n Per of n n e d 勿 ht e d i s itn e t e l e m e nt m e ht o d . hT e of r e e s am o n g ht e Part i e l e s e o n s ide r e d w e r e g a v lty, 衍 e ti o n fo r c e an d 、 恤n D e r Wh a l s I ent ar e t i o n (V D机 ) . hT e m ot i o n o f Part i e l e s i n v o l v e d tr a n s lat i o n an d r o at i on . Fo r 10 0脚 Part i e l e s , if 示 e t i o n an d V D机 w e r e n o t C o n s id e r e d , ht e Pa c ik n g a s s e m b l y Wa s a m i x trU e o f r an d om Pa e k l n g an d or d e r Pa e ik n g w iht Pa e k i n g d e n s iyt o f 0 . 8 6 9 6 , if 创 e t i o n an d V D W I w e r e t a k e n int o a e e o u n t , ht e Pac ikn g d e n s iyt w a s 0 . 8 2 1 3 . oT s t u d y ht e Pr o e e d ur e o f Pa e k i n g , Pa e ik n g d e n s it i e s w iht d i fe r ent t im e , Part i e l e tr aj e c t o yr an d e o o dr i n at i o n n um b e r w e r e e x am i n e d . hT e Pa c k i n g s trU e trU e s o f b i n娜 s y s t e m W iht id fe r e n t e o n d it ion s w e r e a l s o g i v e n . K E Y W O R D S Part i e l e arn d o m p a e ik n g : d i s t i n e t e l e m e in m e ht o d : e o o r d i n a t i o n mnU b e r : p a c k i n g d e n s iyt