砂土地层中不同支护模式下盾构隧道开挖面稳定性研究

防止盾构隧道开挖面失稳的关键是合理设置不同盾构支护平衡模式下的支护压应力。在改进的筒仓楔形体模型计算方法得出的开挖面松动土体对刀盘压力呈近似呈抛物线分布的基础上,研究了气压支护模式、泥水支护模式和土压支护模式下,盾构隧道开挖面分别在地下水位以上和地下水位以下时开挖面的稳定性,研究结果表明:有效支护应力均匀分布时,除粘土开挖面下部失稳外,其余土体均为开挖面中下部失稳;有效支护压应力呈上小下大的梯形分布时,除软粘土开挖面下部失稳外,其余土体均为开挖面上部失稳;有效支护应力呈上大下小的梯形分布时,所有土体开挖面均为下部失稳;在气压、泥水和土压平衡支护模式下,开挖面在未到达筒仓楔形体模型所假设的开挖面整体失稳前,开挖面已经发生了局部失稳,采用筒仓楔形体模型确定的极限稳定支护力是不安全的。最后给出了开挖面松动土体对刀盘压应力公式中计算参数的无量纲化图,以方便实际工程运用。      

A dynamic two‐phase flow model for air sparging

Air sparging (AS) is an in situ soil/groundwater remediation technology, which involves the injection of pressurized air/oxygen through an air sparging well below the zone of contamination. Characterizing the mechanisms governing movement of air through saturated porous media is critical for the design of an effective cleanup treatment system. In this research, micromechanical investigation was performed to understand the physics of air migration and subsequent spatial distribution of air at pore scale during air sparging. The void space in the porous medium was first characterized by pore network consisting of connected pore bodies and bonds. The biconical abscissa asymmetric concentric bond was used to describe the connection between two adjacent pore bodies. Then a rule‐based dynamic two‐phase flow model was developed and applied to the pore network model. A forward integration of time was performed using the Euler scheme. For each time step, the effective viscosity of the fluid was calculated based on fractions of two phases in each bond, and capillary pressures across the menisci was considered to compute the pressure field. The developed dynamic model was used to study the rate‐dependent drainage during air sparging. The effect of the capillary number and geometrical properties of the network on the dynamic flow properties of two‐phase flow including residual saturation, spatial distribution of air and water, dynamic phase transitions, and relative permeability‐capillary pressure curves were systematically investigated. Results showed that all the above information for describing the air water two‐phase flow are not intrinsic properties of the porous medium but are affected by the two‐phase flow dynamics and spatial distribution of each phase, providing new insight to air sparging. Copyright © 2012 John Wiley & Sons, Ltd.     

A numerical study of the effect of soil and grout material properties and cover depth in shield tunnelling

For shield-driven tunnels, the influence of the soil and grout material properties and of the cover depth on the surface settlements, the loading and deformation of the tunnel lining and the steering of the TBM is investigated numerically. To this end, comparative numerical simulations of a mechanised tunnel advance in homogeneous, overconsolidated, soft, cohesive soil below the ground water table are performed and sensitivities are evaluated. The advancement of the step-by-step tunnel construction process is modelled using a three-dimensional finite element model, which takes into account all relevant components of shield tunnelling. The material behaviour of the saturated soil and the tail void grout is modelled by a two-field finite element formulation in conjunction with an elasto-plastic Cam-Clay model for the soil and a hydration-dependent constitutive model for the grout. The analyses provide valuable information with regard to the significance of the investigated parameters and demonstrate the complexity of the various interactions in shield tunnelling.     

Weak discontinuity in porous media: an enriched EFG method for fully coupled layered porous media

In this paper, a series of multimaterial benchmark problems in saturated and partially saturated two‐phase and three‐phase deforming porous media are addressed. To solve the process of fluid flow in partially saturated porous media, a fully coupled three‐phase formulation is developed on the basis of available experimental relations for updating saturation and permeabilities during the analysis. The well‐known element free Galerkin mesh‐free method is adopted. The partition of unity property of MLS shape functions allows for the field variables to be extrinsically enriched by appropriate functions that introduce existing discontinuities in the solution field. Enrichment of the main unknowns including solid displacement, water phase pressure, and gas phase pressure are accounted for, and a suitable enrichment strategy for different discontinuity types are discussed. In the case of weak discontinuity, the enrichment technique previously used by Krongauz and Belytschko [Int. J. Numer. Meth. Engng., 1998; 41:1215–1233] is selected. As these functions possess discontinuity in their first derivatives, they can be used for modeling material interfaces, generating only minor oscillations in derivative fields (strain and pressure gradients for multiphase porous media), as opposed to unenriched and constrained mesh‐free methods. Different problems of multimaterial poro‐elasticity including fully saturated, partially saturated one, and two‐phase flows under the assumption of fully coupled extended formulation of Biot are examined. As a further development, problems involved with both material interface and impermeable discontinuities, where no fluid exchange is permitted across the discontinuity, are considered and numerically discussed. Copyright © 2014 John Wiley & Sons, Ltd.     

Analysis of shield tunnel

This paper proposes a two‐dimensional finite element model for the analysis of shield tunnels by taking into account the construction process which is divided into four stages. The soil is assumed to behave as an elasto‐plastic medium whereas the shield is simulated by beam–joint discontinuous model in which curved beam elements and joint elements are used to model the segments and joints, respectively. As grout is usually injected to fill the gap between the lining and the soil, the property parameters of the grout are chosen in such a way that they can reflect the state of the grout at each stage. Furthermore, the contact condition between the soil and lining will change with the construction stage, and therefore, different stress‐releasing coefficients are used to account for the changes. To assess the accuracy that can be attained by the method in solving practical problems, the shield tunnelling in the No. 7 Subway Line Project in Osaka, Japan, is used as a case history for our study. The numerical results are compared with those measured in the field. The results presented in the paper show that the proposed numerical procedure can be used to effectively estimate the deformation, stresses and moments experienced by the surrounding soils and the concrete lining segments. The analysis and method presented in this paper can be considered to be useful for other subway construction projects involving shield tunnelling in soft soils. Copyright © 2004 John Wiley & Sons, Ltd.     

Computational model for the simulation of the shield tunneling process in cohesive soils

A two-dimensional computational model is developed here in order to simulate the continuous advance of the Earth Pressure Balance (EPB) Shield during the tunneling process in cohesive soils. The model is based on the combination of the plane strain “transverse–longitudinal” sections that can incorporate the three-dimensional deformation of the soil around and ahead of the shield face. This model is capable of prediciting the soil response due to the shield tunneling before the event, especially in soft ground conditions. An elasto-plastic finite element analysis that is based on the coupled theory of mixtures for inelastic porous media for finite deformation is used in this work to describe the time-dependent deformation of the saturated cohesive soils. The results of this model are compared with the in situ field measurements of the N-2 tunnel project excavated in 1981 in San Francisco using the EPB shield tunneling machine. Reasonable agreement is found between the observed field measurements and the predicted deformations of the soil using the proposed numerical simulation. Copyright © 1999 John Wiley & Sons, Ltd.     

The Interaction Between Shield,Ground and Tunnel Support in TBM Tunnelling Through Squeezing Ground

When planning a TBM drive in squeezing ground, the tunnelling engineer faces a complex problem involving a number of conflicting factors. In this respect, numerical analyses represent a helpful decision aid as they provide a quantitative assessment of the effects of key parameters. The present paper investigates the interaction between the shield, ground and tunnel support by means of computational analysis. Emphasis is placed on the boundary condition, which is applied to model the interface between the ground and the shield or tunnel support. The paper also discusses two cases, which illustrate different methodical approaches applied to the assessment of a TBM drive in squeezing ground. The first case history—the Uluabat Tunnel (Turkey)—mainly involves the investigation of TBM design measures aimed at reducing the risk of shield jamming. The second case history—the Faido Section of the Gotthard Base Tunnel (Switzerland)—deals with different types of tunnel support installed behind a gripper TBM.     

非饱和土化学-塑性耦合本构行为的数值模拟

基于Hueckel提出的饱和黏土化学-塑性本构模型和Gallipoli提出的非饱和土弹塑性本构模型,提出了一个新的非饱和多孔介质的化学-塑性本构模型,并建立了该模型的隐式积分算法,算法中考虑了化学软化和非饱和吸力的影响。在已有的非饱和多孔介质有限元分析程序平台上进行了程序研发,对孔隙水中化学污染物浓度变化对非饱和土力学行为的影响进行数值模拟,使所研制的程序能够进行岩土工程问题的化学-力学耦合非线性分析。     

Cracking risk of partially saturated porous media—Part I: Microporoelasticity model

Drying of deformable porous media results in their shrinkage, and it may cause cracking provided that shrinkage deformations are hindered by kinematic constraints. This is the motivation to develop a thermodynamics‐based microporoelasticity model for the assessment of cracking risk in partially saturated porous geomaterials. The study refers to 3D representative volume elements of porous media, including a two‐scale double‐porosity material with a pore network comprising (at the mesoscale) 3D mesocracks in the form of oblate spheroids, and (at the microscale) spherical micropores of different sizes. Surface tensions prevailing in all interfaces between solid, liquid, and gaseous matters are taken into account. To establish a thermodynamics‐based crack propagation criterion for a two‐scale double‐porosity material, the potential energy of the solid is derived, accounting—in particular—for mesocrack geometry changes (main original contribution) and for effective micropore pressures, which depend (due to surface tensions) on the pore radius. Differentiating the potential energy with respect to crack density parameter yields the thermodynamical driving force for crack propagation, which is shown to be governed by an effective macrostrain. It is found that drying‐related stresses in partially saturated mesocracks reduce the cracking risk. The drying‐related effective underpressures in spherical micropores, in turn, result in a tensile eigenstress of the matrix in which the mesocracks are embedded. This way, micropores increase the mesocracking risk. Model application to the assessment of cracking risk during drying of argillite is the topic of the companion paper (Part II). Copyright © 2009 John Wiley & Sons, Ltd.     

A discrete numerical approach for modeling face stability in slurry shield tunnelling in soft soils

In this paper, a numerical simulation method for evaluating tunnelling-induced ground movement is presented. The method involves discrete element simulation of TBM slurry shield advancement and considers explicitly soil excavation from the face, effects of varying face support pressure, and the influence of tunnel cover depth. For the cases studied, it is found that for tunnel cover depths (C/D) between 0.7 and 2.1, ground deformations inducing by the tunnelling can be controlled within a certain extent and tunnel face stability can ensured, provided the support pressure ratio (N) lies between 0.8 and 1.5. The proposed method is reasonably benefited to modeling the face stability in shield-driven tunnels in soft soils.     

非饱和粘性土中气体渗透特征

由于气体的进入形成了非饱和带,饱和土层中的气体渗透实质上是非饱和渗透问题。本文借助基于多孔介质气体渗流理论并考虑气体压缩性的、描述非饱和土中气体渗流运动的本构方程,试验研究了上海地区饱和粘性土层中的气体渗透规律。结果表明,饱水粘性土层中,气体渗透存在起始压力值,数值为100～200kPa, 与上海地区粘土的非饱和进气值较为接近;气体渗透速度与外界所施加的压力平方差梯度之间存在明显的分段特征。后者可能与在充满结合水的微孔隙通道中,只有一部分结合水在气压梯度的作用下发生流动,从而形成气体通道,施加的外力越大,形成的气体通道也越大,气体渗流也越明显。非饱和粘性土中气体渗透规律的研究,对于饱和粘土中的气压法施工,具有重要的理论与工程实践意义。     

A cap model for partially saturated soils

The present paper deals with the extension of a cap model in order to describe the material behavior of partially saturated soils, in particular, of partially saturated sands and silts. The soil model is formulated in terms of two stress state variables, using net stress and matric suction and, alternatively, the average soil skeleton stress and suction, the latter playing the role of a stress‐like plastic internal variable. The yield surface, consisting of a shear failure surface and a hardening cap surface, the plastic potentials for the non‐associated flow rule and the hardening law for the cap are extended by taking into account the effects of matric suction on the material behavior. Furthermore, the third invariant of the deviatoric stress tensor is taken into account in the formulation of the yield surfaces. The developed model is validated by the numerical simulation of an extensive series of suction controlled tests for a silty sand, which were conducted at different constant values of suction. Although both versions of the soil model yield identical results for stress paths at constant values of matric suction, differences are encountered for stress paths involving wetting. Copyright © 2007 John Wiley & Sons, Ltd.     

Analysis of liquefaction susceptibility of nearly saturated sands

A new constitutive formulation for simulating the behaviour of nearly saturated sands under seismic loads is presented. The formulation is based on combining the Henry's law for dissolution of gas in water, the ideal or perfect gas law and the law of conservation of mass. The effects of transient air dissolution in water on the compressibility of partially saturated soils are also taken into account. The model was calibrated based on numerical simulations of isotropically consolidated cyclic triaxial tests conducted on partially saturated samples of Toyoura sand. A multi‐yield plasticity soil constitutive model implemented in the finite element code DYNAFLOW was used for these numerical simulations. It is shown that the formulation proposed here is able to reasonably predict the soil cyclic undrained behaviour at various degrees of saturation (95% and higher). Copyright © 2006 John Wiley & Sons, Ltd.     

One‐dimensional contaminant transport through a deforming porous medium: theory and a solution for a quasi‐steady‐state problem

Contaminant migration through soil is usually modelled mathematically using the dispersion–advection equation. This type of model finds application when planning the remediation of contaminated land, predicting the movement of polluted groundwater and designing engineered landfills. Usually the analysis assumes that the porous media through which the contaminant migrates is stationary. However, the construction of landfills on clay soils means that the soil beneath the landfill will undergo time‐dependent deformation as the soil consolidates. To date, there are no published data on the effect a deforming porous media may have on contaminant transport beneath a landfill; indeed, there appears to be no theory of contaminant migration through a deforming soil. In this paper, a one‐dimensional theory of contaminant migration through a saturated deforming porous media is developed based on a small and large strain analysis of a consolidating soil and conservation of contaminant mass. By selection of suitable parameters, the new transport equation reduces to the familiar one‐dimensional dispersion–advection equation for a saturated soil with linear, reversible, equilibrium controlled sorption of the contaminant onto the soil skeleton. Analytic solutions to a quasi‐steady‐state contaminant transport problem for a deforming media are presented, and a preliminary assessment made of the potential importance of soil deformation on the results of a contaminant migration analysis. Copyright © 2000 John Wiley & Sons, Ltd.     

