Nuwen Xu

The Dynamic Evaluation of Rock Slope Stability Considering the Effects of Microseismic Damage

A state-of-the-art microseismic monitoring system has been implemented at the left bank slope of the Jinping first stage hydropower station since June 2009. The main objectives are to ensure slope safety under continuous excavation at the left slope, and, very recently, the safety of the concrete arch dam. The safety of the excavated slope is investigated through the development of fast and accurate real-time event location techniques aimed at assessing the evolution and migration of the seismic activity, as well as through the development of prediction capabilities for rock slope instability. Myriads of seismic events at the slope have been recorded by the microseismic monitoring system. Regions of damaged rock mass have been identified and delineated on the basis of the tempo-spatial distribution analysis of microseismic activity during the periods of excavation and consolidation grouting. However, how to effectively utilize the abundant microseismic data in order to quantify the stability of the slope remains a challenge. In this paper, a rock mass damage evolutional model based on microseismic data is proposed, combined with a 3D finite element method (FEM) model for feedback analysis of the left bank slope stability. The model elements with microseismic damage are interrogated and the deteriorated mechanical parameters determined accordingly. The relationship between microseismic activities induced by rock mass damage during slope instability, strength degradation, and dynamic instability of the slope are explored, and the slope stability is quantitatively evaluated. The results indicate that a constitutive relation considering microseismic damage is concordant with the simulation results and the influence of rock mass damage can be allowed for its feedback analysis of 3D slope stability. In addition, the safety coefficient of the rock slope considering microseismic damage is reduced by a value of 0.11, in comparison to the virgin rock slope model. Our results demonstrate that microseismic activity induced by construction disturbance only slightly affects the stability of the slope. The proposed feedback analysis technique provides a novel method for dynamically assessing rock slope stability and can be used to assess the slope stability of other similar rock slopes.

Numerical Investigation of Dynamic Rock Fracture Toughness Determination Using a Semi-Circular Bend Specimen in Split Hopkinson Pressure Bar Testing

The International Society for Rock Mechanics (ISRM) has suggested a notched semi-circular bend technique in split Hopkinson pressure bar (SHPB) testing to determine the dynamic mode I fracture toughness of rock. Due to the transient nature of dynamic loading and limited experimental techniques, the dynamic fracture process associated with energy partitions remains far from being fully understood. In this study, the dynamic fracturing of the notched semi-circular bend rock specimen in SHPB testing is numerically simulated for the first time by the discrete element method (DEM) and evaluated in both microlevel and energy points of view. The results confirm the validity of this DEM model to reproduce the dynamic fracturing and the feasibility to simultaneously measure key dynamic rock fracture parameters, including initiation fracture toughness, fracture energy, and propagation fracture toughness. In particular, the force equilibrium of the specimen can be effectively achieved by virtue of a ramped incident pulse, and the fracture onset in the vicinity of the crack tip is found to synchronize with the peak force, both of which guarantee the quasistatic data reduction method employed to determine the dynamic fracture toughness. Moreover, the energy partition analysis indicates that simplifications, including friction energy neglect, can cause an overestimation of the propagation fracture toughness, especially under a higher loading rate.

Numerical Assessment of the Progressive Rock Fracture Mechanism of Cracked Chevron Notched Brazilian Disc Specimens

The International Society of Rock Mechanics (ISRM) suggested cracked chevron notched Brazilian disc method falls into a major testing category of rock fracture toughness measurement by virtue of chevron notched rock samples. A straight through crack front during the whole fracturing process is assumed in the testing principle but is never assessed. In this study, the progressive rock fracture mechanism of cracked chevron notched Brazilian disc rock specimens is numerically simulated for the first time. Two representative sample types with distinct geometry of notch ligaments are modelled. The assumption of a straight through crack front for chevron notched fracture samples is critically assessed. The results show that not only the notch tip but also the saw-cut chevron notch cracks during the experiments. The straight through crack front assumption is never satisfied in the realistic rock fracture progress of chevron notched disc samples. In addition, the crack features prominent curved front, far from being straight. In contrast to the sample type with narrow notch ligament, the acoustic emission (AE) of the simulation on the sample with wide notch ligament depicts obvious biased fracturing of the prescript fracturing route of the notch. The numerically observed progressive fracture mechanism calls for more attention on how to accurately calibrate the critical dimensionless stress intensity factor for a better measurement of Mode I fracture toughness via chevron notched samples.

Static and Dynamic Flexural Strength Anisotropy of Barre Granite

