Study on the Creep and Fracture Evolution Mechanism of Rock Mass with Weak Interlayers

To study the influence of weak interlayers on the creep failure characteristics of rock masses, based on the continuous-discontinuous method (CDEM), the uniaxial compression creep experiments of rock masses containing weak layers were numerically simulated; and the weakened rock masses under different conditions were analyzed in detail. We focused on the final failure mode and creep curve of the rock mass with a weak interlayer (θ� 30°, d� 20mm, c� 1) as examples by selecting the crack distribution state of the model during compression at different time steps. We analyzed the propagation and convergence mode of cracks in a rock mass with weak layers. *e research results show that the existence of weak interlayers affects the integrity of the rock mass and the creep failure mode.With the increase in the inclination of the weak interlayer, the failure mode of the rockmass changes from shear failure through the weak layer to slip along the weak layer. For shear failure, the total strain and steady-state creep rate of the rock mass first decrease and then increase, showing a U-shaped distribution; as the thickness of the weak interlayer increases, the rock mass always follows the shear in the weak layer. Creep failure occurs on the fracture surface, and the total strain and steady-state creep rate of the rock mass are positively correlated with the thickness. If the thickness continues to increase, there is no significant difference in the creep characteristics of the rock mass; the volume occupied by the soft rock in the body increases, the overall rigidity of the rock mass decreases, and the plastic deformation increases. *e form of creep failure of the rock mass changes from sliding shear failure along the weak layer to sliding shear failure through the weak interlayer.*e total strain and steady-state creep rate of the rock mass increase with the increase in the number of weak layers; the greater the distance between the weak layers, the smaller the total strain and steady-state creep rate of the rock mass.*e slower the crack growth rate, the less likely the rock mass to undergo creep damage.


Introduction
e weak layered rock mass is a natural layered material composed of single rocks of different thicknesses, different inclination angles, and different mineral components. e weak layered rock mass is one of the common types of rock masses in rock engineering, along with construction, mines, slopes, and tunnels [1]. e rock mass with a weak interlayer has both the characteristics of strength anisotropy and the transverse isotropy of layered rock mass deformation, so its mechanical properties are complex [2]. Numerous engineering practices have shown that the creep characteristics of rock masses with weak layers are one of the important reasons for the deformation and instability of the engineering. e long-term stability of these rock projects has always been the focus of attention [3][4][5].
erefore, researching the creep characteristics of rock masses with weak interlayers has important theoretical guiding significance.
At present, the research on rock masses with weak layers usually adopts methods such as laboratory experiments [6][7][8][9], numerical simulation [10][11][12], and theoretical analysis [13,14]. In terms of theoretical research, Jaeger [15] first proposed a transverse isotropy criterion with a single weakface rock mass as the research object and then regarded the internal cohesive force in the Mohr-Coulomb criterion as a function of the layer dip angle. e internal friction force is taken as a constant, and the shear strength formula of the layered rock mass is given and revised. Saroglou et al. [16] and Singh et al. [17] applied the revised Mohr-Coulomb criterion and Hoek-Brown failure criterion to layered rock masses based on predecessors. Zienkiewicz et al. studied the mechanical properties of composite rock masses, proposed a simple multilayered rock mass model, and considered the viscoplastic strain rate of the rock mass as the sum of the contributions of each joint group and rock material. e applicability of these methods to the mechanical properties of composite rock mass has not yet been verified. In terms of physical experiments, because of the limitations of on-site sampling of rock masses containing weak layers and sample processing, it is difficult to conduct indoor experiments on real raw rock materials. erefore, many scholars use similar materials to prepare rock masses containing weak layers [18,19] and conduct corresponding research. Because the rock mass with weak layers is an anisotropic, discontinuous, heterogeneous, and random natural aggregate, similar simulation results often differ greatly from the real situation. To better conduct indoor experimental research on layered rock masses, scholars have proposed a new production method [20,21] using epoxy resin glue or marble glue as a binder and using different types of original rocks, cutting, grinding, bonding, and assembling to make a layered rock mass sample with controllable mechanical properties and inclination of the weak layer.
is method is simple to operate for similar simulation experiments, and the sample production cycle is short, which can reflect the macromechanical properties of the layered rock mass but cannot reflect the real interface bonding in the layered rock mass.
Many scholars have studied the strength characteristics and creep characteristics of rock masses containing weak layers. Xu [22] and others used similar materials to make layered rock masses with weak interlayers and conducted triaxial rheological tests on them to study the rheological deformation laws and accelerated rheological characteristics of layered rock masses with weak interlayers. Xiong et al. [23] used PFC2D to numerically simulate the triaxial compression experiment of layered rock mass and studied the influence of the inclination angle of the interlayer rock mass on the strength of the rock mass under triaxial compression load and the law of crack development and failure mode. Li [24] conducted compression and triaxial compression experiments on the salt rock mass with mudstone interlayers and determined that the mudstone interlayer has a significant impact on the deformation and failure characteristics of the salt rock mass. Huang et al. [25] simulated the excavation process of a tunnel with a weak interlayer using a combination of physical test model and numerical analysis. e results indicated that the position, inclination, and thickness of the weak interlayer are important factors affecting the stability of the tunnel.
Most of the above studies simplified the rock mass containing weak layers into an equivalent homogeneous and transversely isotropic composite material for analysis. e research direction focuses primarily on the location and dip angle of weak layers under different confining pressures. e influence of rock mass strength, deformation characteristics, and stability and some new failure criteria and stability evaluation methods for rock masses containing weak layers have been established; and some new progress has been made. However, previous studies mostly analyze the mechanical properties and failure modes of rock masses containing weak layers under uniaxial or triaxial compression from a macro perspective and do not consider the macroand mesoscopic evolution process of the creep fracture of rock masses containing weak layers.
is article is based on the continuous-discontinuous element method (CDEM), and in it, we establish and solve the creep fracture geometric model of a rock mass with weak interlayers. e fracture process of rock mass with different weak interlayer conditions was numerically simulated using crack initiation, propagation, and connection to penetration. We analyzed the evolution of crack propagation in a rock mass with weak interlayers in each creep stage. e effects of different dip angles of weak interlayers, thickness of weak interlayers, and the number of weak interlayers and spacing of weak interlayers on the final failure mode and creep fracture of the rock mass. e creep rupture evolution mechanism of a rock mass with weak interlayers was given. It provided a theoretical basis for preventing and controlling the creep disaster of a rock mass with weak interlayers.

The Evolution Law of Unit Failure
2.1. Basic Concepts of CDEM. CDEM can be defined as a dynamic explicit solution algorithm based on breakable elements under the Lagrangian system. Strict control equations are established through the Lagrangian energy system, and the dynamic relaxation method is used to display an iterative solution, which realizes a continuousdiscontinuous unified description. It can simulate the entire process of material from continuous deformation to fracture to movement, combining continuous and discrete calculations. e advantages of calculation include the finite element, finite volume, and spring element methods for continuous calculation and the discrete element method for discrete calculation. e numerical model of CDEM can be divided into a block part and an interface part, as shown in Figure 1. e block part comprises one or more finite elements, which are used to characterize the continuous deformation characteristics of the material; and the interface part is the contact surface of the block boundary. Introducing a breakable, onedimensional spring on the boundary of the block will lead to cracks in the material.

Discrete Element Calculation of CDEM.
In CDEM, the key to switching from finite element to discrete element calculation is to perform discrete element calculation at the boundary of the block. e spring system model of the interface is shown in Figure 2.
In the failure calculation, the criterion for interface failure is based on the maximum tensile stress and the 2 Advances in Materials Science and Engineering Mohr-Coulomb composite criterion, divided into tensile failure and shear failure. When the interface is under tensile failure, as shown in Figure 3, the interface contact spring adopts the maximum tensile stress criterion; and the interface force state is shown in the following equation. e tensile force after failure and the force on the node are both zero.
where σ t is the tensile strength of the contact surface, k n is the normal stiffness of the spring, and u is the normal displacement. When the interface is under shear failure, as shown in Figure 4, the interface contact spring adopts the Mohr-Coulomb criterion; the interface force state is as shown in the following equation: In the formula, k s is the tangential stiffness of the spring, v is the tangential displacement, c is the contact cohesion, φ is the contact internal friction angle, and σ n is the normal stress.

