Research Article  Open Access
The Influence of Moisture Content on the TimeDependent Characteristics of Rock Material and Its Application to the Construction of a Tunnel Portal
Abstract
Uniaxial compression creep experiments were carried out for lowgrade metamorphic slate samples (located in the southeastern area of Guizhou province, China) with different moisture contents, using an Instron electricfluid servocompression machine. Based on the experimental results, a detailed analysis was made of the effect of moisture content on the strength and deformation behaviour of the slate specimens. The threeparameter generalised Kelvin model was identified to describe the creep behaviour of the lowgrade metamorphic slate with different moisture contents. There is an approximately linear negative correlation between the elastic modulus and the saturation degree, and the viscoelastic modulus and viscosity coefficient show a negative exponent correlation with the saturation degree. The Kelvin creep model considering the moisture degradation effect was established and a threedimensional finite difference model was developed with the software to validate the creep model. A threedimensional numerical analysis was then performed to simulate the tunnel excavation process. The results show that the influence of moisture and creep of the surrounding rocks is needed for estimating the deformation of complex tunnel portals.
1. Introduction
The influence of moisture on the strength and deformability of rock has been analysed by many researchers [1–3]. Dyke and Dobereiner [4] demonstrated the variation in uniaxial compressive with moisture content for three quartz arenites ranging in dry strength from approximately 34 to 70 MPa. They concluded that in general the weaker the rock is, the more sensitive it is to changes in moisture content. Atkinson and Meredith [5] investigated experimentally the effect of water on subcritical microcrack growth parameters of granitic rocks. They found that the increase of humidity does not have a significant influence on the subcritical microcrack growth index but results in a large increase in microcrack growth velocity. Golshani et al. [6] checked the effect of moisture conditions on the timedependent behaviour of Inada granite. The experimental and theoretical results have shown that the timedependent behaviours of rocks are seriously affected by water [7, 8].
Tunnel entrance construction presents challenges and because of its complexity it could easily result in landslides or collapse [9–13]. Generally, the cover depth at the tunnel entrance is shallow as it often moves through a mountain surface layer whose rocks are usually broken and seriously weathered. So, the surrounding rocks at a tunnel portal are easily affected by surface water [14], which makes the rheological properties of surrounding rocks more obvious and involves many different complicated problems due to the uncertainty associated with the response of the rocks to the process of excavation and construction [15].
In this study, uniaxial compression creep experiments were carried out for the lowgrade metamorphic slate samples (the surrounding rock of the Ruipo tunnel portal, which is located in the southeastern area of Guizhou province, China) with different moisture contents using an Instron electricfluid servocompression machine. A detailed analysis was made of the effect of moisture content on the strength and deformation behaviour of slate specimens. Then a threedimensional finite difference analysis was executed on expressway tunnel portal areas while considering tunnel excavation and the moisture degradation effects on creep characteristics of surrounding rocks. In addition, the results were compared with those obtained by common elastoplastic models.
2. Compression and Creep Tests
2.1. Specimen Preparation and Experimental Procedure
Lowgrade metamorphic slate, a quintessential rock of the ShuiDu expressway that was taken from the Ruipo tunnel portal, has a cataclastic texture. Xray diffraction analysis reveals that, apart from the most common quartz, kaolinite and muscovite are present in the collected lowgrade metamorphic slate (see Figure 1). Kaolinite has a rough surface and perfect cleavage, strong water absorption, no water swelling, great plasticity in wet conditions, and a great influence on rock strength. There were 21 cylindrical test specimens, namely, S1 to S21, all 50 mm in diameter and 100 mm in length, which was in accordance with the tested sample size suggested by the ISRM [16]. Three degrees of saturation of specimens were considered for experiments and 6 specimens were prepared at each degree of saturation.
(a)
(b)
3 specimens, S1, S2, and S3, were used to obtain the degree of saturation at different times. Nature dry specimens were considered as the first degree of saturation, specimens submerged in water for 24 hours represented the second degree, and those that remained submerged in water for 240 hours represented the third degree of saturation. The specimen preparation process was as follows. First, all specimens were naturedried for 48 hours, the weights of S1, S2, and S3 were measured, and S4 to S9 were selected as the specimens at the first degree of saturation. Second, all specimens except S4 to S9 were submerged in water, and the weights of S1, S2, and S3 were measured every 24 hours. After having been submerged in water for 24 hours, S10 to S15 were selected as the specimens at the second degree of saturation, and S16 to S21 were specimens at the third degree of saturation after having been submerged in water for 240 hours. In order to avoid water loss from the specimens during compression testing, S10 to S21 were sealed up with wax (see Figure 2). Then, the dried weights of S1, S2, and S3 were obtained after having been ovendried for 24 hours at 105°C. Finally, S1, S2, and S3 were saturated using the vacuum saturation method [17]. The samples were placed in the vacuum saturation apparatus and a sufficient amount of water at approximately 25°C was added to cover the sample. Entrapped air was removed from the samples by applying a residual manometer pressure of 100 kPa. The pressure was maintained for 4 hours and then the saturated weights of samples were measured.
Uniaxial creep experiments were performed on an Instron 1346 electronic hydraulic servocontrolled testing machine, the loading capacity of which is 2 MN. Before the creep tests, uniaxial compression tests were performed on two specimens at each degree of saturation, and triaxial compression experiments were also carried out at each degree of saturation for specimens under three different confining pressures , that is, 2, 5, and 10 MPa. The instantaneous failure strengths and corresponding peak strains were obtained at failure of these specimens, and creep tests were planned from a knowledge of these values. The multilevel test was used: the load was applied at low levels, rising to a high level on the same specimen. The duration of tests was dependent on the displacement rates. In the case when the rate of displacement was less than 0.001 mm/h, the creep under this load level was considered to be basically stable [18]. The next load level could be executed afterwards, and the loading rating was 500 N/s.
2.2. Experimental Results and Discussion
2.2.1. Relationship between Saturation Degree and Time Submerged in Water
The degree of saturation of specimens at different times can be defined aswhere is the weight of the specimen after being soaked for a period of time; and are the dry weight and saturated weight of the specimen, respectively.
The experimental results showed that the degree of saturation increased with the time submerged in water (see Figure 3), and the average degrees of saturation of specimens were 41.5%, 80.1%, and 93.5% when submerged in water for 0, 24, and 240 hours, respectively.
2.2.2. Mechanical Parameters Obtained from Uniaxial and Triaxial Compression Tests
Table 1 lists the mechanical parameters of specimens at three saturation degrees. The uniaxial compression strength of specimens with saturation degrees of 80.1% and 93.5% is, respectively, 69.85 and 58.74 MPa, which decreases by 24.2% and 36.2% compared with nature dry specimens, as well as the fact that elastic modulus has a reduction of 64.4% and 83.8% correspondingly. As a conclusion, a great decrease could be found in elastic modulus when compared with strength, and the increase of moisture causes an increase in Poisson’s ratio.

