Next Article in Journal
A Review on Methods and Systems for Remote Collaboration
Previous Article in Journal
Assessing Nitrate Contamination Risks in Groundwater: A Machine Learning Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hoek-Brown Failure Criterion-Based Creep Constitutive Model and BP Neural Network Parameter Inversion for Soft Surrounding Rock Mass of Tunnels

State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, College of Environment and Civil Engineering, Chengdu University of Technology, Chengdu 610059, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2021, 11(21), 10033; https://doi.org/10.3390/app112110033
Submission received: 16 August 2021 / Revised: 12 October 2021 / Accepted: 19 October 2021 / Published: 26 October 2021
(This article belongs to the Topic Artificial Intelligence (AI) Applied in Civil Engineering)

Abstract

:
This paper summarizes the main factors affecting the large deformation of soft rock tunnels, including the lithology combination, weathering effect, and underground water status, by reviewing the typical cases of largely-deformed soft rock tunnels. The engineering geological properties of the rock mass were quantified using the rock mass block index (RBI) and the absolute weathering index (AWI) to calculate the geological strength index (GSI). Then, the long-term strength σ r and the elastic modulus E 0 of the rock mass were calculated according to the Hoek–Brown failure criterion and substituted into the creep constitutive model based on the Nashihara model. Finally, the creep parameters of the surrounding rock mass of the Ganbao tunnel were inverted and validated by integrating the on-site monitoring and BP neural network. The inversion results were consistent with the measured convergence during monitoring and satisfied the engineering requirements of accuracy. The method proposed in this paper can be used to invert the geological parameters of the surrounding rock mass for a certain point, which can provide important mechanical parameters for the design and construction of tunnels, and ensure the stability of the surrounding rock mass during the period of construction and the safety of the lining structure during operation.

1. Introduction

Tunnel projects in western China often encounter soft rocks with well-developed bedding, such as carbonaceous phyllite, sericite phyllite, schist, carbonaceous slate, sandy slate, and carbonaceous shale. Under high in situ stress, the laminated soft surrounding rock masses are vulnerable to large and rapid deformation and local destruction [1,2,3]. In these conditions, tunnels frequently suffer from large deformation hazards.
Laminated rock mass, frequently seen in engineering practice and presenting oriented grouped bedding, has more heterogeneous mechanical properties than normal rock mass. Many researchers have developed constitutive models for laminated surrounding rock masses. For instance, Jia et al. [4] applied the microscopic element method with the constitutive model based on damage mechanics and statistical theory to the simulation of rock tunnel stability using the finite element method (FEM). Li et al. [5] developed a three-dimensional creep constitutive model for transversely isotropic rock mass, based on the Burgers viscoelastic model. Li et al. [6] proposed three basic creep patterns of shale and a general methodology for developing the anisotropic creep model. These constitutive models have been used to probe the deformation and failure characteristics of laminated rock masses with varied dip angles. However, research on the constitutive theory of laminated surrounding rock masses in tunnels is fairly limited and cannot provide the theoretical guidance for the engineering design of actual tunnels. In terms of the mechanisms of laminated rock mass deformation, previous studies have mainly investigated patterns of stability [7], deformation [8,9,10], mechanical behavior [11], and mechanical properties [12]. The main triggers for large tunnel deformation through highly-dipped laminated soft rock under high in situ stress are high structural stress, unfavorable rock occurrences, and low rock strength [13]. On-site tunnel monitoring has shown that the deformation–failure zone of the surrounding rock mass of laminated soft rock tunnels is concentrated along the direction perpendicular to the rock bedding, instead of the direction of the maximum principal stress; the local deformation of laminated soft rock tunnels is affected by the topography, rock mass structure, and in situ stress [14]. Four mechanisms have been proposed to explain large tunnel deformation, namely, along-bedding sliding, flexural deformation, toppling deformation, and plastic extrusion. For each mechanism, a specific anchoring support plan should be developed [15]. Numerical simulation has been performed to comprehensively analyze the stress environment, the deviator stress distribution [16,17,18,19], and the plastic zone distribution [20,21] of the surrounding rock mass during tunnel excavation. Such simulations, incorporating the rock characteristics of deterioration, stabilization, and accelerated creep [22], reveal the characteristics of the surrounding rock mass failure with varied surrounding rock mass strengths, buried depths, and roof strengths [23], which provide guiding values to support the design of the analogous tunnels.
Due to the complexity of geology, geotechnical materials are typically characterized by discontinuity and heterogeneity, and the rheological parameters of the actual surrounding rock mass of the tunnel are often hard to measure. There are two main methods to obtain creep parameters: one is to calculate creep parameters through laboratory testing of rock combined with a creep constitutive model, but there is a large deviation between the laboratory results for rock and on-site rock mass parameters. Second, to obtain the mechanical parameters of the surrounding rock mass through displacement back analysis based on monitoring data [24]. Wenzheng Cao [25] proposed a novel back analysis program based on a BP neural network, which can realize automatic correction and the adjustment of parameters and adapt to most tunnel projects. Qingdong Wu [26] introduced a support vector machine (SVM) and an artificial neural network (ANN) to predict tunnel surrounding rock mass displacement, and compared and analyzed the results of the two methods. Xianghui Deng [27] established a tunnel risk assessment model by combining a fuzzy method with a BP neural network based on historical data from 50 tunnels. Numerical simulation has been widely combined with back-propagation (BP) neural network analysis to calculate the physical and mechanical parameters of surrounding rock mass, which has achieved good application performance [28,29,30].
This paper proposed a solution for obtaining the creep parameters of the surrounding rock mass in highway tunnels after excavation. First, the long-term strength of the engineering rock mass was obtained using the Hoek–Brown failure criterion, and a creep constitutive model for the rock mass was developed. Then, back-propagation analysis was performed using the on-site monitoring data to obtain the initial values of the creep parameters of the tunnel’s surrounding rock mass. Subsequently, training samples were generated via orthogonal experiments. Then, a BP neural network and numerical simulation were integrated to inversely calculate the tunnel’s surrounding rock mass creep parameters that matched the reality.

2. Factors Affecting the Deformation of Soft Rock Tunnels under High In Situ Stresses

