Study of Simulation Test in Inclusion Rockburst

Based on the existing classification research of rockburst dynamic disasters, a new rockburst model is introduced. -e model specimens with inclusion rock are designed and produced, and the indoor test of inclusion rockburst with or without cementation, different inclusion angles, and prefabricated cracks of different lengths is carried out. -e failure process and displacement variation of rock specimens are recorded and measured by the high-speed camera system and digital image correlation method, and based on this, the characteristics and laws of rockburst are analyzed.-e results show that the failure process of inclusion rock can be divided into several periods, such as crushing and extrusion of the cement body, and the cracks in the upper left tip, the lower right tip, the lower left tip, and the upper tip are expanded and broken in different orders.-e dynamic ejection failure of the rock fragments occurs when the degree of failure is severe, and the static failure occurs when the degree of failure is weak. -e inclusions have cementation, and the better the cementation, the greater the stress loading required in the failure, the greater the energy accumulated, the more severe the failure, and the higher the fractal dimension of the fragments. -e larger the angle of the inclusion tip, the greater the stress loading in failure; the greater the change of the u and v displacements, the more the displacement of the inclusions changed; and the more the energy accumulated, the higher the fractal dimension of the fragment. -e smaller the crack length, the greater the stress loading in failure, and the more the energy accumulated, the higher the fractal dimension of the fragments in failure. With more number of cracks and the more severe expansion before failure, the failure of rock mass can be transformed from the dynamic to the static.


Introduction
Rockburst is a common rock dynamic failure phenomenon in underground engineering and mining engineering [1][2][3], which is featured by huge sound, air waves, and spontaneous rock stripping and quick rock ejection, bringing casualties, damage of construction equipment, and underground engineering scrap.erefore, the rockburst belongs to a dynamically static geological disaster.
e classification of dynamic disasters on rockburst is one of the most fundamental and key scientific issues on rockburst research.Many experts have classified rockburst from the mechanism, location, and failure manifestation.Based on five international authoritative scholars on the mechanism and definition of rockburst, Qian and Zhou [4] proposed that the rockburst should be divided into fault slip or fault slip and rock strain types caused by rock failure.Based on the energy accumulation and transformation characteristics of rockburst, He et al. [5] proposed that rockburst should be divided into the types of single energy induced and composite energy conversion induced, and the single energy type can be divided into the type of solid energy induced, gas energy induced, liquid energy induced, roof falling energy induced, and tectonic energy induced.Based on the caused sudden instability, failure of coal and rock mass from the stress state, the common geological structure, and mining process characteristics, Jiang et al. [6] proposed the rockburst of coal mine can be divided into material instability, slip displacement, and structural instability.Xu and Wang [7] divided rockburst into self-heavy stress, tectonic stress, variable stress, and comprehensive stress from the cause of high ground stress on rock mass of rockburst.
In order to study the mechanism of rockburst, the relevant experts have carried out a large number of laboratory experiments on the basis of the existing rockburst classification and have obtained many valuable research results.Xie et al. [8] pointed out that the energy release of rock structure after energy dissipation was studied from the perspective of the macroscopic energy conservation, which can be used to estimate the intensity of damage in the project or the ejection velocity of the fractured rock and quantitatively analyze the energy transfer process in the rock structure.Li et al. [9] carried out a spallation failure test on granite specimens by the semisinusoidal waveform generated with the modified Hopkinson pressure bar device, and the crack failure process of the specimens is recorded by the high-speed camera.Zhou et al. [10] carried out the buckling rockburst simulation test by a similar rock material and analyzed the failure process and characteristics of the plate buckling rock.Dyskin and Germanovich [11] analyzed the formation causes of slabtype rockburst, it is believed that the stress concentration of surrounding rock caused by excavation promotes the stable expansion of primary fractures, the vacant surface causes the crack to expand unsteadily along the direction parallel to the rock wall, and the rock plate fracture and the surrounding rock are suddenly separated and ejected out.Li et al. [12] carried out a similar material test to simulate the rockburst caused by tunnel excavation, indicating that the tunnel in the field is dominated by spalling and drum-type tensile failure.Gong et al. [13] carried out a three-dimensional loading test on specimens with prefabricated holes by a true triaxial rock test system and reproduced the whole process of plate buckling rockburst in a deep buried hard rock circular tunnel.Su et al. [14] carried out the physical test of the strained rockburst on the cuboid granite by the true triaxial test system and analyzed the characteristics and laws of the failure process on rockburst.Zuo et al. [15] studied the velocity of the failure block in the buckling rockburst under static loading and dynamic-static combined loading from the perspective of energy balance.Zuo et al. [16] studied the deformation and failure characteristics of different coal-rock combinations and the association with dynamic disasters.Sun et al. [17] analyzed the temperature field and displacement field evolution law and model instability failure mechanism during the test of a weak broken surrounding rock tunnel by the noncontact monitoring method.
e above research results have greatly enriched the existing occurrence mechanism of various types of rockburst and have important significance for the prevention and control of rockburst.However, the mechanism of rockburst is extremely complicated, and there is no unified and clear theoretical explanation for rockburst.erefore, based on the engineering site and other existing inclusion materials research [18,19], the model of inclusion rockburst is proposed, and the model specimen with inclusions is designed and processed.e indoor physical test on ejection process of inclusion rockburst is carried out by the high-speed camera and digital image correlation method (DIC), which realizes the real-time monitoring and recording of the whole process of the test, reproduces the rockburst process of the inclusion rockburst, and the test process and results are analyzed in detail so as to study the mechanism of the inclusion rockburst.

