1. Introduction
Rocks are naturally complex geomaterials. Its mechanical properties are significantly influenced by factors such as stress state, lithology and other geological characteristics. Rock strength criteria can quantitatively describe the relationship among rock strength, stress state and field variables. The strength criteria guide the design of ultimate bearing capacity and safety evaluation in practical engineering and can be used to establish the constitutive relationship of rock materials. So far, many rock strength criteria have been proposed, which can be categorized [] into the classical criteria [,,,,,,], the modified criteria [,,,,,,,,,,,,,,], the unified criteria [,,,,,,,,,].
The classical criteria include the Mohr-Coulomb criterion, the Hoke-Brown criterion [,], the Lade-Duncan criterion [,], the spatial mobilized plane (SMP) criterion [,,] and etc. The Mohr-Coulomb (M-C) criterion and the Hoek-Brown (H-B) criterion are two of the most commonly used criteria for describing rock strength characteristics. Numerous researches have demonstrated that the intermediate principal stress has the notable effect on rock strength [,,,,,,,,,,,], which these two criteria fail to characterize. The Lade criterion is suitable for cohesive geomaterials, but the failure curve on deviatoric plane remains unchanged under a given internal friction angle []. The SMP criterion neglects the strength nonlinearity of rock on meridian plane.
Most of the modified criteria is based on the classical criteria, such as the Paul-Mohr-Coulomb criterion [], the 3D Hoke-Brown criteria [,,,,,,,,,,,,], the Zienkiewicz–Pande criterion [] and etc. The modification of the classical criteria for rocks has expanded its application in characterizing the strength characteristics under complex stress state. The Paul-Mohr-Coulomb criterion (P-M-C) proposed by Paul is a 3D linear strength criterion. The two sets of the P-M-C criterion can approximately delineate the nonlinear strength characteristics of rock. Many scholars have adopted the H-B criterion as the foundation for studying 3D rock strength criteria, primarily because it reflects the nonlinear strength of rock on the meridian plane. However, the current 3D H-B criteria have two main limitations. The first is that the strength criteria fail to predict the same rock strength under triaxial compression or triaxial tension as the original H-B criterion [,,,,]. The second is that the failure surfaces of the strength criteria in the principal stress space cannot simultaneously meet smoothness and convexity requirements [12-14, 16, 18, 21, 22]. The Zienkiewicz–Pande criterion is essentially a type of modified M-C criterion, and its failure curve is quadratic on meridian plane and curvilinear on deviatoric plane. For rocks, as the hydrostatic pressure increases, the failure curve on deviatoric plane changes from a curved triangle to a circle, which cannot be characterized by the Zienkiewicz–Pande criterion [].
The unified strength criteria amalgamate various individual criteria. Yu et al. [,,] proposed a unified strength theory (UST) based on the twin-shear model, which can characterize the effect of the intermediate principal stress. The UST has been widely applied in geotechnical engineering [,]. However, the UST predicts the same rock strength under both triaxial tension and triaxial compression, which conflicts with the actual strength of some rocks and limits its application. By introducing a modified parameter χ into the UST to correcting the rock strength under triaxial tension, a generalized unified strength theory (GUST) was proposed []. Although the GUST predicts the strength of various types of rocks more accurately, its failure surface in the principal stress space is non-smooth, leading to inconvenience of the numerical computation. Yao et al. [,] proposed the unified strength criterion (USC) for geomaterials. The USC was established based on nonlinear characteristics of rock on both meridian and deviatoric planes. The stress transformation method [,] was adopted for incorporating the failure functions on both meridian and deviatoric planes. Four parameters of the USC control the shape of failure curve, which can effectively characterize the shape transition of the failure curve from a curved triangle to a circle. Liu et al. [] proposed a 3D nonlinear strength criterion (3D NSC) for rocks by combining the segmented meridian function and the generalized deviatoric function with one parameter. The 3D NSC can accurately describe and predict the strength change of rock in both brittle and ductile domains.
The rock criteria can also be categorized into the linear strength theory and the nonlinear strength theory []. A criterion can be classified as the linear strength theory only if the failure curves are linear on both meridian and deviatoric planes, whereas if the failure curve is nonlinear on either meridian plane or deviatoric plane, it is classified as the nonlinear strength theory. The linear strength theory includes the M-C criterion, the twin-shear stress theory [], etc. The nonlinear strength theory includes the Hoek-Brown (H-B) criterion [,], the SMP criterion [,,], the Lade criterion [,], etc. The failure curves of the M-C criterion and the twin-shear theory are linear on the two planes. The failure curve of H-B criterion is parabolic on meridian plane and linear on deviatoric plane, respectively. The SMP criterion appears as a linear function on meridian plane, and curvilinear shape on deviatoric plane. The Lade criterion also has a linear relationship on meridian plane and a curvilinear relationship on deviatoric plane. The Zienkiewicz–Pande criterion presents failure curve with quadratic forms on meridian plane and the curvilinear form on deviatoric plane.
The strength test methods for rocks under different stress conditions mainly include the conventional triaxial test (CTT) and true triaxial test (TTT). The CTT was perhaps first conducted by Karman under the stress state that the intermediate and minimum principal stresses are equal []. Due to the simplicity of the equipment and ease of sample preparation, the CTT has been widely applied in rock mechanics experiments [,,,]. The CCT is suitable for testing the strength of a type of rock with the weak intermediate principal stress effect. The TTT apparatus designed by Mogi may first enable the application of three mutually independent and uniform loads to the specimen faces []. The TTT conducted by Mogi found that the maximum principal stress σ1 at failure is a function of σ2 with concave curve of σ1 vs. σ2 under constant σ3. The theories and applications based on the TTT have been verified by practical engineering cases such as coal mining [], wellbore stability [] and rock-burst [,]. Therefore, both of the CTT and TTT are the effective methods for measuring rock strength. In particular, the TTT can measure the rock strength under more complex stress state.
In summary, a versatile rock strength criterion needs to meet the following functions: (i) the influence of the intermediate principal stress on rock strength is effectively characterized; (ii) the failure surface of the strength criterion in the principal stress space is smooth and convex; (iii) the parameters of the strength criterion have clear physical meanings; (iv) the criterion can be conveniently integrated with elastoplastic constitutive models. To determine the advantages and disadvantages of the current typical rock strength criteria, three types of rock strength criteria were compared in this paper. These rock strength criteria are the USC [], the GUST [] and 3D H-B []. By comparing the USC with the GUST and the 3D H-B, the advantages of the USC are determined. Finally, based on the conventional triaxial test results for eight types of rock, the accuracy of rock strength predicted by the original H-B criterion and the USC was compared.
2. Three typical rock strength criteria
2.1. 3D H-B criterion
Based on Griffith's strength theory and combined with a large amount of triaxial test for rocks, the original H-B criterion was first proposed by Hoek and Brown in 1980 []. It can be expressed as
where σ1 is the maximum principal stress; σ3 is the minimum principal stress; σc is the uniaxial compressive strength of rock; mi is a constant reflecting the degree of softness and hardness of rock, which ranges from 0.001 to 25.0. The failure curve of the H-B strength criterion in meridian plane is parabolic. It can be expressed as
(2)
where p is the hydrostatic pressure; q is the generalized shear stress.
Based on the SMP criterion and the Drucker-Prager (D-P) criterion, Du et al. [] developed the nonlinear unified strength theory (NUST). A new stress space, β stress space, is adopted in the NUST. On deviatoric plane, the NUST can be expressed as,
where pβ is the hydrostatic pressure in β stress space, and pβ0 is the initial value of pβ; σβi (i=1,2,3) is the principal stress in β stress space.
Considering the effect of the intermediate principal stress, Huang et al. [] integrated the NUST with the H-B criterion by the stress-space transformation method. Therefore, the 3D H-B criterion applicable to describing strength characteristics of rock under triaxial stress state was established. The 3D H-B criterion is defined as,
(4)
where Mf is a material parameter; Mβ is the failure stress ratio; mt is the ratio of the triaxial tension strength and the triaxial compression strength under the same hydrostatic stress and calculated by Mβ, Mf and β.
The 3D H-B criterion constructs the strength curves on meridian and deviatoric planes based on the H-B criterion and the NUST, respectively. By employing a stress-space transformation method, it combines the strength curves on meridian and deviatoric planes, thereby establishing a new 3D rock strength criterion.
2.2. Generalized unified strength theory
Yu and He proposed a unified strength theory (UST) based on a twin-shear model []. It can be expressed as follows:
when τ12 + βσ12 ≥ τ12 + βσ23
(5a)
when τ12 + βσ12 ≤ τ12 + βσ23
(5b)
where b is a parameter reflecting the influence of τ12 and τ23 on the failure of materials; β is the coefficient reflecting the effect of normal stresses on failure; C is a material constant. The values of β and C can be determined by uniaxial tension strength σt and uniaxial compression strength σc respectively.
The UST can be expressed by the principal stresses as follows:
when σ2 ≤ (σ1 + ασ3)/(1 + α)
when σ2 ≥ (σ1 + ασ3)/(1 + α)
where α=σt / σc. The UST expressed by the principal stresses from Eq. 6 incorporate the intermediate principal stress σ2, which implies that the strength characteristics of the rock can be characterized more accurately.
Rock strength predicted by the UST are equal under triaxial compression and triaxial tension, which limits its application. Li et al. [] modified the UST by introducing a parameter χ and proposed the generalized unified strength theory (GUST). The expression of the GUST is as follows:
when σ2 ≥ σ*2
(7a)
when σ2 ≤ σ*2
(7b)
where σ*2 is the turning point of the σ*2 effect corresponding to the maximum strength, σ*2=σ1 [(1-sinφ) + σ3 (1-sinφ)]/[(1+χ) + (1-χ)sinφ]; the factor χ characterizes the mobilization of the weakening effect of σ2 when σ2 ≥ σ*2, and it can be determined from the conventional triaxial tension test,
(8a)
or
(8b)
where φe is the friction angle under triaxial tension. As explained above, the strength characteristics depends on the mobilization of σ2-effect in the complex stress state.
The GUST is based on the twin-shear model. It introduces the parameter χ, which enables the twin-shear model to be applied to some types of rock with higher triaxial tensile strength than triaxial compressive strength.
2.3. Unified strength criterion
Based on experimental law, the rock strength exhibited nonlinear behavior both on deviatoric plane and meridian plane []. The USC was proposed for characterizing the nonlinear behavior of the rock strength. On deviatoric plane, the USC unifies the SMP criterion and the Mises criterion by linear interpolation.
1. The failure curve of the Mises criterion on π-plane in original stress space can be expressed as follows:
2. The failure curve of the SMP criterion on π-plane in original stress space can be expressed as follows:
(10)
where q* is deviator stress; I1, I2 and I3 represent the first, second and third stress invariant, respectively;
A linear interpolation of the SMP criterion and the Mises criterion is conducted by introducing a material parameter α. The expression of the USC can be expressed as follows:
As shown in , α is a factor controlling the shape of the USC on π-plane. α varies within the range of 0 to 1. α=0 and α=1 represent the SMP criterion and the Mises criterion respectively. Therefore, the shape of the USC on deviatoric plane can transmit from a curved triangle to a circle with α increasing.
On meridian plane, a failure function of the USC can be established as a power function [],
where p is mean stress; pr is given reference pressure; Mf and n are material parameters; σ0 is triaxial tension strength.
Notice that if establishing the strength criterion by combining the failure functions on deviatoric and meridian planes, one of the two points should be satisfied at least: (i) the failure function on meridian plane is linear; (ii) the failure curve on deviatoric plane is circular. The failure curves of the USC on deviatoric and meridian planes are a curve triangle and nonlinear respectively. Thus, a transformed method is adopted for the failure function on meridian plane of the USC, which can transform the nonlinear failure curve into the linear failure curve. The failure curves on deviatoric planes in original and transition spaces are shown in . The failure curve passes a fixed point (pr, Mf pr). Transformation relationship is defined as,
(13)
where is mean stress in transformed stress space.
The stress tensor in transition space is expressed as,
(14)
where i (i=1,2,3) represents the principal stress in transformed stress space. In transformed stress space, the new unified strength criterion, the USC, satisfies a condition that is the linear failure function on meridian plane. Therefore, the USC can be established by combing the failure functions on deviatoric and meridian planes (3 > 0),
(16)
The USC constructs the strength curve on meridian and deviatoric planes based on the empirical power function and a linear combination of the SMP and Mises criteria, respectively. It also combines the strength curves on these two planes by using a stress-space transformation method.
3. Rock strength test
3.1. True triaxial test
The true triaxial test results are quoted from the references [,,]. The USC was verified by the true triaxial test from five types of rock including Dunham dolomite [] , KTB amphibolite [], Shirahama sandstone [], Solnhofen limestone [] and Yuubari shale []. The true triaxial tests were conducted under different σ3 levels. The tests data are shown in , , , , .
3.2. Conventional triaxial test
3.2.1. Test equipment
The conventional triaxial tests were conducted on WDT-1500 Rock Test System with an axial load capacity of 1500 kN, confining pressure capacity of 80 MPa (see ). The WDT-1500 controller consists of a hardware component and software applications. It includes four main parts: a digital servo controller, an axial dynamic loading system, a self-balanced pressure chamber and a data acquisition system. The computer is used to control the host to run the experiment.
3.2.2. Test Specimen
Eight types of rock, green sandstone, tuff, mudstone, fine sandstone, cyan sandstone, granite, grey sandstone and micritic limestone are selected to test strength in this paper. The rock samples were obtained from the surrounding rock of a tunnel in Qinlin Mountain, Shaanxi Province, China. The samples were cut into standard cylinders 50 mm in diameter and 100 mm in length according to the ISRM testing procedure and guidelines.
3.2.3. Test method
Uniaxial compression strength tests and conventional triaxial tests were conducted on each type of rock. During the tests, a constant confining pressure was applied to the rock samples, and the axial pressure was gradually increased until the rock failed. The loading rate of axial force was set at 0.1 kN/s. The confining pressures applied to green sandstone and tuff were 2, 7, 10, and 13 MPa, and 4, 8, 12, and 16 MPa, respectively. For the rest rock samples, the confining pressures were 5, 10, 15, and 20 MPa. The test results were summarized in .
4. Comparison among the USC, the GUST and 3D H-B criterion
All parameters of the USC are shown in . The results of true triaxial tests illustrate the strength variation of five types of rock from triaxial compression to triaxial tension. From Figure 4, the following conclusions can be drawn: (1) σ2 significantly affects rock strength; (2) σ2 effect on rock strength exhibits an interval pattern, which refers that as σ2 increases from the lower limit value (σ2=σ3) to the upper limit value (σ1=σ2), rock strength initially increases, then decreases once σ2 reaches a certain value; (3) rock strength increases as σ3 increases; (4) rock strength under triaxial tension is slightly greater than that under triaxial compression. (5) the degree of σ2 effect of Yuubari shale and Shirahama sandstone is smaller than other types of rock, which means the degree of σ2 effect varies among different types of rock.
The 3D H-B criterion and USC describe the strength characteristics of rock in a nonlinear form, while GUST describe those in a piecewise linear form. For Dunham dolomite (see a), σ3 ranges from 25 to 125 MPa with a constant interval of 20 MPa. The value of σ2 corresponding to the maximum strength predicted by the USC is greater than that predicted by the other two criteria, for example, when σ3=125 MPa, σ2USC=550 MPa, σ23D H-B=426 MPa, σ2GUST=370 MPa. The maximum increase ratio of rock strength from triaxial compression to triaxial tension predicted by the USC, the GUST, and the 3D H-B criterion is 34.01%, 17.31%, and 14.42%, respectively.
For KTB amphibolite, σ3 ranges from 0 to 150 MPa with intervals of 30 MPa, 40 MPa, and 50 MPa (see b). The length of the strength increase section predicted by the USC is still greater than that of the other two criteria. Among the three criteria, the rock strength predicted by the 3D H-B under triaxial tension is lower than that under triaxial compression when σ3 is greater than 0 MPa, while the other criteria maintain the prediction that the triaxial tension strength is greater than the triaxial compression strength. The 3D H-B criterion shows a significant drop in rock strength after reaching the maximum value when σ3 is greater than 30 MPa, resulting in an obvious error. Compared with other criteria, the rock strength predicted by the USC and the actual values have good fitting results.
For Shirahama sandstone, σ3 ranges from 8 to 40 MPa with intervals of 5 MPa, 7 MPa, and 10 MPa (see c). When the values of are 8 MPa, 15 MPa, and 20 MPa, the rock strength predicted by the USC criterion is lower than the actual values near triaxial compression. When the values of σ3 are 30 MPa and 40 MPa, the rock strength predicted the 3D H-B criterion is significantly lower than the actual values near the triaxial tension state. The errors between the rock strength predicted by the GUST and the actual values are not obvious. For Solnhofen limestone and Yuubari shale, none of the three criteria shows significant errors in fitting rock strength (see d, e).
In order to compare the USC with other strength criterions quantitatively, the root-mean-square error (RMSE, MPa) between the predicted strength and the tested data is defined as,
(17)
where the unit of RMSE is MPa; σ1itest and σ1icalc are tested data and calculated data of i-th the first principal stress respectively; n is the total number of test samples.
As shown in , the RMSE ranges for the USC, the 3D H-B criterion, and the GUST are 7.35 MPa to 69.12 MPa, 8.74 MPa to 108.32 MPa, and 8.63 MPa to 109.12 MPa respectively, which indicate that USC is the optimal criterion in predicting rock strength. It is known that the effect of σ2 on rock strength can be divided into two parts: a strengthening region at low σ2 values, and a weakening region at high σ2values. In addition, the triaxial tension strength of rock is slightly greater than the triaxial compression strength.
The USC accurately characterizes the effect of σ2 due to integration of the nonlinear characteristics of rock strength both on deviatoric plane and meridian plane. The GUST characterizes the effect of σ2 by introducing the parameter b, and adjusts the fitting results by introducing the mobilized parameter χ to make the triaxial tension strength of the rock greater than the triaxial compression strength. For the 3D H-B criterion, when the stress state of rock is approaching triaxial tension state, the errors between fitting values and actual values become more significant, particularly in KTB amphibolite. The deficiency of 3D H-B criterion is mainly related to the overestimation of the weakening effect at high σ2 values.
5. Rock strength predicted by the USC on meridian plane
To further validate the applicability of the USC for different types of rock, the conventional triaxial tests were conducted on eight types of rock, including green sandstone, tuff, mudstone, fine sandstone, cyan sandstone, granite, grey sandstone and micritic limestone. Mf, n, σ0 and pr, were determined by the methods provided in reference [].
The test results show that the deviatoric stress q increases monotonically with an increase of the hydrostatic pressure p (see ). The failure curve fitted by the USC on meridian plane includes linear and nonlinear forms, which are controlled by the parameter n. The failure curve fitted by the original H-B criterion only has a nonlinear form. The intersection of the failure curve of rock on meridian plane with the p-axis represents the true triaxial tensile strength. In the USC, the true triaxial tensile strength of the rock can be evaluated directly [], whereas the original H-B criterion cannot provide the reasonable true triaxial tensile strength. For example, the true triaxial tensile strength of Fine sandstone calculated by the USC criterion is -3.97 MPa, while the true triaxial tensile strength of Fine sandstone calculated by the original H-B criterion is -26.36 MPa, which is 6.64 times than that calculated by the USC.
The RMSE of q for the USC and the original H-B criterion has been calculated according to Eq. (21). shows the RMSE of q for all types of rock. The RMSE of q for the USC and the original H-B criterion ranges from 0.69 MPa to 4.9 MPa and 0.72 MPa to 7.88 MPa respectively. The RMSE interval of q for the USC falls within that of the original H-B criterion, which suggests that the USC can predict the different types of rock more accurately than the original H-B criterion. On the other hand, the USC overcomes the shortcoming of some 3D H-B criteria, which cannot predict the same rock strength as the H-B criterion [,,,,].
6. Discussion
The GUST, the 3D H-B criterion, and the USC all belong to the unified strength theory. Each of the three criteria can be degenerated into a single strength criterion by changing values of the relevant parameters.
When b=0, the GUST reduces to the Mohr-Coulomb criterion. When b=1/2 and (1-sinφ)/(1+sinφ)=0, the GUST becomes the linear approximation of the Mises criterion. When b=0 and (1-sinφ)/(1+sinφ)=1, the GUST becomes the Tresca criterion. For the 3D H-B criterion, when mt=1, it reduces to the original H-B criterion; it becomes the SMP criterion with β=0 and Drucker-Prager criterion with β=1. For the USC criterion, it reduces to the SMP criterion with α=0 and the Mises criterion with α=1.
The parameters of the three criteria are determined through the true triaxial tests or conventional triaxial tests. The GUST includes four parameters c, φ, χ and b. c and φ can be determined by conventional triaxial tests. b and χ can be determined jointly using the least square method based on the results of true triaxial tests []. The 3D H-B criterion includes six parameters σc, mi, mt, Mβ, Mf and β. These six parameters are not mutually independent. As long as mt is determined, the Mβ, Mf and β will be determined. In general, mt is determined by true triaxial tests. When some types of rock satisfy the negative exponential empirical relationship between mi and mt, mt =[(mi+1)-0.371+1]/2, the five parameters of the 3D H-B can be reduced to two parameters σc and mi []. σc and mi can be determined by conventional triaxial tests. The parameters of USC are determined by true triaxial tests. Therefore, the difficulty level of the method for determining the parameters among these three criteria is similar to the types of rock.
The USC can unify the strength characteristics of different materials by adjusting the values of four independent parameters Mf, σ0, n and α. As shown in , σ0=0 and σ0≠0 in the USC represents cohesiveless and cohesive materials respectively (see a, c, e), while 0<n≤1 in the USC represents frictional materials (see b, d, f). Different combinations of σ0, n and α can be used for describing the strength characteristics of some specific materials. Note that the color of the 3D failure surface serves only for rendering purposes, providing a good visual effect without the specific meaning. The parameter value can be classified into three categories.
The first category for parameter value is that n=1, α=0 (see a, b). When σ0=0, the USC reduces to the SMP criterion, which can be applied to sands and saturated clays (see a). When σ0≠0, cohesion is taken into account in the USC, which is identical to the extended SMP criterion []. In this case, the USC is appropriate for cohesive unsaturated soils and cemented sands (see b).
The second category for parameter value is that n=1, 0<α≤1 (see c, d). The failure surface of the USC changes from curved triangle to circle with α increasing. When σ0=0, the USC can be used to describe the strength characteristics of sands and saturated clays (see c). When σ0≠0, the USC is appropriate for cohesive unsaturated soils and cemented sands (see d).
The third category for parameter value is that 0< n<1, 0<α<1 (see e, f). The coupling effect of the strength on meridian plane and deviatoric plane can be taken into account. When σ0≠0, the USC can be applied to nonlinear cohesiveless materials, such as gravels (see e). When σ0≠0, the USC is in the general form. In this case, the material is cohesive, and the internal friction angle varies nonlinearly with hydrostatic pressure increasing (see f). Thus, the applicable scope of the USC can be expanded under, which means that the strength characteristics of rock and concretes can be described by the USC under as well.
In conclusion, the USC can be applied for describing the strength characteristics of soils, gravels, concretes and rocks by adjusting the value of relevant parameters.
Compared to the original H-B criterion, the GUST, and the 3D H-B criterion, the USC not only clearly reflects the shear yield and failure characteristics of rocks under three-dimensional stress condition but also facilitates the integration with the specific elastoplastic constitutive models [,,]. There are numerous structural planes in rock masse that can affect the stability of engineering projects. To more accurately assess stability of the projects, the primary task is to accurately predict the strength of the rock mass. Although the USC has demonstrated good accuracy in predicting the strength of intact rock, it has not been proven effective in predicting rock mass strength. Therefore, in future research, the focus will be on predicting rock mass strength and establishing the corresponding elastoplastic constitutive model.
7. Conclusion
In this paper, three typical types of rock strength criteria are selected to fit the true triaxial experimental results of five types of rock from the reference. The pros and cons of the USC, 3D H-B, and GUST criteria are compared. Based on the conventional triaxial test results of eight types of rock, the accuracy of the original H-B criterion and the USC in predicting rock strength on meridian plane is compared. Finally, the meaning and determination method of the parameters in each criterion is discussed.
1. The USC is constructed by combining the failure functions on deviatoric and meridian planes. The USC can predict the strength of various geomaterials. The 3D H-B criterion is established in β stress space. The GUST introduces a weighted parameter χ to correct the UST.
2. The predictions of strength on σ1-σ2 plane under the true triaxial test condition for five types of rock are compared among the 3D H-B, the USC and the GUST. The USC effectively reflects the influence of the intermediate principal stress on rock strength and accurately predicts rock strength under both triaxial tension and triaxial compression.
3. The predicted strength on meridian plane under the triaxial test condition for eight types of rock are compared between the original H-B criterion and the USC. The results indicate that the USC criterion is applicable to various types of rock and has a higher degree of accuracy.
4. The difficulty level of the method for determining the parameters among these three criteria is similar for various types of rock. The USC are appropriate for more types of geomaterials.
Author contributions
Qin Yang: Writing – original draft, Validation, Investigation. Mingming He: Writing – original draft, Validation, Investigation, Data curation. Panagiotis G. Asteris: Writing – review & editing, Investigation, Funding acquisition. Qin Zhao: Writing – review & editing, Formal analysis, Data curation.
Acknowledgments
This paper is supported by the National Natural Science Foundation of China (Grants No. 42177158 and U2368203). The authors would like to thank Prof. Yangping Yao, Prof. Dechun Lu and Prof. Jingqi Huang for their help and suggestions in preparing the manuscript.
Conflicts of Interest
All the authors claim that the manuscript is completely original. The authors also declare no conflict of interest.