A semi-spring and semi-edge combined contact model in CDEM and its application to analysis of Jiweishan landslide

Abstract Continuum-based discrete element method (CDEM) is an explicit numerical method used for simulation of progressive failure of geological body. To improve the efficiency of contact detection and simplify the calculation steps for contact forces, semi-spring and semi-edge are introduced in calculation. Semi-spring is derived from block vertex, and formed by indenting the block vertex into each face (24 semi-springs for a hexahedral element). The formation process of semi-edge is the same as that of semi-spring (24 semi-edges for a hexahedral element). Based on the semi-springs and semi-edges, a new type of combined contact model is presented. According to this model, six contact types could be reduced to two, i.e. the semi-spring target face contact and semi-edge target edge contact. By the combined model, the contact force could be calculated directly (the information of contact type is not necessary), and the failure judgment could be executed in a straightforward way (each semi-spring and semi-edge own their characteristic areas). The algorithm has been successfully programmed in C++ program. Some simple numerical cases are presented to show the validity and accuracy of this model. Finally, the failure mode, sliding distance and critical friction angle of Jiweishan landslide are studied with the combined model.

[1]  Dawei Zhao,et al.  A fast contact detection algorithm for 3-D discrete element method , 2004 .

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

[3]  Ove Stephansson,et al.  Fundamentals of discrete element methods for rock engineering , 2007 .

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

[5]  S. H. Li,et al.  Analysis of Critical Excavation Depth for a Jointed Rock Slope Using a Face-to-Face Discrete Element Method , 2007 .

[6]  Ying Rao,et al.  A New Numerical Method For Dem –Block And Particle Model , 2004 .

[7]  Jian-Hong Wu,et al.  A new algorithm to identify contact patterns between convex blocks for three-dimensional discontinuous deformation analysis , 2008 .

[8]  P. Cundall,et al.  FORMULATION OF A THREE-DIMENSIONAL DISTINCT ELEMENT MODEL - PART II. MECHANICAL CALCULATIONS FOR MOTION AND INTERACTION OF A SYSTEM COMPOSED OF MANY POLYHEDRAL BLOCKS , 1988 .

[9]  Kianoosh Hatami,et al.  Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm , 2008 .

[10]  Manhong Zhao,et al.  Stochastic structural model of rock and soil aggregates by continuum-based discrete element method , 2005 .

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

[12]  C. H. Juang,et al.  Vertex‐to‐face contact searching algorithm for three‐dimensional frictionless contact problems , 2005 .

[13]  Jian-Hong Wu.  New edge-to-edge contact calculating algorithm in three-dimensional discrete numerical analysis , 2008, Adv. Eng. Softw..

[14]  W. S. Chen,et al.  New Approach to Determine Three-Dimensional Contacts in Blocks System: Penetration Edges Method , 2006 .

[15]  Lin Gao,et al.  Formulations of the three-dimensional discontinuous deformation analysis method , 2004 .

[16]  Tang Xiaojun Liu Junhua,et al.  Research on dynamic characteristics of multi-sensor system in the case of cross-sensitivity , 2005 .

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

[18]  Chun Feng,et al.  A GPU Accelerated Continuous-Based Discrete Element Method for Elastodynamics Analysis , 2011 .

[19]  李世海,et al.  A New Numerical Method for DEM-Block and Particle Model , 2004 .