Inclusion Rockburst
e author found that the rock fragments on the roadway surface were exfoliated and ejected out in the ventilation roadway of the working face 0092 of Jinggezhuang mine of Kailuan Group, and the rock fragments that were peeled off and ejected to the ground were not the same rock type (siltstone and maroon mudstone), the failure craters and nearby areas formed on the roadway surface are also divided into two color regions, a small part of areas is brownish red (mudstone) and other rocky areas is gray (siltstone), and a small part of the area and the boundary part are more seriously damaged.Based on the geological data of the mine and related literature, the rock mass structure is an inclusion or inclusion layer, and the research on failure of inclusion rock in mining and rock engineering is lesser; however, a few experts have carried out some research on the inclusion problems of other materials.erefore, the inclusion rockburst model is proposed on the study of rockburst in this paper.In other words, the inhomogeneity of stratigraphic sedimentation and complex geological effects result in the nonuniform monolithic rock mass in the complex rock mass structure.Different rock geological bodies are cemented together by cement and compressed by the sedimentary layer to form inclusions, which are characterized by small volume of geological bodies being wrapped by larger volume of geological bodies, and it shows a small volume of geological bodies is wrapped by a larger volume of geological bodies.When the excavation of rock mass is close to the inclusion rock mass, the stress concentration near the rock and coal surface appears.When the rock nearby inclusions continue to be excavated, the pressure is more concentrated under the shear stress, and a part of inclusions and the cement are crushed, which lose the cementation effect so that the inclusions are crushed and sheared and ejected out, resulting in the occurrence of rockburst (Figure 1).
erefore, the inclusion model specimens are designed and processed, considering the angles of the inclusions and prefabricated cracks in rock specimens.e model test is used to simulate the occurrence process and phenomena of inclusion rockburst in rock mass engineering, the digital image correlation method is used to study the deformation field of the surface of inclusion rock, and the high-speed camera is used to obtain the failure law of the surface of the specimens and the fragments in the ejection, which explores the mechanism of inclusion rockburst, triaxial rock test system of State Key Laboratory of Coal Resources and Safe Mining (Figure 2).

Parameters Test of the Rock.
e rock sample is derived from the shale with good uniformity and complete integrity in the diversion tunnel of Longtan Hydropower Station.e rock material is processed into six standard cylindrical 2 Advances in Civil Engineering specimens of φ50 mm × 100 mm.e triaxial rock test system (GCTS) is used to carry out the uniaxial compression experiment on the specimens through displacement loading; the maximum tonnage of the testing machine is 100 t, and the loading speed is 0.05 mm/min.e mechanical parameters of the rock material can be obtained in the stressstrain curves of specimens, and the conventional calculation formulas, such as incremental method, and all of which are shown in Table 1.
Based on the entire stress-strain curve and the loading and unloading curves of the conventional specimens of shale measured by the triaxial testing machine, the relevant parameters obtained by the calculation formula can be used to judge the rockburst proneness of the shale, including elastic deformation energy index W ET , burst energy index K cf , brittleness index B, and maximum elastic strain energy E s [20][21][22][23].

Elastic Deformation Energy Index W et .
When the shale specimen is subjected to 80% of the peak strength in the uniaxial compression test, the stress-strain curve can be obtained.us, the elastic deformation energy E e and the plastic deformation energy E p are obtained, and elastic deformation energy index W et can be obtained from the following formula: Criteria for rockburst proneness: W et : W et ≥ 5.0, strong rockburst proneness; 3.5 ≤ W et < 5.0, medium rockburst proneness; 2.0 ≤ W et < 3.5, weak rockburst proneness; and W et < 2.0, no rockburst proneness.

Burst Energy Index K cf .
Based on the whole process curve of stress and strain of shale, the area F 1 surrounded by the curve and the strain axis before the stress peak and the area F 2 surrounded by the curve and the strain axis after the stress peak are, respectively, calculated, and the burst energy index K cf can be obtained from the following formula: Criteria for rockburst proneness: W et : K cf ≥ 3.0, strong rockburst proneness; 2.0 ≤ K cf < 3.0, medium rockburst proneness; 1.0 ≤ K cf < 2.0, weak rockburst proneness; and K cf < 1.0, no rockburst proneness.

Brittleness Index B.
e brittleness index B is the ratio of the uniaxial compressive strength σ c to the tensile strength σ t , which re ects the brittleness degree of the rock.e larger the brittleness index, the more likely the rockburst occurs:

Maximum Elastic Strain Energy E s .
e maximum elastic strain energy E s stored in the standard specimen of shale is as shown in the following formula, where E is the elastic modulus of the rock: Criteria for rockburst proneness are given as follows: E s : E s < 0.20 MJ/m 3 , no rockburst proneness; 0.2 ≤ E s < 0.5 MJ/ m 3 , weak rockburst proneness; 0.5 ≤ E s < 0.75 MJ/m 3 , medium rockburst proneness; and E s ≥ 0.75 MJ/m 3 , strong rockburst proneness.erefore, based on the four rockburst proneness indicators and the discriminant criteria obtained from the above calculations, the preliminary rockburst proneness of shale is shown in Table 1.

