Stability of long circular tunnels in sloping ground

Assessment of tunnel stability has become increasingly crucial as more and more tunnels are built in difficult terrains such as sloping ground. The required support pressure on the tunnel walls associates both tunnel stability and liner design considerations. The present analysis attempts to find a uniform internal pressure which can support a circular tunnel built in a sloping ground with a particular level of stability in cohesive-frictional soils. The lower bound finite element limit analysis has been applied to find the required minimum uniform internal support pressure presented as a non-dimensional term p/c; where p is the minimum normal internal pressure on the tunnel boundary to avoid collapse and c is the cohesion of soil. The variation of p/c is presented for a range of normalised embedment depth of tunnel (H/D), stability number (γD/c), internal friction angle of soil (?) and slope angle (β); where H is the crown depth of the tunnel, D is the tunnel diameter and γ is the unit weight of soil. Appropriate comparisons have been carried out with available literature. Failure patterns of the tunnel have also been studied to understand the extent and the type of failure zone which may generate during the collapse.     

Upper-bound multi-rigid-block solutions for bearing capacity of two-layered soils

The limit analysis method has been widely used in the stability analysis of geotechnical problems including the bearing capacity of foundations. Two main approaches have been followed in the limit analysis to improve the calculation of the bearing capacity of foundations. One approach is to combine limit analysis with the finite element method and linear/nonlinear programming. The other is to use a multi-rigid-block mechanism to obtain an upper-bound solution. In this paper, the multi-rigid-block upper-bound method with a modified failure mechanism from that of Florkiewicz [Florkiewicz A. Upper bound to bearing capacity of layered soils. Can Geotech J 1989;26(4):730–6.] was employed to calculate the bearing capacity of foundations. Attention was paid particularly to the bearing capacity of strip footings over a two-layered soil. In order to verify the effectiveness of the modified mechanism, comparisons were made with other well-known solutions. The results showed improvements over the best available multi-rigid-block upper-bound solutions given by Michalowski and Shi [Michalowski RL, Shi L. Bearing capacity of footings over two-layer foundation soils. J Geotech Eng ASCE 1995;121(5):421–8.], and fair consistence with the results from the finite element limit analysis in Shiau et al. [Shiau JS, Lyamin AV, Sloan SW. Bearing capacity of a sand layer on clay by finite element limit analysis. Can Geotech J 2003;40(5):900–15.].     

Undrained stability of wide rectangular tunnels

Upper bound rigid block methods and finite element limit analysis are applied to investigate the undrained stability of rectangular tunnels. Solutions for a range of soil properties and tunnel geometries are obtained for situations in which the tunnel is wider than it is high. The upper and lower solutions for tunnel stability, which accurately bracket the true collapse load, are summarised in a series of stability charts for use in tunnel design.     

地表超载作用下隧道失稳破坏的上限有限元分析

针对地表超载作用下隧道稳定性和破坏模式问题,基于刚体平动运动单元上限有限元理论编程并计算分析,获得了浅埋隧道失稳临界超载系数上限解和刚性运动块体体系破坏模式。通过与现有的刚性块体极限分析上限法以及极限分析上、下限有限元法计算结果的对比分析,验证了上限解的可靠性。研究结果表明,（1）临界超载系数 黏聚力c之比 随土体内摩擦角 和隧道埋深C与直径D之比( )的增大而相应增大,随土体重度与黏聚力参数 的增大而减小;（2） 和 对隧道破坏模式的影响较明显; 增大,则隧道破坏范围增加;内摩擦角 增大,刚性运动块体破坏模式相互错动更加显著,相比而言, 对破坏模式的影响并不显著;（3）刚体平动运动单元上限有限元上限解精度高,所得刚性运动块体破坏模式具有滑移线形态,能精细地反映隧道失稳破坏特征。     

软基上沉入式大圆筒结构的水平承载力分析

大圆筒结构是一种适宜于软土地基上的新型港口水工建筑物,但目前对于软基上大圆筒结构的水平极限承载力尚缺乏合理的计算模式和分析方法.结合长江口深水航道治理二期工程,针对大圆筒结构设计方案提出了一种简便而实用的上限极限分析方法,采用有限元数值分析解验证了所提出计算方法的合理性,进而,据此进行了一系列变动参数比较计算,以此对大量的计算结果进行了系统的分析,所得到的数值计算结果与分析结论为工程设计提供参考依据.     

Analysis and Design of Retaining Wall having Reinforced Cohesive Frictional Backfill

The case of a rigid wall with inclined back face retaining reinforced cohesive-frictional backfill subjected to uniformly distributed surcharge load has been analyzed using limit equilibrium approach. The analysis considers the stability of an element of the failure wedge, which is assumed to develop in the reinforced earth mass adjoining the back face of wall. The non-dimensional charts have been developed for computing the lateral earth pressure on wall and the height of its point of application above the base of wall. The theoretical findings have been verified by model tests on a rigid wall retaining a dry cohesive-frictional soil reinforced by geogrid strips. Experimental results are in good agreement with the theoretical predictions. A design example has been included to illustrate the design procedure.     

基于极限超载的无衬砌浅埋土质隧道稳定性研究

为了从理论上判别和度量无衬砌浅埋土质隧道的稳定性,考虑拱形隧道洞身范围斜向滑裂面作用,采用极限平衡法导出隧道发生整体失稳时的极限超载与稳定安全系数理论解,并通过既有典型加载试验进行验证。浅埋隧道极限超载理论解与典型加载试验结果较接近,在此基础上,将抗滑力与滑动力之比定义为整体稳定安全系数,并分析其影响因素。研究表明,稳定安全系数随黏聚力、内摩擦角的增大而增大,黏聚力对其影响显著,而内摩擦角对其影响不显著;稳定安全系数随隧道跨度、覆跨比、超载的增大而减小,但超载越大对稳定安全系数影响越小;随着高跨比的增大,若黏聚力较低则安全系数减小,若黏聚力较高则安全系数先增大后减小。通过典型的加载试验对整体稳定安全系数进行了校验,表明稳定安全系数计算值基本满足工程允许误差要求。     

Probabilistic stability assessment using adaptive limit analysis and random fields

For deterministic scenarios, adaptive finite element limit analysis has been successfully employed to achieve tight bounds on the ultimate load of a geotechnical structure in a much more efficient manner than a dense uniform mesh. However, no probabilistic studies have so far considered finite element limit analysis with adaptive remeshing. Therefore, this research explores the benefits of combining adaptive mesh refinement with finite element limit analysis for probabilistic applications. The outcomes indicate that in order to achieve tight bounds on probabilistic results (such as the probability of failure), the ultimate load in each individual simulation (e.g. factor of safety or bearing capacity) has to be estimated with a very high level of accuracy and this can be achieved more economically using adaptive mesh refinement. The benefits, assessed here for undrained conditions, are expected to be much more pronounced in the case of frictional soils and complex geometries.     

基于方位离散线性化的上限原理有限元法

上限原理有限元法不仅可以得到边坡的安全系数,还可以给出临界滑动面,且具有比极限平衡法更严谨的理论基础,因此,拥有更广阔的应用前景。针对传统的上限有限元法不能考虑强度各向异性的问题,提出了一种新的摩尔-库仑屈服面线性化方法。该方法在对方位角离散化的基础上,建立了线性化的方位离散塑性流动约束方程,丰富了基于线性规划的上限法理论。两个算例结果表明：该方法可以稳定地从极限解的上方收敛;且对边坡进行稳定性分析,若忽略了边坡的强度各向异性,则会高估边坡的稳定性,得到较大的安全系数。     

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.     

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.     

地铁荷载作用下软土地基长期沉降的有限元分析

