Study on Failure Process and Permeation Evolution of Single-Cracked Rock

To evaluate the mechanical properties and permeation evolution of cracked rock mass, failure evolution tests were designed by RFPA software for single-cracked rock mass with (i) different inclination angles under uniaxial compression and (ii) different confining pressures and pore pressures under triaxial compression. 0e results show the following: (1) Angle of the crack significantly affects the crack propagationmode and slightly affects the bearing capacity of rock. During the crack propagation, the peak of permeation is delayed at the peak of stress. 0e stress-strain curve shows a different behavior in the postcritical part of the curve, especially in the case of 45°, where a smooth postcritical curve was clearly observed instead of an abrupt decrease in the stress of other two cases. (2) When the confining pressure is constant, the trend is almost the same when varying pore pressures, and with the increase in pore pressure, crack propagation is accelerated. At a low confining pressure, the crack is extended vertically to the upper and lower ends of the specimen, forming a longitudinal macroscopic crack. At a high confining pressure, the crack gradually extends to the left and right boundaries of the specimen, forming a transverse macroscopic crack. (3) 0e rate of crack initiation and destruction first decreases and then increases with the increase in confining pressure when pore pressure is constant.


Introduction
Safe and efficient green mining of coal resources is important for the national economy and social development in China.However, the geological conditions of coal mine in China are very complex, and water inrush [1][2][3] has been one of the important factors to hinder the production of coal.Water inrush is affected by many factors, for example, water abundance of aquifer, geological structure, and impedance water capacity of floor aquifer.Geological structure, especially the crack structure [4,5], is the main cause of water inrush, closely related to the permeation of rock mass.In nature, rock is a defective body with many tiny cracks, holes, or other structures.erefore, rock has the characteristics of discontinuity, heterogeneity, and inhomogeneity.Under the combined action of load and fluid pressure, preexisting weaknesses can close, open, extend, or induce new cracks, which can in turn change the structure of the rock and alter its fluid flow properties [6][7][8][9][10][11]. Further, skeleton grains of rock will shrink or expand, which will change its permeation characteristics.After the permeability of rock has increased, the confined water of floor is promoted to the area with large permeability.Finally, the channel of water inrush forms.Since the permeation of Westerly granite under a high pressure was investigated through experimental testing by Brace et al. [12], similar laboratory studies in this field have been extensively conducted in different types of rock, for example, by Shiping et al. [13], Wang and Park [14], and Oda et al. [7] in granite; by Wang et al. [15] for coal.
Some theoretical and numerical investigations on the realistic failure process analysis of rock mass have been widely reported.Eftekhari et al. [16] studied the crack propagation in the rock specimen with the inclined crack under uniaxial compression and found that the crack growth direction was deviated from the initial crack direction and propagated from the crack tip toward the loading plates.Yuan and Harrison [17] reviewed previous research and summarized that the following three main methods were employed to model the progressive mechanical breakdown of rock and associated uid ow in rock: discrete models based on crack mechanics.Xu and Yang [18] established the constitutive model for ssured rock mass based on the deformation characteristics of microcracks under compression.
Tang et al. [9] set up a ow-stress-damage (FSD) coupling model for heterogeneous rocks based on a statistical model that considers the growth of existing cracks and the formation of new cracks.is model was coded into a twodimensional realistic failure process analysis (F-RFPA2D) by Yang et al. [19] to simulate hydraulic fracturing in permeable rock, and then it was extended to three dimensions by Li et al. [20].Zhu et al. [21] studied the e ects of axial static stress and dynamic stress on the damage and failure process of rock by RFPA and found that the simulation results were closely related to the rock heterogeneity and were in qualitative agreement with that observed in laboratory SHPB tests.Chen et al. [22] carried out a series of triaxial compression tests with permeation measurements under di erent con ning pressures and pore pressure di erence coupling conditions and found that the shapes of stressstrain curves, permeation evolution curves, and failure patterns are signi cantly a ected by the con ning pressure but are only slightly a ected by the pore pressure di erence.
us, the results of previous research deepen the understanding of the realistic failure process of rock mass and provide several theoretical methods for the numerically simulating corresponding issues.
However, the relationship between failure characteristics and permeation evolution at di erent angles, con ned pressures, and pore pressures has been rarely studied.erefore, in this study, the tests for cracked rock mass at di erent inclination angles under uniaxial compression and di erent con ning pressures and pore pressures under triaxial compression were designed by RFPA software to analyze the relationship between failure characteristics and permeation evolution.e results of this study are expected to provide a valuable reference for excavation in geotechnical engineering and protection of water inrush.