Preparation of Rock Samples.
e shale sample is processed into the cuboid specimen with an overall size of 76 mm × 50 mm × 30 mm. e symmetrical triangular cylinder is cut in the middle of every rock specimen.Since the inclusion rockburst model is proposed for the first time, the tests are carried out twice.
e first batch is a tentative simulation test, and the second test is based on the first test to optimize the test.erefore, in the first test, the angle of the tip of the triangular cylinder (inclusion body) of the rock specimen is 20 °and the tip side is 12 mm from the left and right side.e same has been proceeded in two specimens, which are, respectively, labeled as specimen 1# and specimen 2# (Figure 3(a)).e second batch of specimens is processed with four shale specimens, namely, 3#, 4#, 5#, and 6#.Four specimens are prefabricated symmetric cracks; the outer crack tip is 11 mm from the bottom, and the angle between the crack and the horizontal is 30 °.In order to study the effect of crack length a on the failure of the specimen, the crack length of specimen 3# and specimen 6# is 6.5 mm, and the crack length of specimen 4# and specimen 5# is 5.5 mm.Similarly, in order to analyze the effect of the inclusion angle θ on the failure, the tip angle θ of specimens 3# and 4# is 40 °, and the tip angle θ of specimens 5# and 6# is 30 °, which is shown in Figure 3(b) In order to understand the effect of cement on the inclusion rockburst, the two triangular cylinders (inclusions) are put directly in the triangle areas of specimen 1# without cementation (Figure 4).e two triangular cylinders (inclusions) are cemented in the cut triangle area of specimens 2#, 3#, 4#, 5#, and 6#, which are cured for one week (Figure 3(c)).e cement material is Anda 42.5 grade fast hard sulphoaluminate cement, and the compressive strength of cement is 53.3 MPa, the modulus of elasticity is 0.22 Gpa, Poisson's ratio is 0.16, and the water-cement ratio is 0.6 [24].
e physical and mechanical parameters of the shale are shown in Table 1.

e Layout of Test System.
e uniaxial compression test of specimens 1#, 2#, 3#, 4#, 5#, 6# is carried out by the GCTS triaxial rock test system.Since the diameter φ of the indenter and base of the system is 50 mm, two discs of the same stiffness are placed between the specimen and the indenter and base of the test machine, the diameter φ of the disc is 100 mm, and the stiffness of the disc is 1.75 mN/mm.e rate of displacement loading is 0.05 mm/min.Since the first test is a tentative test, the video recording of the specimens 1# and 2# is performed only by a normal camera.
e uniaxial compression test of specimens 3#, 4#, 5#, and 6# is carried out at the same displacement loading rate in the second test; the high-speed camera FASTCAM1.1 and the digital camera MV-VD500SM/SC are used to synchronously record the deformation process of the specimen in the loading process.e photo interval of the high-speed camera is 0.001333 s, and the time interval of the digital camera is 1 s.
e test system is shown in Figure 5.

e Result of Rockburst Test.
Due to the relatively high rigidity of the disc pad, the pad hardly affects the result of specimens in this test.erefore, Figure 4 shows the whole process stress-strain curve of specimens 1#, 2#, 3#, 4#, 5#, and 6#.Ejection failure occurred in specimens 1#, 2#, 4#, 5#, and 6#, and static brittle failure occurred in specimen 3#. e prepeak curves of the specimens 1#, 2#, 3#, 4#, and 5# have obvious compaction stages, and specimen 6# directly skips the compaction stage on account of the test operation factors in the typical stress-strain curve.e fluctuations of curve in the prepeak stress correspond to the failure phenomenon in the test, respectively.e crack initiation, expansion failure, ejection of the penetration region and inclusions, shear failure of the specimen as a whole from the tip of the inclusion body, and sounds of different sizes are accompanied.
After reaching the peak stress, the curve drops sharply, and Advances in Civil Engineering the curve is steep, which is fully manifested as brittle failure of the specimen.
In order to describe in detail the occurrence process of the inclusion rockburst test, it is analyzed and shown as Figure 4 and Table 2.

Test of Specimen 1#.
Since the test of specimen 1# is tentative, it is only recorded by a normal camera, and the resolution and time interval of the captured image are low, and the typical failure process is as shown in Figure 6. e color of the tip gradually becomes white with the increase of the load, which indicates that the part is concentrated in stress and the damage of the tip is obvious.When the axial stress σ z is 88.39 MPa, the stress-strain curve shows a drop uctuation.e tip of the right inclusion of the specimen is cracked at this time, and the reverse crack with a crack angle of 95 °is extended downward to the bottom of the specimen.
e cracked rock mass and the right inclusion are rapidly ejected (Figure 6(b)), accompanying with a crisp sound.e white part of the left and right inclusions became more and more obvious during the loading process, and the new defects gradually became produced.When the stress σ z is 93.42 MPa, the curve showed small uctuations, and the right tip showed an airfoil crack with a crack angle of 70 °50 0.020 0.025 0.030 0.035 e stress loading continues to last for 6 seconds.When the peak stress σ z of 98.22 MPa is reached, the airfoil crack gradually expands to the bottom of the specimen, and the rock mass formed by the airfoil crack and the reverse crack is ejected again, accompanying with a loud sound.e left tip is cracked in the angle of 90 °at the same time, and when it is extended upward by 11.57mm, the vertical bifurcation appears to the right and extends to the top of the specimen and the main shear fracture zone, respectively.Finally, the overall shear failure of the specimen occurs, which is shown in Figure 6(d).Two large rockbursts occurred on the right side of specimen, and small particles and flaky projectiles occurred locally during the entire failure process.