Granite exhibits anisotropy due to pre-existing microcracks under tectonic loadings; and the mechanical property anisotropy such as flexural/tensile strength is vital to many rock engineering applications. In this paper, Barre Granite is studied to understand the flexural strength anisotropy under a wide range of loading rates using newly proposed semi-circular bend tests. Static tests are conducted with a MTS hydraulic servo-control testing machine and dynamic tests with a split Hopkinson pressure bar (SHPB) system. Six samples groups are fabricated with respect to the three principle directions of Barre granite. Pulse shaping technique is used in all dynamic SHPB tests to facilitate dynamic stress equilibrium. Finite element method is utilized to build up equations calculating the flexural tensile strength. For samples in the same orientation group, a loading rate dependence of the flexural tensile strength is observed. The measured flexural tensile strength is higher than the tensile strength measured using Brazilian disc method at given loading rate and this scenario has been rationalized using a non-local failure theory. The flexural tensile strength anisotropy features obvious dependence on the loading rates, the higher the loading rate, the less the anisotropy and this phenomenon may be explained considering the interaction of the preferentially oriented microcracks.

Fracture Toughness Determination of Cracked Chevron Notched Brazilian Disc Rock Specimen via Griffith Energy Criterion Incorporating Realistic Fracture Profiles

The cracked chevron notched Brazilian disc (CCNBD) specimen has been suggested by the International Society for Rock Mechanics to measure the mode I fracture toughness of rocks, and has been widely adopted in laboratory tests. Nevertheless, a certain discrepancy has been observed in results when compared with those derived from methods using straight through cracked specimens, which might be due to the fact that the fracture profiles of rock specimens cannot match the straight through crack front as assumed in the measuring principle. In this study, the progressive fracturing of the CCNBD specimen is numerically investigated using the discrete element method (DEM), aiming to evaluate the impact of the realistic cracking profiles on the mode I fracture toughness measurements. The obtained results validate the curved fracture fronts throughout the fracture process, as reported in the literature. The fracture toughness is subsequently determined via the proposed G-method originated from Griffith’s energy theory, in which the evolution of the realistic fracture profile as well as the accumulated fracture energy is quantified by DEM simulation. A comparison between the numerical tests and the experimental results derived from both the CCNBD and the semi-circular bend (SCB) specimens verifies that the G-method incorporating realistic fracture profiles can contribute to narrowing down the gap between the fracture toughness values measured via the CCNBD and the SCB method.

Numerical Observation of Three-Dimensional Wing Cracking of Cracked Chevron Notched Brazilian Disc Rock Specimen Subjected to Mixed Mode Loading

The cracked chevron notched Brazilian disc (CCNBD) specimen has been suggested by International Society for Rock Mechanics for measuring mode I fracture toughness of rocks. Subsequently, this specimen geometry has been widely extended to conduct mixed mode fracture tests on rocks as well. A straight through crack front during the fracturing process upon the root of the chevron notch is assumed in the testing principle, but has never been thoroughly evaluated before. In this study, for the first time, the progressive rock fracture mechanism of the CCNBD rock specimen under mixed mode loading is numerically simulated. Specimens under representative mixed mode loading angles are modelled; and the assumption of the straight through crack front growth is critically assessed. The results show that not only the notch tip but also the saw-cut chevron notch cracks during the experiments, yielding a prominent twisted front, far from being straight. The crack front never grows up to the root of the notch ligament and the straight through crack front assumption is never satisfied in the realistic rock fracture progress of this chevron notched specimen subjected to mixed mode loads. In contrast, the fracture progress features typical three-dimensional wing cracking towards the loading ends. The numerically observed progressive fracture mechanism reveals that the measuring principle of mixed mode fracture tests employing CCNBD specimens is significantly violated and the measures of both modes I and II fracture toughness are uncertain.

Experimental Investigation of the Influence of Joint Geometric Configurations on the Mechanical Properties of Intermittent Jointed Rock Models Under Cyclic Uniaxial Compression

Intermittent joints in rock mass are quite sensitive to cyclic loading conditions. Understanding the fatigue mechanical properties of jointed rocks is beneficial for rational design and stability analysis of rock engineering projects. This study experimentally investigated the influences of joint geometry (i.e., dip angle, persistency, density and spacing) on the fatigue mechanism of synthetic jointed rock models. Our results revealed that the stress–strain curve of jointed rock under cyclic loadings is dominated by its curve under monotonic uniaxial loadings; the terminal strain in fatigue curve is equal to the post-peak strain corresponding to the maximum cyclic stress in the monotonic stress–strain curve. The four joint geometrical parameters studied significantly affect the fatigue properties of jointed rocks, including the irreversible strains, the fatigue deformation modulus, the energy evolution, the damage variable and the crack coalescence patterns. The higher the values of the geometrical parameters, the lower the elastic energy stores in this jointed rock, the higher the fatigue damage accumulates in the first few cycles, and the lower the fatigue life. The elastic energy has certain storage limitation, at which the fatigue failure occurs. Two basic micro-cracks, i.e., tensile wing crack and shear crack, are observed in cyclic loading and unloading tests, which are controlled principally by joint dip angle and persistency. In general, shear cracks only occur in the jointed rock with higher dip angle or higher persistency, and the jointed rock is characterized by lower fatigue strength, larger damage variable and lower fatigue life.

