Abstract
This study mainly establishes a novel intelligent assessment model for peak shear strength of rock joints based on the relevance vector machine (RVM). RVM is a state-of-the-art soft computing technique that has been rarely utilized in joint shear strength assessment. To establish the hybrid intelligent model, three-dimensional scanning tests and direct shear tests on 36 granite joint specimens were conducted. The peak shear strength ratio () is perceived as an explanation of four types of influencing factors, including joint surface roughness, strength of rock material, basic friction angle, and normal stress. In particular, the compressive strength and tensile strength of rock material are first considered together. A total of 36 experimental data were used in this study to train the RVM model to predict the peak shear strength of rock joints. The performance of the RVM model was assessed using the direct shear test data of rock joints collected from previous researches. Four different kinds of kernel functions were adopted to obtain the optimal model. Results show that the proposed model is significantly efficient in predicting the peak shear strength of rock joints. The proposed model is also a promising tool for peak shear strength of rock joints and provides a new research approach to research the mechanical properties of rock joints.
1. Introduction
The deformation capacity, strength, and stability of fractured rock mass are strongly influenced by the shear behavior of joints. Therefore, the accurate determination of shear strength is the key to design the safety structure in rock mass or rock masses [1]. Over the past decades, empirical [2–9], semi-theoretical [10], and theoretical methods [11] have been proposed to determine the shear strength of rock joints. Various factors, such as rock type, joint surface roughness, joint size, and infilling materials, exhibit a wide variation of joint shear strength [12, 13]. Considering their strong randomness, fuzziness, and uncertainty, the peak shear strength of rock joints has a complex nonlinear relationship with the aforementioned factors [14, 15]. The traditional methods demonstrate strong usability and dependability in practical engineering, but large computations and complex algorithms limit the use of these methods. Thus, the deterministic analysis methods in assessing the joint shear strength remain difficult. In recent years, data mining and artificial intelligence technologies, such as artificial neural network (ANN) and support vector machine (SVM), have been widely used to estimate the mechanical properties of rocks, such as the compressive strength and tensile strength [16–21]. However, ANN has several disadvantages, such as poor generalizing performance, large training data requirement, slow convergence velocity, and over-fitting problem [22–24]. SVM can effectively solve the aforementioned problems and exhibit excellent generalization capability and high regression accuracy [25, 26]. Despite many advantages, SVM suffers from several limitations. The output of an SVM represents decisions rather than posterior probabilities [27]. In addition, the kernel function of SVM must fulfill Mercer’s condition. The relevance vector machine (RVM) is a probabilistic model based on Bayesian theory for regression and classification; it has high generalization capability, sparse model structure, and low computation complexity, avoiding the principal limitations of SVM [28]. Especially for small samples, RVM can map the limited input data into a high-dimensional feature space and finally use several relevance vectors to deal with the dimensionality problems. The primary purpose of this study is to introduce a reliable estimation method of peak shear strength of rock joints on the basis of RVM. To this end, 36 rock joints under different stress were designed and tested in laboratory. The rest of this article is organized as follows: Section 2 presents the theoretical background of RVM. Section 3 introduces the description of the experimental procedure and results. Section 4 describes the establishment of the RVM model. Section 5 presents the comparison and discussion. The final section presents the conclusions.
2. RVM Method
RVM is a probabilistic learning model based on the Bayesian framework to solve complex regression and classification problems [28]. It generates a mapping between the targets and the corresponding input vectors .where N is the total sample number, n is the nth sample number, y is the output variable, is the weight of the nth sample, is the kernel functions, and is a mean-zero white noise process with variance .
Assuming an independent distribution for , the likelihood of the entire dataset can be expressed as follows:where , , and .
Direct maximum-likelihood estimation of from equation (2) may suffer from overfitting [29]. Therefore, a zero-mean Gaussian prior distribution over every weight with the diagonal covariance is adopted and described as follows:
The posterior distribution for the weights can then be calculated using Bayes’ law based on the Gamma prior distributions [30]:
Integrating the weight parameters, the marginal distribution of hyperparameters and can be obtained:where , , and is a unit vector. Taking its partial derivatives to equation (5) by iteration, the values of and can be obtained.
3. Experimental Procedure and Results
3.1. Preparation of Rock Samples
To investigate the effect of wide-range normal stress on peak shear strength, rectangular granite joints were selected for the direct shear tests because of their high strength. Rectangular blocks of rock specimens with dimension of 100 mm × 100 mm × 200 mm were labeled as SJ1–SJ36. Brazilian splitting tests were conducted by using a TSF-600 shear-flow coupled test system for rock joints at Tongji University. The granite joint surface dimensions were approximately 100 mm wide and 200 mm long (Figure 1). After testing, no apparent defects were found in the specimens. is the uniaxial compressive strength, is the tensile strength, is Young’s modulus, is the rock density, and is Poisson’s ratio, and basic friction angle was determined using the suggested methods of the International Society for Rock Mechanics (ISRM) [31], and the results are listed in Table 1.