Test of Specimen 2#.
e test of specimen 2# is also tentative.Based on the concentrated force and obvious damage, the color gradually becomes white during the loading process.When the axial stress σ z is 83.88 MPa, the tip of the right inclusion is cracked, and the white reverse crack 1# with a crack angle of 85 °is expanded downward (Figure 7(a)).e white part of the left and right inclusions became more and more obvious during the loading process.When the peak stress σ z is 124.46MPa, a second airfoil crack 2# with a cracking angle of 70 °appears on the right tip, and the airfoil crack 3# with a cracking angle of 80 °appears on the left tip (Figure 7(b)).e internal damage zone results in a decrease in the overall load carrying capacity of the specimen; however, the overall failure appears on the specimen.e curve continues to rise from the peak stress 119.76 MPa. e stress σ z is 121.49MPa after 28 seconds of continuous loading; crack 4# with the crack angle of 65 °appears on the left tip (Figure 7(c)), and it gradually expands to the top of the specimen.e video shows that the left inclusion is ejected out, the connect of cracks 1# and 4# makes the upper left part of the specimen eject out, and the connect of cracks 2# and 3# makes the lower right rock and the right inclusion eject out (Figure 7(d)).Finally, the specimen is subjected to overall shear failure based on the shear band.
It can be found that the process of inclusion rockburst can be summarized that the stress concentration of the tip part of the inclusion body results in serious damage, secondary cracks appearing and expanding up and down; inclusions and the nearby cracked rock mass are ejected out   When the axial stress σ z is 111.27MPa, the curve shows a drop state with sound, and the crack 5# at the left tip with a cracking angle of 20 °gradually expands to the bottom of the specimen with an arc line, and the distance from the end of the crack to the end of the left bottom is 14 mm (Figure 8(c)).
e peak stress σ z is 112.23 MPa at loading, the failure of specimen is accompanied by a large violent sound, and the cement between the left and right inclusions and the surrounding rock is crushed (white).Based on cracks 3#, 4#, and 5#, an arc-shaped crack 6# appears.e end of the cutthrough crack at the bottom is 19.24 mm from the end point of the left bottom of the specimen.e cut-through area of the left inclusion body and cracks 3#, 4#, 5#, and 6# only appears as static failure and is not ejected out; the reason is that more cracks in this part appear, the accumulated strain energy is mostly consumed by the surface energy during the fracture process, and the remaining energy is not enough to be converted into the projectile kinetic energy.Vertical upward crack 7# and crack 8# sheared in the specimen appear in the right end of the inclusion body, crack 8# combines with crack 6# and 7#, and the surface energy of the fracture of the right part of the specimen also greatly consumes the strain energy accumulated during the loading process, resulting in specimen 3# only undergoing static brittle failure without rockburst (Figure 8(d)).
Since the energy accumulated in specimen 3# is consumed by the surface energy of more cracks, no violent rock eruption occurs, and it is only accompanied by the ejection of a small number of small rock particles on the surface and the cement between the inclusions and the specimen in initiation and cracking of the crack (Figure 9).In order to assess the threat to the surrounding rock during rockburst, the extent of the failure is characterized by the ejection velocity of the rock fragments.
e continuous photos taken by the high-speed camera are imported into MATLAB to obtain the displacement of the rock fragments ejected in different photos.
e time difference between the photos is known, so the moving displacement of the fragments is divided by the time difference to calculate the ejection velocity of the fragments.
e ejection velocities of the fragments during the ejection has been constantly changing and decreasing, so the extent of the ejection failure is characterized by calculating the initial velocities.Since there are no ejection failures of the rock fragments and only the ejection of the cement fragments in specimen 3#.erefore, based of Figures 9(a e curve continues to rise after it drops to 106.64 MPa.When the peak stress σ z is 116.627MPa, the left front prefabricated crack combines with the secondary crack generated by the tip of the inclusion, and the part of the rock together with the inclusion is broken and ejected out (Figure 10), accompanied by a crisp sound.
e shear failure appears on specimen 4# and grows to the end of left preformed crack.
e left front inclusion is cemented together with the rock in the initial stage of the ejection, the two gradually split, and some white and damaged fragments of the internal part are successively ejected out from the failure area during the ejection process.At the same time, similar failure occurs in the right rear of the specimen; the difference is that the right rear inclusion and the broken rock are separated and ejected out at the initial stage of the rockburst.e crack generated at the tip of the right inclusion is expanded by 11.36 mm at the angle of 85 °and combines with the shear band of the specimen.Many different particle flakes and fragments are ejected out from shear fracture surface of the specimen.
Similarly, the ejector velocity of the left inclusion is 2.4235 m/s through the MATLAB software, the ejection velocity of the rock block cut-through by the left lower fracture and the tip of the inclusion is 3.3373 m/s, the ejection velocity of the upper right cement is 2.3857 m/s, the ejection velocity of the lower right cement is 2.1736 m/ s, the ejection velocity of the right inclusion is 3.9044 m/s, and the ejection velocity of the rock broken at the right rear is 2.5137 m/s.Since the mass of the four blocks is When the peak stress σ z is 116.137MPa, crack 1# extends to the left inclusion, which is 5.45 mm from the left tip.At the same time, crack 4# with a cracking angle of 75 °appears at the tip of the left inclusion and extends upward to the top of the specimen.e secondary crack 2# from the tip of the right inclusion continues to extend to the bottom of the specimen, which is 6.36 mm from the lower right end of the specimen.e upward crack 5# with a cracking angle of 50 °occurs at the end of the lower right prefabricated crack and grows to the surface of the test piece (Figure 11(a)).e shear-pull crack zone 6# appears in the middle of the specimen and combines with other cracks to cause severe failure, accompanying with severe noise.e left inclusion is completely extruded and ejected out, and the half tip of the right inclusion is ejected out in rockburst.A lot of white fragments in different sizes inside slip out, during the opening of the intermediate main shear crack, indicating that a large amount of energy is accumulated inside the specimen and used for rock failure (Figure 11).
Similarly, the ejection velocity of the left inclusion is 7.2873 m/s through the preliminary calculation of MATLAB software, the ejection velocity of the rock cut-through at the lower left crack and the tip of the inclusion is 7.0085 m/s, the ejection velocity of the upper crack of the left inclusion is 6.5037 m/s, the ejection velocity of the nearby internal rock mass is 3.2752 m/s, the ejection velocity of the fractured object on the right side is 6.4197 m/s, and the ejection velocity of the rock damaged at the lower right precut crack is 5.2484 m/s.

