Three-dimensional finite element analysis of spatially variable PVD improved ground

A stochastic approach that investigates the effects of soil spatial variability on stabilisation of soft clay via prefabricated vertical drains (PVDs) is presented and discussed. The approach integrates the local average subdivision of random field theory with the Monte Carlo finite element (FE) technique. A special feature of the current study is the investigation of impact of spatial variability of soil permeability and volume compressibility in the smear zone as compared to that of the undisturbed zone, in conjunction with uncoupled three-dimensional FE analysis. A sensitivity analysis is also performed to identify the random variable that has the major contribution to the uncertainty of the degree of consolidation achieved via PVDs. The results of this study indicate that the spatial variability of soil properties has a significant impact on soil consolidation by PVDs; however, the spatial variability of soil properties in the smear zone has a dominating impact on soil consolidation by PVDs over that of the undisturbed zone. It is also found that soil volume compressibility has insignificant contribution to the degree of consolidation estimated by uncoupled stochastic analysis.     

Analytical solutions for consolidation aided by vertical drains

Soil disturbance caused during the installation of vertical drains reduces the in situ hydraulic conductivity of soft deposits in the immediate vicinity of the drains, resulting in a slower rate of consolidation than would be expected in the absence of disturbance. Experimental investigations have revealed the existence of two distinct zones, a smear zone and a transition zone, within the disturbed zone around the vertical drain. The degree of change in the hydraulic conductivity in the smear and transition zones is difficult to assess without performing of laboratory tests. Based on the available literature, four different profiles of hydraulic conductivity versus distance from the vertical drain were identified. Closed-form solutions for the rate of consolidation for each of these four hydraulic conductivity profiles were developed. It is found that different variations of the hydraulic conductivity profiles in the disturbed zone result in different rates of consolidation.     

涂抹区重叠竖井地基固结特性研究

在软基处理工程中,经常出现竖井打设变密而地基固结效率降低的现象。鉴于此,建立了重叠涂抹区内土体水平向渗透系数的分布函数,给出了涂抹区重叠时竖井地基超静孔压和平均固结度的解析解。通过分析不同工况下竖井地基固结度随竖井间距的变化情况,探究了竖井间距减小而地基固结效率不增反减的成因。最后,探讨了涂抹作用和井阻作用对竖井最小临界间距的影响。结果表明:相邻竖井涂抹区重叠是竖井地基中出现竖井最小临界间距的根本原因。涂抹作用越大,则竖井最小临界间距越大;具体表现为当地基扰动程度增大时或涂抹区半径增大时,竖井最小临界间距随之增大。井阻作用越大,则竖井最小临界间距越小;具体表现为当竖井渗透系数减小时、竖井长度增大时或竖井半径减小时,竖井最小临界间距随之减小。     

真空联合堆载预压下混凝土芯砂石桩复合地基固结特性理论分析

真空联合堆载预压与混凝土芯砂石桩复合地基相结合是处理深厚软土地基的一种新方法。该方法既能利用预制混凝土芯桩提高地基承载力,又能利用砂石外壳缩短排水路径、传递真空负压,加快固结。根据真空联合堆载预压下混凝土芯砂石桩复合地基的固结特点,将堆载预压和真空预压的固结效应分开考虑,然后再进行叠加,推导出了真空联合堆载预压下混凝土芯砂石桩复合地基的平均固结度解析解,并利用三维数值模拟对平均固结度解析解进行了验证。通过对解析解的参数分析,探讨了置换率、涂抹区大小及渗透系数、混凝土芯砂石桩长径比和芯桩率对混凝土芯砂石桩复合地基固结特性的影响。结果表明,混凝土芯砂石桩复合地基的固结速率随着置换率和芯桩率的增大而增大,随涂抹区和砂石外壳直径之比、未扰动区和涂抹区渗透系数之比、长径比的增大而减小。     

基于非等温管道流竖井地基热排水固结模拟

针对新技术热排水固结法,采用非等温管道流模拟竖井中U型导热管的传热过程,考虑温度对竖井扰动区和未扰动区渗透性的影响,在COMSOL Multiphysics有限元软件中进行二次开发,建立了竖井地基热排水固结法的有限元模型。以热排水固结软基处理原型试验为例,重点分析了模型耦合、部分耦合和不耦合情况下软土地基的固结度。结果表明,相对于传统排水固结法的不耦合模型,部分耦合模型下因温度产生的孔压增量延缓了地基固结的发展,固结速率有所减慢;耦合模型下温度虽也产生一定的孔压增量,但温度有效地改善了竖井涂抹区土的渗透特性,地基的固结速率加快,固结周期缩短,与试验结果一致。     

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

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

考虑温度耦合效应的竖井地基固结有限元分析

针对塑料排水板(PVD)安装热源能提升PVD性能、加速竖井地基固结这一工程现象,基于热-水-应力 (T-H-M) 三场全耦合的有限元方法来模拟利用热源进行地基处理新技术(PVTD)。首先,以微分形式与等效弱形式分别给出T-H-M耦合控制方程,并推导出其有限元方程组。然后在多场耦合有限元软件中建立饱和土的T-H-M全耦合模型,并通过与已有解析解比较,验证了模型正确性。最后,对一个经典有涂抹区的竖井地基算例,分不耦合温度(UT)、耦合温度但不考虑其对饱和土物性影响(CT)、耦合温度考虑温度对饱和土渗透性影响(CTP) 3种情况进行固结计算分析。研究结果表明,相对于无热源竖井地基,CT情况下由于热源产生的附加孔隙水压力,固结速度略有下降;CTP情况下,由于热源有效改善涂抹区的渗透性能,竖井地基固结速率明显加快。上述研究结论从理论上较好地阐明了PVTD的作用机制。     

考虑淤堵效应的疏浚淤泥真空固结沉降计算

疏浚淤泥含水率高、强度极低,真空梯度作用下排水板周围会快速形成致密"土柱",导致排水固结能力迅速下降,即出现淤堵现象。由疏浚淤泥真空固结室内模型试验获得了出水量变化情况和土体含水率、渗透系数的时空分布情况,确定了排水体周围淤堵区的形成时间与范围,提出了考虑时间效应与淤堵效应的真空度传递模式,同时考虑疏浚淤泥土体压缩和渗透的非线性,建立了考虑淤堵效应的固结分析模型,获得了相应的解析解,通过与已有数据和现场试验结果对比验证了该解答的有效性。利用该解析解,分析了真空度传递模式、淤堵系数λ和淤堵比c对沉降的影响。结果表明,淤堵效应明显减缓了沉降速率,也严重降低了疏浚淤泥的最终沉降量。     

考虑3种非线性关系的径向排水固结解析解

大部分砂井及排水板固结解析解理论结合了太沙基一维固结理论, 把土体考虑成线弹性问题, 而研究发现, 软土在等向固结条件下平均有效应力和水平渗透系数与孔隙比成半对数线性关系。同时, 实验表明, 在塑料排水板周围的涂抹区中水平渗透系数成抛物线状分布。因此, 为了充分考虑这些非线性的影响, 推导出了该问题的解析解, 并通过退化验证了该解析解的正确性。同时, 还对参数Cc/Ck和对固结速率的影响进行了分析, 结果表明:当Cc/Ck1时, 固结速率较快; 当Cc/Ck1时, 固结速率较慢; 参数增加, 固结速率逐渐变小。     

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

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

考虑井阻及涂抹作用的非饱和土竖井地基固结解析解

以往的非饱和土竖井地基研究中未同时考虑竖井的井阻和涂抹作用,大部分按理想竖井进行研究,然而井阻和涂抹作用是影响非饱和土竖井地基固结的重要因素.针对这种情况,本文基于Fredlund非饱和土一维固结理论及等应变假设,引入变量将超孔隙压力耦合控制方程组转化为等价的线性偏微分方程组,考虑涂抹和井阻条件,并采用分离变量法和待定...     

Analytical solution to axisymmetric consolidation of unsaturated soil stratum under equal strain condition incorporating smear effects