Feasibility Analysis of CDEM.
In recent years, many scholars have used the discrete element method to simulate the propagation process of rock cracks, analyze the rock fracture mechanism, and verify the feasibility of this method for simulating rocks by comparing it with the laboratory experiment data [26][27][28][29]. In this section, CDEM was used to simulate a rock mass with prefabricated fractures under uniaxial compression. e fracture dip angles were 0°, 30°, 45°, and 60°, respectively. We compared the simulation results with the DEM experiment and laboratory experiment results adopted by Jiang [30] to verify the feasibility and applicability of CDEM. Table 1 shows the failure patterns of samples with prefabricated cracks simulated by the CDEM experiment, DEM experiment, and laboratory experiment, respectively. Specific parameter selection and other details have been listed in this article [30]. is section only lists the results to verify the feasibility of CDEM. e CDEM simulated crack growth pattern is similar to the DEM and laboratory experiments, and obvious wing crack growth can be observed. e initial crack develops in the direction of principal stress. e test results prove the applicability of CDEM to simulate rock. In the following section, we used CDEM to simulate the crack propagation mode of a rock mass with weak interlayers under different working conditions.

Numerical Model and Research Plan of Rock Mass with Weak Layers
is numerical model uses a two-dimensional plane stress model with a model size of 50 mm × 100 mm. As shown in Figure 5, the inclination angle of the weak interlayer is θ, the thickness is h, the number of weak layers is c, and the interlayer spacing is d.
e interface contact between layers is completely elastic regardless of the influence of   Based on the comparative analysis of laboratory experiments and the collected data, the basic physical and mechanical parameters in the numerical model are shown in Table 2. e numerical calculation is based on the occurrence of a through-sliding failure surface in the model or the nonconvergence of the calculation as the end condition.  [30] Laboratory experiment [30] 4 Advances in Materials Science and Engineering To avoid the influence of the complex bonding mechanical parameters of the interface on the numerical results, the calculation process does not consider the effect of structural surfaces between adjacent rock masses. It is assumed that the bonding between the rock masses is good [31]. To explore the uniaxial compression values under different conditions and simulate the law of strength characteristics and the creep mechanical behavior of a rock mass with weak layers, we considered the influence of the four main factors of the inclination angle, thickness, number of layers, and the distance between weak interlayers on the creep mechanical behavior of the rock mass, as shown in Table 3, the numerical simulation program.

Crack Propagation Evolution Process.
e numerical model of a rock mass with weak interlayers is composed of many elements, and the element boundary is connected by springs. Under constant pressure, tensile failure and shear failure occur when the interface contact spring reaches the maximum tensile stress and Mohr-Coulomb composite criterion. We took a rock mass with weak interlayers (θ � 30°, d � 20 mm, c � 1) as the research object to analyze the crack propagation path. In the simulation process, the initiation position, quantity, length, and propagation direction of the initial crack can be observed intuitively, and the initial state diagram and final state diagram in the deceleration creep stage, isokinetic creep stage, and accelerated creep stage can be intercepted. By comparing the front and rear pictures, we analyzed the initiation position and propagation path of a tensile crack and shear crack. is method is used to describe the crack evolution process in the creep process of rock mass with weak interlayers, wherein the direction indicated by the red arrow in the figure is the crack propagation direction.
As shown in Figure 6, at the moment the load is applied, the initial pores and cracks inside the model are compressed and closed, the model undergoes instantaneous elastic deformation, and no cracks are generated on the surface of the model. For tensile crack a and shear crack b, the length of crack a and b is about 1 mm. e extension direction of tensile crack a is perpendicular to the interface of soft and hard rock mass, and the shear crack b extends along the direction of the weak layer. When it runs to 2,600, new cracks (c, d, e, f, g, h, j) continue to sprout in the weak layer, and the propagation direction of shear cracks c, e, and h is consistent with the direction of crack b. e propagation direction of extension cracks d, f, and j is the same as the direction of crack a, the propagation direction of shear crack g is opposite to crack a, the extension crack d extends toward the lower interface, and a wing crack develops on the left side of extension crack j. When it gets to time Step 7,600, it enters the stage of constant velocity creep. e lower end of the shear crack c expands in an "S" shape along with the     Figure 7. Under constant load conditions, when the weak layer inclination angle is θ � 0°, the main mode of cracks is shear cracks. After the shear cracks develop to a certain extent, X-type conjugate shear failure occurs through the weak interlayer. When the inclination angle of the weak layer is θ � 15°, the shear cracks first extend along the direction of the weak layer and then extend through the interlayer along the loading direction, finally showing a shear slip failure trend. When the inclination angle of the weak layer is θ � 30°, the sample has a shear fracture at the interface, and sliding shear failure occurs along the interlayer direction. When the inclination angle of the weak layer is θ � 45°, the sample undergoes sliding shear failure along the weak interlayer. With the increase of the inclination angle, the failure form of the rock sample with the weak layer changes from the shear failure through the weak layer to the slip shear failure along with the weak layer.
e creep curves of rock masses with different weak layer dip angles are shown in Figure 8. When the weak layer dip angle is θ � 0°and 15°, the creep curves of the rock mass with weak layers only show the decelerating creep stage and the constant velocity creep stage. When the dip angle of the weak layer continues to increase, the rock mass gradually shows three stages of complete creep. First, in the elastic compaction stage, the pores and cracks inside the rock mass are compressed and closed at the moment of load application, and zero or very few shear cracks appear on the surface of the sample. After entering the deceleration creep stage, the strain of the rock mass increases slowly, but the creep rate gradually decreases. e change in the weak layer inclination has no obvious effect on the deformation of the rock mass. e strain of each rock mass at this stage is about 1.63 × 10 −3 . In the stage of constant velocity creep, the creep rate of each group of rock mass remains stable, and the steady-state creep rate v cp of the rock mass with weak layer dip angles (θ � 0°, 15°, 30°, 45°) is 4.32 × 10 −8 , 9.88 × 10 −9 , 7.12 × 10 −8 , and 1.56 × 10 −7 , respectively. e initial cracks inside the rock mass continue to develop and expand, and many secondary cracks evolve during this stage. Under the long-term action of external force, the cracks continue to evolve, extend, and intersect one another to form a penetrating slip surface. e rock mass has undergone significant macroscopic deformation by this point, and the total strain ε of the rock mass is 0.0024, 0.0020, 0.0034, and 0.0070 in sequence. As shown in Figure 9, the dip angle of the weak layer is one of the factors that affect the creep of the rock mass. In the range of small dip angles, the creep strain and creep rate decrease with the increase in the dip angle of the weak layer. Beyond this range, the strain increases with the increase of the dip angle of the weak layer, which conforms to the U-shaped distribution characteristics.