(a)

(b)
3.2. Surface Topography Measurement
To investigate the influence of the three-dimensional (3D) surface roughness on the mechanical properties of rock joints quantitatively, the joint surface topography of measuring each shear test sample block before using a noncontact direct shear test, 3D stereo-topometric measurement system (positioning accuracy: 20 μm) was developed by Xia et al. [32]. Before the estimation of joint surface roughness, the 3D surface topography was reconstructed on the basis of the interpolation of the scanned point clouds. Tatone and Grasselli indicated that the digitized surface topographies of rock joints with sampling interval <0.55 mm using linear interpolation technique can capture microfeatures of surface morphology [33]. In this study, we set the sampling intervals in length (L) and width (W) directions to be 0.3 mm. Taking upper surface topographies of four joint samples as examples, the digitized surface topographies of rock joints (the height magnified to 3 times its normal size) are shown in Figure 2.

(a)

(b)

(c)

(d)
3.3. Test Procedure and Results
This test was designed to obtain an extensive peak shear strength model of rock joints. We conduct the direct shear tests on the test samples under nine normal stress levels, namely, 1.0, 2.0, 3.0, 5.0, 8.0, 10.0, 12.0, 15.0, and 20.0 MPa. Four direct tests were performed on the test specimen at each level of normal stress. The values of the peak shear strength for each test specimen are listed in Table 2. The rate of shear displacement was set to be 0.5 mm/min, and the maximum displacement was 20 mm which is 10% of the joint length.
4. Establishment of a Peak Shear Strength Assessment Model Based on RVM
4.1. Selecting the Influencing Factors
According to the study by Barton and Choubey [34], the joint surface roughness, rock material strength, basic friction angle, and normal stress exhibit a wide variation of shear behavior. Therefore, normal stress, joint surface roughness, rock strength, and basic friction angle were selected in this work to assess the peak shear strength.
Various parameters, including empirical, statistical, and fractal, have been proposed to quantify the surface roughness of rock joints. Numerous experimental tests demonstrated that partial contact between the lower and upper joint surface was located in the steepest zones facing the shear direction during the shearing [35–37]. Grasselli et al. [38] established a relationship between total potential contact area ratio and the corresponding apparent dip angle , as shown in the following equation:where is the maximum potential contact area under apparent dip angle at 0°; is the maximum apparent dip angle in the shear direction; and C is a dimensionless roughness parameter characterizing the distribution of the apparent dip angle over the surface, which is calculated using equation (6). Considering the possible range of , the surface roughness of rock joints can be measured by the roughness metric [33]. varies significantly in different directions; thus, the 3D roughness metric and are both included to fully describe the roughness anisotropy. The 3D morphology parameters of the granite joints, namely, , , and C, were calculated using the measured data from 36 surfaces on the basis of equation (6). The values of and are listed in Table 3. Figure 3 shows the fitting curves of SJ3, SJ12, SJ21, and SJ36 for calculating the roughness parameter C.

(a)

(b)

(c)