Experimental and Numerical Study on the Cracked Chevron Notched Semi-Circular Bend Method for Characterizing the Mode I Fracture Toughness of Rocks

The cracked chevron notched semi-circular bending (CCNSCB) method for measuring the mode I fracture toughness of rocks combines the merits (e.g., avoidance of tedious pre-cracking of notch tips, ease of sample preparation and loading accommodation) of both methods suggested by the International Society for Rock Mechanics, which are the cracked chevron notched Brazilian disc (CCNBD) method and the notched semi-circular bend (NSCB) method. However, the limited availability of the critical dimensionless stress intensity factor (SIF) values severely hinders the widespread usage of the CCNSCB method. In this study, the critical SIFs are determined for a wide range of CCNSCB specimen geometries via three-dimensional finite element analysis. A relatively large support span in the three point bending configuration was considered because the fracture of the CCNSCB specimen in that situation is finely restricted in the notch ligament, which has been commonly assumed for mode I fracture toughness measurements using chevron notched rock specimens. Both CCNSCB and NSCB tests were conducted to measure the fracture toughness of two different rock types; for each rock type, the two methods produce similar toughness values. Given the reported experimental results, the CCNSCB method can be reliable for characterizing the mode I fracture toughness of rocks.

Microseismic Monitoring and Numerical Simulation of Rock Slope Failure

A numerical code with an elastic-brittle failure model has been developed to simulate the seismic activities in rock failure problems. The feasibility in principle of monitoring and prediction of rock failure was numerically simulated. The heterogeneity was considered to be the main reason for the existence of slope failure precursors. Seismic events could be observed in heterogeneous rock, whereas the homogeneous rocks showed an abrupt fracture mode without any early seismic precursors. The failure process of a slope was numerically investigated by using a gravity increase method (centrifugal loading method), and the application of the microseismic monitoring system in the slope was introduced. The numerical results showed that the fracture of the main faults caused the slope slide, and the microcracking caused by the heterogeneity in the faults prior to the landslide could be considered as the precursors of the slope failure, which were captured by the microseismic monitoring system. The microseismic monitoring technique was proved to be successful in predicting the failure in the slope, and the numerical results will be helpful in interpreting the microseismic monitoring results.

Microseismic Monitoring of the Left Bank Slope at the Baihetan Hydropower Station, China

Layered rock slopes are especially prone to collapse under excavation excitations due to the relatively weak strength characteristics of discontinuities such as joints, bedding planes and foliations (Liu et al. 2004, 2014). The stability analysis of layered rock slopes is of increased importance when subjected to continuous excavation. A systematic study and analysis of high rock slope stability has been performed in related studies (Stacey et al. 2004; Stead et al. 2006; Stead and Wolter 2015). However, approaches are rare that have already formed mature programmes and that have been effectively used in engineering practices. Compared with theoretical and numerical studies, field monitoring provides an actual and visual approach to the study of the deformation failure characteristics of rock slopes, especially for complicated cases (i.e., deep interlayer staggered zones and intraformational faulted zones) where theoretical or numerical solutions are difficult to obtain. Traditional measurement techniques, such as multiple-point extensometers, anchor stress gauges and global positioning system, are useful for monitoring the surface deformation. However, they are unrealistic to effectively monitor the microcracking activities in rock masses prior to the formation of macroscopic fractures on slope surface. In general, macro-fractures and deformation failure in the rock mass often lag behind the initiation, coalescences and propagation of microcracks. Thus, there must be an intrinsic correlation between micro-fractures (microseismicity) and rock slope instability (Xu et al. 2012). Although it is difficult to monitor the full-field stress, the response of the rock slope to stresses (i.e., microcracking) can be monitored. The variation of stress due to engineering excavation can be indirectly obtained through monitoring the microcracking evolution. The microseismic (MS) monitoring technique can help to achieve this goal (Kaiser 2009; Tang et al. 2011). More descriptions about the principle and applications of MS monitoring technique were detailed in related studies (Cai et al. 2001; Hirata et al. 2007; Hudyma and Potvin 2010; Trifu and Shumila 2010; Xu et al. 2014, 2016; Young et al. 2004).