A Review of the Influence of Microscopic Characteristics on the Progressively Brittle Failure of Foliated Rocks Subjected to Compression Loading

Geotechnical and Geological Engineering - The progressively brittle failure of foliated rocks involving the initiation, propagation, and aggregation modes is closely related to the microscopic...     

Numerical Simulation Study on Crack Evolution Characteristics of Coal Specimen Subjected to Conventional Compression Loading

Geotechnical and Geological Engineering - The paper represents a simulation investigation about the crack evolution and acoustic emission characteristics of coal specimen subjected to conventional...     

Dynamic Characteristics of Granite Subjected to Intermediate Loading Rate

Summary A large diameter split Hopkinson pressure bar (SHPB) has been developed. This equipment is briefly described, together with a shaped striker that initiates a half-sine incident waveform to obtain the complete stress-strain relationship of the Bukit Timah granite at medium strain rate. Good constant strain rate was derived, and the dynamic complete stress-strain curves and energy absorption of the granite were measured at a strain rate between 20 and 60 per second. Repeated impact between 60–90% of the static strength of the granite was also conducted. Results from the tests show that the cumulative damage of the granite depends on the peak stress of the dynamic loads with a fixed duration. The dynamic fracture strength of the granite loaded at medium strain rate is directly proportional to the cube root of the strain rate. For the granite loaded at this strain rate, Youngs modulus is unchanged. Energy absorption of the samples loading to fragmentation determined its fragmented size distribution. At high strain rate, the rock possesses large energy absorption and the particle size of the fragments is much smaller.     

Fracture Toughness and Fracture Roughness in Anisotropic Granitic Rocks

In this paper we present an experimental approach aimed at assessing the correlation between fracture toughness (K _IC) and fracture roughness of two granitic rocks (Barre and Stanstead granites) exhibiting significant fracture toughness anisotropy. Roughness values have been estimated for fractured surfaces obtained from Chevron Cracked Notch Brazilian Disc samples failed under mode I along three orthogonal planes with respect to their microstructural fabrics. There exists a clear correlation between roughness and toughness within each rock examined along the three planes. Specific orientation of micro-crack alignment could result in preferred out-of-plane propagation of the test-crack irrespective of grain-size distribution. These experimental observations reinforce the hypothesis of the existence of a link among pre-existing petrofabric anisotropy, fracture toughness, fracture roughness, and the evolution and extent of the associated induced fractures within the process zone of granitic rocks along specific directions. This study also highlights the need for employment of pre-failure and advanced post-failure diagnostic techniques in quantifying these inter-relationships.     

Numerical Modeling of Pore Pressure Influence on Fracture Evolution in Brittle Heterogeneous Rocks

Rock is a heterogeneous geological material. When rock is subjected to internal hydraulic pressure and external mechanical loading, the fluid flow properties will be altered by closing, opening, or other interaction of pre-existing weaknesses or by induced new fractures. Meanwhile, the pore pressure can influence the fracture behavior on both a local and global scale. A finite element model that can consider the coupled effects of seepage, damage and stress field in heterogeneous rock is described. First, two series of numerical tests in relatively homogeneous and heterogeneous rocks were performed to investigate the influence of pore pressure magnitude and gradient on initiation and propagation of tensile fractures. Second, to examine the initiation of hydraulic fractures and their subsequent propagation, a series of numerical simulations of the behavior of two injection holes inside a saturated rock mass are carried out. The rock is subjected to different initial in situ stress ratios and to an internal injection (pore) pressure at the two injection holes. Numerically, simulated results indicate that tensile fracture is strongly influenced by both pore pressure magnitude and pore pressure gradient. In addition, the heterogeneity of rock, the initial in situ stress ratio (K), the distance between two injection holes, and the difference of the pore pressure in the two injection holes all play important roles in the initiation and propagation of hydraulic fractures. At relatively close spacing and when the two principal stresses are of similar magnitude, the proximity of adjacent injection holes can cause fracturing to occur in a direction perpendicular to the maximum principal stress.     

