abstract:
The Monte Carlo simulation (MCS) and First-Order Reliability Method (FORM) provide a reliability analysis in axisymmetric deep tunnels driven in elastoplastic rocks. The Convergence-Confinement method (CV-CF) and Mohr-Coulomb (M-C) criterion are used to model the mechanical interaction between the shotcrete lining and ground through deterministic parameters and random variables. Numerical models synchronize tunnel analytical models and reliability methods, whereas the limit state functions control the failure probability in both ground plastic zone and shotcrete lining. The results showed that a low dispersion of random variables affects the plastic zone's reliability analysis in unsupported tunnels. Moreover, the support pressure generates a significant reduction in the plastic zone's failure, whereas the increase of shotcrete thickness results in great reduction of the lining collapse probability.
Keywords:
axisymmetric deep tunnels; shotcrete lining; CV-CF method; MCS; FORM
resumo:
Este trabalho propõe uma análise de confiabilidade em túneis profundos axissimétricos, escavados em maciços rochosos elasto-plásticos, utilizando o Método de Confiabilidade de Primeira Ordem e simulações de Monte Carlo. O método Convergência-Confinamento (CV-CF) e o critério de Mohr-Coulomb (M-C) são utilizados para modelar analiticamente a interação mecânica entre o revestimento de concreto projetado e o maciço escavado, através de parâmetros determinísticos e variáveis aleatórias. Modelos numéricos são sincronizados com os modelos analíticos de túneis e os métodos de confiabilidade, enquanto as funções estado-limite definem a probabilidade de falha na zona plástica do maciço rochoso e no revestimento de concreto projetado. Os resultados mostram que uma baixa dispersão das variáveis aleatórias afeta a confiabilidade da zona plástica na análise de túneis sem revestimento. Além disso, a pressão de suporte gera uma redução significativa na falha da zona plástica, enquanto o aumento da espessura do concreto projetado resulta em uma grande redução da probabilidade de colapso do revestimento.
Palavras-chave:
túneis profundos axissimétricos; concreto projetado; método CV-CF; simulações de Monte Carlo; FORM
1 INTRODUCTION
Several geotechnical and structural parameters are involved in deep tunnels and underground excavation analysis, which induce high or undefined structural risks. The reliability analysis evaluates structural failure probabilities in tunnels by applying random and deterministic parameters in analytical solutions or numerical methods.
The first serious discussions and analyses of tunnel reliability emerged during the 1990s with [11 S. Kohno, “Reliability-based design of tunnel support systems,” Ph.D. dissertation, Univ. Illinois, Urbana, United States of America, 1989.]-[66 E. Hoek, "Reliability of Hoek-Brown estimates of rock mass properties and their impact on design," Int. J. Rock Mech. Min. Sci., vol. 35, no. 1, pp. 63–68, 1998.]. Low and Tang [55 B. K. Low and W. H. Tang, "Efficient reliability evaluation using spreadsheet," J. Eng. Mech., vol. 123, no. 1, pp. 749–752, 1997.] focused on identifying and evaluating the reliability and risk parameters through numerical procedures, synchronizing reliability codes developed in software like Microsoft Excel, and analytical methods for axisymmetric tunnels. This work used the First Order Reliability Method (FORM) and Monte Carlo Simulation (MCS) in the reliability study. Low and Tang [77 B. K. Low and W. H. Tang, "Reliability analysis using object-oriented constrained optimization," Struct. Saf., vol. 26, no. 1, pp. 69–89, 2004.], [88 B. K. Low and W. H. Tang, "Efficient spreadsheet algorithm for first-order reliability method," J. Eng. Mech., vol. 133, no. 1, pp. 1378–1387, 2007.] developed and improved the FORM numerical procedures. Li and Low [99 H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.] conducted a series of reliability analyses of an axisymmetric tunnel, considering normal and non-normal random variables of the rock mass and the Low and Tang [55 B. K. Low and W. H. Tang, "Efficient reliability evaluation using spreadsheet," J. Eng. Mech., vol. 123, no. 1, pp. 749–752, 1997.] FORM algorithm. The internal pressure produces a reduction of risk; moreover, the authors checked the reliability results using MCS. Lü and Low [1010 Q. Lü and B. K. Low, "Probabilistic analysis of underground rock excavations using response surface method and SORM," Comput. Geotech., vol. 38, no. 8, pp. 1008–1021, 2011.] and Song et al. [1111 L. Song, H.-Z. Li, C. L. Chan, and B. K. Low, "Reliability analysis of underground excavation in elastic-strain-softening rock mass," Tunn. Undergr. Space Technol., vol. 60, no. 1, pp. 66–79, 2016.] studied the reliability of axisymmetric tunnels through the FORM and Second-Order Reliability Method (SORM), synchronizing these methods with analytical solutions based on Mohr-Coulomb (M-C) and Hoek Brown criteria.
Bjureland et al. [1212 W. Bjureland, J. Spross, F. Johansson, A. Prästings, and S. Larsson, "Reliability aspects of rock tunnel design with the observational method," Int. J. Rock Mech. Min. Sci., vol. 1, no. 98, pp. 102–110, 2017.] applied the MCS associated with the Convergence-Confinement method (CV-CF method) to analyze the failure probability in a shotcrete lining of a circular tunnel. The M-C criterion describes the plasticity behavior for the excavated rock mass, and the Observational method checks the failure results, verifying the necessity of putting into action security measures.
The present paper analyses the reliability methods impact on analytical tunneling methods, examining two different tunnels with statistical and deterministic distinct parameters. This study aims to investigate the reliability results of the plastic zone excavation and the shotcrete lining for both tunnels, establishing interesting comparisons between them. For this, the reliability procedures associate the use of the FORM Low and Tang method, developed in Low and Tang [55 B. K. Low and W. H. Tang, "Efficient reliability evaluation using spreadsheet," J. Eng. Mech., vol. 123, no. 1, pp. 749–752, 1997.], and the MCS. The reliability algorithms, developed in VBA and MATLAB codes, employ two analytical methodologies for axisymmetric tunnels: the M-C plasticity criterion for excavated rock masses and the CV-CF method to analyze the interaction between the shotcrete lining and rock mass.
2 RELIABILITY CONCEPTS
2.1 Reliability Index: The Low and Tang FORM Interpretation
Hasofer and Lind [1313 A. M. Hasofer and N. C. Lind, "Exact and invariant second moment cod format," J. Eng. Mech., no. 100, pp. 111–121, 1974.] defined the second-moment reliability index β as a better approach in civil engineering design in comparison with usual safety factors [77 B. K. Low and W. H. Tang, "Reliability analysis using object-oriented constrained optimization," Struct. Saf., vol. 26, no. 1, pp. 69–89, 2004.]. Equation 1, proposed by Low and Tang [77 B. K. Low and W. H. Tang, "Reliability analysis using object-oriented constrained optimization," Struct. Saf., vol. 26, no. 1, pp. 69–89, 2004.], represents the matrix formulation of the Hasofer-Lind reliability index for correlated normal variables:
where is the set of random variables of the vector x, is the set of mean values of the vector , and is the set of standard deviations of the vector . The domain F represents the failure regions, whereas is the correlation matrix, expressed by [1414 A. K. Verma, D. R. Karanki, and S. Ajit, Reliability and Safety Engineering. London: Springer Verlag London Ltd., 2016.]:
where is the correlation coefficient between two random variables and .
For the non-normal distributions, the basic idea is to calculate the equivalent normal mean and equivalent standard deviation , which replaces the values of and in Equation 1. The Rackwitz-Fiessler transformation defines the following analytical solution for and [1515 A. S. Nowak and K. R. Collins, Reliability of Structures. New York: Mc Graw Hill, 2000.]:
In Equation 3, is the inverse of the cumulative distribution function (CDF) (standard normal), and is the original non-normal CDF, evaluated for x. In Equation 4, is the probability density function (pdf) of the standard normal distribution, whereas is the original non-normal pdf, evaluated for x.
In a classical explanation, is the smallest distance from the mean value point to the limit state surface in units of directional standard deviations [77 B. K. Low and W. H. Tang, "Reliability analysis using object-oriented constrained optimization," Struct. Saf., vol. 26, no. 1, pp. 69–89, 2004.]. Low and Tang [55 B. K. Low and W. H. Tang, "Efficient reliability evaluation using spreadsheet," J. Eng. Mech., vol. 123, no. 1, pp. 749–752, 1997.], [77 B. K. Low and W. H. Tang, "Reliability analysis using object-oriented constrained optimization," Struct. Saf., vol. 26, no. 1, pp. 69–89, 2004.] proposed an alternative interpretation of the reliability index, which is a simple method of defining the index, in the original space of the random variables, through the perspective of an expanding ellipsoid. Figure 1 shows this interpretation.
Equivalent dispersion ellipsoids illustrated in the plane [88 B. K. Low and W. H. Tang, "Efficient spreadsheet algorithm for first-order reliability method," J. Eng. Mech., vol. 133, no. 1, pp. 1378–1387, 2007.].
As shown in Figure 1, the design point is a tangency point of the expanding dispersion ellipsoid with the limit state surface, which separates safe values from unsafe values (failure domain). The reliability index is the axis ratio of the ellipse that touches the limit state surface [88 B. K. Low and W. H. Tang, "Efficient spreadsheet algorithm for first-order reliability method," J. Eng. Mech., vol. 133, no. 1, pp. 1378–1387, 2007.]. The dispersion ellipsoid expands from the mean-value point according to the following probability density function of a multivariate normal distribution [88 B. K. Low and W. H. Tang, "Efficient spreadsheet algorithm for first-order reliability method," J. Eng. Mech., vol. 133, no. 1, pp. 1378–1387, 2007.]:
where to minimize , defined by Equation 1, is to maximize the value of the probability density function defined by Equation 5. Hence, it is possible to define the smallest dispersion ellipsoid tangent to the limit state surface (see Figure1) and the most probable failure point (design point).
Low and Tang [55 B. K. Low and W. H. Tang, "Efficient reliability evaluation using spreadsheet," J. Eng. Mech., vol. 123, no. 1, pp. 749–752, 1997.] described that if each axis of the dispersion ellipsoid, in Figure 1, is parallel to the coordinate axis, the random variables are uncorrelated. So, in correlated normal random variables, each axis of the dispersion ellipsoid is inclined about the coordinate axis.
The index is particularly useful in evaluating the failure probability through the application of the CDF of , , and considering standard normal variables [99 H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.]:
Low and Tang [88 B. K. Low and W. H. Tang, "Efficient spreadsheet algorithm for first-order reliability method," J. Eng. Mech., vol. 133, no. 1, pp. 1378–1387, 2007.] performed a new efficient algorithm using FORM, where the first and third terms under the square root in Equation 1 are equivalent standard vectors . Thus, Equation 1 is defined as:
The optimization routine Solver of Microsoft Excel solves Equation 7. Autonomously, the routine changes the values of components (computing the parameters and , when necessary) and subjecting the limit state function to a zero-value restriction. The components in the original space, , are program-calculated from:
2.2 Monte Carlo Simulation
The Monte Carlo Simulation (MCS) consists of developing an analytical and numerical model and using many simulations or cycles to check the real behavior of a system. Each cycle provides output variables (results) based on random input variables, using statistical methods and probability distributions [1616 B. M. Ayyub and R. H. McCuen, Probability, Statistics, and Reliability for Engineers and Scientist, 1st ed. New York: Taylor & Francis Group, 2011.].
Considering a random variable X, with a CDF , to generate sample values included in X, the MCS code needs to accomplish some steps. First, MATLAB numerical procedures generate a sample value of between 0 and 1 for a random variable with uniform distribution. Next, the inverse of computes a sample value of ().
For two correlated normal random variables, the numerical procedure generates two random numbers and , used to compute the random variables x and y from the following equations [1717 A. H.-S. Ang and W. H. Tang, Probability Concepts in Engineering Planning and Design, vol. II: Decision, Risk, and Reliability, 1st ed. New York: John Wiley & Sons, Inc., 1984.]:
where the mean and standard deviation of the variable y, which are dependent on x value, are evaluated through the following equations:
Each MCS cycle consists of introducing random input variables in an appropriated mathematical model: the limit state or performance functions . A negative performance function implies random variables located in a failure domain. After the implementation of all cycles, the following relation, between the total failure cycles and the total cycles of MCS N, checks the failure probability [1515 A. S. Nowak and K. R. Collins, Reliability of Structures. New York: Mc Graw Hill, 2000.]:
Figure 2 shows a flowchart with the steps of the MCS numerical procedure developed in MATLAB. In summary, the numerical procedure considers a predefined number of cycles and checking the failure cycles, . Equation 14 checks the failure probability estimation error through a normal distribution approximation, with a reliability interval of 95%.
3 ANALYTICAL FRAMEWORK FOR ELASTOPLASTIC DEFORMATION IN AXISYMMETRIC TUNNELS
3.1 Elastoplastic fields in the surrounding rock mass
Figure 3 shows the axisymmetric tunnel characteristics, excavated in a homogeneous elastic-plastic rock mass and with a radius . The cylindrical coordinate system , , and represents the stress and displacement components. The stresses and are the in-plane normal components of the rock stress tensor , whereas and is the far-field pressure and uniform support pressure into the excavation. It is assumed that the stress component does not vary, which implicitly amounts to disregarding the confinement losses along the z-axis direction.
In the fully elastic regime, the radial displacement at any a point of the rock mass located at distance r from the center of the excavation is given by ([1818 J. Salençon, "Contraction quasi-statique d’une cavite a symetrie spherique ou cylindrique dans un milieu elasto-plastique," Ann. Ponts Chaussees, vol. 4, pp. 231–236, 1969.], [1919 F. Corbetta, “Nouvelles méthodes d’étude des tunnels profonds: calculs analytiques et numériques,” Ph.D. dissertation, Éc. Natl. Super. Mines, Paris, France, 1990.]):
where and are the Poisson's rate and Young's modulus, respectively.
In the elastoplastic regime, the mechanical analysis of the tunnel deformation indicates that two regions should be distinguished, as illustrated in Figure 4: the elastoplastic region undergoing irreversible strains, and the elastic behavior region . The plastic radius defines the limit between the elastic and elastoplastic regions of the rock mass.
The analytical methods use the M-C criterion concepts for the elastoplastic behavior, defined according to Salençon [1818 J. Salençon, "Contraction quasi-statique d’une cavite a symetrie spherique ou cylindrique dans un milieu elasto-plastique," Ann. Ponts Chaussees, vol. 4, pp. 231–236, 1969.] and Corbetta [1919 F. Corbetta, “Nouvelles méthodes d’étude des tunnels profonds: calculs analytiques et numériques,” Ph.D. dissertation, Éc. Natl. Super. Mines, Paris, France, 1990.]. The plastic radius is:
where:
in which c is the cohesion and is the friction angle. The limit pressure defines a rock mass pressure when as follows:
Referring to Figure 4, the radial displacement in the elastoplastic region is defined by:
where Equation 21 defines the dilatancy coefficient , representing the dilatancy angle, whereas Equations 22, 23 and 24 characterize the parameters A, B, and C.
3.2 Convergence-confinement method
The CV-CF method allows modeling in a simplified manner the evolution of normalized radial displacement of the rock mass and support, or the convergence at the tunnel wall, associated with the internal pressure in a specific cross-section of the tunnel wall. The convergence is given by:
where is the radial displacement at the tunnel wall, defined by Equation 15 (elastic behavior) or Equation 20 (elastoplastic behavior). The principle of the CV-CF method is schematized in Figure 5.
Schematized description of the Convergence-confinement method for a specific cross-section of a tunnel.
Figure 5 shows the Ground response curve, or Convergence curve (CV), which describes only the mechanical behavior of the unsupported excavation wall, while the Support curve, or Confinement curve (CF), describes the mechanical behavior of the elastic shotcrete lining.
When the excavation front is close to a specific cross-section, changes occur in the ground stress state. At some distance, before the excavation reaches the cross-section, there are small changes of the stress state, and elastic convergence occurs as decreases. The internal pressure continues to decrease where, eventually, the decrement of stresses in the rock mass will reach a limit , and the stresses cause plastic behavior of the ground in a zone with a radius surrounding the tunnel periphery [1212 W. Bjureland, J. Spross, F. Johansson, A. Prästings, and S. Larsson, "Reliability aspects of rock tunnel design with the observational method," Int. J. Rock Mech. Min. Sci., vol. 1, no. 98, pp. 102–110, 2017.]. Thus, considering an interval , the CV curve shows an elastic behavior (blue curve in Figure 5), and for the interval , the CV curve shows an elastoplastic behavior. Equations 16 and 19 defines and
The CF curve avoids the ground confinement loss, which causes normalized radial displacements of the rock mass. According to the red curve in Figure 5, a ground confinement loss generates a growing internal pressure and convergence on the cross-section lining. Moreover, numerical or analytical approaches estimate the convergence at the beginning of the support installation. Panet and Guenot [2020 M. Panet and A. Guenot, “Analysis of convergence behind the face of a tunnel”, in Tunnelling 82, Proc. 3rd Int. Symp., Brighton, 1982, pp. 197–204.] proposed an analytical approximation for :
where is the convergence for elastoplastic behavior, at a larger distance between the excavation face and the cross-section (Equation 25 with ).
The pressure at the shotcrete lining is evaluated as [2121 D. Bernaud, “Tunnels profonds dans les milieux viscoplastiques: approaches experimentale et numerique,” Ph.D. dissertation, L’Ec. Natl. Ponts et Chaussées, Paris, France, 1991.]:
In Equation 27, is the shotcrete convergence (determined when varies between 0 and in Equation 27), and is the support stiffness coefficient, calculated according to the relation ([2222 M. Panet and P. Guellec, “Contribution à l’etude du Soutènement d’un Tunnel à l’Arrière du Front de Taille,” in Proc. 3rd Int. Cong. Rock Mechanics, Denver, 1982.], [2323 E. Hoek and E. T. Brown, Underground Excavations in Rock, 1st ed. London: E&FN Spon, 1980.]):
where is the shotcrete lining thickness, while and are the Young modulus and the Poisson Coefficient of the support, respectively. The maximum support pressure is defined by the following Equation:
standing for the uniaxial compressive strength of the support lining.
The maximum convergence is evaluated by setting into Equation 27. Referring to Figure 5, the equilibrium state of the structure ( and ) is defined as the intersection point between the CV and CF curves. The design value can be computed combining Equations 16, 20 and 27:
By using Equation 16 one can calculate with . Then Equation 30 is solved numerically by the bisection method inside the MCS MATLAB code.
4 PERFORMANCE FUNCTIONS AND ADOPTED PARAMETERS
The limit state function, or performance function, refers to the plastic zone's reliability analyses, whereas refers to the shotcrete lining reliability. The functions and consider the relationship between an accepted resistance and a demand (loading), according to the equations below ([99 H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.], [1313 A. M. Hasofer and N. C. Lind, "Exact and invariant second moment cod format," J. Eng. Mech., no. 100, pp. 111–121, 1974.]):
The ratios L and are the allowed largest values (resistances) for the following demands (loadings): plastic zone ratio () and design convergence ratio (). Equation 31 evaluates using Equation 16, while Equation 32 evaluates by fixing in Equation 27, with . When takes a negative value, an unacceptable large plastic zone occurs, and it exceeds the L ratio. If takes a negative value, the radial displacement exceeds the largest support displacement .
The geotechnical random variables were obtained in the research of Hoek [66 E. Hoek, "Reliability of Hoek-Brown estimates of rock mass properties and their impact on design," Int. J. Rock Mech. Min. Sci., vol. 35, no. 1, pp. 63–68, 1998.] (tunnel A) and Brantmark [2424 J. Brantmark, “Rock support in weak rock – a case study based on the Uri Project,” Ph.D. dissertation, KTH Roy. Inst. Technol., Stockholm, Sweden, 1998.] (tunnel B). Tables 1 and 2 show the geotechnical statistical parameters for tunnels A and B.
Geotechnical statistical parameters for tunnel A [66 E. Hoek, "Reliability of Hoek-Brown estimates of rock mass properties and their impact on design," Int. J. Rock Mech. Min. Sci., vol. 35, no. 1, pp. 63–68, 1998.].
Geotechnical statistical parameters for tunnel B [2424 J. Brantmark, “Rock support in weak rock – a case study based on the Uri Project,” Ph.D. dissertation, KTH Roy. Inst. Technol., Stockholm, Sweden, 1998.].
All statistical parameters described in Table 1 and 2 are normally distributed. However, it is considered a negative correlation between c and , with the following correlation coefficients: for tunnel A and for tunnel B. It should be emphasized that the correlation between some of the design parameters significantly affects the reliability predictions. In that respect, the assumption of negatively correlated shear strength parameters was found to be conservative with respect to uncorrelated variables [2525 A. Hamrouni, D. Dias, and B. Sbartai, "Reliability analysis of shallow tunnels using the response surface methodology," Undergr. Space, vol. 1, no. 2, pp. 246–258, 2017.]. The analysis developed, for instance, in Laso et al. [44 E. Laso, M. S. G. Lera, and E. Alarcón, "A level II reliability approach to tunnel support design," Appl. Math. Model., vol. 19, no. 1, pp. 371–382, 1995.] has included the correlation between some of the relevant problem parameters, such as correlation between the ground Young modulus and ground density, shotcrete average resistance and shotcrete layer thickness, as well as between tunnel radius and shotcrete layer thickness. Although the subject of correlation is not fully addressed in the present reliability study, this important issue will be foreseen in the continuation line of this research.
The rock mass failure criterion features a non-associated flow rule in both tunnels, with the following deterministic parameters: and (Equation 21) for tunnel A; and for tunnel B. The Poisson coefficients are , and for tunnels A and B, respectively. For both tunnels, a radius of is adopted.
Table 3 shows the statistical parameters for the shotcrete lining. Both tunnels use the same parameters, which are statistically independent and normally distributed. The Poisson coefficient is .
Statistical parameters of shotcrete support [2424 J. Brantmark, “Rock support in weak rock – a case study based on the Uri Project,” Ph.D. dissertation, KTH Roy. Inst. Technol., Stockholm, Sweden, 1998.].
Li and Low [99 H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.] suggested applying the lognormal or beta distribution instead of the normal distribution when the coefficient of variation a random variable (ratio between standard deviation and mean) is 0.25 or higher. This is necessary to avoid negative geotechnical parameters, which is irrational. All coefficients of variation are smaller than 0.25 in tunnels A and B, except for the cohesion c in tunnel A, where, alternatively, the cohesion is constrained to be greater than zero.
5 RELIABILITY ANALYSIS REGARDING THE EXTENT OF PLASTIC ZONE
This section analyzes the reliability of tunnels A and B, contemplating the performance function of Equation 31. The first reliability analysis performs the characterization of the fixed resistances in regarding unsupported tunnels, whereas the second analysis of the reliability methods will contemplate the shotcrete lining. The 2007 FORM algorithm of Low and Tang and MCS compute the reliability indexes and failure probabilities, employing the analytical theories explained in section 2.
In summary, the Low and Tang algorithm is a code developed in Visual Basic (VBA) language and Microsoft Excel. A constrained optimization routine, denominated Solver, assists in the FORM procedure development, where the equations and parameters described in sections 2.1 and 4 have been employed directly in the Excel worksheets. Low and Tang [77 B. K. Low and W. H. Tang, "Reliability analysis using object-oriented constrained optimization," Struct. Saf., vol. 26, no. 1, pp. 69–89, 2004.], [88 B. K. Low and W. H. Tang, "Efficient spreadsheet algorithm for first-order reliability method," J. Eng. Mech., vol. 133, no. 1, pp. 1378–1387, 2007.], and Li and Low [99 H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.] give more details about the FORM algorithm of Low and Tang.
First, it is necessary to verify the results obtained by FORM and MCS to check the algorithms' proper operation. Li and Low [99 H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.] analyzed the plastic zone reliability, employing the FORM Low and Tang code, the performance function , and tunnel A parameters. Thus, Table 4 shows a comparison between results of this work and Li and Low [99 H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.], with the following deterministic values: the parameters of tunnel A, , , and 5E+06 Monte Carlo simulations, which generates an error between when .
The failure probabilities evaluated in Tab. 4 have shown a good agreement with the Li and Low [99 H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.] results, considering both FORM and MCS results of this work. Thus, considering random variables with normal distribution, the reliability algorithms used in this work have shown efficiency when applied in unsupported and supported tunnels analysis.
The first analysis intends to verify the influence of resistance L in the reliability index , on the plastic zone reliability in tunnels A and B. Figure 6 shows the values evaluated for the statistical parameters of unsupported tunnels A and B, employing FORM Low and Tang algorithm, and varying the resistance L, in Equation 31, considering the following intervals: for tunnel A and for tunnel B.
Reliability index results with the variation of the resistance L, regarding unsupported tunnels A and B ().
Comparing the curves in Figure 6, a great sensibility on the variation with the increase of L occurs for tunnel B (red curve). For example, the index in tunnel B increases from to if the resistance increases from to . On the other hand, tunnel A (blue curve) does not show a considerable increase in , i.e., the index increases from to when increases from to . Table 5 illustrates the values of and according to the curves of Figure 6.
The results of Tab. 5 highlight the reduction in both studies. In tunnel A, decreases from 44.45% to 25.86% between and , while in tunnel B decreases from 18.85% to 8.57E-07% between and . In contrast, the reliability index increases from 0.14 to 0.65 for tunnel A and from 0.88 to 5.64 for tunnel B, with L variation. The greater variation of reliability results in tunnel B is a consequence of the correlation between the random variables and the plastic zone results , performed by MCS.
Figure 7 illustrates the correlation between the random variables and through dispersion charts, using 30E+03 Monte Carlo simulations. The number of simulations was varied to assess the accuracy of the MCS results. The errors in estimating the failure probability were 1.28% and 2.40%, checked through Equation 14, using the parameters of unsupported tunnels A and B ().
Dispersion results of the plastic zone () regarding cohesion (c), friction angle (), and far-field stress (): (a) Tunnel A; (b) Tunnel B.
There is a strong correlation between the plastic zone and the parameters (, and in tunnel B (Figures 7b). It does not occur with tunnel A parameters, which are more dispersed than tunnel B parameters (Figures 7a). About these figures, when the resistance L is close to the plastic zone average, the results are more dispersed in the failure zone . Due to this, the reliability results of tunnel A, in Figure 6 and Table 5, shows a low variation when adopting . In contrast, when for tunnel B, the strong correlation causes a little dispersion in the failure zone , and a great variation in reliability results when adopting
In an overall analysis, these results say that the designers must be careful when adopting elevated values for L, in reliability analysis that involves highly correlated parameters. For this purpose, this study adopted the average resistance values in the performance function : for tunnel B, and, for comparison purposes, for tunnel A.
Turning now to the shotcrete lining tunnels, the Low and Tang FORM algorithm and MCS estimates the plastic zone reliability through Equation 31. Tables 6 and 7 summarizes the and results of the plastic zone. The following interval constraints the shotcrete lining pressure : . The MCS results for and consider 30E+03 simulations in both tunnels, whereas the MCS results for and consider 20E+06 simulations in both tunnels. Regarding , the MCS results refer to 10E+07 simulations for tunnel A and 20E+06 simulations for tunnel B. The number of simulations increases with to keep the error below 10%.
The results in Tables 6 and 7 show that is greater in tunnel A than in tunnel B, when the support pressure , whereas the opposite occurs when . This shows the importance of adopting an appropriate lining for tunnel safe. The and results obtained by the Low and Tang Algorithm and MCS show good concordance. The variation in both tunnels reported reductions of and increases in .
The MCS code provides the plastic zone histograms in Figures 8a and 8b, where the MATLAB application Distribution Fitter adjusts probability distribution functions with these histograms.
Probability Distribution Functions (PDF) of the plastic zone through Monte Carlo simulations: (a) Tunnel A; (b) Tunnel B.
The plastic zone intervals of PDF curves of Tunnel B (Figure 8b) is smaller than the intervals of tunnel A (Figure 8a). It occurs because the standard deviations for tunnel B are smaller than the values of tunnel A. Also, the mean values of tunnel B are smaller than of tunnel A, due to the higher geotechnical parameters of tunnel B.
Figures 8a and 8b emphasize the histogram distribution types. Tunnel A fits the Lognormal distribution for all histograms, whereas tunnel B fits the Normal distribution for all histograms.
All PDF curves exhibit reductions in the statistical parameters when increases. In Figure 8a, the mean and standard deviation reduces from and to and , whereas in Figure 8b, the parameters reduce from and to and . To summarize, the increasing causes the reliability index increasing and failure probability reduction during the reliability analyses.
6 RELIABILITY ANALYSIS REGARDING THE TUNNEL CONVERGENCE
In this section, the reliability analysis evaluates the probability of and exceeding the allowed parameters and . The MCS employs the CV-CF method linked with the performance function (Equation 32). The shotcrete lining parameters are those described in section 4 and Table 3.
As evaluated in the previous section, Equation 14 checks the MCS accuracy considering some number of simulations and a shotcrete lining thickness . This work adopts 5E+06 simulations, which generate the following results for tunnels A and B: with an error of 0.43% and with an error of 0.18%.
Figure 9 shows the variation of the failure probability with the increase of . The results were obtained using 5E+06 MCS cycles for each , where in tunnel A the error varies between regarding , and in tunnel B the error varies between regarding . The increase of produces a considerable reduction in tunnel A. The increase of just one centimeter in of tunnel A reduces from 4.06% to 0.93% , arrives at 0.02% when . With the tunnel B parameters, Figure 9 shows a reduction from to , arrives at 0.08% when . These results say that tunnel B needs more shotcrete than tunnel A to get acceptable reliability levels due to the far-field pressure in tunnel B is much higher than in tunnel A.
Failure probability reduction with the increasing of , concerning the performance function .
Table 8 exhibits the design and maximum support parameters, considering for tunnel A and for tunnel B, which keep an adequate structural safety level for the support . The convergences and have been calculated using the convergence of the support installation as the origin, whereas the radial displacements and have been calculated through Equation 25, i.e., and .
The results in Table 8 reveal a difference of 72% between the mean values for tunnels A and B, whereas the difference between the mean values is only 9%. Concerning and , the thickness controls the support resistance, which increases the stiffness value and , but does not change in both tunnels.
Figures 10a and 10b illustrate the failure probability through the and histograms obtained with the parameters of Table 8.
Demand and resistance histograms, generated using 5E+06 Monte Carlo simulations: (a) Tunnel A; (b) Tunnel B.
The histograms in Figures 10a and 10b show the support failure region, between 6.00 and 7.00 mm approximately, where . The area of this region is proportional to the failure probability, . From Figure 9, it can be seen that on (tunnel A) and (tunnel B) have close values. This occurs due to standard deviations since the histogram in Figure 10a is sparser than the corresponding histogram in Figure 10b. As well as in Figure 10, it is easy to evaluate the failure probability through , according to Figures 11a and 11b.
Histograms of performance function , generated using Monte Carlo simulations: (a) Tunnel A; (b) Tunnel B.
The negative values of , located on the left side of the dashed line in Figures 11a and 11b represents a failure point of the shotcrete lining. That is to say that the area below the histograms and on the left side of the dashed line (in Figure 11) has the same value of the failure probability, . It is possible to see in Figures 11a and 11b, that this area reduces with the increase of , and, therefore, also reduces . Similar to earlier findings, the standard deviation in tunnel A is greater than tunnel B in histograms of Figure 11.
Figures 12a and 12b show an analysis through the CV-CF method, using the mean values of tunnels A and B to verify the accordance with the results defined in Table 8. In these cases, it has been considered and for tunnels A and B, respectively. The results in Figure 12 must be interpreted with caution due to referring to the shotcrete lining displacement, which starts from the displacement .
There are similarities between the design parameters defined in Figures 12a and 12b and the design parameters of Table 8, where it can be seen that the interaction between the Convergence and Confinement curves has been influencing the reliability analysis. The Confinement curve of tunnel A tends to perform better than tunnel B due to the internal pressure being higher in tunnel B.
These results suggest that the shotcrete lining of tunnel A, considering , produces acceptable stiffness results , which increases the inclination of the Convergence curve in Figure 12. Tunnel B, likewise, produces acceptable results using . The designer must consider adopting support parameters higher than the before adopted, aiming to increase , thus keeping the failure probability low.
7 SUMMARY AND CONCLUSIONS
This study performed a reliability analysis in deep axisymmetric tunnels, integrating FORM and MCS methods with tunnel analytical methodologies. VBA and MATLAB codes verified the reliability in rock masses and shotcrete lining, examining two different tunnels.
The first aim was to check the plastic zone reliability through VBA FORM and MATLAB MCS algorithms. The reliability results verification has shown an excellent agreement between the reliability indexes and failure probabilities, obtained by both algorithms, and with the Li and Low [99 H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.] results.
The reliability study of the variation in the resistance L, considering unsupported tunnels, suggests that, in general, it is necessary to pay attention to the adopted L value in rock masses with low dispersion. This work recommends adopting a resistance value L around the mean value to avoid higher values.
The increase of internal pressure produces great changes in the plastic zone reliability. Both tunnels develop significant reductions when : reduces 43.80% in tunnel A and 18.40% in tunnel B. The support pressure maintains reliability percentages above 99% in both tunnels; moreover, for , results are close to 0%.
The MCS code used the CV-CF method equations to provide the convergence reliability. A shotcrete thickness results in tunnel A, whereas results in tunnel B. When increases for 13 cm in tunnel A, the failure probability reduces to ; whereas for tunnel B, reduces to 0.02% when increases to 20 cm. These findings suggest that tunnel A must have a shotcrete thickness between and . Tunnel B must have a shotcrete thickness between and . In summary, the failure probabilities are smaller than 0.1%, considering in tunnel A and for tunnel B. Moreover, the internal design pressure for these thicknesses is greater than , which keeps plastic zone failures close to zero.
This paper contributes to check safety or failure probabilities in tunnel design stages. The reliability codes provided statistical data for different rock mass parameters. So, designers and engineers can use these methodologies to ensure security in underground design.
A natural progression of this work is to link two-dimensional and three-dimensional numerical models with the reliability codes. Moreover, this study suggests using reliability analysis in different cross-section geometries and supports (steel lining, precast concrete lining, bolts).
ACKNOWLEDGEMENTS
The authors acknowledge UFRGS and FURG for their support and CNPq and CAPES for their research grants.
-
Financial support: This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brazil (CAPES) - Finance Code 001.
-
How to cite: C. C. C. Silva, M. V. Real, and S. Maghous, “A simplified approach to reliability evaluation of deep rock tunnel deformation using First-Order Reliability Method and Monte Carlo simulations,” Rev. IBRACON Estrut. Mater., vol. 15, no. 1, e15104, 2022, https://doi.org/10.1590/S1983-41952022000100004
CITATIONS
-
1S. Kohno, “Reliability-based design of tunnel support systems,” Ph.D. dissertation, Univ. Illinois, Urbana, United States of America, 1989.
-
2E. Laso, “Application of Level II reliability methods to the design of tunnels liners,” Ph.D. dissertation, Univ. Politéc. Madrid, Madrid, Spain, 1991. [in Spanish].
-
3S. Kohno, A. H.-S. Ang, and W. H. Tang, "Reliability evaluation of idealized tunnel system," Struct. Saf., vol. 11, no. 1, pp. 81–93, 1992.
-
4E. Laso, M. S. G. Lera, and E. Alarcón, "A level II reliability approach to tunnel support design," Appl. Math. Model., vol. 19, no. 1, pp. 371–382, 1995.
-
5B. K. Low and W. H. Tang, "Efficient reliability evaluation using spreadsheet," J. Eng. Mech., vol. 123, no. 1, pp. 749–752, 1997.
-
6E. Hoek, "Reliability of Hoek-Brown estimates of rock mass properties and their impact on design," Int. J. Rock Mech. Min. Sci., vol. 35, no. 1, pp. 63–68, 1998.
-
7B. K. Low and W. H. Tang, "Reliability analysis using object-oriented constrained optimization," Struct. Saf., vol. 26, no. 1, pp. 69–89, 2004.
-
8B. K. Low and W. H. Tang, "Efficient spreadsheet algorithm for first-order reliability method," J. Eng. Mech., vol. 133, no. 1, pp. 1378–1387, 2007.
-
9H. Li and B. K. Low, "Reliability analysis of circular tunnel under hydrostatic stress field," Comput. Geotech., vol. 37, no. 1-2, pp. 50–58, 2010.
-
10Q. Lü and B. K. Low, "Probabilistic analysis of underground rock excavations using response surface method and SORM," Comput. Geotech., vol. 38, no. 8, pp. 1008–1021, 2011.
-
11L. Song, H.-Z. Li, C. L. Chan, and B. K. Low, "Reliability analysis of underground excavation in elastic-strain-softening rock mass," Tunn. Undergr. Space Technol., vol. 60, no. 1, pp. 66–79, 2016.
-
12W. Bjureland, J. Spross, F. Johansson, A. Prästings, and S. Larsson, "Reliability aspects of rock tunnel design with the observational method," Int. J. Rock Mech. Min. Sci., vol. 1, no. 98, pp. 102–110, 2017.
-
13A. M. Hasofer and N. C. Lind, "Exact and invariant second moment cod format," J. Eng. Mech., no. 100, pp. 111–121, 1974.
-
14A. K. Verma, D. R. Karanki, and S. Ajit, Reliability and Safety Engineering London: Springer Verlag London Ltd., 2016.
-
15A. S. Nowak and K. R. Collins, Reliability of Structures. New York: Mc Graw Hill, 2000.
-
16B. M. Ayyub and R. H. McCuen, Probability, Statistics, and Reliability for Engineers and Scientist, 1st ed. New York: Taylor & Francis Group, 2011.
-
17A. H.-S. Ang and W. H. Tang, Probability Concepts in Engineering Planning and Design, vol. II: Decision, Risk, and Reliability, 1st ed. New York: John Wiley & Sons, Inc., 1984.
-
18J. Salençon, "Contraction quasi-statique d’une cavite a symetrie spherique ou cylindrique dans un milieu elasto-plastique," Ann. Ponts Chaussees, vol. 4, pp. 231–236, 1969.
-
19F. Corbetta, “Nouvelles méthodes d’étude des tunnels profonds: calculs analytiques et numériques,” Ph.D. dissertation, Éc. Natl. Super. Mines, Paris, France, 1990.
-
20M. Panet and A. Guenot, “Analysis of convergence behind the face of a tunnel”, in Tunnelling 82, Proc. 3rd Int. Symp., Brighton, 1982, pp. 197–204.
-
21D. Bernaud, “Tunnels profonds dans les milieux viscoplastiques: approaches experimentale et numerique,” Ph.D. dissertation, L’Ec. Natl. Ponts et Chaussées, Paris, France, 1991.
-
22M. Panet and P. Guellec, “Contribution à l’etude du Soutènement d’un Tunnel à l’Arrière du Front de Taille,” in Proc. 3rd Int. Cong. Rock Mechanics, Denver, 1982.
-
23E. Hoek and E. T. Brown, Underground Excavations in Rock, 1st ed. London: E&FN Spon, 1980.
-
24J. Brantmark, “Rock support in weak rock – a case study based on the Uri Project,” Ph.D. dissertation, KTH Roy. Inst. Technol., Stockholm, Sweden, 1998.
-
25A. Hamrouni, D. Dias, and B. Sbartai, "Reliability analysis of shallow tunnels using the response surface methodology," Undergr. Space, vol. 1, no. 2, pp. 246–258, 2017.
Edited by
Publication Dates
-
Publication in this collection
28 July 2021 -
Date of issue
2022
History
-
Received
13 Mar 2021 -
Accepted
19 Apr 2021