(d)
Joint wall compressive strength (JCS) is commonly used for representing the influence of rock material strength on the shear behavior of rock joints [3, 8, 25, 39–41]. For the fresh rock, JCS can be replaced by the uniaxial compressive strength . Based on experimental observation and numerical simulations, several researchers have also found that the failure of asperities in rock joints under direct shear loading was possibly attributed to shear-induced tensile fracture [5, 6, 42–44]. However, the failure mechanism of asperities has not been clearly stated, and controversies still exist. Therefore, uniaxial compressive strength and tensile strength of the rock are both selected in this study to describe the effect of rock material strength on the peak shear strength of rock joints.
Basic friction angle is crucial to determine the shear strength of rock joints. Barton suggested a fundamental form of joint shear strength [3]:where is the normal stress, is the peak dilation angle, and is the shear component angle. On the basis of equation (7), various 2D or 3D morphology parameters have been used for representing the sum of peak dilation angle and shear component angle . No uniform approach for obtaining the basic angle is currently available. In this work, laboratory direct shear testing method for two smooth joint surfaces is adopted to obtain the basic angle and the extensive shear strength model, as shown in Table 1.
4.2. Training and Testing of the RVM Model
In this study, the aforementioned original influencing factors, namely, , , and , can be dimensionless preprocessed as two factors: and . Therefore, the five influencing factors, which include , , , , and , form the input vector of the regression model. All these index values can be measured in the laboratory (see Tables 1–3). Mapping the above five influencing factors to the five-dimensional feature space according to the definition of RVM, the appropriate kernel function was chosen to simplify the nonlinear transformation. The linear kernel function was first used due to its parameter-free and rapid calculation as follows:where are the training dataset of five influencing factors. The results of direct shear tests on 36 rocks in this study were used to cultivate and advance the RVM model. The results of direct shear tests by Grasselli and Egger [5] and Yang et al. [8] were used for testing the RVM model, in which a total of 50 tests were conducted on five different material types of rock joints (granite, sandstone, limestone, marble, and serpentinite). The 3D morphology parameters of 50 groups of rock joints are represented in Table 4. After the analysis was completed, the quality of the regression model was estimated by the average error formula and correlation coefficient [45]:where is the average estimation error, n is the sample size, is the measured peak shear strength, is the estimated peak shear strength, is the average value of the measured peak shear strength, is the comprehensive estimation error of the RVM model, and is the correlation coefficient of the train and test model; train_n is the number of the training sample, and test_n is the number of the testing sample. In this work, train_n and test_n are 36 and 50, respectively. The average estimation errors of the RVM model for the training and testing data are 6.46% and 21.12%, respectively. The correlation coefficient of the RVM model for the training and testing data is 0.977 and 0.769, respectively. The comprehensive estimation error and correlation coefficient of the RVM model for all the data are 14.98% and 0.856, respectively. This finding indicates that the RVM model exhibits good training accuracy and high generalization capacity. Therefore, the RVM model is feasible and effective for estimating the peak shear strength of rock joints and can be used in practical applications.
5. Comparison and Discussion
5.1. RVM with Different Kernel Functions
The inner product of kernel functions is usually substituted for the complex calculation of high-dimensional feature space to deal with the dimensionality problem under high-dimensional feature space. The nonlinear shining mode from the raw feature space to high-dimensional space relays on the type of kernel functions. The approximate kernel function and model parameters have a considerable effect on the model stability and regression precision [46].
To further enhance the accuracy of the proposed model, three other common kernel functions were introduced: polynomial kernel function, Gaussian kernel function, and sigmoid kernel function.where a, b, c, d, k, and p are kernel function parameters.
In the RVM model, the kernel function parameters must be defined. The most relevant papers which include the determination of the optimal kernel function parameters for RVM is lacking. Several different values of kernel function parameters were selected as determining which one provides better regression accuracy. The four kernel function parameters and the accompanying average estimation errors and correlation coefficients are represented in Table 5. Synthetically considering the regression performance and the kernel function parameter number, a peak shear strength of RVM model based on a Gaussian kernel function is established. The optimal value of the kernel function parameter d is 52.43.
5.2. Comparison with Existing Criteria
To certify the rationality of the proposed model, the two models developed by Grasselli and Egger [5] and Yang et al. [8] have been chosen for comparison.
The 3D morphology parameters of the 50 rock joints are represented in Table 4. The estimated results using equations (16) and (17), the RVM model results, and the measured peak shear strength are listed in Table 4. The prediction results of the peak shear strength of rock joints are compared with the measured data, as shown in Figures 4 and 5.


