Fracturing Behaviors and Mechanism of Serial Coal Pillar Specimens with Different Strength

The fracturing behaviors of serial coal pillars is significant for understanding their failure mechanism. To reveal this, the bearing stress, acoustic emission, electrical resistivity, local strain, force chain distribution, and cracks evolution of serial coal pillars under uniaxial compression were evaluated by experiment and numerical simulation. The results show that four bearing stages are observed during the fracturing process (i.e., nonlinear growth, linear growth, yielding growth, and weakening stages). The acoustic emission features, electrical resistivity responses, strain develops, force chain distributions, cracks evolutions, and local displacement are highly consistent to illustrate the fracturing behaviors. System fracturing of serial coal pillar specimens is appeared along with the collapse of lower uniaxial compressive strength coal pillar specimen. The limit bearing capacity of serial coal pillar specimens is almost equal to the strength of lower uniaxial compressive strength coal pillar specimen. The unbalanced deformation characteristics of serial coal pillar specimens are presented due to the strength differences. The evolution of the key deformation element is the rooted reason for the overall fracturing mechanism of serial coal pillar specimens. For serial coal pillar specimens with different strengths, the critical condition of system fracturing is that the sum of secant modulus of upper and bottom coal pillars is zero, which is expected to predict the system fracturing of serial pillars in the underground coal mining.


Introduction
In recent years, an increasing number of ultra-closed mining activities have been conducted in China [1,2], America [3,4], Australia [5,6], India [7], and other famous coalmining countries [8,9]. According to different mining methods at the upper and bottom levels, various multi-seam mining configurations are defined. Multi-seam pillar mining is one of the main configurations in practice [10]. Generally, several residual coal pillars are arranged and left permanently in ultra-closed mining to support the weight of overburden strata, which are overlapped or staggered at the different levels [11,12]. Overlapped residual coal pillars can be considered as a whole system [13], which can also be termed as serial coal pillars.
Long-term stability of serial coal pillars is one of the significant issues to be focused on in the ultra-closed goaf [14]. With the coupled effects of environmental weathering, results. Berisavljević et al. [31] inspected the influence of lithological heterogeneity on serial sandstone-siltstone specimens' strength. Furthermore, Lu et al. [32], Huang et al. [33], and Liu et al. [34] investigated the mechanical behavior and failure mode of serial coal-rock specimens by a group of uniaxial compression tests.
These laboratory tests provided a deeper understanding about the fracturing of serial coal/rock masses. However, all of the experimental tests were just concentrated on the serial rock-rock specimens or serial coal-rock specimens. Due to the differences of combining size and contacting form, the above laboratory results could not be directly applied in the serial coal-coal specimen testing. The problem to be solved was also different. The serial rock-rock specimens were prepared to reveal the mechanism of rock burst. The serial coal-rock specimens were proposed to investigate the interaction of roof, coal seam and floor at the single mining level. However, the serial coal-coal specimens were designed to explore the interaction of upper and bottom coal pillars at the ultra-closed goafs. Therefore, further laboratory observations about serial coal pillar specimens with different mechanical properties are necessary to make a throughout insight about the fracturing behaviors and mechanism. Meanwhile, it should also be verified by corresponding numerical simulation and explained by the related theorical method.
In this study, the laboratory uniaxial compression tests and numerical simulations are conducted for serial coal pillar specimens with different strength. During the loading process, the fracturing behaviors are investigated comprehensively based on bearing stress, acoustic emission (AE) features, electrical resistivity responses, strain develops, force chain distributions, cracks evolutions, and local displacement. A schematic model of serial coal pillar specimens is also established. Then, the key deformation element is defined and the fracturing evolution is analyzed to reveal the unbalanced deformation characteristics. Finally, the progressive fracturing mechanism is revealed theoretically. Figure 1 shows a flowchart summarizing the main content of this research. uniaxial compressive experiments to determine the influence of weak rock thickness and moisture content on the strength of serial sandstone-shale specimens. Tziallas [30] proposed a methodology to predict the strength of serial sandstone and siltstone specimens based on the laboratory results. Berisavljević et al. [31] inspected the influence of lithological heterogeneity on serial sandstone-siltstone specimens' strength. Furthermore, Lu et al. [32], Huang et al. [33], and Liu et al. [34] investigated the mechanical behavior and failure mode of serial coal-rock specimens by a group of uniaxial compression tests. These laboratory tests provided a deeper understanding about the fracturing of serial coal/rock masses. However, all of the experimental tests were just concentrated on the serial rock-rock specimens or serial coal-rock specimens. Due to the differences of combining size and contacting form, the above laboratory results could not be directly applied in the serial coal-coal specimen testing. The problem to be solved was also different. The serial rock-rock specimens were prepared to reveal the mechanism of rock burst. The serial coal-rock specimens were proposed to investigate the interaction of roof, coal seam and floor at the single mining level. However, the serial coal-coal specimens were designed to explore the interaction of upper and bo om coal pillars at the ultra-closed goafs. Therefore, further laboratory observations about serial coal pillar specimens with different mechanical properties are necessary to make a throughout insight about the fracturing behaviors and mechanism. Meanwhile, it should also be verified by corresponding numerical simulation and explained by the related theorical method.
In this study, the laboratory uniaxial compression tests and numerical simulations are conducted for serial coal pillar specimens with different strength. During the loading process, the fracturing behaviors are investigated comprehensively based on bearing stress, acoustic emission (AE) features, electrical resistivity responses, strain develops, force chain distributions, cracks evolutions, and local displacement. A schematic model of serial coal pillar specimens is also established. Then, the key deformation element is defined and the fracturing evolution is analyzed to reveal the unbalanced deformation characteristics. Finally, the progressive fracturing mechanism is revealed theoretically. Figure  1 shows a flowchart summarizing the main content of this research.

Experimental Test Setup
In this section, the experimental specimens, experimental testing system, and experimental scheme will be introduced in detail.

