非球体填充的组合球模型及松弛算法-NSFC.PDFVIP

第 26 卷 第 3 期 计 　　算 　　物 　　理 Vol . 26 ,No . 3 2009 年 5 月 　CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 　 May , 2009 　 ( ) [文章编号] 　1001246X 2009 非球体填充的组合球模型及松弛算法 李水乡 , 　赵 　健 (北京大学工学院力学与空天技术系 湍流与复杂系统国家重点实验室 , 北京 　100871) [摘 　要] 　现有的松弛算法由于仅用于球填充而只考虑颗粒的平动 , 故提出考虑非球体转动的改进松弛算法并 采用组合球模型 , 使其能够模拟任意形状非球体的随机填充以及多种非球体的混合填充. 用多个球体的外包络 面近似一个非球体外形的组合球模型 , 将非球体之间的接触转化为球体之间的接触 , 从而简化并统一非球体接 触判断算法. 通过引入非球体的转矩和转角松弛机制 , 使改进松弛算法克服了 “自锁”现象 , 并能生成非球体 的随机密填充. 算例表明 , 填充结果与现有的数值模拟及实验结果相符. [ 关键词] 　填充 ; 堆积 ; 非球形颗粒 ; 松弛算法 ; 随机密填充 ; 组合球 [ 中图分类号] 　O411 ; TP391 [文献标识码] 　A 0 　引言 填充问题的研究在理论和实践中都具有极其重大的意义. 自上世纪 50 年代末以来 ,数值模拟逐渐成为 研究填充问题的主要手段 ,并已形成了序列添加 、集合重排 、前缘法和优化法四类方法. 在当前的填充研究 [ 1] 中 ,球填充的研究占主导地位且已较成熟. 球填充中填充率的上下界等参数都有广泛认同的结果 ,不同球 径的组合填充 ,如双径球 、三径球 、四径球 、直径正态分布球等 ,都有比较系统的研究. 为简化模型 , 以往的填 充研究往往以等量直径把非球体简化为球体. 然而近年来的研究表明 ,球填充的填充率与非球体填充的填充 率之间差距 明显. Torquato 和 Stillinger 等对椭球填充 的数值模拟[2 ] 及实验[3 ] 表 明 , 当椭球 的三轴 比为 ( ) [4 ] 125 :1 :08时 ,椭球随机密填充 random close packing 的填充率可达 74 % ,平均配位数可达 12 . Williams 等 ( ) ( 对球柱体 spherocylinder 填充的研究表明 ,其随机密填充率可达 70 %. 这些结果均远高于球填充 随机密填 充率 64 % ,平均配位数 6) . 这些研究表明 , 以往将所有非球体都简化为球体是不适当的 ,非球体填充的性质 与球填充的性质有较大的差别. 在非球体模型和接触判断方面 ,现有的方法大致可分为两类 :一类是解析或近似解析方法 ,这类研究均 针对一种特定的非球体来建立其几何模型和接触判断算法 ,如椭球[5 ]