Physical Modeling of Displacement and Failure Monitoring of Underground Roadway in Horizontal Strata

Construction Engineering Department, North China Institute of Aerospace Engineering, Langfang 065000, China School of Civil and Environmental Engineering, University of Science and Technology Beijing, Beijing 100083, China State Key Laboratory for Geo-Mechanics and Deep Underground Engineering, China University of Mining and Technology, Beijing 100083, China School of Mechanics and Civil Engineering, China University of Mining and Technology, Beijing 100083, China


Introduction
Coal resources accounted for about 70% of China's energy consumption.More than 90% of the total output is obtained by underground mining, but the mining condition is bad, and the coal resources mining depth increases 8-12 meters every year; over the next 20 years, most of the coal mine in China will enter the 1000-1500 meters mining depth.
e environment of deep mining is very complex; this caused the coal mine safety accident rate to be higher and accurately grasped the status and development trend of the surrounding rock through monitoring and early warning, which is very important for the security of coal mine production.
e stability problem of the roadway surrounding rock is studied by many experts: Ma et al. [1], Shen et al. [2], Zhang et al. [3], Fekete et al. [4], Cohe et al. [5], Fraldi and Guarracino [6], Eberhardt [7], Li et al. [8], Lai et al. [9], and Małkowski [10], by adopting the method of physical model for the study of the issue.e influence of the structural plane and the main structural surface has not been considered in these existing experiments, but these are abundant in the engineering rock mass; at the same time, the displacement and strain of the large-scale physical model are measured by the method of arrangement of the strain monitoring point, and it is difficult to realize the large deformation measurement of the physical model.
At present, the main technology for roadway surrounding rock condition monitoring is deformation monitoring, microseismic monitoring, electromagnetic radiation monitoring, bolt stress monitoring, and so on which are studied by Trifu et al. [11], Brady and Leighton [12], Kent et al. [13], Phillips et al. [14], Tan et al. [15], Yang et al. [16], Lian et al. [17], and so on, by adopting the method of field application for the study of the issue, but the accuracy of the monitoring and early warning of roadway surrounding rock stability needs to be improved.
e physical model is constructed by the "physically finite elemental slab assemblage (PFESA)," the structural plane in which the rock stratum is considered.
e numerical 2D digital image correlation (DIC) technology is used to carry out the real-time monitoring of the surface displacement of the model.Aiming at the problems existing in the stability monitoring of the surrounding rock of the roadway, a small constant resistance bolt (SCRB) is designed with the mechanical and geometric similarity ratio of the constant resistance bolt, and the SCRB is built in the physical model for the monitoring of the mechanical information during the experiment.

Physical Model Construction
2.1.Rock Model Material.Da Tai coal mine is the experimental engineering geological background, and the physical and mechanical parameters of the surrounding rock of the tunnel can be obtained from the field and indoor physical mechanics experiment, as shown in Table 1.
Conventionally, the construction of a large-scale geological model is very expensive and difficult to operate.To overcome these problems, a new approach, that is, the PFESA, initially proposed by He Manchao, was used for the construction of large-scale physical models to mimic stratified layers with different inclinations.is PFESA method uses a specimen-sized rectangular plate, that is, the elemental slab, made of artificial materials similar to the real rocks in terms of the physical and mechanical properties.e elemental slabs can be put into mass production at a low cost.Large numbers of these elemental slabs are arranged, making layers at specified orientations in the simulation of these strata.e elemental slabs were made of similar materials, that is, gypsum and water.e ratio of water to gypsum was obtained according to the test of rock specimens.After testing, the elemental slabs can be recycled and repeatedly used.Hence, large-scale physical models can be constructed conveniently at lower costs [18,19].
Samples were taken from the surrounding rock masses in the roadway in the Da Tai coal mine.
e geometric and physical parameters for the elemental slabs are presented in Table 2.
e dimensions of the elemental slabs are 400 mm × 400 mm × 10/20/30 mm.A large number of these elemental slabs were assembled in layers with different mechanical properties, simulating the rock strata of the coal seam, diabase, and green schist of the real rock masses.e rock layers were arranged at a zero dip angle with respect to the horizontal structure in order to mimic the horizontally distributed strata in the real rock in the test conducted in our investigation.