Experimental Specimens
The serial coal pillars were generally distributed at different goaf levels. They played a combining bearing role with the surrounding rock strata [6,10,35]. That is, the roof strata, upper coal pillar, parting strata, bottom coal pillar, and floor strata should be considered as a whole system. Previous researches had confirmed that the coal pillars were easier to be damaged than the rock strata in the combining bearing system. That is, the stability of serial coal pillars should be paid much more attention [36]. Once the initial fracturing of serial coal pillars was controlled with proper measurements, the stability of the whole bearing system could then be maintained. Therefore, two overlapped cylindrical specimens were designed to model the serial coal pillars without consideration of parting rock in this study.
Obtained from the Datong coal basins of Shanxi in China, the experimental coal pillar specimens were prepared from natural coal mass. According to the methods suggested by the International Society for Rock Mechanics [37], the coal pillar specimens were cored, shaped and polished carefully. The diameter of experimental coal specimens was 50 mm, while the height was 100 mm. Afterwards, the polished coal pillar specimens were kept in a constant temperature and humidity incubator for 24 h before the experiment. The temperature and humidity were set to 20 • and 70%, respectively. Figure 2 shows the photograph of experimental coal pillar specimens.

Experimental Test Setup
In this section, the experimental specimens, experimental testing system, and experimental scheme will be introduced in detail.

Experimental Specimens
The serial coal pillars were generally distributed at different goaf levels. They played a combining bearing role with the surrounding rock strata [6,10,35]. That is, the roof strata, upper coal pillar, parting strata, bo om coal pillar, and floor strata should be considered as a whole system. Previous researches had confirmed that the coal pillars were easier to be damaged than the rock strata in the combining bearing system. That is, the stability of serial coal pillars should be paid much more a ention [36]. Once the initial fracturing of serial coal pillars was controlled with proper measurements, the stability of the whole bearing system could then be maintained. Therefore, two overlapped cylindrical specimens were designed to model the serial coal pillars without consideration of parting rock in this study.
Obtained from the Datong coal basins of Shanxi in China, the experimental coal pillar specimens were prepared from natural coal mass. According to the methods suggested by the International Society for Rock Mechanics [37], the coal pillar specimens were cored, shaped and polished carefully. The diameter of experimental coal specimens was 50 mm, while the height was 100 mm. Afterwards, the polished coal pillar specimens were kept in a constant temperature and humidity incubator for 24 h before the experiment. The temperature and humidity were set to 20° and 70%, respectively. Figure 2 shows the photograph of experimental coal pillar specimens. Two types of coal pillar specimens were prepared in this experiment, which were collected from different mining level in Datong coal field. Table 1 lists the basic physical and mechanical parameters of testing coal pillar specimens, which are determined averagely from five series of tests. Two types of coal pillar specimens were prepared in this experiment, which were collected from different mining level in Datong coal field. Table 1 lists the basic physical and mechanical parameters of testing coal pillar specimens, which are determined averagely from five series of tests.

Experimental Testing System
The experimental testing system consisted of uniaxial loading, AE monitoring, electrical resistance testing, and strain testing systems, as shown in Figure 3. During the fracturing process of specimens, the mechanical behaviors, AE features, electrical resistivity responses and local strain develops could be monitored simultaneously by this testing system.

Experimental Testing System
The experimental testing system consisted of uniaxial loading, AE monitoring, electrical resistance testing, and strain testing systems, as shown in Figure 3. During the fracturing process of specimens, the mechanical behaviors, AE features, electrical resistivity responses and local strain develops could be monitored simultaneously by this testing system.

Uniaxial Loading System
In the present experiments, the confining stress of serial coal pillar specimens was regarded to be very small since the adjacent coal resources were excavated. Furthermore, the axial stress of serial coal specimen with a deeper buried condition was relatively large. Therefore, the axial stress was viewed as the main loading factor associated to the fracturing of serial coal pillar specimens. The axial loading condition was thus selected in the following laboratory experiments.
A servo-controlled testing machine was applied to achieve the uniaxial loading condition. As shown in Figures 3 and 4, two rigid steel loading plates, with dimensions of a × a × h = 400 mm × 400 mm × 80 mm, were placed between the loading frame and the testing coal specimens.

Uniaxial Loading System
In the present experiments, the confining stress of serial coal pillar specimens was regarded to be very small since the adjacent coal resources were excavated. Furthermore, the axial stress of serial coal specimen with a deeper buried condition was relatively large. Therefore, the axial stress was viewed as the main loading factor associated to the fracturing of serial coal pillar specimens. The axial loading condition was thus selected in the following laboratory experiments.
A servo-controlled testing machine was applied to achieve the uniaxial loading condition. As shown in Figures 3 and 4, two rigid steel loading plates, with dimensions of a × a × h = 400 mm × 400 mm × 80 mm, were placed between the loading frame and the testing coal specimens.

Experimental Testing System
The experimental testing system consisted of uniaxial loading, AE monitoring, electrical resistance testing, and strain testing systems, as shown in Figure 3. During the fracturing process of specimens, the mechanical behaviors, AE features, electrical resistivity responses and local strain develops could be monitored simultaneously by this testing system.

Uniaxial Loading System
In the present experiments, the confining stress of serial coal pillar specimens was regarded to be very small since the adjacent coal resources were excavated. Furthermore, the axial stress of serial coal specimen with a deeper buried condition was relatively large. Therefore, the axial stress was viewed as the main loading factor associated to the fracturing of serial coal pillar specimens. The axial loading condition was thus selected in the following laboratory experiments.
A servo-controlled testing machine was applied to achieve the uniaxial loading condition. As shown in Figures 3 and 4, two rigid steel loading plates, with dimensions of a × a × h = 400 mm × 400 mm × 80 mm, were placed between the loading frame and the testing coal specimens.  In this study, the displacement control model was selected according to similar papers related to the fracturing of coal/rock pillars [14,21,23]. Each displacement growth meant the mining disturbance to the serial coal pillar specimens [14,38]. Along with the growth of axial displacement, the stress state at the defect areas in the serial coal pillar specimens was satisfied to the strength criterion. Then, the damage would be induced at such areas to reach a new equilibrium state. At this situation, the stress redistribution was induced, which would give birth to further ruptures of adjacent areas in the serial coal pillar specimens [14]. That is, the stress redistribution was constantly repeated during the loading process. It would finally lead to the macroscopic fracturing of serial coal pillar specimens. Therefore, the displacement rate was 0.12 mm/min in the control model.
The upper loading plate was stationary during the loading process, while the bottom loading plate was movable. The displacement boundary conditions were shown as follows: (1) the axial displacement was zero at the upper boundary of serial coal pillar specimens; (2) for the bottom boundary, the axial displacement was equal to the rising distance of bottom loading plate. It was the product of displacement-control rate (0.12 mm/min) and the loading time. Besides, the stress boundary conditions were also presented: the horizontal stress was zero at the upper and bottom boundaries of serial coal pillar specimens since the end frictional effect in this testing was neglected [39].

