Elastic analysis of ground movements around a tunnel considering a buoyant lining moving upwards

A complex variable solution is derived for the ground movements around a tunnel considering the lining at different positions in the drilled hole in a linear elastic half-plane. The soil displacements for three different convergence patterns are presented and compared. It is shown that the surface settlement trough will become flatter and wider when the lining moves upwards; the horizontal displacements are more concentrated in the region below the lining; and the region influenced by displacements is wider. Most calculated settlement troughs can be well represented by Gaussian curves with simple parameters, but the obtained width parameters are larger than the traditional empirical values. The rigorous solution can be used as a simple tool to understand the fundamental mechanism of the lining buoyancy problem and as a benchmark for numerical procedures based on more sophisticated models. The importance of buoyancy effect on the displacement field should be taken into account.     

Complex variable analysis for stress distribution of an underwater tunnel in an elastic half plane

When an underwater tunnel is excavated, the groundwater may flow into the tunnel. The seepage forces consequently induced can have important effects on the effective stresses around the tunnel. Moreover, the influences of the free surface of a shallow underwater tunnel should also be considered. In this research, an analytical solution is presented to calculate the seepage‐induced effective stresses around a shallow underwater tunnel in an elastic half plane. The solution uses the complex variable method and consists of conformally mapping the half plane with a hole onto a transformed circular ring. The coefficients of the various terms in the Laurent series expansions of the stress functions in the transformed region can be obtained from the boundary conditions. The total stress distribution around a shallow underwater tunnel can be calculated by the potentials in the half plane. The effective stress can be obtained by subtracting the pore pressure from the total stress. The analytical solution is validated by numerical simulations and can be used to perform both the short‐term and long‐term analyses. By using the proposed solution, it is found that the circumferential effective stresses around the tunnel increase greatly because of seepage, and they increase with the increase of water depth in both the undrained and drained conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Deformation of a poroelastic layer overlying an elastic half‐space due to dip‐slip faulting

An analytical solution of the plane strain problem of the deformation of a homogeneous, isotropic, poroelastic layer of uniform thickness overlying a homogeneous, isotropic, elastic half‐space due to two‐dimensional seismic sources buried in the elastic half‐space has been obtained. The integral expressions for the displacements, stresses and pore pressure have been obtained using the stress function approach by applying suitable boundary conditions at the free surface and the interface. The solution obtained is in the Laplace–Fourier transform domain. The case of a vertical dip‐slip line dislocation for the oceanic crust model of Earth is studied in detail. Schapery's formula is used for the Laplace inversion and the extended Simpson's formula for the Fourier inversion. Diffusion of pore pressure in the layer is studied numerically. Contour maps showing the pore pressure in the poroelastic layer have been plotted. The effect of the compressibility of the solid and fluid constituents on pore pressure has also been studied. Copyright © 2015 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.     

Forced vertical and horizontal movements of a rectangular rigid foundation on a transversely isotropic half‐space

An analytical investigation of a half‐space containing transversely isotropic material under forced vertical and horizontal displacements applied on a rectangular rigid foundation is presented in this paper. With the goal of a rigorous solution to the shape‐ and rigidity‐ induced singular mixed boundary value problem, the formulation employs scalar potential representation, the Fourier expansion and the Hankel integral transforms method to obtain the surface arbitrary point‐load solution in cylindrical coordinate system. The obtained Green's functions are rewritten in rectangular coordinate system, allowing the response of the half‐space because of an arbitrary distributed load on a rectangular surface area be given in terms of a double integral. The numerical evaluations of stresses are done with the use of an element, which is singular at the edge and the corner of the rectangle. Upon the imposition of the rigidity displacement boundary condition for a rigid foundation and the use of a set of two‐dimensional adaptive‐gradient elements, which can capture the singular behavior in the contact stress effectively, a set of new numerical results are presented to illustrate the effect of transverse isotropy on the foundation response. Copyright © 2012 John Wiley & Sons, Ltd.     

An analytical solution for a jointed shield‐driven tunnel lining

This paper presents an analytical solution for the prediction of internal forces and displacements of a jointed segmental precast circular tunnel lining. The effects of joint stiffness on the performance of the tunnel lining are discussed. The ‘force method’ is used to determine the internal forces and displacements of jointed tunnel lining. Five shield‐driven tunnel cases are adopted to study the effects of joint stiffness, soil resistance, joint distribution and joint number on the internal forces and displacements of circular tunnels. Laboratory model tests are conducted to verify the proposed analytical solution. Copyright © 2001 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.     

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.     

Response of a tunnel deeply embedded in a viscoelastic medium

This paper presents a three‐dimensional energy‐based solution for the time‐dependent response of a deeply embedded and unsupported semi‐infinite tunnel of circular cross‐section. The tunnel is taken to be excavated quasi‐instantaneously from an infinite rock body that initially exhibits an isotropic stress state and that is made up of a homogeneous, isotropic and viscoelastic material. The viscoelastic behaviour is modelled by means of Burger's model, and the rock is taken to behave volumetrically linear elastic and to exhibit exclusively deviatoric creep. This viscoelastic problem is transformed into the Laplace domain, where it represents a quasi‐elastic problem. The displacement fields in the new solution are taken to be the products of independent functions that vary in the radial and longitudinal directions. The differential equations governing the displacements of the system and appropriate boundary conditions are obtained using the principle of minimum potential energy. The solutions for these governing equations in the Laplace domain are then obtained analytically and numerically using a one‐dimensional finite difference technique. The results are then transformed back into the time domain using an efficient numerical scheme. The accuracy of the new solution is comparable with that of a finite element analysis but requires much less computation effort. Copyright © 2011 John Wiley & Sons, Ltd.     

深埋圆形隧道弹塑性位移统一解及其比较分析

基于统一强度理论和非关联线性流动法则,考虑塑性区的真实弹性应变,推导了深埋圆形隧道弹塑性位移统一解。位移统一解合理地考虑了中间主应力效应和围岩剪胀特性的影响,具有广泛的理论意义。经与文献简化位移解析解比较分析知,围岩塑性区的弹性应变对隧道位移影响显著,简化解析解所得结果偏小,简化解析解和文中统一解之间的差异与中间主应力效应和围岩剪胀特性密切相关;考虑中间主应力的影响,能更加充分发挥围岩的强度潜能;隧道支护设计应考虑围岩剪胀特性的重要影响     

Bayesian Estimation of Boundary Conditions with Application to Deep Tunneling

An iterative procedure is proposed to determine the far-field state of stress that exists in a rock mass with non-linear behavior around a tunnel. Absolute displacements of the tunnel wall obtained by means of geodetic measurements are used to identify the boundary stress conditions. Previous information is accounted for, so that boundary conditions can be updated at the various stages of a project, as soon as new information becomes available. The application to the 2-D plane strain synthetic model of a tunnel in a yielding rock mass shows the fast convergence of the procedure and the need to use absolute displacements (rather than relative displacements) in order for the identification problem to be well-posed.     

半空间饱和土中圆形衬砌对弹性压缩波的散射

借助于Biot 波动理论和弹性波的传播理论,采用复变函数和多级坐标法,对半空间饱和土中圆形衬砌结构对弹性稳态压缩波的散射问题进行求解和分析。利用一个半径很大的圆弧来逼近半空间直边界,将待解问题转化为稳态弹性压缩波在一个大圆孔和一个弹性衬砌结构的散射问题。通过引入势函数,将饱和土的Biot波动方程和衬砌的弹性波动方程解耦成Helmholtz 方程,借助复变函数级数展开便可以预先写出该组Helmholtz方程的通解。然后,通过引用复变量,把饱和土和衬砌结构中的应力、位移及孔压用设定的势函数表示出来,再利用半空间饱和土和衬砌结构的连续性条件和近似直边界的圆弧边界和衬砌内边界的边界条件求解出该组势函数的特解。最后,利用势函数的特解,得到饱和土中的位移,应力和孔压及衬砌结构的位移和应力;变换不同的参数求解衬砌结构内外边界的动应力和孔压的集中系数,通过对算例结果的分析得出一系列有益的结论。     

An analytical solution for axisymmetric tunnel problems in elasto-viscoplastic media

An analytical solution is presented for the time-dependent stresses and displacements in plane strain around a circular hole when it is loaded by an axisymmetric internal and far-field pressure. The material is assumed to be elasto-viscoplastic with dilatant plastic deformations according to a non-associated flow rule. Strain softening is considered by a modified St. Venant slider which is characterized by Mohr-Coulomb yield conditions for both the peak and the residual strengths.     

