盾构施工引起地表固结沉降问题的研究

将隧道周围土体视为均质连续各向同性的饱和弹性介质,采用保角变换的方法将含有隧道的半无限平面映射为同心圆环计算域。根据Terzaghi-Rendulic二维固结理论,建立隧道在不透水的情况下周围土体超孔隙水压力分布的控制方程。然后,采用分离变量法计算得到土体超孔隙水压力分布的解析解,最后,根据弹性理论计算得出隧道中线上方地表固结沉降的计算公式。结合算例,分析盾构施工扰动程度、土体渗透系数、土体弹性模量及隧道埋深等因素对隧道中心上方地表处固结沉降的影响。研究结果表明,地表固结沉降的增加值与隧道外侧初始超孔隙水压力值C0的变化量成正比例关系,施工扰动程度越大所引起的固结沉降越大;土体的渗透系数越大固结沉降速度越快,但土体的渗透系数与最终的地表固结沉降量无关;土体的弹性模量越大,最终的地表固结沉降量越小;隧道埋深越深,地表固结沉降所需时间越长,最终的地表固结沉降量也越大。     

考虑涂抹区渗透系数变化的非饱和土砂井地基固结特性分析

在以往对非饱和土砂井地基固结理论研究中,均将涂抹区与非涂抹区土体渗透系数视为相等,这与实际工程并不相符。本文将考虑涂抹区土体渗透系数的变化,分析其对超孔隙气、水压力消散规律的影响。基于Fredlund一维固结理论以及Darcy定律和Fick定律,对有限厚度线弹性非饱和土砂井地基,在大面积均布瞬时荷载作用下,考虑涂抹区土体渗透系数的变化,利用Laplace变换并引入Bessel函数推导出Laplace变换下的解,再通过Crump方法编程实现Laplace逆变换得到超孔隙气压力、超孔隙水压力的半解析解。利用典型算例进行计算,分别得到在不同半径、不同涂抹区半径和不同涂抹程度的情况下,超孔隙气压力、超孔隙水压力随时间的变化规律。得出考虑涂抹作用时,超孔隙气、水压力的消散速度降低;涂抹区半径越大、涂抹程度越高速度越慢,反之消散越快。本研究丰富了非饱和土砂井固结理论,对非饱和土砂井固结特性的研究具有一定的工程参考价值。     

基于一维和三维固结挤土桩沉降时效计算分析

饱和软黏土地基中的挤土桩,施工引起的桩侧土体超孔隙水压力消散和竖向荷载引起的桩底土体固结是导致其沉降的主要因素。基于三维固结理论,研究了桩侧地基土再固结沉降的变化规律,并从桩-土相互作用的原理出发,结合边界条件,推导出了在桩侧地基土再固结沉降的作用下挤土桩沉降的计算公式。运用太沙基一维固结理论并采用分层总和法计算竖向荷载作用下桩底土体的沉降。提出了饱和软土中挤土桩沉降时效的计算方法。并对工程实例进行了计算分析,其结果符合工程实际沉降规律。     

Foundation settlement statistics via finite element analysis

The dispersion observed in soil data comes both from the spatial variability which greatly influences the behavior of large structures and from errors in testing. Thus, the geotechnical engineering deals with uncertainties for which deterministic approaches are not suitable. The resort to probabilistic techniques, enables modeling uncertainties by analyzing their dispersion effect on the global behavior of the structure. The scope of this paper is analyzing settlement and differential settlement variability of a pair of foundations on random heterogeneous medium. The random soil properties of interest are the elastic modulus, and the Poisson ratio. The elastic modulus is modeled as a spatially random field by adopting the lognormal distribution, which enables analyzing its large variability. Because soil Poisson ratio is bounded in practice between two extreme values, its random field is obtained by using the Beta distribution. In this study, one proposes for the Beta field determination, a mapping technique on the probability distribution function diagram, by solving a non-linear equation. However, the mean and variance are unchanged through the mapping operation. Because the soil Poisson ratio is a positive parameter, one prefers to perform the mapping operation with the probability function of the lognormal distribution. Also, the proposed technique can be used for other bounded soil properties such as the porosity. In this paper, settlement and differential settlement statistics prediction are carried out using Monte Carlo simulations combined with deterministic finite element method (DFEM). A performed parametric study shows the following: (i) as the variability of the elastic modulus increases as settlement and differential settlement statistics are important, also, settlement statistics decreases as the Poisson ratio variability increases, and differential settlement statistics do not seem be affected by its variability. (ii) settlement and differential settlement statistics are important for positive inter-property correlation. (iii) a great influence of the correlation lengths on settlement and differential settlement statistics.     

A closed form analytical solution for two‐dimensional plane strain consolidation of unsaturated soil stratum

This paper discusses the excess pore‐air and pore‐water pressure dissipations and the average degree of consolidation in the 2D plane strain consolidation of an unsaturated soil stratum using eigenfunction expansion and Laplace transformation techniques. In this study, the application of a constant external loading on a soil surface is assumed to immediately generate uniformly or linearly distributed initial excess pore pressures. The general solutions consisting of eigenfunctions and eigenvalues are first proposed. The Laplace transform is then applied to convert the time variable t in partial differential equations into the Laplace complex argument s. Once the domain is obtained, a simplified set of equations with variable s can be achieved. The final analytical solutions can be computed by taking a Laplace inverse. The proposed equations predict the two‐dimensional consolidation behaviour of an unsaturated soil stratum capturing the uniformly and linearly distributed initial excess pore pressures. This study investigates the effects of isotropic and anisotropic permeability conditions on variations of excess pore pressures and the average degree of consolidation. Additionally, isochrones of excess pore pressures along vertical and horizontal directions are presented. It is found that the initial distribution of pore pressures, varying with depth, results in considerable effects on the pore‐water pressure dissipation rate whilst it has insignificant effects on the excess pore‐air pressure dissipation rate. Copyright © 2015 John Wiley & Sons, Ltd.     

盾构施工引起的固结沉降分析

盾构在低渗透性土层中开挖,常常伴随着地表下沉,究其原因为超孔隙水压力消散引起固结沉降的结果。为对盾构施工引起土体的固结沉降进行研究,首先,根据隧道施工前后土体应力的变化值,应用Henkel超孔隙水压力理论,推导了隧道开挖引起的初始超孔隙水压力的计算公式,并采用数值分析方法,考虑了由于土体的固结引起的沉降变形。研究成果应用到上海地铁2号线,根据具体的地质条件进行理论计算分析。结果表明,隧道开挖引起的初始超孔隙水压力最大值位于隧道起拱线处,地表固结沉降预测值与实测值吻合较好。     

Prediction of Probabilistic Settlements by the Perturbation-Based Spectral Stochastic Meshless Local Petrov–Galerkin Method

Originating an attempt of understanding the reliability and serviceability of foundations, an interest of comparing the difference between settlements predicted with and without considering the uncertainty in such as the spatial variability of soil properties is born. This study selectively compares between settlements predicted with and without considering the uncertainty in the spatial variability of Young’s modulus. The tool is a coupling of perturbation expansions of Young moduli and a two-dimensional meshfree weak-strong form in elastostatics. Two further examples show that the spatial variability of Young’s modulus causes apparent difference between probabilistic and deterministic settlement components along the direction of a surcharge. We can derive an autocorrelation function to describe the spatial variability of Young’s modulus and understand how it affects predicted settlements depending upon autocorrelation function values. In addition, the spectral stochastic meshless local Petrov–Galerkin method is a time-saving tool for predicting probabilistic settlements with the uncertainty in the spatial variability of soil properties.     

循环荷载下双曲线模型修正土体一维固结解

假定土体的应力-应变关系满足双曲线模型、渗透系数的降低与压缩系数的减小成正比以及初始有效应力沿深度均匀分布,推导出任意荷载下单层地基的非线性一维固结方程。结合具体的定解条件采用分离变量法求出土体有效应力的通解,并进一步求出超静孔隙水压力、按沉降定义的平均固结度Us以及按孔隙水压力定义的平均固结度Up的通解。根据通解求出梯形低频循环荷载作用下的解,并找出循环荷载作用下土体固结性状的影响因素。基于循环荷载作用下的解,采用FORTRAN语言编制了计算程序,计算并分析考虑土体本构关系非线性时的固结性状,以及反映土体本构特征的参数E和n对固结性状的影响。计算结果表明：①Us稍大于Up;②距排水平面越近,附加有效应力受加载方式的影响越明显,距排水平面越远则相反;③E和n对Us的影响很小,可以忽略不计;④E和n对Up和附加有效应力的影响较大且不容忽略。     

考虑幂函数渗透压缩变化的双层砂井地基固结性状

采用Barron轴对称固结及大变形固结问题的某些简化与假定,推导建立了砂井地基大变形固结控制方程,利用建立的双层砂井地基大变形固结方程及编制的计算程序,通过引入软土渗透系数、有效应力与孔隙比之间的幂函数关系k =ced与e=a( )b,对瞬时加载下双层砂井地基固结性状进行算例计算。结果表明：（1）双层软土幂函数渗透关系及压缩关系中诸参数对双层砂井地基固结性状有重要影响：随着两层软土幂函数渗透关系中参数c1、c2的增加（渗透性增加）、或幂函数压缩关系中参数a1、a2的增加,各土层水平径向与竖向孔隙比减小更快,沉降发展速率与超静孔压消散速率也相应增加,且沉降发展速率快于孔压消散速率。（2）两层土在分界面处的孔隙比及平均超静孔压均出现明显的突变,将沿深度分布曲线分成形状不同的两段,表现出不同的固结性状。     

Semi‐analytical solution to one‐dimensional consolidation for unsaturated soils with semi‐permeable drainage boundary under time‐dependent loading

This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

非饱和土一维大变形固结模型

基于分段线性差分法,建立了一种非饱和土一维大变形固结模型。该模型可考虑土性参数非线性变化,可计算与分析大变形问题,并编制了Fortran计算程序。在现有解答和试验数据的基础上,对该模型进行了验证,瞬时加载情况下模型数值解与现有解答基本吻合,考虑加载过程下的数值解与试验数据吻合。进行了大变形算例分析,对比了加荷压密与消散固结阶段土层变形,探讨了孔隙气、水渗透系数比对土层沉降量、饱和度和不同应变情况下固结度的影响规律,分析了非饱和土大、小变形固结理论计算孔隙水（气）压和沉降量的差异。     

强夯荷载作用瞬间饱和土层固结变形计算

在已有强夯荷载作用下饱和土层内孔隙水压力分布理论研究成果的基础上,提出了一个预估强夯荷载作用瞬间由于孔隙水压力消散所引起的饱和土层固结变形量的计算方法;讨论了固结系数、卸荷固结比等土性参数对强夯冲击荷载作用下固结变形量大小的影响,进一步解释了强夯法加固饱和土层的机理。     

不同温度应力路径下饱和黏土剪切特性

在精准温控动三轴试验系统上开展了不同温度及不同升温路径饱和黏土剪切试验研究,探讨了不同温度对饱和软黏土不排水剪切特性的影响,分析不同升温固结方式对饱和软黏土孔压发展、体变、强度以及模量的影响规律。试验结果显示：在4～76 ℃试验研究范围内,环境温度升高导致饱和软黏土的不排水剪切强度有所减少,但温度升高对土体模量增加影响明显,温度T和模量ET关系可用ET = 2.69T 0.3表达;升温变化时正常固结黏土产生超孔隙水压力并随着温度增大而增大,升温热固结后土的剪切强度将明显提高,且排水状态下升温固结对土剪切强度增长小于升温完成后再固结情况;土体从26 ℃分别升高20、40 ℃时,升温引起的超孔压比分别为0.41、0.61,剪切峰值强度分别增加8.23%、22.37%。研究表明：升温幅值增大会使土体热固结程度越大,升温分级越多,热固结也越充分,其对应的体变、强度增长率则越大;同时最终温度及热固结路径对其剪切相转换特征存在影响,升温越高、热固结路径越多其剪胀性越明显,但温度变化范围、固结分级、热固结路径总体上对孔隙水压力的发展基本不产生影响。     

饱和软土成层地基一维非线性固结解析解

假定土体在固结过程中压缩性和渗透性的变化成正比,基于 - 及 - 关系,推导出饱和软土成层地基一维非线性固结解析解,分别给出了按沉降定义和按有效应力定义的每层土平均固结度及整个土层总固结度的计算公式。采用Fortran语言编制了相应的计算程序,将计算得到的结果与已有双层地基一维非线性固结解析解计算结果进行比较,验证该解析解的正确性。利用该程序分析成层地基一维非线性固结性状,分别讨论了初始竖向渗透系数、初始体积压缩系数、荷载值及土层厚度对地基固结性状的影响。分析结果表明：在成层地基一维非线性固结过程中,初始竖向渗透系数对超静孔压的影响较为复杂,对上下层地基固结速率影响不同;初始体积压缩系数增大,超静孔压增大,固结速率变小;所加荷载值越大,超静孔压消散越慢,固结发展越慢;超静孔压消散速率不仅取决于土层厚度,同时取决于各层土渗透性的相对大小。     

竖井地基的大应变固结分析

简单介绍了真空-堆载联合预压下饱和软土竖井地基大变形固结沉降模型（VRCS1模型）,该模型能考虑大应变、材料非线性、非达西流、真空折损和竖井未打穿等因素。使用该模型对澳大利亚巴利纳支路某现场沉降进行预测,结果与现场实测数据吻合良好。基于该工程案例,讨论了分级堆载、循环荷载、真空联合堆载、上边界为等应力/等应变条件以及排水板打入深度等因素对竖井地基固结过程的影响规律。计算结果表明：分级堆载可降低土体超静孔压峰值进而改善土体稳定性;循环荷载使土体固结过程发生震荡,且土体沉降峰值迟于超静孔压峰值和荷载峰值出现;真空荷载和堆载作用机制有本质区别,真空荷载不可简单以堆载替代;上边界等应变条件下的固结速率一般快于等应力条件,工程真实条件一般介于两者之间;排水板打入深度超过0.9倍土层厚度时,再加大打入深度加速固结效果不明显。     

考虑渗流与变形耦合作用的基坑工程空间效应

基于三维比奥固结理论,编制了考虑地下水渗流与土骨架变形耦合效应的有限元程序,分析了基坑变形以及坑内外超静孔压和土水势的空间分布规律,并与二维分析的结果进行了对比。研究表明,由于空间效应的影响,基坑中部截面处的围护结构水平位移和坑后地表沉降较基坑边缘截面处的围护结构水平位移和坑后地表沉降大,而基坑中部截面处的超静孔压和土水势却较小;与三维有限元分析的结果相比,二维分析时围护结构的水平位移和坑后地表沉降都更大。     

海相淤泥软土地基强夯置换砂桩试验分析

针对某海相淤泥软土的处理工程,选取了一强夯砂桩置换试验区,监测了强夯置换过程中夯坑及其周围土体的变形,超孔隙水压力的增长和消散情况以及深层土体水平位移。分析结果表明,单个置换砂桩夯击8次较为合理;强夯置换后静置一段时间,可以有效减小强夯过程中的超静孔隙水压力,且一遍强夯后静置7 d后超孔隙水压力基本消散完;深层土体的水平位移主要发生在地表以下3～6.5 m深度范围内。     

Study of the 1-D consolidation behavior of two-layered soft soils: Parametric studies using a rheological model with viscoplastic body

The governing equations for the coupled processes of consolidation and creep of two-layered soft soils are established. The Nishihara rheological model is adopted to simulate the elasto-viscoplastic characteristics of soft soils, disregarding the effects of the soil self-gravity. A semi-analytical theory combined with numerical and analytical methods is introduced to solve the governing equations of the one-dimensional rheological model. The computational procedure and the approximate solutions for two-layered soft soils subjected to surface loading are obtained for two drainage conditions. The solutions and the computational procedure are used to study the effects of the two layers and constitutive parameters on rheological consolidation behavior of soft soils. It can be concluded that two layers affect the rate of excess pore water pressure dissipation and settlement development. The parametric studies show that when the parameters of the upper layer remain constant, increases in the permeability and elastic modulus in the lower layer accelerate the dissipation of the excess pore water pressure, and meanwhile increases in the viscosity coefficient and viscoplastic limit slows down the dissipation of the excess water pressure.     

考虑变渗透系数的饱和黏土一维流变固结分析

为进一步探讨饱和黏性土的弹黏塑性和变渗透性对一维固结进程的影响,引入考虑时间效应的统一硬化（UH）本构模型描述饱和黏土的弹黏塑性,并用Taylor的经验关系式表示其渗透系数与孔隙比的关系,修正了太沙基一维固结理论,并利用有限差分法对方程进行求解。通过与一维流变固结试验结果对比,验证本文计算方法的有效性以及UH模型的适用性。在此基础上,讨论了UH模型参数、压缩指数与渗透指数的比值以及荷载强度等对固结过程的影响。计算结果表明:饱和黏土的黏滞性使得在固结初期地基内部出现了孔压升高现象,并减缓了固结中后期地基孔压的整体消散,增大了地基的沉降量,并且这种现象随着初始超固结参数的增大而更明显;另外,随着压缩指数与渗透指数比值的减小,地基固结进程会加快,但地基最终沉降量不受影响;同时,在固结后期,当压缩指数不大于渗透指数时,较小荷载强度对应的地基内孔压消散较慢,然而,当压缩指数大于渗透指数时,荷载强度对孔压消散的影响则相对很小。     

盾构法隧道施工工后横向地表总沉降研究

盾构法隧道施工引起的土体超孔隙水压力消散,会导致工后地表继续沉降。在前人研究基础上,提出了隧道周围土体的初始超孔隙水压力计算方法。采用分层总和法,计算由于扰动区内土体超孔隙水压力消散引起的地表固结沉降量,结合施工阶段地表沉降量,叠加得到总的工后横向地表沉降。提出了固结开始t时刻的地表总沉降量计算方法,研究了地表沉降速率随时间的变化。算例分析结果表明：本文方法计算得到的横向地表沉降曲线符合正态分布规律,与实测值吻合;地表总沉降速率的发展也与实测沉降速率曲线相符。