Model of Single-Crack Propagation
At present, crack propagation of brittle materials occurs through (i) tension crack and (ii) slide crack.e sliding crack model indicating the crack mechanism of microcracks is gradually taking the dominant position in explaining the microcrack propagation of brittle materials and accepted by the academic community.Brace and Bombolakis [23] rst proposed the sliding crack model; then, it was popularized by Horii and Nemat-Nasser [24] to study the mechanics of crack brittle materials under compression conditions.

Single-Crack Propagation Model under Triaxial Stress.
As shown in Figure 1, for a crack with a length of 2a inside an in nite plate in the plane, the model was subjected to an axial pressure of σ 1 and a con ning pressure of σ 3 ; the angle between crack surface and axial stress is β; the internal cohesion of the crack is c w ; the internal friction angle is ψ w ; the internal cohesion of rock is c; and the internal friction angle of rock is ψ [25].e stress and shear stress of the crack can be shown as follows: According to the theory of elastic crack mechanics and judgment basis of shear crack under a pressure condition (k Π k Πc ) [26], the initial crack strength of the crack can be expressed as follows: where k Πc is the shear toughness under compression conditions, f w is the friction coe cient of crack, c w is very small when the uniaxial compression, c w 0; therefore, the initial crack strength can be expressed as follows: erefore, k Πc is the critical value of stress intensity factor for crack propagation.Tests show that when k Π > k Πc , the crack starts to expand.

Single-Crack Propagation Model under Pore Pressure.
To simplify the analysis, Song [27] simpli ed the single-crack propagation model as follows: a standard specimen with a height of 100 mm and diameter of 50 mm was used.e inside of the model has a crack with a length of l and an angle of a, and the internal crack space of the specimen is lled with pore pressure with a size of P, as shown in Figure 2.

Advances in Materials Science and Engineering
According to rock mechanics, crack tip effect and single-crack rock mass pressure shear crack criterion (λK I + K II � K IIC ) and the crack strength can be expressed as follows: where σ crack is the crack strength; when σ x ≥ σ crack , the crack began to expand.It can be seen from ( 3) and ( 4), under the action of compressive stress and pore pressure, the microcracks can close, open, extend, or induce new cracks, and the prefabricated cracks are relatively sliding under the action of shear stress.When the crack is activated, the tip is extended and connected through the concentrated effect, changing the permeation and strength of rock mass.erefore, it is necessary to analyze the stress model and crack strength of cracked rock mass under stress and pore pressure.

Experimental
In this study, the numerical simulation tests were divided into two groups: (i) the crack characteristics and permeation evolution test of the single crack in different inclination angles under uniaxial compression and (ii) the relationship between failure characteristics and different confining pressures and pore pressures in single-crack specimens.
In the first numerical simulation, a numerical model was established based on the constitutive model of the Mohr-Coulomb mechanics.e lengths of the model's X and Y axes were 0.05 and 0.10 m, respectively, with a total of 45,000 units.e model was calculated by plane strain.e inside of the model had a crack with a length of 0.02 m and width of 0.001 m. e angle between the crack and horizontal direction is known as the crack angle β, and its values were 15 °, 30 °, 45 °, 60 °, and 75 °, as shown in Figure 3. e model used the standard model uniaxial compression base load.e loading type was the stress control.e initial value was 0, and the single step increment was 0.6 MPa.A pore pressure of 0.8 MPa was applied to the hole of crack cavity, and the water head of model's top, bottom, left, and right was 0. e mechanical parameters of materials used in the model are shown in Table 1.
In the second numerical simulation, to study the influence of different confining pressures and different pore pressures on the failure process and permeation evolution of single-cracked rock, avoiding the impact of the crack angle on the failure process, the angle of the crack significantly affects the crack propagation mode and slightly affects the bearing capacity of rock; therefore, β � 45 °was selected for simulation in this experiment.
is simulation uses the confining pressure compression under the standard mode.To verify the effect of different pore pressures and confining pressures on rock crack characteristics, pore pressures and confining pressures of this test were set as 1.0, 1.5, 1.8 MPa and 2, 3, 5 MPa, respectively, as shown in Figure 4. e other boundary conditions and mechanical parameters were the same as the first simulation experiment.5 shows that the angle of the crack significantly affects the crack propagation mode with the increase in the crack angle, and the main crack is divided into two cracks and finally intersects with the main crack during the crack propagation.e main trends are as follows:

