A new closed‐form solution for analysis of unlined pressure tunnels under seepage forces

In this study, a simplified analytical closed‐form solution, considering plane strain and axial symmetry conditions, for analysis of a circular pressure tunnel excavated underwater table, is developed. The method accounts for the seepage forces with the steady‐state flow and is based on the generalized effective stress law. To examine the effect of pore pressure variations and also the boundary conditions at the ground surface, the formulations are derived for different directions around the tunnel. The proposed method can be applied for analysis and design of pressure tunnels. Illustrative examples are given to demonstrate the performance of the proposed solution and also to examine the effect of seepage forces on the stability of tunnels. The simplified analytical solution derived in this study is compared with numerical analyses. It is concluded that the classic solutions (Lame's thick‐walled solution), considering the internal pressure as a mechanical load applied to the tunnel surface, are not applicable to pervious media and can result in an unsafe design. Copyright © 2012 John Wiley & Sons, Ltd.     

Analytical stress and displacement due to twin tunneling in an elastic semi‐infinite ground subjected to surcharge loads

The construction of twin tunnels at shallow depth has become increasingly common in urban areas. In general, twin tunnels are usually near each other, in which the interaction between tunnels is too significant to be ignored on their stability. The equivalent arbitrarily distributed loads imposed on ground surface were considered in this study, and a new analytical approach was provided to efficiently predict the elastic stresses and displacements around the twin tunnels. The interaction between 2 tunnels of different radii with various arrangements was taken into account in the analysis. We used the Schwartz alternating method in this study to reduce the twin‐tunnel problem to a series of problems where only 1 tunnel was contained in half‐plane. The convergent and highly accurate analytical solutions were achieved by superposing the solutions of the reduced single‐tunnel problems. The analytical solutions were then verified by the good agreement between analytical and numerical results. Furthermore, by the comparison on initial plastic zone and surface settlement between analytical solution and numerical/measured results of elastoplastic cases, it was proven that the analytical solution can accurately predict the initial plastic zone and its propagation direction and can qualitatively provide the reliable ground settlements. A parametric study was finally performed to investigate the influence of locations of surcharge load and the tunnel arrangement on the ground stresses and displacements. The new solution proposed in this study provides an insight into the interaction of shallow twin tunnels under surcharge loads, and it can be used as an alternative approach for the preliminary design of future shallow tunnels excavated in rock or medium/stiff clay.     

A complex variable solution for a deforming buoyant tunnel in a heavy elastic half‐plane

An analytical solution is presented for a buoyant tunnel in an elastic half‐plane. The tunnel undergoes a prescribed displacement along its boundary, and the surface of the half‐plane is stress‐free. The solution uses the complex variable method and consists of conformally mapping the hole and half‐plane to an annular region. The complex potentials each contain two logarithmic terms due to the resultant buoyancy force acting on the tunnel, and as a result the displacements at infinity are unbounded. An example is presented for the case of a rigid buoyant tunnel in a half‐plane with gravity loading. Even though the tunnel does not deform, the buoyancy effect can be clearly seen in the contours of the displacements and stresses. Copyright © 2002 John Wiley & Sons, Ltd.     

基于复变函数方法的水下隧道围岩弹性分析

水下隧道属于浅埋隧道,故其围岩的应力位移分析从力学角度可简化为含圆孔半平面弹性体在水平边界上受任意分布荷载的问题。与深埋隧道不同,浅埋隧道围岩分析在数学处理上历来存在较大的困难。借助Verruijt提供的共形映射函数,把含圆孔半无限平面映射为圆环域,然后将像平面上的解析函数展成Laurent级数,利用Muskhelishvili 的复变函数解法,求得问题的应力场和位移场,最后用Fortran语言编写了计算程序。利用该程序给出一个水下隧道算例的围岩应力、位移结果,并分析了其受力变形特点。     

An analytical solution for deforming twin‐parallel tunnels in an elastic half plane

The interaction between twin‐parallel tunnels affects the tunnelling‐induced ground deformation, which may endanger the nearby structures. In this paper, an analytical solution is presented for problems in determining displacements and stresses around deforming twin‐parallel tunnels in an elastic half plane, on the basis of complex variable theory. As an example, a uniform radial displacement was assumed as the boundary condition for each of the two tunnels. Special attention was paid to the effects of tunnel depth and spacing between the two tunnels on the surface movement to gain deep insight into the effect of the interaction between twin‐parallel tunnels using the proposed analytical approach. It is revealed that the influence of twin tunnel interaction on surface movements diminishes with both the increase of the tunnel depth and the spacing between the two tunnels. The presented analytical solution manifests that, similar to most of the existing numerical results, the principle of superposition can be applied to determine ground deformation of twin‐parallel tunnels with a certain large depth and spacing; otherwise, the interaction effect between the two tunnels should be taken into account for predicting reliable ground movement. Copyright © 2014 John Wiley & Sons, Ltd.     

Elasto‐plastic seepage‐induced stresses due to tunneling

When tunneling is carried out beneath the groundwater table, hydraulic boundary is altered, resulting in seepage entering into the tunnel. The development of flow into the tunnel induces seepage stresses in the ground and the lining is subjected to additional loads. This can often cause fine particles to move, which clog the filter resulting in the long‐term hydraulic deterioration of the drainage system. However, the effect of seepage force is generally not considered in the analysis of tunnel. While several elastic solutions have been proposed by assuming seepage in an elastic medium, stress solutions have not been considered for the seepage force in a porous elasto‐plastic medium. This paper documents a study that investigates the stress behavior, caused by seepage, of a tunnel in an elasto‐plastic ground and its effects on the tunnel and ground. New elasto‐plastic solutions that adopt the Mohr–Coulomb failure criterion are proposed for a circular tunnel under radial flow conditions. A simple solution based on the hydraulic gradient obtained from a numerical parametric study is also proposed for practical use. It should be noted that the simple equation is useful for acquiring additional insight into a problem on a tunnel under drainage, because only a minimal computational effort is needed and considerable economic benefits can be gained by using it in the preliminary stage of tunnel design. The proposed equations were partly validated by numerical analysis, and their applicability is illustrated and discussed using an example problem. Comments on the tunnel analysis are also provided. Copyright © 2010 John Wiley & Sons, Ltd.     

软土地层浅埋圆形隧道非达西渗流场解析研究

软土在低水力坡降下的渗流会偏离达西定律,即为非达西渗流模式。假设孔隙水渗透服从指数渗流模式,采用镜像法原理推导了浅埋单孔和双孔圆形隧道非达西渗流场的解析解。结合算例,对浅埋圆形隧道非达西渗流解析解与达西渗流解析解进行了对比分析与验证,并对非达西渗流指数、隧道周围土体与衬砌渗流系数比值对隧道渗流场的影响进行了讨论。结果表明：非达西渗流指数、渗流系数比值对隧道渗流量和周围土体孔压均有较大的影响;随着渗流指数逐渐增大,土体内水头损失加快,隧道周围土体孔压及渗流量逐渐减小;随着土体与衬砌渗流系数比值逐渐增大,衬砌排水能力增强,隧道渗流量逐渐增大,隧道周围土体孔压减小更大。     

考虑不完全土拱效应的浅层地基竖向应力计算

浅层岩溶土洞塌陷和浅埋隧道施工等时常会引发浅层地基出现局部沉陷,导致地基可能承受不完全土拱效应作用。如何定量分析不完全土拱效应对浅层地基竖向应力的影响尤为重要。统计了国内外浅层活动门试验,将浅层地基滑移面概化为塔形,同时考虑了浅层地基不同深度处土层差异沉降及主应力偏转过程。通过建立主应力偏转角与活动门相对位移之间的数量关系,量化了浅层地基不同深度对应的不完全土拱效应发挥程度,优化了考虑不完全土拱效应的浅层地基竖向应力计算方法。分析了主要参数对不完全土拱效应的影响,结果表明,随着浅层活动门高宽比及相对位移的增大,应力转移量增加,土体有效内摩擦角及滑移面倾角则相反。可为局部沉降作用下的浅层地基竖向应力计算提供理论指导。     

Response of circular tunnel with imperfect interface bonding in layered ground subjected to obliquely incident plane P or SV wave

A semianalytical procedure is proposed for evaluating the internal forces of circular tunnel with imperfect interface bonding in layered ground subjected to an obliquely incident plane P or SV wave. In this study, the hoop bending moment and hoop axial force are related to the free‐field responses of the ground. A one‐dimensional numerical approach is firstly presented to obtain the free‐field responses of a layered half‐space with an obliquely incident plane P or SV wave propagation. Then, the free‐field stress state is transformed and decomposed in the polar coordinate system. The internal forces of tunnel caused by the isotropic stress state and the pure shear state are calculated and then summed up to obtain the overall analytical solutions. Finally, the validity of the proposed semianalytical procedure is demonstrated by numerical examples.     

半无限含水层中带衬砌隧洞渗流解析研究

现有关于半无限含水层中隧洞渗流的解析解不能考虑衬砌部分排水的特性。假定土体和衬砌均为饱和均匀连续介质,采用土体与衬砌分算;利用复变函数解决孔口问题基本方法,对土体采用共形映射坐标变换,将土体中的渗流方程转换为与衬砌部分类似的二维Laplace方程圆环域的Dirichlet问题;通过Fourier解法,根据边界条件并利用流量连续条件,求得半无限空间中带衬砌隧洞渗流量和衬砌周边水头解析解答。分析了径深比、土体与衬砌渗透系数相对值、内壁水压力等因素对渗流量和衬砌周边水头的影响,并与数值算例进行了比较。结果表明,对浅埋隧洞该解析解仍具有较好的精度,是简便、实用的计算方法。     

A finite strain closed‐form solution for the elastoplastic ground response curve in tunnelling

The ground response to tunnel excavation is usually described in terms of the characteristic line of the ground (also called ‘ground response curve’, GRC), which relates the support pressure to the displacement of the tunnel wall. Under heavily squeezing conditions, very large convergences may take place, sometimes exceeding 10–20% of the excavated tunnel radius, whereas most of the existing formulations for the GRC are based on the infinitesimal deformation theory. This paper presents an exact closed‐form analytical solution for the ground response around cylindrical and spherical openings unloaded from isotropic and uniform stress states, incorporating finite deformations and linearly elastic‐perfectly plastic rock behaviour obeying the Mohr–Coulomb failure criterion with a non‐associated flow rule. Additionally, the influence of out‐of‐plane stress in the case of cylindrical cavities under plane‐strain conditions is examined. The solution is presented in the form of dimensionless design charts covering the practically relevant parameter range. Finally, an application example is included with reference to a section of the Gotthard Base tunnel crossing heavily squeezing ground. The expressions derived can be used for preliminary convergence assessments and as valuable benchmarks for finite strain numerical analyses. Copyright © 2014 John Wiley & Sons, Ltd.     

Generalized Hoek-Brown solution of circular tunnel considering effects of axial stress and seepage forceEI北大核心CSCD

This paper studies the joint effect of seepage force and axial stress on the stress and displacement of circular tunnel. The circular tunnel is simplified as an axisymmetric model and the seepage field is simplified as volumetric force in the stress field. The excavation cross-section of surrounding rock is assumed as a plane as well, and an axial stress perpendicular to the plane is further introduced. Nonlinear solutions for the stress and displacement of circular tunnel are deduced considering the joint effect of axial stress and seepage force, based on the generalized Hoek-Brown failure criterion and the non-associated flow rule in elastic-brittle-plastic rock mass. Numerical simulations are also employed to analyze the distribution of stress field and displacement field in plastic zone of a circular tunnel under the joint effect of axial stress and seepage force. The calculated results show that the displacement in plastic zone increases significantly with the gradient increment of the seepage pressure, compared with the situation without seepage force. The radius and stress of surrounding rock in plastic zone increase when axial stress is the intermediate principal stress, while the radius and stress have less change when axial stress is the major or minor principal stress. It can be concluded that the seepage force has negative effects on the stability of circular tunnel, and the axial stress significantly influences the stress and displacement of the circular tunnel, especially in water-rich areas. Therefore, it is necessary to consider the joint effects of axial stress and seepage force to ensure the stability of circular tunnel in water-rich area.     

考虑轴向力和渗透力时圆形隧道广义Hoek-Brown解

为了研究渗透水压力和轴向应力共同作用时隧道围岩的应力和位移变化趋势,将圆形隧道简化为轴对称模型,假定渗透场以渗透体积力作用在原应力场,以围岩开挖断面为假定平面,引入垂直于该平面的轴向应力。基于广义Hoek-Brown强度准则和非关联流动法则,推导出考虑轴向应力和渗透场共同作用时弹-脆-塑性围岩的应力和位移非线性解,采用数值算例分析了轴向应力和渗透力共同作用时隧道围岩塑性区应力场和位移场的分布规律。计算结果表明：与无渗透水压力作用下的模型相比,渗透力作用使得围岩塑性区各点位移增大,并且内外水头差越大,位移增大越明显。当轴向应力为中主应力时,围岩塑性区半径和塑性区各点应力增大,轴向应力为大主应力和小主应力时,围岩塑性区半径和塑性区应力变化较小。因此,渗透力的存在不利于隧道围岩的稳定性,并且轴向应力的大小对于富水地区隧道围岩的应力和位移分布具有较大影响。在施工设计时考虑渗透力以及轴向应力的共同影响对于保证隧道围岩稳定性具有重要意义。     

Analytical and Numerical Study of the Mechanics of Rockbolt Reinforcement around Tunnels in Rock Masses

Summary  This paper addresses the problem of quantifying the mechanical contribution of rockbolts installed systematically around tunnels excavated in rock masses. The mechanical contribution referred to here is that of increased stress confinement and decreased tunnel convergences as compared with corresponding stresses and displacements obtained for non-reinforced tunnels. The problem is treated analytically first by presenting a closed-form solution for stress and displacement distributions around a circular tunnel excavated in elastic material and reinforced by grouted or anchored rockbolts. The analytical solution assumes that rockbolts are regularly spaced around the tunnel and that axi-symmetry conditions of geometry and loading apply. The results obtained with the closed-form solution are shown to be equivalent to the results of the same problem solved with traditional numerical methods. Based on the analytical and numerical results and by introducing dimensionless ratios that allow to quantify the increase of radial stresses and the decrease of radial displacements in the reinforced region of the tunnel, the paper shows that reinforcement can have a significant mechanical effect (i.e., increasing the confinement and decreasing the convergences) in tunnels excavated in rock masses of poor to very poor quality. The paper analyzes then the mechanical contribution of rockbolt reinforcement when the rock mass is assumed to behave elasto-plastically. For this case, it is shown that rockbolt reinforcement can also have a critical effect in controlling the extent of the plastic failure zone and the convergences of the tunnel. Correspondence: C. Carranza-Torres, Department of Civil Engineering, University of Minnesota, Duluth Campus, 1305 Ordean Court, Duluth, USA     

Radial Deformations Induced by Groundwater Flow on Deep Circular Tunnels

Summary. The magnitude and distribution of ground deformations around a tunnel are often monitored during construction and provide key information about ground-support interaction and ground behavior. Thus it is important to determine the effects of different parameters on ground deformations to accurately and effectively evaluate what contributes to ground and support behavior observed during excavation. This paper investigates one such relation: the effects of seepage on radial deformations. A number of numerical analyses have been conducted with the following assumptions: deep circular unsupported tunnel, elastic ground, isotropic far field stresses, dry ground or saturated ground with steady-state water seepage. The analyses cover a wide range of tunnel sizes, effective stresses, and pore pressures. Results from the numerical simulations confirm previous analytical solutions for normalized radial deformations behind the face (i.e. on the tunnel side of the face) of a tunnel excavated in dry ground, and have been used to propose a new analytical formulation for normalized radial displacements ahead and behind the tunnel face for both dry and saturated ground with water flow. Water seepage substantially increases the magnitude and distribution of the normalized radial deformations ahead of the face and at the tunnel face, but does not change much the displacement distribution behind the tunnel face.     

Analytical study of ground responses induced by the excavation of quasirectangular tunnels at shallow depths

The construction of quasirectangular tunnels at shallow depths is becoming increasingly common in urban areas to efficiently utilize underground space and reduce the need for backfilling. To clarify the mechanical mechanism of the stresses and displacements around the tunnels, this study proposes analytical solutions that precisely account for quasirectangular tunnel shapes, the ground surface, the tunnel depth, and the ground's elastic/viscoelastic properties. The Schwarz alternating method combined with complex variable theory is employed to derive the elastic solution, and convergent and highly accurate solutions are obtained by superposing the solutions in the alternating iterations. Based on the solution and the extended corresponding principle for the viscoelastic problem, the time-dependent analytical solutions for the displacement are obtained for the ground assuming any viscoelastic model. The analytical solutions agree well with the finite element method (FEM) numerical results for models that are completely consistent, and qualitatively agree with field data. Furthermore, based on the stress solution combined with the Mohr-Coulomb failure criterion, the predicted initial plastic zone and propagation directions around the tunnels are qualitatively consistent with those determined by the limit analysis. A parametric study is performed to investigate the influences of the rectangular/quasirectangular tunnel shape, burial depth, and supporting pressure on the ground stresses and displacements.     

Plasticity solutions for soil behaviour around contracting cavities and tunnels

The action of tunnel excavation reduces the in-situ stresses along the excavated circumference and can therefore be simulated by unloading of cavities from the in-situ stress state. Increasing evidence suggests that soil behavior in the plane perpendicular to the tunnel axis can be modelled reasonably by a contracting cylindrical cavity, while movements ahead of an advancing tunnel heading can be better predicted by spherical cavity contraction theory. In the past, solutions for unloading of cavities from in-situ stresses in cohesive-frictional soils have mainly concentrated on the small strain, cylindrical cavity model. Large strain spherical cavity contraction solutions with a non-associated Mohr–Coulomb model do not seem to be widely available for tunnel applications. Also, cavity unloading solutions in undrained clays have been developed only in terms of total stresses with a linear elastic-perfectly plastic soil model. The total stress analyses do not account for the effects of strain hardening/softening, variable soil stiffness, and soil stress history (OCR). The effect of these simplifying assumptions on the predicted soil behavior around tunnels is not known. In this paper, analytical and semi-analytical solutions are presented for unloading of both cylindrical and spherical cavities from in-situ state of stresses under both drained and undrained conditions. The non-associated Mohr-Coulomb model and various critical state theories are used respectively to describe the drained and undrained stress-strain behaviors of the soils. The analytical solutions presented in this paper are developed in terms of large strain formulations. These solutions can be used to serve two main purposes: (1) to provide models for predicting soil behavior around tunnels; (2) to provide valuable benchmark solutions for verifying various numerical methods involving both Mohr–Coulomb and critical state plasticity models. Copyright © 1999 John Wiley & Sons, Ltd.     

Stresses around Pressure Tunnels with Semi-permeable Liners

Summary. The liner of a pressure tunnel needs to be designed such that it can withstand the loads from the ground, the internal pressure, and minimize the development of significant pore pressures at the liner-ground interface. Pore pressures behind the liner reduce the effective stresses in the ground immediately in contact with the liner and can ultimately produce loss of support from the ground. Deformations and loads of the liner are intimately connected to the interplay that exists between liner, ground, and pore pressures in the ground. A closed-form analytical solution has been derived that accounts for the inter-relation between liner, ground, and pore pressures. Elastic response of the liner and ground, and plane strain conditions at any cross-section of the tunnel are assumed. The solution shows that stresses in the ground depend on the following dimensionless factors: relative stiffness of the ground and liner, ground Poisson’s ratio, surface slope angle, coefficient of earth pressure at rest, relative tunnel depth, and magnitude of the pore pressure behind the liner relative to the internal pressure. The minimum ground effective tangential stresses at the ground-liner interface increase with the relative stiffness of the liner, with the coefficient of earth pressure at rest, and with tunnel depth. They decrease with increasing surface slope angle and pore pressures behind the liner. As leakage through the liner increases, the pore pressures in the ground increase. This results in a decrease of effective radial and tangential stresses in the ground while displacements and loads of the liner are relatively less affected.     

马蹄形隧洞考虑支护滞后过程的应力分析

基于平面弹性复变函数中的保角变换方法,推导出带有衬砌的非圆形隧洞在原始地应力作用下的应力解析解。根据衬砌内边界的应力边界条件及围岩衬砌接触面上的应力和位移连续条件,获得求解围岩和衬砌解析函数的基本方程,计算了围岩和衬砌中的应力和位移。在求解过程中,考虑了支护滞后于开挖的力学过程,并认为围岩和衬砌之间紧密接触,不会相互分开和相对滑动。以马蹄形隧洞为例,获得了围岩开挖边界和衬砌内外边界的切向应力及围岩与衬砌接触面上的接触应力分布规律,并与ANSYS数值方法结果对比,算例表明两种方法的计算结果吻合很好。     

Analytical‐numerical solution for stress distribution around tunnel reinforced by radial fully grouted rockbolts

This paper presents an analytical‐numerical approach to obtain the distribution of stresses and deformations around a reinforced tunnel. The increase in the radial stress of the reinforced tunnel, based on the performance of a bolt, is modeled by a function, which its maximum value is in the vicinity of the bolt periphery and it exponentially decreases in the far distance from the bolt. On the basis of this approach, the shear stiffness between the bolt and the rock mass and the shear stress distribution around the bolt within the rock mass are also analytically obtained. The results are compared with those obtained by the assumption of ‘uniform increase of radial stress’ method, which is made by the previous studies. The analyses show when the bolts' spacing is large, the safety factor must be increased if the ‘uniform increase of radial stress’ method is used for the design. Finally, a procedure is introduced to calculate the non‐equal deformation of the rock mass between the bolts at any radius that can be useful to compute the bending moment in shotcrete layer in New Austrian Tunnelling Method (NATM) approach. Copyright © 2016 John Wiley & Sons, Ltd.