The engineering data from some typical largely-deformed soft rock tunnels in China are summarized in Table 1, including the buried depth, formation lithology, in-situ stress, deformation magnitude, and characteristics of deformation and failure.
Surrounding rock masses that are found with large deformations are mostly soft rock with a highly developed bedding texture, such as shale, slate, and phyllite (Table 1). The bedding plane profoundly affects the strength and deformation characteristics of tunnel-surrounding rock mass and is considered the key factor determining the mechanical behavior and deformation-failure characteristics of anisotropic rock mass. In general, the compressional deformation of the surrounding rock mass of tunnels in laminated soft rock occurs mainly at the two sides along the normal direction of the weakest bedding plane. Moreover, the smaller the bedding spacing is, the worse the integrity and self-stabilizing ability of the surrounding rock mass are, which leads to increased proneness to local large deformation and thus a relatively large compressional load on the support structure and highly-deformed lining and surrounding rock mass in local areas. The factors affecting the large deformation of the surrounding rock mass can be summarized as follows:
Based on lithology, the surrounding rock mass can be grouped into two types, namely, plastic and brittle rock. The former consists of soft rock and is typically characterized by low mechanical strength, softening and swelling on wetting, and low slaking resistance, which are all unfavorable for the load-bearing and stability of the tunnel’s surrounding rock mass. On the contrary, the latter is of hard rock and generally has a high load-bearing capacity, low softening and swelling tendencies on wetting, and high slaking resistance, which are favorable for the load-bearing and stability of the tunnel’s surrounding rock mass. The lithology combination can be generally summarized into four cases, namely, the consistent lithology, hard interbeds (in much softer surrounding rock mass), soft interbeds (in much harder surrounding rock mass), and soft–hard alternating lithology. Different lithology combinations are associated with varied deformation strengths. A soft interbed in the surrounding rock mass is equivalent to a weak part in the surrounding rock mass. In most cases, deformation-failure occurs first in the soft interbed during tunnel excavation. If a hard interbed exists in the surrounding rock mass, it will often suppress the rock deformation, due to its more rigid mechanical properties than those of the surrounding rock mass. Therefore, the deformation magnitude is often relatively small at the hard interbed.
The weathering effect is a common geological phenomenon in nature, which can weaken the cementation between rock particles, form damage cracks, reduce surface roughness, and worsen the physical and mechanical properties of rock. Rapid weathering of soft rock in tunnels will reduce the stability of the surrounding rock mass and affect the safety of the tunnel construction. After the tunnel excavation, the weak layered surrounding rock mass will form internal cracks due to the unloading effect; then, an effective seepage channels are formed, groundwater seeps through these channels, reducing the effective stresses and thus the shear strength along discontinuities and therefore the strength of the rock mass. Therefore, the presence of underground water makes the rock mass more prone to large deformation. The varied underground water seepage pathways or inconsistent thickness of the broken rock zone in the surrounding rock mass after excavation may result in local large deformation of the surrounding rock mass.
Due to the complex geological conditions outlined above, obvious asymmetry is observed for the deformation occurring at the two sides of the central axis, crown, and inverted arch, and local deformation is a common phenomenon for the practice of tunnel engineering. Therefore, it is necessary to quantitatively analyze the engineering geology of each part of the tunnel’s surrounding rock mass and accurately invert the creep parameters at each local point from the on-site deformation monitoring data, so as to guide the design and construction of tunnels.

3. Rock Strength and Creep Constitutive Model Based on the Hoek–Brown Failure Criterion

