多向荷载作用下层状地基刚性桩筏基础分析

提出一种多向荷载作用下层状地基中刚性桩筏基础的计算方法。基于剪切位移法,采用传递矩阵形式分析了竖向荷载下桩顶面-桩顶面相互作用;引入修正桩侧地基模量,采用有限差分法分析了水平荷载下桩顶面-桩顶面相互作用;基于层状弹性半空间理论,分析了多向荷载下桩顶面-土表面、土表面-桩顶面、土表面-土表面的相互作用关系。建立了桩土体系柔度矩阵,得到了多向荷载下层状地基中刚性桩筏基础的受力和变形的关系以及桩的内力和变形沿桩身分布规律。通过与有限元对比,验证了该方法的合理性和修正地基模量的优越性,并对多向荷载作用下的桩筏基础进行了计算分析,计算结果表明,水平力将会引起桩筏基础的倾斜。     

Analysis of pile‐raft foundations under complex loads in layered soils

Considering there is hardly any concerted effort to analyze the pile‐raft foundations under complex loads (combined with vertical loads, horizontal loads and moments), an analysis method is proposed in this paper to estimate the responses of pile‐raft foundations which are subjected to vertical loads, horizontal loads and moments in layered soils based on solutions for stresses and displacements in layered elastic half space. Pile to pile, pile to soil surface, soil surface to pile and soil surface to soil surface interactions are key ingredients for calculating the responses of pile‐raft foundations accurately. Those interactions are fully taken into account to estimate the responses of pile‐raft foundations subject to vertical loads, horizontal loads and moments in layered soils. The constraints of the raft on vertical movements, horizontal movements and rotations of the piles as well as the constraints of the raft on vertical movements and horizontal movements of the soils are considered to reflect the coupled effect on the raft. The method is verified through comparisons with the published methods and FEM. Then, the method is adopted to investigate the influence of soil stratigraphy on pile responses. The study shows that it is necessary to consider the soil non‐homogeneity when estimating the responses of pile‐raft foundations in layered soils, especially when estimating the horizontal responses of pile‐raft foundations. The horizontal loads and the moments have a significant impact on vertical responses of piles in pile‐raft foundations, while vertical loads have little influence on horizontal responses of piles in pile‐raft foundations in the cases of small deformations. The proposed method can provide a simple and useful tool for engineering design. Copyright © 2013 John Wiley & Sons, Ltd.     

A simplified analysis method for piled raft foundations in non‐homogeneous soils

A simplified method of numerical analysis based on elasticity theory has been developed for the analysis of axially and laterally loaded piled raft foundations embedded in non‐homogeneous soils 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 and the piles as elastic beams and the soil is treated as springs. The interactions between structural members, pile–soil–pile, pile–soil–raft and raft–soil–raft interactions, are approximated based on Mindlin's solutions for both vertical and lateral forces with consideration of non‐homogeneous soils. The validity of the proposed method is verified through comparisons with some published solutions for single piles, pile groups and capped pile groups in non‐homogeneous soils. Thereafter, the solutions from this approach for the analysis of axially and laterally loaded 4‐pile pile groups and 4‐pile piled rafts embedded in finite homogeneous and non‐homogeneous soil layers are compared with those from three‐dimensional finite element analysis. Good agreement between the present approach and the more rigorous finite element approach is demonstrated. Copyright © 2002 John Wiley & Sons, Ltd.     

Elastoplastic multiphase model for simulating the response of piled raft foundations subject to combined loadings

A multiphase model and corresponding computational time‐saving finite element code is proposed in this paper for predicting the settlements experienced by a piled raft foundation when subject to the combined action of vertical and lateral loadings. This model, which is formulated in the framework of an elastoplastic behaviour for the soil and the reinforcing piles as well, explicitly accounts for the shear and flexural behaviour of the latter. Starting from a simple analytical example where all the concepts attached to this model are clearly illustrated, the main stages leading to its finite element implementation are then presented. The numerical tool thus elaborated, is applied to the simulation of a pile‐reinforced strip foundation submitted to a horizontally applied seismic load in addition to a permanent vertical load. One of the key results of such a simulation in terms of design recommendation, lies in the conclusion that, while the shear and flexural contributions of the reinforcement play quite a negligible role in the case of a vertical load (as compared with their axial resistance), they remain absolutely essential for withstanding the seismic lateral loading. Copyright © 2006 John Wiley & Sons, Ltd.     

Vertical response of pile raft foundations subjected to tunneling‐induced ground movements in layered soil

A simplified analysis method has been developed to estimate the vertical movement and load distribution of pile raft foundations subjected to ground movements induced by tunneling based on a two‐stage method. In this method, the Loganathan–Polous analytical solution is used to estimate the free soil movement induced by tunneling in the first stage. In the second stage, composing the soil movement to the pile, the governing equilibrium equations of piles are solved by the finite difference method. The interactions between structural members (such as pile–soil, pile–raft, raft–soil, and pile–pile) are modeled based on the elastic theory method of a layered half‐space. The validity of the proposed method is verified through comparisons with some published solutions for single piles, pile groups, and pile rafts subjected to ground movements induced by tunneling. Good agreements between these solutions are demonstrated. The method is also used for a parametric study to develop a better understanding of the behavior of pile rafts influenced by tunneling operation in layered soil foundations. Copyright © 2011 John Wiley & Sons, Ltd.     

Multiphase model for inclusion‐reinforced geostructures: Application to rock‐bolted tunnels and piled raft foundations

A multiphase model is proposed to describe the mechanical behaviour of geomaterials reinforced by linear inclusions. This macroscopic approach considers the reinforced soil or rock mass as the superposition of continuous media. Equations of motion and constitutive laws of the model are first derived. Its implementation in a finite element computer code is then detailed. A modified implicit algorithm for elastoplastic problems is proposed. The model and its implementation are fully validated for rock‐bolted tunnels (comparison with scale model experiments) and piled raft foundations (comparison with the classical ‘hybrid method’). The Messeturm case history is finally presented to assess the handiness of the approach for real structures. Copyright © 2001 John Wiley & Sons, Ltd.     

软土地基超高层建筑补偿桩筏基础案例再分析

坐落于上海的金茂大厦、上海环球金融中心和上海中心大厦三幢闻名世界的超高层建筑的基础均为补偿深埋桩筏基础。在补偿深埋桩筏基础设计中如何充分考虑补偿基础的优越性,如何考虑深埋基础的抗风能力,是值得研究的问题。利用高层建筑与地基基础共同作用的分析方法,结合统计-经验公式对该三幢大楼的桩筏深基础进行了再分析。通过再分析在节约相当数量的桩之后,其地基承载力和变形仍然能够满足设计要求。同时在再分析中进一步探讨优化设计的能力,论证地下连续墙可以分担相当数量的荷载。相关案例的再分析对进一步完善补偿深埋桩筏基础设计理论提出建议,可为国家节省大量投资。     

群桩与筏板的竖向动力相互作用研究

对沿深度离散形成一系列均质的土层,应用传递矩阵法推导出单桩的动力阻抗。借助推求的桩-桩动力相互作用因子,研究承台板与群桩的竖向动力相互作用,通常假设承台板是刚性的,并充分考虑了承台板具有一定的刚度、可以变形的柔性承台板,提出承台板-群桩动力相互作用模型,然后推出其动力相互作用方程。通过算例说明了承台板的刚度对承台板的动力作用是有影响的。研究结果表明,振动频率越低,柔性承台板的振幅与刚性承台板的振幅相差越大;随着频率的增大,二者差别也减小,并逐渐趋于相等。     

Consolidation analysis of saturated multi‐layered soils with anisotropic permeability caused by a point sink

This paper presents the analytical layer‐element method to analyze the consolidation of saturated multi‐layered soils caused by a point sink by considering the anisotropy of permeability. Starting from the governing equations of the problem, the solutions of displacements and stresses for a single soil layer are obtained in the Laplace–Hankel transformed domain. Then, the analytical layer‐element method is utilized to further derive the solutions for the saturated multi‐layered soils in the transformed domain by combining with the boundary conditions of the soil system and continuity conditions between adjacent layers. The actual solutions in the physical domain can be acquired by the inversion of Laplace–Hankel transform. Numerical results are carried out to show the accuracy and stability of the proposed method and evaluate the influence of sink depth and anisotropic permeability on excess pore pressure and surface settlement. Copyright © 2011 John Wiley & Sons, Ltd.     

Deformation analysis of tunnel excavation below existing pipelines in multi‐layered soils based on displacement controlled coupling numerical method

The effect of tunneling on surrounding environments, especially on existing buried pipelines is a problem that engineers designing and practicing in urban geotechnical environments encounter more frequently than in the past. However, previous studies are usually based on the assumption that the soil is homogeneous. How to reflect soil stratification is the main focus for the problem of tunneling in multi‐layered soils. A displacement controlled coupling numerical method is presented for the displacement analysis of tunnel excavation below existing pipelines in multi‐layered soils. On the basis of the layered soil model, to consider the soil nonhomogeneous characteristic, the finite element method and boundary element method are coupled to simulate the deformation of existing pipelines induced by tunneling. The solutions indicate that good agreements are obtained between the proposed coupling numerical method and the commercial software. The accuracy of the proposed numerical method is better than the two stages method based on the existing closed‐form solutions. Moreover, the results discussed in this paper show that the error obtained by the previous method of weighted average on the basis of homogeneous half space converted from layered soils is not negligible for the obvious difference of elastic parameters among successive layers. Copyright © 2012 John Wiley & Sons, Ltd.     

A quasistatic analysis of a plate on consolidating layered soils by analytical layer‐element/finite element method coupling

In this paper, a coupling method between finite element and analytical layer‐elements is utilized to analyze the time‐dependent behavior of a plate of any shape and finite rigidity resting on layered saturated soils. Based on the integral transform techniques together with the aid of an order reduction method, an analytical layer‐element solution is derived from the governing equations for three‐dimensional Biot consolidation with respect to a Cartesian coordinate system and then extended to be the fundamental solution for the layered saturated soil under a point load. The Mindlin plate is modeled by eight‐noded isoparametric elements. The governing equations of the interaction between soil and plate in the Laplace‐Fourier transformed domain are deduced by referring to the coupling theory of FEM/BEM, and the final solution is obtained by applying numerical inversion. Numerical examples concerned with the time‐dependent response of a plate are performed to demonstrate the influence of soil and plate properties on the interaction process. Copyright © 2014 John Wiley & Sons, Ltd.     

Impedances of rigid cylindrical foundations embedded in transversely isotropic soils

A complete formulation and implementation for assessment of the response to dynamic loads of cylindrical rigid structures embedded in transversely isotropic elastic half‐spaces is presented. The analysis is performed in the frequency domain and the steady‐state structure response is obtained. The method is based on a non‐singular version of the indirect boundary element method which uses influence functions, instead of Green's functions, as fundamental solutions. These influence functions are the response of an elastic half‐space to distributed, internally applied loads. The proposed method imposes full bonding contact between the foundation and the surrounding soil. Numerical results for displacement (vertical and horizontal) and rotation (twisting and rocking) impedances, showing the influence of the soil anisotropy, are presented. Results for the soil–structure interface tractions and for the displacement field throughout the half‐space are also shown. Copyright © 2006 John Wiley & Sons, Ltd.     

Simple superposition approach for dynamic analysis of piled embedded footings

The effectiveness and accuracy of the superposition method in assessing the dynamic stiffness and damping coefficients (impedance functions) of embedded footings supported by vertical piles in homogeneous viscoelastic soil is addressed. To this end, the impedances of piled embedded footings are compared to those obtained by superposing the impedance functions of the corresponding pile groups and embedded footings treated separately, with the magnitude of the relative average differences being around 10–30%. The results are presented in a set of dimensionless graphs and simple expressions that can be used to estimate the dynamic stiffness and damping of piled embedded footings, provided that the impedance functions of the two individual components are known. This is precisely the reason why the superposition approach studied here is appealing, because such impedance functions for both embedded footings and pile groups are available for a wide range of cases. How to estimate the kinematic response functions of the system when those of the individual components are known is also discussed. To address the problem, parametric analyses performed using a 3D frequency‐domain elastodynamic BEM‐FEM formulation are presented for different pile–soil stiffness contrasts, embedment depths, pile‐to‐pile separations and excitation frequencies. Vertical, horizontal, rocking, and cross‐coupled horizontal‐rocking impedance functions, together with translational and rotational kinematic response functions, are discussed. The results suggest that the superposition concept, in conjunction with a correction strategy as that presented herein, can be employed in geotechnical design. For kinematic effects, the response functions of the embedded footing are found to provide reasonable estimates of the system's behaviour. Copyright © 2011 John Wiley & Sons, Ltd.     

下卧基岩饱和地基中埋置刚性基础的竖向振动问题研究

基于Biot提出的饱和土波动方程,研究了埋置于单层饱和地基中的有质量的刚性圆柱基础的竖向振动问题。运用Hankel积分变换求解饱和土基本动力方程,并考虑基础与地基接触面的混合边值条件,求得了基础底面和基础侧面的动反力,结合基础振动的动力平衡方程,得到了基础竖向振动振幅和地基等效动力刚度的表达式。数值分析结果表明：下卧基岩饱和地基的等效动力刚度和基础振动振幅都呈现出明显的波动现象,土层厚度、基础埋深比、基础质量比对埋置基础的竖向振动有很大影响。     

Importance of footing–soil separation on dynamic stiffness of piled embedded footings

Different phenomena such as soil consolidation, erosion, and scour beneath an embedded footing supported on piles may lead to loss of contact between soil and the pile cap underside. The importance of this separation on the dynamic stiffness and damping of the foundation is assessed in this work. To this end, a numerical parametric analysis in the frequency domain is performed using a rigorous three‐dimensional elastodynamic boundary element–finite element coupling scheme. Dimensionless plots relating dynamic stiffness functions computed with and without separation effects are presented for different pile–soil configurations. Vertical, horizontal and rocking modes of oscillation are analyzed for a wide range of dimensionless frequencies. It is shown that the importance of separation is negligible for frequencies below those for which dynamic pile group effects start to become apparent. Redistribution of stiffness contributions between piles and footing is also addressed. Copyright © 2009 John Wiley & Sons, Ltd.     

Effective block diagonal preconditioners for Biot's consolidation equations in piled‐raft foundations

The finite element (FE) simulation of large‐scale soil–structure interaction problems (e.g. piled‐raft, tunnelling, and excavation) typically involves structural and geomaterials with significant differences in stiffness and permeability. The symmetric quasi‐minimal residual solver coupled with recently developed generalized Jacobi, modified symmetric successive over‐relaxation (MSSOR), or standard incomplete LU factorization (ILU) preconditioners can be ineffective for this class of problems. Inexact block diagonal preconditioners that are inexpensive approximations of the theoretical form are systematically evaluated for mitigating the coupled adverse effects because of such heterogeneous material properties (stiffness and permeability) and because of the percentage of the structural component in the system in piled‐raft foundations. Such mitigation led the proposed preconditioners to offer a significant saving in runtime (up to more than 10 times faster) in comparison with generalized Jacobi, modified symmetric successive over‐relaxation, and ILU preconditioners in simulating piled‐raft foundations. Copyright © 2012 John Wiley & Sons, Ltd.     

Finite element simulations of the behavior of piled‐raft foundations using a multiphase model

Finite element simulations of the behavior of a piled raft foundation have been carried out using a multiphase model conceived as an improved homogenization approach. According to this model, the ground reinforced by a group of piles is treated as a homogeneous continuous medium. In this approach, no specific interface elements are necessary to account for the mechanical interaction between the piles and the ground: this interaction is described by means of two scalar parameters, one stiffness parameter and one which can easily be derived from the maximum ground‐pile friction. The implementation of the model into a finite element code provides an efficient tool for the analysis of the influence of the pile number or length on the settlement and bearing capacity of a square piled raft foundation and of the way the total applied load is shared between the raft and the piles. Results are compared with a standard 3D finite element analysis. The comparison highlights the fact that the proposed approach remains to be improved to account for tip resistance. Copyright © 2012 John Wiley & Sons, Ltd.     

层状地基中瑞利波随深度的衰减特性

层状地基中瑞利波（Rayleigh wave）的传播规律以及随深度的衰减规律是工程环境微振动分析与控制研究的关键问题之一。采用解析法（矩阵传递法）分析了层状地基中瑞利波弥散曲线,并在此基础上研究了各种频率下不同类型土层（上软下硬、软夹层、硬夹层）中瑞利波位移随深度衰减规律,通过与弹性半空间解的对比发现,瑞利波位移随深度的衰减规律与频率、土层分布密切相关：频率越高,传播深度越浅;特殊土层（软夹层、硬夹层）的存在会局部显著地改变其衰减规律,并且位移峰值的分布受到频率和土层分布的共同影响。最后,将该方法用于分析上海光源工程地基条件下瑞利波随深度的衰减规律,并与现场深孔微振动测试数据进行了比较。结果证明,与弹性半空间解相比,层状地基中瑞利波随深度的衰减规律可以更为准确地反映真实条件下瑞利波随深度的衰减规律。     

Study on interaction between dissimilar piles in layered soils

This paper, under the assumption that the pile–pile interaction showed elasticity, analyzed the pile–pile interaction between two dissimilar piles in layered soils using the shear displacement method and taking the sheltering effect into account to modify the conventional interaction factor between two dissimilar piles. The linear stiffness of single pile was adopted herein to analyze the influence of the pile stiffness factor, the pile length factor, and the diameter of shorter pile on the pile–pile interaction factor between two dissimilar piles. The interaction factor between two dissimilar piles drew from the present method was generally in good agreement with the interaction factor computed by Wong and Poulos. However, this paper intends to provide a different but more convenient means of estimating the interaction factor between piles that have dissimilar diameters, lengths and pile tip conditions, and the settlement of dissimilar piles in pile groups. Copyright © 2010 John Wiley & Sons, Ltd.     

A coupled two‐phase fluid flow and elastoplastic deformation model for unsaturated soils: theory,implementation, and application

Although numerous numerical models have been proposed for simulating the coupled hydromechanical behaviors in unsaturated soils, few studies satisfactorily reproduced the soil–water–air three‐phase coupling processes. Particularly, the impacts of deformation dependence of water retention curve, bonding stress, and gas flow on the coupled processes were less examined within a coupled soil–water–air model. Based on our newly developed constitutive models (Hu et al., 2013, 2014, 2015) in which the soil–water–air couplings have been appropriately captured, this study develops a computer code named F²Mus3D to investigate the coupled processes with a focus on the above impacts. In the numerical implementation, the generalized‐α time integration scheme was adopted to solve the equations, and a return‐mapping implicit stress integration scheme was used to update the state variables. The numerical model was verified by two well‐designed laboratory tests and was applied for modeling the coupled elastoplastic deformation and two‐phase fluid flow processes in a homogenous soil slope induced by rainfall infiltration. The simulation results demonstrated that the numerical model well reproduces the initiation of a sheared zone at the toe of the slope and its propagation toward the crest as the rain infiltration proceeds, which manifests a typical mechanism for rainfall‐induced shallow landslides. The simulated plastic strain and deformation would be remarkably underestimated when the bonding stress and/or the deformation‐dependent nature of hydraulic properties are ignored in the coupled model. But on the contrary, the negligence of gas flow in the slope soil results in an overestimation of the rainfall‐induced deformation. Copyright © 2015 John Wiley & Sons, Ltd.