For the triaxial test, the linear MohrCoulomb criterion, which can be expressed simply as a linear relationship between peak axial stress and confining pressure , that is, , is selected here to determine the peak strength parameters of cohesion and internal friction angle [19]. The values of and for the specimens in these tests decreased with the increase of saturation degree (see Figure 4). Several studies [20, 21] concluded that there was usually an exponential relation between the strength of rock material and moisture. The values of and of the slate obtained by Li et al. [20] also had an exponential dependency on the rock moisture, and so the exponential function was chosen to describe the relationship between strength parameters and degree of saturation:where is the saturation degree of specimens (%).
(a)
(b)
2.2.3. Creep Behaviour
Uniaxial experimental creep results for specimens are shown in Figure 5. Overall, transient and creep axial strains of the specimens increase with degree of saturation. At MPa, the transient axial strain of the three specimens was, respectively, 0.073% , 0.074% , and 0.413% , and the axial strain of specimens increased by 0.003%, 0.026%, and 0.066% correspondingly after creeping for 12 hours. At MPa, the transient axial strain of the three specimens was, respectively, 0.092% , 0.137% , and 0.582% , and the axial strain of specimens increased by 0.003%, 0.080%, and 0.041% correspondingly after creeping for 12 hours. When reaching the failure stress level, the specimens showed tertiary creep behaviour and succumbed to creep failure. The failure stress level of specimens was the highest while the saturation degree was 41.5% (66.2 MPa) and then decreased when the moisture content was higher (56.1 and 38.2 MPa for specimens with % and 93.5%).
(a)
(b)
(c)
From the compression and creep tests results, it can be seen that the moisture has a great effect on the strength and deformability of the lowgrade metamorphic slate. Hawkins and McConnell [3] studied the sensitivity of sandstone strength and deformability to changes in moisture content. They concluded that the sensitivity increases progressively with higher proportions of clay minerals and rock fragments. The lowgrade metamorphic slate studied in this paper has 31% sericite and 10% siliceous cement, which have similar characteristics with clay minerals. Another important reason lies in the fact that the welldeveloped fissures might make it much easier for water to get into the rock specimens.
3. Creep Model and Validation
3.1. Creep Model considering Moisture Degradation Effect
From Figure 5, the deformation of specimens reached stability after a certain period of time and the strain rate reduced to zero. Therefore, the threeparameter generalised Kelvin model (see Figure 6), which connects in series an elastic component (Hookean body) and a Kelvin model, can be used to describe the deformation characteristics [22]. The onedimensional creep equation of the Kelvin model iswhere , , and are the elastic modulus, the viscoelastic modulus, and the viscosity coefficient, respectively.
For the creep test, in general, the least square method is used to statistically process the test data [23]. The estimated creep parameters using the least square method are shown in Table 2. There is an approximately linear negative correlation between the elastic modulus and the saturation degree (see Figure 7(a)):

