1 Introduction
In coal mining operations, numerous coal pillars are typically left to stabilize the overlying rock layer in the extraction zone, forming a coal-rock composite structure. When neighboring coal seams have narrow layer spacing, cracks in the surrounding rock above the lower coal seam increase the likelihood of collapses or toppling, significantly influencing the development of impact ground pressure. Thus, studying the mechanical properties of coal-rock composites requires careful consideration of the roof rock layer's initial damage state. This consideration is critical for designing mine support structures and provides a robust theoretical foundation for engineering applications [1–4].
Numerous scholars have examined the mechanical properties and crack evolution of prefabricated cracks in coal-rock masses. Chen et al. [5] observed that penetrating cracks between coal and rock layers significantly influence the mechanical properties of coal-rock systems. The experiments revealed that crack length and angle play a critical role in strength, energy evolution, and destabilization mechanisms. Zhao et al. [6] studied the strength-deformation characteristics and crack extension laws in cracked rock masses under combined water and stress conditions. They found that increasing water pressure transformed crack-bearing rock masses from ductile to brittle damage, with crack angle significantly affecting mechanical properties and damage morphology. Wang et al. [7] used CT scanning and PFC simulations to analyze crack evolution in coal masses containing prefabricated cracks under varying perimeter pressures. They observed that increasing perimeter pressure reduced crack initiation angles and shifted damage modes from tensile-shear to shear failure. Their proposed stress-strength factor formula closely aligned with experimental results. Li et al. [8] examined how prefabricated crack angles affect mechanical properties, damage modes, and acoustic-thermal characteristics of coal specimens through uniaxial compression tests. They found that crack angle significantly influenced peak stress and strain, revealing crack expansion through infrared thermal imaging and acoustic emission analysis. Wang et al. [9] analyzed the mechanical properties of granite with single penetrating cracks. They reported that increasing crack angle significantly reduced compressive strength, damage stress, and elastic modulus, with damage modes transitioning from crack-surface failure to sliding failure along the crack. Li et al. [10] analyzed the effects of various parallel pre-cracking forms on the mechanical properties of coal-rock composites using PFC2D numerical simulations. They found that parallel pre-cracking significantly reduced peak strength, peak strain, elastic modulus, and crack initiation stress, identifying three primary crack initiation modes. Ma et al. [11] examined the mechanical behavior and failure mechanisms of coal-rock composites with different crack angles using uniaxial compression tests and numerical simulations. They found that crack angle significantly influenced strength, deformation, and damage modes, and proposed a failure criterion consistent with experimental results. Chen et al. [12] discovered that water-rock interactions accelerated crack expansion between coal and rock layers, resulting in mechanical property degradation and structural instability. Crack width and soaking time significantly influenced strength, energy characteristics, and damage mechanisms, posing risks to coal mining safety. Zhang et al. [13] employed RFPA2D software to simulate coal-rock composites with coal-sustained joints. They found that joint angle and height ratio significantly influenced uniaxial compressive strength, acoustic emission characteristics, and damage modes. The study proposed that acoustic emission energy accumulation curves could serve as precursors for failure. Li et al. [14] performed dynamic impact tests on pre-cracked coal-rock composites using a split Hopkinson pressure bar. They observed that crack location and angle significantly influenced strength and failure modes, with the most pronounced effects occurring at 30°. They further noted that cracks altered stress concentration and stability. Lv et al. [15] examined the mechanical properties, acoustic emission behavior, and damage modes of coal-rock composites with a single prefabricated non-penetrating crack through theoretical analysis, experiments, and numerical simulations. They concluded that crack angle significantly influenced unconfined compressive strength and failure modes, with strength reaching its lowest at 45°. Wang et al. [16] investigated how layering and prefabricated cracks influence cracking characteristics and breakage mechanisms in black shale using Brazilian splitting tests. They determined that prefabricated cracks weakened specimen strength, while crack initiation and expansion paths were significantly affected by layering angles and crack orientations. The study summarized five failure modes and two acoustic emission counting features. Li et al. [17] investigated how various defect types influence crack initiation, propagation, and damage modes in rock specimens through experiments and numerical simulations. They observed that wing-like cracks originated at defect tips and that crack agglomeration intensified as defect numbers increased. The study summarized the characteristics of cleavage and shear damage in rocks. Zhao et al. [18] explored the cracking and stress-strain behavior of rock-like materials with two defects using uniaxial compression tests and strain analysis. They found that crack merging modes were strongly correlated with defect geometries and that local strain concentration revealed crack development processes. Zeng et al. [19] developed a stable water pressure testing system to study fracture characteristics and acoustic emission behavior of fractured rock under 1 MPa water pressure coupled with stress. Their results indicated that water pressure significantly influenced peak strength, crack extension paths, and acoustic emission characteristics. Wang et al. [20,21] investigated the effect of CO2 fracturing on rocks and proposed a fractal damage calculation model. The experiments show that the static load guides the fracture direction, and the dynamic load affects the number and morphology of fracture, which provides a new support for the design of drilling vibration prevention and the study of the damage intrinsic properties of coal rock. Li et al. [22] proposed a two-dimensional eigenstructural model to effectively describe the nonlinear stress-strain behavior of rocks, which provides theoretical support for the damage analysis of coal-rock composites.
In summary, numerous studies on coal-rock composites have primarily focused on the mechanical properties and crack evolution characteristics influenced by prefabricated crack angles and sizes in single materials or interfaces [23,24]. However, the effects of internal cracks within the rock mass have received less attention. This study integrates laboratory experiments with DEM numerical simulations to model the initial damage state of the top rock layer by varying the geometrical parameters of prefabricated cracks within the rock. The results examine the damage modes and evolution processes of cracked coal-rock composites and establish a crack damage constitutive equation, offering theoretical support for stability assessments of coal-rock composites in coal mining.
1 Uniaxial compression experiment
1.1 Specimen preparation
In this study, pure coal, sandstone, and coal-rock composites were chosen as the experimental materials. The specimens were machined into standard cylindrical shapes with a diameter of 50 mm and a height of 100 mm. The pure coal specimens were collected from the Hepan coal mine located in Shenmu County, Shanxi Province. After on-site sampling, the specimens were cut and polished to prepare standard samples. The sandstone specimens were prepared by first extracting a core from a rock slab using a core drilling machine, followed by processing into a standard cylinder using a stone grinding machine. To maintain uniformity in structure and composition, all sandstone specimens were obtained from the same batch of rock slabs.
The preparation steps for the coal-rock composites are as follows: the treated sandstone and pure coal were cut transversely along the center and divided equally into two specimens of the same size. The cutting surfaces were sandpapered to ensure smoothness and uniformity. Subsequently, high-strength mica glue was used to bond the polished coal-rock interface, ensuring a tight fit, and the specimens were processed into a φ50 mm × 100 mm composite cylinder. After bonding, the composite specimens were placed in a cool environment for 24 hours to ensure adequate curing. The end faces of all specimens were finely polished to control non-parallelism and non-perpendicularity within 0.02 mm [25].
In addition, prefabricated cracks in the coal-rock composites were generated by wire-cutting. According to Yuan Chao [26] and Zhu et al. [27], the crack angle of approximately 45° has a greater effect on the strength of the specimen, while the effect is less when it is close to 90°. Therefore, the angles of the prefabricated cracks were chosen to be 45° and 90° for this experiment, and the size of the cracks was set to be 25 mm in length and 1.5 mm in width. To minimize the variability in experimental results, each specimen was subjected to three repetitions of the experiment, and the dimensions and mass of the specimens were accurately measured using vernier calipers and electronic scales. The physical parameters obtained are detailed in Table 1.
[Figure omitted. See PDF.]
1.2 Experimental results
The uniaxial compression experiments were conducted using the WDW-300 microcomputer-controlled universal testing machine produced by Changchun Kexin Co., Ltd. The displacement control method was employed for loading, with the loading rate of 0.12 mm/min [13].
The stress-strain curves obtained from the experiments are shown in Fig 1. One stress-strain curve from each group of specimens was selected to classify the stages accordingly. The curves can generally be divided into four stages: the initial compaction stage, the linear elastic deformation stage, the yield stage, and the failure stage. The curve behaves differently at each stage: ① Initial compaction stage (I stage): Due to the presence of primary cracks within the specimen, these cracks close under axial compression. This is represented in the graph as the “concave” portion of the curve, where the stress increases rapidly, but the deformation is relatively small. ② Linear elastic deformation stage (L stage): At this stage, the relationship between stress and strain follows Hooke’s law, meaning stress is proportional to strain, and the figure shows a linear relationship. ③ Yield stage (Y stage): When the stress reaches a certain value, the material enters the yield stage. At this stage, the specimen begins to undergo plastic deformation, and the deformation becomes irreversible. The figure shows the yield stage up to the peak stress. ④ Failure stage (F stage): At this point, the stress reaches its maximum value, fracture occurs within the material, the specimen is no longer intact, and the stress begins to decrease. The curve in the figure shows the post-peak stage.
[Figure omitted. See PDF.]
The results indicate that sandstone exhibits the highest peak strength and elastic modulus, with average values of 39.63 MPa and 4.82 GPa, respectively. The peak strength and elastic modulus of pure coal rank second, with an average peak strength of 16.28 MPa and an average elastic modulus of 1.82 GPa. The mechanical parameters of the coal-rock composites are between those of sandstone and pure coal, which is consistent with the conclusion obtained in the literature [28]. The average peak strength of the non-crack coal-rock composites is 23.30 MPa, and the average elastic modulus is 2.64 GPa. The analysis shows that the different material properties of coal and sandstone, as well as the weakening effect of the contact interface between coal and rock, may affect the stress transfer and failure mode, leading to differences in the stress-strain curve. The data from the uniaxial compression tests are summarized in Table 1. The letters CC in the table represent pure coal specimens, YS represents sandstone specimens, and MY represents coal-rock composite specimens without prefabricated cracks, e.g., MYS-45-1 represents coal-rock composite specimen No. 1 with a 45° crack, and so on.
2 Establishment of numerical model of fractured coal-rock mass
2.1 Numerical simulation
The particle flow code (PFC) is a software tool for modeling and analyzing the motion and interactions of particulate media using the discrete element method. The software effectively simulates the microscopic structural characteristics of coal rock mass, thus accurately reflecting their macroscopic mechanical behavior. In PFC2D, a variety of rigid particle contact models are provided. In this study, the parallel bond model (PBM) and the smooth joint model (SJM) are primarily used as microscopic contact models to analyze the mechanical properties of coal-rock mass.
Fig 2(a) illustrates the PBM model, which is capable of transmitting both force and moment. When the tensile or shear stress exceeds the tensile or shear strength, the particle connections in the PBM model break, thereby more realistically reflecting the mechanical behavior of the coal rock mass. In PFC simulations, the conventional approach to joint simulation is typically achieved by removing the contact on both sides of the joint. However, this method may not adequately model the natural roughness of the crack surface [29]. To address this issue, Deisman [30] proposed the SJM model, as shown in Fig 2(b). The model allows particles to interpenetrate and overlap each other in the direction of the cracks, which effectively improves the simulation accuracy of the mechanical properties of the joints. The model has been widely applied by scholars both domestically and internationally in the simulation study of joint cracks, with results highly consistent with experimental data [31–33].
[Figure omitted. See PDF.]
2.2 Calibration of microscopic parameters
PFC simulation uses local contact behavior to describe the macroscopic mechanical properties of materials, allowing numerical experiments to be conducted by setting the geometric and mechanical parameters of particles and bonds [34,35]. Fig 3 shows the comparison of uniaxial compression stress-strain curves obtained from both laboratory experiments and PFC simulations, revealing that the simulation curve closely matches the experimental data. However, in the PFC simulation, due to the rigidity of the particles in the numerical model, the ‘concave’ phenomenon observed in the actual crack compaction process during the initial loading stage is not present. In addition, the peak strength of coal is lower than that of sandstone, indicating that coal fails before sandstone under the same pressure, which is consistent with the results in the literature by Yang et al. [36,37]. During axial compression, the cracks generated inside the coal-rock composites not only expand from the coal body to the sandstone but also cause stress concentration at the crack tip. This stress concentration may cause the crack to further expand, connecting with other primary cracks and eventually leading to shear failure and scattered damage to the coal body. Therefore, the composite structure of sandstone and coal can effectively delay the formation and expansion of coal cracks by increasing the overall structural stability, thereby reducing the risk of composite failure.
[Figure omitted. See PDF.]
Therefore, it is necessary to consider three mesoscopic parameters: pure coal, sandstone, and the joint contact interface. The determination of these three parameters is typically calibrated using the ‘trial and error’ method. After repeated calibration, the simulation data become consistent with the experimental data, and the obtained microscopic parameters are shown in Table 2. The introduction of prefabricated cracks may lead to internal instability, which can eventually cause discrepancies in the stress-strain curves. Therefore, the simulation is based on the average value of the experimental data.
[Figure omitted. See PDF.]
In Table 3, the differences in macroscopic mechanical parameters between the average values obtained from the numerical model and the experimental tests are compared in detail. The results show that the simulated uniaxial compressive strength and elastic modulus closely match the experimental results. The relative errors between the experimental and simulated values were less than 10%. In addition, the failure modes of the two fracture numerical models are essentially the same as the actual failure modes observed in the experimental tests, indicating that the mesoscopic parameters and modeling methods used are suitable for related research.
[Figure omitted. See PDF.]
2.3 Simulation scheme
To study the effect of initial roof damage on the impact hazard of coal mining, this paper simplifies and quantifies the roof and coal seam, combining them into a coal-rock composite with a ratio of 1:1. Based on the actual working conditions of the upper rock roof and the lower coal seam at the site, a numerical model of the coal-rock composites with dimensions of φ50 mm × 100 mm was established using the numerical simulation software PFC2D.
In order to ensure that the crack features in the model are consistent with the crack characteristics observed in the laboratory experiments, the prefabricated cracks were generated by removing the particles. The crack length was fixed at L = 25 mm, and the crack widths were set at K = 1 mm, 1.5 mm, and 2.0 mm. Additionally, the crack angles were selected as 0°, 30°, 45°, 60°, and 90°, based on their common usage in the literature. These angles have been widely adopted in similar studies, such as those by Liu et al., Luo et al., Chen et al., Zhang et al., and Wen et al. [5,38–41], to investigate the mechanical behavior of fractured rock masses under various conditions. As an example, the coal-rock composite with 90° cracks is modeled, and the relevant model is shown in Fig 4.
[Figure omitted. See PDF.]
Uniaxial compression numerical simulation experiment was conducted on the coal-rock composite using a displacement loading control method, with a loading rate set at 0.12 mm/min. The aim was to explore the damage mechanism of the coal-rock composite in the presence of cracks and to analyze the effects of the crack inclination angle and crack width on the composite. The detailed simulation scheme is shown in Table 4. The number L25-K1-α0° indicates that the crack length of the sample is 25 mm, the crack width is 1 mm, and the inclination angle is 0°.
[Figure omitted. See PDF.]
3 Macro-micro mechanical characteristics of fractured coal-rock mass
3.1 Macroscopic mechanical properties
The stress-strain curves under uniaxial compression for the specimens of coal-rock mass with different crack types are shown in Fig 5(a), (b) and (c). Specimens with crack angles of 0°, 30°, 45°, 60°, and 90° exhibited slight oscillations near the peak strength due to crack propagation, and produced jagged horizontal fluctuations in the vicinity of the peak, along with some fluctuations both before and after the peak. This behavior is primarily attributed to local micro-cracks and closure caused by stress concentration at the crack tip, as well as dynamic changes in the microstructure and the redistribution of stress as the specimen approaches failure. These microscopic mechanisms have a certain influence on the macroscopic mechanical properties of the coal and rock mass.
[Figure omitted. See PDF.]
Fig 6(a) and (b) show the comparison of macroscopic mechanical parameters of coal and rock specimens with different crack widths under uniaxial compression. It is observed that the macroscopic parameters of the cracked coal-rock composites are lower than those of the non-fractured composites. This is because the fracture, as a defect within the material, reduces the overall continuity and uniformity of the coal-rock mass, leading to a decrease in its mechanical properties. As the fracture angle increases, coal and rock mass specimens with different fracture widths exhibit different mechanical behaviors.
[Figure omitted. See PDF.]
It can be observed that the compressive strength and modulus of elasticity of the cracked coal-rock composites exhibit significant fluctuating trends as the crack width and angle change. Firstly, the compressive strength shows a decreasing trend followed by an increasing trend with the crack angle. For the specimen with K = 1 mm, the compressive strength gradually decreases from 15.39 MPa at α = 0° to 15.17 MPa at α = 30°, then reaches a peak of 16.36 MPa at α = 45°, and finally increases to 17.10 MPa at α = 90°. The fluctuation in compressive strength intensifies with increasing crack width. For instance, the maximum value of 20.98 MPa is reached at α = 60° for K = 2 mm. Overall, larger crack widths and angles significantly influence the fluctuation of compressive strength, particularly with an increase in strength at larger angles, which may be related to the change in crack propagation pathways. The modulus of elasticity exhibits a ‘W’-shaped fluctuation trend, which becomes more pronounced, especially as the crack width increases. For the specimen with K = 1 mm, the modulus of elasticity decreases from 1.084 GPa at α = 0° to 1.148 GPa at α = 30°, peaks at 1.628 GPa at α = 45°, decreases to 1.346 GPa at α = 60°, and finally rises to 1.376 GPa at α = 90°. The fluctuation in the modulus of elasticity is more pronounced with K = 2 mm, reflecting the complex effect of prefabricated cracks on the mechanical properties at different angles.
Overall, the crack angle and width have a significant effect on the compressive strength and modulus of elasticity of the cracked coal-rock composites, and the presence of cracks alters the internal stress distribution and crack propagation paths, leading to periodic fluctuations in mechanical properties as the crack characteristics change. This provides an important theoretical foundation and guidance for geotechnical engineering design and rock crack prediction. The numerical simulation results are summarized in Table 5.
[Figure omitted. See PDF.]
3.2 Micromechanical response
In order to further explore the failure characteristics of coal-rock composites with prefabricated cracks, the meso-cracks are studied in greater detail. According to the failure mode, the meso-cracks can be classified into: ① wing-shaped crack (W); ② coplanar secondary cracks (C); ③ inclined secondary cracks (O) [42–44]. Based on the mode of force failure, the meso-cracks can be divided into: ① tensile crack (T); ② shear crack (S).
The failure diagram of the preset crack is shown in Fig 7. The failure mode, overall displacement vector diagram, and crack sketch of different specimens are presented in Fig 8. When the tensile stress and shear stress on the particles exceed the maximum tensile stress and maximum shear stress, respectively, the parallel bond is destroyed, resulting in tensile and shear cracks in the material. This is the crack generation mechanism in the PFC numerical simulation tests. In Fig 8, blue represents shear failure cracks, and red represents tensile failure cracks. Tensile-shear composite cracking occurs both within the rock material and at the contact interface. Tensile cracking predominates within the pure coal material, but as axial compression increases, crack evolution becomes more complex within the specimen. During the simulation, due to differences in the cracking mechanism and the shape of the prefabricated cracks, the cracks within the specimen are generally categorized into three types: S crack, T crack, and S-T mixed crack. By comparing the displacement vector diagram with the crack sketch, it is easier to observe that the tensile crack propagates along the stress direction in a relatively uniform manner in most cases, while secondary cracks form near the main crack zone.
[Figure omitted. See PDF.]
[Figure omitted. See PDF.]
As shown in Fig 8, it can be observed that the damage to coal-rock composites is primarily concentrated in the coal body and the interface, which aligns with the findings in the literature [45]. When the width of the prefabricated crack is kept constant, the upper rock material is damaged along the crack tip when the crack angle is 0° and 30°. With the transfer of stress, the lower coal body also experiences minor tensile failure. As the crack angle increases to 45°, the cracks formed at this point create a through-face as the prefabricated crack tip propagates through the coal-rock interface into the coal body. As the crack angle increases to 60° and 90°, the number of cracks in the upper rock material decreases, and the degree of damage gradually diminishes. However, the lower coal material experiences the opposite effect. The analysis suggests that increasing the crack angle alters the path of crack downward expansion, affecting the lower coal material and potentially inducing new crack formation or increasing the activity of primary cracks in the coal body during the extension process.
Combined with the analysis of the stress-strain curve in Fig 5, it can be seen that under the action of external load, not only tensile cracks are formed within the composite, but shear cracks also develop. This indicates that the instability and failure of the sample are caused by the combined action of tensile and shear forces, rather than a single failure mechanism. Therefore, when evaluating the structural stability of the sample, the interaction between shear and tensile failure must be considered. Focusing solely on one failure mechanism may result in an overestimation of the structural safety factor, thereby posing a greater safety risk.
4 Analysis of crack initiation stress and damage evolution of fractured coal-rock mass
4.1. Analysis of crack initiation stress of fractured coal-rock mass
In the process of axial loading, the numerical model is subjected to force in both directions, and the force is transmitted through particles. As the load gradually increases, when the ultimate strength of contact bonding is reached, the contact bonds between some particles break, generating cracks. As the load continues to increase, the fractures also expand synchronously, and the cumulative number of cracks reaches a new peak. The deterioration of the internal microstructure of the particles leads to further damage in the coal-rock composites, resulting in significant macroscopic damage. Fig 9, Fig 10 and Fig 11 show the crack number evolution curve and crack orientation accumulation diagram of the specimen under axial compression.
[Figure omitted. See PDF.]
[Figure omitted. See PDF.]
[Figure omitted. See PDF.]
Through the crack azimuth accumulation diagram in Fig 9, Fig 10 and Fig 11, it is observed that the cracks generated in the sample are primarily concentrated at azimuth angles of 90° and 270°, with a high consistency in the distribution density of cracks at these angles. This indicates that the cracks predominantly propagate along the axial compression direction.
In the study of the crack growth process of coal and rock mass, the development of cracks can be compared with the four stages of the uniaxial experimental stress-strain curve, providing a detailed understanding of the changes in the number of cracks and their influence: ① In the initial compaction stage, the crack increment is 0. This is because the natural primary cracks in the coal rock mass are compacted under the action of external pressure, and no new cracks are generated. ② With the increase in loading, the coal-rock mass transitions from elastic deformation to the linear elastic deformation stage. At this point, new secondary cracks begin to form, based on the existing primary cracks. The total number of cracks reaches about 100, indicating the initiation of crack activity. ③ During the yield stage, stress within the coal-rock mass reaches the yield strength of the material, leading to accelerated crack growth. The primary fractures and the newly formed secondary fractures further propagate under stress, beginning to converge in specific regions, forming one or more main fracture zones. The number of cracks increases sharply, rising from hundreds to between 400 and 1000, and the crack distribution begins to exhibit localized characteristics. ④ In the failure stage, the number of cracks reaches the peak, and further crack expansion leads to the complete destruction of the coal-rock composites. Localized cracks are formed, and the structural integrity of coal and rock mass is lost.
4.2. Establishment of damage evolution constitutive equation
In order to further study the influence of fracture width and fracture dip angle on coal-rock composites, a damage constitutive model is constructed to accurately characterize the relationship between these parameters. During the loading process of the composite, the development of damage is closely related to the generation and propagation of cracks, which can be quantified by introducing the damage factor D. As a key index, the number of cracks effectively reflects the distribution of micro-cracks within the sample and the corresponding damage level [46]. Based on this, assuming that the specimen is in good condition before compression, the damage factor D can be defined as [47,48]:
(1)
where D represents the damage factor under axial compression, ranging from 0 to 1, where 0 indicates no damage and 1 indicates complete failure; N is the number of cracks generated due to specimen damage, and N0 is the total number of cracks that can potentially be generated in the specimen.
It is assumed that the relationship between the probability of failure of micro-units and the cumulative number of cracks under axial compression follows the Weibull distribution statistical probability model [49,50]. As the crack area size increase, the number of micro-units in the composite specimen under axial compression is:
(2)
where N is the number of cracks in the specimen; p1 is the shape parameter, a dimensionless constant; p2 and p3 are the scale parameters in N-0.5, respectively.
It is assumed that there is the following relationship between the damage factor D and the number of cracks N generated at any time:
(3)
According to Eq. (4), the damage factor D can be expressed as:
(4)
If the cumulative number of cracks in the whole compression process of the specimen is Nt, N(t) represents the cumulative number of cracks at any time in the axial compression process of the specimen, then the relationship between Nt and N(t) is:
(5)
Therefore, the composites of Eq. (4) and Eq. (5) is transformed, and the damage factor expression of coal-rock composites specimen can be obtained:
(6)
Considering that the damage factor D, calculated based on the number of cracks, is highly consistent with the prediction results of the Weibull model, it is reasonable to use the number of cracks to characterize the damage evolution. Based on the effective stress theory presented in references [51, 52], the stress-strain relationship of coal rock damage can be expressed as follows:
(7)
where ε is the strain of the samples and E is the elastic modulus of the samples.
By substituting Eq. (6) into Eq. (7), the damage evolution constitutive equation can be calculated as follows:
(8)
The stress-strain, elastic modulus and fracture number are substituted into Eq. (8), and the parameters p1, p2 and p3 are solved by nonlinear fitting.
The specific parameters in the damage evolution constitutive equation of the coal-rock composites are listed, and the verification results of the constitutive model are shown in Fig 12. The solid line in the figure represents the simulated value, while the dashed line represents the theoretical value. It can be observed that, although the simulated curve exhibits fluctuates near the peak, the stress-strain curve calculated by the constitutive equation is generally consistent with the simulated curve in terms of trend, and maintains a good fit in the post-peak stage.
[Figure omitted. See PDF.]
According to the equations, the fracture area S and the fracture dip angle α can be used as independent variables, while the parameters p1, p2 and p3 are dependent variables that can be used to fit the coal-rock composites. The relationship between these variables is expressed in the following form:
(9)
4.3 Damage evolution analysis of cracked coal-rock composites
According to the constitutive relation in Eq. (9), the damage evolution equation of the cracked coal-rock composites, considering sensitivity factors such as the fracture inclination angle and fracture area, can be effectively characterized, as shown in Eq. (10). Thus, the damage evolution curve of the cracked coal-rock composites is obtained. Taking the K = 1.5 mm model as an example, the evolution curve is shown in Fig 13.
[Figure omitted. See PDF.]
It can be observed that: ① The damage factor D increases nonlinearly with the strain, and the fitting curve's R2 value is close to 1, indicating a good fitting effect. ② Under different inclination angles, the damage curve of coal-rock composites is consistent. The increase in axial pressure leads to the formation of new cracks, aggravating the damage, which then attenuates after the peak value. ③ In the initial loading stage, the slope of the damage curve is 0, indicating no cracks. As the loading progresses, the damage factor D increases slowly, and the cracks expand rapidly. In the peak stress stage, the damage becomes significant, accounting for 60%-80% of the failure process. The cracks connect, causing a sharp drop in stress. ④ Taking the specimen with a 45° inclination as a demarcation, it can be observed that the strain at the peak stress damage stage for the specimens with 0° and 30° inclination angles is higher than that for specimens with 60° and 90° inclination angles. This indicates that the specimens with larger inclination angles experience more severe damage and faster degradation. This also shows that increasing the fracture inclination angle can enhance the compressive performance of the coal-rock composites, but it accelerates the instability process of the coal-rock composite structure.
(10)
5 Conclusions and discussion
This study systematically analyzed the mechanical properties and crack damage evolution of coal-rock composites containing cracks under uniaxial compression, through both laboratory uniaxial compression experiments and numerical simulations. A constitutive equation for crack damage evolution was developed. The results demonstrate that the presence of fractures significantly influences the mechanical behavior of coal-rock composites, with particular emphasis on the role of crack inclination in regulating the compressive strength and failure mode of the composite material.
1. In uniaxial compression experiments on different specimens, sandstone exhibited the highest peak strength and elastic modulus, followed by pure coal, while the macro-mechanical parameters of the coal-rock composites were intermediate. The strength and stiffness of coal-rock composites containing cracks were significantly lower than those without cracks due to the presence of cracks. The crack angle and width exerted a significant influence on the mechanical behavior of the coal-rock composites: the strength and elastic modulus of K1 and K1.5 gradually increased with the increase of the inclination angle, while K2 exhibited a nonlinear change of increasing and then decreasing.
2. The failure modes of the coal-rock composites are primarily dominated by the coal component, exhibiting bulk failure, which indicates that the failure process is governed by the combined action of tension and shear. This emphasizes the need to consider the interaction of these two damage mechanisms in the assessment of structural stability. Analysis of the crack evolution curves and radial accumulation diagrams shows that the cracks are mainly concentrated at azimuthal angles of 90° and 270°, indicating that they primarily propagate along the axial compression direction. With the increase in crack inclination angle, the crack propagation path changes significantly, resulting in a reduced damage degree in the rock component and an intensified damage degree in the coal component, further highlighting the crucial role of the crack inclination angle in regulating the damage process of coal-rock composites.
3. The ontological equations and damage evolution curves based on the number of cracks show that there is a nonlinear relationship between the damage factor and the strain, and the shape of the damage curves remains unaffected by changes in the crack inclination angle. Taking the 45° inclination specimen as the demarcation, it can be observed that the strain at the peak stress stage for the specimens with 0° and 30° inclination angles is lower than that of the specimens with 60° and 90° inclination angles, indicating that while an increase in the crack inclination angle can improve the compressive performance of the coal-rock composites, it accelerates the process of structural instability.
These findings provide valuable insights for mine crack management. In practical applications, mining engineers should pay particular attention to the distribution and inclination of cracks, especially in regions with larger crack inclinations. More precise crack control and reinforcement measures should be implemented in these areas. For instance, enhanced crack monitoring and analysis, as well as techniques like grouting, could be employed to mitigate crack propagation. Moreover, adaptive support designs should be adopted in high-risk areas to improve the overall stability of the mine structure.
We recommend that when designing mine support structures, the characteristics of fractures particularly crack inclination, distribution, and failure modes be comprehensively considered. Suitable support materials, tailored to the fracture characteristics of coal-rock composites, should be selected, with adjustments made according to variations in crack inclination to ensure that the support system can effectively address different failure modes. Additionally, the monitoring and maintenance of the support structure should be strengthened to mitigate the risks associated with crack propagation and structural instability.
Supporting information
S1 Data. Raw data for figures.
https://doi.org/10.1371/journal.pone.0316586.s001
(ZIP)
S1 Fig. All figures presented in the study.
https://doi.org/10.1371/journal.pone.0316586.s002
(ZIP)
References
1. 1. Xie HP, Zhang ZT, Gao F, Zhang R, Gao MZ, Liu JF. Stress-fracture-seepage field behavior of coal under different mining layouts. Journal of China Coal Society. 2016;41(10):2405–17.
* View Article
* Google Scholar
2. 2. Xie H, Zhou H, Xue D, Gao F. Theory, technology and engineering of simultaneous exploitation of coal and gas in China. Journal of China Coal Society. 2014;39(08):1391–7.
* View Article
* Google Scholar
3. 3. Wang W, Chen YP, Yuan L, Chen RZ, Wang HF, Du K. Floor fracture evolution and relief gas drainage timeliness in deeper underground short-distance upper protective coal seam extraction. Journal of China Coal Society. 2016;41(01):138–48.
* View Article
* Google Scholar
4. 4. Ding K, Gu QH, Wang LG, Guo JX. Study on the coal rib burst mechanism in deep roadway and energy release. Journal of China Coal Society. 2024;49(09):3759–80.
* View Article
* Google Scholar
5. 5. Chen GB, Tang W, Li T, Wang CY, Wang EY, Zhang GH. Mechanical response and instability model of fractured coal-rock combined body. Rock and Soil Mechanics. 2024;45(09):2633–52.
* View Article
* Google Scholar
6. 6. Zhao C, Xin JQ, Niu JL, Ma CC. Experimental study on crack propagation of precrack rock-like specimens under hydro-mechanical coupling. Chinese Journal of Rock Mechanics and Engineering. 2019;38(S1):2823–30.
* View Article
* Google Scholar
7. 7. Wang L, Liu HQ, Li SB, Chen LP, Liu HQ. Fracture evolution characteristics of prefabricated crack coal under different confining pressures. Journal of Mining & Safety Engineering. 2023;40(04):786–97.
* View Article
* Google Scholar
8. 8. Li Z, Tian H, Niu Y, Wang EY, Zhang X, He S, et al. Study on the acoustic and thermal response characteristics of coal samples with various prefabricated crack angles during loaded failure under uniaxial compression. Journal of Applied Geophysics. 2022;200:104618.
* View Article
* Google Scholar
9. 9. Wang C, Wang L, Liu JF, Liu J. Mechanical properties and failure characteristics of granite intersected with single fractures under uniaxial compression. Journal of Southwest Jiaotong University. n.d.;59(02):369–76.
* View Article
* Google Scholar
10. 10. Li Y, Liu J, Yu Q. Patterns of influence of parallel rock fractures on the mechanical properties of the rock–coal combined body. Sustainability. 2022;14(20):13555.
* View Article
* Google Scholar
11. 11. Ma S, Liu K, Guo T, Yang J, Li X, Yan Z. Experimental and numerical investigation on the mechanical characteristics and failure mechanism of cracked coal & rock-like combined sample under uniaxial compression. Theoretical and Applied Fracture Mechanics. 2022;122:103583.
* View Article
* Google Scholar
12. 12. Chen GB, Tang W, Chen SJ, Wang EY, Wang CY, Li T, et al. Damage effect and deterioration mechanism of mechanical properties of fractured coal–rock combined body under water-rock interaction. Rock Mechanics and Rock Engineering. 20241–25.
* View Article
* Google Scholar
13. 13. Zhang D, Liu CX, Zhou Y, Ma DP. Characteristics of fracture mechanism and acoustic emission of rock-coal combined body with prefabricated fissure. Geotechnical and Geological Engineering. 2020;38:6245–54.
* View Article
* Google Scholar
14. 14. Li CJ, Xu Y, Chen PY, Li HL, Lou PJ. Dynamic mechanical properties and fragment fractal characteristics of fractured coal-rock-like combined bodies in split Hopkinson pressure bar tests. Natural Resources Research. 2020; 29: 3179–3195.
* View Article
* Google Scholar
15. 15. Lv H, Wang D, Cheng Z, Zhang Y, Zhou T. Study on mechanical characteristics and failure modes of coal-mudstone combined body with prefabricated crack. Mathematics. 2022;10(2):177.
* View Article
* Google Scholar
16. 16. Wang H, Li Y, Cao SG, Pan KR, Yang HY, Zhang KW. Brazilian splitting test study on crack propagation process and macroscopic failure mode of pre-cracked black shale. Chinese Journal of Rock Mechanics and Engineering. 2020;39(05):912–26.
* View Article
* Google Scholar
17. 17. Li X, Bai Y, Chen X, Zhao X, Lv M. Experimental and numerical study on crack propagation and coalescence in rock-like materials under compression. The Journal of Strain Analysis for Engineering Design. 2021;56(8):548–62.
* View Article
* Google Scholar
18. 18. Zhao YL, Zhang LY, Wang WJ, Pu CZ, Wan W, Tang JZ. Cracking and stress–strain behavior of rock-like material containing two flaws under uniaxial compression. Rock Mechanics and Rock Engineering. 2016;49:2665–87.
* View Article
* Google Scholar
19. 19. Zeng J, Pu C, Wang Q, Shen Q, Zeng Q, Yang Z. Tunnel and underground engineering rock mass water inrush damage and acoustic emission characteristics. PLoS One. 2024;19(9):e0307700. pmid:39303002
* View Article
* PubMed/NCBI
* Google Scholar
20. 20. Wang K, Pan H, Fujii Y. Study on energy distribution and attenuation of CO2 fracturing vibration from coal-like material in a new test platform. Fuel. 2024;356:129584.
* View Article
* Google Scholar
21. 21. Wang K, Chang C. Study on the characteristics of CO2 fracturing rock damage based on fractal theory. Theoretical and Applied Fracture Mechanics. 2024;134:104691.
* View Article
* Google Scholar
22. 22. Li Z, Fujii Y, Alam A, Li ZH, Du F, Wei WJ. Implementing a simple 2D constitutive model for rocks into finite element method. Computers and Geotechnics. 2024;167:106095.
* View Article
* Google Scholar
23. 23. Pu CZ, Yang SJ, Zhang CY. Fracture mechanism of pre-cracked specimens influenced by opening width. Chinese Journal of Geotechnical Engineering. 2019;41(10):1836–1844.
* View Article
* Google Scholar
24. 24. Tang SB, Sun K, Zhang YL, Wu F. Fracture growth of Brazilian disk with open/closed cracks. Journal of China Coal Society. 2021;46(11):3459–69.
* View Article
* Google Scholar
25. 25. Wang L, Yuan QP, Xie GX, Gu SH, Jiao ZH, Liu HQ, et al. Length-diameter ratio effect of energy dissipation and fractals of coal samples under impact loading. Journal of China Coal Society. 2022; 47(04):1534–1546.
* View Article
* Google Scholar
26. 26. Yuan C. Research on Mechanical Properties and Strength Prediction of Rock Mass with Different Fracture Geometry Characteristics. Xi`an University of Science and Technology; 2021.
27. 27. Zhu QQ, Li DY. Experimental investigation on crack behavior and stress thresholds of sandstone containing a square inclusion under uniaxial compression. Applied Sciences. 2020;10(21):7621.
* View Article
* Google Scholar
28. 28. Cai YB, Wang K, Xu C. Comparative experimental study on deformation and damage characteristics of single coal rock and primary coal-rock combination. Journal of Mining Science and Technology. 2020;5(03):278–83.
* View Article
* Google Scholar
29. 29. Kulatilake P, Malama B, Wang JL. Physical and particle flow modeling of jointed rock block behavior under uniaxial loading. International Journal of Rock Mechanics and Mining Sciences. 2001;38(5):641–57.
* View Article
* Google Scholar
30. 30. Deisman N, Mas Ivars D, Pierce M. PFC2D smooth joint contact model numerical experiments. Proceedings of the GeoEdmonton. 2008;8.
* View Article
* Google Scholar
31. 31. Zhen Z, Yang SQ, Chen CP, Tian WL, Qian JH, Li X. Simulation study on mechanical behavior of bedding yellow sandstone under triaxial cycling loading-unloading. Journal of Basic Science and Engineering. 2023;31(03):731–40.
* View Article
* Google Scholar
32. 32. Li W, Si L, Lu Z, Yi K, Wu L, . Determination of coal cracking initiation strength under uniaxial compression and analysis of its key factors. Journal of China Coal Society. n.d.;46(S2):670–80.
* View Article
* Google Scholar
33. 33. Liu X, Deng Z, Liu Y, Lu Y, Liu S, Han Y. Macroscopic and microscopic analysis of particle flow in pre-peak cyclic direct shear test of rock joint. Journal of China Coal Society. 2019;44(07):2103–15.
* View Article
* Google Scholar
34. 34. Wang XQ. Jointed coal mass model reconstruction and its mechanical response characteristics. China University of Mining & Technology, Beijing. 2017.
35. 35. Shi C, Zhang Q, Wang SN. Numerical simulation technology and application of particle flow. Beijing: China Architecture & Building Press; 2018.
36. 36. Yang L, Gao FQ, Wang XQ. Mechanical response and energy partition evolution of coal-rock combinations with different strength ratios. Chinese Journal of Rock Mechanics and Engineering. 2020;39(S2):3297–305.
* View Article
* Google Scholar
37. 37. Yang L, Gao FQ, Wang XQ, Li JZ. Energy evolution law and failure mechanism of coal-rock combined specimen. Journal of China Coal Society. 2019;44(12):3894–902.
* View Article
* Google Scholar
38. 38. Liu X, Li Y, Wang W, Zhou Y, Cheng L, Fan Z. Research on mechanical properties and strength criterion of carbonaceous shale with pre-existing fissures under drying-wetting cycles. Chinese Journal of Rock Mechanics and Engineering. 2022;41(02):228–39.
* View Article
* Google Scholar
39. 39. Luo D, Lu SH, Su GS, Tao HH. Experimental study on rock burst of granite with prefabricated single crack under true-triaxial stress condition with a free face. Rock and Soil Mechanics. 2023;44(01):75–87.
* View Article
* Google Scholar
40. 40. Zhang W, Zhou GQ, Zhang HB, Zhang Y. Experimental research on the influence of obliquity on the mechanical characteristics of a fractured rock mass. Journal of China University of Mining & Technology. 2009;38(01):30–3.
* View Article
* Google Scholar
41. 41. Wen L, Feng WJ, Li MY, Kou ZL, Wang L, Yu JH. Strain rate effect on crack propagation and fragmentation characteristics of red sandstone containing pre-cracks. Explosion and Shock Waves. 2023;43(11):88–107.
* View Article
* Google Scholar
42. 42. Zhou X, Long Y, Ye W. Experimental investigations on the cracking and mechanical responses of PMMA samples with two 3D embedded elliptic flaws under uniaxial compression. Geohazard Mechanics. 2023;1(1):77–85.
* View Article
* Google Scholar
43. 43. Qi J, Zhou L, Zhang HD, Chen JX, Ma LJ, Shi TY. Research on crack evolution law and mechanical analysis of three cracked rock masses subjected to compression load. Theoretical and Applied Fracture Mechanics. 2023;127:104035.
* View Article
* Google Scholar
44. 44. Li Q, Cai ST, Li HJ, Zhong YQ, Liu YW. Particle flow simulation of fracture characteristics of rock-concrete combination with single crack. Chinese Journal of High Pressure Physics. 2024;38(05):190–205.
* View Article
* Google Scholar
45. 45. Xu C, Wang W, Wang K, Hu K, Cao Z, Zhang Y. Influence of coal-rock interface inclination on the damage and failure law of original coal-rock combination. Engineering Failure Analysis. 2024;161:108275.
* View Article
* Google Scholar
46. 46. Li Q, Nie F, Pan C, Li L, Zhong Y, Yu M, et al. Energy dissipation damage constitutive relation of CFRP passively confined coal sample. Heliyon. 2024;10(18):e37586. pmid:39318801
* View Article
* PubMed/NCBI
* Google Scholar
47. 47. Huang D, Qiao S, Chang X, Wang X, Lu H, Pan X. Study on macro–micro mechanical behavior of rock like samples with hole and cracks. Comp Part Mech. 2023;11(4):1579–98.
* View Article
* Google Scholar
48. 48. Ge XR, Ren JX, Pu YB. Macro-Microscopic damage mechanics of rock and earth mass. Beijing: Science Press; 2004.
49. 49. Zhang L, Ren T, Li XC, Tan LH. Acoustic emission, damage and cracking evolution of intact coal under compressive loads: Experimental and discrete element modelling. Engineering Fracture Mechanics. 2021;252:107690.
* View Article
* Google Scholar
50. 50. Yi M, Zhao T, Ma FF, Wang L, Yuan C. Study of constitutive model of frozen sandstone damage based on Weibull distribution. Coal Geology & Exploration. n.d.;50(8):116–24.
* View Article
* Google Scholar
51. 51. Li Q, Gao SL, Hu LL, Yu MM, Liu YW, Zeng XG, et al. Constitutive relation of energy dissipation damage of heterogeneous coal samples under different loading rates. Journal of China Coal Society. 2022;47(S1):90–102.
* View Article
* Google Scholar
52. 52. Zhang L, Qi QX, REN T, Li XC, Gao K, Li CY. Study on the damage and fracture characteristics of coal rock based on the X-ray micro-CT scanning technology and statistical strength theory. Coal Science and Technology. 2023;51(S2):1–12.
* View Article
* Google Scholar
Citation: Li Q, Pan C, Zhong Y, Li W, Li L, Nie F, et al. (2025) Study on the damage constitutive characteristics of coal-rock composites under uniaxial compression: Influence of prefabricated crack angle and geometric dimensions. PLoS ONE 20(3): e0316586. https://doi.org/10.1371/journal.pone.0316586
About the Authors:
Qingwen Li
Roles: Conceptualization, Methodology, Writing – original draft, Writing – review & editing
E-mail: [email protected]
Affiliation: School of Civil and Architectural Engineering, Liaoning University of Technology, Jinzhou, China
ORICD: https://orcid.org/0009-0004-5467-362X
Chuangchuang Pan
Roles: Conceptualization, Data curation, Methodology, Writing – original draft, Writing – review & editing
Affiliation: School of Civil and Architectural Engineering, Liaoning University of Technology, Jinzhou, China
ORICD: https://orcid.org/0009-0000-1533-014X
Yuqi Zhong
Roles: Data curation
Affiliation: School of Civil and Architectural Engineering, Liaoning University of Technology, Jinzhou, China
Wenxia Li
Roles: Formal analysis
Affiliation: School of Civil and Architectural Engineering, Liaoning University of Technology, Jinzhou, China
Ling Li
Roles: Validation
Affiliation: School of Civil and Architectural Engineering, Liaoning University of Technology, Jinzhou, China
Fanfan Nie
Roles: Formal analysis
Affiliation: School of Civil and Architectural Engineering, Liaoning University of Technology, Jinzhou, China
Jiabo Chen
Roles: Investigation, Validation
Affiliation: School of Civil and Architectural Engineering, Liaoning University of Technology, Jinzhou, China
[/RAW_REF_TEXT]
[/RAW_REF_TEXT]
[/RAW_REF_TEXT]
[/RAW_REF_TEXT]
1. Xie HP, Zhang ZT, Gao F, Zhang R, Gao MZ, Liu JF. Stress-fracture-seepage field behavior of coal under different mining layouts. Journal of China Coal Society. 2016;41(10):2405–17.
2. Xie H, Zhou H, Xue D, Gao F. Theory, technology and engineering of simultaneous exploitation of coal and gas in China. Journal of China Coal Society. 2014;39(08):1391–7.
3. Wang W, Chen YP, Yuan L, Chen RZ, Wang HF, Du K. Floor fracture evolution and relief gas drainage timeliness in deeper underground short-distance upper protective coal seam extraction. Journal of China Coal Society. 2016;41(01):138–48.
4. Ding K, Gu QH, Wang LG, Guo JX. Study on the coal rib burst mechanism in deep roadway and energy release. Journal of China Coal Society. 2024;49(09):3759–80.
5. Chen GB, Tang W, Li T, Wang CY, Wang EY, Zhang GH. Mechanical response and instability model of fractured coal-rock combined body. Rock and Soil Mechanics. 2024;45(09):2633–52.
6. Zhao C, Xin JQ, Niu JL, Ma CC. Experimental study on crack propagation of precrack rock-like specimens under hydro-mechanical coupling. Chinese Journal of Rock Mechanics and Engineering. 2019;38(S1):2823–30.
7. Wang L, Liu HQ, Li SB, Chen LP, Liu HQ. Fracture evolution characteristics of prefabricated crack coal under different confining pressures. Journal of Mining & Safety Engineering. 2023;40(04):786–97.
8. Li Z, Tian H, Niu Y, Wang EY, Zhang X, He S, et al. Study on the acoustic and thermal response characteristics of coal samples with various prefabricated crack angles during loaded failure under uniaxial compression. Journal of Applied Geophysics. 2022;200:104618.
9. Wang C, Wang L, Liu JF, Liu J. Mechanical properties and failure characteristics of granite intersected with single fractures under uniaxial compression. Journal of Southwest Jiaotong University. n.d.;59(02):369–76.
10. Li Y, Liu J, Yu Q. Patterns of influence of parallel rock fractures on the mechanical properties of the rock–coal combined body. Sustainability. 2022;14(20):13555.
11. Ma S, Liu K, Guo T, Yang J, Li X, Yan Z. Experimental and numerical investigation on the mechanical characteristics and failure mechanism of cracked coal & rock-like combined sample under uniaxial compression. Theoretical and Applied Fracture Mechanics. 2022;122:103583.
12. Chen GB, Tang W, Chen SJ, Wang EY, Wang CY, Li T, et al. Damage effect and deterioration mechanism of mechanical properties of fractured coal–rock combined body under water-rock interaction. Rock Mechanics and Rock Engineering. 20241–25.
13. Zhang D, Liu CX, Zhou Y, Ma DP. Characteristics of fracture mechanism and acoustic emission of rock-coal combined body with prefabricated fissure. Geotechnical and Geological Engineering. 2020;38:6245–54.
14. Li CJ, Xu Y, Chen PY, Li HL, Lou PJ. Dynamic mechanical properties and fragment fractal characteristics of fractured coal-rock-like combined bodies in split Hopkinson pressure bar tests. Natural Resources Research. 2020; 29: 3179–3195.
15. Lv H, Wang D, Cheng Z, Zhang Y, Zhou T. Study on mechanical characteristics and failure modes of coal-mudstone combined body with prefabricated crack. Mathematics. 2022;10(2):177.
16. Wang H, Li Y, Cao SG, Pan KR, Yang HY, Zhang KW. Brazilian splitting test study on crack propagation process and macroscopic failure mode of pre-cracked black shale. Chinese Journal of Rock Mechanics and Engineering. 2020;39(05):912–26.
17. Li X, Bai Y, Chen X, Zhao X, Lv M. Experimental and numerical study on crack propagation and coalescence in rock-like materials under compression. The Journal of Strain Analysis for Engineering Design. 2021;56(8):548–62.
18. Zhao YL, Zhang LY, Wang WJ, Pu CZ, Wan W, Tang JZ. Cracking and stress–strain behavior of rock-like material containing two flaws under uniaxial compression. Rock Mechanics and Rock Engineering. 2016;49:2665–87.
19. Zeng J, Pu C, Wang Q, Shen Q, Zeng Q, Yang Z. Tunnel and underground engineering rock mass water inrush damage and acoustic emission characteristics. PLoS One. 2024;19(9):e0307700. pmid:39303002
20. Wang K, Pan H, Fujii Y. Study on energy distribution and attenuation of CO2 fracturing vibration from coal-like material in a new test platform. Fuel. 2024;356:129584.
21. Wang K, Chang C. Study on the characteristics of CO2 fracturing rock damage based on fractal theory. Theoretical and Applied Fracture Mechanics. 2024;134:104691.
22. Li Z, Fujii Y, Alam A, Li ZH, Du F, Wei WJ. Implementing a simple 2D constitutive model for rocks into finite element method. Computers and Geotechnics. 2024;167:106095.
23. Pu CZ, Yang SJ, Zhang CY. Fracture mechanism of pre-cracked specimens influenced by opening width. Chinese Journal of Geotechnical Engineering. 2019;41(10):1836–1844.
24. Tang SB, Sun K, Zhang YL, Wu F. Fracture growth of Brazilian disk with open/closed cracks. Journal of China Coal Society. 2021;46(11):3459–69.
25. Wang L, Yuan QP, Xie GX, Gu SH, Jiao ZH, Liu HQ, et al. Length-diameter ratio effect of energy dissipation and fractals of coal samples under impact loading. Journal of China Coal Society. 2022; 47(04):1534–1546.
26. Yuan C. Research on Mechanical Properties and Strength Prediction of Rock Mass with Different Fracture Geometry Characteristics. Xi`an University of Science and Technology; 2021.
27. Zhu QQ, Li DY. Experimental investigation on crack behavior and stress thresholds of sandstone containing a square inclusion under uniaxial compression. Applied Sciences. 2020;10(21):7621.
28. Cai YB, Wang K, Xu C. Comparative experimental study on deformation and damage characteristics of single coal rock and primary coal-rock combination. Journal of Mining Science and Technology. 2020;5(03):278–83.
29. Kulatilake P, Malama B, Wang JL. Physical and particle flow modeling of jointed rock block behavior under uniaxial loading. International Journal of Rock Mechanics and Mining Sciences. 2001;38(5):641–57.
30. Deisman N, Mas Ivars D, Pierce M. PFC2D smooth joint contact model numerical experiments. Proceedings of the GeoEdmonton. 2008;8.
31. Zhen Z, Yang SQ, Chen CP, Tian WL, Qian JH, Li X. Simulation study on mechanical behavior of bedding yellow sandstone under triaxial cycling loading-unloading. Journal of Basic Science and Engineering. 2023;31(03):731–40.
32. Li W, Si L, Lu Z, Yi K, Wu L, . Determination of coal cracking initiation strength under uniaxial compression and analysis of its key factors. Journal of China Coal Society. n.d.;46(S2):670–80.
33. Liu X, Deng Z, Liu Y, Lu Y, Liu S, Han Y. Macroscopic and microscopic analysis of particle flow in pre-peak cyclic direct shear test of rock joint. Journal of China Coal Society. 2019;44(07):2103–15.
34. Wang XQ. Jointed coal mass model reconstruction and its mechanical response characteristics. China University of Mining & Technology, Beijing. 2017.
35. Shi C, Zhang Q, Wang SN. Numerical simulation technology and application of particle flow. Beijing: China Architecture & Building Press; 2018.
36. Yang L, Gao FQ, Wang XQ. Mechanical response and energy partition evolution of coal-rock combinations with different strength ratios. Chinese Journal of Rock Mechanics and Engineering. 2020;39(S2):3297–305.
37. Yang L, Gao FQ, Wang XQ, Li JZ. Energy evolution law and failure mechanism of coal-rock combined specimen. Journal of China Coal Society. 2019;44(12):3894–902.
38. Liu X, Li Y, Wang W, Zhou Y, Cheng L, Fan Z. Research on mechanical properties and strength criterion of carbonaceous shale with pre-existing fissures under drying-wetting cycles. Chinese Journal of Rock Mechanics and Engineering. 2022;41(02):228–39.
39. Luo D, Lu SH, Su GS, Tao HH. Experimental study on rock burst of granite with prefabricated single crack under true-triaxial stress condition with a free face. Rock and Soil Mechanics. 2023;44(01):75–87.
40. Zhang W, Zhou GQ, Zhang HB, Zhang Y. Experimental research on the influence of obliquity on the mechanical characteristics of a fractured rock mass. Journal of China University of Mining & Technology. 2009;38(01):30–3.
41. Wen L, Feng WJ, Li MY, Kou ZL, Wang L, Yu JH. Strain rate effect on crack propagation and fragmentation characteristics of red sandstone containing pre-cracks. Explosion and Shock Waves. 2023;43(11):88–107.
42. Zhou X, Long Y, Ye W. Experimental investigations on the cracking and mechanical responses of PMMA samples with two 3D embedded elliptic flaws under uniaxial compression. Geohazard Mechanics. 2023;1(1):77–85.
43. Qi J, Zhou L, Zhang HD, Chen JX, Ma LJ, Shi TY. Research on crack evolution law and mechanical analysis of three cracked rock masses subjected to compression load. Theoretical and Applied Fracture Mechanics. 2023;127:104035.
44. Li Q, Cai ST, Li HJ, Zhong YQ, Liu YW. Particle flow simulation of fracture characteristics of rock-concrete combination with single crack. Chinese Journal of High Pressure Physics. 2024;38(05):190–205.
45. Xu C, Wang W, Wang K, Hu K, Cao Z, Zhang Y. Influence of coal-rock interface inclination on the damage and failure law of original coal-rock combination. Engineering Failure Analysis. 2024;161:108275.
46. Li Q, Nie F, Pan C, Li L, Zhong Y, Yu M, et al. Energy dissipation damage constitutive relation of CFRP passively confined coal sample. Heliyon. 2024;10(18):e37586. pmid:39318801
47. Huang D, Qiao S, Chang X, Wang X, Lu H, Pan X. Study on macro–micro mechanical behavior of rock like samples with hole and cracks. Comp Part Mech. 2023;11(4):1579–98.
48. Ge XR, Ren JX, Pu YB. Macro-Microscopic damage mechanics of rock and earth mass. Beijing: Science Press; 2004.
49. Zhang L, Ren T, Li XC, Tan LH. Acoustic emission, damage and cracking evolution of intact coal under compressive loads: Experimental and discrete element modelling. Engineering Fracture Mechanics. 2021;252:107690.
50. Yi M, Zhao T, Ma FF, Wang L, Yuan C. Study of constitutive model of frozen sandstone damage based on Weibull distribution. Coal Geology & Exploration. n.d.;50(8):116–24.
51. Li Q, Gao SL, Hu LL, Yu MM, Liu YW, Zeng XG, et al. Constitutive relation of energy dissipation damage of heterogeneous coal samples under different loading rates. Journal of China Coal Society. 2022;47(S1):90–102.
52. Zhang L, Qi QX, REN T, Li XC, Gao K, Li CY. Study on the damage and fracture characteristics of coal rock based on the X-ray micro-CT scanning technology and statistical strength theory. Coal Science and Technology. 2023;51(S2):1–12.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2025 Li et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
In coal mining environments with complex topographic and geological conditions, the presence of primary cracks in the rock strata of the upper mining airspace is critical to mine safety, especially when roof slabs are broken. Cracked roof slabs not only increase risks during mining but also make mining operations more challenging. Therefore, studying the initial damage state of the roof rock formation is great significance. In this study, the effects of different prefabricated crack sizes and inclination angles on the mechanical properties of the coal-rock composite containing cracks were analyzed through the uniaxial compression experiments and PFC2D numerical simulations. The results show that the peak strength and elastic modulus of the coal-rock composites fall between those of pure coal and sandstone, while the macroscopic mechanical parameters of crack-containing composites are significantly lower than those of non-crack-containing composites. Coal-rock composites with different crack characteristics exhibited different mechanical properties, with their damage modes were caused by the combined effects of tensile and shear damage. The increase in crack inclination altered the crack extension path, and the final damage of the specimen manifested first in the upper part, then the middle part, and ultimately in the lower part of the coal body, with tension-induced bulk damage being the dominant failure mode. Analysis of the radial cumulative map revealed that cracks primarily extended along 90° and 270° directions, indicating a strong tendency for crack propagation under axial pressure. The damage evolution curves indicate a nonlinear relationship between the damage factor and strain. While increased crack inclination enhances the compressive performance of coal-rock composites, it simultaneously accelerates structural destabilization. These findings offer theoretical insights into the damage mechanisms of coal-rock composites with cracks, serving as valuable references for coal mining safety.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details
; Pan, Chuangchuang
; Zhong, Yuqi; Li, Wenxia; Li, Ling; Nie, Fanfan; Chen, Jiabo