Dynamic Analysis of Tunnel in Soil Subjected to Internal Blast Loading

The present work deals with the three-dimensional nonlinear finite element (FE) analyses of the tunnel in soil subjected to internal blast loading. The analyses are performed using the coupled Eulerian–Lagrangian analysis tool in FE software Abaqus/Explicit. The soil and reinforced concrete lining are modeled using the Lagrangian elements. The explosive Trinitrotoluene (TNT) is modeled using the Eulerian elements. The stress–strain response of soil, concrete, and reinforcement are simulated using strain rate dependent Drucker–Prager plasticity, concrete damaged plasticity and Johnson–Cook (J–C) plasticity models, respectively. The pressure–volume relationship of the TNT explosive is simulated using the Jones-Wilkins-Lee equation of state. Parametric sensitivity studies have been performed for different (1) tunnel lining thicknesses, (2) explosive charge weights and (3) angles of internal friction of soil. It is observed from the results that blast induced pressure on the tunnel lining increases with the increase in charge weight. Both the lining and the surrounding soil undergo significant deformation. The deformation of the tunnel lining increases with increasing charge weight and decreases with increasing lining thickness and increasing the angle of internal friction of soil. Blast-induced velocity in soil attenuates with increasing distance from the source of the blast.     

Application of Stochastic Fracture Network with Numerical Fluid Flow Simulations to Groundwater Flow Modeling in Fractured Rocks

ＩＮＴＲＯＤＵＣＴＩＯＮＧｒｏｕｎｄｗａｔｅｒｏｒｆｌｕｉｄｆｌｏｗｍｏｄｅｌｉｎｇｉｎｆｒａｃｔｕｒｅｄｒｏｃｋｓｉｓａｃｏｍｐｌｉｃａｔｅｄｔｈｅｏｒｅｔｉｃａｌａｎｄａｐｐｌｉｅｄｔｏｐｉｃ.Ｂｏｔｈｔｈｅｏｒｅｔｉｃａｌｌｙａｎｄｏｐｅｒａｔｉｏｎａｌｌｙ ,ｉｔｉｓｉｍｐｏｒｔａｎｔｉｎｍａｎｙｆｉｅｌｄｓｓｕｃｈａｓｇｅｏｌｏｇｉｃａｌａｎｄｈｙｄｒｏｇｅｏｌｏｇｉｃａｌｅｎｇｉｎｅｅｒｉｎｇ ,ｅｎｖｉｒｏｎｍｅｎｔａｌｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｅｔｒｏｌｅｕｍｅｎｇｉｎｅｅｒｉｎ…     

Linear Elastic and Cohesive Fracture Analysis to Model Hydraulic Fracture in Brittle and Ductile Rocks

Hydraulic fracturing technology is being widely used within the oil and gas industry for both waste injection and unconventional gas production wells. It is essential to predict the behavior of hydraulic fractures accurately based on understanding the fundamental mechanism(s). The prevailing approach for hydraulic fracture modeling continues to rely on computational methods based on Linear Elastic Fracture Mechanics (LEFM). Generally, these methods give reasonable predictions for hard rock hydraulic fracture processes, but still have inherent limitations, especially when fluid injection is performed in soft rock/sand or other non-conventional formations. These methods typically give very conservative predictions on fracture geometry and inaccurate estimation of required fracture pressure. One of the reasons the LEFM-based methods fail to give accurate predictions for these materials is that the fracture process zone ahead of the crack tip and softening effect should not be neglected in ductile rock fracture analysis. A 3D pore pressure cohesive zone model has been developed and applied to predict hydraulic fracturing under fluid injection. The cohesive zone method is a numerical tool developed to model crack initiation and growth in quasi-brittle materials considering the material softening effect. The pore pressure cohesive zone model has been applied to investigate the hydraulic fracture with different rock properties. The hydraulic fracture predictions of a three-layer water injection case have been compared using the pore pressure cohesive zone model with revised parameters, LEFM-based pseudo 3D model, a Perkins-Kern–Nordgren (PKN) model, and an analytical solution. Based on the size of the fracture process zone and its effect on crack extension in ductile rock, the fundamental mechanical difference of LEFM and cohesive fracture mechanics-based methods is discussed. An effective fracture toughness method has been proposed to consider the fracture process zone effect on the ductile rock fracture.     

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.     