Test of Specimen 6#.
When the axial stress σ z is 90.97 MPa in the loading of the specimen 6#, the stressstrain curve drops to 89.32 MPa, and crack 1# appears at the end of the lower left preformed crack of the specimen and gradually expands upward at a crack angle of 50 °.When the peak stress σ z is 111.074MPa, the curve shows a wave-fold increase, and the cement between the left and right inclusions and the rock is crushed and extruded, and crack 2# with an upwards cracking angle of 64 °appears from the tip of the right inclusion.At the same time, the reverse crack 3# with the downward angle of 91 °appears at the tip and combines with the lower right preformed crack.e crack 4# with a cracking angle of 80 °upward appears at the tip of the left inclusion.e secondary crack 5# with a cracking angle of 40 °appears at the end of the left prefabricated crack, extending to the tip of the left and left inclusions.e left and right inclusions are, respectively, ejected out, and the rock blocks formed after each crack penetration are also ejected out.Due to the shear tensile stress, the shear cracks 2# and 3# are torn (Figure 12), and the damaged fragments in the cracks gradually collapse.Similarly, the ejection velocity of the left inclusion of the specimen 6# is 7.4125 m/s, the velocity of the small rock near the left inclusion is 7.934 m/s, the ejection velocity of the rock block cut-through by the lower left crack and the tip of the inclusion is 5.3595 m/s, the ejection velocity of the right inclusion is 3.6242 m/s, the velocity of the rock block 1 at the lower right prefabricated crack is 5.8814 m/s, and the ejection velocity of the rock mass 2 is 6.6238 m/s.