This paper proposes closed‐form analytical solutions to the axisymmetric consolidation of an unsaturated soil stratum using the equal strain hypothesis. Following the 1‐dimensional (1D) consolidation theory for unsaturated soil mechanics, polar governing equations describing the air and water flows are first presented on the basis of Fick's law and Darcy's law, respectively. The current study takes into account the peripheral smear caused by an installation of vertical drain. Separation of variables and Laplace transformation are mainly adopted in the analytical derivation to obtain final solutions. Then, the hydraulic conductivity ratio, the radius of influence zone and smear parameters influencing time‐dependent excess pore pressures, and the average degree of consolidation are graphically interpreted. In this study, a comparison made between the proposed equal strain results and the existing free strain results suggests that both hypotheses would deliver similar predictions. Moreover, it is found that the smear zone resulting from vertical drain installations would hinder the consolidation rate considerably.     

土−水特征曲线测量过程中孔隙分布的演化规律探讨

测量了在宽广吸力范围内原状样和压实样的脱湿持水曲线,对比分析了单双峰结构持水性能的差异;并利用压汞试验测试两种土样在脱湿过程的孔隙分布,分析了两者的差异并探讨了脱湿过程孔隙的演化规律;在考虑收缩变形的基础上,基于孔隙分布曲线确定了土?水特征曲线的基本参数。试验结果表明：原状样在宽广吸力范围内基本上呈单峰孔隙结构;饱和压实样具有单峰孔隙结构,随着吸力的增加,双峰结构越来越明显,当吸力达到很大时,演化成完全双峰孔隙结构。原状样的持水曲线为经典的S形,而压实样的持水曲线在过渡段出现了水平台阶状;低吸力段,压实样的持水曲线低于原状样,而高吸力段,两者的持水曲线基本重合。基于孔隙分布曲线确定了控制持水曲线进气值和残余值的孔径,并计算出对应的吸力值,其值更符合实际物理意义。     

考虑涂抹区压缩及渗透系数变化的堆载预压固结解

为评估涂抹区土体压缩和渗透系数变化对含竖向排水体地基固结的影响,采用等体积应变假设,考虑涂抹区土体的压缩变形及其水平向渗透系数沿径向分别呈线性和抛物线分布,并考虑井阻作用以及地基附加球应力沿深度任意分布,推导了随时间线性堆载预压条件下固结微分方程的显式解析解答,分析了涂抹区半径、水平向渗透系数的分布模式、以及体积压缩系数对地基整体平均固结度的影响。结果表明,涂抹区土体采用均匀折减的水平向渗透系数明显低估了地基的固结速率,而当涂抹区半径较大时,不考虑涂抹区土体的压缩变形将会高估地基的固结速率。在含竖向排水体地基固结问题的分析中,这些影响不可忽视。     

考虑加载过程及桩体固结变形的碎石桩复合地基固结解析解

为了完善碎石桩复合地基固结理论,通过假设从桩体排出的水量等于流入桩体的水量与桩体体积变化之和以及地基扰动区土体水平渗透系数呈线性变化,并考虑上部荷载逐渐施加,推导了考虑桩体体积变化的碎石桩复合地基超静孔压及固结度解析解。当加载时间趋于零时,本文解可退化为瞬时加载情况下的解;当加载时间及桩径同时趋于零时,本文解可进一步退化为Terzaghi一维固结解,这证明了本文解的正确性。通过与已有解的比较,对地基固结性状进行了分析。结果表明,加载过程对地基固结度影响显著,加载历时越长,固结越慢;在各种条件下,不考虑桩体固结变形时地基固结始终比考虑桩体变形时快,并且其影响随着加载历时变小、桩径比变小、桩土模量比变小、桩土渗透系数变小而逐渐增大,这说明在实际工程固结计算中不考虑桩体固结变形是偏于不安全的。     

A general theoretical solution for the consolidation of a composite foundation

A theoretical study is made of the consolidation of a composite foundation reinforced by columns. A more reasonable initial condition for the consolidation of a composite foundation is derived from the equilibrium equation and the equal strain assumption. Then, the general analytical solutions for excess pore water pressures in the column and in the surrounding soil are obtained using this new initial condition under the influence of several important factors. These factors include: the gradual variation of horizontal permeability in the disturbed soil, the change in the total average stress in a composite foundation with depth, and the effect of time-dependent external loads. The general expressions for the overall average degree of consolidation of a composite foundation, both in terms of stress and in terms of deformation, are then derived. On the basis of the above theory, the explicit solutions for a particular case, which considers ramp loading and a linear change of total stress with depth, are given in detailed forms. Finally, a parametric study and a comparison of some available solutions are made. The results show that the average degree of consolidation obtained in terms of stress and in terms of deformation is the same; the increase in construction time, the disturbance intensity during column construction and the size of the disturbed zone reduce the consolidation rate of a composite foundation. For the case of PTIB (pervious top surface and impervious bottom surface), the increase in the ratio of top total stress to bottom total stress in the composite foundation accelerates consolidation.     

吹填泥浆中土颗粒沉降-固结规律研究

吹填泥浆中土颗粒沉降-固结过程十分复杂,其中土颗粒的沉降固结规律并不明确。室内沉降试验采用一种改进的方法--分层抽取法,先得到土颗粒体积分布规律,再利用体积通量函数法计算不同粒组颗粒平均沉速。计算结果表明,吹填泥浆在沉降过程中会出现3个区域,分别是干涉沉降区、絮体压缩沉降区和自重固结区;在沉降初期0～1 h和1～7 h内,干涉沉降区内土颗粒沉速远小于Stokes公式沉速计算值,颗粒沉速分别与颗粒直径的0.488 1和0.111 7次方有关。絮体压缩沉降区和自重固结区内颗粒基本不出现分选现象,均属于泥浆颗粒的沉积物;该试验中絮体压缩沉降区的密度约为1.04 g/cm3。     

Consolidation of vertical drain with depth-varying stress induced by multi-stage loading

An explicit analytical solution is developed for the consolidation of vertical drain with both radial and vertical drainage by adopting a depth-varying stress induced by multi-stage loading. The well resistance and smear effect are also considered. The smear effect is described by three decay patterns of horizontal permeability towards drains within the smeared zone, including a reduced constant pattern, a linear decay pattern and a parabolic decay pattern. A parameter analysis is performed to investigate the consolidation behavior of the vertical drain. The convergence of the proposed series solution is discussed.     

卸荷损伤原状膨胀土剪切力学特性试验研究

通过GDS三轴试验系统对经历3种卸荷速率损伤后的原状南阳膨胀土样进行再加荷不排水三轴伸长剪切试验,同时考虑了超固结比与固结状态的影响。试验结果表明,膨胀土再加荷剪切力学特性与初始卸荷速率有关。在相同的轴向应变下,初始卸荷速率越小,其偏应力单调越小。在主应力方向改变前后,其应力?应变关系曲线斜率显著变化。相同固结方式与超固结比状态下,孔隙水压力均表现为先增大后减小趋势,孔隙水压力峰值应变随卸荷速率增大而减小。无论是等压固结还是K0固结,初始卸荷速率越大,不排水剪切强度越大。膨胀土样经历了初始卸荷损伤后,再加荷常规三轴伸长试验所得剪切强度均低于无损伤时的强度。以膨胀土破坏强度所得损伤度SD低估了卸荷速率对膨胀土的损伤程度,建议采用孔隙水压力峰值强度进行膨胀土边坡设计计算。原状膨胀土力学性状随卸荷速率损伤的演化规律受卸荷阶段轴向应变大小及裂隙性综合影响。     

Trial Embankment Analysis to Predict Smear Zone Characteristics Induced by Prefabricated Vertical Drain Installation

In this study, FLAC finite difference software has been adopted to simulate the performance of the ground improved using prefabricated vertical drains assisted preloading, considering smear zone characteristics. The numerical code has been applied to predict smear zone properties employing a back calculation procedure using the results of several case studies. The construction of a trial embankment is proposed as a reliable method to predict the smear zone characteristics. The proposed back calculation method is applied to estimate the minimum required degree of consolidation and consequently the minimum required preloading time, resulting in a reliable estimation of the smear zone permeability and extent. Three preloading case studies considering both conventional preloading and vacuum assisted preloading have been simulated to verify the numerical code and to conduct the parametric study using the back calculation procedure. According to the results, the properties of the smear zone can be back-calculated reliably, when at least 33 % degree of consolidation due to trial embankment construction is achieved.