Abstract

The complex mechanical properties of deep surrounding rocks during excavation and unloading have always been a concern in engineering communities. Based on deep roadway engineering, the mechanical properties of granite strain-softening were investigated by laboratory tests. An exponential relationship between granite peak softening modulus and confining pressure was obtained using a nonlinear fitting method. A strain-softening model was developed based on plasticity theory for granite which took into account the nonlinear dilatancy angle and confining pressure. A mathematical model was developed using FLAC3D as platform. By developing a numerical model for a deep roadway, the local fracture characteristics of deep surrounding rocks under strain-softening conditions were evaluated. It has been found that the postpeak failure of granite had a tendency for brittle-ductile transformation. Under high confining pressure conditions, granite exhibited brittle failure characteristics during the postpeak period, and the postpeak softening modulus was decreased with the increase of confining pressure. From FLAC3D numerical calculations, it was found that the numerical models of different mesh densities had basically the same characteristic curves of surrounding rocks, which showed that the local cracking phenomenon had little effect on calculation results when the convergence constraint method was applied for the calculation of the stability of supporting structures and surrounding rocks. It was seen from the numerical simulation analyses of surrounding rock local fracture properties in deep roadways that plastic shear strain appeared in the local areas of roadway vault and arch foot, which basically coincided with the damage location and depth of surrounding roadway rocks.

1. Introduction

Local rupture of rock masses, which refers to successive ruptures and unfractured areas in surrounding rocks on roadway sides and working faces, is common in deep mining [1, 2]. The local rupture phenomena of rock masses decrease the bearing capacity of surrounding rocks. This process is considered as a major reason for the progressive failure of surrounding rocks after the peak excavation process [3]. In actual engineering, brittle failure characteristics are generally observed after the strength of rock reaches its peak value. That is to say, rocks generate and develop cracks under stress and show complex mechanical behaviors such as strength attenuation and expansion in the areas with the maximum strength [4, 5]. Laboratory tests revealed a clear correlation between the mechanical properties of rocks after the peak and confining pressure [6]. In other words, rocks are transformed from brittle failure under low confining pressures to plastic flow failure under high confining pressures. Postpeak brittleness and ductile failure processes could be attributed to the rock strain-softening behavior. In deep rock mining, the strain-softening behavior of rocks can highly affect the local fracture characteristics and stability of surrounding rocks.

Many research works had been conducted on the local fracture of deep rock masses. Ma et al. [7] showed that the creation of zonal disintegrations is mainly dependent on in-suite stress levels and material heterogeneities. The fracture patterns formed during zonal disintegration depend on the size and shape of the tunnel as well as local weakness distribution in surrounding rocks. Makarov et al. [8] established a novel non-Euclidean mathematical model rock masses under high stress values and developed methods to determine model parameters and major relationships between rock mass strength and cracking zone width. Bi and Zhou [9] developed a GPD method-based constitutive model for the simulation of zonal disintegration in deep rock masses. Varas et al. [10] proposed a theoretical analytical method for strain local rupture and investigated the distribution of local fracture zones by establishing numerical models with different fine levels. Egger [11] comprehensively investigated the strain localization of tunnels and found that the phenomenon of strain localization was common during tunnel excavation and unloading. This phenomenon was found to be related to the lithology and ground stress of surrounding rocks and had a certain distribution regularity. The developed theoretical analytical method was used to explain and analyze strain localization phenomenon, and valuable theoretical results were obtained. Li et al. [12] applied numerical methods to simulate the occurrence and evolution of regional rock failure in deep rocks. In relevant findings of recent years, certain randomness and empiricism are witnessed in the peak and postpeak strength parameters of rock masses and the degradation law of rock mass postpeak strength when calculating local fractures in rock masses, especially when numerical analysis software is applied to simulate the calculations.