AE Monitoring System
Real-time AE features can make a clear explanation about the internal damage of testing materials [40,41]. The AE monitoring system was utilized to observe the AE characteristics of serial coal pillar specimens. As shown in Figure 3, it consisted of AE sensors, AE preamplifiers, an AE data acquisition instrument and data processing system.
Four sensors (RS2A-1, RS2A-2, RS2A-5, RS2A-6) were employed to record the AE parameters. Two sensors (RS2A-2, RS2A-6) were attached on the surface of upper coal pillar specimen. The other ones (RS2A-1, RS2A-5) were located on the surface of bottom coal pillar specimen. The arrangement position of AE sensors is shown symmetrically in Figure 4. To improve the signal transmission, silicon grease was applied as coupling agent between AE sensors and coal pillar specimens. The data acquisition frequency and amplitude threshold were set at values of 2.5 MHz and 40 dB, respectively.

Electrical Resistance Testing System
Electrical resistivity is an inherent physical property that defines the resistance characteristics of testing specimen [42,43]. The responses of electrical resistivity can reveal the fracturing behaviors of coal/rock at different loading stages [44] which are helpful for predicting the dynamic collapse induced by pillars' instability. In the current experiments, serial coal pillar specimens as a whole were regarded as an electrical resistance element.
Generally, the electrical resistivity of a single coal pillar specimen is determined [45] by Equation (1).
where ρ is the electrical resistivity of a single coal pillar specimen (Ω·m); R is the electrical resistance of a single coal pillar specimen (Ω); S is the effective cross-sectional area of a single coal pillar specimen (m 2 ); H is the height of a single coal pillar specimen (m); r is the radius of a single coal pillar specimen (m) For the serial coal pillar specimens, the electrical resistivity can be calculated by the Equation (2).
where ρ p is the electrical resistivity of serial coal pillar specimens (Ω·m); R p is the electrical resistance of serial coal pillar specimens (Ω); S p is the effective cross-sectional area of serial coal pillar specimens (m 2 ); H p is the total height of serial coal pillar specimens (m); R u and R b are the electrical resistance of upper and bottom coal pillar specimen (Ω), respectively. A two-electrode insulation resistance testing system was adopted for measuring the electrical resistance of serial coal pillar specimens. It consisted of resistance testing holders, copper electrodes, insulation boards, a data acquisition instrument and an analysis software. The maximum testing resistance was 1012 Ω, and the measuring deviation was within 1%.
The data acquisition frequency was 5 Hz. As shown in Figures 3 and 4, copper electrodes, with dimensions of a 2 × a 2 × h 2 = 200 mm × 200 mm × 1 mm, were fixed symmetrically at both sides of the serial coal pillar specimens. They were sandwiched between insulation boards (a 1 × a 1 × h 1 = 300 mm × 200 mm × 1 mm) and testing coal pillar specimens. Resistance testing holders were connected with upper and bottom copper electrodes, which were also linked with the data acquisition instrument.

Strain Testing System
A strain testing system was utilized to record the axial and lateral strains of serial coal pillar specimens. In this experiment, two pairs of strain gauge rosettes were attached by silastic glue at the central surface of upper and bottom coal pillar specimens, respectively (see Figure 4). Each strain gauge rosette was composed of two strain gauges with a sensitivity of 10 −6 ε. They were arranged along the horizontal and longitudinal orientations of the specimens. During the fracturing process, the strain gauge rosettes were connected with the data acquisition system by the lead wires. The local strain of serial coal pillar specimens could be recorded continuously.

Experimental Scheme
Three groups of overlapped coal pillar specimens were prepared and tested to obtain the fracturing information under the uniaxial compression condition. Table 2   It should be noted that the upper and bottom coal specimen didn't fail simultaneously due to the different mechanical parameters. That is, two possible conditions may be observed after the local collapsing of upper or bottom coal specimen. One condition was that the upper coal specimen may lose the support of bottom collapsed specimen and become unstable. The other one was that a relatively long distance was formed between the upper loading plate and the bottom coal specimen after the collapse of upper coal specimen. At both conditions, the loading stress could not directly transfer to the non-collapsing coal pillar specimen. The deviation and dropping of AE sensors, copper electrodes and strain gauges, attached on the surface of collapsed coal specimen, may be also observed. Therefore, the data monitoring had to be stopped after the local collapse of upper or bottom coal specimen.

Experimental Results
In this section, fracturing behaviors of serial coal pillar specimens are presented, which are both described by the mechanical behaviors, AE features, electrical resistivity responses and local strain develops. For different groups of serial coal pillar specimens, the overall evolution of fracturing behaviors was very similar. Therefore, only a group of testing specimens was selected as the representative one. The peak load of the selected serial specimens was very close to the average UCS.

