Development of contact algorithm for three‐dimensional numerical manifold method

SUMMARY This paper customizes a contact detection and enforcing scheme to fit the three-dimensional (3-D) numerical manifold method (NMM). A hierarchical contact system is established for efficient contact detection. The mathematical mesh, a unique component in the NMM, is utilized for global searching of possible contact blocks and elements, followed by the local searching to identify primitive hierarchies. All the potential contact pairs are then transformed into one of the two essential entrance modes: point-to-plane and crossing-lines modes, among which real contact pairs are detected through a unified formula. The penalty method is selected to enforce the contact constraints, and a general contact solution procedure in the 3-D NMM is established. Because of the implicit framework, an open-close iteration is performed within each time step to determine the correct number of contact pairs among multi-bodies and to achieve complete convergence of imposed contact force at corresponding position. The proposed contact algorithm extensively utilizes most of the original components of the NMM, namely, the mathematical mesh/cells and the manifold elements, as well as the external components associated with contacts, such as the contact body, the contact facet and the contact vertex. In particular, the utilization of two mutually approaching mathematical cells is efficient in detecting contacting territory, which makes this method particularly effective for both convex and non-convex bodies. The validity and accuracy of the proposed contact algorithm are verified and demonstrated through three benchmark problems. Copyright © 2013 John Wiley & Sons, Ltd.

[1]  Jae-Joon Song,et al.  NEW CONTACT-DEFINITION ALGORITHM USING INSCRIBED SPHERES FOR 3D DISCONTINUOUS DEFORMATION ANALYSIS , 2011 .

[2]  P. A. Cundall,et al.  FORMULATION OF A THREE-DIMENSIONAL DISTINCT ELEMENT MODEL - PART I. A SCHEME TO DETECT AND REPRESENT CONTACTS IN A SYSTEM COMPOSED OF MANY POLYHEDRAL BLOCKS , 1988 .

[3]  Guoxin Zhang,et al.  SIMULATION OF TOPPLING FAILURE OF ROCK SLOPE BY NUMERICAL MANIFOLD METHOD , 2010 .

[4]  Guowei Ma,et al.  Numerical analysis of 2-D crack propagation problems using the numerical manifold method , 2010 .

[5]  J. Bonet,et al.  An alternating digital tree (ADT) algorithm for 3D geometric searching and intersection problems , 1991 .

[6]  Gen-Hua Shi,et al.  Manifold Method of Material Analysis , 1992 .

[7]  Yu-Fei Wu,et al.  Numerical manifold method based on the method of weighted residuals , 2005 .

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

[9]  Y. H. Zhang,et al.  Formulation of a Three-dimensional Numerical Manifold Method with Tetrahedron and Hexahedron Elements , 2008 .

[10]  Xian-Zhong Mao,et al.  Least‐square‐based radial basis collocation method for solving inverse problems of Laplace equation from noisy data , 2010 .

[11]  Chuangbing Zhou,et al.  Modeling Unconfined Seepage Flow Using Three-Dimensional Numerical Manifold Method , 2010 .

[12]  P. Wriggers,et al.  On contact between three-dimensional beams undergoing large deflections , 1997 .

[13]  Larsgunnar Nilsson,et al.  A unified contact algorithm based on the territory concept , 1996 .

[14]  Guowei Ma,et al.  Modeling complex crack problems using the numerical manifold method , 2009 .

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

[16]  Mao Kurumatani,et al.  Finite cover method with multi‐cover layers for the analysis of evolving discontinuities in heterogeneous media , 2009 .

[17]  John R. Williams,et al.  A linear complexity intersection algorithm for discrete element simulation of arbitrary geometries , 1995 .

[18]  Guowei Ma,et al.  A NUMERICAL MANIFOLD METHOD FOR PLANE MICROPOLAR ELASTICITY , 2010 .

[19]  Louis Ngai Yuen Wong,et al.  Frictional crack initiation and propagation analysis using the numerical manifold method , 2012 .

[20]  Dian-Qing Li,et al.  A three-dimensional numerical manifold method based on tetrahedral meshes , 2009 .

[21]  Guowei Ma,et al.  Footwall slope stability analysis with the numerical manifold method , 2011 .

[22]  Kenjiro Terada,et al.  An integrated procedure for three‐dimensional structural analysis with the finite cover method , 2005 .

[23]  Bernard Amadei,et al.  Extensions of discontinuous deformation analysis for jointed rock masses , 1996 .

[24]  G. Ma,et al.  A new way to treat material discontinuities in the numerical manifold method , 2011 .

[25]  Kenjiro Terada,et al.  Finite cover method for linear and non‐linear analyses of heterogeneous solids , 2003 .

[26]  Kenjiro Terada,et al.  Eulerian finite cover method for quasi‐static equilibrium problems of hyperelastic bodies , 2006 .

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

[28]  Guowei Ma,et al.  DEVELOPMENT OF 3D NUMERICAL MANIFOLD METHOD , 2010 .

[29]  Q. H. Jiang,et al.  A Model of Point-to-Face Contact for Three-Dimensional Discontinuous Deformation Analysis , 2004 .

[30]  M. R. Yeung,et al.  A model of edge-to-edge contact for three-dimensional discontinuous deformation analysis , 2007 .

[31]  Ted Belytschko,et al.  The splitting pinball method for contact-impact problems , 1993 .

[32]  Lin Gao,et al.  Static and dynamic stability analysis using 3D-DDA with incision body scheme , 2006 .

[33]  Qinghui Jiang,et al.  Validation of block theory and three-dimensional discontinuous deformation analysis as wedge stability analysis methods , 2003 .