Investigation of the rock strain-softening behavior can be performed by three approaches. The first approach is to start with laboratory tests and the analysis of postpeak failure characteristics of rock samples. For example, Perez-Rey. et al. [13] studied the total stress-strain responses of several types of intact rock samples such as sedimentary, metamorphic, and igneous rocks through more than 300 triaxial and uniaxial strength tests. Lu et al. [14] performed routine triaxial compression tests on soft mudstone samples to achieve complete stress-strain relationship curves under different confining pressures and obtained a yield surface model for soft rocks having the two state parameters of generalized cohesion and generalized internal friction angle. Yu et al. [15] evaluated the brittle-extended-plastic transition characteristics of Jinping deep-buried marble peaks after experimental deformation. The second approach was the theoretical analysis. For example, Yu et al. [16] studied a time-dependent deformation mechanism in depth and obtained the viscoelastoplastic deformation using a nonassociated flow rule and the modified Nishihara model. Zareifard [17] proposed a numerical-analytical solution for the ground response curve of a tunnel excavated in an elastoplastic strain-softening rock mass by taking into account the gravitational loading. Alejano et al. [18, 19] developed a postpeak strain-softening model on the basis of the Mohr–Coulomb criterion and Hoek–Brown criterion. Sharan et al. [20] provided an analytical expression for the strain-softening model in the Hoek–Brown elastic-brittle-plastic model during the stress-deformation of surrounding rock masses. In the third approach, numerical calculation methods are proposed for reasonably expressing the strain-softening behavior of rocks based on experimental research and theoretical analyses. For example, Ghorbani and Hasanzadehshooiili [21] studied all influential parameters on strain-softening, including novel and previously studied features, and proposed a complete solution for calculating the GRC of circular tunnels. Wang et al. [22] developed a numerical calculation method for the simulation of the evolution characteristics of parameters during the strain-softening of geotechnical media. Also, the stress-strain curves of geotechnical materials during the evolution of different strength parameters were calculated and analyzed based on the Mohr–Coulomb criterion. Zhang et al. [23] presented a large strain numerical solution of a circular tunnel for the ground reaction curve (GRC) in strain-softening rock masses. Sun et al. [24] investigated the interaction between the surrounding rocks of the chamber and the supporting structures by numerical simulation using a plasticity theory-based rock strain-softening model. The above research results are very valuable in engineering applications. However, developing a postpeak strain-softening model for rocks is very complicated, which limits its wide application in practical engineering. The theoretical calculation was oversimplified, and the calculation results had large errors.

In view of the above problems, in this paper, we have studied the mechanical characteristics of strain-softening of granite samples in deep roadway engineering using laboratory tests. An exponential relationship was found between the peak softening modulus of granite and confining pressure using a nonlinear fitting method. Assuming a constant rock dilatancy angle, a plasticity theory based on the strain-softening model of granite was developed by taking into account the influences of dilatancy angle and confining pressure. A mathematical model was developed using FLAC3D as platform. By developing a numerical model for a deep roadway, the local fracture characteristics of deep surrounding rocks under strain-softening conditions were evaluated. The research results enriched the research system of postpeak mechanical behavior and provided references for practical engineering.

2. Rock Strain Localization and Mechanical Model

2.1. Stress Bifurcation and Localization

In a mathematical point of view, stress bifurcation behavior is due to different combinations of boundary conditions and equilibrium equations. In other words, when the critical value of a parameter is reached, stress bifurcation or strain localization occurs [25]. In compression deformation of rocks, local shear fracture zones are created leading to heterogeneous deformation failure. When surrounding rocks on roadway surface reach a specific unloading level pi during the excavation and unloading of roadways, immediate stress bifurcation and strain localization take place, as shown in Figure 1.

2.2. Mechanical Model of Rock Strain-Softening

In the plastic strain-softening model for rock materials, plastic potential energy equation and yield criterion are expressed by stress tensor and include softening parameters, as shown in the following equation [26]:where is the maximum principal stress and is the minimum principal stress.

