Mechanical responses of anchoring structure under triaxial cyclic loading

Dynamic load on anchoring structures(AS)within deep roadways can result in cumulative damage and failure.This study develops an experimental device designed to test AS under triaxial loads.The device enables the investigation of the mechanical response,failure mode,instability assessment criteria,and anchorage effect of AS subjected to combined cyclic dynamic-static triaxial stress paths.The results show that the peak bearing strength is positively correlated with the anchoring matrix strength,anchorage length,and edgewise compressive strength.The bearing capacity decreases significantly when the anchorage direction is severely inclined.The free face failure modes are typically transverse cracking,concave fracturing,V-shaped slipping and detachment,and spallation detachment.Besides,when the anchoring matrix strength and the anchorage length decrease while the edgewise compressive strength,loading rate,and anchorage inclination angle increase,the failure intensity rises.Instability is determined by a negative tangent modulus of the displacement-strength curve or the continued deformation increase against the general downward trend.Under cyclic loads,the driving force that breaks the rock mass along the normal vector and the rigidity of the AS are the two factors that determine roadway stability.Finally,a control measure for surrounding rock stability is proposed to reduce the internal driving force via a pressure relief method and improve the rigidity of the AS by full-length anchorage and grouting modification.

Failure characterization of fully grouted rock bolts under triaxial testing

Confining stresses serve as a pivotal determinant in shaping the behavior of grouted rock bolts.Nonetheless,prior investigations have oversimplified the three-dimensional stress state,primarily assuming hydrostatic stress conditions.Under these conditions,it is assumed that the intermediate principal stress(σ_(2))equals the minimum principal stress(σ_(3)).This assumption overlooks the potential variations in magnitudes of in situ stress conditions along all three directions near an underground opening where a rock bolt is installed.In this study,a series of push tests was meticulously conducted under triaxial conditions.These tests involved applying non-uniform confining stresses(σ_(2)≠σ_(3))to cubic specimens,aiming to unveil the previously overlooked influence of intermediate principal stresses on the strength properties of rock bolts.The results show that as the confining stresses increase from zero to higher levels,the pre-failure behavior changes from linear to nonlinear forms,resulting in an increase in initial stiffness from 2.08 kN/mm to 32.51 kN/mm.The load-displacement curves further illuminate distinct post-failure behavior at elevated levels of confining stresses,characterized by enhanced stiffness.Notably,the peak load capacity ranged from 27.9 kN to 46.5 kN as confining stresses advanced from σ_(2)=σ_(3)=0 to σ_(2)=20 MPa and σ_(3)=10 MPa.Additionally,the outcomes highlight an influence of confining stress on the lateral deformation of samples.Lower levels of confinement prompt overall dilation in lateral deformation,while higher confinements maintain a state of shrinkage.Furthermore,diverse failure modes have been identified,intricately tied to the arrangement of confining stresses.Lower confinements tend to induce a splitting mode of failure,whereas higher loads bring about a shift towards a pure interfacial shear-off and shear-crushed failure mechanism.

Damage evolution of rock-encased-backfill structure under stepwise cyclic triaxial loading

Rock-encased-backfill(RB)structures are common in underground mining,for example in the cut-andfill and stoping methods.To understand the effects of cyclic excavation and blasting activities on the damage of these RB structures,a series of triaxial stepwise-increasing-amplitude cyclic loading experiments was conducted with cylindrical RB specimens(rock on outside,backfill on inside)with different volume fractions of rock(VF=0.48,0.61,0.73,and 0.84),confining pressures(0,6,9,and 12 MPa),and cyclic loading rates(200,300,400,and 500 N/s).The damage evolution and meso-crack formation during the cyclic tests were analyzed with results from stress-strain hysteresis loops,acoustic emission events,and post-failure X-ray 3D fracture morphology.The results showed significant differences between cyclic and monotonic loadings of RB specimens,particularly with regard to the generation of shear microcracks,the development of stress memory and strain hardening,and the contact forces and associated friction that develops along the rock-backfill interface.One important finding is that as a function of the number of cycles,the elastic strain increases linearly and the dissipated energy increases exponentially.Also,compared with monotonic loading,the cyclic strain hardening characteristics are more sensitive to rising confining pressures during the initial compaction stage.Another finding is that compared with monotonic loading,more shear microcracks are generated during every reloading stage,but these microcracks tend to be dispersed and lessen the likelihood of large shear fracture formation.The transition from elastic to plastic behavior varies depending on the parameters of each test(confinement,volume fraction,and cyclic rate),and an interesting finding was that the transformation to plastic behavior is significantly lower under the conditions of 0.73 rock volume fraction,400 N/s cyclic loading rate,and 9 MPa confinement.All the findings have important practical implications on the ability of backfill to support underground excavations.

