Search

Overhanging rock instability is a serious geological problem in mountainous regions. The initiation and propagation of a controlling crack are the fundamental causes of overhanging rock failure. In the view of fracture mechanics, the instability of toppling overhanging rock can be equivalent to the effect of pure shear and pure bending moment. This paper analyzes the evolution of displacement-load curves and failure modes of overhanging rock models at different crack lengths, crack angles, and load distances based on fracture tests of the toppling overhanging rock model. Maximum circumferential stress theory is used as the fracture criterion of mixed mode I-II fractures to calculate the stress intensity factor (SIF) of the crack tip and the theoretical failure load of the overhanging rock model. The fracture mode mechanism is analyzed by the calculated SIF. We use the extended finite element method to simulate the crack propagation path and analyze the change process from equivalent stresses of the numerical model using propagation steps. Comparisons between the physical tests, theoretical calculations, and numerical simulations verify the rationality of the experimental results. [ABSTRACT FROM AUTHOR]

In this paper, the embedded discontinuity approach is applied in finite element modeling of rock in compression and tension. For this end, a rate-dependent constitutive model based on (strong) embedded displacement discontinuity model is developed to describe the mode I, mode II and mixed mode fracture of rock. The constitutive model describes the bulk material as linear elastic until reaching the elastic limit. Beyond the elastic limit, the rate-dependent exponential softening law governs the evolution of the displacement jump. Rock heterogeneity is incorporated in the present approach by random description of the mineral texture of rock. Moreover, initial microcrack population always present in natural rocks is accounted for as randomly-oriented embedded discontinuities. In the numerical examples, the model properties are extensively studied in uniaxial compression. The effect of loading rate and confining pressure is also tested in the 2D (plane strain) numerical simulations. These simulations demonstrate that the model captures the salient features of rock in confined compression and uniaxial tension. The developed method has the computational efficiency of continuum plasticity models. However, it also has the advantage, over these models, of accounting for the orientation of introduced microcracks. This feature is crucial with respect to the fracture behavior of rock in compression as shown in this paper. [ABSTRACT FROM AUTHOR]

Grouting plays a crucial role in dam safety. Due to the concealment of grouting activities, complexity of fracture distribution in rock masses and rheological properties of cement grout, it is difficult to analyze the effects of grouting. In this paper, a computational fluid dynamics (CFD) simulation approach of dam foundation grouting based on a 3D fracture network model is proposed. In this approach, the 3D fracture network model, which is based on an improved bootstrap sampling method and established by VisualGeo software, can provide a reliable and accurate geometric model for CFD simulation of dam foundation grouting. Based on the model, a CFD simulation is performed, in which the Papanastasiou regularized model is used to express the grout rheological properties, and the volume of fluid technique is utilized to capture the grout fronts. Two sets of tests are performed to verify the effectiveness of the Papanastasiou regularized model. When applying the CFD simulation approach for dam foundation grouting, three technical issues can be solved: (1) collapsing potential of the fracture samples, (2) inconsistencies in the geometric model in actual fractures under complex geological conditions, and (3) inappropriate method of characterizing the rheological properties of cement grout. The applicability of the proposed approach is demonstrated by an illustrative case study—a hydropower station dam foundation in southwestern China. [ABSTRACT FROM AUTHOR]

Fracture network modelling plays an important role in many application areas in which the behaviour of a rock mass is of interest. These areas include mining, civil, petroleum, water and environmental engineering and geothermal systems modelling. The aim is to model the fractured rock to assess fluid flow or the stability of rock blocks. One important step in fracture network modelling is to estimate the number of fractures and the properties of individual fractures such as their size and orientation. Due to the lack of data and the complexity of the problem, there are significant uncertainties associated with fracture network modelling in practice. Our primary interest is the modelling of fracture networks in geothermal systems and, in this paper, we propose a general stochastic approach to fracture network modelling for this application. We focus on using the seismic point cloud detected during the fracture stimulation of a hot dry rock reservoir to create an enhanced geothermal system; these seismic points are the conditioning data in the modelling process. The seismic points can be used to estimate the geographical extent of the reservoir, the amount of fracturing and the detailed geometries of fractures within the reservoir. The objective is to determine a fracture model from the conditioning data by minimizing the sum of the distances of the points from the fitted fracture model. Fractures are represented as line segments connecting two points in two-dimensional applications or as ellipses in three-dimensional (3D) cases. The novelty of our model is twofold: (1) it comprises a comprehensive fracture modification scheme based on simulated annealing and (2) it introduces new spatial approaches, a goodness-of-fit measure for the fitted fracture model, a measure for fracture similarity and a clustering technique for proposing a locally optimal solution for fracture parameters. We use a simulated dataset to demonstrate the application of the proposed approach followed by a real 3D case study of the Habanero reservoir in the Cooper Basin, Australia. [ABSTRACT FROM AUTHOR]

Different failure modes during fracture shearing have been introduced including dilation, sliding, asperity cut-off and degradation. Several laboratory studies have reported the complexity of these failure modes during shear tests performed under either constant normal load (CNL) or constant normal stiffness (CNS) conditions. This paper is concerned with the mechanical behaviour of synthetic fractures during direct shear tests using a modified shear cell and related numerical simulation studies. The modifications made to an existing true triaxial stress cell (TTSC) in order to use it for performing shear tests under CNL conditions are presented. The large loading capacity and the use of accurate hydraulic pumps capable of applying a constant shear velocity are the main elements of this cell. Synthetic mortar specimens with different fracture surface geometries are tested to study the failure modes, including fracture sliding, asperity degradation, and to understand failure during shearing. A bonded particle model of the direct shear test with the PFC2D particle flow code is used to mimic the tests performed. The results of a number of tests are presented and compared with PFC2D simulations. The satisfactory results obtained both qualitatively and quantitatively are discussed. [ABSTRACT FROM AUTHOR]

Different failure modes during fracture shearing have been introduced including normal dilation or sliding, asperity cut-off and degradation. Attempts have been made to study these mechanisms using analytical, experimental and numerical methods. However, the majority of the existing models simplify the problem, which leads to unrealistic results. With this in mind, the aim of this paper is to simulate the mechanical behaviour of synthetic and rock fracture profiles during direct shear tests by using the two-dimensional particle flow computer code PFC2D. Correlations between the simulated peak shear strength and the fracture roughness parameter D recently proposed by Rasouli and Harrison () are developed. Shear test simulations are carried out with PFC2D and the effects of the geometrical features as well as the model micro-properties on the fracture shear behaviour are studied. The shear strength and asperity degradation processes of synthetic profiles including triangular, sinusoidal and randomly generated profiles are analysed. Different failure modes including asperity sliding, cut-off, and asperity degradation are explicitly observed and compared with the available models. The D parameter is applied to the analysis of synthetic and rock fracture profiles. Accordingly, correlations are developed between D and the peak shear strength obtained from simulations and by using analytical solutions. The results are shown to be in good agreement with the basic understanding of rock fracture shear behaviour and asperity contact degradation. [ABSTRACT FROM AUTHOR]