Cracking Behavior and Acoustic Emission Characteristics of Rock Containing a Single Preexisting Flaw

The existence of flaws in brittle rocks or rock-like materials has an obvious influence on the material mechanical properties and cracking behavior of civil engineering projects. In this work, the two-dimensional particle flow code PFC2D was used to study the deformation and strength properties, failure processes, and acoustic emission (AE) characteristics of mudstone with a single preexisting flaw. First, the procedure to construct a parallel bond model is introduced. The Weibull distribution is used to reflect the mechanical heterogeneity of rocks. Then, the microscopic parameters used in PFC2D are calibrated to the macroproperties of mudstone obtained from laboratory tests under the uniaxial compression. The results indicate that the increases of the flaw inclination lead to the increasing uniaxial compressive strength and elastic modulus. In terms of microcrack evolution, the initiation, propagation, and coalescence of microcracks are closely related to the force chain. Specifically, an “X” shaped tension force chain concentrated area around the preexisting flaw is founded, which is the most prone area for microcracks to initiate. With an increase in flaw inclination, the b value of AE also shows an increasing trend. By incorporating the AE event numbers into a damage variable, this paper derives a constitutive model, which is verified by numerical results on brittle rocks with a single preexisting flaw under uniaxial compression.

Compressive Behavior of Oil Shale with Calcareous Concretion: Parametric Study

The non-uniformly distributed calcareous concretion among the oil shale in the Junggar basin of China has led to the difficulty in achieving the slope stability. This paper presents the numerical simulation of the behavior of oil shale with calcareous concretion via the Particle Flow Code (PFC2D) program based on the trial experimental test results. The critical parameters investigated in this research covered the size, distribution, strength, and number of the calcareous concretion. The following conclusions can be drawn based on the discussions and analysis: (1) the hard concretion always results in the high compressive strength of the specimen compared with that without concretion; (2) when the radius of the concretion size raised from 2.5 mm to 20 mm, the peak strength of tested specimens is approximately 50 MPa, whereas, the specimen with large concretion is much more ductile under compression; (3) the compressive behavior of tested specimens is similar even when the position of the concretion is variable; and (4) different from the specimens with only one concretion, these specimens contained two concretions featured with the double “X” failure mode. Meanwhile, the peak strength of the specimens with two hard concretions is about 2.5 times that of its counterparts with two soft concretions. The numerical simulation results are meaningful in guiding the design and analysis of the oil shale slope with the concretion.

Experimental and Numerical Investigation of Failure Characteristics and Mechanism of Granites with Different Joint Angles

The uniaxial compression tests were conducted on granite samples with different joint dip angles to more favorably explore the influences of the nonconsecutive joint on mechanical properties and deformation characteristics of the rock mass. The stress-strain curves, deformation and strength characteristics, and energy evolution process of the samples were analyzed. Numerical simulation using particle flow code (PFC) is employed to study the crack propagation process. The mode of jointed and fractured rock was investigated. The research results showed a significant reduction in both the peak strength and elastic modulus of jointed samples compared with intact ones: the peak strength and elastic modulus drop to the minimum at the joint dip angle of about 45°, especially for the peak strength, which takes up about 55% of the intact samples. The fractured samples’ total energy, elastic strain energy, and dissipated energy during the uniaxial compression drop significantly relative to intact samples. The proportion of the fracture modes varies with different joint dip angles, in which the ratio of shear cracks grows at first and then declines, with the highest balance at the dip angle of 45°. The damage stress’s sensitivity to the dip angle change is greater than that of the peak stress, with reduction amplitude more extensive than the latter.

Scenario Simulation of The Geohazard Dynamic Process of Large-Scale Landslides: A Case Study of The Xiaomojiu Landslide Along The Jinsha River

Large-scale landslides often cause severe damage due to their long run-out distances and having disaster chain effects. Scenario simulation has been adopted in the current work to analyze the Xiaomojiu landslide dynamic processes, such as sliding velocity, deposition characteristics, and flood outburst after a landslide-dam failure using Particle Flow Code (PFC-3D) which introduced the changeable friction coefficient and the HEC-RAS software. The landslide characteristics and topography data were obtained via field investigation, whereas high-resolution topographic data (0.17 m) was obtained using an Unmanned Aerial Vehicle (UAV). The results showed that: 1. The landslide presents a scallop shape with a length of 1566 m, a width ranging from 809~1124 m, and an area of 1.34×106 m2. The average thickness and volume of the sliding body is approximately 40 m, 5.1×107 m3. The InSAR deformation analysis showed that the Xiaomojiu landslide has a maximum annual displacement rate of 60 mm/y, and a maximum accumulation deformation of 180 mm since November 25, 2017. 2. From the landslide simulation results, the failure process of the Xiaomojiu landslide lasted for 65 s with a maximum velocity of 78.2 m/s. The deposited area is approximately 2023 m long, 900 m wide, with a maximum height of approximately 149 m. 3. After the landslide blocks the Jinsha River, a landslide-dammed lake with an elevation of 2940 m and a storage capacity of 4.13×109 m3 is formed. The maximum peak flow rate of the breach is 12051.7 m3/s, 43451.4 m3/s, 148635.6 m3/s, and 304544.7 m3/s for the landslide-dammed failure degrees of 15%, 25%, 50%, and 75%, respectively. These results provide a scientific reference for the risk analysis and mitigation of the landslide.

Simulating Crack Development and Failure Characteristic of Toppling Rock Slope under Seismic Loading on Lancang River in China

Toppling rock slopes, induced by rapid and continuous downcutting of Lancang River, are widely distributed in the mountainous area of southwest China. To investigate the instability mechanism of 1# toppling rock slope of Huangdeng Hydropower Station under seismic loading, particle flow code (PFC) is applied to simulate the dynamic response and failure mode. The study considers the particle characteristics of displacement, velocity, energy, and cracks. According to numerical results, the potential failure mechanism of toppling rock slope is identified: multisliding surfaces form at the interfaces between the highly and moderately toppled rock mass and between the highly/moderately and weak toppled-crept rock mass; intersecting faults cut rock mass at the toe, leading to shear-toppling deformation; tension cracks develop, penetrate, and coalesce in the weak toppled-crept rock mass, resulting in tension-toppling-bending deformation. During the 2 to 5 s of strong seismic intensity, crack increases sharply and energy of particles fluctuates greatly. The impacts of the amplitude of seismic loading and loading method in PFC are investigated. This study will provide a practically useful reference for seismic design of rock slopes.

Numerical Simulations of the Soil–Rock Mixture Mechanical Properties Considering the Influence of Rock Block Proportions by PFC2D

Soil–rock mixtures (S-RMs), as a kind of special engineering geological material, need to be studied because of the special structure and complex movement mechanism of their rock blocks, their physical and mechanical properties, and the factors underlying rock block movement in the process of their deformation and failure. In this paper, a series of discrete-element numerical models are constructed in particle flow code software (PFC2D). First, the random structure numerical models of S-RMs with different rock block proportions are established. Then, the parameters of the soil meso-structure are inversed by the biaxial simulation test, and a series of biaxial compressive tests are performed. The characteristics of stress and strain, deformation and failure, and rock block rotation and energy evolution are systematically investigated. The results show the following. (1) As the rock block proportion (confining pressure 0.5 MPa) increases, the peak strength of increases, the fluctuations of the post-peak become more obvious, and the dilatancy of the sample increases. (2) As the rock block proportion increases, the width of the shear band increases, the distribution of cracks becomes more complex and dispersed, and the range of the shear zone increases. (3) The number of rock blocks with rotation also increases significantly as rock block proportion increases, and the rotation angles are mostly between −5° and 5°. (4) The strain energy of S-RMs with different rock block proportions follows the same change rule as axial strain, showing a trend of first increasing and then decreasing, like the stress–strain curve.

Negative Stiffness, Incompressibility and Strain Localisation in Particulate Materials

In this paper, we consider two mechanisms capable of inducing strain localisation in particulate geomaterials in compression: the apparent negative stiffness and the incremental incompressibility caused by dilatancy. It is demonstrated that the apparent negative stiffness can be produced by the rotation of clusters of particles in the presence of compression. The clusters are formed by connecting the particles by the bonds that still remain intact in the process of bond breakage in compression. We developed a 2D isotropic model of incremental incompressibility showing that a single strain localisation zone is formed inclined at 45° to the direction of axial compressive loading. This mechanism of localisation was analysed through Particle Flow Code (PFC) 2D and 3D simulations. It is shown that, in the simulations, the peak stress (the onset of localisation) does correspond to the incremental Poisson’s ratio, reaching the critical values of 1 (in 2D) and 0.5 (in 3D).

The meso-fracturing mechanism of marble under unloading confining pressure paths and constant axial stress

Failure tests on marble during unloading confining-pressure under constant axial stress and simulations with the particle flow code were performed. The influence mechanism of the unloading rate of the confining pressure, initial unloading stress, and confining pressure on the failure characteristics of, and crack propagation in, marble was studied. By using the trial-and-error method, the conversion relationship between the unloading rates of confining pressures in laboratory tests and numerical simulations was ascertained. Micro-cracks formed in the unloading process of confining pressure are dominated by tension cracks, accompanied by shear cracks. The propagation of shear cracks lags that of tension cracks. As the confining pressure is increased, more cracks occur upon failure of the samples. The proportion of shear cracks increases while that of tension cracks decreases. The failure mode of samples undergoes a transition from shear-dominated failure to conjugated shear failure.