Intermittent disturbance mechanical behavior and fractional deterioration mechanical model of rock under complex true triaxial stress paths

Mechanical excavation,blasting,adjacent rockburst and fracture slip that occur during mining excavation impose dynamic loads on the rock mass,leading to further fracture of damaged surrounding rock in three-dimensional high-stress and even causing disasters.Therefore,a novel complex true triaxial static-dynamic combined loading method reflecting underground excavation damage and then frequent intermittent disturbance failure is proposed.True triaxial static compression and intermittent disturbance tests are carried out on monzogabbro.The effects of intermediate principal stress and amplitude on the strength characteristics,deformation characteristics,failure characteristics,and precursors of monzogabbro are analyzed,intermediate principal stress and amplitude increase monzogabbro strength and tensile fracture mechanism.Rapid increases in microseismic parameters during rock loading can be precursors for intermittent rock disturbance.Based on the experimental result,the new damage fractional elements and method with considering crack initiation stress and crack unstable stress as initiation and acceleration condition of intermittent disturbance irreversible deformation are proposed.A novel three-dimensional disturbance fractional deterioration model considering the intermediate principal stress effect and intermittent disturbance damage effect is established,and the model predicted results align well with the experimental results.The sensitivity of stress states and model parameters is further explored,and the intermittent disturbance behaviors at different f are predicted.This study provides valuable theoretical bases for the stability analysis of deep mining engineering under dynamic loads.

Theoretical insights into thermal transport and structural stability mechanisms of triaxial compressed methane hydrate

The heat transfer and stability of methane hydrate in reservoirs have a direct impact on the drilling and production efficiency of hydrate resources,especially in complex stress environments caused by formation subsidence.In this study,we investigated the thermal transport and structural stability of methane hydrate under triaxial compression using molecular dynamics simulations.The results suggest that the thermal conductivity of methane hydrate increases with increasing compression strain.Two phonon transport mechanisms were identified as factors enhancing thermal conductivity.At low compressive strains,a low-frequency phonon transport channel was established due to the overlap of phonon vibration peaks between methane and water molecules.At high compressive strains,the filling of larger phonon bandgaps facilitated the opening of more phonon transport channels.Additionally,we found that a strain of0.04 is a watershed point,where methane hydrate transitions from stable to unstable.Furthermore,a strain of0.06 marks the threshold at which the diffusion capacities of methane and water molecules are at their peaks.At a higher strain of0.08,the increased volume compression reduces the available space,limiting the diffusion ability of water and methane molecules within the hydrate.The synergistic effect of the strong diffusion ability and high probability of collision between atoms increases the thermal conductivity of hydrates during the unstable period compared to the stable period.Our findings offer valuable theoretical insights into the thermal conductivity and stability of methane hydrates in reservoir stress environments.

Effects of pressure relief holes on coal burst:Insights from true-triaxial unloading tests

Drilling pressure relief is one of the methods to reduce the risk of coal bursts in deep mines.However,the effect of the drill hole orientations has not been studied well enough to understand their impact on the burst failure mechanism.In this study,we investigated two designs of drill hole orientations.The first design includes drill holes located on the upper free face of the rectangular samples and labelled as upper hole(UH)and centre hole(CH)e the long axes of the drill holes are aligned with minor principal stress,s3,direction.The second design includes drill holes at the top(TH)and the side(SH)of the rectangular samples in which the long axes of the drill holes are aligned with the maximum,s1,and intermediate principal stress,s2,directions,respectively.The coal samples with the proposed drill hole orientations were subjected to the true-triaxial unloading coal burst tests.The results show that the drill holes reduce the risk of coal bursts.However,we found that the intensity of coal burst was significantly reduced with the SH-type,followed by the CH-types.We also observed that the coal burst intensity is reduced better for the CH,UH,TH,and SH-type drilling patterns.However,it was found that the orientations of drill holes have little influence on the failure mode(splitting).The acoustic emission(AE)activities for coal with drill holes noticeably decreased,especially for the UH and CH layouts.The drill holes reduced the upper limit of the AE entropy(chaos of microcracks generation).However,regarding reducing the coal burst risk,the TH and SH are less effective than UH and CH.