The Mohr–Coulomb yield criterion could be stated aswhere is the uniaxial compressive strength of rock and is defined as

When it is assumed that c and in the Mohr–Coulomb constant decay linearly with , the strain-softening relationship of Mohr–Coulomb could be stated aswhere is the peak parameter, is the residual parameter, is the critical softening parameter of rock from strain-softening to residual stage, and can replace the cohesion c and friction angle in the Mohr–Coulomb model and the constants m and s in the Hoek–Brown criterion. In other words, when  = 0, rock mass is in elastic deformation stage, when , it is in the strain-softening stage, and when , it is in the residual stage. The rock softening process is calculated from the postpeak stress-strain curve slope, i.e., the softening modulus M of rock, as shown in Figure 2. When M approaches infinity, elastic and brittle deformation of the material takes place, while when M approaches zero, an ideal elastoplastic deformation is observed. Both cases could be considered as special cases of rock strain-softening behavior.

Plasticity parameter can be stated in the form of intrinsic variables. In other words, plasticity parameter is defined as plastic shear strain, which is calculated from the difference between the minimum and maximum principal plastic strains. That iswhere is the plastic shear strain and and are the maximum and minimum plastic strains, respectively.

Based on the simplified strain-softening curve after the rock peak shown in Figure 2, the maximum principal plastic strain could be stated aswhere is the maximum principal strain of elasticity before the peak, is the softening strain after the peak, and is the maximum principal strain of elasticity. Under constant confining pressure, the parameters in equation (7) could be stated aswhere and are peak and residual principal stress, respectively.

Considering dilatancy angle and based on the form of strain increment, dilatancy angle could be expressed as [27]where (j = 1, 3) is plastic principal strain rate. equation (8) can be transformed intowhere

Under static conditions, the relationship between and at constant dilatation angle can be calculated from equation (10) as

Considering the variation of nonlinear dilatancy angle, the dilatancy angle can be calculated using the plastic strain trajectory, as shown in Figure 3. At any point a, b, and c in the plastic strain trajectory, the plastic strain of two discontinuities (ab and bc intervals) can be expressed aswhere

Finally, the mathematical relation of constant dilatancy angle and confining pressure could be obtained by combining equations (1) and (6)–(11):

In FLAC3D, plastic shear strain of the Mohr–Coulomb strain-softening model is expressed by plastic parameter , and its incremental form iswherewhere (j = 1, 3) is the main increment of plastic shear strain.

When considering rock dilatancy behavior, the nonassociated flow law corresponding to plastic shear potential function in FLAC3D could be expressed as

The relationship between plastic parameters and plastic shear strain in FLAC3D is

3. Experimental Study on Strain-Softening Characteristics of Deep Granite

3.1. Granite Laboratory Test

The granite rock samples selected for experimental tests were collected from the roadway excavation face at a depth of 1,000 meters in Xincheng Gold Mine of Shandong Province in China. Rock blocks had a uniform texture with good overall integrity. Rock blocks were cut into φ 50 mm × 100 mm cylindrical and φ 50 mm × 25 mm disc samples. The processing accuracy of the samples met the requirements of the test specifications recommended by the International Rock Mechanics Society. Brazilian splitting and triaxial compression tests of rocks were conducted on a TAW-2000 microcomputer controlled electrohydraulic servo rigid pressure tester (Chaoyang Instrument Co., Ltd). The test system had the advantages of high accuracy, fast data processing, and high operability and reliability.

In triaxial compression tests, the confining pressure was set at 0 MPa, 2 MPa, 5 MPa, 10 MPa, 15 MPa, and 30 MPa. First, a small hydrostatic pressure was applied to the sample, and then a constant confining pressure was applied. Then, an axial load was applied at a constant strain rate. Axial strain , radial strain , axial stress , and radial stress were obtained from the experimental data obtained by the automatic data acquisition system. It was found from repeated experiments on a large number of rock samples that, because of the high homogeneity of granite samples, the dispersions of elastic modulus and peak strength of the rocks were low. The complete stress-strain curves were experimentally obtained for granite samples from typical triaxial compression, as shown in Figure 4(a).

