Compared with two-dimensional (2D) contact cases, three-dimensional (3D) contact problems are more complicated in both contact geometry representation and contact detection. The entrance block theory proposed by Shi lays the foundation for solving contact problems in a uniform format. This paper presents a new contact detection approach for 3D polyhedral blocks with the use of entrance block theory. In this approach, contacts between 3D polyhedra were classified into three categories: finite-contact covers (face to face, edge to face, and parallel edge to edge), single-contact covers (vertex to face, crossing edge to edge, vertex to edge, and vertex to vertex), and basic-contact covers (vertex to face and crossing edge to edge). A finite-contact cover consists of several single-contact covers, and all single-contact covers can be identified as basic-contact covers. Thus, the contact detection process can be fulfilled in three phases: a rough search phase to form a potential contacting block list f...

For the stability analysis of rock slope, it is very critical to obtain the spatial geometric characteristics of the structural surfaces of the rock mass accurately and effectively. As for a high-steep rock slope of an iron ore mine, in order to solve the problems of inefficiency and high risk of traditional manual geological survey, the geological survey and stability evaluation of the slope was carried out by adopting unmanned aerial vehicle digital photogrammetry (UAV-DP) technology. Firstly, a large number of high-resolution images of the slope were obtained by UAV-DP. Then, the structure from motion (SFM) method was used to construct the fine 3D point cloud model of the slope, which was subjected to coplanarity detection and K-means clustering for identifying the structural surfaces. Finally, the stability and failure model of the slope cut by the structural surfaces are analyzed by using the stereo-projection and discrete element methods. The research results show that the error between UAV-DP and manual measurement is within the acceptable range, which demonstrates the reliability of UAV-DP used in the geological investigation. Furthermore, the stability state and failure model of the slope is also consistent well with the field observation.

In order to avoid the risk of tailing pond failures and to minimize the post-failure losses, it is necessary to analyze the current operation status of tailings ponds, to explore the evolution law of their failure process, to grasp their post-failure impact range, and to propose corresponding effective prevention and control measures. Based on a tailings pond in China, this paper establishes a 1:200 scale indoor model to explore the evolution law of post-failure tailings discharge in a tailings pond under flooded roof conditions; secondly, the finite element difference method and smooth particle fluid dynamics are combined to compare and analyze the post-failure impact area and to delineate the risk prevention and control area. The results of the study show that, during the dam break, the lower tailing sand in the breach is the first to slip, and after forming a steep can, the upper tailing sand in the steep can is pulled to slip, so that the erosion trench mainly develops vertically first, and then laterally. The velocity of the discharged tailing sand will quickly reach its peak value in a short period of time and then decrease to the creeping stage; the front edge of the sand flow is the first to stop moving, and the trailing edge of the tailing sand accumulation depth continues to increase until the end of the dam failure, at which point the initial bottom dam area of the discharge tailing sand flow velocity is the largest. The further the tailings are released from the initial dam, the smaller the accumulation depth and the larger the particle size, and the larger the elevation of the foundation in the same section, the smaller the accumulation depth and the larger the particle size; further, the presence of blocking materials will increase the local tailings accumulation depth. Based on the maximum flow velocity of the discharged tailings and the accumulation depth, the risk area downstream of the tailings pond is divided, so that relocation measures can be formulated. The test results can provide an important reference for the operation and management of similar tailings ponds.

Majiagou landslide, a major ancient landslide in Three Gorges Reservoir region, is located in the high earthquake area of southwest China. The 2013 Badong earthquake caused an obvious deformation of landslide monitored by the sliding inclinometer. A strong earthquake may induce the reactivation of ancient landslide. So, it is necessary to research the seismic dynamic response of Majiagou landslide. For this purpose, discontinuous deformation analysis (DDA), improved by introducing the artificial joint and viscous boundary, is applied in this study. The displacements at monitoring points caused by Badong earthquake are calculated and compared with the field data, verifying the numerical method and model. Further, a strong earthquake with the peak acceleration of 1 g is assumed to act on the landside, the initiation and evolution process of landslide is simulated, and the movement features of landslide are discussed. The dynamic failure of landslide and the local amplification of seismic wave can be embodied, indicating that the improved DDA provides an alternative approach for analyzing the seismic dynamic response of jointed rock.

Evaluation of blasting‐induced rock damage and fragmentation is very important for safety control of construction in the jointed rock mass. The discontinuous numerical models are commonly applied due to the advantages in modeling fragmentation and treating discontinuities. In this paper, the rock fracturing algorithm and rate dependent strength law are incorporated into the discontinuous deformation analysis (DDA) to study the wave propagation and rock fragmentation phenomena in dynamic problems. By reproducing Hopkinson pressure bar tests under different loading rates, the improved method is validated to be capable of solving dynamic failure problem. Finally, taking the Xiucun tunnel as an example, its failure process under the action of the explosive wave is simulated, and some failure features are captured. In addition, the explosion wave propagation and its induced particle vibration in surrounding rock are studied. The reasonable simulation results indicate that the modified DDA method can effectively model the stress wave propagation and joint growth process in the jointed rock under blasting load.

In view of the fact that the field monitoring information can more intuitively and accurately reflect the stability state of slopes, this paper takes the warning threshold of slope deformation rate monitoring as the slope stability evaluation standard, and puts forward a method for solving the safety coefficient of rocky slopes. The discrete element method 3DEC, which is suitable for rocky slopes, is selected as the numerical calculation tool, and the convergence criterion of its strength reduction method is modified to the slope deformation rate threshold, and the method is realized by the bifurcation method through its built-in FISH programming language. The results of the classical case show that, by selecting the slope deformation rate threshold in the appropriate interval, the results of this paper's method are very close to those of the finite unit stress method and the limit equilibrium method, which verifies the reliability of this paper's method. Further, the method of this paper is applied to an open-pit mine slope project, based on the slope deformation on-site monitoring data, through the time series prediction method to determine the slope deformation rate warning threshold, using the threshold as an evaluation criterion to solve the slope's coefficient of safety, the calculation results show that: the slope's coefficient of safety in natural working conditions is 1.086, basically stable, suffered from continuous rainfall, the slope body is gradually saturated, the slope's coefficient of safety is reduced to 0.987, the slope's safety is reduced to 0.987. After continuous rainfall and gradual saturation of the slope, the coefficient of safety decreased to 0.987, destabilizing and destroying, which is consistent with the site conditions.