According to equation (9), the mean errors of the estimated value results calculated using Grasselli’s model, Yang’s model, and RVM model are 16.3%, 8.9%, and 15.4%, respectively. Based on equation (11), the correlation values of the estimated value results calculated using Grasselli’s model, Yang’s model, and RVM model are 0.889, 0.940, and 0.938, respectively. From the mean error and correlation coefficient analysis, the average estimation error values of the three models are less than 20% and the correlation coefficients are greater than 0.85, which indicate that the three models were quite efficient to be used to predict the shear strength of rock joints. Moreover, the proposed model is more accurate than that of Grasselli’s model and has almost equivalent precision to that of Yang’s model. After analysis of the above average errors, the existent limitations of the three models should be pointed out. In Grasselli’s and Yang’s models, the function forms are complex. According to Grasselli’s model (i.e., equation (14)), the ratio () was taken as a measure of joint surface roughness, but cannot be zero ( for saw-tooth joints). Yang’s model is not fully explained on how to use the two parameters and in combination, and is not even considered. Furthermore, Grasselli and Egger [5] and Yang et al. [8] deemed that the actual contact area was closely related with and , respectively. Compared with Grasselli’s and Yang’s models, the proposed model was less complex for large samples and comprehensively considered many factors; in particular, the effect of compressive strength and tensile strength on the peak shear strength of rock joints was first investigated. However, the determination of kernel function parameter needs to be further improved and optimized.
In this study, the investigation is limited to the shear behavior of unfilled joints. Influencing factors, such as size of joint (scale effect), sampling interval, joint matching coefficient, and joint filling, should be further investigated. Although the proposed model has limitations, it is believed to be a novel research idea for assessing peak shear strength of rock joints.
6. Conclusions
The hybrid probabilistic regression assessment algorithm RVM was accustomed to identify the peak shear strength of the rock joint, with Gaussian kernel function being introduced to improve the prediction accuracy. The direct shear tests on joint specimens were conducted to train the new model, and its rationality on assessing the peak shear strength was also validated. Some concluding remarks are presented as follows:(1)A novel intelligent peak shear strength model combined with RVM was supposed to assess the peak shear strength of rock joints, which can comprehensively consider the elements influencing the peak shear strength, such as normal stress, joint surface roughness, rock strength, and basic friction angle. The 36 natural granite joint specimens were accustomed to examine the suggested model, and the result reflects that the RVM model exhibits good training accuracy and high generalization capacity.(2)Four types of kernel functions were investigated for their performance in regression accuracy in the RVM model. According to the regression performance and the kernel function parameter number, the Gaussian kernel function was determined to be the better one combined with the RVM model. The optimal value of the kernel function parameter d is 52.43.(3)The peak shear model of Grasselli’s and Yang’s was chosen for comparison with the proposed model. Based on the 3D morphology parameters of 50 rock joints, the prediction results show that the proposed model can meet the accuracy requirements. In addition, the new model has less complexity and can take into account the effects of compressive strength and tensile strength on the peak shear strength of rock joints.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
The study was funded by the Natural Science Foundation of Zhejiang Province (no. LY18D020003) and the National Natural Science Foundation of China (nos. 41327001, 41572299, 41427802). This support was gratefully acknowledged.
Supplementary Materials
1. The file named topograph_3D is the code for the 3D parameter of . 2. The file named Data contains two sheets, one of which is the prediction of the shear strength, and the second is the comparison of the kernel functions. 3. The file named 3D roughness parameter is the fitting curves of SJ3, SJ12, SJ21, and SJ30 for calculating the roughness parameter C. 4. The file named Comparation is the display of the shear strength prediction comparison with the model of Grasselli’s and Yang’s. (Supplementary Materials)