The classic rock mass failure criterion considers rock as a continuous homogeneous medium, which greatly simplifies the reality. A constitutive model for the rock mass with multiple structural planes was developed via superimposition of the constitutive equations of the rock mass with a single joint, which cannot accurately capture the complex mechanical characteristics of the tiny fissures and multiple groups of structural planes in the rock mass [31,32]. Hence, the semi-empirical semi-theoretical approach was used to characterize the laminated rock mass in engineering practice.
In 1980, E. Hoek and E. T. Brown derived the correlation expression among the limiting principal stresses for the failure of rock (mass) (namely, the Hoek–Brown failure criterion) via a trial-and-error process, based on the statistical analysis of massive data from rock triaxial tests and field testing of rock mass, and the Griffith theory. Then, E. Hoek further proposed the generalized Hoek–Brown empirical failure criterion. The Hoek–Brown empirical failure criterion can reflect the effects of various factors of the rock mass—such as the rock strength, the number of structural planes, and the in situ stress—on the rock mass strength, and overcome the disadvantages of the conventional theoretical equation. The generalized Hoek–Brown failure criterion can be expressed as below:
σ r = σ 3 + σ c ( m b σ 3 σ c + s ) a
where σ r is the maximum principal stress when the rock fails, namely, the rock strength; σ 3 is the minimum principal stress when the rock fails; σ c is the uniaxial compressive strength of the intact rock; m b , s, and a are all empirical parameters; m b represents the hardness of the rock, with a value of 0.0000001–25 (0.0000001 for severely-disturbed rock and 25 for intact hard rock); s represents the broken degree of the rock mass, with a value of 0–1 (zero for completely broken rock and one for intact rock); a is related to the rock mass quality.
For rock mass with good quality, its strength characteristics are mainly controlled by the rock particle strength, due to the tight packing of rock particles. The restricted Hoek–Brown empirical failure criterion is more applicable to this case and a = 0.5. On the contrary, for rock mass with poor quality, the packing of fragments in the rock mass is loosened by shearing or weathering, which results in loss of the tensile strength of the rock mass (zero cohesion). Under such circumstances, if there is no confinement, the rock mass will collapse and a should be assigned to other values. The expressions of the parameters in Equation (1) are listed below:
m b = m i exp ( G S I 100 28 14 D )
s = exp ( G S I 100 9 3 D )
a = 1 2 + 1 6 ( e G S I / 15 e 20 / 3 )
0 < σ 3 < σ c 4
where D is the disturbance factor (zero for the undisturbed rock mass, and one for completely disturbed rock mass). mi indicates the values of constants for intact rock, which can be determined by laboratory testing (uniaxial compression test and conventional triaxial compression test). Based on laboratory data and engineering experience, E. Hoek et al. [33,34] developed a comprehensive and detailed mi value table that was able to cover a variety of rocks; the mi index of intact phyllite in this paper refers to this mi value table.
Afterwards, E. Hoek and E. J. Brown further derived the estimations of the relevant mechanical parameters of the rock mass [35], based on Equation (1):
σ cr = s σ c
σ tr = 1 2 σ c ( m b m b 2 + 4 s )
E rm = 10 5 1 D / 2 1 + exp [ ( 75 + 25 D G S I ) / 11 ]
Equation (1) shows that the GSI value is the key to determining the mechanical parameters of the Hoek–Brown failure criterion for the rock mass. However, E. Hoek proposed only a general range for the texture type and weathering conditions of the rock mass and offered no quantitative calculation method.
Rock quality designation (RQD) is the most commonly used index to describe the structural characteristics of rock mass, which is defined as the percentage of the sum of the lengths of intact core pieces longer than 0.1 m relative to the chosen length of an evaluated drill core ( R Q D = ( the   lengths   of   intact   core   pieces     10   cm   per   footage     the   length   of   an   evaluated   drill   core ( % ) ); however, the definition of RQD requires that the quality of core pieces must be “hard and sound”, meaning it can only calculate the percentage of the core with a length beyond 10 cm in the total cored length, which is not suitable for integrity evaluations when characterizing the multi-joint rock mass.
This paper introduces the rock mass block index (RBI) and the absolute weathering index (AWI) to quantify the GSI. The RBI, proposed by Hu et al. [36], can thoroughly characterize the block dimension, texture type, and structural packing of blocks for the rock mass. The RBI defines the percentages of drilling cores with the measured lengths of 3–10 cm, 10–30 cm, 30–50 cm, 50–100 cm, and >100 cm as the weights, and the value of the RBI is the sum of the weights multiplied by the corresponding coefficients.
Although the rock mass has the same RQD value, the RBI values can be different. For example, the rock mass of RQD = 90% can be 10~30 cm with a mosaic structure, or 30~50 cm with a block structure, or 50~100 cm with a block structure, or even more than 100 cm with intact structure. Under the same RQD, the larger the RBI is (from a mosaic structure to an intact structure), the more intact the rock mass is, and the RBI can be regarded as an extension of RQD. The calculation formula of RBI is shown below:
R B I = 3 C r 3 + 10 C r 10 + 30 C r 30 + 50 C r 50 + 100 C r 100
where C r 3 ,   C r 10 ,   C r 30 ,   C r 50 ,   and   C r 100 are the acquisition rates (weights in percentages) of the cores with the lengths of 3–10 cm, 10–30 cm, 30–50 cm, 50–100 cm, and >100 cm, respectively. The details of the rock mass texture representation by RBI, associated with the field description, are shown in Table 2.
The mechanical properties of rock are also related to the weathering level. Chemical weathering can alter the mineral composition of rock and thus alter its physical and mechanical performance.
Parker [37] proposed an index of weathering for silicate rocks, based on the proportions of the alkali and alkaline earth metals present. The bond strengths of the most mobile of the major elements with oxygen are used as weighting factors in the index; the proposed index is defined by the following expression:
W I = [ a ( Na ) 0.35 + a ( Mg ) 0.9 + a ( K ) 0.25 + a ( Ca ) 0.7 ] × 100
where a(X) (X = Na, Mg, K, Ca) indicates the atomic proportion of element X, defined as the atomic percentage divided by the atomic weight, and the denominator is the bond strength of element X with oxygen, which represents the stability of the element in the weathering process. From fresh rock to weathered rock, the index WI is gradually reduced.
The weathering index of rock mass A W I 0 [38] can be expressed as below:
A W I 0 = W I / W I
where WI stands for the weathering index of the weathered rock, whereas WI′ represents that of the fresh rock. Thus, a higher value of AWI indicates fresher rock with a lower level of chemical weathering and correspondingly better mechanical properties. On the contrary, a higher AWI value denotes a high level of chemical weathering. Obviously, AWI quantitatively characterizes the weathering conditions of the engineering rock mass. Based on the research on the relationship between the weathering index and the weathering level of the rock [39], the weathering level of rock based on AWI is shown in Table 3/Figure 1. In order to conveniently determine the absolute weathering index A W I 0 of the rock mass, can use the ratio of uniaxial compressive strength of fresh rock to weathered rock.
A correction factor λ is introduced to describe the effects of underground water on the weathering conditions of the jointed rock mass (Table 4). Accordingly, the rock mass weathering index, considering the effects of underground water, can be expressed as A W I = λ A W I 0 .
By integrating the quantitative indexes in Table 2, Table 3 and Table 4 for the texture and weathering conditions of the rock mass, the geological strength index (GSI) quantification method (Figure 2) can be obtained.
The widely applied Nishihara model can characterize the decelerating, steady-state, and accelerating creep phases of rock. Extensive experiments have shown that notable damage of the rock occurs only in the accelerating creep phase. The rock strength σ r , presented upon the failure of the rock, is the long-term strength for rock creep. The creep elements of the rock mass model are shown in Figure 3.
The creep equation can be expressed as below:
ε = { σ E 0 + σ E 1 [ 1 exp ( E 1 η 1 t ) ] ( σ σ r ) σ E 0 + σ E 1 [ 1 exp ( E 1 η 1 t ) ] + σ σ r η 2 t ( σ > σ r )
where E 0 is the elastic modulus; E 1 is the viscoelastic modulus; η 1 and η 2 are the viscosity coefficients; σ is the load stress; ε is the total strain; σ r is the long-term strength calculated using Equation (1) (the Hoek–Brown failure criterion); and t is the creep time.
The strength parameters of the laminated jointed rock mass at each position (Figure 4) can be obtained according to the Hoek–Brown failure criterion. Then, they are input into the numerical simulator via the creep constitutive model (Equation (12)) to calculate the deformation of the surrounding rock mass (Figure 5).

4. Field Application

The stress state of the surrounding rock mass changes after starting the excavation of the tunnel. Hence, the support design based on the surrounding rock mass parameters obtained in the previous investigation is inappropriate, in particular for a tunnel engineering design with a high buried depth, large span, and highly-developed structural planes. To provide a basis for theoretical reference and design, this paper inverts and validates the surrounding rock mass parameters by integrating the on-site monitoring with a BP neural network.
The BP neural network is a typical nonlinear algorithm, composed of the input, output, and several (one or more) hidden layers (Figure 6), and each layer has several nodes. The connection of nodes between layers is represented by the weight. The BP neural network with one hidden layer is the traditional shallow neural network, whereas that containing multiple hidden layers is the deep learning neural network (Figure 7).
The core steps for training the BP neural network are illustrated below (Figure 8), in which the solid line represents forward propagation and the dashed line represents back propagation. The forward propagation means that the data (information or signal) are first imported into the input end, delivered along the network direction, and multiplied by the corresponding weight. The products of the input data and the corresponding weights are then summed, and the results are input into the activation function for calculation. Then calculation results are delivered to the next node as the input. The calculation is performed successively until the final output is obtained.
Here we summarize the inversion workflow of the creep parameters of the surrounding rock mass, as shown below (Figure 9).
The Ganbao tunnel project of the Wenchuan–Maerkang expressway was taken as an engineering case study. It is an ultra-long left-right separated tunnel project (the left tunnel is 4777 m long and the right one is 4796 m long). The longitudinal profile of the tunnel is illustrated in Figure 10, in which the red surrounding rock mass mainly represents the sericite phyllite; the yellow part represents the sericite phyllite interbedded by carbonaceous phyllite and metasandstone; the blue part shows the carbonaceous phyllite and sericite phyllite interbedded by metasandstone; the green part represents the sericite phyllite and carbonaceous phyllite; and the black part indicates the carbonaceous phyllite. The tunnel was excavated using the drilling and blasting method. A large deformation occurs at the mile sections K141 + 540~K141 + 860 and K143 + 485~K143 + 650, which results in the circumferential cracking of multiple points along the initially supported arch and local cracking and spalling of the sidewall and tunnel crown (Figure 11). The monitoring section is located at K143 + 621 and starts monitoring after 12 h of excavation. The overburden near the monitoring section consists of collapsed debris deposits, sericite phyllite and carbonaceous phyllite (from top to bottom, Figure 12),the lithology presented on the face of the surrounding rock mass of the largely-deformed tunnel section (Figure 13) is dominated by the carbonaceous phyllite and sericite phyllite, mainly the Grade-V surrounding rock mass. The rock mass is found with a laminated texture and has a buried depth of 646 m, with a vertical stress of 21.7 MPa and a horizontal stress of 8.3 MPa. Field observations indicate minimal effects of underground water, and thus the correction factor λ for the effects of underground water on the weathering conditions of the rock mass is set as one.
The pre-excavation boreholes are arranged on the tunnel face (Figure 14), the rock cores are collected via pre-excavation boreholes for investigation, and the RBI value of rock mass (Table 5) for each borehole is calculated using Equation (9) and the core length measured on site. For the monitoring of surrounding rock mass deformation after excavation, the on-site monitoring point placement is illustrated in Figure 14.
Taking Point 4 as an example, the field investigation shows that the RBI of the surrounding rock mass cores at Point 4 is about 2.8; the AWI is about 0.65. According to Table 4, the GSI value is 33. For the intact phyllite rock, the Hoek–Brown constant m i = 10 and the uniaxial compressive strength σ c = 35 MPa. From Equation (5), we have σ 3 m a x = σ c 4 = 8.75   MPa and from Equations (2)–(4), we have m b = 0.916, s = 0.000587, and a = 0.518.
Finally, according to Equation (1), the uniaxial comprehensive strength of the rock mass σ r = 16.33 MPa and the elastic modulus E 0 = 1.11 GPa. The basic parameters at the other points are summarized in Table 6.
The fitting of the creep constitutive model (Equation (11)) to the surrounding rock mass monitoring data (Figure 15) was performed using the least square method to calculate the creep parameters of the rock mass at each monitoring point. The calculation results are shown in Table 7.
Since the factors affecting the deformation of the surrounding rock mass are varied in the different positions of the tunnel, a neural network was constructed and trained for each point to improve the inversion accuracy. The single-variable method was used for training, which means when the neural network of one point is being trained, the parameters of the other points are constant. Here, Point 4 was taken as an example to demonstrate the training of the neural network. As stated above, the parameters of the other points were fixed, and the learning sample dataset of the surrounding rock mass mechanical properties was built based on the orthogonal test design (Figure 16). The sample parameters in the training dataset were substituted into the numerical simulation model developed using Midas (Figure 17), and the forward calculation was performed to obtain the surrounding rock mass deformation magnitude at each point (Table 8). We used 80% of the data in Table 8 as the training dataset, whereas we used the other 20% as the testing dataset.
Section K143 + 621 was chosen for numerical simulation. Its buried depth was 646 m, we set the size of the model as long × wide = 360 m × 150 m, the fixed end restraint was set at the bottom, the overlying pressure (γH = 33.6 KN/m3 × (646 − 150) m = 16.67 MPa) was added above the model, the horizontal stress (8.3 MPa) was added on both sides of the model, and the surrounding rock masses all adopted the M-C failure criterion, as shown in Figure 17. The physical parameters of the surrounding rock mass, dominated by sericite and carbonaceous phyllite, were: average density ρ = 3.42 g/cm3, average elastic modulus E = 5.52 GPa, Poisson ratio ν = 0.38, internal cohesion c = 175 kPa, internal friction angle φ = 27°.
Using the BP neural network toolkit in MATLAB, the forward training of the rock mass parameters and deformation at Point 4 (Table 7) was performed, which produced the neural network model of Point 4, referred to as net4. Similarly, the training samples for the other point were generated and so were their neural network models, namely net1, net2, net3, and net5, respectively. Finally, the deformation measured at each point was input into the corresponding neural network model and the surrounding rock mass parameters were inverted. For the purpose of validation, the inverted mechanical parameters of the surrounding rock mass were substituted into the Midas numerical model to determine the deformation at each point of the surrounding rock mass (Table 9).
The relative error between the measured displacement and the numerically simulated displacement obtained through inversion was within 15%, indicating the good consistency and high applicability of the proposed method. Therefore, one can build a monitoring system for the largely-deformed surrounding rock mass section, according to the on-site information (such as the lithology combination, underground water status, bedding thickness, rock mass texture, surrounding rock mass integrity, maximum principal stress, and on-site measured rock deformation) for the purposes of gaining real-time information of the tunnel’s surrounding rock mass parameter, adjusting the support plan in a timely manner, and providing references for design and construction.

5. Conclusions

Due to the concealment of tunnel engineering, it is difficult to obtain the parameters of the surrounding rock mass after excavation. In this study, we combined the traditional Nashihara model with the Hoek–Brown failure criterion and developed the creep constitutive model for the jointed rock mass. Then, the creep parameters of the surrounding rock mass of the Ganbao tunnel were inverted and validated by integrating the on-site monitoring and the BP neural network. The following conclusions were drawn:
(1)
By reviewing the typical cases of large deformation in soft rock tunnels, the main influential factors can be summarized as the lithology combination, weathering effect, and underground water status. With the classical rock mass failure criterion, it is hard to thoroughly incorporate the geological characteristics of the actual rock mass and therefore the semi-empirical semi-theoretical Hoek–Brown approach is more fit-for-purpose.
(2)
The geological characteristics of the engineering rock mass were quantitatively characterized using two indexes, namely, the rock mass block index (RBI) and the absolute weathering index (AWI). Following the Hoek–Brown criterion, the long-term strength σ r and elastic modulus E 0 of the rock were obtained and then substituted into the rock mass creep constitutive model based on the Nashihara model and the Hoek–Brown failure criterion. By doing so, the original five creep parameters that needed to be determined in the creep equation were reduced to three, which simplifies the calculation of the constitutive model, while well reflecting the engineering, geological, and creep characteristics of the rock mass on site.
(3)
Considering the fact that the actual engineering geology varies at different positions of the tunnel’s surrounding rock mass, a specific BP neural network model was built for each monitoring point. Then, the rock mass parameters at each point were inverted from the on-site measured deformation. At last, the inverted parameters were input into the numerical model to calculate the deformation at each point, which was compared with the corresponding measured deformation. The resultant errors were all within 15%, satisfying the engineering requirements and demonstrating the reliability of the proposed method.
(4)
The values of the rock mass GSI (Figure 2) were all determined according to the engineering geology handbook, relevant standards, and existing literature. Therefore, these values can be adjusted as per the field condition. The inversion of the surrounding rock mass parameters is highly affected by the basic parameters of the engineering geological characteristics and the in-situ stress field of the tunnel’s surrounding rock mass. Hence, during applications of the proposed method, the basic parameters need to be accurately measured to improve the accuracy of the inversion.
(5)
This method is suitable for tunneling in unsupported rock mass or plainly supported tunnels after excavation.

Author Contributions

Writing—original draft preparation, C.C.; writing—review and editing, T.L.; methodology, C.M.; investigation, H.Z.; software, J.T.; data curation, Y.Z.; funding acquisition, T.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number U19A20111.

Institutional Review Board Statement

Not involved.

Informed Consent Statement

Not involved.

Data Availability Statement

The study did not report any data.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (NO. U19A20111). The financial support is gratefully acknowledged. Special thanks are due to Tianbin Li and Chunchi Ma for their help in improving this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, M.; Jing, H.; Sun, X. Mechanics of Soft Rock Engineering; Science Press: Beijing, China, 2002. [Google Scholar]
  2. Kong, P.; Jiang, L.S.; Shu, J.M.; Sainoki, A.; Wang, Q.B. Effect of fracture heterogeneity on rock mass stability in a highly heterogeneous underground roadway. Rock Mech. Rock Eng. 2019, 52, 4547–4564. [Google Scholar] [CrossRef]
  3. Chen, Z.; He, C.; Xu, G.; Ma, G.; Wu, D. A case study on the asymmetric deformation characteristics and mechanical behavior of deep-buried tunnel in phyllite. Rock Mech. Rock Eng. 2019, 52, 4527–4545. [Google Scholar] [CrossRef]
  4. Jia, P.; Tang, C.; Yang, T.; Wang, S. Numerical Stability Analysis of Surrounding rock mass Mass Layered by Structural Planes with Different Obliquities. J. Northeast Univ. 2006, 27, 1275–1278. [Google Scholar]
  5. Li, L.; Guan, J.; Xiao, M.; Liu, H.; Tang, K. A creep constitutive model for transversely isotropic rocks. Rock Soil Mech. 2020, 41, 2922–2930, 2942. [Google Scholar]
  6. Li, C.; Wang, J.; Xie, H. Anisotropic creep characteristics and mechanism of shale under elevated deviatoric stress. J. Pet. Sci. Eng. 2020, 185, 106670. [Google Scholar] [CrossRef]
  7. Li, A.; Xu, N.; Dai, F.; Gu, G.; Hu, Z.; Liu, Y. Stability analysis and failure mechanism of the steeply inclined bedded rock masses surrounding a large underground opening. Tunn. Undergr. Space Technol. 2018, 77, 45–58. [Google Scholar] [CrossRef]
  8. Song, D.; Chen, J.; Cai, J. Deformation monitoring of rock slope with weak bedding structural plane subject to tunnel excavation. Arab. J. Geosci. 2018, 11, 251. [Google Scholar] [CrossRef]
  9. Gao, F.; Guo, J. Stress mechanism and deformation monitoring of bias tunnel. Int. J. Simul. Syst. Sci. Technol. 2016, 17, 14.1–14.4. [Google Scholar] [CrossRef]
  10. Yin, C.; Li, H.; Che, F.; Li, Y.; Hu, Z.; Liu, D. Susceptibility mapping and zoning of highway landslide disasters in China. PLoS ONE 2020, 15, e0235780. [Google Scholar] [CrossRef]
  11. Liu, C.H.; Li, Y.Z. Analytical study of the mechanical behavior of fully grouted bolts in bedding rock slopes. Rock Mech. Rock Eng. 2017, 50, 2413–2423. [Google Scholar] [CrossRef]
  12. Vergara, M.R.; Kudella, P.; Triantafyllidis, T. Large scale tests on jointed and bedded rocks under multi-stage triaxial compression and direct shear. Rock Mech. Rock Eng. 2015, 48, 75–92. [Google Scholar] [CrossRef]
  13. Li, L.; Tan, Z.S.; Guo, X.L.; Wu, Y.; Luo, N. Large deformation of tunnels in steep dip strata of interbedding phyllite under high geostresses. Chin. J. Rock Mech. Eng. 2017, 36, 1611–1622. [Google Scholar]
  14. Li, X.H.; Xia, B.W.; Li, D.; Han, C.R. Deformation characteristics analysis of layered rockmass in deep buried tunnel. Rock Soil Mech. 2010, 31, 1163–1167. [Google Scholar]
  15. Fu, X. Numerical Simulation of Asymmetric Large-Deformation Energy-Releasing Bolt Support for Layered Soft Rock Tunnel; Chengdu Univerisity of Technology: Chengdu, China, 2020. [Google Scholar]
  16. Li, Z.; Shan, R.; Wang, C.; Yuan, H.; Wei, Y. Study on the distribution law of stress deviator below the floor of a goaf. Geomech. Eng. 2020, 21, 301–313. [Google Scholar] [CrossRef]
  17. Tian, M.; Han, L.; Meng, Q.; Ma, C.; Zong, Y.; Mao, P. Physical model experiment of surrounding rock mass failure mechanism for the roadway under deviatoric pressure form mining disturbance. KSCE J. Civ. Eng. 2020, 24, 1103–1115. [Google Scholar] [CrossRef]
  18. Wang, Z.-J.; Luo, Y.-S.; Guo, H.; Tian, H. Effects of initial deviatoric stress ratios on dynamic shear modulus and damping ratio of undisturbed loess in China. Eng. Geol. 2012, 143–144, 43–50. [Google Scholar] [CrossRef]
  19. Kroon, M.; Faleskog, J. Numerical implementation of a J2 and J3 dependent plasticity model based on a spectral decomposition of the stress deviator. Comput. Mech. 2013, 52, 1059–1070. [Google Scholar] [CrossRef]
  20. Wang, Q.; Pan, R.; Jiang, B.; Li, S.; He, M.; Sun, H.; Wang, L.; Qin, Q.; Yu, H.; Luan, Y. Study on failure mechanism of roadway with soft rock in deep coal mine and confined concrete support system. Eng. Fail. Anal. 2017, 81, 155–177. [Google Scholar] [CrossRef]
  21. Wang, W.; Zhang, C.; Wei, S.; Zhang, X.; Guo, S. Whole section anchor–grouting reinforcement technology and its application in underground roadways with loose and fractured surrounding rock mass. Tunn. Undergr. Space Technol. 2016, 51, 133–143. [Google Scholar] [CrossRef]
  22. Xu, G.; He, C.; Chen, Z.; Yang, Q. Transversely isotropic creep behavior of phyllite and its influence on the long-term safety of the secondary lining of tunnels. Eng. Geol. 2020, 278, 105834. [Google Scholar] [CrossRef]
  23. Li, X.; Ju, M.; Yao, Q.; Zhou, J.; Chong, Z. Numerical investigation of the effect of the location of critical rock block fracture on crack evolution in a gob-side filling wall. Rock Mech. Rock Eng. 2016, 49, 1041–1058. [Google Scholar] [CrossRef]
  24. Zhou, J.; Wei, Q.; Liu, G. Back analysis on rock mechanics parameters for highway tunnel by BP neural network method. Chin. J. Rock Mech. Eng. 2004, 23, 941–945. [Google Scholar]
  25. Cao, W.; Jiang, Y.; Sakaguchi, O.; Li, N.; Han, W. Predication of Displacement of Tunnel Rock Mass Based on the Back-Analysis Method-BP Neural Network. Geotech. Geol. Eng. 2021, 2021, 1–14. [Google Scholar] [CrossRef]
  26. Wu, Q.; Yan, B.; Zhang, C.; Wang, L.; Ning, G.; Yu, B. Displacement prediction of tunnel surrounding rock mass: A comparison of support vector machine and artificial neural network. Math. Probl. Eng. 2014, 2014, 351496. [Google Scholar] [CrossRef]
  27. Deng, X.; Xu, T.; Wang, R. Risk evaluation model of highway tunnel portal construction based on BP fuzzy neural network. Comput. Intell. Neurosci. 2018, 2018, 8547313. [Google Scholar] [CrossRef] [Green Version]
  28. Wen, H.; Yin, J.; Qin, Z.; Xie, R. Application of BP Neural Network to the Back Analysis of Mechanical Parameters of Tunnel Surrounding rock mass. J. Yangtze River Sci. Res. Inst. 2013, 30, 47–51, 56. [Google Scholar]
  29. He, Y.; Sun, X.; Zhang, Y.; Guo, H.; Li, Q. Intelligent fusion model and analysis method for rock parameter inversion of water diversion tunnel. J. Hydroelectr. Eng. 2021, 40, 114–126. [Google Scholar]
  30. Ma, F.; Jia, S. Back analysis of elastoplastic parameters of surrounding rock mass for roadway in mudstone and its long-term stability prediction. Rock Soil Mech. 2014, 7, 1987–1994. [Google Scholar]
  31. Amadei, B. Strength of a regularly jointed rock mass under biaxial and axisymmetric loading conditions. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1988, 25, 3–13. [Google Scholar] [CrossRef]
  32. Xiao, S.; Yang, S. Rock Mass Mechanics; Geological Publishing House: Beijing, China, 1987. [Google Scholar]
  33. Hoek, E.; Brown, E.T. Practical estimates of rock mass strength. Int. J. Rock Mech. Min. Sci. 1997, 34, 1165–1186. [Google Scholar] [CrossRef]
  34. Marinos, P.; Hoek, E. Estimating the geotechnical properties of heterogeneous rock masses such as flysch. Bull. Eng. Geol. Environ. 2001, 60, 85–92. [Google Scholar] [CrossRef]
  35. Hoek, E.; Brown, E.T. The Hoek–Brown failure criterion and GSI–2018 edition. J. Rock Mech. Geotech. Eng. 2019, 11, 445–463. [Google Scholar] [CrossRef]
  36. Hu, X.; Zhong, P.; Ren, Z. Rock-mass block index and its engineering practice significance. J. Hydraul. Eng. 2002, 33, 80–83. [Google Scholar]
  37. Parker, A. An index of weathering for silicate rocks. Geol. Mag. 1970, 107, 501–504. [Google Scholar] [CrossRef]
  38. Su, Y.; Feng, L.; Li, Z.; Zhao, M. Quantification of elements for geological strength index in Hoek-Brown criterion. Chin. J. Rock Mech. Eng. 2009, 28, 36–43. [Google Scholar]
  39. Li, R.; Wu, L. Research on characteristic indexes of weathering intensity of rocks. Chin. J. Rock Mech. Eng. 2004, 23, 3830. [Google Scholar]
Figure 1. Schematic diagram of rock weathering.
Figure 1. Schematic diagram of rock weathering.
Applsci 11 10033 g001
Figure 2. GSI quantification of rock mass [31]. Note: The following example is used to demonstrate how to use the above table. In the case of AWI = 0.52 and RBI = 9.7, two vertical and horizontal lines are drawn according to the index values (the dashed lines in the table); the intersection point is found between 40 and 45; the GSI value at the point is calculated via linear interpolation (the result is 41).
Figure 2. GSI quantification of rock mass [31]. Note: The following example is used to demonstrate how to use the above table. In the case of AWI = 0.52 and RBI = 9.7, two vertical and horizontal lines are drawn according to the index values (the dashed lines in the table); the intersection point is found between 40 and 45; the GSI value at the point is calculated via linear interpolation (the result is 41).
Applsci 11 10033 g002
Figure 3. Combination of the creep elements for rock mass.
Figure 3. Combination of the creep elements for rock mass.
Applsci 11 10033 g003
Figure 4. Surrounding rock mass strength parameters based on the Hoek–Brown failure criterion.
Figure 4. Surrounding rock mass strength parameters based on the Hoek–Brown failure criterion.
Applsci 11 10033 g004
Figure 5. Schematic diagram of the surrounding rock mass deformation.
Figure 5. Schematic diagram of the surrounding rock mass deformation.
Applsci 11 10033 g005
Figure 6. Schematic diagram of the BP neural network structure. Note: The weight w represents the connection strength; the bias b indicates whether or not the node is easy to activate (threshold).
Figure 6. Schematic diagram of the BP neural network structure. Note: The weight w represents the connection strength; the bias b indicates whether or not the node is easy to activate (threshold).
Applsci 11 10033 g006
Figure 7. Deep BP neural network structure.
Figure 7. Deep BP neural network structure.
Applsci 11 10033 g007
Figure 8. Steps for training the BP neural network.
Figure 8. Steps for training the BP neural network.
Applsci 11 10033 g008
Figure 9. Inversion workflow for the creep parameters of the surrounding rock mass.
Figure 9. Inversion workflow for the creep parameters of the surrounding rock mass.
Applsci 11 10033 g009
Figure 10. Longitudinal profile of the tunnel.
Figure 10. Longitudinal profile of the tunnel.
Applsci 11 10033 g010
Figure 11. The twisted steel arch and back haunch fillet of the arch, and rockfall of the initial support.
Figure 11. The twisted steel arch and back haunch fillet of the arch, and rockfall of the initial support.
Applsci 11 10033 g011
Figure 12. The overburden near monitoring section K143 + 621.
Figure 12. The overburden near monitoring section K143 + 621.
Applsci 11 10033 g012
Figure 13. Tunnel face (unexcavated) at K143 + 630. (a) Geological sketch picture. (b) Photo of the tunnel face.
Figure 13. Tunnel face (unexcavated) at K143 + 630. (a) Geological sketch picture. (b) Photo of the tunnel face.
Applsci 11 10033 g013
Figure 14. Monitoring point and pre-excavation borehole placement in the tunnel (K143 + 621).
Figure 14. Monitoring point and pre-excavation borehole placement in the tunnel (K143 + 621).
Applsci 11 10033 g014
Figure 15. Measured displacement vs. time at each monitoring point.
Figure 15. Measured displacement vs. time at each monitoring point.
Applsci 11 10033 g015
Figure 16. Training the dataset at a certain point in a single-variable approach.
Figure 16. Training the dataset at a certain point in a single-variable approach.
Applsci 11 10033 g016
Figure 17. Numerical simulation model of the tunnel.
Figure 17. Numerical simulation model of the tunnel.
Applsci 11 10033 g017
Table 1. Deformation-failure characteristics of tunnels’ surrounding rock masses and support structures under complex geological conditions in western China.
Table 1. Deformation-failure characteristics of tunnels’ surrounding rock masses and support structures under complex geological conditions in western China.
Tunnel NameTunnel Length
/m
Buried Depth
/m
Lithology Max   Principal   Stress   σ 1 m a x
/MPa
Strength of Rock
/MPa
Maximum
Deformation
/mm
Deformation CharacteristicsStrain ε = Tunnel Closure/Tunnel Diameter × 100
Maoxian Tunnel,
Chenglan Railway
25,000675Sericite phyllite, carbonaceous phyllite27.521.95510Compressive large deformation occurs, with prolonged deformation growth and notable time-dependency; horizontal convergence exceeds crown settlement, due to steep inclination of the surrounding rock mass8.05
Wuqiaoling Tunne,
Lanxin Railway
20,0501100Phyllite, slate32.80.7~2.51209It penetrates the compressive fault; the overall stability of the surrounding rock mass is low; large deformation, early rapid deformation growth, and prolonged duration of deformation are observed due to intensive compression.9.69
Zhegushan Tunnel,
317 National Rd.
44231000Thin layers of carbonaceous phyllite17~2012300The relatively large magnitude and prolonged duration of surrounding rock mass deformation are manifested as support breakage, steel arch twisting, and their intrusion into tunnel clearance; the tunnel is prone to collapse.9.31
Maoyushan Tunnel,
Lanyu Railway
8503700Thin bedded slate21.285.63540Large rapid deformation, with notable rheological effects; severe twisting and fracturing of the steel arch; horizontal convergence far larger than crown settlement8.94
Gonghe Tunnel,
Yusha Expressway
47791000Sandy shale29.8611.4200Longitudinal cracking and steel frame bending occurs at the initial support of the right spandrel and left arch foot, indicating severe biased compression of the surrounding rock mass7.52
Table 2. Representation of rock mass textures by RBI.
Table 2. Representation of rock mass textures by RBI.
Texture TypeRBIRock Mass Characteristics
Laminated mosaic texture30–10Relatively intact, barely, or partially disturbed, often developing 3 groups of structural planes with the spacing of 30–50 cm
Mosaic texture10–3Less intact, mostly disturbed, broken yet with tightly packed fragments, generally developing 3–4 groups of structural planes with the spacing of 10–30 cm
Broken texture3–1Broken rock mass, sufficiently disturbed, composed of fragments or thin layers, with extensive structural planes presenting spacing generally smaller than 10 cm
Loose texture1–0Extremely crushed rock mass, extremely disturbed, composed of loose rock blocks, and angular fragments with crushed debris
Table 3. Representation of the rock mass weathering characteristics by AWI.
Table 3. Representation of the rock mass weathering characteristics by AWI.
Weathering Condition A W I 0 Weathering Characteristics
Non-weathered>0.90Very good fracture surface: very coarse, fresh, indicating well-sealed fresh rock matrix and no weathering
Slightly weathered0.90–0.75Good fracture surface: coarse, relatively fresh, with the presence of rush and the slight alteration of minerals with low weathering resistance, indicating slight weathering
Weakly weathered0.75–0.55Ordinary fracture surface: smooth with no filling, partial alteration of minerals with low weathering resistance, indicating weak weathering
Intensively weathered0.55–0.35Poor fracture surface: the presence of slickensides, covering of tight films or filling of angular debris on the surface, high alteration of minerals with low weathering resistance, indicating intensive weathering
Extremely weathered≤0.35Very poor fracture surface: the presence of slickensides, and soft clay films or clay filling; the vast majority of minerals with low weathering resistance are altered; indicating extreme weathering
Table 4. Recommended values for the correction factor λ to capture the effects of underground water on the weathering conditions of rock mass.
Table 4. Recommended values for the correction factor λ to capture the effects of underground water on the weathering conditions of rock mass.
Production Status of Underground WaterRBI Values of Jointed Rock Mass
30–1010–33–11–0
Humid or dripping0.950.95–0.890.89–0.830.83–0.76
Rain-like or spring-like production with water pressure < 0.1 MPa; or unit water production rate < 10 L/min·m0.95–0.890.89–0.830.83–0.760.76–0.71
Rain-like or spring-like production with water pressure > 0.1 MPa; or unit water production rate > 10 L/min·m0.89–0.830.83–0.760.76–0.710.71–0.67
Table 5. Cores collected via pre-excavation boring at each monitoring point and the RBI values.
Table 5. Cores collected via pre-excavation boring at each monitoring point and the RBI values.
Monitoring PointCoresRBI Values
1 Applsci 11 10033 i0016.763
2 Applsci 11 10033 i0025.626
3 Applsci 11 10033 i0033.748
4 Applsci 11 10033 i0042.8
5 Applsci 11 10033 i0054.331
Table 6. Basic parameters of the surrounding rock mass at each monitoring position.
Table 6. Basic parameters of the surrounding rock mass at each monitoring position.
Monitoring PointRBIAWIGSI m b s a σ r
(MPa)
σ 3 = σ c 4
(MPa)
E 0
(GPa)
16.7631.57079.7082.2120.0014180.2736.845.292.69
25.6261.30666.3051.8400.0011790.2232.2837.672.24
33.7480.87044.1691.2260.0007860.1521.5725.101.49
42.80.65330.9160.0005870.1116.3318.751.11
54.3311.00651.0501.4170.0009080.1722.829.011.73
Table 7. Creep parameters of the rock mass at each monitoring point.
Table 7. Creep parameters of the rock mass at each monitoring point.
Point No. E 0 / GPa σ r   ( MPa ) E 1 / GPa η 1 / ( GPa · d ) η 2 / ( GPa · d )
12.6936.80.816.3893.71
22.2432.280.7111.2288.63
31.4921.570.58.786.55
41.1116.330.47.875.62
51.7322.80.610.1178.52
Table 8. The orthogonal test design plan and calculation results for Point 4.
Table 8. The orthogonal test design plan and calculation results for Point 4.
No. E 0
/GPa
σ r
/MPa
E 1
/GPa
η 1
/(GPa·d)
η 2
/(GPa·d)
Displacement at Point 1
/mm
Displacement at Point 2
/mm
Displacement at Point 3
/mm
Displacement at Point 4
/mm
Displacement at Point 5
/mm
12.6936.800.8016.3893.717.319.1113.1127.9411.50
22.6941.801.3021.3898.7112.5115.5922.4347.8019.67
32.6946.801.8026.38103.7111.0213.7419.7642.1117.33
42.6951.802.3031.38108.719.8412.2717.6437.6015.47
52.6956.802.8037.38113.718.9011.0915.9534.0013.99
62.2436.801.3026.38108.7115.8319.7325.8533.6622.41
72.2441.801.8031.38113.7114.3117.8423.3730.4320.26
82.2446.802.3037.3893.7111.7614.6619.2125.0116.65
92.2451.802.8016.3898.7117.7322.1028.9537.7025.10
102.2456.800.8021.38103.7110.9113.6017.8223.2015.45
111.4936.801.8037.3898.7110.6013.2118.7430.8916.44
121.4941.802.3016.38103.718.7110.8515.4025.3813.51
131.4946.802.8021.38108.7114.9018.5726.3543.4323.11
141.4951.800.8026.38113.7113.1316.3623.2238.2620.36
151.4956.801.3031.3893.7111.7214.6120.7334.1618.18
161.1136.802.3021.38113.7123.7328.0540.7651.4934.56
171.1141.802.8026.3893.7120.9124.7135.9145.3730.45
181.1146.800.8031.3898.7118.6722.0732.0640.5127.19
191.1151.801.3037.38103.7116.8819.9528.9936.6224.58
201.1156.801.8016.38108.7113.8716.3923.8230.0920.20
211.7336.802.8031.38103.7118.2621.8525.0831.6921.27
221.7341.800.8037.38108.7116.3119.5118.1222.8915.37
231.7346.801.3016.38113.7114.7417.6416.3920.7013.89
241.7351.801.8021.3893.7112.1214.4913.4617.0111.42
251.7356.802.3026.3898.7120.7324.8033.2341.4628.48
Table 9. Parameter inversion of surrounding rock mass and error analysis.
Table 9. Parameter inversion of surrounding rock mass and error analysis.
Monitoring
Point
Measured Deformation
/mm
Output: Inverted Mechanical Parameters of the Surrounding Rock MassSimulated Deformation
/mm
Relative Error
/%
E 0
/GPa
σ r
/MPa
E 1
/GPa
η 1
/(GPa·d)
η 2
/(GPa·d)
123.733.4435.350.9216.85104.4926.6712.37
228.052.9131.890.8112.3997.2030.699.41
340.762.1221.670.5910.019.1345.5711.81
451.501.7116.910.489.198.1158.5413.67
534.562.3922.630.7111.1588.0037.518.53
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, C.; Li, T.; Ma, C.; Zhang, H.; Tang, J.; Zhang, Y. Hoek-Brown Failure Criterion-Based Creep Constitutive Model and BP Neural Network Parameter Inversion for Soft Surrounding Rock Mass of Tunnels. Appl. Sci. 2021, 11, 10033. https://doi.org/10.3390/app112110033

AMA Style

Chen C, Li T, Ma C, Zhang H, Tang J, Zhang Y. Hoek-Brown Failure Criterion-Based Creep Constitutive Model and BP Neural Network Parameter Inversion for Soft Surrounding Rock Mass of Tunnels. Applied Sciences. 2021; 11(21):10033. https://doi.org/10.3390/app112110033

Chicago/Turabian Style

Chen, Chao, Tianbin Li, Chunchi Ma, Hang Zhang, Jieling Tang, and Yin Zhang. 2021. "Hoek-Brown Failure Criterion-Based Creep Constitutive Model and BP Neural Network Parameter Inversion for Soft Surrounding Rock Mass of Tunnels" Applied Sciences 11, no. 21: 10033. https://doi.org/10.3390/app112110033

APA Style

Chen, C., Li, T., Ma, C., Zhang, H., Tang, J., & Zhang, Y. (2021). Hoek-Brown Failure Criterion-Based Creep Constitutive Model and BP Neural Network Parameter Inversion for Soft Surrounding Rock Mass of Tunnels. Applied Sciences, 11(21), 10033. https://doi.org/10.3390/app112110033

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop