A Model of Point-to-Face Contact for Three-Dimensional Discontinuous Deformation Analysis

SummaryThe key to three-dimensional discontinuous deformation analysis (3D DDA) is a rigorous contact theory that governs the interaction of many three-dimensional blocks. This theory must provide algorithms to judge contact types and locations and the appropriate state of each contact, which can be open, sliding or locked. This paper presents a point-to-face contact model, which forms a part of the contact theory, to be used in 3D DDA. Normal spring, shear spring and frictional force submatrices are derived by vector analysis and the penalty method. Also given are the “open-close” iteration criteria and operations performed for different changes in contact state. Sliding at a contact can occur in any direction parallel to the contact face, as opposed to one of two directions in two-dimensional DDA. This point-to-face contact model has been implemented into a 3D DDA computer program, and numerical results from several test cases demonstrate the validity of the model and the capability of the program.

[1]  O. C. Zienkiewicz,et al.  A note on numerical computation of elastic contact problems , 1975 .

[2]  Noboru Kikuchi,et al.  Finite element methods for constrained problems in elasticity , 1982 .

[3]  Robert D. Cook,et al.  An iterative procedure for finite-element stress analysis of frictional contact problems , 1984 .

[4]  Anil Chaudhary,et al.  A SOLUTION METHOD FOR PLANAR AND AXISYMMETRIC CONTACT PROBLEMS , 1985 .

[5]  A. Klarbring A mathematical programming approach to three-dimensional contact problems with friction , 1986 .

[6]  Anil Chaudhary,et al.  A solution method for static and dynamic analysis of three-dimensional contact problems with friction , 1986 .

[7]  R. J. Melosh,et al.  Solving discretized contact problems using linear programming , 1987 .

[8]  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 .

[9]  Richard E. Goodman,et al.  Generalization of two‐dimensional discontinuous deformation analysis for forward modelling , 1989 .

[10]  Genki Yagawa,et al.  A dynamic contact buckling analysis by the penalty finite element method , 1990 .

[11]  M. R. Yeung Analysis of a mine roof using the DDA method , 1993 .

[12]  K. Yamazaki,et al.  Penalty method for three-dimensional elastic contact problems by boundary element method , 1994 .

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

[14]  Nicholas Sitar,et al.  Kinematics and Discontinuous Deformation Analysis of Landslide Movement , 1997 .

[15]  Hong Xiao,et al.  Three dimensional elasto-plastic contact boundary element analysis for rolling with consideration of friction , 1998 .

[16]  J. Chern,et al.  Modification of the DDA method for rigid block problems , 1998 .

[17]  Chen Wanji,et al.  Smoothing Newton method for solving two‐ and three‐dimensional frictional contact problems , 1998 .

[18]  Jonathan D. Bray,et al.  CAPTURING NONSPHERICAL SHAPE OF GRANULAR MEDIA WITH DISK CLUSTERS , 1999 .

[19]  Yong-Il Kim,et al.  Modeling the effect of water, excavation sequence and rock reinforcement with discontinuous deformation analysis , 1999 .

[20]  Chen Wanji,et al.  A new non-smooth model for three dimensional frictional contact problems , 2000 .

[21]  Antonio Munjiza,et al.  Penalty function method for combined finite–discrete element systems comprising large number of separate bodies , 2000 .

[22]  Y. H. Hatzor,et al.  The validity of dynamic block displacement prediction using DDA , 2001 .

[23]  Lanru Jing,et al.  Modeling of fluid flow and solid deformation for fractured rocks with discontinuous deformation analysis (DDA) method , 2001 .

[24]  David M. Doolin,et al.  Displacement Accuracy of Discontinuous Deformation Analysis Method Applied to Sliding Block , 2002 .