Stress-AE Energy-Electrical Resistivity Response
The curves of axial stress, AE energy, accumulated AE energy and electrical resistivity for the representative specimens are shown in Figure 5. The black curve in Figure 5 presents the response of axial stress. The red curve in Figure 5 shows the evolution of AE energy. The accumulated AE energy is presented as the blue curve in Figure 5. The green curve in Figure 5 exhibits the evolution of electrical resistivity. It should be noted that the electrical resistivity of serial coal pillar specimens in different groups was calculated based on the real-time electrical resistance monitored during the fracturing process.
ity for the representative specimens are shown in Figure 5. The black curve in Figure 5 presents the response of axial stress. The red curve in Figure 5 shows the evolution of AE energy. The accumulated AE energy is presented as the blue curve in Figure 5. The green curve in Figure 5 exhibits the evolution of electrical resistivity. It should be noted that the electrical resistivity of serial coal pillar specimens in different groups was calculated based on the real-time electrical resistance monitored during the fracturing process.
According to the bearing behaviors, it could be roughly divided into four stages. Stage I was the nonlinear growth stage of bearing capacity. Stage II was the linear growth stage of bearing capacity. Stage III was the yielding growth stage of bearing capacity. Stage IV was the weakening stage of bearing capacity. A thorough description of stress evolution is presented as follows (see Figure 5). (1) At stage I, axial stress of serial coal pillar specimens was growing with an upconcaved form. Original pores and cracks inside the testing coal pillar specimens were gradually enclosed. That is, the serial coal pillar specimens mainly experienced the compaction deformation at stage I. AE energy was primarily released due to the gradual closure of initial pores and cracks at stage I. A slight fluctuation tendency was detected. It maintained a relatively small value. As a result, a weak growth of accumulated AE energy was observed. The initial electrical resistivity of serial coal pillar specimens was 2.48 According to the bearing behaviors, it could be roughly divided into four stages. Stage I was the nonlinear growth stage of bearing capacity. Stage II was the linear growth stage of bearing capacity. Stage III was the yielding growth stage of bearing capacity. Stage IV was the weakening stage of bearing capacity. A thorough description of stress evolution is presented as follows (see Figure 5).
(1) At stage I, axial stress of serial coal pillar specimens was growing with an upconcaved form. Original pores and cracks inside the testing coal pillar specimens were gradually enclosed. That is, the serial coal pillar specimens mainly experienced the compaction deformation at stage I. AE energy was primarily released due to the gradual closure of initial pores and cracks at stage I. A slight fluctuation tendency was detected. It maintained a relatively small value. As a result, a weak growth of accumulated AE energy was observed. The initial electrical resistivity of serial coal pillar specimens was 2.48 MΩ·m before the axial loading, respectively. With the evolution of compaction deformation, original and newborn electrical channels were combined to improve the conductive properties [45]. The electrical resistivity of serial coal pillar specimens was cutting down gradually as a whole. The fluctuation characteristics of electrical resistivity were also revealed at stage I. It could be explained by the heterogeneity of the coal pillar specimen and the uncertainty of initial damage's expansion direction.
(2) Stage II had the longest evolution time during the entire loading fracturing process. The duration of stage II was approximately 1070s. At this stage, the linear growth of bearing capacity was monitored. That is, elastic deformation was dominant for the serial coal pillar specimens. With the linear growth of bearing capacity, AE behavior of serial coal pillar specimens was activated little by little. The accumulated AE energy presented a gentle smooth increasing tendency at stage II. Moreover, a continuous descending tendency of electrical resistivity was detected at stage II. Due to local damage of conductive channels in the testing specimens, the electrical resistivity was decreasing with a weakening rate. The fluctuation at stage II almost disappeared.
(3) The duration of stage III is shortest for the testing serial specimens. The carrying capacity began to gradually deviate from linear growth, which shown the yield growth trend. The growth rate was weakened slightly at stage III. At the end of this stage, the maximum value of loading capacity developed to 7.50 MPa. Obviously, the limit bearing capacity of serial coal pillar specimens is similar to the UCS of lower coal pillar specimen. The activity of AE behavior became more intense at stage III. Massive AE energies released and rapid growth of accumulated AE energy was observed. The maximum value of AE energy was recorded as soon as the appearance of peak stress. During the stage III, the electrical resistivity maintained a very low level. It still basically presented a decreasing tendency as a whole. At the end of stage III, the electrical resistivity of serial coal pillar specimens increased to 2.09 MΩ·m sharply. It was attributed to the formation and expansion of micro-fracture in the testing coal pillar specimens [43,44]. The dilatation phenomenon was also observed at this moment. The electrical conductive channels then suffered serious damaged and broken [36,45]. As a result, the electrical resistivity presented an accelerated rising tendency.
(4) At stage IV, the carrying capacity of serial coal pillar specimens turned to decrease. At the initial period of stage IV, the weakening rate was relatively slow. The bearing capacity had a weak fluctuation change. Afterwards, a large drop of bearing capacity occurred within a short time. It indicated that the serial coal pillar specimens experienced intense damage and the system stability was seriously destructed. The output of AE energy continued to be intense at stage IV. It was primarily induced by the propagation of macrocracks. An increasing tendency of accumulated AE energy was monitored at this stage. As shown in Figure 5, the slope of accumulated AE energy, representing the damage rate, was larger enough. It was an illustration about the violent damage. Along with the weakening of bearing capacity, the change of electrical resistivity was relatively weak. This was due to the fact that irreversible damage occurred in the conductive channel of testing specimens after the formation of macro-cracks. It made the overall conductivity challenging to recover.

AE Events Distribution
The stress-AE energy-electrical resistivity response could reflect the overall fracturing behaviors of serial coal pillar specimens. However, it was difficult to explain the local fracturing behaviors of upper specimen or bottom specimen during the loading process. The interaction relationship between the upper specimen or bottom specimen was even more difficult to be revealed. AE features could not only be reflected by the evolution of AE energy quantitatively, but also be revealed by the distribution of AE events qualitatively [40,41]. Therefore, the AE events distribution of serial coal pillar specimens at different loading stages is shown in Figure 6. It is presented to reflect the micro-cracks and internal fractures in the upper specimen or bottom specimen [46,47]. The letters showed in Figure 6 are corresponding to the purple markers and time slices in Figure 5.
Only a sporadic amount of AE events appeared inside the upper coal pillar specimen. No AE events were distributed inside the bottom coal pillar specimen at stage I, as shown in Figure 6a,b. At stage II, earlier AE events appeared inside the left area of bottom coal pillar specimen (lower UCS). The number of AE ruptures inside the specimen increased step by step, which gradually extended to the upper-right area of bottom coal pillar specimen (lower UCS). In this process, the AE events also increased weakly inside the upper coal pillar specimen (higher UCS). The main distribution area was at the right surface area, as shown from Figure 6c-f. At stage III, the AE events distribution density and area increased rapidly inside the bottom coal pillar specimen. It spread with a relatively fast speed. Internal fracturing was still weak for the upper coal pillar specimen (higher UCS), as shown in Figure 6g,h. Almost no new-born AE events appeared [see Figure 6h]. At stage IV, the number and density of AE events scattered inside the bottom coal pillar specimen (lower UCS). The AE events were penetrated throughout the two end surfaces of bottom coal pillar specimen (lower UCS), as shown in Figure 6i,j. This showed that the lower UCS coal pillar specimen became to be failure. At this time, the changes of AE fracturing point in the upper higher UCS coal pillar specimen were still very weak. Only a small number of new rupture points appeared on the right-top edge of the specimen, indicating that the upper coal pillar specimen was hardly fracturing during the entire loading process. After that, the lower coal pillar specimen gradually lost the ability to support the upper coal pillar specimen. This would cause the overall destruction of serial coal pillars. internal fractures in the upper specimen or bo om specimen [46,47]. The le ers showed in Figure 6 are corresponding to the purple markers and time slices in Figure 5. Only a sporadic amount of AE events appeared inside the upper coal pillar specimen. No AE events were distributed inside the bo om coal pillar specimen at stage I, as shown in Figure 6a,b. At stage II, earlier AE events appeared inside the left area of bo om coal pillar specimen (lower UCS). The number of AE ruptures inside the specimen increased step by step, which gradually extended to the upper-right area of bo om coal pillar specimen (lower UCS). In this process, the AE events also increased weakly inside the upper The above distribution characteristics of AE events showed that: the bottom coal pillar specimen was earlier and more prone to damage. Then, it lost the support of upper coal pillar specimen, which also triggered the destruction of serial coal pillar system. The overall limit bearing capacity of serial coal pillar specimens was 7.50 MPa. It was close to the strength (7.62 MPa) of bottom coal specimen. The above phenomenon was appeared since the UCS of bottom coal pillar specimen was less than that of upper testing specimen. As a result, the systemic fracturing of serial coal pillar specimens was accompanied by the fracturing of lower UCS specimen. It could be concluded that the UCS difference between the upper and bottom specimens was one of the important factors influencing the overall bearing capacity and system stability.

Local Strain Develops
To make a further explanation about the fracturing of serial coal pillar specimens, the local strain was measured by the strain gauge rosettes. It was widely applied to reflect the local deformation of coal and rock [48,49]. Local strain of testing coal pillar specimens is plotted as Figure 7. The positive value of local strain represents the compressive strain (CS) measured by the longitudinal gauge, while the negative value of local strain means the tensile strain (TS) measured by the horizontal gauge.

Local Strain Develops
To make a further explanation about the fracturing of serial coal pillar specimens, the local strain was measured by the strain gauge rose es. It was widely applied to reflect the local deformation of coal and rock [48,49]. Local strain of testing coal pillar specimens is plo ed as Figure 7. The positive value of local strain represents the compressive strain (CS) measured by the longitudinal gauge, while the negative value of local strain means the tensile strain (TS) measured by the horizontal gauge.
(1) At stage I, CS and TS were very small for the testing coal pillar specimens. Nonlinear growth of local strain was observed as a whole.
(2) At stage II, the growth rate of local strain was more rapid for the bo om coal pillar specimen (lower UCS). It revealed that internal fracturing inside the lower UCS coal pillar specimen was more intense. As a result, both the CS and TS of lower UCS coal pillar specimen were larger in value.
(3) At stage III, the stable growth of local strain was monitored for the upper coal pillar specimen (higher UCS). At the peak stress point, the maximum values of CS and TS were captured in the higher UCS coal pillar specimen. (1) At stage I, CS and TS were very small for the testing coal pillar specimens. Nonlinear growth of local strain was observed as a whole.
(2) At stage II, the growth rate of local strain was more rapid for the bottom coal pillar specimen (lower UCS). It revealed that internal fracturing inside the lower UCS coal pillar specimen was more intense. As a result, both the CS and TS of lower UCS coal pillar specimen were larger in value.
(3) At stage III, the stable growth of local strain was monitored for the upper coal pillar specimen (higher UCS). At the peak stress point, the maximum values of CS and TS were captured in the higher UCS coal pillar specimen.
(4) At stage IV, a rapid growth of CS and TS was observed for the bottom coal pillar specimen. Then, the maximum CS and TS values of above specimen were reached. Afterwards, CS and TS began to decline for the upper coal pillar specimen (higher UCS). It was worth to noted that the weakening rate was basically contrary to the growth rate at stage III. That is, the rebounding of deformation was appeared in the corresponding specimens. At the end of stage IV, local strain of higher UCS coal pillar specimen dropped suddenly along with the losing of bearing stress.
In summary, the uniaxial mechanical behaviors were mainly divided into four stages for the serial coal pillar specimens. The AE features, electrical resistivity responses and local strain develops were highly consistent to illustrate the mechanical behaviors. The lower UCS coal pillar specimen was more likely to suffer internal fracturing throughout the loading progress. System fracturing of serial coal pillar specimens was appeared along with the collapse of lower UCS coal pillar specimen. The limit bearing capacity of serial coal pillar specimens with different UCS was almost equal to the strength of lower UCS coal pillar specimen.

Numerical Simulation
To make a further demonstration about experimental results' accuracy, the numerical fracturing behaviors were investigated by particle flow code (PFC) for the serial coal pillars with different strength.

Establishment of Numerical Model
PFC is a discrete element numerical simulation software that is widely used in the study of the mechanical behavior of coal and rock. The numerical model of serial coal pillar specimens was shown in Figure 8. In accordance with laboratory tests, a model of the specimen was constructed with a diameter of 50 mm and a height of 200 mm, using the same scale as the experimental specimen. The particle size was set to 0.5-0.75 mm. The specimens constituted 7217 particles. The model consisted of bottom coal pillar specimen C1 and upper coal pillar specimen C2. Parallel bonding model was used for the contact model of the specimens to reflect the mechanical properties. Table 3 showed the meso-parameters used in the PFC model for serial coal pillar specimens in this research. coal pillar specimen.

Numerical Simulation
To make a further demonstration about experimental results' accuracy, the numerical fracturing behaviors were investigated by particle flow code (PFC) for the serial coal pillars with different strength.

Establishment of Numerical Model
PFC is a discrete element numerical simulation software that is widely used in the study of the mechanical behavior of coal and rock. The numerical model of serial coal pillar specimens was shown in Figure 8. In accordance with laboratory tests, a model of the specimen was constructed with a diameter of 50 mm and a height of 200 mm, using the same scale as the experimental specimen. The particle size was set to 0.5-0.75 mm. The specimens constituted 7217 particles. The model consisted of bo om coal pillar specimen C1 and upper coal pillar specimen C2. Parallel bonding model was used for the contact model of the specimens to reflect the mechanical properties. Table 3 showed the mesoparameters used in the PFC model for serial coal pillar specimens in this research.   Table 3. Meso-Parameters of serial coal pillar specimens with different strength.

Meso-Mechanical Parameters Values
Particle minimum radius, (mm) 0.5 Particle radius ratio 1.5 The density of the particle, (kg/m 3 ) 2500 The damp of the particle 0. Three types of contact bonding models were applied to simulate the serial coal pillar specimens: the contact bonding in specimen C1, the contact bonding in specimen C2, and the contact bonding at the interface of C1 and C2. Figure 9 shows the comparison of stress-strain curves obtained from experimental and numerical serial coal pillar specimens under the uniaxial compression. It could be seen that the numerical results agreed well with the experimental results.
Three types of contact bonding models were applied to simulate the serial coal pillar specimens: the contact bonding in specimen C1, the contact bonding in specimen C2, and the contact bonding at the interface of C1 and C2. Figure 9 shows the comparison of stressstrain curves obtained from experimental and numerical serial coal pillar specimens under the uniaxial compression. It could be seen that the numerical results agreed well with the experimental results.

Stress Distribution Pa ern
To obtain the stress evolution pa ern, 1200 measurement circles were regularly arranged inside the serial specimens during the loading process. The measurement circles had a certain amount of overlapping to ensure the particles in the range were monitored completely. In the loading direction, the numerical stress evolution pa ern of serial coal pillar specimens was shown in Figure 10. Six representative stress state were presented Figure 10. The stress of point A state was 20% of maximum bearing capacity. The stress of point B state was 40% of maximum bearing capacity. The stress of point C state was 60% of maximum bearing capacity. The stress of point D state was 80% of maximum bearing capacity. The stress of point E state was equal to the maximum bearing capacity. The stress of point F state was 35% of maximum bearing capacity after the peak stress. It should be

Stress Distribution Pattern
To obtain the stress evolution pattern, 1200 measurement circles were regularly arranged inside the serial specimens during the loading process. The measurement circles had a certain amount of overlapping to ensure the particles in the range were monitored completely. In the loading direction, the numerical stress evolution pattern of serial coal pillar specimens was shown in Figure 10. Six representative stress state were presented Figure 10. The stress of point A state was 20% of maximum bearing capacity. The stress of point B state was 40% of maximum bearing capacity. The stress of point C state was 60% of maximum bearing capacity. The stress of point D state was 80% of maximum bearing capacity. The stress of point E state was equal to the maximum bearing capacity. The stress of point F state was 35% of maximum bearing capacity after the peak stress. It should be noted that the positive stress values indicated the tensile stress, while the negative stress values meant the compressive stress.  It could be seen from Figure 10 that the compressive stress existed within the serial coal pillar specimens before the peak stress. At point A, the compressive stress was relatively low and the stress concentrations were found on the interface of upper-bo om specimens. From point B to point E, the stress concentration degree decreased. With the It could be seen from Figure 10 that the compressive stress existed within the serial coal pillar specimens before the peak stress. At point A, the compressive stress was relatively low and the stress concentrations were found on the interface of upper-bottom specimens. From point B to point E, the stress concentration degree decreased. With the growth of loading, the stress became progressively uniform within the serial specimens. At point F, the non-uniform tensile stress zone appeared inside the serial coal pillar specimens. It indicated that the tensile stress was accompanied with the post-peak stage. However, no significant differences were observed in the stress distribution between the upper and lower coal pillar specimens. That is, the serial coal pillar specimens had a balanced bearing characteristic, which could also be confirmed by the following stress data (see Figure 11). terials 2023, 16, x FOR PEER REVIEW Figure 11. Quantitative relationship between the stress and timestep of serial Obviously, the carrying capacity of the upper testing specimen specimen was basically equal in value. The overall bearing capacity of system was also equal to the independent carrying capacity of testing s there were many differences in the contact force chain distribution bet specimen and bo om specimen. It could be seen from Figure 12 that t tion, extension and penetration of contact force chain mainly occurre pillar specimen (lower UCS). When the macro-contact force chain frac in the bo om coal pillar specimen, local instability was caused. It wou for the upper coal pillar specimen (higher UCS). Eventually, the sy serial coal pillar specimens was induced.  Figure 11 presented the quantitative stress-time relationship of upper testing specimen, bottom testing specimen and serial coal pillar specimens, which were expressed as the blue curve, red curve and black curve. Figure 12 represents the contact force chain distribution of serial coal pillar specimens during the loading process. The step in Figure 12 was corresponding to the purple time circle in Figure 11.
Obviously, the carrying capacity of the upper testing specimen and bottom testing specimen was basically equal in value. The overall bearing capacity of the series coal pillar system was also equal to the independent carrying capacity of testing specimen. However, there were many differences in the contact force chain distribution between the upper coal specimen and bottom specimen. It could be seen from Figure 12 that the rupture, bifurcation, extension and penetration of contact force chain mainly occurred in the bottom coal pillar specimen (lower UCS). When the macro-contact force chain fracture zone appeared in the bottom coal pillar specimen, local instability was caused. It would lose the support for the upper coal pillar specimen (higher UCS). Eventually, the systemic fracturing of serial coal pillar specimens was induced.
there were many differences in the contact force chain distribution between the upper coal specimen and bo om specimen. It could be seen from Figure 12 that the rupture, bifurcation, extension and penetration of contact force chain mainly occurred in the bo om coal pillar specimen (lower UCS). When the macro-contact force chain fracture zone appeared in the bo om coal pillar specimen, local instability was caused. It would lose the support for the upper coal pillar specimen (higher UCS). Eventually, the systemic fracturing of serial coal pillar specimens was induced.

Crack Evolution
From a qualitative and quantitative perspective, the crack number, form, and location of serial coal pillar specimens were investigated in the following section. Figure 13 and Table 4 showed the cracks number evolution of upper specimen, bo om specimen and serial coal pillar specimens. Figure 14 represented the cracks form and location of upper specimen, bo om specimen and serial coal pillar specimens at different loading timestep. It should be pointed out that the blue lines in the Figure 14 indicated the tensile cracks, and the red lines meant the shear cracks.

Crack Evolution
From a qualitative and quantitative perspective, the crack number, form, and location of serial coal pillar specimens were investigated in the following section. Figure 13 and Table 4 showed the cracks number evolution of upper specimen, bottom specimen and serial coal pillar specimens. Figure 14 represented the cracks form and location of upper specimen, bottom specimen and serial coal pillar specimens at different loading timestep. It should be pointed out that the blue lines in the Figure 14 indicated the tensile cracks, and the red lines meant the shear cracks.
From a qualitative and quantitative perspective, the crack number, form, and location of serial coal pillar specimens were investigated in the following section. Figure 13 and Table 4 showed the cracks number evolution of upper specimen, bo om specimen and serial coal pillar specimens. Figure 14 represented the cracks form and location of upper specimen, bo om specimen and serial coal pillar specimens at different loading timestep. It should be pointed out that the blue lines in the Figure 14 indicated the tensile cracks, and the red lines meant the shear cracks. Figure 13. Evolution curve of crack number of serial coal pillar specimens.
When the loading timestep was 15,000~20,400, the number of micro-cracks in the upper specimen and bottom specimen increased slowly, but it was mainly sprouted and expanded in the bottom specimen, as shown in Figure 14. Specifically, the number of tensile cracks in the upper specimen grew from 6 to 22, and the number of shear cracks increased from 1 to 2. In the bottom specimen, the tensile cracks number increased to 58 from 3, while the shear crack number grew to 31 from 1. Moreover, the micro-cracks were mainly distributed in the lower-left corner and upper-right corner area in the bottom specimen, as shown in Figure 14b-d.
When the loading timestep was 20,400~22,800, the micro-cracks in the upper specimen still grew slowly. The number of tensile cracks and shear cracks increased to 50 and 3, respectively. However, the micro-cracks in the bottom specimen showed a rapid growth trend, as shown in Figure 13 and Table 4. The number of tensile cracks in the bottom specimen increased from 58 to 869, and the number of shear cracks rose from 31 to 220. Meanwhile, the tensile and shear cracks were still mainly distributed lower-left corner and upper-right corner area in the bottom specimen. The distribution area also gradually increased. That is, the micro-cracks in the above two distribution areas gradually extended and gathered. However, the two areas did not penetrate each other. In the end, a W shaped cracks presented in the lower-left corner, and a triangular shaped cracks gathered in the upper-right corner, as shown in Figure 14d- was accompanied. Then, the overall fracturing of serial coal pillar specimens was also induced.  Figure 15 represented locally the displacement-timestep relationship of upper and bo om specimens. It can be seen from Figure 15 that: the vertical displacement of bo om specimen was always larger than that of upper specimen during the loading process. There were large differences in the displacement generated by the upper and the bo om specimen. That is, the serial coal pillar specimens showed obvious non-balanced deformation characteristics under the axial loading condition. When the loading timestep was 22,800~23,487, the number of micro-cracks remained unchanged inside the upper specimen, as shown in Figure 13 and Table 4. The micro-cracks continued to increase sharply in the bottom specimen. The number of tensile cracks and shear cracks increased to 1156 and 266, respectively. It should be pointed out that the micro-cracks gradually expanded at the lower-left corner and upper-right corner. As a result, the cracks were penetrated throughout each other and formed several macro cracks, as shown in Figure 14j-l. It also indicated the macro-failure of bottom specimen, which lost the support to the upper specimen. It was believed that the system stability could not be maintained for the serial coal pillar stability.

Local Displacement
In summarize, the germination, extension and penetration of micro-cracks occurred mainly in the specimen with lower UCS. Local instability appeared of lower UCS specimen when the macro-cracks generated. Then, the chain fracturing of higher UCS specimen was accompanied. Then, the overall fracturing of serial coal pillar specimens was also induced. Figure 15 represented locally the displacement-timestep relationship of upper and bottom specimens. It can be seen from Figure 15 that: the vertical displacement of bottom specimen was always larger than that of upper specimen during the loading process. There were large differences in the displacement generated by the upper and the bottom specimen. That is, the serial coal pillar specimens showed obvious non-balanced deformation characteristics under the axial loading condition.

Discussion
The related numerical results and experimental results were consistent as a The common conclusion could be obtained that the unbalanced deformation cha tics of serial coal pillar specimens was presented. The balanced bearing charac were also exhibited. In order to make a deeper illustration about the related pheno a schematic model of serial coal pillar specimens was established in the following Then, the progressive fracturing mechanism of serial coal pillar specimens with d UCS was investigated. Figure 16 showed the schematic model of serial coal pillar specimens with d strength. In the one-degree-of-freedom system, two nonlinear and inelastic co specimens were connected in serial [20]. Furtherly, the upper specimen and bo o imen were simplified to two springs, which was represented by Su and Sb, respe That is, the serial coal pillar specimens were regarded to be composed by two in nected springs. Both of them were possessed with the softening and hardening ch istics. For the serial coal pillar specimen model, the UCS of upper pillar specim more than the bo om one, while the elastic modulus was same in value.

Discussion
The related numerical results and experimental results were consistent as a whole. The common conclusion could be obtained that the unbalanced deformation characteristics of serial coal pillar specimens was presented. The balanced bearing characteristics were also exhibited. In order to make a deeper illustration about the related phenomenon, a schematic model of serial coal pillar specimens was established in the following section. Then, the progressive fracturing mechanism of serial coal pillar specimens with different UCS was investigated. Figure 16 showed the schematic model of serial coal pillar specimens with different strength. In the one-degree-of-freedom system, two nonlinear and inelastic coal pillar specimens were connected in serial [20]. Furtherly, the upper specimen and bottom specimen were simplified to two springs, which was represented by S u and S b , respectively. That is, the serial coal pillar specimens were regarded to be composed by two interconnected springs. Both of them were possessed with the softening and hardening characteristics. For the serial coal pillar specimen model, the UCS of upper pillar specimen was more than the bottom one, while the elastic modulus was same in value.

Schematic Model of Serial Coal Pillar Specimens
strength. In the one-degree-of-freedom system, two nonlinear and inelastic coal pillar specimens were connected in serial [20]. Furtherly, the upper specimen and bo om specimen were simplified to two springs, which was represented by Su and Sb, respectively. That is, the serial coal pillar specimens were regarded to be composed by two interconnected springs. Both of them were possessed with the softening and hardening characteristics. For the serial coal pillar specimen model, the UCS of upper pillar specimen was more than the bo om one, while the elastic modulus was same in value. Figure 16. Schematic model of serial coal pillar specimens before and after the uniaxial loading.
For the upper and bo om coal pillar specimen, the relationship between the bearing load and displacement could be expressed by Equation (3) [14,21]. For the upper and bottom coal pillar specimen, the relationship between the bearing load and displacement could be expressed by Equation (3) [14,21].
where F u and F b meant the bearing load of upper and bottom coal pillar specimen, respectively. u u and u b represented the vertical displacement of upper and bottom coal pillar specimen, respectively. When the interaction between upper and bottom pillar specimens was regarded to be quasi-static, the constitutive relation of serial coal pillar system could be described by Equation (4).
where F s was the global bearing load of serial coal pillar specimens. u s was the total displacement of serial coal pillar specimens. Furthermore, the relationship between bearing load and displacement for serial coal pillar specimens could be derived as the Equation (5) [14]: where f s (u s ) was the secant modulus of serial coal pillar specimens, f u (u u ) and f b (u b ) were the secant modulus of upper and bottom coal pillars, respectively.