e Influence of the ickness of the Weak Layer on the Creep Characteristics of Rock Mass.
To deeply study the influence of the thickness of the weak interlayer (h � 10 mm, 20 mm, 30 mm, 40 mm) on the creep characteristics of the weak layered rock mass under a constant load and a given working condition (θ � 30°, c � 1), numerical simulation calculations were conducted on the samples. e crack propagation characteristics and equivalent displacement cloud diagrams of the rock masses with different weak layer thicknesses are shown in Figure 10.
Under constant load conditions, the failure of the rock mass with weak layers starts from the weak layer and then gradually expands to the hard rock on both sides. When the thickness of the weak layer is h � 10 mm, the cracks in the weak layer are mainly shear cracks, and the rock mass tends to slide along the rock formation. When the thickness of the weak layer is h � 20 mm, shear fracture occurs in the weak layer, and the sample runs along the weak layer. e layer has shear slip failure when the thickness of the weak layer is h � 30 mm and 40 mm. Local shear failure occurs on the lower left side of the weak layer, and the sample slips along the shear surface. e four failure modes of weak layercontaining rock masses with weak layer thickness always result in sliding shear failure along the weak layer shear failure surface.
Under the same working conditions, the creep curves of each weak layer of thick rock mass are shown in Figure 11. e figure indicates that the creep curves of the weak layer thickness, h � 30 mm and h � 40 mm, have similar trends. When a constant load is applied, the instantaneous strains of the rock mass with weak layers (h � 10 mm, 20 mm, 30 mm, 40 mm) are 5.02 × 10 −4 , 5.08 × 10 −4 , 5.11 × 10 −4 , and 5.29 × 10 −4 , respectively. At this stage, the initial pores and cracks inside the rock mass are closed. After entering the deceleration creep stage, the internal cracks of the rock mass gradually begin to develop and expand. e creep curve indicates that the thickness of the weak layer gradually influences the deformation and failure trend of the rock mass. In the constant velocity creep stage, the steady-state creep rate v cp of each rock mass is 1.12 × 10 −8 , 7.27 × 10 −8 , 1.99 × 10 −7 , and 2.01 × 10 −7 in sequence. After entering the accelerated creep stage, the cracks are connected and penetrated; and a shear fracture surface appears inside the weak layer. e rock eventually creeps and fails. e total strain ε of the rock mass with different weak layer thicknesses is 0.0019, 0.0034, 0.0054, and 0.0055. Figure 12 indicates that as the thickness of the weak interlayer increases, the creep rate and creep strain of the rock mass gradually increase as well, and the rock mass is more prone to instability failure. Under long-term load, the thickness of the weak interlayer has an important influence on the creep deformation of the rock mass.  Advances in Materials Science and Engineering h � 10 mm). e influence law of rupture, the rupture mode, and equivalent displacement cloud diagram are shown in Figure 13. Under constant load conditions, the cracks in the weak layer of the rock with one weak layer (c � (1) are mainly shear cracks; and the rock tends to slide and fail along the rock. e displacement value of the left side of the upper interface of the lower weak layer in the rock mass with two weak layers (c � (2) is about 0.001186 mm. e displacement change is the largest.

e Influence of the Number of Weak Layers on the Creep
ere is an obvious shear fracture surface near this spot, and the rock mass along the sliding shear failure is on the fractured surface of the rock mass. e cracks in the rock mass with three weak layers (c � (3) penetrate the weak interlayer, and the displacement value on the left side of the upper interface of the middle weak layer is 0.001613 mm, which is the largest value; thus, a local shear fracture surface is generated at the middle weak layer. e bearing capacity of the rock mass decreases, and the rock mass slips and fails along the shear plane. As the number of interlayers increases, the form of creep failure of the rock mass with weak layers changes from sliding shear failure along the weak layer to sliding shear failure through the weak interlayer. Figure 14 shows the creep curves of rock masses with different numbers of weak layers. e comparative analysis shows that when the number of weak layers is c � 1, the rock mass mainly presents a decelerating creep stage in which the creep rate gradually decreases with time. e characteristics of the constant velocity creep phase remain approximately constant with the creep rate. When the number of weak layers is c � 2 or 3, the rock mass presents a complete threestage creep process. Figure 14 shows that the steady-state creep rate v cp of each group of rock masses (c � 1, 2, 3) is 3.35 × 10 −8 , 5.86 × 10 −8 , and 6.19 × 10 −8 , respectively, and the total strain ε is 0.00162, 0.00302, and 0.00345, respectively. For rock masses with different weak layers, the total amount of instantaneous strain and creep under the same stress level is different, and both increase as the number of layers increases. According to Figure 15, the total strain and steadystate creep rate of the rock mass with weak layers increase as the number of weak layers increases. is is because of the increase in weak interlayers: the larger the volume of soft rock in the rock mass with weak layers is, the lower the overall rigidity of the rock mass will be, the greater the increase in plastic deformation will be, and the more likely it is that the rock mass will be broken. e numerical  simulation result is the same as the simulation result from the literature [33].

Influence of Weak Layer Spacing on Creep Characteristics of Rock Mass.
Considering the influence of the distance between weak layers on the creep failure of a rock mass with double weak layers under a constant load, keep other conditions unchanged (θ � 30°, h � 10 mm, c � 1), and take d � 10, 20, and 30 mm for the distance between weak layers. e three groups of rock samples with double weak layers were numerically simulated, and the simulation results are shown in Figure 16.
Under the given working conditions, the creep rupture mode of the rock mass with the weak layer spacing (d � 10 mm) is shown in this figure. e displacement value on the left side of the upper interface of the lower weak layer is about 0.001186 mm, and the displacement change is the largest. ere is an obvious shear fracture surface near this point, and the rock mass has sliding shear failure along the fracture surface. e displacement of the left unit at the upper end of the weak layer spacing (d � 20 mm) shows the largest change, with a value of 0.001544; and the upper rock mass of the lower weak layer moves downward as a whole and slips along the shear fracture surface of the destroyed lower weak layer. e rock mass with the weak layer spacing (d � 30 mm) also has the largest displacement on the left side of the upper end. When the calculation is completed, it is clear that the cracks in the lower weak layer do not intersect each other to form a through-plane, so there is no creep failure overall. e existence of a weak interlayer destroys the integrity of the rock mass. Under the long-term action of the load, cracks first occur at the weak layer. As time passes, the cracks gradually expand, connect, and penetrate the weak interlayer, forming a macroscopic shear. e fracture surface and the crack propagation speed at the hard rock are lower than those at the soft rock. erefore, when the distance between the weak layers is larger, the time required for the crack to expand and connect at the hard rock in the two weak layers is longer; and the cracks at the two interlayers are more difficult to penetrate. erefore, the larger the distance between the weak layers is, the less likely it is to undergo creep instability failure. Figure 17 shows the creep curve with different weak layer spacing. Under the same stress level, when the weak layer spacing is d � 10 mm, the instantaneous strain of the weak layer rock mass is 1.34 × 10 −4 , the steady-state creep rate v cp is 5.8 × 10 −8 , and the creep strain is 0.00358. When the distance between weak layers is d � 20 mm, the instantaneous strain of the rock mass with weak layers is 1.32 × 10 −4 , the steady-state creep rate v cp is 2.8 × 10 −8 , and the total strain ε is 0.00238. When the distance between weak layers is d � 30 mm, the instantaneous strain of the rock mass with weak layers is 1.89 × 10 −4 , the steady-state creep rate v cp is 1.6 × 10 −8 , and the total strain ε is 0.00161. As shown in Figure 18, the creep strain and steady-state creep rate of the rock mass with weak layers decrease with the increase in the interlayer spacing; and the change rule of the instantaneous strain is not obvious.

Conclusion
is article mainly describes the evolution process of local element deformation and failure combined with the CDEM method, and it numerically simulates the influence of different weak layer inclinations, thicknesses, numbers of layers, and layer spacing on the creep characteristics of rock mass under the conditions of the uniaxial compression creep test. e conclusions are as follows: (1) In CDEM, when calculating the interface failure, the judgment standard for the interface failure is based on the maximum tensile stress and the Mohr-Coulomb composite criterion. e failure modes are divided into tensile failure and shear failure.  (2) With the increase in the dip angle of the weak layer, its creep failure mode changes from shear failure through the weak layer to slip shear failure along the weak layer. When the dip angle is within a certain range, as the dip angle of the weak layer increases, the total strain and steady-state creep rate of the rock mass decrease. After exceeding the threshold, the total strain and steady-state creep rate increase with the dip angle of the weak layer, showing a U-shaped distribution. (3) Under long-term load, regardless of the thickness, the creep failure mode of the rock mass with weak layers always occurs along the fracture surface; as the thickness of the weak layer increases, the total strain of the rock mass with weak layers increases. e steady-state creep rate also increases. (4) As the number of interlayers increases, the form of creep failure changes from sliding shear failure along the weak layer to sliding shear failure through the weak interlayer. e total strain and steady-state creep rate increase with the number of weak layers.
is is because of the increase in weak interlayers. e volume of soft rock in the weak layered rock mass increases, the overall rigidity of the rock mass decreases, the plastic deformation increases, and the rock mass is more prone to damage. (5) e distance between the weak interlayers affects the crack propagation speed during the creep process: e larger the interlayer spacing, the slower the crack propagation speed. e steady-state creep rate and total strain gradually decrease, and creep failure occurs along the shear fracture surface.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.