Results Analysis and Discussion
(1) When the crack angle is small (β �15 °), the crack propagation does not occur in the crack tip but takes place near the crack tip.e crack propagation is relatively easy before the stress reaches the peak, but the crack initiation requires a larger stress.( 2) With the increase in the crack angle, the crack propagation mainly occurs at the crack tip expansion.e stress required for crack propagation first decreases and then increases with the increase in the crack angle.e direction of crack propagation is approximately perpendicular to the extension line of the main crack.As the axial pressure continues to increase, the energy starts to transfer to the main crack.e cracks are then compacted, and a stress concentration occurs at the tip of the crack parallel to the axial stress direction.
e new crack is extended to the tip of the specimen, resulting in macroscopic damage to the rock specimen.(3) When β � 75 °, the crack angle is almost parallel to the axial loading direction, and the direction of crack propagation is parallel to the direction of main cracks.At this time, the crack initiation is easy, but the crack propagation is strongly inhibited, increasing its peak strength.In conclusion, with the increase in the crack angle, the strength of rock first decreases and then increases, consistent with that reported by Liang et al. [28].According to the preliminary analysis of the numerical simulation results, when β � 45 °, the strength of rock may be the lowest.
According to the evolution of pore pressure, the trend of permeation evolution is consistent with that of mechanical properties.Pore pressure of model instability was 0.79 MPa; therefore, the crack angle slightly affects the bearing capacity of rock.Pore pressure in vertical (Y) evolution degree is higher than the horizontal (X) evolution degree, and the permeation of model in the longitudinal direction is higher than the horizontal permeation.

Permeation Evolution under Uniaxial Compression.
e analysis indicates that with the increase in the crack angle, the strength of rock first decreases and then increases, ranging from 30 °to 60 °.erefore, to study this process, the stress-strain and flow-strain curves of 30 °, 45 °, and 60 °were selected as the research object, as shown in Figure 6.

Advances in Materials Science and Engineering
Figure 6 shows that the stress-strain curves were subjected to compaction, elastic, plastic, and residual load stages.In the early stage of axial stress loading, the original crack did not develop, and the permeation of the specimen was small.As the axial stress continued to increase, the specimen entered the plastic failure stage.e specimens are constantly destroyed, and the crack constantly sprouted, expanded, and penetrated; therefore, the permeation constantly increased.e stress-strain curve shows that after entering the plastic zone, pressure drop occurred in the three di erent angles, and the crack developed.Small cracks joined, and a big crack formed.e ow-strain curve shows that the permeation evolution and stress-strain trends are consistent, but the peak of permeation delayed at the peak of stress.During the residual load stage, the specimens can still bear axial stress, and the permeation rst decreases and then maintains a higher level or continues to increase.
In addition, when β 45 °, the stress-strain curve is di erent from that of β 30 °and 60 °, and the stress-strain curve shows a di erent behavior in the postcritical part of the curve, especially in the case of 45 °, where a smooth postcritical curve was clearly observed instead of an abrupt decrease in the stress of other two cases (β 30 °and 60 °).
is is because when the strength of rock decreases with the increase in strain, additional load must be delivered to the sample, as the rocks show residual strength in the postdestructive region.In the other case, the destruction process is uncontrollable, and the piston of the material testing machine must be withdrawn to stop it [29,30].In addition, according to the e ect of the single joint angle on the strength of the single-crack rock specimen, when β satis es the following (5), sliding failure along the joint plane direction occurs in the single-crack rock, and the corresponding strength is the lowest [31].At this time, the failure of the specimen is mainly caused by the shear failure.
erefore, the stress-strain curve shows a di erent behavior in the postcritical part of the curve, especially in the case of 45 °, where a smooth postcritical curve was clearly observed instead of an abrupt decrease in the stress of other two cases: where φ is the internal friction angle of the joint plane and μ is the friction coe cient of the joint plane, which is 0 in this study.us when β is 45 °, the uniaxial compressive strength of the single-crack rock specimen is the lowest.Subsequently, in this study, the uniaxial compressive strength of the rock is the lowest, and the sliding failure occurs along the joint plane within the rock body.