Petrofabric Analyses of Rhum and Skaergaard Layered Rocks

Petrofabric analyses of layered rocks from Rhum have revealeda preferred orientation for felspar in the allivalites and forolivine in the peridotites; a regional petrofabric map of felsparorientation contains a radial pattern which suggests the presenceof convection currents during crystal settling. An orientedspecimen of Skaergaard ferrogabbro from the margin of a troughband has allowed comparison to be made between a known magmacurrent direction and the preferred orientation of felspar,olivine, clinopyroxene, and apatite crystals in the rock.     

Dynamic Study on Fracture Problems in Viscoelastic Sedimentary Rocks Using the Numerical Manifold Method

The viscoelastic deformation behavior of a sedimentary rock under different loading rates is numerically modeled and investigated by the numerical manifold method (NMM). By incorporating a modified 3-element viscoelastic constitutive mode in the NMM, crack initiation and propagation criteria, and crack identification and evolution techniques, the effects of the loading rates on the cracking behavior of a sedimentary rock, such as crack open displacement, crack sliding displacement, crack initiation, crack propagation and final failure mode, are successfully modeled. The numerical results reveal that under a high loading rate (>1,000 MPa/s), due to the viscoelastic property of the sedimentary rock, not only the structural behavior deviates from that of elastic model, but also different cracking processes and final failure modes are obtained.     

Fatigue Behavior of Granite Subjected to Cyclic Loading Under Triaxial Compression Condition

A series of laboratory tests were performed to examine the fatigue behavior of granite subjected to cyclic loading under triaxial compression condition. In these tests, the influences of volumetric change and residual strain on the deformation modulus of granite under triaxial cyclic compression were investigated. It is shown that the fatigue behavior of granite varies with the tendency for volumetric change in triaxial cyclic compression tests. In the stress–strain space, there are three domains for fatigue behavior of rock subjected to cyclic loading, namely the volumetric compaction, volumetric dilation with strain-hardening behavior, and volumetric dilation with strain-softening behavior domains. In the different domains, the microscopic mechanisms for rock deformation are different. It was also found that the stress level corresponding to the transition from volumetric compaction to volumetric dilation could be considered as the threshold for fatigue failure. The potential of fatigue deformation was compared with that of plastic deformation. The comparison shows that rocks exhibit higher resistances to volumetric deformation under cyclic loading than under plastic loading. The influence of residual strain on the fatigue behavior of rock was also investigated. It was found that the axial residual strain could be a better option to describe the fatigue behavior of rock than the loading cycle number. A constitutive model for the fatigue behavior of rock subjected to cyclic loading is proposed according to the test results and discussion. In the model, the axial residual strain is considered as an internal state variable. The influences of confining pressure and peak deviatoric stress on the deformation modulus are considered in a term named the equivalent stress. Comparison of test results with model predictions shows that the proposed model is capable of describing the prepeak fatigue behavior of rock subjected to cyclic loading.     

Mixed-Mode Fracturing of Rocks Under Static and Cyclic Loading

Static diametrical compression tests conducted on inclined cracked chevron notched Brazilian disc (CCNBD) Brisbane tuff specimens showed that the notched cracks at the centre of the specimens opened (Mode I) up to 30° crack inclination angle (β), whereas crack closure (Mode II) started for β > 33°, and closure became more pronounced at even higher β of 45° and 70°. Both the experimental and numerical results showed that the crack initiation angle (θ) was a function of the β. Scanning electron microscope (SEM) images showed that fatigue damage on cyclic loading of Brisbane tuff is strongly influenced by the failure of the matrix due to both intergranular and transgranular fracturing.