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

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

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

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

非等温分布条件下考虑半透水边界时饱和黏土的一维固结解析解

温度的变化会导致土体的物理力学性质改变,且在一些实际工程中,饱和黏土会处于非等温分布状态。为此,针对非等温分布条件下饱和黏土的一维固结问题,考虑了更具普遍性的半透水边界,通过某些假定推导了单级线性加荷形式下饱和黏土一维固结控制方程,并利用分离变量法求解得到了控制方程的解析解。通过将所提解析解分别与已有解析解和有限差分解展开对比分析,验证了所提解答的正确性。基于所提解析解,利用某一算例分析了温度梯度、半透水边界参数及加荷时间对固结性状的影响。结果表明：温度梯度 M 越大,土体的渗透性越大,土体的固结速率越快;半透水边界参数 R1和 R2越大,相同时间内土体的超孔隙水压力越小,土体的平均固结度越大;土体的平均固结度随加荷时间 tc 的增大而减小,这主要是由于加荷阶段所施加的外荷载小于最终荷载,但加荷时间tc的延长可一定程度减小土体中产生的最大超孔隙水压力。     

Stress history effects on 1‐D consolidation of soft soils: a rheological model

Stress history plays an important role in controlling the consolidation behavior of soft clays, but few models exist that can provide quantitative estimate of its influence. In this paper, the Gibson–Lo rheological model is used to simulate the coupled processes of drainage and creep of soft soils that takes stress history into account. A hybrid combination of analytical and numerical methods is adopted to solve the governing equations of consolidation with the nonlinear rheological model. The methodology is applied to a saturated soft soil subjected to surface loading. The soil profile is separated into normally consolidated and overconsolidated layers by a boundary that is allowed to move. Comparisons of the model predictions and its simulations are used to evaluate the effects of stress history, model parameters, and loading pattern on consolidation behavior. It is shown that stress history influences the location of the moving boundary, variations of the profiles of excess pore water pressure dissipation, stress and deformation‐based average degrees of consolidation. Parametric studies conducted show that when soil is stiffer, the excess pore water pressure dissipates much more quickly, and thus the soil consolidates much faster especially at the early stages. The results also show that soil viscosity influences the deformation‐based average degree of consolidation at the latter stages. The consolidation process of soil layer under linear loading is shown to lag behind those under instantaneous loading: the longer the loading period is, the smaller the average degrees of consolidation are no matter how they are defined. Copyright © 2012 John Wiley & Sons, Ltd.     

考虑起始坡降双层地基的大应变非线性固结

天然软土成层分布特性及土中渗流存在起始水力坡降的现象已被人们熟知。但变荷载下能同时考虑黏土中起始水力坡降、软土非线性压缩渗透特性及大应变特性的双层地基固结理论还鲜见报道。在拉格朗日坐标系中建立以超静孔压为变量的双层软土地基大应变非线性固结模型并给出其有限差分解。通过与考虑起始水力坡降的单层地基大应变非线性固结数值计算结果对比,验证了差分解的可靠性。着重分析了上、下土层起始坡降无量纲参数R1、R2对双层地基固结性状的影响,分析在大应变与小应变假定下双层地基超静孔压消散及固结沉降变形的异同。结果表明：上层土无量纲参数R1对双层地基固结性状的影响程度较下层土无量纲参数R2显著;大应变假定下双层地基渗流前锋的下移速度要快于小应变假定下的移动速度;大应变假定下考虑起始水力坡降的双层软土地基超静孔压消散速率要比小应变假定下快,且大应变假定下考虑起始水力坡降的双层地基最终沉降量要比小应变假定下大。     

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.     

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

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

起始水力梯度对饱和黏土一维固结的影响