考虑动态泥膜效应的盾构开挖面稳定理论研究

泥水盾构工法被广泛应用于高压富水复杂地质环境中修筑隧道,其中开挖面稳定控制是工程成败的关键,在盾构掘进时,由于刀盘不断旋转切削地层,开挖面上的泥膜处于“形成―破坏―再形成”的动态循环过程,泥膜透水性增强,这种动态泥膜作用下开挖面的稳定如何评价有待于进一步探索。对此,对泥水盾构掘进过程中的动态泥膜进行了分类,提出了考虑盾构运动和泥浆渗滤特征的动态泥膜理论。在此基础上,采用空间离散法将泥水压力引入旋转体模型,建立了考虑动态泥膜效应的开挖面稳定力学模型,提出了泥水盾构动态掘进过程中开挖面临界泥浆压力及泥浆支护效率评估方法,研究指出：当盾构在强渗透性地层中掘进时,在尽量提高泥浆成膜率的同时,应注重泥浆黏度调节,充分发挥渗透力支护作用,同时可适当提高盾构掘进速度并降低刀盘转速,有利于开挖面稳定控制。以泥浆成膜率、等效渗透系数为依据给出了泥浆支护效率设计图,研究工作对盾构掘进支护压力控制和泥浆配置有一定的借鉴意义。     

Correction to the gas phase pressure term in the continuum model for partially saturated granular media presented by Pietruszczak and co‐workers

Pietruszczak and coworkers (Internat. J. Numer. Anal. Methods Geomech. 1994; 18 (2):93–105; Comput. Geotech. 1991; 12 ( ):55–71) have presented a continuum‐based model for predicting the dynamic mechanical response of partially saturated granular media with viscous interstitial liquids. In their model they assume that the gas phase is distributed uniformly throughout the medium as discrete spherical air bubbles occupying the voids between the particles. However, their derivation of the air pressure inside these gas bubbles is inconsistent with their stated assumptions. In addition the resultant dependence of gas pressure on liquid saturation lies outside of the plausible range of possible values for discrete air bubbles. This results in an over‐prediction of the average bulk modulus of the void phase. Corrected equations are presented. Copyright © 2003 John Wiley & Sons, Ltd.     

大直径盾构进出洞加固体稳定性判别方法研究

从监测指标和监测范围两个方面对钱江隧道盾构进出洞加固体稳定性的实时判别进行了研究,即首先在分析板块理论和弹塑性静力学理论的基础上,提出能实时反映盾构进出洞加固体稳定性的监测项目,包括水土压力P及其相应的稳定性评价指标 、 ,而后采用理论分析和试验研究加固体的破坏模式,在此基础上,提出加固体实时监测点布置方案。通过监测点的合理布置,监测加固体所受水土压力P,计算 、 ,可实现加固体稳定性的实时评价     

The effective stress concept and the evaluation of changes in pore air pressure under jacketed isotropic compression tests for dry sand based on the 2‐phase porous theory

The effective stress concept for solid‐fluid 2‐phase media was revisited in this work. In particular, the effects of the compressibility of both the pore fluid and the soil particles were studied under 3 different conditions, i.e., undrained, drained, and unjacketed conditions based on a Biot‐type theory for 2‐phase porous media. It was confirmed that Terzaghi effective stress holds at the moment when soil grains are assumed to be incompressible and when the compressibility of the pore fluid is small enough compared to that of the soil skeleton. Then, isotropic compression tests for dry sand under undrained conditions were conducted within the triaxial apparatus in which the changes in the pore air pressure could be measured. The ratio of the increment in the cell pressure to the increment in the pore air pressure, m, corresponds to the inverse of the B value by Bishop and was obtained during the step loading of the cell pressure. In addition, the m values were evaluated by comparing them with theoretically obtained values based on the solid‐fluid 2‐phase mixture theory. The experimental m values were close to the theoretical values, as they were in the range of approximately 40 to 185, depending on the cell pressure. Finally, it was found that the soil material with a highly compressible pore fluid, such as air, must be analyzed with the multi‐phase porous mixture theory. However, Terzaghi effective stress is practically applicable when the compressibilities of both the soil particles and the pore fluid are small enough compared to that of the soil skeleton.     

Frictionless contact formulation for dynamic analysis of nonlinear saturated porous media based on the mortar method

A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd.