An improved 3D DEM-FEM contact detection algorithm for the interaction simulations between particles and structures

The interaction simulations between particles and structures are often performed in the context of the combined discrete-finite element (DEM-FEM) method, where an efficient, robust and accurate contact algorithm for challenging contact problems is essential. A three-dimensional (3D) discrete and finite element contact algorithm, named ZGL here, has been proposed in our research group (Zang, M.Y. et al., 2011. A contact algorithm for 3D discrete and finite element contact problems based on penalty function method. Comput. Mech. 48, 541–550.). Despite being quite well-established, ZGL algorithm may lack efficiency, robustness and accuracy in certain situations, e.g., when special mesh pattern is applied. However, since such cases are readily to appear in practice engineering application, an algorithm named DZCell is developed to resolve these issues and thus significantly applicable to challenging contact problems. The proposed DZCell algorithm includes an improved global phase to directly find the potential segments rather than the nearest node for discrete elements as contact counterparts. Furthermore, both the memory cost and time consumption of this algorithm are linear, and the algorithm is readily to extend to parallel computing. Several numerical examples demonstrate the achievable improvements in terms of efficiency, robustness and accuracy for 3D contact analysis and validate the capability of the DZCell algorithm in the granular science and mechanical engineering.

[1]  David J. Benson,et al.  Sliding interfaces with contact-impact in large-scale Lagrangian computations , 1985 .

[2]  H. Shen,et al.  Comparisons of physical experiment and discrete element simulations of sheared granular materials in an annular shear cell , 2009 .

[3]  Jianqiao Ye,et al.  Discrete element modeling of the microbond test of fiber reinforced composite , 2010 .

[4]  Fei Wang,et al.  The modelling of multi‐fracturing solids and particulate media , 2004 .

[5]  Mengyan Zang,et al.  Investigation of impact fracture behavior of automobile laminated glass by 3D discrete element method , 2007 .

[6]  Aibing Yu,et al.  DEM study of the flow of cohesive particles in a screw feeder , 2014 .

[7]  Larsgunnar Nilsson,et al.  A contact searching algorithm for general contact problems , 1989 .

[8]  Jeffrey Lawrence Slade.  DEVELOPMENT OF A NEW OFF-ROAD RIGID RING MODEL FOR TRUCK TIRES USING FINITE ELEMENT ANALYSIS TECHNIQUES , 2009 .

[9]  Mengyan Zang,et al.  LC-Grid: a linear global contact search algorithm for finite element analysis , 2014 .

[10]  Eric Perkins,et al.  A fast contact detection algorithm insensitive to object sizes , 2001 .

[11]  Zhenhan Yao,et al.  FFS contact searching algorithm for dynamic finite element analysis , 2001 .

[12]  Hiroshi Shimizu,et al.  2D FE–DEM analysis of tractive performance of an elastic wheel for planetary rovers , 2016 .

[13]  Hiroshi Nakashima,et al.  Algorithm and implementation of soil–tire contact analysis code based on dynamic FE–DE method , 2004 .

[14]  P. Hazell,et al.  An accurate and robust contact detection algorithm for particle‐solid interaction in combined finite‐discrete element analysis , 2015 .

[15]  Mark O. Neal,et al.  Contact‐impact by the pinball algorithm with penalty and Lagrangian methods , 1991 .

[16]  Zhaolin Gu,et al.  Discrete element simulation of particle flow in arbitrarily complex geometries , 2011 .

[17]  Jidong Zhao,et al.  Realistic Generation and Packing of DEM Sand Samples , 2012 .

[18]  Mengyan Zang,et al.  A three-dimensional computational framework for impact fracture analysis of automotive laminated glass , 2015 .

[19]  Hiroshi Nakashima,et al.  FE-DEM Analysis of the Effect of Tread Pattern on the Tractive Performance of Tires Operating on Sand , 2009 .

[20]  Mengyan Zang,et al.  Analysis of rigid tire traction performance on a sandy soil by 3D finite element–discrete element method , 2014 .

[21]  Mengyan Zang,et al.  A contact algorithm for 3D discrete and finite element contact problems based on penalty function method , 2011 .

[22]  David R. Owen,et al.  Performance comparisons of tree‐based and cell‐based contact detection algorithms , 2007 .

[23]  Susan Fisher,et al.  An improved finite-element contact model for anatomical simulations , 2003, The Visual Computer.

[24]  Tan Yuanqiang,et al.  A virtual-surface contact algorithm for the interaction between FE and spherical DE , 2016 .

[25]  Eric Perkins,et al.  A contact algorithm for partitioning N arbitrary sized objects , 2004 .

[26]  Bernhard Peters,et al.  DEM–FEM coupling simulations of the interactions between a tire tread and granular terrain , 2015 .

[27]  Xiao Yun Xu,et al.  An improved version of RIGID for discrete element simulation of particle flows with arbitrarily complex geometries , 2014 .

[28]  M. Shaw,et al.  Investigation of comminution in a Wiley Mill: Experiments and DEM Simulations , 2013 .

[29]  A. Munjiza,et al.  NBS contact detection algorithm for bodies of similar size , 1998 .

[30]  Hu Guoming,et al.  A new algorithm for contact detection between spherical particle and triangulated mesh boundary in discrete element method simulations , 2013 .

[31]  David J. Benson,et al.  A single surface contact algorithm for the post-buckling analysis of shell structures , 1990 .

[32]  Runyu Yang,et al.  Discrete particle simulation of particulate systems: A review of major applications and findings , 2008 .

[33]  Kaiming Xia.  Finite element modeling of tire/terrain interaction: Application to predicting soil compaction and tire mobility , 2011 .

[34]  David R. Owen,et al.  A combined finite/discrete element simulation of shot peening processes – Part II: 3D interaction laws , 2000 .

[35]  Mica Grujicic,et al.  Computational Analysis of Mine Blast on a Commercial Vehicle Structure , 2007 .

[36]  Yuanqiang Tan,et al.  Discrete element method (DEM) modeling of fracture and damage in the machining process of polycrystalline SiC , 2009 .

[37]  Rimantas Kačianauskas,et al.  Discrete element method and its application to the analysis of penetration into granular media , 2004 .

[38]  S. Wang,et al.  The inside–outside contact search algorithm for finite element analysis , 1997 .

[39]  Mengyan Zang,et al.  An approach to combining 3D discrete and finite element methods based on penalty function method , 2010 .

[40]  Jidong Zhao,et al.  3D generation of realistic granular samples based on random fields theory and Fourier shape descriptors , 2014 .

[41]  Yunming Yang,et al.  Three‐dimensional finite element modeling of tire/ground interaction , 2012 .

[42]  David R. Owen,et al.  Parallelised finite/discrete element simulation of multi‐fracturing solids and discrete systems , 2001 .

[43]  Siegmar Wirtz,et al.  Experimental and numerical investigation on the influence of particle shape and shape approximation on hopper discharge using the discrete element method , 2013 .

[44]  Jingwei Zheng,et al.  GPU-based parallel algorithm for particle contact detection and its application in self-compacting concrete flow simulations , 2012 .

[45]  Eugenio Oñate,et al.  The Double Hierarchy Method. A parallel 3D contact method for the interaction of spherical particles with rigid FE boundaries using the DEM , 2016 .

[46]  Youssef M A Hashash,et al.  Three‐dimensional discrete element simulation for granular materials , 2006 .