采用考虑起始水力梯度的非Darcy渗流方程,修正了Terzaghi饱和黏土一维固结理论,给出了相应的有限差分格式,并据此探讨了起始水力梯度对渗流前锋、孔隙水压力和土层固结度的影响。计算结果表明,如果起始水力梯度大于0,离排水面较远土层的渗流就存在滞后现象,这样土层的固结要慢于Terzaghi固结理论。此时土层的最终平均固结度小于1,这可归结为土层中存在一定的残余孔隙水压力无法完全消散的缘故。     

变荷载下基于指数渗流双层地基一维固结分析

在土中渗流服从指数形式的前提下,建立了变荷载作用下双层地基的一维固结控制微分方程。利用有限差分法求得孔隙水压力的数值解,并通过与解析解对比对其可靠性进行了验证。对双层地基在指数形式渗流时不同参数下的固结性状进行分析,结果表明：单面排水条件下,双层地基中上层土渗流指数的大小对固结速率起决定性作用,而下层土渗流指数大小对固结速率的影响很小;如果上、下两层土体的压缩性不同,则地基按变形定义的平均固结度和按孔压定义的平均固结度不再相等;地基中下层土与上层土的相对压缩性越低、相对渗透性越高,则地基的固结速率越快;增大压缩性小、渗透性高的土层相对厚度,会加快双层地基的固结速率。     

基于黏弹性理论的沉桩再固结问题分析

假设桩周土体为饱和黏弹性介质,采用Burgers流变模型进行描述,同时考虑竖向和径向固结,建立了固结控制方程。根据不排水和自由排水情况,将边界条件分为3类并分别得到超孔隙水压力消散的级数解答,该解答能够为孔压静力触探反求固结系数提供一定的理论依据。在此基础上编制了应用程序,对Burgers流变模型中主要参数进行了分析。结果表明,地基表面自由排水、桩端地基不排水条件下,在一定深度以内的桩周土体的固结速度随深度降低,但超过某一范围后固结速度趋于稳定;上、下边界均自由排水条件下,固结速度随深度增加呈现下降、稳定、升高;上、下边界均不排水条件时,孔压消散速度不随深度变化,可简化为本解答仅考虑径向固结的特例。同时土体的流变特性对超孔隙水压力消散的影响比较显著,流变参数G1/1的变化使超孔隙水压力趋于某不为0的定值,且该值随G1/1比值的增大而增加;其他参数不变时,土体剪切刚度比G1/G2的增大会引起孔压消散速度的下降。     

Consolidation of sensitive clay with vertical drain

The coefficient of consolidation is one of the most important parameters that control the rate of consolidation. Conventional consolidation theories assume that the coefficient of consolidation is constant during the whole consolidation process. In the case of sensitive clay, the coefficient of consolidation is strongly dependent on the level of effective stress of clay. With the dissipation of pore water pressure and the increase of effective stress, the soil structure of the upper subsoil is gradually destroyed downwards and its coefficient of consolidation becomes smaller. At the same time, the coefficient of permeability of the vertical drains drops down because of the kinking or bending effect. The destructured upper subsoil and the kinking of the vertical drain hinder the dissipation of the pore pressure in the lower subsoil. This paper presents a model to describe the above important phenomena during the consolidation of sensitive clay with vertical drain. The solution for proposed model can be obtained by extending the closed‐form solution of the consolidation of double‐layered ground with vertical drain by the interactive method introducing the concept of the moving boundary and the reduction of discharge capacity of vertical drain. The numerical results obtained from the finite element method package PLAXIS (Ver. 7.2) are adopted to compare those obtained from the present algorithm. Moreover, the rationality of the moving boundary is explained by the distributions of the excess pore water pressure in natural soil zone along the radial direction. Wenzhou airport project is taken as a case study in this paper. The results for the sensitive soil with decaying sand drain agree very well with the in situ measured data. The calculated results can properly explain two general phenomena observed in the consolidation of soft sensitive soil ground with vertical drains: one is that the time to achieve the same consolidation degree is much longer under heavy loading than that under light loading; the other is that the consolidation speed is much slower in the lower subsoil than in the upper subsoil. Finally, it is indicated that the vertical drains can decrease the hindrance effect of the destructured subsoil. Copyright © 2007 John Wiley & Sons, Ltd.     