Effect of fractures on mechanical behavior of sand powder 3D printing rock analogue under triaxial compression

In practical engineering applications,rock mass are often found to be subjected to a triaxial stress state.Concurrently,defects like joints and fractures have a notable impact on the mechanical behavior of rock mass.Such defects are identified as crucial contributors to the failure and instability of the surrounding rock,subsequently impacting the engineering stability.The study aimed to investigate the impact of fracture geometry and confining pressure on the deformation,failure characteristics,and strength of specimens using sand powder 3D printing technology and conventional triaxial compression tests.The results indicate that the number of fractures present considerably influences the peak strength,axial peak strain and elastic modulus of the specimens.Confining pressure is an important factor affecting the failure pattern of the specimen,under which the specimen is more prone to shear failure,but the initiation,expansion and penetration processes of secondary cracks in different fracture specimens are different.This study confirmed the feasibility of using sand powder 3D printing specimens as soft rock analogs for triaxial compression research.The insights from this research are deemed essential for a deeper understanding of the mechanical behavior of fractured surrounding rocks when under triaxial stress state.

Microdynamic mechanical properties and fracture evolution mechanism of monzogabbro with a true triaxial multilevel disturbance method

The far-field microdynamic disturbance caused by the excavation of deep mineral resources and underground engineering can induce surrounding rock damage in high-stress conditions and even lead to disasters.However,the mechanical properties and damage/fracture evolution mechanisms of deep rock induced by microdynamic disturbance under three-dimensional stress states are unclear.Therefore,a true triaxial multilevel disturbance test method is proposed,which can completely simulate natural geostress,excavation stress redistribution(such as stress unloading,concentration and rotation),and subsequently the microdynamic disturbance triggering damaged rock failure.Based on a dynamic true triaxial test platform,true triaxial microdynamic disturbance tests under different frequency and amplitudes were carried out on monzogabbro.The results show that increasing amplitude or decreasing frequency diminishes the failure strength of monzogabbro.Deformation modulus gradually decreases during disturbance failure.As frequency and amplitude increase,the degradation rate of deformation modulus decreases slightly,disturbance dissipated energy increases significantly,and disturbance deformation anisotropy strengthens obviously.A damage model has been proposed to quantitatively characterize the disturbance-induced damage evolution at different frequency and amplitude under true triaxial stress.Before disturbance failure,the micro-tensile crack mechanism is dominant,and the micro-shear crack mechanism increases significantly at failure.With the increase of amplitude and frequency,the micro-shear crack mechanism increases.When approaching disturbance failure,the acoustic emission fractal dimension changes from a stable value to local large oscillation,and finally increases sharply to a high value at failure.Finally,the disturbance-induced failure mechanism of surrounding rock in deep engineering is clearly elucidated.

Anisotropic strength,deformation and failure of gneiss granite under high stress and temperature coupled true triaxial compression

The anisotropic mechanical behavior of rocks under high-stress and high-temperature coupled conditions is crucial for analyzing the stability of surrounding rocks in deep underground engineering.This paper is devoted to studying the anisotropic strength,deformation and failure behavior of gneiss granite from the deep boreholes of a railway tunnel that suffers from high tectonic stress and ground temperature in the eastern tectonic knot in the Tibet Plateau.High-temperature true triaxial compression tests are performed on the samples using a self-developed testing device with five different loading directions and three temperature values that are representative of the geological conditions of the deep underground tunnels in the region.Effect of temperature and loading direction on the strength,elastic modulus,Poisson’s ratio,and failure mode are analyzed.The method for quantitative identification of anisotropic failure is also proposed.The anisotropic mechanical behaviors of the gneiss granite are very sensitive to the changes in loading direction and temperature under true triaxial compression,and the high temperature seems to weaken the inherent anisotropy and stress-induced deformation anisotropy.The strength and deformation show obvious thermal degradation at 200℃due to the weakening of friction between failure surfaces and the transition of the failure pattern in rock grains.In the range of 25℃ 200℃,the failure is mainly governed by the loading direction due to the inherent anisotropy.This study is helpful to the in-depth understanding of the thermal-mechanical behavior of anisotropic rocks in deep underground projects.

Experimental and numerical study on dynamic mechanical behaviors of shale under true triaxial compression at high strain rate