Failure Analysis of Specimens. Figures 6∼12 show the
failure and ejection process of the specimens in the uniaxial compression test.e stress-strain curves, the shape of the fragments, and the fractures of the specimens are, respectively, shown in Figures 5 and 13.
Rock fragments in different sizes are produced in failure of inclusion rock model failure, which are accompanied by a large amount of energy consumption in the process of generating rock fragments.Since the lumpiness of rock fragments has statistical self-similarity, fractal dimension is an index to measure the characteristics of rockburst.e  2), and the fractal dimension is calculated by the mass M, the equivalent length L eq of the fragment [25][26][27], which is given as follows: where k is the slope of M L eq /M − L eq in double logarithmic coordinates, M L eq /M is the cumulative percentage of rock fragments mass in the equivalent length <L eq , M L eq is the mass of rock fragments in equivalent length L eq , and M is the total mass of the rock fragments in the calculation scale.e equivalent side length L eq is the side length converted into a cube by the length, width, and thickness of the rock measured.e relationship curve of lgM L eq /M − lg L eq is plotted, and the linear slopes of the six inclusions specimens are obtained by linear regression fitting, and the fractal dimensions D � 3 − k are obtained.e calculation results are shown in Figure 14 and Table 3.
A large amount of energy is consumed in failure of the inclusion rock, which causes a large degree of failure, and the excess energy is provided to the ejection of rock fragments.
e fractal dimension of rock fragments formed is related to the complexity of failure.erefore, the fractal dimension can measure the complexity of the rock failure process; the larger the fractal dimension, the more the number of fragments, the higher the failure degree of the rock sample, and the more the energy consumption before the severe failure.Advances in Civil Engineering 12 Advances in Civil Engineering From the results of the above-mentioned inclusion rockburst model test, it can be seen that the pinch-o portion (tip) of the inclusions is severely damaged in uniaxial compression test, and cracks appear and expand continuously, causing the rock fragments formed by the cracks to eject out, accompanied by cracking sounds.Rock fragments in the surface and inside failure area of the sample are broken and ejected out (white) in the process, and the entire sample is completely sheared and broken.
(1) Based on the failure process and characteristics of specimens 1# and 2#, it can be found that the tip of the inclusion body is seriously damaged (turns white), and cracks appear on the basis of this and expand to the periphery, and the area cut-through is large.e peak stress, strain, and fractal dimension of specimen 2# with the cementation of inclusions and surrounding rock are larger than that of specimen 1# without cement, indicating the inclusions in the rock mass are cemented with the surrounding rock and the better the cementation, so the higher the strength, the larger the deformation before failure, the more the energy accumulated, the larger the energy consumption before the sudden failure, the more the formation of rock fragments, and the greater the degree of failure.(2) For specimens 3#, 4#, 5#, and 6# with symmetric preformed cracks and cementation, the failure strength and fractal dimension values of the specimens are lower than that of specimen 2#; the preformed cracks can reduce the failure strength of the inclusions and the generation of fragments.(3) e failure strength and fractal dimension of specimens 4#, 5#, and 6# with the cemented inclusions and preformed cracks are higher than that of specimen 1# Advances in Civil Engineering without the cementation and preformed cracks, indicating no cementation (reducing the cementation between inclusions and the surrounding rock) is more conducive to reducing the failure strength of the rock than by pressure relief of prefabrication cracks, and it is not easy to form more fragments, and the consumption of surface energy is also less.(4) e larger the tip angle of the inclusions, the greater the stress loading required in failure for specimens 3# and 6# and specimens 4# and 5# in the same length of prefabricated crack; in other words, the failure strength of specimen 3# is greater than that of specimen 6#.However, due to the heterogeneity of the rock, many large cracks appear before failure of specimen 3#, so that a large amount of energy is consumed here, and fewer large cracks appear and extend, that is, the fractal dimension of specimen 3# is much lower than that of specimen 6#. e failure strength and fractal dimension of specimen 4# are larger than that of specimen 5#; the reason is that the larger the tip angle of the inclusions, the larger the volume of the inclusions, the more inclined the force surface, the greater the support effect of the inclusions on the surrounding rock mass, and the greater the force required to be overcome in failure.Moreover, it is easy to form more small cracks under high stress intensity to generate more small-sized fragments.( 5) e larger the length of the preformed crack, the more prone the crack propagation to failure, and the smaller the peak stress in failure of the specimen, the smaller the number and scale of fragments for specimens 3# and 4# and specimens 5# and 6# in the same angle of inclusion tip.
In summary, in order to reduce the threat to surrounding people, equipment, and engineering in inclusion rockburst, it is possible to reduce the degree of cementation between the inclusions and the surrounding rock mass, increase the length of the prefabricated cracks, and increase the angle of the inclusions tip, which can reduce the failure degree of the rock mass.

Analysis of Digital Image Correlation Method
e digital image correlation method (DIC) [28] is a fullfield, noncontact optical measurement method that directly extracts deformation information from randomly distributed spots or pseudorandom distribution artificial speckle fields.e basic principle is to match the geometric points on the digitized speckle image in different states of the surface of the object and track the movement of the point to obtain the deformation information of the surface of the object.Given f(x, y) before image deformation and g (x, y) after image deformation, the image subregion corresponding to the image before deformation is identified in the deformed image, and correlation matching operation is performed.
e MV-VD500SM/SC camera is used to record the image of every moment (the speckle movement) in the uniaxial compression process of the inclusion model specimen, and the shooting speed is 1 piece per second.e digital image processing program compiled by MATLAB is used to process the acquired image and convert it into a DIC image with a resolution of 700 pixels × 1200 pixels.e digital marker points are set in the image, and the distance between the adjacent points is 20 pixel × 20 pixel (Figure 15), and the displacement of every point is calculated in the speckle field; the graph is drawn based on the obtained data to measure the horizontal displacement field (u) and the vertical displacement field (v) in the compression process of the inclusion model specimen, and the displacement field before failure immediately is selected, which is shown in Figure 16.
Figure 16 shows the instantaneous u and v displacement fields (magnitude and direction of displacement in the surface of the specimen) of specimens 3#, 4#, 5#, and 6#   inclusion and the lower cut-through part of specimen 4# cause the upper part of the area to have no support, so that the left and lower parts are displaced.It can be seen from Figures 10,16(b), and 16(f ) that the inclusion body and the leftward area of the specimen are 4 to 6 pixels in horizontal displacement (u), and the vertical displacement (v) is downward 12 to 14 pixels.ere is no crack propagation failure in the front region of the right inclusion; only the cement in the right inclusion and the specimen is crushed and ejected out, so that the change of the displacement to the right is great in the right inclusion; the displacement field u is 1 to 3 pixels, and the displacement field u is 13 to 15 pixels.
It can be seen from Figures 11,16(c), and 16(g) in specimen 5# that the local area of the specimen in the deformation and failure corresponds to the u and v displacements obtained by the digital image, that is, the specimen is subjected to tensile shear failure in the lower left and upper right.e horizontal displacement of the lower left rock mass of the specimen is 0 to 1 pixels, and the horizontal displacement of the lower right rock block is 1.5 to 3 pixels.e cut-through cracking failure occurs in the left inclusion and the upper and lower local regions, the corresponding displacement field (Figure 16(c)) is the tip and strip region in the lower left, and the displacement changing to the left is 0.5∼1 pixels.Due to the blocking effect of the inclusions on the vertical displacement of the specimen, the v displacement field changes greatly in the vicinity of the inclusion, and the value is 6.5 to 10 pixels.e right inclusion and its lower region are cut-through, the u displacement field is 2 to 3 pixels, and the v displacement field is 6 to 8 pixels.
It can be seen from Figures 12,16(d), and 16(h) in specimen 6# that the left and right inclusions and its upper and lower regions are cut-through cracking, and the change in the displacement field is relatively uniform.
e horizontal change of displacement in the left failure region is 0 to 0.2 pixels, and the v field displacement is 4.5 to 7 pixels.e change of the u displacement field in the right failure region is 0.6 to 1.8 pixels, and the v field displacement is 3.5 to 7 pixels.
e change of the u and v displacements of specimen 5# is greater than that of specimen 6# in deformation and failure; the reason is that the length of preformed crack of specimen 5# is smaller than that of specimen 6#, and the shorter the length of the crack, the more difficult the rock to be destroyed, the higher the peak stress, and the greater the displacement change before failure.e change of the u and v displacements of specimen 4# is greater than that of specimen 5# in deformation and failure, which indicates that the larger the angle of the tip of the inclusion is, the greater the stress is, and the larger the displacement of the specimen surface is.It is consistent with the law of the stress-strain curve in the test.
e variation law of the four specimens in the direction of the v displacement field is consistent.A large change of displacement occurs above the horizontal band between the inclusions, and the changed displacement in vertical direction near the horizontal band of the inclusion and the area below is reduced; the stratification is obvious, indicating that the inclusions suppress the change of the vertical displacement, the middle area in the vertical displacement field is elliptical, and the displacement of the middle region is larger than the change of the v displacement near the left and right inclusions in the vicinity of the horizontal band, which causes the inclusions to withstand a large concentration of stress and becomes prone to damage.

Analysis of Finite Element Simulation
e finite element software ANSYS 17.0 is used to simulate and analyze the six model specimens in inclusion rockburst test.Based on the size of inclusion specimens (Figure 3), the element type of mud shale rock and cement of models are both three-dimensional solid structural unit SOLID185 with 8 nodes, and every node includes three degrees of freedom to translate along the x-y-z direction.
e mesh size of six models is 2 mm, and there is no cement between the inclusions and the surrounding rock in specimen 1#. e cement, the inclusions, and the specimens are bonded by the command of Glue in specimens 2# to 6#.Local fine division (level 3) is in tip of the model inclusions and the prefabricated cracks (Figure 17).e material parameters are shown in Table 1.e displacement of the bottom surface of the model is constrained in u y direction, and the four sides of the bottom of the model are subjected to the displacement constraint in u xyz direction.e displacement constraint of the surface is applied to the top surface of the model in the linear elastic condition; the displacement is the instantaneous displacement loading of the specimen before failure (Table 2).e distributions of strain energy density, von Mises stress, and x-y shear stress obtained after the model calculation are shown in Figure 18.
It can be seen from Figure 18 that the distribution of the stress and strain energy density obtained at the loading under the elastic condition is roughly consistent with the failure characteristics of the specimen and the displacement change obtained by the digital image method.In the finite element analysis of the specimen, the stress concentration region in the distribution of x-y shear stress is obliquely symmetrically distributed at the tip of the left and right inclusions, that is, light blue (upper left tip and lower right tip), yellow (lower left tip and right tip).e stress concentration gradually weakens from the tip end of the inclusions to the top and bottom of the model, which is consistent with the crack propagation failure from the tip of the inclusion in the test.However, the rock in the test is a heterogeneous brittle material, and there are defects inside, which leads to nonuniformity failure in the process of stress concentration, defects or serious defects near the tip of the inclusions, and the crack appear firstly in failure, and the chronological order of failure is right tip down-left tiptensile failure in specimen 1#; the chronological order of failure is right tip down-left tip up-left tip up-right tip in specimen 2#.Similarly, the stress concentration region of the x-y shear stress distribution in finite element simulation of specimens 3# to 6# is also symmetrically distributed at the tip of the left and right inclusions.Due to the prefabricated 16 Advances in Civil Engineering cracks in the lower part of the specimen, the stress concentration appears near the tips of the cracks.In other words, a large stress concentration band appears between the tip of the inclusion body and the end of the fracture, which is consistent with the fracture cut-through failure region in the lower portion of the tip in the test.e stress band in the tip of Inclusion body tip, the end of prefabricated crack, and the strip bottom area of specimens 3# to 6# is more concentrated    Advances in Civil Engineering than those of specimens 1# to 2#, so the former is easy to break in this area.It can be seen that the results of finite element analysis agrees well with the experimental results in failure mode of the inclusion model specimens.Figure 18(a) shows the energy density of every element of the six specimens before failure, and the strain energy accumulated during the load is calculated.Since the area under the stress-strain curve of the specimen is the work done by an external force on a rock mass unit [29,30], the work (energy) of the specimen obtained by the test system before failure is calculated by the data obtained by the test, which is shown in Table 4.
It can be seen from the above table that the strain energy calculated in the finite element under elastic is greater than the energy calculated in the test; the reason is that the rock in the test is a heterogeneous material, and there are natural defects, resulting in a decrease in energy accumulation.e energy difference affected by the length of the crack in the finite element calculation of specimen 3# and specimen 4# is 0.17 J, and the energy obtained in the test differs greatly.Cracks in the static failure of specimen 3# appear and extend continually during the loading before the peak stress, and the fracture consumed a part of the energy.It can be seen from Tables 2 and 4 that the energy change and the stress changes are consistent in specimens 3# and 6# and specimens 4# and 5#, that is, the larger the angle of the inclusions, the greater the failure stress, the larger the interaction between the inclusions and the rock, the more the stress concentration of the inclusions, and the more the energy accumulated; the distribution of the strain energy density is shown in Figure 18(a).e law of energy calculated is consistent with the results of experimental failure, fractal dimension calculated, and digital image method.

