盾构隧道施工引起地表下土体变位的分析评估

目前对盾构隧道施工引起的地表下土体变位进行预测分析的研究还相对较少,尤其是缺乏简单实用的工程估算方法。在深埋隧道周围土体弹性位移计算方法及盾构间隙参数研究成果的基础上,通过相应的分析和假设提出了一种预测盾构隧道施工引起地表下土体水平变位的简便估算方法。另外,在已有研究成果基础上,提出了一种新的盾构隧道沉降槽的描述方法,并结合Mair等提出的计算公式对盾构隧道施工引起的地表下土体沉降变位进行预测。通过与有限元计算结果及一些典型盾构隧道监测数据的对比分析,证明提出的估算方法能够较好地预测实际工程中盾构隧道施工引起的地表下土体的变位情况。     

A simplified analysis method for piled raft foundations subjected to ground movements induced by tunnelling

A simplified method of numerical analysis has been developed to estimate the deformation and load distribution of piled raft foundations subjected to ground movements induced by tunnelling and incorporated into a computer program ‘PRAB’. In this method, a hybrid model is employed in which the flexible raft is modelled as thin plates, the piles as elastic beams, and the soil is treated as interactive springs. The interactions between structural members, pile–soil–pile, pile–soil–raft and raft–soil–raft interactions, are modelled based on Mindlin's solutions for both vertical and lateral forces. The validity of the proposed method is verified through comparisons with some published solutions for single piles and pile groups subjected to ground movements induced by tunnelling. Thereafter, the solutions from this approach for the analysis of a pile group and a piled raft subjected to ground movements induced by tunnelling are compared with those from three‐dimensional finite difference program. Good agreements between these solutions are demonstrated. The method is then used for a parametric study of single piles, pile groups and piled rafts subjected to ground movements induced by tunnelling. Copyright © 2005 John Wiley & Sons, Ltd.     

Axial kinematic response of end‐bearing piles to P waves

Kinematic pile–soil interaction under vertically impinging seismic P waves is revisited through a novel continuum elastodynamic solution of the Tajimi type. The proposed model simulates the steady‐state kinematic response of a cylindrical end‐bearing pile embedded in a homogeneous viscoelastic soil stratum over a rigid base, subjected to vertically propagating harmonic compressional waves. Closed‐form solutions are obtained for the following: (i) the displacement field in the soil and along the pile; (ii) the kinematic Winkler moduli (i.e., distributed springs and dashpots) along the pile; (iii) equivalent, depth‐independent, Winkler moduli to match the motion at the pile head. The solution for displacements is expressed in terms of dimensionless transfer functions relating the motion of the pile head to the free‐field surface motion and the rock motion. It is shown that (i) a pile foundation may significantly alter (possibly amplify) the vertical seismic excitation transmitted to the base of a structure and (ii) Winkler moduli pertaining to kinematic loading differ from those for inertial loading. Simple approximate expressions for kinematic Winkler moduli are derived for use in applications. Copyright © 2013 John Wiley & Sons, Ltd.     

Analysis of soil–pile–structure interaction in a two‐layer ground during earthquakes considering liquefaction

This study is conducted with a numerical method to investigate the seismic behaviour among certain soils, single piles, and a structure. A series of numerical simulations of the seismic behaviour of a single‐pile foundation constructed in a two‐layer ground is carried out. Various sandy soils, namely, dense sand, medium dense sand, reclaimed soil, and loose sand, are employed for the upper layer, while one type of clayey soil is used for the lower layer. The results reveal that when a structure is built in a non‐liquefiable ground, an amplification of the seismic waves is seen on the ground surface and in the upper structure, and large bending moments are generated at the pile heads. When a structure is built in a liquefiable ground, a de‐amplification of the seismic waves is seen on the ground surface and in the upper structure, and large bending moments are generated firstly at the pile heads and then in the lower segment at the boundary between the soil layers when liquefaction takes place. Copyright © 2007 John Wiley & Sons, Ltd.     

Beam–solid contact formulation for finite element analysis of pile–soil interaction with arbitrary discretization

This paper presents an embedded beam formulation for discretization independent finite element (FE) analyses of interactions between pile foundations or rock anchors and the surrounding soil in geotechnical and tunneling engineering. Piles are represented by means of finite beam elements embedded within FEs for the soil represented by 3D solid elements. The proposed formulation allows consideration of piles and pile groups with arbitrary orientation independently from the FE discretization of the surrounding soil. The interface behavior between piles and the surrounding soil is represented numerically by means of a contact formulation considering skin friction as well as pile tip resistance. The pile–soil interaction along the pile skin is considered by means of a 3D frictional point‐to‐point contact formulation using the integration points of the beam elements and reference points arbitrarily located within the solid elements as control points. The ability of the proposed embedded pile model to represent groups of piles objected to combined axial and shear loading and their interactions with the surrounding soil is demonstrated by selected benchmark examples. The pile model is applied to the numerical simulation of shield driven tunnel construction in the vicinity of an existing building resting upon pile foundation to demonstrate the performance of the proposed model in complex simulation environments. Copyright © 2014 John Wiley & Sons, Ltd.     

基于间隙参数模型的盾构隧道周围土体位移分析

基于一种盾构加载模型并引进间隙参数的概念,对盾构隧道周围土体位移的机理及预测方法进行了总结。盾构相对于设计轴线的偏移是影响土体变形的主要原因之一。确定间隙参数是应用解析方法的关键步骤,提出了一个半经验、半解析方法,并针对一个工程实例把现场测量值和理论计算值进行了对比分析。     

挡墙水平变位诱发地表沉降的显式解析解

将基坑看作平面应变边值问题,借助挡墙刚性平移诱发地表沉降的基本解析解,利用微积分思想,推导得到挡墙刚性变位和柔性变位下的地表沉降显式解析解。其中,在挡墙刚性平移影响范围无穷大时,显示解析解与基本解析解形式一致。通过与挡墙刚性平移解析解和绕墙趾转动、绕墙顶转动、三角形变位、抛物线柔性变位模式等几种已有的积分形式解析解作对比,验证了显式解析解的正确性。将显式解析解应用于实际工程,通过与实测数据对比分析,对解析解的适用范围进行了探讨。当围护结构水平位移较小时,显式解析解可较好地预测墙后地表沉降;当围护结构水平位移较大时,由显式解析解计算得到的墙后地表沉降归一化曲线可偏安全地估计墙后地表沉降,说明了该显式解析解的工程实用性。     

A non‐linear coupled finite element–boundary element model for the prediction of vibrations due to vibratory and impact pile driving

This paper presents a non‐linear coupled finite element–boundary element approach for the prediction of free field vibrations due to vibratory and impact pile driving. Both the non‐linear constitutive behavior of the soil in the vicinity of the pile and the dynamic interaction between the pile and the soil are accounted for. A subdomain approach is used, defining a generalized structure consisting of the pile and a bounded region of soil around the pile, and an unbounded exterior linear soil domain. The soil around the pile may exhibit non‐linear constitutive behavior and is modelled with a time‐domain finite element method. The dynamic stiffness matrix of the exterior unbounded soil domain is calculated using a boundary element formulation in the frequency domain based on a limited number of modes defined on the interface between the generalized structure and the unbounded soil. The soil–structure interaction forces are evaluated as a convolution of the displacement history and the soil flexibility matrices, which are obtained by an inverse Fourier transformation from the frequency to the time domain. This results in a hybrid frequency–time domain formulation of the non‐linear dynamic soil–structure interaction problem, which is solved in the time domain using Newmark's time integration method; the interaction force time history is evaluated using the θ‐scheme in order to obtain stable solutions. The proposed hybrid formulation is validated for linear problems of vibratory and impact pile driving, showing very good agreement with the results obtained with a frequency‐domain solution. Linear predictions, however, overestimate the free field peak particle velocities as observed in reported field experiments during vibratory and impact pile driving at comparable levels of the transferred energy. This is mainly due to energy dissipation related to plastic deformations in the soil around the pile. Ground vibrations due to vibratory and impact pile driving are, therefore, also computed with a non‐linear model where the soil is modelled as an isotropic elastic, perfectly plastic solid, which yields according to the Drucker–Prager failure criterion. This results in lower predicted free field vibrations with respect to linear predictions, which are also in much better agreement with experimental results recorded during vibratory and impact pile driving. Copyright © 2008 John Wiley & Sons, Ltd.     

Elastoplastic medium for foundation settlements and monotonic soil–structure interaction under combined loadings

Foundation settlements and soil–structure interaction are important problems to structural and geotechnical engineers. This study introduces a novel elastoplastic three‐degree‐of‐freedom medium which models foundations settlements under combined loadings. A soil–structure interaction problem can then be solved by replacing the soil mass with this three‐degree‐of‐freedom elastoplastic medium, thus reducing significantly the size of the problem. The model was developed by extending the classical plasticity concepts to the force‐deformation level. Its ability to predict foundation deformations was evaluated using finite element solutions of a typical shallow foundation problem and was found reasonably accurate while producing significant time savings. Copyright © 2006 John Wiley & Sons, Ltd.     

盾构偏航引起的地表位移预测

对于盾构隧道的设计和施工，解析方法预测地表沉降，具有使用简便、物理意义明确的特点。在Sagaseta提出的基于由地层损失引起的地表沉降理论基础上，结合工程实际情况，提出了模拟盾构推进过程的三维地层损失模式，并推导了相应的地表位移计算解析公式，对实际工程具有参考价值。     

Frictional contact algorithms in SPH for the simulation of soil–structure interaction

Simulation of frictional contact between soils and rigid or deformable structure in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two algorithms are implemented into the SPH code to describe contact behavior, where the contact forces are calculated using the law of conservation of momentum based on ideal plastic collision or using the criteria of partial penetrating. In both algorithms, the problem of boundary deficiency inherited from SPH is properly handled so that the particles located at contact boundary can have precise acceleration, which is critical for contact detection. And the movement and rotation of the rigid structure are taken into account so that it is easy to simulate the process of pile driving or movement of a retaining wall in geotechnical engineering analysis. Furthermore, the capability of modeling deformability of a structure during frictional contact simulations broadens the fields of SPH application. In contrast to previous work dealing with contact in SPH, which usually use particle‐to‐particle contact or ignoring sliding between particles and solid structure, the method proposed here is more efficient and accurate, and it is suitable to simulate interaction between soft materials and rigid or deformable structures, which are very common in geotechnical engineering. A number of numerical tests are carried out to verify the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions or results from finite element method analysis. Good agreement obtained from these comparisons suggests that the proposed algorithms are robust and can be applied to extend the capability of SPH in solving geotechnical problems. Copyright © 2013 John Wiley & Sons, Ltd.     

Analytical and numerical analyses of the load‐bearing capacity of retaining walls laterally supported at both ends

This paper investigates the load‐bearing capacity of a perfectly smooth retaining wall laterally supported at both ends assuming that the wall fails by the development of three plastic hinges. The study considers the case of a cohesionless elastic–perfectly plastic backfill with a Mohr–Coulomb yield criterion and an associative flow rule in drained conditions. A kinematically admissible soil–structure failure mechanism is proposed and compared with the conventional solutions and with results from a numerical finite element modelling. The study shows that the proposed solution and the numerical solution are in good agreement. These solutions are found to be much more favourable for the wall than the conventional solutions. Copyright © 2010 John Wiley & Sons, Ltd.     

Implementation of DSC model and application for analysis of field pile tests under cyclic loading

The disturbed state concept (DSC) model, and a new and simplified procedure for unloading and reloading behavior are implemented in a nonlinear finite element procedure for dynamic analysis for coupled response of saturated porous materials. The DSC model is used to characterize the cyclic behavior of saturated clays and clay–steel interfaces. In the DSC, the relative intact (RI) behavior is characterized by using the hierarchical single surface (HISS) plasticity model; and the fully adjusted (FA) behavior is modeled by using the critical state concept. The DSC model is validated with respect to laboratory triaxial tests for clay and shear tests for clay‐steel interfaces. The computer procedure is used to predict field behavior of an instrumented pile subjected to cyclic loading. The predictions provide very good correlation with the field data. They also yield improved results compared to those from a HISS model with anisotropic hardening, partly because the DSC model allows for degradation or softening and interface response. Copyright © 2000 John Wiley & Sons, Ltd.     

Performance of seismically loaded shearwalls on nonlinear shallow foundations

If utilized, the energy dissipative capability of seismically loaded shallow foundations due to inelastic behavior can result in more economic design, provided the consequences, such as excessive deformations are accounted for. In this article, a Beam‐on‐Nonlinear‐Winkler‐Foundation (BNWF) model is used to assess the performance of shearwall‐foundation systems with different attributes, when subjected to ground motions of varied hazard levels. The numerical study indicates that the force and drift demands to the shearwall reduce significantly, when nonlinear foundation behavior is realized, while permanent settlement is well below the permissible limit. These results support the concept of shallow foundation capacity mobilization in future design. Copyright © 2010 John Wiley & Sons, Ltd.     

盾构施工中土体损失引起的地面沉降预测

土质软硬决定了隧道周围土体的移动方向,移动焦点在隧道中心点与隧道底部位置之间变动。采用两圆相切的土体损失模型,通过引入移动焦点的坐标参数,建立了统一的土体移动模型,该模型能将Park模型与Loganathan模型包括在内。假定土体不排水,利用源汇法推导了由土体损失引起的地面沉降通用计算公式,该方法适用于施工阶段。算例分析表明该方法的计算结果与实测值非常吻合,适用于各种土质条件。Loganathan公式只适用于土质较差的情况,当土质较好时计算得到的地面沉降量要比实测值偏小。     

金属矿山崩落采矿法引起的岩层移动规律分析

以典型陡倾结构面条件下的金属矿山--程潮铁矿西区为例,通过对矿区的地表变形监测资料及宏观破坏特征分析,认为矿区的岩层移动分为2个阶段,第1阶段为采空区顶板岩体破坏扩展至地表引起塌陷阶段;第2阶段为采空区周边围岩向采空区的倾倒破坏阶段,并得出了倾倒滑移区的地表岩体变形规律：岩体主要发生水平移动,水平移动值大于沉降值;变形先以缓慢变形为主,然后进入一个快速变形阶段,存在明显的转折点;开采沉陷和地形引起的应力同向叠加作用,使得地表岩体沿下坡方向的变形值增大,特别是水平移动值。同时揭示了矿区岩层移动角的分布特征;南部岩层移动角大于北部,究其原因是北部受最为发育的NNW、NNE结构面影响,倾倒破坏较为严重。所得成果为其他类似的金属矿山工程提供可借鉴的规律。     

抽水地面沉降预计的随机介质模型

本文视抽水引起的地面沉降为一随机过程,应用随机介质理论和土力学基本原理,建立了随机介质模型,对抽水地面沉降及变形分布进行了分析预测,计算实例表明,该模型是有铲的。     

Transverse seismic kinematics of single piles by a modified Vlasov model

Within the framework of soil–pile interaction, a novel displacement scheme for the transverse kinematic response of single piles to vertically propagating S waves is proposed on the basis of the modified Vlasov foundation model. The displacement model contains a displacement function along the pile axis and an attenuation function along the radial direction. The governing equations and boundary conditions of the two undetermined functions are obtained in a coupled form by using Hamilton's principle. An iterative algorithm is adopted to decouple and solve the two unknown functions. In light of the governing equation of the pile kinematics, a mechanical model is proposed to evaluate the present method on a physical basis considering material damping. The coefficient of the equivalent Winkler spring is derived explicitly as function of the displacement decay parameter γ and soil Poisson's ratio. A parametric study is performed to investigate the effects of the soil–pile system properties on the kinematic response of single piles. The results show that the dimensionless pile length controls the transverse kinematics of piles. In terms of the theory of beams on elastic foundation, the classification limits of the dimensionless pile length may be π ∕ 4 and π, respectively. Copyright © 2014 John Wiley & Sons, Ltd.     

Load transfer approach for laterally loaded piles

A two‐parameter model has been proposed previously for predicting the response of laterally loaded single piles in homogenous soil. A disadvantage of the model is that at high Poisson's ratio, unreliable results may be obtained. In this paper, a new load transfer approach is developed to simulate the response of laterally loaded single piles embedded in a homogeneous medium, by introducing a rational stress field. The approach can overcome the inherent disadvantage of the two‐parameter model, although developed in a similar way. Generalized solutions for a single pile and the surrounding soil under various pile‐head and base conditions were established and presented in compact forms. With the solutions, a load transfer factor, correlating the displacements of the pile and the soil, was estimated and expressed as a simple equation. Expressions were developed for the modulus of subgrade reaction for a Winkler model as a unique function of the load transfer factor. Simple expressions were developed for estimating critical pile length, maximum bending moment, and the depth at which the maximum moment occurs. All the newly established solutions and/or expressions, using the load transfer factor, offer satisfactory predictions in comparison with the available, more rigorous numerical approaches. The current solutions are applicable to various boundary conditions, and any pile–soil relative stiffness. Copyright © 2001 John Wiley & Sons, Ltd.     

Elastic models for nonlinear response of rigid passive piles

Recent study indicates that the response of rigid passive piles is dominated by elastic pile–soil interaction and may be estimated using theory for lateral piles. The difference lies in that passive piles normally are associated with a large scatter of the ratio of maximum bending moment over maximum shear force and induce a limiting pressure that is ~1/3 that on laterally loaded piles. This disparity prompts this study. This paper proposes pressure‐based pile–soil models and develops their associated solutions to capture response of rigid piles subjected to soil movement. The impact of soil movement was encapsulated into a power‐law distributed loading over a sliding depth, and load transfer model was adopted to mimic the pile–soil interaction. The solutions are presented in explicit expressions and can be readily obtained. They are capable of capturing responses of model piles in a sliding soil owing to the impact of sliding depth and relative strength between sliding and stable layer on limiting force prior to ultimate state. In comparison with available solutions for ultimate state, this study reveals the 1/3 limiting pressure (of the active piles) on passive piles was induced by elastic interaction. The current models employing distributed pressure for moving soil are more pertinent to passive piles (rather than plastic soil flow). An example calculation against instrumented model piles is provided, which demonstrates the accuracy of the current solutions for design slope stabilising piles. Copyright © 2014 John Wiley & Sons, Ltd.