Numerical and experimental investigation of rock breaking method under free surface by TBM disc cutter

To study the rock breaking method under the free surface induced by disc cutter,the rock breaking simulations were first conducted based on the discrete element method,and the dynamic process of rock breaking under the free surface was studied including stressed zone,crush zone,crack initiation and propagation.Then the crack propagation conditions,specific energy,etc.under different free surface distance(S)were also investigated combined with linear cutting experiments.The results show that the rock breaking process under the free surface induced by disc cutter is dominated by tension failure mode.There exists a critical S to promote crack propagation to free surface effectively.And this rock breaking method can improve the rock breaking force and breaking efficiency significantly when proper.

Distinct element modelling of fracture plan control in continuum and jointed rock mass in presplitting method of surface mining

Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surface blast design. The purpose of presplitting is to form a fracture plane across which the radial cracks from the production blast cannot travel. The purpose of this study is to investigate of effect of presplitting on the generation of a smooth wall in continuum and jointed rock mass. The 2D distinct element code was used to simulate the presplitting in a rock slope. The blast load history as a function of time was applied to the inner wall of each blasthole. Important parameters that were considered in the analysis were stress tensor and fracturing pattern. The blast loading magnitude and blasthole spacing and jointing pattern were found to be very significant in the final results.

Analysis profile of the fully grouted rock bolt in jointed rock using analytical and numerical methods

The purpose of this study was to investigate the effect of bolt profile on load transfer mechanism of fully grouted bolts in jointed rocks using analytical and numerical methods. Based on the analytical method with development of methods, a new model is presented. To validate the analytical model, five different profiles modeled by ANSYS software. The profile of rock bolts T3 and T4with load transfer capacity,respectively 180 and 195 kN in the jointed rocks was selected as the optimum profiles. Finally, the selected profiles were examined in Tabas Coal Mine. FLAC analysis indicates that patterns 6+7 with2 NO flexi bolt 4 m better than other patterns within the faulted zone.

Modeling rock failure using the numerical manifold method followed bythe discontinuous deformation analysis

A complete rock failure process usually involves opening/sliding of preexisting discontinuities as well as frac- turing in intact rock bridges to form persistent failure sur- faces and subsequent motions of the generated rock blocks. The recently developed numerical manifold method （NMM） has potential for modelling such a complete failure process. However, the NMM suffers one limitation, i.e., unexpected material domain area change occurs in rotation modelling. This problem can not be easily solved because the rigid body rotation is not represented explicitly in the NMM. The discontinuous deformation analysis （DDA） is specially de- veloped for modelling discrete block systems. The rotation- induced material area change in the DDA modelling can be avoided conveniently because the rigid body rotation is represented in an explicit form. In this paper, a transition technique is proposed and implemented to convert a NMMmodelling to a DDA modelling so as to simulate a complete rock failure process entirely by means of the two methods, in which the NMM is adopted to model the early fracturing as well as the transition from continua to discontinua, while the DDA is adopted to model the subsequent motion of the generated rock blocks. Such a numerical approach also im- proves the simulation efficiency greatly as compared with a complete NMM modelling approach. The fracturing of a rock slab with pre-existing non-persistent joints located on a slope crest and the induced rockfall process are simulated. The validity of the modelling transition from the NMM to the DDA is verified and the applicability of the proposed nu- merical approach is investigated.

A review of discrete modeling techniques for fracturing processes in discontinuous rock masses

The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.

Numerical evaluation of strength and deformability of fractured rocks