Effect of Pore Pressure and Confining Pressure on Failure Characteristics and Permeation Evolution of Single Crack.
To simplify the analysis of crack propagation trend, the crack initiation and final failure of the specimen were studied in this test.Figures 7-9 show that with the increase in confining pressure, the velocity of crack initiation to failure did not increase monotonically but first decreased and then increased.However, with the increase in pore pressure, the velocity of crack initiation to failure clearly increased.A detailed analysis of this phenomenon is given below: (1) When the confining pressure is constant, as the pore pressure increases, the velocity of crack initiation to failure clearly increased.is is because an increase in pore pressure leads to severe damage to the integrity of material inside the specimen.e crack remains on the base of the original damage and continues to grow until the specimen is eventually destroyed.e higher the pore pressure, the higher the damage to the specimen, shortening the time of the specimen being destroyed.An increase in pore pressure contributes to the development of cracks.In addition, a large difference was observed between the low and high confining pressures in the crack propagation mode.In the low confining pressure stage (σ � 2.0 MPa and 3.0 MPa), the mode of crack propagation is the same as that in the previous study, that is, with the increase in axial compression, the crack tip wings start to crack and propagate in reverse direction, and the crack propagation line is at an angle and then gradually extended.Axial stress is parallel to the axial load applied to the main crack.A "wing crack" forms and vertically extends to the specimen on the lower end surface, forming a longitudinal crack.e process of crack propagation is accompanied by secondary antiwing and coplanar cracks, but the direction is approximately parallel to the axial direction.In the high confining pressure stage (σ � 5.0 MPa), the crack propagates at a certain angle; then, with the increase in axial pressure, the propagation direction is perpendicular to the axial stress.Finally, the crack gradually expands to the left and right boundaries of the specimen with the form of "antiwing crack," forming transverse macrocracks.
is is because when the confining pressure plays a major role in the coupling process, the cracks are under confining pressure.Homogeneous specimens in the horizontal direction are significantly lower than the axial direction, and the specimens are damaged in the horizontal direction.erefore, the direction of crack propagation changes.e cracks are gradually transformed to "antiwing cracks" form the level extended to the specimen around the border.
(2) When the pore pressure is constant, the velocity of crack initiation to failure does not increase monotonically but first decreases and then increases.is is because during the crack development, (i) new cracks are produced, (ii) original cracks disappear, (iii) cracks propagate, and (iv) cracks interpenetrate.e four stages do not exist independently but exist together and undergo mutual transformation.Under the action of con ning pressure, the small cracks are compacted, leading to the disappearance of the original crack.
e homogeneity of specimens increased, and the integrity of specimens improved, leading to crack propagation.When σ 3.0 MPa, the crack propagation is more than that when σ 2.0 MPa.With the increase in con ning pressure, the specimen severely damaged.e number of original cracks that disappeared is much smaller than the number of new cracks, and the speed of crack propagation is higher than that of disappearance of original cracks, leading to a much higher rate of crack propagation.
(3) e trend of permeation evolution is consistent with the direction of crack propagation.In the low con ning pressure stage (σ 2.0 MPa and 3 MPa), the pore pressure produces cracks in the "elliptical" longitudinal direction, extending to the specimen under the surface, and the longitudinal evolution degree is greater than the lateral evolution degree.ese two factors rapidly increase the permeation of the specimen in the longitudinal direction (Y direction), much higher than transverse permeation (X direction).In the high con ning pressure stage (σ 5.0 MPa), pore pressure extends at the level of feather shape, and the degree of lateral evolution is higher than that of longitudinal evolution.is leads to a sharp increase in the permeation of specimens in the transverse direction (X direction), much higher than that in the longitudinal direction (Y direction).

Conclusions
In this study, tests were designed by RFPA software for cracked rock mass with (i) di erent inclination angles under (1) Crack angle significantly affects the crack propagation mode with the increase in the crack angle and slightly affects the bearing capacity of rock; it can be divided into three modes: (i) the crack propagation does not appear in the crack tip, (ii) the direction of  crack propagation is approximately perpendicular to the extension line of the main crack, and (iii) the direction of crack propagation is parallel to the direction of the main cracks.In addition, the trend of pore pressure evolution is consistent with that of crack propagation.(2) With the increase in the crack angle, the strength of rock first decreases and then increases.e peak of permeation is delayed at the peak of stress.e stress-strain curve shows a different behavior in the postcritical part of the curve, especially in the case of 45 °, where a smooth postcritical curve was clearly observed instead of an abrupt decrease in the stress of other two cases (β � 30 °and 60 °).
(3) When the confining pressure is constant, as pore pressure increases, the velocity of crack initiation to failure clearly increased.ere is a significant difference between the low and high confining pressures in the crack propagation mode.(4) When the pore pressure is constant, the velocity of crack initiation to failure does not increase monotonically but first decreases and then increases.e trend of permeation evolution is consistent with the direction of crack propagation.
erefore, the results of this study are expected to provide a valuable reference for excavation in geotechnical engineering and protection of water inrush.In addition, we will test the simulation results in the future and expect to find more discoveries on this basis.

Figure 3 :
Figure 3: Model schematic diagram and RFPA model diagram under uniaxial compression.

Table 1 :Figure 2 :
Figure 2: Compressive shear model of rock mass with the single crack under pore pressure.