Rock Strata.
Figure 1 shows the principle for producing rock strata.A large number of the elementary slabs (shown in Figure 1(b)) with the same model rock property were used to ensemble a rock layer of the same property, and a certain number of the rock layers with the same rock property were used to construct a rock stratum.As shown schematically in Figure 1(a), the green schist stratum and coal seam stratum are assembled by placing the elementary slabs in layers with perfectly mated interfaces.A weak surface is formed between the two different model rock strata.
is weak surface is homogeneous along its striking at a macroscopic scale, while it has many minor asperities at a "microscopic scale" constituted by many small gaps between the parallel-placed elemental slabs in the rock layer.e gaps and asperities formed by the layered slabs also exist within the rock strata.
ese minor flaws can be viewed as the origin of the heterogeneity of the rock under external loading.ree classes of the model rock strata were constructed in our test, that is, diabase stratum, green schist stratum, and coal stratum which create two types of interfaces or weak faces, that is, green schist-coal interface and diabase-coal interface.

Geological Model.
e large-scale physical model was built with a total of six strata including one diabase, three green schists, and two coal seams, and all the strata are inclined at an angle of 0 °with respect to the horizontal structure, to simulate near horizontal geological structures.
e constructed geological model, as shown schematically in   3.

Digital Image Correlation.
e digital image correlation (DIC) is a no-contact displacement-strain measurement technique, which is also called the digital speckle correlation method (DSCM), and the image measurement is realized based on the computer vision technology.rough noncontact means, the full-eld displacement and strain measurement can be achieved.
e DIC technique has obvious advantages compared with the contact-type displacement and strain measurement.
rough real-time monitoring of surface displacement, the displacement and strain information of the measured object surface can be obtained, collected, and extracted.
e MTI-2D digital image correlation measurement system is shown in Figure 3. e system mainly includes the image acquisition system, lighting system, and digital image   processing system.In the experiment, the common white light can be used as the light source, and the illumination system can be added to ensure that the uniform light eld is formed on the surface of the test piece.In order to obtain the accurate deformation image, the CCD camera is placed on the precision three-dimensional support and is made parallel to the surface of the specimen.Image acquired is directly deposited onto the image card and imported into the computer.Image saving format can be of BMP or ti format.
In the image analysis step, the image extracted will be saved for further analysis.

SCRB Monitoring.
With the increase of mining depth, the failure of roadways often has the characteristics of large nonlinear deformation.Under the mining conditions such as mining disturbance and high pressure, the traditional small deformation rock bolts can only bear the deformation of the roadway surrounding rock within 200 mm and cannot adapt to the large deformation and failure characteristics of the surrounding rock of the deep roadway; therefore, in roadway support projects, it often appears that the bolt cannot adapt to the large deformation and tensile failure of the surrounding rock, resulting in the tunnel collapse and roof fall accident.
In view of the disadvantage of the small deformation of the bolt in the deep roadway support, domestic and foreign scholars have conducted studies on the tensile strength of bolts and achieved fruitful results.e energy-absorbing bolt can provide constant support resistance in the roadway creep stage and provide a large amount of extension.e mechanical characteristics, structural characteristics, and working mechanism of several typical bolts are analyzed by Professor He MC, and a new type of an energy-absorbing anchor rod is developed which is called the constant resistance and large deformation bolt (CRLB).
rough a large number of laboratory experiments and the engineering eld experimental studies, the CRLB can also provide a constant support resistance in the process of deformation.e structure of the CRLB includes the nut, the ball pad, the tray, the constant resistance device, and the rod body.
e constant resistance device can provide deformation and constant resistance; the tray and the nut are arranged at the tail part of the constant resistance device, and the constant resistance device is connected with the nut through the screw thread.CRLB is shown in Figure 4.
e constant resistance device is the main component of the CRLB to ensure the constant resistance of the device, and its design strength is 80%-90% of the anchor rod body.When the surrounding rock deformation is large, CRLB resistance can be kept constant and the deformation energy of the surrounding rock can be absorbed by CRLD; hence, CRLD can still play a good supporting e ect when the surrounding rock is under large deformation state.
A small constant resistance large deformation bolt (SCRB) is developed, which is suitable for indoor physical model experiments.e supporting structure used in the experiment only considers the similarity of mechanics and geometry.According to the similarity theory of the physical model, the geometric similarity ratio of the bolt is C l 15, which is used to calculate the length and diameter of the anchor rod in the physical model.e mechanical similarity  4. e SCRB is composed of rod and constant resistance device; the constant resistance device meet the mechanical and geometric similar ratio, and it can provide 40 mm deformation.
e constant resistance device is assembled by tray and round tube; the outer diameter of the round tube is 60 mm, and the thickness is 2.5 mm. e deformation length of the SCRB can be adjusted by setting the length of the constant resistance device; in this paper, the constant resistance value of the constant resistance is 90% of the strength of the rod body.e diagram of SCRB is shown in Figure 5.
With reference to the situation of roadway support in Da Tai coal mine, the arrangement of SCRB and SB monitoring devices is shown in Figure 6. e constant resistance device and anchor rod body and pull pressure sensor are assembled into the SCRB monitoring device, and the layout and number of the bolts are shown in Table 5.