Knowledge of the strength and deformability of fractured rocks is important for design, construction and stability evaluation of slopes, foundations and underground excavations in civil and mining engineering. However, laboratory tests of intact rock samples cannot provide information about the strength and deformation behaviors of fractured rock masses that include many fractures of varying sizes, orientations and locations. On the other hand, large-scale in situ tests of fractured rock masses are economically costly and often not practical in reality at present. Therefore, numerical modeling becomes necessary. Numerical predicting using discrete element methods(DEM) is a suitable approach for such modeling because of their advantages of explicit representations of both fractures system geometry and their constitutive behaviors of fractures, besides that of intact rock matrix. In this study, to generically determine the compressive strength of fractured rock masses, a series of numerical experiments were performed on two-dimensional discrete fracture network models based on the realistic geometrical and mechanical data of fracture systems from feld mapping. We used the UDEC code and a numerical servo-controlled program for controlling the progressive compressive loading process to avoid sudden violent failure of the models. The two loading conditions applied are similar to the standard laboratory testing for intact rock samples in order to check possible differences caused by such loading conditions. Numerical results show that the strength of fractured rocks increases with the increasing confning pressure, and that deformation behavior of fractured rocks follows elasto-plastic model with a trend of strain hardening. The stresses and strains obtained from these numerical experiments were used to ft the well-known Mohr-Coulomb(MC) and Hoek-Brown(H-B) failure criteria, represented by equivalent material properties defning these two criteria. The results show that both criteria can provide fair estimates of the compressive strengths for all tested numerical models. Parameters of the elastic deformability of fractured models during elastic deformation stages were also evaluated, and represented as equivalent Young’s modulus and Poisson’s ratio as functions of lateral confning pressure. It is the frst time that such systematic numerical predicting for strength of fractured rocks was performed considering different loading conditions, with important fndings for different behaviors of fractured rock masses, compared with testing intact rock samples under similar loading conditions.

Fluid–solid coupling analysis of rock pillar stability for concealed karst cave ahead of a roadway based on catastrophic theory

In order to study the mechanism of water inrush from a concealed, confined karst cave, we established a fluid–solid coupling model of water inrush from a concealed karst cave ahead of a roadway and a strength reduction method in a rock pillar for preventing water inrush based on catastrophic theory. Fluid–solid coupling effects and safety margins in a rock pillar were studied. Analysis shows that rock pillar instability, exerted by disturbance stress and seepage stress, is the process of rock pillar catastrophic destabilization induced by nonlinear extension of plastic zones in the rock pillar. Seepage flow emerges in the rock pillar for preventing water inrush, accompanied by mechanical instability of the rock pillar. Taking the accident of a confined karst cave water-inrush of Qiyi Mine as an example, by studying the safety factor of the rock pillar and the relationship between karst cave water pressure and thickness of the rock pillar,it is proposed that rock pillar thickness with a safety factor equal to 1.5 is regarded as the calculated safety thickness of the rock pillar, which should be equal to the sum of the blasthole depth, blasting disturbance depth and the calculated safety thickness of the rock pillar. The cause of the karst water inrush at Qiyi Mine is that the rock pillar was so small that it did not possess a safety margin. Combining fluid–solid coupling theory, catastrophic theory and strength reduction method to study the nonlinear mechanical response of complicated rock engineering, new avenues for quantitative analysis of rock engineering stability evaluation should be forthcoming.

Determination and applications of rock quality designation(RQD)

Characterization of rock masses and evaluation of their mechanical properties are important and challenging tasks in rock mechanics and rock engineering. Since in many cases rock quality designation （RQD） is the only rock mass classification index available, this paper outlines the key aspects on determination of RQD and evaluates the empirical methods based on RQD for determining the deformation modulus and unconfined compressive strength of rock masses. First, various methods for determining RQD are presented and the effects of different factors on determination of RQD are highlighted. Then, the empirical methods based on RQD for determining the deformation modulus and unconfined compressive strength of rock masses are briefly reviewed. Finally, the empirical methods based on RQD are used to determine the deformation modulus and unconfined compressive strength of rock masses at five different sites including 13 cases, and the results are compared with those obtained by other empirical methods based on rock mass classification indices such as rock mass rating （RMR）, Q-system （Q） and geological strength index （GSI）. It is shown that the empirical methods based on RQD tend to give deformation modulus values close to the lower bound （conservative） and unconfined compressive strength values in the middle of the corresponding values from different empirical methods based on RMR, Q and GSI. The empirical methods based on RQD provide a convenient way for estimating the mechanical properties of rock masses but, whenever possible, they should be used together with other empirical methods based on RMR, Qand GSI.