Fracturing Evolution Analysis
According to the fracturing behaviors of serial coal pillar specimens, the deformation of upper and bottom pillar specimens was not always coordinated during the loading process. Figure 17 presented the schematic relationship of load-vertical displacement for serial coal pillar specimens (different strength) during the whole loading process. F u and F b mean the limit bearing load of upper and bottom coal specimen, respectively. The load-displacement relation of lower UCS coal pillar specimen was expressed by the red curve, while it was plotted by the purple curve for the higher UCS coal pillar specimen. For the dashed purple curve in Figure 17, it stood for the hypothetical variation path of load-displacement for higher UCS coal pillar specimen when the bearing load was more than F b . Without the effect of lower UCS coal pillar specimen, it would experience the strain hardening and softening stage when the stress was more than F b . The solid purple curve in Figure 17 represented the actual evolution path of higher UCS coal pillar specimen during the loading process. Obviously, there was a significant rebound deformation when the stress was more than F b . Moreover, the blue curve represented the global load-displacement relation of serial coal pillar specimens.
Materials 2023, 16, x FOR PEER REVIEW Figure 17. Relationship of load-displacement for serial coal pillar specimens (different stre From Figure 17, the unbalanced deformation characteristics were illustrate again for the serial coal pillar specimens. Therefore, the key deformation element fined to reveal the unbalanced deformation characteristics in this study. It was dete by comparing the values of uu and ub. Both the testing coal pillar specimens were re as the key deformation elements when uu = ub. The upper coal pillar specimen, gen more deformation than the bo om one, was considered as the key deformation e Contrarily, the key deformation element was the bo om coal pillar specimen wh ub. Specifically, the fracturing evolution of serial coal pillar specimens was analyz deformation of upper and bo om testing specimens at the initial loading period sically coordinated. Both the upper and bo om coal pillar specimens were regarde key deformation elements. Then, the bo om coal pillar specimen (lower UCS) tu be the key deformation element. Owing to the lower strength, the axial stress wa to reach the limit bearing load of bo om coal pillar specimen. As a result, the fra evolution of lower UCS coal pillar specimen was presented as the solid path: O → Bb → Cb → Db' → Eb. At that same time, the higher UCS coal pillar specimen was sti elastic deformation period. Afterward, the bearing load of lower UCS coal pillar sp was gradually weakened. It would induce the elastic rebounding deformation of UCS coal pillar specimen. Consequently, the fracturing evolution of higher UCS co specimen was presented as the solid path: O → Au → Bu' → Cu' → Du' → Eu'.

Progressive Fracturing Mechanism
The fracturing of serial coal pillar specimens was progressive. Table 5 summ the key deformation element evolution and unbalanced bearing characteristics at d the loading stages. The following results could be obtained: Table 5. Progressive fracturing evolution of serial coal pillar specimens (different strength) From Figure 17, the unbalanced deformation characteristics were illustrated once again for the serial coal pillar specimens. Therefore, the key deformation element was defined to reveal the unbalanced deformation characteristics in this study. It was determined by comparing the values of u u and u b . Both the testing coal pillar specimens were regarded as the key deformation elements when u u = u b . The upper coal pillar specimen, generating more deformation than the bottom one, was considered as the key deformation element. Contrarily, the key deformation element was the bottom coal pillar specimen when u u < u b .
Specifically, the fracturing evolution of serial coal pillar specimens was analyzed. The deformation of upper and bottom testing specimens at the initial loading period was basically coordinated. Both the upper and bottom coal pillar specimens were regarded as the key deformation elements. Then, the bottom coal pillar specimen (lower UCS) turned to be the key deformation element. Owing to the lower strength, the axial stress was earlier to reach the limit bearing load of bottom coal pillar specimen. As a result, the fracturing evolution of lower UCS coal pillar specimen was presented as the solid path: At that same time, the higher UCS coal pillar specimen was still at the elastic deformation period. Afterward, the bearing load of lower UCS coal pillar specimen was gradually weakened. It would induce the elastic rebounding deformation of higher UCS coal pillar specimen. Consequently, the fracturing evolution of higher UCS coal pillar specimen was presented as the solid path: O → A u → B u ' → C u ' → D u ' → E u '.

Progressive Fracturing Mechanism
The fracturing of serial coal pillar specimens was progressive. Table 5 summarized the key deformation element evolution and unbalanced bearing characteristics at different the loading stages. The following results could be obtained:

Conclusions
In this paper, uniaxial compression tests of serial coal pillar specimens with different strengths were conducted. The "stress-acoustic-electric" response characteristics of specimens during the failure process were analyzed. The PFC2D software was used to analyze the force chain evolution and crack development characteristics. The main conclusions are as follows: (1) The bearing stress of serial coal pillar specimens are mainly divided into four stages during the loading process: nonlinear growth stage of bearing capacity, linear growth stage of bearing capacity, yielding growth stage of bearing capacity, weakening stage of bearing capacity. The AE features, electrical resistivity responses, strain develops, force chain distributions, cracks evolutions, and local displacement are highly consistent to illustrate the fracturing behaviors.
(2) The lower UCS coal pillar specimen is more likely to suffer internal damage throughout the loading progress. System fracturing of serial coal pillar specimens is appeared along with the collapse of lower UCS coal pillar specimen. The limit bearing capacity of serial coal pillar specimens is almost equal to the strength of lower UCS coal pillar specimen.
(3) Due to the differences of UCS, the unbalanced deformation characteristics of serial coal pillar specimens are presented. The balanced bearing characteristics are also exhibited. The key deformation element is defined to reflect the unbalanced deformation characteristics. The evolution of key deformation element is the rooted reason to account for the overall fracturing mechanism of serial coal pillar specimens.
(4) For the serial coal pillar specimens with different UCS, both the lower and higher UCS coal pillar specimens are the initial key deformation element. It turned out to be the lower UCS coal pillar specimen at the later loading stages. The limit bearing capacity of lower UCS coal pillar specimen reached earlier, while higher UCS coal pillar specimen is still at the elastic deformation stage. Rebounding deformation of higher UCS coal pillar specimen occurred after the peak stress. It would accelerate the system fracturing of lower UCS coal pillar specimen. The limit carrying capacity is almost equal to the strength of lower UCS coal pillar specimen.
(5) The critical condition of system fracturing is f u (u u ) + f b (u b ) = 0 for serial coal pillar specimens with different strength, which is expect to predict the system fracturing of serial pillar in the underground coal mining.