Rock Mass Engineering Research Articles

Probabilistic evaluation method for the stability of large underground cavern considering the uncertainty of rock mass mechanical parameters: A case study of Baihetan underground powerhouse project

For deep large underground cavern projects under complex geological conditions, the determination of rock mass mechanical parameters and stability analysis is fraught with a great deal of uncertainty. This paper proposes a set of methods for estimating rock mass mechanical parameters and probabilistic stability assessment suitable for construction characteristics of large underground caverns. On the one hand, according to the Hoek-Brown criterion and Monte Carlo simulation, a dynamic estimation method for rock mass mechanical parameters is proposed by using Rock Mass Rating (RMR), the uniaxial compressive strength (UCS) of intact rock, and the material constant (mb) as input parameters. On the other hand, considering the uncertainty of rock mass mechanical parameters, a probabilistic evaluation method applicable to the unloading response characteristics of rock mass (including displacement and fracturing zones) around the large cavern is put forward, combining the point estimation method and numerical simulation. The proposed methods were applied to an underground powerhouse, with an excavation span of 34 m and height of 88.7 m, at the Baihetan hydropower station in the southwest of China. Based on extensive field investigations and tests, the probability distributions of key rock mass mechanical parameters were obtained. Furthermore, the high sensitivity of input parameters in the H-B criterion for estimating rock mass mechanical parameters were revealed. Through dynamic simulation, the probability distributions of the displacement and fracturing zones in surrounding rock during the layer-by-layer excavation of the cavern were presented. Comprehensive on-site tests showed that the simulated results were in basic agreement with the actual unloading response behavior during the cavern excavation, validating the correctness and applicability of the proposed methods. The study provides a relatively comprehensive approach for estimating rock mass mechanical parameters and stability evaluation in similar underground engineering projects, and also has guiding significance for predicting and preventing engineering rock mass hazards.

Scale effect of rock discontinuity considering all morphological information

Having an accurate understanding of the scale effect of surface morphology characteristics is crucial to examining the mechanical behavior of rock structural plane. At present, the quantification and sampling methods of surface morphology show diversity, which is the potential reason for the inconsistent research conclusions on scale effect. Firstly, based on mathematical statistics and correlation analysis, the most representative parameter is proposed from hundreds of morphological parameters. Then, the previous scale effect sampling methods are analyzed. In order to ensure that the selected samples are representative, a novel sampling method, considering all morphological information, is proposed. By means of the novel quantification and sampling methods, the size effect characteristics are systematically analyzed. Under the conditions of different rock types, shear directions and sampling locations, etc., discontinuity roughness does not change significantly with sampling scale. As sampling scale increases, the distribution range of representative samples is gradually concentrated, the total amount decreases, and the proportion increases. However, the distribution of representative samples on the initial structural plane does not show obvious regularity. These findings would provide theoretical support for the deformation control and stability analysis of rock mass in engineering.

Rock strength prediction based on machine learning: A study from prediction model to mechanism explanation

Rational and effective assessment of rock strength, especially uniaxial compressive strength (UCS), is the crucial basis of rock engineering practices. As computer technology advances, the application of machine learning techniques in predicting rock strength has become increasingly prevalent. However, due to the black-box effect of machine learning algorithms, the prediction process is always indescribable. Therefore, an innovative approach for the interpretable predicting uniaxial compressive strength (UCS) is initially introduced by using Explainable AI to crack the black-box effect of machine learning algorithms and explain the prediction mechanism. The findings reveal that the Xgboost model, optimized with interpretable hyper-parameters, achieves notable success in estimating the uniaxial compressive strength, with an RMSE of 14.69 and an R2 that can reach 0.9117. Additionally, multi-angle comparisons have verified that the robustness of the AOA-Xgboost model surpasses that of the other four hybrid models. The model employs SHAP analysis to measure how much each input factor contributes to the outcomes of the predictions. This work has potential for the prediction of crucial parameters and analysis of internal mechanisms of rock mass in geotechnical engineering.

Research on triaxial compressive mechanical properties and damage constructive model of post-thawing double-fractured quasi-sandstones with different angles