Characterization of transient groundwater flow through a high arch dam foundation during reservoir impounding

Numerous deep underground projects have been designed and constructed in China, which are beyond the current specifications in terms of scale and construction difficulty. The severe failure problems induced by high in situ stress, such as rockburst, spalling, damage of deep surrounding rocks, and timedependent damage, were observed during construction of these projects. To address these problems, the dynamic design method for deep hard rock tunnels is proposed based on the disintegration process of surrounding rocks using associated dynamic control theories and technologies. Seven steps are basically employed：（i） determination of design objective,（ii） characteristics of site, rock mass and project, and identification of constraint conditions,（iii） selection or development of global design strategy,（iv）determination of modeling method and software,（v） preliminary design,（vi） comprehensive integrated method and dynamic feedback analysis, and（vii） final design. This dynamic method was applied to the construction of the headrace tunnels at Jinping II hydropower station. The key technical issues encountered during the construction of deep hard rock tunnels, such as in situ stress distribution along the tunnels, mechanical properties and constitutive model of deep hard rocks, determination of mechanical parameters of surrounding rocks, stability evaluation of surrounding rocks, and optimization design of rock support and lining, have been adequately addressed. The proposed method and its application can provide guidance for deep underground projects characterized with similar geological conditions.

Geotechnical investigations and remediation design for failure of tunnel portal section： a case study in northern Turkey

Mass movements are very common problems in the eastern Black Sea region of Turkey due to its climate conditions, geological, and geomorphological characteristics. High slope angle, weathering, dense rainfalls, and anthropogenic impacts are generally reported as the most important triggering factors in the region. Following the portal slope excavations in the entrance section of Cankurtaran tunnel, located in the region, where the highly weathered andesitic tuff crops out, a circular toe failure occurred. The main target of the present study is to investigate the causes and occurrence mechanism of this failure and to determine the feasible remedial measures against it using finite element method(FEM) in four stages. These stages are slope stability analyses for pre-and postexcavation cases, and remediation design assessments for slope and tunnel. The results of the FEM-SSR analyses indicated that the insufficient initial support design and weathering of the andesitic tuffs are the main factors that caused the portal failure. After installing a rock retaining wall with jet grout columns and reinforced slope benching applications, the factor of safety increased from 0.83 to 2.80. In addition toslope stability evaluation, the Rock Mass Rating(RMR), Rock Mass Quality(Q) and New Austrian Tunneling Method(NATM) systems were also utilized as empirical methods to characterize the tunnel ground and to determine the tunnel support design. The performance of the suggested empirical support design, induced stress distributions and deformations were analyzed by means of numerical modelling. Finally, it was concluded that the recommended stabilization technique was essential for the dynamic long-term stability and prevents the effects of failure. Additionally, the FEM method gives useful and reasonably reliable results in evaluating the stability of cut slopes and tunnels excavated both in continuous and discontinuous rock masses.

Anisotropy of strength and deformability of fractured rocks