(a)
(b)
(c)
The viscoelastic modulus and viscosity coefficient show a negative exponent correlation with saturation degree (see Figures 7(b)7(c)):
It can be seen that all the Kelvin model parameters, , , and , decrease with the saturation degree. Therefore, the damage to the rock induced by water can be considered as the damage of each creep parameter. The damage variable of the elastic modulus is defined as follows:where is the elastic modulus of rock with saturation degree and is the elastic modulus of rock with saturation degree . Substituting (4) into (6), we have
In the same way, the damage variables of the viscoelastic modulus and viscosity coefficient can be described as follows:Substituting (7)(8) into (3), the Kelvin creep model considering the moisture degradation effect can be expressed as
3.2. Validation of the Creep Model
To validate the Kelvin creep model considering the moisture degradation effect, a threedimensional finite difference model is developed with [24]. In classical fluid mechanics there are the following three hypotheses [25]: (1) the creep deformation of rock material results from the partial stress tensor, but the spherical stress tensor does not cause the creep deformation; that is, no volume flow occurs during the creep deformation; (2) rock is a kind of isotropic material, and the shortterm stressstrain curve and creep curve in the tensile and compressive stress states are very similar; and (3) the Poisson’s ratio of rock material is not dependent on time during the creep deformation. Based on the above hypotheses, we can derive the threeparameter generalised Kelvin model creep model in the threedimensional stress state, which is described as follows:where is the strain tensor, is the spherical stress tensor, is the Kronecker delta, is the partial stress tensor, is the bulk modulus, is the elastic shear modulus, and is the viscoelastic shear modulus. In the threedimensional space, we can obtain the following equations:
Due to the limitations of the Kelvin model in dealing with the plastic characteristic of rock material, a MohrCoulomb model is then used. The viscoelastic and plastic strain rate components are assumed to act in series. The viscoelastic constitutive law corresponds to the Kelvin model, and the plastic constitutive law corresponds to the MohrCoulomb model. The creep model is developed by using a user subroutine, which is an option provided by . The shape and size of the finite difference model are the same as for the tested specimens, and the mesh contains 1,000 elements and 1,111 nodes. The parameter values in the model are obtained from the laboratory test.
The calculated results for the specimens are shown in Figure 8. Although there is a certain deviation, the calculated curves coincide basically with the experimental curves with the same general rules, which demonstrates the validity of the creep model when considering the moisture degradation effect.
(a)
(b)
(c)
4. An Application of ThreeDimensional Analysis of the Ruipo Tunnel Portal
4.1. Project Background
The ShuiDu (Shuikou to Duyun city in Guizhou province, China) expressway is one of the most difficult sections in the construction of the Xiamen to Chengdu highway, which is the number 16 eastwest expressway of the National Trunk Highway System. 58 tunnels were designed on the ShuiDu expressway, with a total length of 71,988 m (34.6% of the whole line) and more than 200 tunnel portals. Almost all of the rocks surrounding the tunnel portals are broken and have poor stability. About half of the portals are shallow buried and unsymmetrically loaded, such as the Laozhai tunnel portal [26] and the Ruipo tunnel portal (see Figure 9). Another challenge was the rainy season that set in during the excavation of tunnel portals. The heavy and constant rainfall would lead to a serious decrease in surrounding rock quality. In order to overcome the difficulties in excavation, many prereinforcements (such as retaining wall and antislide pile) were built before excavation. Also, many complicated excavation methods, for example, the cross diaphragm (CRD) method [27] and ringlike drift heading method [28], were adopted. The excavation was slow (less than 2 m/d), but even then collapse occurred in many tunnel portals.
The Ruipo tunnel is located in AT11 section of the ShuiDu expressway, with starting and ending milestones of ZK49+865 and ZK51+395. The surrounding rocks of the tunnel portal section (ZK49+865 and ZK50+000) were composed mainly of strongweathered slate and completely weathered slate, which have very poor stability. The tunnel portal section was seriously shallow buried (the depth was 0–38 m) and unsymmetrically loaded (the transversal gradient was about 40°). According to design, the excavation in the tunnel portal section was to be performed by the CRD method. Before excavation the surrounding rocks were reinforced by a grouting method and a 108 mm diameter steel pipeline shed was installed along the tunnel outline, for which the intervals of steel were 50 cm. In order to mitigate the influence of an unsymmetrical load, seven antislide piles were installed in the shallow buried side. The support system and excavation sequence are shown in Figure 10.
(a)
(b)
During excavation rain began to fall and lasted for nearly 3 days, and then the portal suffered severe damage. Cracks first appeared in the slope, initial support, and landfills behind the antislide piles, and more cracks appeared and spread in the next few days (see Figure 11). For the reason that damage happened shortly after the tunnel portal excavation, only the forces of antislide piles were measured. The monitoring results are shown in Figure 12 (many measurements were not picked up for the reason that the earth pressure cells were destroyed). As can be seen in Figure 12, the maximum value of the forces was only 55 kPa; the low forces showed that antislide piles had less reinforcement effect. The reason for this is that the rainfalls reduced the loading capacity of the surrounding rocks and landfills in the shallow buried side, which could not transmit the unsymmetrical pressure from the deep buried side to the antislide piles.
(a)
(b)
(c)
(a)
(b)
(c)
4.2. Elements and Boundary Conditions
In order to analyse the moisture degradation effect on creep characteristics of surrounding rocks, a threedimensional finite difference analysis was performed by using the program . The analysis model includes information from the horizontal direction, 5.5 times the tunnel span from the tunnel centre to the model boundary, and from the vertical direction, 3 times the tunnel height from the tunnel bottom to the base boundary, and from the height above the tunnel from the tunnel crown to the nature ground surface. The model dimensions (, , and ) were approximately 143 m, 71 m, and 72 m. The boundary conditions include a free ground surface, fixed left and right boundaries in the axis direction, fixed front and rear boundaries in the axis direction, and a fixed base boundary in the axis direction. The tunnel liners, steel pipeline shed, horizontal support, and median septum are modelled with shell elements, and the antislide piles are modelled with solid elements. There are 63,519 elements in this numerical model (see Figure 13).
4.3. Simulation Phases and Rock Mass Physical and Mechanical Properties
Before excavation, the system was brought to elastic equilibrium (time 0 of modelling) under gravity. Boundary and initial conditions were set. Referring to Figure 9, the excavation sequence can be defined for the CRD method as follows.
Stage 1. Part A was excavated once with a length of 2 m.
Stage 2. Parts A and B were excavated once with lengths of 2 m.
Stage 3. Parts A, B, and C were excavated once with lengths of 2 m.
Stage 4. Parts A, B, C, and D were excavated once with lengths of 2 m, until the tunnel excavation was completed.
Based on the actual situation of tunnel excavation and support implementation, no matter which part was excavated, the lining and anchors corresponding to that part should be activated before the creep time of 12 hours had passed. And the next excavation stage was performed 24 hours after the support implementation.
For comparison, another condition, where the influences of rainfall and creep were not considered, was also simulated. The lining and anchors corresponding to the excavated part should be activated simultaneously to simulate the support of the tunnel.
The constitutive models used in the numerical calculation are as follows: elastic model for all of the support system, including the tunnel liners, steel pipeline shed, horizontal support, median septum, and antislide piles; the KelvinMC model established in Section 4.1 for surrounding rocks; and the MohrCoulomb model for surrounding rocks under the condition that the influences of rainfall and creep were not considered.
The evaluation of the rock mass quality has been carried out through detailed geomechanical field surveys according to the ISRM suggested method [29]. The geological strength indexes of the strongweathered slate and completely weathered slate are 25 and 10, respectively. The results obtained by field surveys and laboratory tests allowed us to apply the HoekBrown failure criterion and estimate the rock mass strength and elastic properties [30]. The decreases in the rock mass uniaxial compression strength were used to estimate the creep properties of rock mass [31]. The rock mass and support physical and mechanical properties are listed in Table 3.

During excavation, a heavy rainfall lasted for 3 days. Based on the relationship between saturation degree and time submerged in water (Section 2.2.1), the saturation degree of rock mass was estimated to be 84.0%.
4.4. Simulation Results
The excavation section ZK49+883, which was in the middle portion of the freshly constructed tunnel, was chosen as the subject for the research. Large surface subsidence happened after ZK49+883 was excavated (see Figure 14); the maximum value for the vertical displacement was −791 mm. The unsymmetrical surface subsidence squeezed the tunnel support structure, which made downward vertical displacements on the deep buried side and upward vertical displacements in the arch foot and haunch of shallow buried side. The maximum value for the horizontal displacement of surrounding rocks near the tunnel was 36 mm, which lies in the haunch of the deep buried side. The horizontal displacement of landfills behind the antislide piles was also very large. The simulation results show an acceptable agreement with the real situation.
(a)
(b)
The vertical and horizontal displacements at section ZK49+883, where the influences of rainfall and creep were not considered, are shown in Figure 15. The maximum values for the vertical and horizontal displacements were 7.2 mm and 2.1 mm, respectively. The low displacements indicate that the tunnel excavation is safe under this condition.
(a)
(b)
5. Conclusions
In this study, laboratory experiments were performed to understand the compression and creep characteristics of the lowgrade metamorphic slate with different moisture states. Then a threedimensional analysis on the Ruipo tunnel portal was performed while considering tunnel excavation and the moisture degradation effect on creep characteristics of surrounding rocks. The laboratory tests showed that the moisture has a great effect on the strength and deformability of the lowgrade metamorphic slate, and the mathematical relation between the creep parameters and degree of saturation was established. The simulation results show that the tunnel excavation is safe without considering the influence of rainfall and creep, whereas with these being considered large surface subsidence and displacement would appear. Therefore, a threedimensional analysis that considers the influence of moisture and creep of the surrounding rocks is needed for estimating the deformation of complex tunnel portals.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (no. 51408464), the Foundation of Shaanxi Educational Committee (no. 14JK1413), and the Foundation for the Talents of Xian University of Architecture and Technology (no. RC1365).
References
 L. S. Burshtein, “Effect of moisture on the strength and deformability of sandstone,” Soviet Mining Science, vol. 5, no. 5, pp. 573–576, 1969. View at: Publisher Site  Google Scholar
 E. M. van Eeckhout and S. S. Peng, “The effect of humidity on the compliances of coal mine shales,” International Journal of Rock Mechanics and Mining Sciences and, vol. 12, no. 11, pp. 335–340, 1975. View at: Publisher Site  Google Scholar
 A. B. Hawkins and B. J. McConnell, “Sensitivity of sandstone strength and deformability to changes in moisture content,” Quarterly Journal of Engineering Geology, vol. 25, no. 2, pp. 115–130, 1992. View at: Publisher Site  Google Scholar
 C. G. Dyke and L. Dobereiner, “Evaluating the strength and deformability of sandstones,” Quarterly Journal of Engineering Geology, vol. 24, no. 1, pp. 123–134, 1991. View at: Publisher Site  Google Scholar
 B. K. Atkinson and P. G. Meredith, “The theory of subcritical crack growth with applications to mineral and rocks,” in Fracture Mechanics of Rock, Academic Press, London, UK, 1987. View at: Google Scholar
 A. Golshani, M. Oda, Y. Okui, T. Takemura, and E. Munkhtogoo, “Numerical simulation of the excavation damaged zone around an opening in brittle rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 44, no. 6, pp. 835–845, 2007. View at: Publisher Site  Google Scholar
 E. Z. Lajtai, R. H. Schmidtke, and L. P. Bielus, “The effect of water on the timedependent deformation and fracture of a granite,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 24, no. 4, pp. 247–255, 1987. View at: Publisher Site  Google Scholar
 P. G. Meredith and B. K. Atkinson, “Stress corrosion and acoustic emission during tensile crack propagation in Whin Sill dolerite and other basic rocks,” Geophysical Journal International, vol. 75, no. 1, pp. 1–21, 1983. View at: Publisher Site  Google Scholar
 P. A. Vermeer, S. C. Möller, and N. Ruse, “On the application of numerical analysis in tunnelling,” in Proceedings of the 12th Asian Regional Conference on Soil Mechanics and Geotechnical Engineering, Singapore, August 2003. View at: Google Scholar
 M. Karakus and R. J. Fowell, “An insight into the New Austrian Tunneling Method (NATM),” in Proceedings of the 7th Regional Rock Mechanics Symposium (ROCKMEC '04), Sivas, Turkey, October 2004. View at: Google Scholar
 H. Gercek, “Underground space use: analysis of the past and lessons for the future,” in ThreeDimensional Regions in Tunnels, S. Erdem, Ed., Taylor & Francis Group, London, UK, 2005. View at: Google Scholar
 Z. Z. Wang, B. Gao, Y. J. Jiang, and S. Yuan, “Investigation and assessment on mountain tunnels and geotechnical damage after the Wenchuan earthquake,” Science in China Series E: Technological Sciences, vol. 52, no. 2, pp. 546–558, 2009. View at: Publisher Site  Google Scholar
 T. B. Li, “Damage to mountain tunnels related to the Wenchuan earthquake and some suggestions for aseismic tunnel construction,” Bulletin of Engineering Geology and the Environment, vol. 71, no. 2, pp. 297–308, 2012. View at: Publisher Site  Google Scholar
 F. Ye, C. He, S.M. Wang, and J.L. Zhang, “Landscape design of mountain highway tunnel portals in China,” Tunnelling and Underground Space Technology, vol. 29, pp. 52–68, 2012. View at: Publisher Site  Google Scholar
 Y.W. Pan and J.J. Dong, “Timedependent tunnel convergence—I. Formulation of the model,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 28, no. 6, pp. 469–475, 1991. View at: Publisher Site  Google Scholar
 C. E. Fairhurst and J. A. Hudson, “Draft ISRM suggested method for the complete stressstrain curve for intact rock in uniaxial compression,” International Journal of Rock Mechanics and Mining Sciences, vol. 36, no. 3, pp. 281–289, 1999. View at: Google Scholar
 J. MillsBeale, Z. You, R. C. Williams, and Q. L. Dai, “Determining the specific gravities of coarse aggregates utilizing vacuum saturation approach,” Construction and Building Materials, vol. 23, no. 3, pp. 1316–1322, 2009. View at: Publisher Site  Google Scholar
 Y. S. Li and C. C. Xia, “Timedependent tests on intact rocks in uniaxial compression,” International Journal of Rock Mechanics and Mining Sciences, vol. 37, no. 3, pp. 467–475, 2000. View at: Publisher Site  Google Scholar
 S.Q. Yang and H.W. Jing, “Strength failure and crack coalescence behavior of brittle sandstone samples containing a single fissure under uniaxial compression,” International Journal of Fracture, vol. 168, no. 2, pp. 227–250, 2011. View at: Publisher Site  Google Scholar
 C. Y. Li, H. L. Fu, H. L. Cai, X. N. She, and C. L. Nie, “Water character of flowering sheet stone,” Journal of Railway Science and Engineering, vol. 6, no. 1, pp. 74–77, 2009 (Chinese). View at: Google Scholar
 H. Lin, J. G. Deng, L. B. Hu, J. L. Yuan, Z. J. Chen, and X. D. Lai, “Study on the effect of water content on rock strength and sand production,” Science Technology and Engineering, vol. 13, no. 13, pp. 3710–3713, 2013 (Chinese). View at: Google Scholar
 W.D. Yang, Q.Y. Zhang, S.C. Li, and S.G. Wang, “Estimation of in situ viscoelastic parameters of a weak rock layer by timedependent plateloading tests,” International Journal of Rock Mechanics and Mining Sciences, vol. 66, pp. 169–176, 2014. View at: Publisher Site  Google Scholar
 W. Korzeniowski, “Rheological model of hard rock pillar,” Rock Mechanics and Rock Engineering, vol. 24, no. 3, pp. 155–166, 1991. View at: Publisher Site  Google Scholar
 X. R. Liu, X. Yang, and J. B. Wang, “A nonlinear creep model of rock salt and its numerical implement in ${\text{FLAC}}^{3\text{D}}$,” Advances in Materials Science and Engineering, vol. 2015, Article ID 285158, 8 pages, 2015. View at: Publisher Site  Google Scholar
 S. Q. Yang, P. Xu, P. G. Ranjith, G. F. Chen, and H. Jing, “Evaluation of creep mechanical behavior of deepburied marble under triaxial cyclic loading,” Arabian Journal of Geosciences, vol. 8, no. 9, pp. 6567–6582, 2015. View at: Publisher Site  Google Scholar
 X. J. Liu and Y. X. Zhang, “Analysis of reasonable excavation sequence and stress characteristics of portal section of shallow tunnel with unsymmetrical loadings,” Chinese Journal of Rock Mechanics and Engineering, vol. 30, supplement 1, pp. 3066–3073, 2011 (Chinese). View at: Google Scholar
 F. C. Xue, J. L. Ma, L. P. Yan, and Y. M. Zhao, “Threedimension FEM analysis of large crosssection tunnel in collapsible loess constructed by CRD method,” in Proceedings of the GeoFlorida, Advances in Analysis, Modeling, and Design, pp. 2349–2358, ASCE, West Palm Beach, Fla, USA, February 2010. View at: Google Scholar
 H. J. Li, “Study on construction technology for tunneling in the plateau permafrost region,” in Proceedings of the 14th Conference on Cold Regions Engineering, pp. 16–21, ASCE, Duluth, Minn, USA, AugustSeptember 2009. View at: Google Scholar
 ISRM International Society for Rock Mechanics, Rock Characterization, Testing and Monitoring. ISRM Suggest Methods, Pergamon Press, London, UK, 1981.
 E. Hoek, C. CarranzaTorres, and B. Corkum, “The HoekBrown failure criterion—2002 edition,” in Proceedings of the 5th North American Rock Mechanics Symposium and 17th Tunneling Association of Canada Conference (NARMSTAC '02), pp. 267–271, Toronto, Canada, 2002. View at: Google Scholar
 M. Huang, Study on the creep properties of waterbearing siltite and its application in soft rock tunnel engineering [Ph.D. thesis], Chongqing University, Chongqing, China, 2010 (Chinese).
Copyright
Copyright © 2015 Xiaojun Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.