Loading Path and Experimental Process.
e loading path is derived from the measurement and calculation results of in situ stress, and the load and stress similarity ratio is also set by in situ detection and measurement.e vertical and horizontal stress of the roadway can be obtained.e horizontal and vertical direction of the model tunnel is loaded according to the stress similarity ratio.e di erent levels of model boundary load calculations are as follows: where σ x is the horizontal stress boundary, σ y is the vertical stress boundary, σ v is the weight of the di erent depth of stress, σ h is the horizontal stress, and C σ is the stress similarity ratio.e loading process of the model experiment is designed for four stages: A-D, by applying σ x and σ y realization.e loading process and the loading path are shown in Figure 7.
e horizontal stress and vertical stress increase at the same time during the A stage, and the horizontal stress remains unchanged; the vertical stress increases at the B stage, and the vertical stress remains unchanged; the horizontal stress increases at the C stage; the horizontal stress and vertical stress increase at the same time at the D stage.
At stage A, the horizontal and vertical stresses increase to 0.2 MPa.After 4 times of loading, 0.8 MPa stress is achieved, and each loading time is for 30 minutes; at this stage, the gaps between PFESAs are closed.At stage B, the horizontal stress (0.8 MPa) remains unchanged.After 10 times of loading, the vertical stress of 2.8 MPa is achieved, and each loading time is 30 minutes.e loading increases 0.2 MPa each time.en, at the C loading stage, the vertical stress remains unchanged, the horizontal stress is increased from

Vertical Displacement Field.
rough the digital image acquisition and processing in the model experiment, the evolution characteristics of the whole vertical displacement eld of the physical model can be obtained as shown in Figure 8.
During stage A, the downward vertical displacement value has been in a state of growth, and in the A4 stage, the maximum vertical downward displacement value is 27.65 mm.With the increase of the vertical load, the vertical displacement value increases; when the value of the vertical load reaches 1.4 MPa, the maximum vertical displacement value is about 31.50 mm.At the C1 stage, the vertical load (2.8 MPa) remains unchanged and the horizontal load increases; there is a rapid increase of the vertical displacement, and the maximum vertical displacement is about 32 mm to 40 mm, which remains basically in the stable state in the C1∼C8 stage vertical displacement.
During the C9 stage, the horizontal load is 2.8 MPa, the vertical load is 4.4 MPa, the model of this stage begins to break down, the bottom drum and roof fall o , the bottom heave value is 111.8 mm, and the roof falls down 75.89 mm.

Horizontal Displacement Field.
rough the DIC in the model experiment, the evolution characteristics of the whole horizontal displacement eld of the physical model can be obtained as shown in Figure 9.
From the graph, we can see that, during the A1-A2 loading stage, the horizontal displacement rapidly generates after loading and the value is about 4 mm-12 mm; the displacement region is concentrated on two sides of the roadway model at the A3-B10 stage, mainly when the vertical load is applied, and the horizontal displacement is smaller in magnitude.
When the vertical load of 2.8 MPa is applied, at the C1 stage, the horizontal load began to increase, and the model produces about 10 mm-22 mm left displacement, until the C8 stage, the distance of the roadway on the right side is about 2.69 mm and the left side is 8.72 mm.
At the C9 loading stage, the horizontal load value is 4.4 MPa, the vertical load value is 2.8 MPa, the nal distance  Advances in Civil Engineering of the roadway on the left side is about 36.37 mm, and on the right side, it moves 37.84 mm away.

Stress Monitoring of the Rock Bolt.
In the course of the experiment, curves of bolt axial force with time are obtained as shown in Figure 10.e roof and two sides monitoring bolts are arranged as shown in Figure 6, in which R1 and R3 are SCRB monitoring devices and R2 and R4 are SB monitoring devices.From the curves, it can be obtained that, during the A1-A4 loading stages, the bolt axial force monitoring curve of the roof continued to rise.During B1-B10 loading stages, the monitoring curve slows down, and from the C loading stage, the monitoring curve rises rapidly.R2 and R4 monitoring curves in the axial force peak decreased dramatically after that in the roadway deformation and failure process because the force has exceeded its capacity, the rod body occurrence pulled, and failure occurred.R1 and R3 monitoring curves, respectively, maintain 220 N and 230 N stress, until the complete model failure.e SCRB can withstand the constant resistance monitoring of the surrounding rock and will not break until the roadway is completely destructed, and the SCRB monitoring device can monitor the whole failure process of the roadway.
e monitoring curve of bolt axial force on two sides of the roadway is shown in Figure 11, in which S1 and S3 are SCRB monitoring curves and S2 and S4 are SB monitoring curves.In the deformation and failure mode, the roadway is dominated by the roof oor heave and roof falling; therefore, the two-side roadway has less stress, resulting in the small bolt force value, and the bolt force did not exceed the load limit.
rough analysis of the monitoring curve, we conclude that the SCRB monitoring device can realize the whole  8 Advances in Civil Engineering process monitoring of the roadway deformation and failure.
In the experiment, deformation and failure of the roadway mainly occurred on the roof and oor; two sides destroy model is overall launch, so the force monitoring e ect of the two sides is poor.

Stress and Deformation
Monitoring. Figure 12 is the axial force monitoring curve of R1 and R3 bolts and is the roof horizontal displacement and vertical displacement monitoring curve at the same time period.During the loading phase B, the roof displacement remained unchanged, and the bolt stress kept slowly rising.At C5 loading stage, the roadway cracks of the model began to develop, the bolt axial force also began to accelerate, and the roof displacement is still unchanged.From the C9 loading stage, the roadway began to destroy; at this time, the horizontal displacement and vertical displacement of the moment changed.By the analysis of the acceleration of the axial force of the bolt, it is found that it can be used as the precursor information of the failure of the model roadway.
Figures 13 and 14 are the S1 and S3 axial force monitoring curves and the displacement monitoring curves at the corresponding position.From the gures, we can see that the force value before and after the deformation and failure of the roadway has not reached its extreme value, which is caused by the failure of the roadway; after the model failure, the force curves have a rising trend, but the force monitoring e ect is not obvious compared to the roof bolt.
R1 and R3 have the best warning e ect, the bolt axial force monitoring curves have experienced the slow rising stage and the accelerated phase, and the acceleration stage   Advances in Civil Engineering can be used as the information of the early warning stage for the deformation and failure of the roadway.

Numerical Simulation Verification
To validate the model experiment, the discrete element method is used to simulate the same conditions but di erent rock inclination angles.e numerical model of the dip rock stratum at 15, 30, 45, 60, 75, and 90 degrees is established, and the force of the bolt and the displacement of the surrounding rock are monitored.e SCRB monitoring devices are numbered R1, R2, S1, and S2, and the numerical model is shown in Figure 15.
e displacement eld of di erent rock inclination angles is simulated as shown in Figure 16; therefore, the deformation and failure of the roadway surrounding rock are signi cantly a ected by rock inclination.e failure of the roadway surrounding rock is caused by large asymmetric deformation, and the instability model is dominated by shear slip in the key part.
e monitoring e ect of the monitoring bolt is a ected by the distribution of the roadway failure area, and the monitoring e ect in key parts of the surrounding rock is better.When arranging the monitoring bolt, the failure eld distribution of the surrounding rock should be taken into full consideration.
rough the physical model experiment and numerical simulation of di erent rock inclination angles, the monitoring value of the SCRB monitoring device is 0.7-0.8times of its constant resistance.e characteristic changes of the mechanical monitoring curve of the SCRB monitoring device is shown in Figure 17.
e monitoring curve has a signi cant ascending stage; when it reaches the warning value, the corresponding time is the warning time; when the monitoring curve exceeds the warning value to the maximum value, the roadway surrounding rock failure occurs.

Conclusion
Aiming at the problem of monitoring and early warning of the stability of the surrounding rock of the layered roadway, the monitoring technology of surrounding rock stability of the roadway based on the CRLB is put forward.According  Advances in Civil Engineering to the similarity theory, the SCRB is developed which can be applied to the indoor physical model experiment, and SB and SCRB monitoring devices are built into the physical model body.rough the physical model experiment and numerical simulation, the following conclusions are obtained: (1) e SCRB can withstand large deformation of roadway failure, and the whole process of monitoring roadway deformation destruction can be achieved, while the SB found it difficult to realize the roadway deformation and failure, which is the drawback of the whole process of monitoring and early warning.(2) e force monitoring curves of the SCRB are in an accelerating stage before the large deformation and failure of the roadway, and then the curve rises to a certain value after the roadway deformation and failure; so the accelerated phase of the monitoring curve can be used as the early warning information of the roadway surrounding rock.(3) Numerical simulations verify the correctness of the physical model experimental results, and the monitoring value of the SCRB monitoring device is 0.7-0.8times its constant resistance.

Figure 2 (
Figure 2(a), is 1600 mm high, 1600 mm wide, and 40 mm thick.According to the size of the surrounding rock loose circle and the model loading system, the geometric similarity ratio is C l 15 and the size of the roadway of the model has a length 300 mm and the height 200 mm.According to the maximum load set degree of the model and the mechanical parameters of the engineering rock mass, the similarity coe cient of stress is C σ 9, the similarity coe cient of density is C ρ 9/15 0.6, and the physical model is shown in Figure 2. e alternating strata are indexed with 1-6 from top to bottom with the geological pro le, and geometric parameters of the physical model are shown in Table3.

Figure 2 :
Figure 2: Physical model design: (a) schematic of the model; (b) photograph of the model; (c) rock reality of the deep underground roadway simulated by the gypsum board.

Figure 1 :
Figure 1: Weak bedding surface of the rock strata simulated by the similar materials; (a) schematic of the rock strata and their interface; (b) photograph of the gypsum boards.

Figure 4 :
Figure 4: Sketch map of the constant resistance and large deformation bolt.

Figure 5 :
Figure 5: Diagram of the small-sized constant resistance and large deformation bolt.

Figure 6 :
Figure 6: Schematic diagram of the physical model of the monitoring bolt device arrangement.

Figure 7 :
Figure 7: Loading scheme of the physical model.

6
Advances in Civil Engineering e vertical load during the B3∼B10 stages is increased from 1.4 MPa to 2.8 MPa, and the vertical displacement is basically unchanged and is between −1.85 mm and 32.74 mm.

Table 1 :
Mechanical parameters of the engineering rock mass.

Table 2 :
Mechanical parameters of the physically finite elemental slab assemblage.

Table 3 :
Geological section, material, and structural parameters of the physical model from top to bottom.

Table 4 :
Prototype and model physical and mechanical parameters of the bolt.

Table 5 :
Layout and number of monitoring bolt devices.
Figure 10: Bolt axial force-time curve of the roadway roof.