Anisotropy of the strength and deformation behaviors of fractured rock masses is a crucial issue for design and stability assessments of rock engineering structures, due mainly to the non-uniform and non- regular geometries of the fracture systems. However, no adequate efforts have been made to study this issue due to the current practical impossibility of laboratory tests with samples of large volumes con- taining many fractures, and the difficulty for controlling reliable initial and boundary conditions for large-scale in situ tests. Therefore, a reliable numerical predicting approach for evaluating anisotropy of fractured rock masses is needed. The objective of this study is to systematically investigate anisotropy of strength and deformability of fractured rocks, which has not been conducted in the past, using a nu- merical modeling method. A series of realistic two-dimensional （2D） discrete fracture network （DFN） models were established based on site investigation data, which were then loaded in different directions, using the code UDEC of discrete element method （DEM）, with changing confining pressures. Numerical results show that strength envelopes and elastic deformability parameters of tested numerical models are significantly anisotropic, and vary with changing axial loading and confining pressures. The results indicate that for design and safety assessments of rock engineering projects, the directional variations of strength and deformability of the fractured rock mass concerned must be treated properly with respect to the directions of in situ stresses. Traditional practice for simply positioning axial orientation of tunnels in association with principal stress directions only may not be adequate for safety requirements. Outstanding issues of the present study and su^zestions for future study are also oresented.

Fracture development around deep underground excavations: Insights from FDEM modelling

Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage development around underground excavations represents a key issue in several rock engineeringapplications, including tunnelling, mining, drilling, hydroelectric power generation, and the deepgeological disposal of nuclear waste. The goal of this paper is to show the effectiveness of a hybrid finitediscreteelement method （FDEM） code to simulate the fracturing mechanisms associated with theexcavation of underground openings in brittle rock formations. A brief review of the current state-of-theartmodelling approaches is initially provided, including the description of selecting continuum- anddiscontinuum-based techniques. Then, the influence of a number of factors, including mechanical and insitu stress anisotropy, as well as excavation geometry, on the simulated damage is analysed for threedifferent geomechanical scenarios. Firstly, the fracture nucleation and growth process under isotropicrock mass conditions is simulated for a circular shaft. Secondly, the influence of mechanical anisotropy onthe development of an excavation damaged zone （EDZ） around a tunnel excavated in a layered rockformation is considered. Finally, the interaction mechanisms between two large caverns of an undergroundhydroelectric power station are investigated, with particular emphasis on the rock mass responsesensitivity to the pillar width and excavation sequence. Overall, the numerical results indicate that FDEMsimulations can provide unique geomechanical insights in cases where an explicit consideration offracture and fragmentation processes is of paramount importance. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Safety assessment of waste rock dump built on existing tailings ponds

The construction of waste rock dumps on existing tailing ponds has been put into practice in China to save precious land resources. This work focuses on the safety assessment of the Daheishan molybdenum mine waste rock dump under construction on two adjoining tailings ponds. The consolidation of the tailings foundation and the filling quality of the waste rock are investigated by the transient electromagnetic method through detecting water-rich areas and loose packing areas, from which, the depth of phreatic line is also estimated. With such information and the material parameters, the numerical method based on shear strength reduction is applied to analyzing the overall stability of the waste rock dump and the tailings ponds over a number of typical cross sections under both current and designed conditions, where the complex geological profiles exposed by site investigation are considered. Through numerical experiments, the influence of soft lenses in the tailings and possible loose packing areas in the waste rock is examined. Although large displacements may develop due to the soft tailings foundation, the results show that the waste rock dump satisfies the safety requirements under both present and designed conditions.

DEM investigation of weathered rocks using a novel bond contact model

The distinct element method（DEM） incorporated with a novel bond contact model was applied in this paper to shed light on the microscopic physical origin of macroscopic behaviors of weathered rock, and to achieve the changing laws of microscopic parameters from observed decaying properties of rocks during weathering. The changing laws of macroscopic mechanical properties of typical rocks were summarized based on the existing research achievements. Parametric simulations were then conducted to analyze the relationships between macroscopic and microscopic parameters, and to derive the changing laws of microscopic parameters for the DEM model. Equipped with the microscopic weathering laws, a series of DEM simulations of basic laboratory tests on weathered rock samples was performed in comparison with analytical solutions. The results reveal that the relationships between macroscopic and microscopic parameters of rocks against the weathering period can be successfully attained by parametric simulations. In addition, weathering has a significant impact on both stressestrain relationship and failure pattern of rocks.