High-energy gas fracturing of shale is a novel,high efficacy and eco-friendly mining technique,which is a typical dynamic perturbing behavior.To effectively extract shale gas,it is important to understand the dynamic mechanical properties of shale.Dynamic experiments on shale subjected to true triaxial compression at different strain rates are first conducted in this research.The dynamic stress-strain curves,peak strain,peak stress and failure modes of shale are investigated.The results of the study indicate that the intermediate principal stress and the minor principal stress have the significant influence on the dynamic mechanical behaviors,although this effect decreases as the strain rate increases.The characteristics of compression-shear failure primarily occur in shale subjected to triaxial compression at high strain rates,which distinguishes it from the fragmentation characteristics observed in shale under dynamic uniaxial compression.Additionally,a numerical three-dimensional Split Hopkinson Pressure Bar(3D-SHPB),which is established by coupling PFC3D and FLAC3D methods,is validated to replicate the laboratory characteristics of shale.The dynamic mechanical characteristics of shale subjected to different confining stresses are systematically investigated by the coupling PFC3D and FLAC3D method.The numerical results are in good agreement with the experimental data.

Simulation of ductile fracture initiation in steels using a stress triaxiality-shear stress coupled model

Micromechanics-based models provide powerful tools to predict initiation of ductile fracture in steels. A new criterion is presented herein to study the process of ductile fracture when the effects of both stress triaxiality and shear stress on void growth and coalescence are considered. Finite-element analyses of two different kinds of steel, viz. ASTM A992 and AISI 1045, were carried out to monitor the history of stress and strain states and study the methodology for determining fracture initiation. Both the new model and void growth model (VGM) were calibrated for both kinds of steel and their accuracy for predicting fracture initiation evaluated. The results indicated that both models offer good accuracy for predicting fracture of A992 steel. However, use of the VGM leads to a significant deviation for 1045 steel, while the new model presents good performance for predicting fracture over a wide range of stress triaxiality while capturing the effect of shear stress on fracture initiation.

Numerical investigation of the mechanical behavior of the backfill–rock composite structure under triaxial compression

To ensure safe and economical backfill mining,the mechanical response of the backfill–rock interaction system needs to be understood.The numerical investigation of the mechanical behavior of backfill–rock composite structure(BRCS)under triaxial compression,which includes deformation,failure patterns,strength characteristics,and acoustic emission(AE)evolution,was proposed.The models used in the tests have one rough interface,two cement–iron tailings ratios(CTRs),four interface angles(IAs),and three confining pressures(CPs).Results showed that the deformation,strength characteristics,and failure patterns of BRCS under triaxial compression depend on IA,CP,and CTR.The stress–strain curves of BRCS under triaxial compression could be divided into five stages,namely,compaction,elasticity,yield,strain softening,and residual stress.The relevant AE counts have corresponding relationships with different stages.The triaxial compressive strengths of composites increase linearly with the increase of the CP.Furthermore,the CP stress strengthening effect occurs.When the IAs are45°and 60°,the failure areas of composites appear in the interface and backfill.When the IAs are 75°and 90°,the failure areas of composites appear in the backfill,interface,and rock.Moreover,the corresponding failure modes yield the combined shear failure.The research results provide the basis for further understanding of the stability of the BRCS.

Experimental study of the dynamic mechanical responses and failure characteristics of coal under true triaxial confinements

Investigations on the dynamic mechanical properties and failure mechanisms of coal under in-situ stress is essential for the prevention of dynamic disasters in deep coal mines.Thus,a modified true triaxial Hopkinson bar was employed to explore the dynamic mechanical behaviors of coal at different confining pressures(0–20 MPa)and strain rates(40–220 s^(-1)).The results show that the dynamic peak stress is positively correlated with lateral static pre-stressσy andσz,but negatively correlated with axial static prestressσx.At approximate strain rates,increasing the lateral static pre-stress facilitates increasing the dynamic peak stress,but the minimum lateral static pre-stress is the primary factor limiting a significant increase in dynamic peak stress of coal.Furthermore,the dynamic differential stress is linearly related to the logarithm of strain rate,and the peak strain varies linearly with strain rate.However,there is no significant correlation between confining pressure and peak strain.Moreover,X-ray CT images and photographic fracture observations of coal samples show the failure patterns under uniaxial and triaxial conditions are splitting failure and shear failure,respectively.The device provides a viable approach for fully comprehending the dynamic mechanical behaviors of rock-like material in complex stress conditions.

The Influence of Triaxiality Parameter γ on the Chiral Doublet Bands with (πg_(9/2))～(-1)(νh_(11/2))～2 Configuration

The chiral doublet bands with three-quasiparticle configuration （πg9/2）-1 （vh11/2）2 are studied by the fully quantal triaxial particle rotor model. The energy spectra and B（M1）/B（E2） ratios of the doublet bands with different triaxiality parameter γ are systematically analyzed. It is found that γ is a sensitive parameter for the properties of these doublet bands.

Liquefaction proneness of stratified sand-silt layers based on cyclic triaxial tests

Most studies on liquefaction have addressed homogeneous soil strata using sand or sand with fine content without considering soil stratification.In this study,cyclic triaxial tests were conducted on the stratified sand specimens embedded with the silt layers to investigate the liquefaction failures and void-redistribution at confining stress of 100 kPa under stress-controlled mode.The loosening of underlying sand mass and hindrance to pore-water flow caused localized bulging at the sand-silt interface.It is observed that at a silt thickness of 0.2H(H is the height of the specimen),nearly 187 load cycles were required to attain liquefaction,which was the highest among all the silt thicknesses with a single silt layer.Therefore,0.2H is assumed as the optimum silt thickness(t_(opt)).The silt was placed at the top,middle and bottom of the specimen to understand the effect of silt layer location.Due to the increase in depth of the silt layer from the top position(capped soil state)to the bottom,the cycles to reach liquefaction(N_(cyc,L))increased 2.18 times.Also,when the number of silt layers increased from single to triple,there was an increase of about 880%in N_(cyc,L).The micro-characterization analysis of the soil specimens indicated silty materials transported in upper sections of the specimen due to the dissipated pore pressure.The main parameters,including thickness(t),location(z),cyclic stress ratio(CSR),number of silt layers(n)and modified relative density(D_(r,m)),performed significantly in governing the lique-faction resistance.For this,a multilinear regression model is developed based on critical parameters for prediction of N_(cyc,L).Furthermore,the developed constitutive model has been validated using the data from the present study and earlier findings.

Mutual impact of true triaxial stress, borehole orientation and bedding inclination on laboratory hydraulic fracturing of Lushan shale

Unconventional resources like shale gas has been the focus of intense research and development for two decades. Apart from intrinsic geologic factors that control the gas shale productivity (e.g. organic matter content, bedding planes, natural fractures, porosity and stress regime among others), external factors like wellbore orientation and stimulation design play a role. In this study, we present a series of true triaxial hydraulic fracturing experiments conducted on Lushan shale to investigate the interplay of internal factors (bedding, natural fractures and in situ stress) and external factors (wellbore orientation) on the growth process of fracture networks in cubic specimens of 200 mm in length. We observe relatively low breakdown pressure and fracture propagation pressure as the wellbore orientation and/or the maximum in situ stress is subparallel to the shale bedding plane. The wellbore orientation has a more prominent effect on the breakdown pressure, but its effect is tapered with increasing angle of bedding inclination. The shale breakdown is followed by an abrupt response in sample displacement, which reflects the stimulated fracture volume. Based on fluid tracer analysis, the morphology of hydraulic fractures (HF) is divided into four categories. Among the categories, activation of bedding planes (bedding failure, BF) and natural fractures (NF) significantly increase bifurcation and fractured areas. Under the same stress regime, a horizontal wellbore is more favorable to enhance the complexity of hydraulic fracture networks. This is attributed to the relatively large surface area in contact with the bedding plane for the horizontal borehole compared to the case with a vertical wellbore. These findings provide important references for hydraulic fracturing design in shale reservoirs.

A testing system to understand rock fracturing processes induced by different dynamic disturbances under true triaxial compression

In this context,a testing system to understand rock fracturing processes induced by different dynamic disturbances under true triaxial compression was developed.The system is mainly composed of a static loading subsystem,a dynamic loading subsystem,a specimen box subsystem,and a data measurement subsystem.The static loading subsystem uses low stiffness loss frame structure technology,which greatly improves the frame stiffness in the three principal stress directions(up to 20 GN/m)and ensures the demand of the disturbance experiment in both the prepeak and postpeak stages.The disturbance loads with frequency of 0e20 Hz and stress level of 0e30 MPa were applied using large flow parallel oil source technology characterized with high heat dissipation efficiency.For the disturbance loads with frequency of 100e500 Hz and stress level of 0e30 MPa,they were realized by using high-frequency and centimeter-per-second-scale low-speed disturbance rod technology.Three rigid self-stabilizing specimen boxes were utilized to provide support for the specimen and deformation sensors,ensuring the stability and accuracy of the data obtained.To verify the performance of the true triaxial test system,disturbance experiments were conducted on granite specimens.The results show that the experimental device satisfies the requirements of original design,with an excellent repeatability and reliable testing results.

Effects of intermediate stress on deep rock strainbursts under true triaxial stresses

The effect of intermediate stress(in situ tunnel axial)on a strainburst is studied with a threedimensional(3D)bonded block distinct element method(DEM).A series of simulations of strainbursts under true triaxial in situ stress conditions(i.e.high tangential stress,moderate intermediate stress and low radial stress)of near-boundary rock masses are performed.Compared with the experimental results,the DEM model is able to capture the stress-strain response,failure pattern and energy balance of strainbursts.The fracturing processes of strainbursts are also numerically reproduced.Numerical results show that,as the intermediate stress increases:(1)The peak strain of strainbursts increases,the yield stress increases,the rock strength increases linearly,and the ratio of yield stress to rock strength decreases,indicating that the precursory information on strainbursts is enhanced;(2)Tensile and shear cracks increase significantly,and slabbing and bending of rock plates are more pronounced;and(3)The stored elastic strain energy and dissipated energy increase linearly,whereas the kinetic energy of the ejected rock fragments increases approximately exponentially,implying an increase in strainburst intensity.By comparing the experimental and numerical results,the effect of intermediate stress on the rock strength of strainbursts is discussed in order to address three key issues.Then,the Mogi criterion is applied to construct new strength criteria for strainbursts by converting the one-face free true triaxial stress state of a strainburst to its equivalent true triaxial stress state.In summary,the effect of intermediate stress on strainbursts is a double-edged sword that can enhance the rock strength and the precursory information of a strainburst,but also increase its intensity.

Predicting triaxial compressive strength of high-temperature treated rock using machine learning techniques

The accurate prediction of the strength of rocks after high-temperature treatment is important for the safety maintenance of rock in deep underground engineering.Five machine learning(ML)techniques were adopted in this study,i.e.back propagation neural network(BPNN),AdaBoost-based classification and regression tree(AdaBoost-CART),support vector machine(SVM),K-nearest neighbor(KNN),and radial basis function neural network(RBFNN).A total of 351 data points with seven input parameters(i.e.diameter and height of specimen,density,temperature,confining pressure,crack damage stress and elastic modulus)and one output parameter(triaxial compressive strength)were utilized.The root mean square error(RMSE),mean absolute error(MAE)and correlation coefficient(R)were used to evaluate the prediction performance of the five ML models.The results demonstrated that the BPNN shows a better prediction performance than the other models with RMSE,MAE and R values on the testing dataset of 15.4 MPa,11.03 MPa and 0.9921,respectively.The results indicated that the ML techniques are effective for accurately predicting the triaxial compressive strength of rocks after different high-temperature treatments.

A true triaxial strength criterion for rocks by gene expression programming

Rock strength is a crucial factor to consider when designing and constructing underground projects.This study utilizes a gene expression programming(GEP)algorithm-based model to predict the true triaxial strength of rocks,taking into account the influence of rock genesis on their mechanical behavior during the model building process.A true triaxial strength criterion based on the GEP model for igneous,metamorphic and magmatic rocks was obtained by training the model using collected data.Compared to the modified Weibols-Cook criterion,the modified Mohr-Coulomb criterion,and the modified Lade criterion,the strength criterion based on the GEP model exhibits superior prediction accuracy performance.The strength criterion based on the GEP model has better performance in R2,RMSE and MAPE for the data set used in this study.Furthermore,the strength criterion based on the GEP model shows greater stability in predicting the true triaxial strength of rocks across different types.Compared to the existing strength criterion based on the genetic programming(GP)model,the proposed criterion based on GEP model achieves more accurate predictions of the variation of true triaxial strength(s1)with intermediate principal stress(s2).Finally,based on the Sobol sensitivity analysis technique,the effects of the parameters of the three obtained strength criteria on the true triaxial strength of the rock are analysed.In general,the proposed strength criterion exhibits superior performance in terms of both accuracy and stability of prediction results.