修正等应变假定下刚性桩复合地基固结特性分析

刚性桩复合地基固结特性是影响工程进度的重要因素。为研究外荷载作用下刚性桩复合地基固结特性，结合刚性桩和桩间土应力?应变特点对等应变假定进行修正。在此基础之上，考虑涂抹效应、桩间土和下卧层土体固结特性的影响，推导得到桩间土和下卧层土体固结微分方程。基于双层地基固结理论，考虑单级匀速施加荷载和附加应力沿深度变化等边界条件，推导得到其固结度解析解。最后，利用数值模拟方法对解析解进行验证，并对影响刚性桩复合地基固结特性的因素进行分析。结果表明：理论计算与数值模拟结果较为吻合；刚性桩复合地基固结速率受到贯入比影响较大，当贯入比较小时，与天然地基相比，其固结速率较小，贯入比较大且桩端阻力系数较小时，复合地基固结速率大于天然地基；随着贯入比、桩?土相对刚度和垫层?土体相对刚度的增大，复合地基整体固结速率加快；随着置换率的增大，复合地基整体固结速率降低。     

循环荷载作用下考虑非达西渗流的软黏土一维流变固结分析

基于Hansob渗流理论和Merchant流变模型,将微分型流变方程耦合到一维固结理论中,建立了循环荷载作用下双变量表述的软黏土一维流变固结耦合方程组。采用FlexPDE软件对方程组进行了求解,通过与已有文献结果对比,验证了算法的可靠性,在该基础上研究了循环荷载作用下饱和软黏土的一维流变固结特性,分析了Hansob渗流参数、Merchant流变模型参数以及循环荷载类型对饱和软黏土固结特性的影响。研究结果表明,按超孔压定义的平均固结度 和按变形定义的平均固结度 均呈现循环变化特征,且 > ;当循环周次较大时, 进行稳定循环状态,由于流变效应, 峰值仍呈逐渐增大趋势;随着Hansob渗流参数 和 的逐渐增大,地基固结速率逐渐降低,同一循环周次内 、 的峰值也随之降低;随着流变参数F的逐渐增大,同一循环周次内 的峰值逐渐降低, 的峰值显著提高;随着 的逐渐增大,同一循环周次内 、 的峰值逐渐降低。循环荷载作用下地基的平均固结度始终处于循环状态,无法达到固结完成状态,其最大平均固结度也远小于线性荷载的情况。     

基于非牛顿指数描述的非达西渗流一维固结分析

基于非牛顿指数描述的非达西渗流定律,同时考虑地基内部竖向附加应力随深度线性变化以及变荷载的影响,建立了一维固结控制方程并应用有限差分法进行数值求解,同时对不同参数单级加荷下的固结性状进行分析。结果表明：基于非达西渗流比达西渗流下固结速率要慢,且渗流模型中非牛顿指数越大,土层的固结速率越慢;土层厚度越厚,固结速率越慢,因此,传统固结理论中室内土样与地基土层之间的相似关系不再成立;作用于土层的平均附加应力越大,土层的固结速率越快;在单面排水情况下,附加应力分布对土层固结速率有较大影响;相反,双面排水条件下土层固结速率与附加应力的分布是无关的;荷载的加荷速率越快,则土层的固结速率越快。最后,讨论了达西渗流计算固结变形的适用范围。     

非均匀地基与均匀地基一维固结解答之间的转换关系

基于一维固结理论,研究了土层物理力学性质沿深度非均匀连续变化的固结问题。首先,利用分离变量法,获得了渗透系数和压缩系数随深度按指数函数变化的非均匀地基一维固结方程的解析表达式; 然后与Terzaghi固结理论的经典解答进行了比较,获得了两种解答之间的相似转换关系。这样,非均匀地基固结问题就可以用相同荷载以及边界条件下的均匀地基固结的经典解线性表示。因此,非均匀地基固结问题的求解转化为对相似转换系数的计算。该系数集中反映土层非均匀性对地基固结的影响,从而为解决非均匀地基一维固结问题提供了便捷途径。     