Joints and fractures lead to different failure mechanisms in rock masses under different environments. The mechanical properties and failure mechanisms of rocks with fissures are key problems in rock mass engineering. Parallel double-fracture quasi-sandstone specimens with different dip angles were prepared and subjected to triaxial compression tests after a single freeze-thaw cycle. Pore development, crack propagation, damage evolution, and failure characteristics were analysed. Combined with the strength distribution theory of microelements and the static elastic modulus theory, a damage constitutive model of double-fracture quasi-sandstone under freeze-thaw cycles and loads was established. This study explored the pore development, fracture propagation, damage evolution, and failure characteristics of fractured sandstone after thawing. The results showed that the compression wave velocity of the thawed specimens decreased, the nuclear magnetic resonance (NMR) T2 curve shifted to the right, and the frost heave force promoted the development of the internal porosity in the specimens. With an increase in the crack dip angle, peak stress, expansion stress, cohesion and internal friction angle, the specimen showed a ‘U’ shaped change trend, compression cracks, and rock bridge penetration rate after failure decreased, and mixed failure of tension and shear gradually changed into shear failure. When the dip angles were 30° and 60°, the double fractured quasi-sandstone had larger total damage and more obvious brittle failure characteristics.

Drilling Process Monitoring for Predicting Mechanical Properties of Jointed Rock Mass: A Review

Reliably assessing the quality and mechanical properties of rock masses is crucial in underground engineering. However, existing methods have significant limitations in terms of applicability and accuracy. Therefore, a field measurement method that meets the real-time monitoring and safety requirements for the quality of engineering rock masses is needed. Firstly, the research findings of domestic and international scholars on the application of drilling process monitoring technology are comprehensively analyzed. Rotary cutting penetration tests are conducted on tuff rock masses containing fractures and joints. Various rock mass classification and evaluation standards are integrated with rotary penetration tests. Rotary cutting penetration tests are used to determine the residual strength of rock, based on this review. The rationality of the calculated mi parameter values is validated. The peak strength, residual strength, and errors of the rock are obtained based on the penetration method. The rock quality index rock quality designation from drilling (RQDd) is redefined, based on the drilling process monitoring apparatus (DPMA). Rock mass classification is conducted, based on the correlation between the standard deviation of rotary drilling energy and the rock quality designation (RQD). Additionally, a new relational formula is introduced to determine the RQD from variations in drilling energy, based on discontinuity frequency. This field measurement method undoubtedly provides a crucial scientific basis for rock design and construction, ensuring long-term safety in engineering applications.

Mechanical responses of chemically corroded sandstone under cyclic disturbance: Insights from fatigue properties and macro-micro fracturing mechanism

Engineering rock masses are likely to be subjected to coupled chemical corrosion and cyclic load disturbance, deteriorating the internal structure of rocks and further threatening the long-term stability of engineering structures. In this study, to investigate the combined effect of cyclic loading and chemical corrosion, the cyclic uniaxial compression tests are conducted on sandstone specimens chemically corroded by different corrosion times (0, 20, 40, 60 days) and pH values (1, 3, 7, 11, 13). Our results indicate that chemically corroded sandstone specimens featured by an increase in irreversible strain and hysteresis energy density compared with natural rocks. A coupling damage variable is proposed to reflect the initial damage induced by chemical corrosion and the mechanical damage induced by cyclic disturbance, and the accumulation rate of coupled chemical-fatigue damage variables significantly increases with increasing corrosion time or enhancing acidity and alkalinity. Fatigue failure modes of tested specimens transform from splitting failure to unidirectional shear failure as corrosion time increases, and chemically corroded specimens are featured by ductile failure characteristics compared to natural state. Moreover, X-ray diffraction (XRD) and Scanning Electron Microscope (SEM) are employed to analyze the micro-fracturing mechanism of sandstone. Cement minerals dominate the reactions in acidic solutions, whereas quartz takes the lead in alkaline conditions. Due to pronounced influence of cement minerals on the macroscopic mechanical properties of sandstone, acid corrosion specimens tend to experience more severe damage. Micro-cracks inside rocks under cyclic loading have sufficient time to pass through weakly bonded particles, and chemical reactions markedly weaken inter-particle bonding degree, and thus the fatigue fracture of chemically corroded specimens predominantly exhibits intergranular (IG) fracture.

Experimental study on the bearing mechanical characteristics and deformation failure mechanism of damaged sandy mudstone