为讨论地铁荷载作用下上海软土地基的长期沉降问题,本文利用边界面本构模型结合弹塑性动力固结有限元理论,通过ABAQUS数值软件分别建立了轨道-隧道-土体相互作用模型,考察单孔和平行双孔两种隧道布置情形。分析了单次、多次地铁列车动荷载加载引起的隧道地基沉降规律,并结合3种经验模型预测了隧道的长期沉降,得到3种不同加载工况下由列车动荷载引起的动应力幅值与沉降;并对比分析了单双孔隧道的沉降特性,将模型计算结果与沉降实测值进行比较。结果表明,采用模型的预测值与沉降实测值接近;双孔隧道的土体沉降不仅是单孔隧道沉降值的简单线性叠加,而且存在一定程度的增大。     

浅埋隧道掌子面稳定性二维自适应上限有限元分析

为研究浅埋隧道掌子面稳定性及获取精细化的破坏模式,提出了一种上限有限元非结构化网格自适应加密策略。以单元耗散能权重指标作为网格自适应加密评判准则,该策略同时兼顾了单元尺度与塑性应变。应用高阶的6节点三角形单元并建立上限有限元线性规划模型,以多次反复计算和网格加密的方式实现了二维自适应上限有限元分析并编制了计算程序。利用条形基础地基极限承载力课题,从上限解精度和网格加密形态方面验证了该程序的有效性。针对浅埋隧道掌子面稳定性问题,展开多参数条件下的自适应上限有限元计算,分析了网格加密过程中单元总数与上限解精度的关系,列出不同隧道埋深和内摩擦角对应的隧道掌子面稳定性临界值的上限解,揭示出掌子面稳定性变化规律及精细化的破坏模式。     

考虑桩土侧移的被动桩中土拱效应数值分析

被动桩对侧向位移的土层起到遮拦作用的机制主要是土拱效应。采用土工有限元软件Plaxis Tunnel 3D 1.2,对堆载荷载作用下邻近桩基中的土拱效应产生机制和性状进行三维数值分析,指出目前被动桩中土拱效应二维有限元分析存在的问题。考虑桩土侧移与相对位移,再利用土工有限元软件Plaxis2D 8.2详细地研究了侧向土体位移大小、桩身水平位移大小、土体性质以及桩土接触面性质等影响因素对土拱效应性态和桩土荷载分担比的影响。     

深埋盾构隧道开挖面三维对数螺旋破坏模式的上限分析

相对于盾构隧道施工的大量需求与快速发展的状况,国内在盾构工法特别是大型深埋盾构隧道施工技术和理论研究方面还存在不足,特别是水压条件下深埋盾构隧道开挖面稳定问题。基于极限分析上限法和水土压力统一参数,对考虑水压影响的均质土深埋隧道开挖面稳定性计算方法进行研究,建立了考虑水压影响的深埋盾构隧道开挖面三维对数螺旋破坏模式模型,并推导了其极限支护压力计算公式。然后利用土层厚度加权平均法,可将上述方法应用于多层土深埋盾构隧道开挖面稳定性的评价中。最后,以上海长江盾构隧道实际工程为例,采用本文推导的极限分析上限三维对数螺旋破坏模式方法计算并分析其极限支护压力,并将计算结果与前人研究和规范方法计算的结果进行对比分析。通过该研究可改进与完善水压条件下深埋盾构隧道极限支护压力确定方法,从而为考虑水压条件下盾构隧道施工支护压力的合理确定提供理论依据。     

Analytical solutions of linear finite‐ and small‐strain one‐dimensional consolidation

Analytical solutions are presented for linear finite‐strain one‐dimensional consolidation of initially unconsolidated soil layers with surcharge loading for both one‐ and two‐way drainage. These solutions complement earlier solutions for initially unconsolidated soil layers without surcharge and initially normally consolidated soil layers with surcharge. Small‐strain solutions for the consolidation of initially unconsolidated soil layers with surcharge loading are also presented, and the relationship between the earlier solutions for initially unconsolidated soil without surcharge and the corresponding small‐strain solutions, which was not addressed in the earlier work, is clarified. The new solutions for initially unconsolidated soil with surcharge loading can be applied to the analysis of low stress consolidation tests and to the partial validation of numerical solutions of non‐linear finite‐strain consolidation. They also clarify a formerly perplexing aspect of finite‐strain solution charts first noted in numerical solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Appropriate Uses and Practical Limitations of 2D Numerical Analysis of Tunnels and Tunnel Support Response

In spite of the gradual development of three-dimensional analysis packages utilizing finite element models or finite difference algorithms for stress–strain calculations, two-dimensional (2D) analysis is still used as the primary engineering tool for practical analysis of tunnel behavior and tunnel support performance for design—particularly at the preliminary stage of a project. The applicability of 2D finite element analysis or analytical convergence confinement solutions to staged support installation depend on the application of an assumed or validated longitudinal displacement profile. Convergence in commonly applied 2D staged models is controlled by boundary displacement or internal pressure relaxation. While there have been developments to improve this methodology, this often assumes independence between the ground reaction curve and the support resistance, independence between the longitudinal displacement profile to support application, and the assumption that non-isotropic stresses and non-circular geometries can be handled in the same way as circular tunnels in isotropic conditions. This paper examines the validity of these assumptions and the error inherent in these extensions to 2D tunnel analysis. Anisotropic stresses and lagged (staged) excavation present a particular problem. Practical solutions are proposed for support longitudinal displacement LDPs in simplified conditions.     

极限平衡条件的拓展与土工结构稳定分析

与土体稳定有关的工程问题诸如挡土结构物上的土压力、地基承载力以及土坡稳定性都与土的极限平衡状态相关 联。将土体一点（单元）的极限平衡条件推广,给出土体沿滑动面的极限平衡条件,证明了该条件的充分必要性。依据该条件,提出了有限元应力分析与极限平衡结合的有限元极限平衡法及土体稳定的安全系数,并明确了安全系数的物理意义,以评价土工结构的稳定性。随后,列举了该方法在土工结构稳定分析中的一些具有代表性的应用,证明了该法是有效的,可以用于分析岩土工程中的稳定问题。     

Analytical analysis of working face passive stability during shield tunneling in frictional soils

This paper develops the 2D and 3D kinematically admissible mechanisms for analyzing the passive face stability during shield tunneling using upper-bound analysis. The mechanisms consider trapezoidal distribution of support pressure along tunnel face and partial failure originated at tunnel face above invert. For cohesionless soils, the support pressure is a function of soil effective frictional angle φ′ which determines the inclination of failure block and the normalized soil cover depth C/D (soil cover depth/tunnel diameter) which affects the origination of the passive failure. For cohesive soils, the support pressure is a function of φ′, C/D, and the effective cohesion c′. The cohesion c′ has a relatively smaller impact on the support pressure than φ′ and C/D have. The mechanisms are verified by comparing the current solutions with a previous upper-bound solution. The comparison shows that the current solutions are a general solution which is capable of predicting the passive face failure originated at any depth along tunnel face and the previous solution is a particular solution with the assumption that the face failure originated at tunnel invert. The mechanisms are validated through application to a practical project of shallowly buried, large diameter underwater tunnel. The validation shows that the mechanisms are capable of assessing the tunnel face passive instability rationally.

     

正常固结黏土中平板锚基础的吸力和抗拉力

平板锚是新近出现的一种系泊深海浮式结构的基础型式。当黏土地基中的平板锚承受上拔力时,平板上、下表面超静孔压差形成的吸力使其抗拉承载力显著增加,对于风浪等快速加载条件尤其如此。利用有限元软件ADINA建立有效应力形式的轴对称动力有限元模型,研究圆形平板锚在缓慢加载与快速加载时的超静孔压分布与地基破坏型式。快速加载时重黏土和高岭土两组典型正常固结土样所得极限承载力系数与塑性极限分析解一致。进而通过变动参数分析,讨论加载速率和埋深对吸力和总抗拉力极限值的影响,并给出排水和不排水加载条件对应的临界加载速率。结果表明,不排水加载条件的总抗拉力可能达到排水总抗拉力的3倍。