One-dimensional consolidation of soil under multistage load based on continuous drainage boundary

In engineering practice, a rapid loading rate can result in ground failure when the strength of soft soils is relatively low, and a multistage loading scheme is always utilized to deal with this situation. Firstly, under a multistage load and the continuous drainage boundary, an analytical solution of excess pore-water pressure and consolidation degree is obtained by virtue of the superposition formula of excess pore-water pressure, and a more general continuous drainage boundary under arbitrary time-dependent load is developed. Then, a comparison with existing analytical solutions is conducted to verify the present solution. A preliminary attempt on applying the continuous drainage boundary into the finite element model is made, and the feasibility of the numerical model for the one-dimensional consolidation under the continuous drainage boundary is verified by comparing the results calculated by FEM with that from present analytical solution. Finally, the consolidation behavior of soil is investigated in detail for different int erface parameters or loading scheme. The results show that, in land reclamation projects, a horizontal drain should be placed close to the boundary with a smaller interface parameter to improve the consolidation efficiency. The degree of consolidation is also related to the applied time-dependent load and interface parameters.     

复杂荷载作用下考虑下卧层三维渗流的未打穿竖井地基固结分析

针对任意复杂变荷载作用下未打穿的竖井地基,通过在竖井底面以下土层中设置虚拟竖井,使其能够合理考虑下卧层土体三维渗流问题,运用Laplace变换,求得频域内竖井地基的固结解。通过Laplace逆变换,得到了任意荷载作用下竖井地基的平均固结度、孔压消散曲线、沉降曲线。结合具体算例,对影响竖井地基固结的主要影响因素进行了详细分析,并将现有未打穿竖井地基平均固结度近似计算方法的精度和适用的范围进行了对比,得到了一些有益于工程实践的结论。     

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.     

A simple analytical solution to one‐dimensional consolidation for unsaturated soils

This paper presents a simple analytical solution to Fredlund and Hasan's one‐dimensional (1‐D) consolidation theory for unsaturated soils. The coefficients of permeability and volume change for unsaturated soils are assumed to remain constant throughout the consolidation process. The mathematical expression of the present solution is much simpler compared with the previous available solutions in the literature. Two new variables are introduced to transform the two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved with standard mathematical formulas. It is shown that the present analytical solution can be degenerated into that of Terzaghi consolidation for fully saturated condition. The analytical solutions to 1‐D consolidation of an unsaturated soil subjected to instantaneous loading, ramp loading, and exponential loading, for different drainage conditions and initial pore pressure conditions, are summarized in tables for ease of use by practical engineers. In the case studies, the analytical results show good agreement with the available analytical solution in the literature. The consolidation behaviors of unsaturated soils are investigated. The average degree of consolidation at different loading patterns and drainage conditions is presented. The pore‐water pressure isochrones for two different drainage conditions and three initial pore pressure distributions are presented and discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

考虑非Darcy渗流时循环荷载下饱和黏土一维固结分析

考虑到土体处于超固结状态下的压缩性一般要比正常固结状态下的低,引入描述非Darcy渗流的Hansbo方程,修正了Terzaghi饱和黏土一维固结方程,并将其推广到低频循环荷载作用的情况。采用有限体积法对该方程进行了求解,并讨论了非Darcy渗流参数、循环荷载周期以及超固结状态下土体压缩性对固结进程的影响。计算结果表明,在矩形波载作用下,按孔压定义和按变形定义的固结度都随时间增长而震荡增加,且按孔压定义的固结度的震荡幅度明显大于按变形定义的固结度。另外,非Darcy渗流减缓了地基的沉降速率,且循环荷载周期越短,或超固结状态下压缩性越高,地基的沉降速率就越慢。