In the Brazilian splitting test, a splitting fixture was placed in the middle of the upper and lower pressure plates of the testing machine. Once the test was started, the clamping screws on both sides of the fixture were loosened and a uniform load at a loading speed of 0.5 MPa/s was applied until failure. Tensile strength was obtained using the automatic data acquisition system.

It was seen from Figure 4(a) that the rock samples showed elastic and brittle failure during uniaxial compression and were dropped rapidly after reaching peak strength . The residual strength of rock samples was about 0 MPa, and they were broken, as shown in Figure 5(a). The failure of rock samples under low confining pressure () still showed elastic and brittle failure characteristics. By further increase of confining pressure, the stress-strain curve slope after rock peak was gradually decreased. As the confining pressure was increased, the failure characteristics of rocks tended to change from brittleness to ductility, and the rock failure mode was mainly shear failure after the peak. Based on the peak strength curves of rocks under different confining pressures, it was found that confining pressure also affected the elastic parameters of rocks. In this paper, we have focused on the postpeak mechanical characteristics and connected the peak strength point to the residual strength point. The slope of this line can be defined as the softening modulus of the rock sample, which is denoted by M in this paper, as shown in Figure 4(b). Softening modulus M could be obtained from elastic modulus E as

The peak elastic parameters of rock were simply calculated, and the values of elastic modulus E = 34.4 GPa and Poisson’s ratio  = 0.28 were obtained.

Characteristics of typical rock destruction are shown in Figure 5.

3.2. Study on the Peak and Postpeak Strength Characteristics of Granite

In the analysis of engineering stability, linear Mohr–Coulomb and nonlinear Hoek–Brown criteria are generally considered as the yield criteria of rock masses. The basic parameters of the two models were still obtained from laboratory tests. When the Hoek–Brown criterion was adopted, the uniaxial compressive strength and rock hardness m of the complete rock samples were obtained from laboratory tests and the mechanical parameters of surrounding rocks were obtained by combining with the GSI surrounding rock rating system. When the Mohr–Coulomb criterion was adopted, the uniaxial compressive strength of intact rock and the peak and residual friction angles and cohesions c of the rock were determined. The fitting method proposed by Hoek and Brown [26] was applied for the calculation of peak and residual strength principal stress curves of the two criteria in this paper, as shown in Figure 6.

It can be seen from Figure 6 that high fitting accuracies were obtained for both Hoek–Brown strength curve and Mohr–11Coulomb intensity curve. The residual strength and peak value were obtained by the fitting method, as shown in Table 1. The Mohr–Coulomb strength parameters are easy to obtain, and the influence of rock dilatancy is considered in the practical engineering applications. Therefore, in this paper, the postpeak mechanical model of rock was investigated according to the Mohr–Coulomb strength criterion.

It can be seen from the analysis above that softening modulus coefficient is the key to build a strain-softening model for rock. Optimal exponential fitting was performed in this work using the indoor triaxial compression test data. The equation of under different confining pressures could be obtained as

Compared with other fitting equations, the exponential fitting method is more suitable for the secondary development and application of FLAC3D. The average correlation coefficient of the fitting curve is more than 0.9, indicating that the nonlinear fitting effect of the model is remarkable. The fitting curve is shown in Figure 7.

4. Numerical Simulation of Local Fracture in Deep Granite

4.1. Numerical Model and Parameters

The local fracture characteristics of deep granite under strain-softening condition were evaluated using Mohr–Coulomb and Hoek–Brown strain-softening models. The model size was adopted to be 10 R, where R = 2 m to avoid boundary effects. The boundaries were normal fixed constraints. The initial ground stress was 25 MPa, and the lateral pressure coefficient was 1. Boundary conditions and calculation model are shown in Figure 8. In the process of numerical calculation, the strain-softening model of granite was written into the FLAC3D strain-softening model by FISH language. The cohesion and internal friction angle were set as the bilinear function of the plastic parameter , and the same plastic strain was experienced from the softening to residual stage. The nonlinear dilatancy angle of rock varied with the increase of plastic parameter . After excavating the deep roadway surrounding rock, the virtual supporting pressure Pi was gradually reduced to realize the stress release process of the surrounding rock and analyze the evolution characteristics of the local rock local fracture. The calculation parameters of surrounding rock are shown in Table 1.

4.2. Local Fracture Characteristics of Granite

Calculations of the nonlinear dilatancy angle variation characteristics of the surrounding rock revealed that, with the continuous release of the surrounding rock stress, the dilatancy angle distribution moved from the surrounding rock surface to the deep layer, as shown in Figure 9. When the stress of the surrounding rock was released by 20%, the surface of the surrounding rock was the main area of capacity expansion. With the continued release of the stress of the surrounding rock, the dilatancy angle decreased and gradually shifted to the depth of the surrounding rock. It can be seen that the main capacity expansion phenomenon occurred in the surrounding rock near the working face with the excavation of the tunnel. When the working face continued to advance, the influence of the dilatancy angle on the surrounding rock surface gradually decreased and the dilatancy angle gradually decreased and moved to the deep rock mass. The nonlinear dilatancy angle was not concentrated in a certain area but changed with the increase of shear strain.

Calculations of the local fracture characteristics of surrounding rock revealed that the localized distribution ranges of the two yield criteria were relatively close, as shown in Figures 10(a) and 10(b). When the distortion of surrounding rock reached a threshold value, the inner wall section of the surrounding rock was first broken down and the rock entered residual strength stage. By further release of stress from surrounding rock, stress continued to expand into the surrounding rock. In addition, unbroken rock masses were distributed among banded spirals. Due to higher strain-softening parameters of granite, the ranges of stress bifurcation and strain localization generated by surrounding rocks were larger. The spiral fracture zone inside surrounding rock was very obvious. The main reason was that the strength of surrounding rock was decreased and its initial stress was high. During stress release, damage to surrounding rock was more serious such that the damage degree of surrounding rock surface in the cloud diagram of the Mohr–Coulomb model was more serious than that calculated by the Hoek–Brown model.

The bolt elements were installed when the load was released by 60% after surrounding rock excavation (length = 2.0 m, diameter = 18 mm, spacing = 1.2 m). By analyzing the local fracture characteristics of the surrounding rock under the condition of bolt support, it was known that the bolts had a controlling effect on the local fracture range of the surrounding rock. The calculation results are shown in Figures 10(c) and 10(d). The plastic strain rate of the bolts reinforcement zone decreased obviously, but the surrounding rock still showed the characteristics of local fracture. The calculated axial force distribution characteristics of bolts are shown in Figure 11. There was a certain difference in the axial force of bolts calculated by the two constitutive models. The main reason was that the Mohr–Coulomb model took into account the change of nonlinear dilatancy angle, which led to the increase of axial force of bolt. This was more in line with the actual situation of the project.

4.3. Analysis of Surrounding Rock Deformation Characteristic Curve

Many designers use numerical simulation methods for calculating the support structure of deep-buried roadways in engineering practices and usually establish rough numerical models. Numerical models with different grid densities were applied to investigate the effect of local cracking on the calculation of the stability of surrounding rock and supporting structure in this paper. The models are shown in Figure 12, where R = 2 m. The Mohr–Coulomb model and Hoek–Brown model were used in the calculation. It was assumed that the virtual internal support stress Pi of surrounding rocks was gradually decreased by 2% of the original rock stress each time, for 50 times in a cycle. Thus, the relation curve between the displacement and pressure of surrounding rock was obtained.

It can be seen from Figure 13 that the Hoek–Brown model and Mohr–Coulomb model were used to calculate the three different mesh accuracy models. The results showed that surrounding rock displacement in superfine model (a) was larger than that of the other two models, but the difference was small. It was seen that the obtained calculation accuracy met the requirements when calculating the stability of surrounding rocks and supporting the structures by numerical analysis methods, even when local cracking was not considered. However, the superfine mesh model could better study the distribution characteristics of local rupture and the interaction between supporting structure and surrounding rock, which could provide guidance for deep mining projects.

5. Study on the Local Fracture of Surrounding Rock in Deep Roadway

The design production level of Xincheng Gold Mine of Shandong Province in China is 1505 m (ground elevation + 50 m). Surrounding rocks in this section are mainly composed of granite with almost uniform and complete distribution of rock mass. The maximum horizontal ground stress value is 30 MPa, and the vertical ground stress value is 25 MPa. This is a high geostress environment. For accurate prediction of potential damage area of roadway surrounding rock, a two-dimensional fine numerical grid element model of roadway was developed in this paper. Based on the strain-softening model of granite, the location of post-peak failure area of deep surrounding rock was analyzed. The size and boundary conditions of the numerical model are shown in Figure 14, where horizontal ground stress value Px = 30 MPa and vertical ground stress value Py = 25 MPa. The surrounding rock parameters in calculation are shown in the Mohr–Coulomb strength parameters in Table 1.

According to calculations and analyses, under the action of high geostress, a large area of plastic shear failure area still appeared in roadway vault and arch foot although the lithology of roadway surrounding rocks was relatively hard, as shown in Figure 15. The damage area of roadway vault was still a low confining pressure area. Surrounding rock failure mode was mainly brittle failure, and failure area did not extend to deep surrounding rocks. During actual construction and for a period of time after construction, some damage took place in vault surrounding rocks, as shown in Figure 16. Damage depth was basically the same as the calculation result. It can be seen that the granite strain-softening model developed in this work could reasonably reflect the postpeak failure characteristics of rocks, and calculation results could provide relevant suggestions for the construction site.

6. Conclusion

Based on the deep mining project, a strain-softening model for granite was established in this work. The local fracture characteristics of deep granite were studied, and the obtained results were applied to engineering practice. The following conclusions are drawn:(1)By analyzing the properties of triaxial complete stress-strain curves obtained for granite, it was found that the postpeak failure characteristics of granite obviously correlated with the confining pressure. The postpeak failure tended to brittle-ductile transformation. The increase of the confining pressure decreased the postpeak softening moduli of rocks, but it was not obvious under low and medium confining pressures. However, rocks still exhibited brittle failure characteristics after the strength of rock reached its peak value under high confining pressure conditions.(2)Based on the theory of plastic mechanics, the functional relationship between postpeak softening modulus of granite and surrounding rock was obtained by the exponential fitting method. The strain-softening model of granite considering nonlinear dilatancy angle was deduced, and the model was solved by FLAC3D. Calculations revealed that the deep granite had obvious local fracture characteristics, and the surrounding rock characteristic curves obtained by the numerical models with different grid densities were basically the same. The results showed that, when the convergence constraint method or numerical model was used to calculate the stability of surrounding rock and supporting structure, the local cracking phenomenon would have little influence on the calculation results.(3)A wide plastic shear strain area was seen in the roadway vault and arch foot by numerical analyses on the failure characteristics of surrounding rocks in the deep granite roadway. It was basically the same location and depth of the surrounding rock failure of the onsite roadway. The postpeak strain-softening model of granite based on the Mohr–Coulomb model can take into account the change of dilatancy angle, which was more suitable for deep rock mass engineering.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors would like to thank the support from the National Key R&D Projects (Grant no. 2017YFC1503101), National Natural Science Foundation of China (Grant no. 51704144), and Talent Project of Revitalizing Liaoning Funding Project (Grant no. XLYC1807107).