Abstract
A rock failure criterion is very important for the prediction of the failure of rocks or rock masses in rock mechanics and engineering. Least squares support vector machines (LSSVM) are a powerful tool for addressing complex nonlinear problems. This paper describes a LSSVM-based rock failure criterion for analyzing the deformation of a circular tunnel under different in situ stresses without assuming a function form. First, LSSVM was used to represent the nonlinear relationship between the mechanical properties of rock and the failure behavior of the rock in order to construct a rock failure criterion based on experimental data. Then, this was used in a hypothetical numerical analysis of a circular tunnel to analyze the mechanical behavior of the rock mass surrounding the tunnel. The Mohr-Coulomb and Hoek-Brown failure criteria were also used to analyze the same case, and the results were compared; these clearly indicate that LSSVM can be used to establish a rock failure criterion and to predict the failure of a rock mass during excavation of a circular tunnel.
1. Introduction
The failure of rock has been one of the most active research areas in the engineering of rock structures since the establishment of the International Society for Rock Mechanics (ISRM) in 1962. Ulusay and Hudson gave a detailed description of the importance of the ability to predict rock failure in practical engineering [1]. Rock failure criteria include the theoretical basis of methods of predicting when and how rock materials fail under the action of external loads and are therefore essential for optimization of engineering design and construction in rock masses. Over the past several decades, a number of failure criteria have been proposed (e.g., Mohr-Coulomb, Hoek-Brown, and Griffith’s and Drucker-Prager criteria) and developed to predict the failure of rock [2–11]. Recently, ISRM suggested methods have been published as criteria for the prediction of rock failure [12, 13]. Both the Mohr-Coulomb and Hoek-Brown criteria are generally considered to be reliable predictors of rock failure. In these two models, the major principal stress is seen as function of the uniaxial compressive strength (), the minor principal stress (), and some constants, that is, . These constants can often be determined by laboratory experiment, in situ tests, and back-analysis.
Empirical failure criteria established under specific conditions such as rock type, stress state, or functional form are not generic and cannot be adopted uncritically for any set of conditions [10] but must be established experimentally for any given set of rock formation and stress conditions.
For this reason, artificial neural network (ANN) computational models are widely used to establish rock failure criteria for geotechnical engineering problems. This is due to the fact that ANN perform function approximation without requiring additional assumptions of the function form. More detailed descriptions of the ANN-based failure criterion method have been reported by, for example, Meulenkamp and Grima, Singh et al., Çanakci and Pala, Tiryaki, Zorlu et al., and Rafiai and Jafari [14–20]. However, a disadvantage of ANN for establishing rock failure criteria is its slow convergence and overfitting. To overcome these limitations, in the present study a LSSVM-based rock failure criterion was developed based on experimental data. Support vector machines (SVM) have been widely used in geotechnical engineering for nonlinear mapping [21–23].
This paper is organized as follows. Section 2 focuses on the formulation of LSSVM. A LSSVM-based rock failure criterion is described in detail in Section 3. In Section 4, a hypothetical numerical experiment is used to verify the LSSVM-based rock failure criterion by comparing with earlier failure criteria. Finally, conclusions are given in Section 5.
2. Least-Squares Support Vector Machine (LSSVM)
LSSVM models are an alternative formulation of SVM regression proposed by Suykens and Vandewalle [24]. Suppose there is a given training set of data points with input data and output , where is an -dimensional vector space and is a one-dimensional vector space; LSSVM models in feature space are given by where the nonlinear mapping represents the input data in a higher-dimensional feature space; ; ; is an adjustable weight vector; and is the scalar threshold. For function estimation in LSSVM, the following optimization problem is formulated: where is the regularization parameter for determining the trade-off between the fitting error minimization and smoothness, and is an error variable.
The Lagrangian for the above optimization problem (2) is where is the Lagrange multiplier. The conditions for optimality are given by without regard to and ; the solution is given by the following set of linear equations: where , , , and Mercer’s theorem is applied within the matrix, , , where is the kernel function, and the matrix is invertible by choosing and can be written as Then the analytical solution of and is given by The LSSVM model can then be expressed as In (8), is the kernel function, composed of the following three forms.(1)Polynomial kernel: (2)Radial kernel: (3)Sigmoidal kernel: The brief procedure LSSVM is the following.
Step 1. Build the training samples based on rock mechanical experiment.
Step 2. Get (6) according to the LSSVM algorithm.
Step 3. Get the analytical solution of and in (7).
Step 4. Get the LSSVM model (8).
3. LSSVM-Based Rock Failure Criteria
In order to better predict the behavior failure of rock masses, LSSVM was used in this work to represent the nonlinear relationship between rock properties and rock failure and thus to establish a LSSVM-based rock failure criterion. Similar to traditional failure criteria in form, LSSVM-based rock failure criteria can be defined as where and are the major and minor principal stresses, respectively. These can be obtained by laboratory experiment and field testing. The LSSVM-based rock failure criterion can be rewritten in the form: where is the series of experimental data and and are calculated using the algorithm introduced in Section 2. Once a LSSVM-based rock failure criterion is established, it can be used in numerical analysis to predict rock failure. In this work, Excel VBA software was used to execute the code of the LSSVM-based rock failure criterion.
To implement the LSSVM-based failure criteria in a numerical model, it is useful to express it as an equivalent Mohr-Coulomb criterion with variable parameters. The parameters of the equivalent criterion are obtained as The RBF kernel function (Equation (10)) was adopted: Substituting (16) into (15) gives the following form: and the equivalent Mohr-Coulomb criterion may be calculated from the following equations: where and are the instantaneous values of cohesion and friction angle, respectively.
LSSVM-based rock failure criterion presents the mechanical character of rock through combining L-SVM and rock failure criterion. The brief procedure is showed in Figure 1.