AbstractStudying the mechanical properties and failure mechanisms of damaged rocks holds practical significance in guiding tunnel or roadway support and assessing the stability of engineering rock masses. To uncover the failure mechanisms of sandy mudstone with varying degrees of damage, we conducted cyclic loading‐unloading tests to create specimens at different damage levels. Subsequently, uniaxial loading tests were performed to investigate the fundamental mechanical properties, wave velocity, and acoustic emission (AE) parameter evolution characteristics of the damaged specimens. The classification of damaged rocks was based on the energy evolution characteristics of damaged sandy mudstone. Finally, we discuss the damage and deterioration mechanisms of different types of sandy mudstone. The findings contribute to establishing a theoretical foundation for evaluating the degree of rock damage.

Investigation on uniaxial compression and fracture damage mode of prefabricated parallel double-jointed red sandstone.

As a special type of joint fracture, the fracture evolution characteristics of parallel double joints have important engineering significance for the stability analysis of fractured rock mass. In this work, a new method for calculating stress intensity factor of parallel double-jointed fractures was importantly proposed. Physical uniaxial compression tests were carried out on parallel double jointed red sandstone filled with cement mortar under different geometric parameters, and the macroscopic mechanical properties and failure characteristics of red sandstone are deeply analyzed. The results show that the larger the connectivity rate is, the smaller the peak stress and strain are. The increase of connectivity rate will affect the change rate of transverse strain in the center of rock bridge. The closer the dip angle of the joint is, the lower the peak stress is and the shorter the failure time is. The damage mode of joint tip encroachment affects the lateral displacement of the rock bridge center, and the displacement is always close to the first damage section. The closer the joint tip is to the load, the easier the end-face penetrating cracks occur. The research content can provide basic support for guaranteeing the stability of underground engineering rock mass.

A complex-network-based estimation of the representative elementary volume and equivalent permeability coefficient for fractured rock masses

The determination of the representative elementary volume (REV) and equivalent permeability coefficient (EPC) of fractured rock masses is crucial for studying the hydraulic behaviors of engineering rock masses. The complex network originated from graph theory and represents a rapidly developing interdisciplinary field in the 21st century. In this paper, we propose a novel method that utilizes complex-network-based metrics to estimate the REV and EPC of fractured rock masses. In this method, fractures and their intersections are treated as nodes and edges of a complex network, respectively. We propose a transformation procedure of a discrete fracture network (DFN) into a complex network, and the corresponding algorithms are conducted in MATLAB. Subsequently, we obtain a complex network, and several complex-network-based metrics, such as average degree, average clustering coefficient, and graph density, are calculated using Gephi software. We employ the average degree to estimate the REV size of DFN models, while the coefficient of variation is used to quantitatively assess changes in the average degree of multiple random models. Furthermore, we establish an empirical formula utilizing the average degree, average clustering coefficient, and graph density to estimate the EPC of a large-scale DFN model based on the EPC of a small-scale DFN model. To verify the developed method, we conduct numerical experiments using the MAFIC (Matrix/Fracture Interaction Code) solver within the FracMan software. This enabled us to calculate the EPC, which was then used to determine the REV. The results demonstrate that: (a) the REV size estimated by the average degree aligns with the REV size obtained through numerical simulations for both isotropic and anisotropic DFN models, and (b) the error rates between the average EPC calculated by the empirical formula and that obtained by numerical simulations are within 10%. Thus, the proposed seepage parameter estimation methods are deemed valid.

Context-aware semantic segmentation network for tunnel face feature identification

The automated interpretation of tunnel face geological information is significant to the construction decision-making of rock mass engineering. An intelligent recognition algorithm, named the Transformer and Convolution neural networks Semantic segmentation Network (TCSeNet), is introduced to overcome the low interpretation accuracy caused by certain limitations of existing automated interpretation methods. Firstly, a hybrid encoder is constructed to extract global and local context information using Transformer and Convolution Neural Networks, respectively. Furthermore, a context feature adaptive selection module is designed to cross-fusion the information. Secondly, a feature pyramid-like decoder is constructed for multi-scale object problems. Furthermore, a multi-scale feature refinement module is designed to improve the ability of multi-scale feature expression. A tunnel face dataset is constructed for training and testing. The research shows that the mIoU of TCSeNet reaches 92.57%. This indicates that the TCSeNet is suitable for automatic interpretation of tunnel face geological information.

Experimental and numerical study on the damage evolution and acoustic emission multi-parameter responses of single flaw sandstone under uniaxial compression

Cracks are generally the weakest position in engineering rock mass. It is of great practical significance to study the influence of cracks on the fracture process of rock mass and the precursor criterion of failure for disaster prevention and mitigation projects. In this paper, particle flow program (PFC2D) is used to analyze the mechanical strength and damage evolution of fractured rock. The multi-parameter stage characteristics of acoustic emission (AE) and the difference of fracture precursor information of fractured rock mass under uniaxial compression (UC) are explored. The results are as follows: (1) The strength and modulus of elasticity of fractured rock samples are smaller than that of intact rock mass, and increase in quadratic function with the increase of fracture angle, and the strength of B-90 sample is closest to intact rock samples; (2) The PFC analysis shows that the cracks in the defective rock samples all originate from the tip of the prefabricated crack, and the crack initiation angle is basically in the range of 60°∼90°. The loading process mainly produces tensile damage, and the number of tensile cracks are far more than that of shear cracks. Compressive stress and tensile stress increase with the increase of load, but the force chain of failure surface is the smallest after complete failure of rock; (3) The acoustic emission ringing number, RA/AF value, b value and entropy value of sandstone have obvious stage characteristics, the main frequency signal is mainly low frequency signal, and large-scale cracks are mainly produced in the failure process. The sudden increase point of AE ringing times and RA/AF value accumulation curve can be used as critical failure point. The fluctuation of B value and entropy value reflects the instability of fracture development, and “cliff-like” descending point and ascending point can be regarded as the precursor of critical failure of rock; (4) The failure precursor of acoustic emission parameters of sandstone appears in the peak load period of 85 %–98 %. The order of response time of different parameters is: cumulative RA/AF value, cumulative event number, entropy value, B value. The response time of B-45 and B-15 is earlier, while that of B-60 and B-30 is later. The research results are expected to deepen the understanding of the relevant standards of damage evolution and failure precursors of defective rocks in engineering applications.

Acoustic emission characteristics and fracture mechanism of sandstone in open-pit mines under different types of cyclic loads

Rock mass is one of the most important load-bearing media in geotechnical engineering. It has been continually vulnerable to geological tectonic movements, natural calamities, and human excavation activities. Its inherent weak surfaces such as primary pores, joints, and fissures have resulted in varying damage degrees. In mining operations, the damaged rock mass has a variety of negative impacts on the stability of its overlying structures and is frequently disturbed by the load. To study the damage law of rock mass under cyclic loading, in this paper, an acoustic emission (AE) device was employed to monitor the rock under the action of two types of cyclic loads: the variable upper and lower pre-loads, and the fixed upper and lower pre-loads. The damage of the loaded rock was split into three stages in this research, based on the features of the AE signals of the rock under uniaxial load, and the damage evolution of the loaded rock was analyzed in distinct stages. The AE signals of the rock under cyclic loading were mainly emitted in the first loading stage. When the stress did not exceed the maximum stress value in the stress history of the loaded rock, few new AE event was generated in the loaded rock. After the low-frequency cyclic static load, the AE signals varied with the load-bearing stress of the rock during the whole process from initial loading to failure, which was consistent with the characteristics of the AE signals of the loaded rock. The research results can be adapted to rock mass in open-pit mines stability analysis and risk prediction while providing some references for the early warning and danger relief of rock masses in engineering.

Sulfoaluminate cement-based double-liquid material grouting in fractured rock

Grouting is a very effective means to repair modified fractured rock. However, problems such as severe grout leakage, slow increase in stone body strength, uncontrolled solidification time and poor permeability are often encountered in the process of traditional grouting reinforcement of fractured surrounding rock. Sulfoaluminate cement-based double-liquid material has the characteristics of quick setting, early strength, micro-expansion and high strength, which shows good applicability in grouting reinforcement of fractured surrounding rock. Currently, there are few reports on this material. Therefore, this study mainly introduces the performance of sulfoaluminate cement-based grouting material, grouting technical process and engineering application case. The sulfoaluminate cement-based double-liquid material grouting technology proposed in this study will provide a new treatment method for the reinforcement of fractured rock masses in underground engineering.

Cross-scale analysis for the thermo-hydro-mechanical (THM) effects on the mechanical behaviors of fractured rock: Integrating mesostructure-based DEM modeling and machine learning

