Abstract

This paper aims at the phenomenon of the fractured rock column in underground engineering which is prone to collapse when subjected to ground pressure. Uniaxial compression test and particle flow code (PFC2D) are used to analyze the influence of crack dip angle on the mechanical properties, crack propagation, and the failure mode of red sandstone. From the results, it is observed that the stress–strain curve of the precracked red sandstone can be divided into five stages, and there is a critical stress value in the stage of accelerated crack propagation and unstable fracture of the rock sample. Further, the peak strength, maximum strain, and the elastic modulus of the precracked red sandstone increase with the increase of crack dip angle, and the ultimate failure mode of rock sample changes from the “ladder” type failure to slope uneven failure. Furthermore, from PFC2D simulation, it is found that the tensile microcracks contribute more towards the failure of rock samples than the shear cracks. The contact force chain is very weak at the places where the precracks and macroshear planes are formed. This indicates that the original contact force is weakened due to particle fracture. Therefore, the bearing capacity of the precracked rock samples decreases with the increase in load. From the simulation results, it is found that the displacement at the shear plane of the rock sample is large, and the shear dilatation occurs. With the increase in load, the specimen falls off and is ejected. This is due to the weakening of the contact force between the internal particles. Thereafter, it fractures to produce microcracks, which gradually converge, thus providing a prerequisite for the transformation of elastic strain energy into kinetic energy.

1. Introduction

The underground pillars are under uniaxial compression. In general, the underground rock column has fracture development and some initial damage. Further damage occurs under loading. When the underground pillars are affected by dynamic loads such as blasting and rock burst, instability occurs. This leads to a mine disaster. The stability of the rock column is related to the stability and safety of the whole underground structure [1]. It is often difficult to make closed cracks in laboratories and closed cracks in practical engineering are also rare. Therefore, it is of general significance to study the problem of precast nonclosed cracks in rock samples that are under loading [25]. It is difficult to effectively monitor the microcracks in laboratory rock samples. The discrete element particle flow code (PFC) can be used to simulate the propagation of microcracks. It can simulate the problems from point to surface, step by step to final failure. It is more suitable for explaining the damage and failure evolution process of rock [68].

Some scientific research has been carried out on precracked rocks. Cao et al. observed three types of coalescence patterns in between the intermittent fractures [9]. Yang et al. studied the fracture characteristics of brittle sandstone with two noncoplanar fractures by using experiments and numerical simulation [10]. Huang et al. summarized the tensile failure mode and three-dimensional damage evolution characteristics of sandstone through acoustic emission signal [11]. Chen et al. studied the mechanical parameters and the failure characteristics of T-shaped cross fractures [12].

The geometric parameters and the filling properties of preexisting flaws have an important influence on the mechanical properties, crack propagation, and the deformation failure mode of rock. Extensive studies were conducted under uniaxial loading in order to investigate the damage evolution of a rock or rock-like materials containing prefabricated crack. Wong et al. studied the crack propagation mechanism, crack initiation angle, and the failure mode of the precracked specimens using the uniaxial compression test [1316]. Further, crack initiation and propagation for samples with a single nonstraight fissure were investigated [17, 18]. Zhang [19] simulated the cracking process in a rock-like material containing a single flaw using a numerical approach based on the particle element model. Lin et al. derived the relationship between the crack initiation angle of the open flaw tip and the dip angle of the prefabricated flaw under compression conditions [20]. Niu identified the change from microcracking to macrocracking using the acoustic emission system in the uniaxial compression tests [21, 22]. Zhao et al. analyzed the influence of the flaw dip angle and the filling material on the crack initiation stress and the peak strength [23, 24]. Moreover, the full-field strain evolution, fracture mode and the acoustic emission (AE) fractal characteristics of samples with a single preexisting flaw were comprehensively studied [2527].

From the literature, it is observed that the crack initiation, propagation, and mechanical properties of precracked rocks were extensively investigated. However, limited studies were conducted on the failure mode, damage effect, and the fracture characteristic stress of soft rocks with a single flaw using experiments and the mesoscopic simulation. Therefore, in this work, the red sandstone with a prefabricated single crack is considered for the study, and the uniaxial compression experiment is carried out. Simulation is carried out using PFC2D in order to explore its mechanical properties and study the crack propagation law from the macro- and microview. It is also studied to summarize the failure mode of instability in order to reveal the failure characteristics of precracked sandstone. Moreover, it has important reference significance for the evaluation of the stability of the underground rock column.

2. Experimental Details of the Uniaxial Compression Experiment on Dip Angle Effect of the Single-Cracked Red Sandstone

The materials, experimental details, and the experimental plan of the uniaxial compression experiment are briefly described below.

2.1. Experimental Materials

The material that is used for the experiments is red sandstone. For conducting the experiment, red sandstone with high homogeneity, a few natural joints, and good integrity is selected. The physical and mechanical properties of the rock are shown in Table 1. The “Rock Experiment Specification for Water Conservancy and Hydropower Engineering” (sl264-2001) is taken as the standard of experimental processing [28]. The red sandstone is processed according to the standard into a cylinder with dimensions φ 50 mm × 100 mm. Further, the specimen is flattened on both ends and controlled within ±0.05 mm. The precast cracked red sandstone specimen is shown in Figure 1.

2.2. Experimental Details

In this work, RYL-600 which is a shear rheological tester developed by the Provincial Key Laboratory of Safe Mining Technology of Coal Mine, Hunan University of Science and Technology, is used (Figure 2). The uniaxial compression tests of the precracked red sandstone under different dip angles are carried out. The instrument is used for conducting strength tests, relaxation tests, triaxial rheological tests, and shear rheological tests, with the characteristics of high stiffness, high control precision, and good stability. The technical parameters considered for curing the experiment is as follows. The maximum axial test force and the maximum load are taken as 600 kN and the error of deformation measurement is considered to be ≤± 0.5%. The axial control modes that are used for conducting the experiment are displacement control mode and experimental force control mode. Further, the control accuracy is considered as ±0.1% and ±0.2% for displacement error and experimental force measurement error, respectively.

2.3. Experimental Plan

Initially, the red sandstone is prefabricated into an unclosed single-crack rock sample with crack dip angles α of 0°, 15°, 30°, 45°, 60°, 75°, and 90°, respectively. Further, the crack is located at the center of the specimen. The length and width of the prefabricated crack are 15 mm and 1 mm, respectively. According to the crack angle, the rock is divided into seven groups, three rock samples in each group which total up to twenty-one rock samples. The axial loading system of the RYL-600 shear rheometer is used in the experiment, and the loading speed is controlled at 100 N/s. The plan of the uniaxial compression test for the precracked red sandstone is shown in Table 2. According to the loading rate (0.5 MPa/s) which is recommended by the International Society of Rock Mechanics, and considering the need for an experimental observation in the cracking process, the set loading rate must be lower than 0.5 MPa/s. Therefore, the loading rate is set as 100 N/s in this experiment.

3. Analysis of the Experimental Results

The analysis of the experimental results such as the stress–strain curves, deformation parameters, and the failure modes are discussed in detail below.

3.1. Stress-Strain Curve

The stress-strain curve reflects the progressive process of the internal pore compaction, microcrack closure, new crack initiation, propagation, coalescence, and the instability of the failure. The stress-strain curves of intact red sandstone and single precracked red sandstone (α = 45°) under uniaxial compression loading are shown in Figure 3. According to the loading characteristics of the experimental curve, the stress-strain (σ-ε) curve of red sandstone is divided into five stages for analysis [2931]. The first stage is the microfracture and pore compaction stage. During the initial stage of axial stress loading, the original open structural planes or the microcracks close gradually. The pore volume decreases with increasing axial force, and the rock sample is compacted in order to form nonlinear deformation at the early stage. The σ-ε curve is approximately a “concave” nonlinear segment, and the axial force at the end of this stage is the closure stress σcc of the rock sample [32, 33]. The second stage is the elastic deformation stage. During this stage, the σ-ε curve is approximately inclined “straight line.” The pore volume strain tends to be stable and the axial force at the end of this stage is the crack initiation stress σci. The third stage is crack initiation and stable growth stage. In this stage, with the increase in axial force, the microcracks further crack to form new cracks, and the cracks continue to expand at a stable speed. Further, the loading time is longer, and the axial force at the end is the damage strength of rock sample σcd. [32, 33] The fourth stage is accelerated crack propagation and unstable fracture of the rock sample. During this stage, the turning point of the rock sample changes from elastic to plastic, and the yield point is generated. Further, the cracks accelerate to form a certain cross fracture surface, and the rock sample surface is accompanied by an unstable fracture. When the bearing capacity of the rock sample reaches the peak strength, the internal structure of the rock sample is dislocated. Further, it is accompanied by a weak cracking sound with a short duration. The phenomenon of sliding off the macroshear plane does not occur, and the overall shape of the rock sample is maintained. The axial force at the end of this stage is the peak strength σc. The final stage is the postpeak fracture instability stage. During this stage, the bearing capacity of the rock sample decreases slightly with the increase in deformation. Further, the intersecting fracture surfaces become macroscopic fracture surfaces. Thereafter, the bearing capacity of the rock sample is suddenly lost with a sudden loud noise. The axial stress drops to zero instantaneously, and finally, the rock sample becomes unstable. This phase takes place in a very short duration of time.

In Figure 3, the results of the precracked sample are compared with the results of the intact sandstone sample. It is observed that the prefabricated crack causes initial damage to the rock sample, resulting in a lower σc value and shorter duration of each stage. The obvious stress drop before the peak value may be due to the initiation and expansion of microcracks in the rock sample. The macrofracture on the surface of the rock sample and the slip and dislocation of the rock structure may lead to the decrease of the rock bearing capacity and the drop of the axial stress of the rock sample.

3.2. The Strength and the Deformation Parameters

Figure 4 shows the σ-ε curves of precracked red sandstone samples and intact rock samples under uniaxial compression. Figure 5 shows the comparison of the peak strength changes and marks the instability fracture point in the figure. The color from blue to red indicates the value of σc. From Figures 4 and 5, it is observed that the peak strength σc increases, and the strain ε increases with the increase in crack dip angle α. Further, from Figure 5, it is observed that when the rock sample reaches σc, the color changes from red to blue instantly. This indicates that the axial stress decreases rapidly, and finally the rock sample develops instability fracture. The instability fracture occurs after the decrease of σc. Therefore, σc can be taken as the critical stress value for instability failure. The peak strength and elastic modulus increase with the increase in the crack dip angle.

By combining the obtained experimental data, it is found that the relationship between the peak strength and the crack dip angle is fitted linearly (Figure 6(a)). Moreover, it is observed that the peak strength and the elastic modulus increase as the crack dip angle increases. Among them, the elastic modulus increases nonlinearly as the crack dip angle increases (Figure 6(b)).

Figure 7 shows the relationship between the maximum strain, peak strength, elastic modulus, and the crack dip angle. The results show that the maximum strain (i.e., the ultimate strain of the rock specimen after failure and the instability) increases step by step as the crack dip angle increases, and thus it further increases the amplitude in the later stages. Further, the variation of the peak strength and the elastic modulus with the crack dip angle is found to be similar. Moreover, the peak strength and the elastic modulus increase slowly in the early stage but increase sharply after the crack dip angle exceeds 45°. Further, it is observed that they continue to increase almost linearly until the maximum value at the intact rock sample. Furthermore, from Figures 57, it is observed that the peak strength of the red sandstone specimens with various crack dip angles decreases to varying degrees when compared with the intact specimens. Therefore, in order to characterize the damage effect of the prefabricated crack with different dip angles on the rock specimens, the damage coefficient Kd is defined as shown in the following equation [34]:

In the previous equation, σcw refers to the peak strength of the intact rock specimen and σc refers to the peak strength of the rock specimen with a single preexisting crack.

From Figure 8, it is observed that the average damage coefficient of the rock samples is 0.533 when the crack dip angle is 0°, and the average damage coefficient of the rock samples is found to be 0.056 when the crack dip angle is 90°. Further, the damage coefficient of the red sandstone specimen with a single prefabricated crack generally shows a nonlinear decreasing trend with the increase of the crack dip angle, and the reduction gradually becomes flat. Thus, the damage coefficient Kd can reflect the damage effect of the crack dip angle on the peak strength of the rock sample in a better way.

3.3. Failure Mode of Red Sandstone

The uniaxial compression tests of precracked sandstone and intact sandstone are carried out by using a shear rheometer. The final failure mode is shown in Figure 9. The results show that the prefabricated cracks of the rock sample are through fractures, and the failure modes of both sides of sandstone are similar. Therefore, only one side is selected to analyze the influence of the crack dip angle α on the failure mode. Tensile crack (T) and shear crack (S) is produced in rock sample under loading, and the angle between shear crack and prefabricated crack is β.

When the dip angle of the prefabricated crack is α = 0°, the inclined shear cracks are formed at both the ends of the prefabricated cracks which are accompanied by macrotensile cracks. The airfoil crack forms between the shear crack and the prefabricated crack and finally penetrates through the rock sample to form a stepped failure. When the dip angle of the prefabricated crack is α = 15°, a shear crack inclined to the upper left is formed at the left end of the prefabricated crack. A shear crack inclined to the right and down is formed at the right end of the prefabricated crack which is accompanied by a macrotensile crack. With the action of load, some fragments are formed, and finally, the rock sample is destroyed during the antistep failure. When the dip angle of the prefabricated crack is α = 30°, the fracture mode of the rock sample is similar to that of α = 0°. The angle β between the shear crack and the prefabricated crack at both the ends of the prefabricated crack increases and the increased angle is approximately equal to α.

When the dip angle of the prefabricated crack is α = 45°, under loading condition, the rock sample forms a small angle between the shear crack and the airfoil crack at both the ends and thereafter exhibits the fragment peeling phenomenon. The final shear mode crack penetrates through the specimen and causes failure. When the dip angle of the prefabricated crack is α = 60°, with the increase in the dip angle, the peak strength of the rock sample increases gradually, and its bearing capacity is stronger. Under the combined action of composite stress, the ultimate failure mode of the rock sample is two cones with one block on each side.

When the dip angle of the prefabricated crack is α = 75°, the inclined shear cracks are formed at both the ends of the prefabricated crack, which are approximately coplanar with the prefabricated crack. Finally, the rock sample is destroyed in the form of oblique shear, and the volume expansion of the rock sample is large. When the dip angle of the prefabricated crack is α = 90°, the rock samples show uneven failure. Small fragment peeling and block shedding exist simultaneously. The peak strength is close to the intact rock sample, and the accumulated energy is stronger than the other precracked sandstone. Therefore, the energy release will lead to more serious damage and uneven failure. The damage of intact sandstone is more serious than that of the precracked sandstone. Debris, small fragments, and blocks are more serious than the precracked sandstone. Although the peak strength is slightly higher, the accumulated energy is much higher than that of the precracked sandstone, and the uneven failure along the shear plane is more serious.

With the increase in the dip angle α of the prefabricated crack, the crack and the prefabricated crack is subjected to combined stress. The secondary cracks around the prefabricated cracks initiate and gradually expand to both the ends of the specimen, which eventually leads to failure. The prefabricated crack will change the failure mode of the sandstone, crack initiation, and propagation to ultimate instability under uniaxial compression. When 0° ≤ α ≤ 60°, the crack propagation path is regular, the dip angle α of the prefabricated crack increases approximately with the angle β between the shear crack and the prefabricated crack, and αβ. When α > 60°, heterogeneous failure of rock sample occurs. The failure mode of the precracked sandstone under uniaxial compression is as follows: when 0° ≤ α ≤ 30°, the shear failure of the precracked sandstone transits to tension shear composite failure; when 45° ≤ α ≤ 60°, there is a transition from tensile shear composite failure to tensile failure of precracked sandstone; and when 75° ≤ α ≤ 90°, the tensile failure of precracked sandstone is transitional to heterogeneous failure.

4. Numerical Experiment Using PFC

The details of the numerical experiment conducted using PFC is discussed below.

4.1. Basic Theory of PFC

When particle flow code PFC is used for simulation, the particle shape and contact state should be correspondingly assumed. The basic assumptions of PFC include the following. The basic particle element of the numerical model is assumed to be a disk or a sphere, and the particle element is considered as a rigid body. In a very small space area, the particle elements are assumed to be almost in the form of point contact, and the contact parts are assumed to have certain bonding characteristics. Further, it is assumed that there is a certain amount of overlap at the contact between the particle elements, and the relationship between the particle element and the contact force is established by the force-displacement theory [3537].

In the numerical calculation process of PFC, Newton’s second law and force-displacement law are the basic theories, in which the force-displacement law constitutes the physical equation and Newton’s second law constitutes the motion equation [35]. Ten contact models are provided in PFC5.0. In this work, the parallel bonding model in PFC is mainly used, and the schematic diagram of the parallel bonding model is shown in Figure 10.

4.2. Construction of the Numerical Model

The discrete element model of red sandstone with a single precrack is established by the following method. The rock particles in the model are simulated by using a single sphere. PFC simulation of uniaxial compression includes the following four steps. In the first step, the numerical experiment container is generated, and in the second step, the parallel bonding parameters are set. During the third step, the model of a single crack is established by using AutoCAD, and the contour of a single crack is generated by importing into the PFC program. Further, by deleting the particles in the crack contour, an unclosed crack is formed. Finally, in the last step, the upper and lower walls are controlled in order to perform uniaxial loading with a loading speed of 0.1 mm/s.

According to the abovementioned establishment of a discrete element uniaxial compression experiment, the parameters of the red sandstone uniaxial test are calibrated. The parallel bonding model is used to simulate and the mesocontact model parameters of the red sandstone are obtained (Table 3).

Figure 11 shows the comparison results of the uniaxial tests and numerical experiments for red sandstone. From Figure 11(a), it is observed that the stress–strain curves obtained by numerical experiments are in good agreement with the laboratory experiments. From Figure 11(b), it is observed that the shear plane formed by the failure of the numerical test specimen is consistent with the uniaxial test results.

4.3. Failure Evolution of Red Sandstone with a Single Prefabricated Crack

Using the calibrated mesoparameters, the uniaxial compression experiments of the single precracked sandstone under different dip angles are simulated (Figure 12). Figure 12 shows the experimental failure diagram, PFC simulation failure diagram (red for tensile crack and black for shear crack), contact force chain diagram, and the displacement field diagram. From Figure 12, it is observed that, during the PFC simulation, microcracks initially appear in the sample. Further, with the increase in load, the microcracks converge together, and cracks and shear planes appear in the macroscopic view. However, the red tensile cracks are more in the final failure of rock samples, which indirectly indicates that tensile microcracks contribute more towards the failure of the rock samples than the shear cracks. With the increase in crack dip angle α, the nonuniformity of the specimen failure is enhanced, and the nonuniformity of the PFC simulation specimen is also enhanced. However, there are some differences between the simulated and the actual failure samples. The results show that the contact force chain is very weak at the precrack and macroshear plane, which indicates that the original contact force is weakened due to particle fracture. Therefore, with the increase of load, the bearing capacity of the specimen will decrease. The displacement at the shear plane of the PFC simulation sample is larger and is consistent with the experiment. Thereafter, the bulging effect of the sample is simulated. With the increase in the load, the fragments fall off and the fragments eject from the rock sample. This is due to the weakened contact force between the particles. Thereafter, fracture occurs with the formation of microcracks and the cracks converge. Thus, it provides a prerequisite for the transformation of elastic strain energy into kinetic energy, and, finally, the rock sample shows debris shedding and ejection.

4.4. Analysis of the Damage and the Energy Change for the Rock Sample

In order to reveal the damage and the failure mechanism of the rock with a single precrack and to understand the process of crack initiation and propagation, based on the fracture mechanism of particle bonding in PFC, the mesodamage variable D is defined as D = Cd/C, where C is the total number of bonds between the particles and Cd is the number of broken bonds.

The energy can be traced and recorded in PFC. This includes the boundary energy (Eboundary), strain energy of the parallel bond (Epbstrain), particle sliding energy (Eslip), and the particle strain energy (Estrrain). The accumulated elastic strain energy of the rock sample is composed of the Epbstrain and the Estrain. In this paper, the elastic strain energy is defined as We = Wpb + Ws, where Wpb is the strain energy of parallel bond and Ws is the particle strain energy. The damage statistics of the simulated samples are shown in Table 4.

From Table 4, it is observed that there is a positive correlation between the number of microcracks and the mesodamage variable D. Therefore, the damage degree of the rock sample is described by D. It is also observed that, as the crack dip angle increases, the peak strength, the number of microcracks, and the mesodamage variable increase. Thus, the damage degree of the rock increases accordingly.

The relationship between the stress, energy, and axial strain during loading is shown in Figure 13. From Figure 13, it is observed that the elastic strain energy reaches the maximum value on loading when the rock sample reaches the peak strength. Further, as the bearing capacity of the rock sample decreases, it is observed that the accumulated elastic strain energy is transformed into other forms of energy. At the same time, a certain number of microcracks occur in the region of the compressive stress concentration, and macrocracks are formed.

From the energy point of view, it is found that the elastic strain energy We is transformed into the sliding energy Wslip. Further, a part of the energy of Wslip is used to drive the rock sample in order to produce microcracks, and the macrofragments fall off in the crack cross regions. Another part of the energy of Wslip makes the rock sample unstable, which may drive the rock sample to eject. Finally, the microcracks coalesce in order to make the specimen fail.

5. Conclusions

In this paper, studies are carried using the single precracked red sandstone. By conducting the uniaxial compression test and PFC2D simulation, the mechanical properties, crack propagation, failure mode, and the failure evolution characteristics of the single precracked red sandstone under different crack dip angles are studied. From the study, the following conclusions are obtained. The stress–strain curve of the precracked sandstone can be divided into five stages, namely, microfracture and pore compaction stage, elastic deformation stage, crack initiation and stable growth stage, crack acceleration and unstable fracture stage, and postpeak fracture instability stage.

Further, the peak strength, peak strain, and elastic modulus parameters of the precracked rock samples increase with the increase in crack dip angle, and the ultimate failure mode of the precracked rock samples transits from the “ladder” type failure to slope uneven failure.

PFC2D is used to study the crack propagation, force chain change, and displacement field of the single precracked red sandstone. The results show that tensile microcracks contribute more to the failure of rock samples than the shear cracks. Further, the contact force chain is very weak at the prefabricated cracks and the macroshear plane, which indicates that the original contact force is weakened due to particle fracture, so the bearing capacity of the specimen decreases with the increase in load.

In addition, from the results of the PFC simulation, it was observed that, as the uniaxial load increases, the contact force between the particles of the specimen weakens and the microcracks gradually propagate, resulting in the phenomenon of fragments falling off and the rock sample fragment ejection. Further, from the viewpoint of dynamics, the gradual convergence of the microcracks provides a prerequisite for the transformation of the elastic strain energy into kinetic energy.

Data Availability

The data used to support the findings of this study are available within this article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Natural Science Foundation of Hunan Province, China (no. 2018JJ3166), and the Natural Science Foundation of China (nos. 51504091, 51774131, and 51774130).