4. Numerical Experiment
4.1. Matching the Experimental Data
In this work, the RBF kernel function was adopted for pattern analysis or recognition, with the following parameters: . The value of is 14.5061. The values of are listed in Table 1. As in any other use of LSSVM, the SVM must be trained and tested; 98 data sets from triaxial compression tests, consisting of 79 training samples and 19 testing samples, were used as the experimental data in this study (seen in Table 1). Once the performance of the SVM model was satisfactory, it was used to represent the nonlinear relationship between rock properties and the mechanical behavior of the rock to establish the LSSVM-based rock failure criterion.
As shown in Figures 2 and 3, the predicted LSSVM values and experimental values of were almost identical, indicating that the LSSVM-based failure criterion model effectively simulated the failure behavior. Figure 4 shows a comparison between the failure envelopes determined by the LSSVM, Mohr-Coulomb, and Hoek-Brown models and those from the experimental data. It can be seen from Figure 4 that predicted by the LSSVM model agreed well both with the experimental data and the Hoek-Brown model at large stresses but differed a little from the Mohr-Coulomb model; however, differed a little from the Hoek-Brown model at small stresses. These results indicate that the LSSVM-based rock failure criterion trained on experimental data predicts the failure behavior of rock with reasonable accuracy.



4.2. Numerical Analysis Using the LSSVM-Based Rock Failure Criterion
To verify the feasibility of the LSSVM-based rock failure criterion in numerical analysis, it was combined with FLAC3D modeling code to simulate the failure behavior of a 5 m radius circular rock tunnel as a hypothetical numerical case. An initial in situ stress and gravity was presupposed in the numerical model, and the experimental data from Section 4.1 was used. The properties of the rock mass listed in Figure 5 are based on the Mohr-Coulomb model. The numerical model is built using FLAC3D together with the algorithm in Section 3. The horizontal displacements are almost the same as the value of Mohr-Coulomb by FLAC3D (seen in Figure 6). The stress of surrounding rock is in well agreement with the law of Mohr-Coulomb by FLAC3D (seen in Figure 7). The horizontal displacements and the horizontal, vertical, and shear stresses induced in the surrounding rock by the excavation show that the LSSVM model results were almost identical to those obtained by the Mohr-Coulomb model. The induced horizontal and vertical stresses calculated by the two models are also shown as contour plots in Figures 8 and 9. It shows the mechanical character of rock was presented by LSSVM-based rock failure criteria. Overall, the results show that the LSSVM-based rock failure criteria present well the mechanical behavior and character and can be used for numerical analysis.