The excavation process induces the emergence of damage zones in the on-site engineering rock masses, significantly affecting the parameter evaluation of surrounding rock masses. Therefore, an in-depth understanding of the thermo-hydro-mechanical (THM) effects on the cross-scale failure process of fractured rock becomes progressively crucial in natural energy exploitation projects. In this study, a coupled thermo-hydro-mesostructure-based DEM (T-H-MSBM) model was developed to reconstruct rock microstructures and distinguish the THM responses of varying mineral grains, micro-defects and cracks in fractured granite. Five sets of numerical fractured granite in terms of different damage degrees were generated by restoring the form of excavation damage on site, and a series of compression simulations were conducted on the T-H-MSBMs under coupled temperature (25–250 °C), pore pressure (0–15 MPa) and confining pressure (15 MPa) conditions. Based on the cross-scale failure analysis on the fractured granite during the THM-compression loading process, the interplay between THM treatment and damage degree on the mechanical properties of fractured granite was revealed, and the main mechanisms affecting the varied macro mechanical properties were further discussed. Results indicate that both temperature and pore pressure exert the amplified deteriorating effect on the macro mechanical properties of fractured granite with increasing damage degree, while the temperature dependence becomes significantly more pronounced in the fractured granite with low damage degree. With an increase in the damage degree of fractured rock, the abundance of intra-mineral tensile cracks accelerates the initiation of pore pressure-induced cracks, and the coupled increase in pore pressure and temperature enhances the role of quartz-feldspar cracks in promoting the cracks in quartz. Finally, utilizing the optimized kernel ridge regression (KRR) model employed the grid search method, the correlation information database of rock macro-mechanical strength and THM factors were generated. The entire training time was less than 3 s, with the relative average errors of 0.7 MPa, which demonstrated the robust capability of machine learning method in predicting the coupled THM phenomena spanning multiple scales in geological materials and systems.

Study on energy damage evolution of multi-flaw sandstone with different flaw lengths

Flaws are common in rock mass engineering. The length and number of flaws have a great influence on the stability of rock mass. In this paper, uniaxial compression tests were carried out on six kinds of multi-flawed sandstone with different flaw lengths. Based on damage mechanics and energy theory, the mechanical variation characteristics, total strain energy, elastic strain energy and dissipation energy evolution law of multi-flawed sandstone with different flaw lengths under uniaxial action were analyzed. Based on the energy dissipation, the damage variable is redefined, and the sandstone damage process is divided into three stages: the initial damage weakening stage I, the damage stability stage II, and the damage accelerated development stage III. Using this conclusion, considering the initial compaction weakening stage of rock samples, Duncan model is introduced to modify the constitutive model of equivalent strain hypothesis. The results show that the theoretical curve is in good agreement with the experimental curve, indicating that the damage constitutive model can accurately reflect the damage and failure process of multi-flaw rock samples with different flaw lengths.

Experimental investigation on the influence of temperature on the fracture surface variations of granite after Brazilian splitting tests

Understanding the variation in the morphology of rock fracture surfaces after high-temperature is of great significance for evaluating the strength, deformation, and seepage characteristics of deep engineering rock masses. In this study, the variations in the granite fracture surface roughness after exposure to different temperatures (25–800 °C) were investigated, and the damage deterioration mechanism was revealed. The results showed that the maximum fluctuation (Z) and root mean square of the first derivative (Z2) of the granite fracture surface tended to increase with increasing temperature, whereas the fractal dimension (D) decreased, suggesting that the higher temperature contributed to greater roughness of the fracture surface. The anisotropy of the fracture surface was more evident at higher temperatures and had greater roughness in the 45°–75°, 105°–135°, and 155°–170° directions after different temperatures. Physicochemical reactions, such as dehydration, quartz phase transformation, and mineral oxidation, occurred in the granite owing to high temperatures. These reactions changed the crystalline state of the minerals and facilitated the increase and expansion of microcracks, thereby weakening the cementing property between mineral particles. Fracture surfaces with higher roughness were formed after tensile failure of the rock. The research results may provide theoretical references for the construction and operation of deep engineering.

Study on the mechanism of energy evolution and bearing degradation in pre-cracked sandstone under non uniform cyclic loading

The engineering rock mass is affected by the original cracks, which increases the uncertainty of instability and failure mechanisms. Under the action of excavation unloading or other disturbance loads, its bearing properties are more complex and the deformation characteristics is more special, which has always been a difficult issue in the analysis of rock mechanics. In order to explore the propagation and failure mechanism of spatial micro cracks in fractured rock masses, acoustic emission (AE) experiment of sandstone was carried out under pre-static loading and low-frequency stepwise cyclic loading–unloading. The research results indicate that the peak strain and crack initiation stress of rocks with different pre-cracked positions have the highest differences of 12.62% and 18.65%, respectively. The difference between peak stress and crack initiation strain is within 4%. The pre-cracked position mainly affects the peak strain and crack initiation stress of pre-cracked rocks, and its effect on peak stress and crack initiation strain is limited. The crack initiation angle of prefabricated cracks is mainly distributed within the range of 45° − 90°, and it decreases with the increase of quantities of prefabricated cracks. The crack initiation angle of rocks with prefabricated three cracks is 30% lower than that of single crack. The research results have deepened the analysis of energy evolution and deformation damage of fractured rock masses.

Development of correlations between various engineering rockmass classification systems using railway tunnel data in Garhwal Himalaya, India

Engineering rockmass classifications are an integral part of design, support and excavation procedures of tunnels, mines, and other underground structures. These classifications are directly linked to ground reaction and support requirements. Various classification systems are in practice and are still evolving. As different classifications serve different purposes, it is imperative to establish inter-correlatability between them. The rating systems and engineering judgements influence the assignment of ratings owing to cognition. To understand the existing correlation between different classification systems, the existing correlations were evaluated with the help of data of 34 locations along a 618-m-long railway tunnel in the Garhwal Himalaya of India and new correlations were developed between different rock classifications. The analysis indicates that certain correlations, such as RMR-Q, RMR-RMi, RMi-Q, and RSR-Q, are comparable to the previously established relationships, while others, such as RSR-RMR, RCR-Qn, and GSI-RMR, show weak correlations. These deviations in published correlations may be due to individual parameters of estimation or measurement errors. Further, incompatible classification systems exhibited low correlations. Thus, the study highlights a need to revisit existing correlations, particularly for rockmass conditions that are extremely complex, and the predictability of existing correlations exhibit high variations. In addition to augmenting the existing database, new correlations for metamorphic rocks in the Himalayan region have been developed and presented that can serve as a guide for future rock engineering projects in such formations and aid in developing appropriate excavation and rock support methodologies.

Study on energy evolution and fractal characteristics of sandstone with different fracture dip angles under uniaxial compression

In order to investigate the failure modes and instability mechanism of fractured rock. Uniaxial compression tests were conducted on sandstone specimens with different dip angles. Based on rock energy dissipation theory and fractal theory, the energy evolution characteristics and fragmentation fractal characteristics in the process of deformation and failure of specimens were analyzed. The results show that the peak strength and elastic modulus of fractured rock mass are lower than those of intact samples, and both show an exponential increase with the increase of fracture dip angle. The energy evolution laws of different fracture specimens are roughly similar and can be classified into four stages based on the stress–strain curve: pressure-tight, elastic, plastic, and post-destructive. The total strain energy, elastic strain energy, and dissipated strain energy of the specimen at the peak stress point increased exponentially with crack inclination, and the dissipated strain energy and compressive strength conformed to a power function growth relationship. The distribution of the fragments after the failure of the fracture sample has fractal characteristics, and the fractal dimension increases with the increase of the fracture dip angle. In addition, the higher the compressive strength of the specimen, the greater the energy dissipation, the more serious the degree of fragmentation, and the greater the fractal dimension. The data fitting further shows that there is a power function relationship between the dissipated strain energy and the fractal dimension. The research results can provide a theoretical basis for the stability of rock mass engineering and structural deformation control.

Investigation on the damage model of rock subjected to freeze-thaw cycles considering its fracture voids compaction

The investigation on the mechanical properties and damage model of rock subjected to freeze-thaw cycles (FTCs) holds significant theoretical for assessing the deformation of rock mass engineering in cold regions. When the rocks subjected to FTCs, the fracture voids of rocks increase with increasing number of FTCs, and the rocks exhibit nonlinear deformation characteristics as its fracture voids compaction. Based on the dissipation energy evolution law and distribution characteristics of rocks under loading, the critical point of rock compaction deformation was determined. A constitutive model for rock subjected to FTCs considering its fracture voids compaction was established based on macroscopic and microscopic damage levels. The rationality of the model was validated using experimental results, and a comparison was made between the calculation results of proposed model and model that not consider fracture voids compaction. The evolution properties of rock damage variable and stress under different FTCs were analyzed. The variation of model parameters with the FTCs and the influence of model parameters on the deformation characteristics of rock was analyzed.