第 35 卷 第 1 期 力 学 与 实 践 离散元法及其耦合算法的研究综述 徐 爽 2) 朱浮声 2013 年 2 月 1) 张 俊 (东北大学资源与土木工程学院 沈阳 110819) 徐爽 1978 年 12 月生 2001 年毕业于东北大学土木工程专业 2005 年获硕士学位, 东北大学结构工程专业博士研究生 师从于朱浮声教授 从事于非均匀脆性介质破 坏过程的数值模拟方面的研究. 摘要 介绍了离散单元法的基本理论及其研究现状 以及离散单元法与有限单元法 边界单元法 界面单 元法等数值计算方法耦合的研究现状和最新进展 并讨论了离散单元法今后的发展趋势及亟待解决的问题. 关键词 离散单元法 数值方法 耦合 中图分类号: O482 文献标识码: A DOI 10.6052/1000-0879-12-269 A OVERVIEW OF THE DISCRETE ELEMENT METHOD AND ITS COUPLING ALGORITHMS1) XU Shuang2) ZHU Fusheng ZHANG Jun (School of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China) Abstract The basic theory and the current development of the discrete element method are reviewed in this paper, as well as those of the coupling algorithms of the discrete element method with the finite element method, the boundary element method, and the interface element method. The future development of the discrete element method is also discussed. Key words discrete element, numerical method, coupling 引 言 脆性无序介质具有各项异性 非均匀性 不连续 性等特点 因此 对于脆性无序介质 (如混凝土 岩 石 玻璃和陶瓷等) 的变形 微细观损伤乃至宏观破 坏全过程的研究 已成为各界十分关注的研究热点 之一. 目前 在对脆性无序介质的研究中 常用的数 值分析方法有有限差分法 有限单元法和边界单元 法等. 但是 这几种方法都是以连续介质为出发点 2012 07 04 收到第 1 稿 2012 09 26 收到修改稿. 1) 国家自然科学基金资助项目 (51008053). 2) E-mail: phd.xushuang@gmail.com 而且往往是基于小变形的假设 因此所得的结果与 实际相差甚远 为了解决这个问题 离散单元法就应 运而生了. 刘凯欣等 [1] 曾对离散单元法给出过明确定义 通过物理元的单元离散方式并构成具有明确物理意 义的节点关系建立有限离散模型的数值计算方法. 经过三十多年的发展 离散单元法在理论和应 用研究方面取得了长足的进步. 在理论方面 从刚性 元发展到可变形元 从单纯的离散元模拟发展到该
1 9... 1 ( 1(a)). ( 1(b)).. (a) (b) 1... 2 Cundall 1971 [2-3].. 1974 [4].. 1978 FORTAN [5]. Cundall Strack [6-7] BALL Drescher [8] BALL. Cundall [9] 1980. Cundall UDEC (universal distinct element code) 1985 [10-11]. UDEC.. 3DEC (3-dimensinal distince element code) Cundall ITASCA 1986 [12]. UDEC. BALL TRUBAL [13] BALL.
10 2013 35 [14-15] 3DEC [16-18].... [19-20] 1986. [21].. [22-30]. [31-33]. [34-35]. [36]. Liu [37-41],. Mohr--Coulomb. ( ). [42].. 3... 3.1 Dowding [43-44]... 1991 Pan Reed [45]. [46]..... ( 2)...
1 11 2 [40].. [47-48] --.... 3 ( 3) ( ) ( ). ( ). -- ( ). 3 [47-48] 1996 Kun Herrmann [49] --. 1997 Muhlhaus [50]. ( ) --. Mohammadi [51] Ransing [52]. 3.2 Lorig. Lorig Brady [53]. 1983 Lemos [54] [55]. 1984 Lorig [56].. HYDEBE (hybrid discrete element boundary element) CRE-
12 2013 35 ATE BOUND BLOCK. [57]... [58]... Γ1.. Γ2.. [59].. 3.3 -- Zubelewicz Bazant [60].. 4............ 1,.., 2003,33(4): 483-490 (Liu Kaixin, Gao Lingtian. A review on the discrete element method. Advances in Mechanics, 2003, 33(4): 483-490 (in Chinese)) 2 Cundall PA. A computer model for simulating progressive large scale movements in blocky systems. In: Proceedings of the Symposium of the International Society of Rock Mechanics, Nancy, France, 1971 3 Cundall PA. The measurement and analysis of acceleration
1 13 in rock slopes. [PhD Thesis]. London: University of London, 1971 4 Cundall PA. Rational design of tunnel supports: A computer model for rock mass behavior using interactive graphics for the input and output of geometrical data. Tech Report MRD-2-74, Missouri River Division, US Army Corps of Engineers, 1974 5 Cundall PA, Marti J, Beresford PJ, et al. Computer modeling of jointed rock masses. Tech Report N-78-4, US Army Engineer Waterways Experiment Station, Vicksburg, Mississippi, 1978 6 Cundall PA, BALL A program to model granular media using the distinct element method. London: Dames & Moore Advanced Technology Group, 1978 7 Strack ODL, Cundall PA. The Distinct Element Method as a Tool for Research in Granular Media, Part. Minnesota: University of Minnesota, 1978 8 Drescher A. Photoelastic verification of a mechanical modol for the flow of a granular material. J Mech Physics Solids, 1972, 20: 337-351 9 Cundall PA. UDEC A generalized distinct element program for modeling jointed rock. Report PCAR-1-80, Peter Cundall Associates, European Research Office, US Army, 1980 10 Cundall PA, Hart RD. Development of generalized 2-D and 3-D distinct element programs for modeling jointed rock. ITASCA Consulting Group, Misc Paper SL-85-1, US Army Corps of Engineers, 1985 11 UDEC News and World Report, Official Publication of the UDEC Users Group, ITASCA Consulting Group, Inc., 1985 12 3DEC 3-D Distinct Element Code, ITASCA Consulting Group, 1987 13 Cundall PA. Formulation of a three-dimensional distinct element model-part. A scheme to detect and represent contacts in a system composed of many polyhedral blocks. Int J Rock Mech Min Sci & Geomch Abstr, 1988, 25(3): 107-116 14 Lorig LJ. A hybrid computational model for excavation and support desing in jointed media. [PhD Thesis]. Minnesota: University of Minnesota, 1984 15 Hart R, Cundall PA, Lemos J. Formultaion of a threedimensional distinct element model-parte II. Mechanical calculations for motion and interaction of a system composed of many polyhedral blocks. Int J Rock Mech Min Sci & Geomech Abstr, 1988, 25(3): 117-125 16 Shen B, Stephansson O. Cyclic loading characteristics of joints and rock bridges in a jointed rock specimen. In: Proe of Int Conference on Rock Joints, Leon Norway, 1990 17 Stephansson O, Shen B, Lemos J. Modeling of excavation, thermal loading and bentonite swelling pressure for a waste repository. In: Int High-level Radioactive Waste Management Conference, Las Vegas, Nevada, USA, 1991 18 Stephansson O, Shen B. Modeling of rock masses for site location of a nuclear waste repository. In: 7th Int. Congress of Int. Society for Rock Mechanics, Aachen, Germany, 1991 19.. :, 1986. 32-37 (Wang Yongjia. Discrete element method a numerical method suitable for analysis in jointed rocks. In: First National Rock Mechanics Numerical Calculations and Model Tests Seminar Proceedings, 1986. 32-37 (in Chinese)) 20.. :, 1986. 43-46 (Chou Wanxi. Principle of the discrete element method and its application in rock mechanics and engineering. In: First National Rock Mechanics Numerical Calculations and Model Tests Seminar Proceedings. 1986. 43-46 (in Chinese)) 21.. :, 1986 (Zhang Qing. Fundamentals of Rock Mechanics. Beijing: China Railway Publishing House. 1986 (in Chinese)) 22,.. ( ), 1986. 88-94 (Wang Yongjia, Xing Jibo. A new method of ore research discrete element method. In: Second National Mining Conference Proceedings (Metallic Volume), 1986. 88-94 (in Chinese)) 23,,.., 1987, 39(2): 20-26 (Wang Yongjia, Liu Xingguo, Xing Jibo. Research on applications of discrete element method in ore drawing for caving system. Nonferrous Metals, 1987, 39(2): 20-26 (in Chinese)) 24.. [ ]. :, 1987 (Dai Qing. The study on stability of electric rake bottom structure in natural caving method. [PhD Thesis]. Beijing: University of Science and Technology Beijing, 1987 (in Chinese)) 25,.., 1988, (2): 148-153 (Xing Jibo, Wang Yongjia. Discrete element simulation of ore drawing by caving method in mining. Journal of Northeast University of Technology, 1988, (2): 148-153 (in Chinese)) 26.. [ ]. :, 1988 (Lai Senhua. The research on dynamics problems of caving method with bottom structure. [PhD Thesis]. Beijing: University of Science and Technology Beijing, 1988 (in Chinses)) 27,,.. ( ), 1989, (2): 3-8 (Zheng Yongxue, Wang Yongjia, Yao Zanxun. Study on the caving mechanism of natural caving. Nonferrous Metals (Mining Section), 1989, (2): 3-8 (in Chinese)) 28,.. ( ), 1988. 22-26 (Wang Yongjia, Xing jibo. Computer simulation of caving method ore-drawing. In: The National Non-metallic Minerals Conference Proceedings(b), 1988. 22-26 (in Chinese))
14 2013 35 29 Wang Yongjia, Xing Jibo. Computerized simulation of ore drawing of caving system. In: Proe. 2nd Wold Congress on Non-Metallic Minerals, 1989 30.. [ ]. :, 1989 (Chen Ronghua. The study on stability of bottom structure in natural caving method. [Master Thesis]. Shenyang: Northeast University of Technology, 1989 (in Chinese)) 31.. [ ]. :, 1987 (Shen Baotang. Discrete element analysis of slope stability. [Master Thesis]. Shenyang: Northeast University of Technology, 1987 (in Chinese)) 32,...,, 1988 (Wang Yongjia, Shen Baotang. Discrete element analysis of slope stability. Structural Model Destruction Test Symposium Proceeding. Nanjing, Houhai University, 1988 (in Chinese)) 33,.., 1989, (4): 340-354 (Shen Baotang, Wang Yongjia. Discrete element method for slope failure mechanism. Journal of Northeast University of Technology, 1989, (4): 340-354 (in Chinese)) 34.. [ ]. :, 1988 (Chang Laishan. The analysis of discrete element method on failure mechanism of guangdong yunfu pyrite slope. [Master Thesis]. Ma Anshan: Ministry of Metallurgical Industry Maanshan Mining Research Institute (in Chinese)) 35,.., 1987, 6(4): 337-346 (Da He, Xing Jibo. Analysing panshi nickle mine slope stability by distinct element method. Chinese Journal of Rock Mechanics and Engineering, 1987, 6(4): 337-346 (in Chinese)) 36.. : ( ),, 1987 (Yao Jianguo. The numerical simulation of block element on strata movement and deformation in caving influence. First National calculation Rock and Soil Mechanics Symposium (a), E mei, 1987 (in Chinese)) 37 Liu KX, Gao LT, Tanimura S. Application of discrete element method in impact problems. JSME Int J Ser A, 2004, 47: 138-145 38 Liu KX, Gao LT. The application of discrete element method in solving three dimensional impact dynamics problems. Acta Mech Sol, 2003, 16: 256-261 39 Liu KX, Liu WF. Application of discrete element method for continuum dynamic problems. Arch Appl Mech, 2006, 76: 229-243 40 Cheng M, Liu WF, Liu KX. New discrete element models for elastoplastic problems. Acta Mech Sin, 2009, 25: 629-637 41 Shan L, Cheng M, Liu KX, et al. New discrete element models for three-dimensional impact problems. Chin Phys Lett, 2009, 26(12): 120202 42 Shan L, Cui N, Liu KX, et al. A new combined scheme of discrete element method and meshless method for numerical simulation of impact problems. Computer Modeling in Engineering & Science, 2012, 83(4): 353-383 43 Dowding CH, Gilbert C. Dynamic stability of rock slopes and high frequency traveling waves. Journal of Geotechnical Engineering, 1988, 114(10): 1069-1088 44 Dowding CH, Belytschko TB, Yen HJ. A coupled finite element-rigid block method for transient analysis of rock caverns. Int J Num Analy Methods Geomech, 1983, 7: 117-127 45 Pan XD, Reed MB. A coupled distinct element finite element method for large deformation analysis of rock masses. International Journal of Rock Mechanics and Mining Science & Geomechanics Abstract, 1991, 28(1): 93-99 46,.., 2003, 20(6): 476-482 (Xu Jianlong, Tang Zhiping. Combined discrete/finite element multiscale numerical method and its application. Chinese Journal of Computational Physics, 2003, 20(6): 476-482 (in Chinese)) 47. --., 1999 (Xing Jibo, Introduction of the Beam-Particle Model. Beijing: Seismological Press, 1999 (in Chinese)) 48,,... 2005, 24(4): 570-574 (Zhang Dehai, Zhu Fusheng, Xing Jibo, et al. Numerical simulation of fracture propagation of anisotropic rock-like materials. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(4): 570-574 (in Chinese)) 49 Kun F, Herrmann HJ. Fragmentation of colliding discs. Int J of Mod Phys C, 1996, 7(6): 837-855 50 Muhlhaus HB, Sakaguchi H, Wei Y. Particle based modeling of dynamic fracture in jointed rock. In: Yuan J X ed., Proc. 9th Int. Conf. on Numerical Methods in Geomechanics, Balkema, Rotterdam, 1997 51 Mohammadi S, Owen DRJ, Peric D. A combined finite/discrete element algorithm for delamination analysis of composites. Finite Elements in Analysis and Design, 1998, 28: 321-336 52 Ransing RS, Gethin DT, Khoei AR, et al. Powder compaction modeling via the discrete and finite element method. Materials and Design, 2000, 21: 26-269 53 Lorig LJ, Brady BHG. A hybrid discrete element-boundary element method of stress analysis. In: Proc. 23rd U.S. Symposium on Rock Mechanics, California, 1982. 628-636 54 Lemos JV. A hybrid distinct element computational model for the half-plane. [Master Thesis]. Minnesota: University of Minnesota, 1983 55 Lemos JV, Brady BHG. Stress distribution in a jointed and fractured medium. In: the 24th U.S. Symposium on Rock Mechanics, New York, 1983 ( 19 )
1 19 17 Butterfield G R. Application of thermal anemometry and high-frequency measurement of mass flux to aeolian sediment transport research. Geomorphology, 1999, 29 (1-2): 31-58 18 Arnold S. Development of the saltation system under controlled environmental conditions. Earth Surface Processes and Landforms, 2002, 27: 817-829 19 Wiggs GFS, Atherton RJ, Baird AJ. Thresholds of aeolian sand transport: Establishing suitable values. Sedimentology, 2004, 51: 95-108 20 Anderson RS, Haff PK. Simulation of eolian saltation. Science, 1988, 241: 820-823 21 Spies PJ, McEwan IK, Butterfield GR. One-dimensional transition behaviour in saltation. Earth Surface Processes and Landforms, 2000, 25: 505-518 22 McEwan IK, Willetts BB. Sand transport by wind: A review of the current conceptual model. In: Pye K, editor. The Dynamics and Environmental Context of Aeolian Sedimentary Systems. London, 1993 23 Spies PJ. The Transport of Sand in Unsteady Winds. [PhD Thesis]. Scotland: University of Aberdeen, 1996 24 Spies PJ, McEwan IK. Equilibration of saltation. Earth Surface Processes and Landforms, 2000, 25: 437-453 25 Sørensen M. On the effect of time variability of the wind on rates of aeolian sand transport. Aarhus: University of Aarhus, 1997 26,,.., 2005, 3(1): 82-87 (Yan Min Sun Baoping Lu Zhenyou et al. Discrete particle dynamics model and dynamic process for blown sand sport. Science of Soil and Water Conservation 2005, 3(1): 82-87 (in Chinese)) 27 Huang N, Zheng XJ, Zhou YH, et al. Simulation of windblown sand movement and probability density function of liftoff velocities of sand particles. Journal of Geophysical Research, 2006, 111: D20201 28 Zhou YH LI WQ Zheng XJ. PDM simulations of stochastic collisions of sandy grain-bed with mixed size in aeolian sand saltation. Journal of Geophysical Research 2006 111: D15108 29,,.., 2008, 26(6): 66-69 (Yan Min Yin Jianye Sun Baoping. DEM modeling of sand group dynamics. Science & Technology Review 2008, 26(6): 66-69 (in Chinese)) 30 Zheng XJ. Mechanics of Wind-blow Sand Movemnts. Berlin Springer, 2009 31,,.., 2006, 36(4): 580-590 (Yang Bin, Wang Yuan Wang Dawei. Development of wind-sand flow measurement technique. Advances in Mechanics, 2006, 36(4): 580-590 (in Chinese)) 32,,.., 2010, 24(1): 55-64 (Wang Yuan Yang Bin, Wang Dawei. Advances in wind-blown sand flow optical measurement and image processing techniques. Journal of Experiments in Fluid Mechanics, 2010, 24(1): 55-64 (in Chinese)) ( : ) 14 ) 56 Lorig LJ. A hybrid computational model for excavation and support design in jointed media. [PhD Thesis]. Minnesota: University of Minnesota, 1984 57,.., 1987 (Wang Yongjia, Xing Jibo. The coupling of discrete element method and boundary element method. In: The Application of the Boundary Element Method in Rock Mechanics and Engineering Proceedings, 1987 (in Chinese)) 58,,.., 1996, 15: 231-235 (Zhou Xiaoqing, Wang Yuanhan, Ha Qiuling. External coupling computation between distinct element method and boundary element method. Chinese Journal of Rock Mechanics and Engineering, 1996, 15: 231-235 (in Chinese)) 59,,. -., 2001, (2): 24-28 (Jin Feng, Jia Weiwei, Wang Guanglun. Coupling model of distinct element-boundary element. Journal of Hydraulic Engineering, 2001, (2): 24-28 (in Chinese)) 60 Zubelewicz A, Bazant ZP. Interface element modeling of fracture in aggregate composites. Journal of Engineering Mechanics, 1987, 113(11): 1619-1630 ( : )