Surface movements in a rock massif induced by drainage associated to tunnel excavation

Surface movements were measured in the Gotthard massif as the Gotthard Base Tunnel was excavated. These movements might damage concrete dams constructed on the surface valleys. The leading assumption of this work is that deformation is induced by the dissipation of pore pressures in the massif caused by tunnel drainage. Deformations induce both horizontal and vertical surface displacements. Horizontal displacements, may lead to valley closures if they are in opposite direction, which would induce negative effects on arch dams. An analytical solution is derived using the method of images and an approximated integration of deformations to calculate the movements and the flow rate collected in the tunnel. Numerical calculations were carried out in 2D (vertical cross section) and 3D to investigate the problem under different conditions and to study the effect of parameters. The 3D models permit to incorporate the presence of a vertical fracture perpendicular to the tunnel that increases the drainage and pressure drop as it is hit by the tunnel. It was also possible to simulate the impermeabilization works in the tunnel to reduce drainage and consequently, movements. Copyright © 2012 John Wiley & Sons, Ltd.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Source signature and elastic waves in a half‐space under a momentary shear line impulse

The transient deformation of an elastic half‐space under a line‐concentrated impulsive vector shear load applied momentarily is disclosed in this paper. While in an earlier work, the author gave an analytical–numerical method for the solution to this transient boundary‐value problem, here, the resultant response of the half‐space is presented and interpreted. In particular, a probe is set up for the kinematics of the source signature and wave fronts, both explicitly revealed in the strained half‐space by the solution method. The source signature is the imprint of the spatiotemporal configuration of the excitation source in the resultant response. Fourteen wave fronts exist behind the precursor shear wave S: four concentric cylindrical, eight plane, and two relativistic cylindrical initiated at propagating centres that are located on the stationary boundaries of the solution domain. A snapshot of the stressed half‐space reveals that none of the 14 wave fronts fully extend laterally. Instead, each is enclosed within point bounds. These wave arresting points and the two propagating centres of the relativistic waves constitute the source signature. The obtained 14 wave fronts are further combined into 11 disparate wave fronts that are grouped into four categories: an axis of symmetry wave—so named here by reason of being a wave front that is contiguous to the axis of symmetry, three body waves, five surface waves and two inhibitor waves—so named here by reason that beyond them the material motion dies out. Of the three body waves, the first is an unloading shear wave, the second is a diffracted wave and the third is a reflected longitudinal two‐branch wave. Of the two inhibitor waves, the first is a two‐joint relativistic wave, while the second is a two‐branch wave. The wave system, however, is not the same for all the dependent variables; a wave front that appears in the behaviour of one dependent variable may not exist in the behaviour of another. It is evident from this work that Saint–Venant's principle for wave propagation problems cannot be formulated. Therefore, the above results are valid for the particular proposed model for the momentary line‐concentrated shear load. The formulation of the source signature, the wave system, and their role in the half‐space transient deformation are presented here. Copyright © 2003 John Wiley & Sons, Ltd.     

Torsional dynamic response of a large‐diameter pipe pile in viscoelastic saturated soil

An analytical solution is developed in this paper to investigate the dynamic response of a large‐diameter end‐bearing pipe pile subjected to torsional loading in viscoelastic saturated soil. The wave propagation in saturated soil and pile are simulated by Biot's two‐phased linear theory and one‐dimensional elastic theory, respectively. The dynamic equilibrium equations of the outer soil, inner soil, and pile are established. The solutions for the outer and inner soils in frequency domain are obtained by Laplace transform technique and the separation of variables method. Then, the dynamic response of the pile is obtained on the basis of the perfect contacts between the pile and the outer soil as well as the inner soil. The results in this paper are compared with that of a solid pile in elastic saturated soil to verify the validity of the solution. Furthermore, the solution in this paper is compared with the classic plane strain solution to verify the solution further and check the accuracy of the plane strain solution. Numerical results are presented to analyze the vibration characteristics and illustrate the effect of the soil parameters and the geometry size of the pile on the complex impedance and velocity admittance of the pile head. Finally, the displacement of the soil at different depth and frequency is analyzed. Copyright © 2014 John Wiley & Sons, Ltd.