(a)

(b)

(c)

(a) Mohr-Coulomb model

(b) LSSVM model

(a) Mohr-Coulomb model

(b) LSSVM model
To verify the LSSVM-based rock failure criterion for different initial in situ stress states, the circular tunnel was investigated at different vertical stress , 10, and 15 MPa and corresponding horizontal stress in each case. The maximum error and maximum relative error of horizontal displacements are about 0.08 and 2.8%, respectively. The maximum relative error of vertical stress and shear stress is about 18% at the beginning of excavation (seen Figure 11(b)) and then it will be less than 5% with the excavation. Using the LSSVM and Mohr-Coulomb models, the calculated horizontal displacements in the rock surrounding the tunnel are shown in Figure 10 and the major principal stresses are shown in Figure 11. Those show it is feasible to combine LSSVM rock failure criterion with numerical analysis.


(a)

(b)

(c)
The displacements for both models agree well for the three initial in situ stress states (Figure 10), and the principal stresses are also similar for both models (Figure 11), with and showing small differences near the wall of tunnel. These results generally reflect the above results (in which MPa).
The LSSVM and Hoek-Brown models were also compared. Figures 12–14 show the results for the same initial stress states as above, and it is seen that displacements in the surrounding rock agreed well with both the Hoek-Brown and Mohr-Coulomb results. Some small differences are evident in and compared to the Mohr-Coulomb model, but the results are consistent with the Hoek-Brown model because, as shown in Figure 4, the LSSVM and Hoek-Brown failure envelopes are similar. Thus, the proposed failure criterion can be used in numerical analysis to effectively reflect the mechanical behavior of rock.

(a) Displacement

(b)

(c)

(d)

(a) Displacement

(b)

(c)

(d)

(a) Displacement

(b)

(c)

(d)
5. Conclusions
The results of comparisons between the proposed LSSVM-based rock failure criterion and the Mohr-Coulomb and Hoek-Brown criteria and experimental data demonstrated clearly that LSSVM provided an effective rock failure criterion for the purpose of numerical analysis. Comparisons of the displacements in the rock surrounding a circular tunnel from the LSSVM, Mohr-Coulomb and Hoek-Brown models, and experimental data showed that the LSSVM model mapped the nonlinear relationship between the mechanical properties of the rock and its failure behavior.
List of Symbols
| Number of samples | |
| : | -dimensional vector space |
| : | One-dimensional vector space |
| : | Input and output of training samples |
| : | Nonlinear mapping in a higher-dimensional feature space |
| : | Adjustable weight vector |
| : | Scalar threshold |
| : | Regularization parameter |
| : | Error variable |
| : | Lagrange multiplier |
| : | Kernel function |
| : | Matrix of kernel function |
| : | matrix |
| : | Major principal stress |
| : | Minor principal stress |
| : | Vertical in situ stress |
| : | Horizontal in situ stress |
| : | Horizontal stress of surrounding rock mass |
| : | Vertical stress of surrounding rock mass |
| : | Shear stress |
| : | Experimental data |
| : | Instantaneous values of cohesion |
| : | Instantaneous values of friction angle |
| : | LSSVM-based rock failure criteria |
| : | Parameter of RBF kernel function. |
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
Support by the Program for Innovative Research Team (in Science and Technology) in University of Henan Province (no. 15IRTSTHN029) and National Fund of Science in China (nos. 51104057, 41172244) are gratefully acknowledged.