Conclusions
e classification of rockburst has always been one of the foundations for the study of rockburst problems.e model of inclusion rockburst is given in this paper.e simulation test on rockburst is carried out on specimens with inclusions by high-speed photography, digital image correlation method, and finite element simulation to study the failure mechanism of inclusion rockburst, and the effect of the angle of the tip of the inclusion and the preformed crack on the degree of rockburst is analyzed.e following conclusions are drawn: (1) e cracks appear from the tip of the inclusion and expand to the up and down direction in inclusion rockburst.e cracks expand in different orders of the left tip, the lower tip, the left tip, and the right tip; the cut-through rock mass would have been sharply ejected out when the accumulated energy is more, and the fractal dimension is larger.e static failure occurs in the low energy accumulated, and the fractal dimension is low.e lumpiness of rockburst fragments is lower in static load failure.
(2) e analysis from the inclusions shows that the greater the compressive strength of the rock mass, the higher the energy accumulated, and the rockburst ejection damage is severe in the existence of cementation.erefore, the degree of cementation between different rocks should be reduced in the prevention and control of inclusion rockburst.
(3) e angle of the tip of the inclusion has a significant effect on the failure of the rock mass.e larger the angle of the inclusions is, the greater the strength is, the more difficult the failure is, and the larger the fractal dimension of the fragments formed after failure is.When the angle of inclusion is small, the displacement in horizontal and vertical is small.e inclusions have a barrier effect on the change of the rock sample in the vertical direction.e larger the angle of the inclusion is, the greater the barrier effect in the vertical direction is, and the more likely the stress concentration is, thus causing a large accumulation of energy.(4) Prefabrication crack near the inclusions helps to reduce the strength of failure and the accumulation of energy.e larger the length of the prefabricated crack is, the smaller the stress loading required in failure is, the smaller the energy accumulated is, and the smaller the fractal dimension of the fragment is, which can reduce the degree of rockburst.
(5) From specimens 3# and 4#, the increase of the number of cracks and the extended length before failure of rock mass can contribute to the transformation of the dynamic failure of rock mass to the static failure, thereby reducing the occurrence of rockburst.

3. 3 . 3 .
Test of Specimen 3#.When the axial stress σ z is 102.79MPa in the loading of the specimen 3#, secondary cracks 1# and 2# appear at the tip of the lower left prefabricated crack, and the reverse crack 1# with a cracking angle of 123 °extends straightway from the tip to the bottom of the specimen; the airfoil crack 2# with a cracking angle of 20 °gradually extends from the tip end downward to the bottom of the specimen at an arc angle (Figure8(a)).When the axial stress σ z is 106.01MPa, crack 3# and crack 4# appear at the tip of the left inclusion, and crack 3# appears from the tip of the inclusion at an angle of 103 °and continues to expand downward at different angles, which connects to crack 2#. e connecting crack at the bottom of the specimen is 4.52 mm from the left bottom end point.Crack 4# with the cracking angle of 30 °connects to the end of lower left crack (Figure8(b)).When σ z � 111.27MPa, the curve shows a drop state with sound, and the crack at the left tip has a crack angle of 20 °, which gradually expands to the bottom of the test piece with a curved line, and the distance from the end of the crack to the end of the left bottom is 14 mm (Figure8(c)).

Figure 17 :
Figure 17: Nodes distribution and mesh generation.(a) Nodes distribution in no precast cracks.(b) Nodes distribution in precast cracks.(c) Mesh generation in no precast cracks.(d) Mesh generation in precast cracks.

Figure 18 :
Figure 18: e distribution of strain energy density and stress.(a) Strain energy density.(b) vons Mises stress.(c) x-y shear stress.

Table 1 :
Mechanical parameters of shale and rockburst proneness.

Table 2 :
Test results of inclusion rockburst

Table 3 :
Fractal dimension of fragments of rock.Specimen number e division equivalent side length/mm Fit the linear equation Coefficient of association R Fractal dimension D

Table 4 :
Energy calculation results before failure.