沉积泥砂非线性大变形固结沉降计算模型

沉积形成的水底黏性泥砂自重固结过程表现出显著非线性大变形固结特征,应采用大变形固结理论进行泥砂沉积固结计算。基于软黏土一维非线性大应变固结理论,应用有效应力、渗透系数与孔隙比间扩展幂次函数固结本构关系,由达西定律、有效应力原理、连续介质方程等建立大变形固结控制方程,根据固结单元孔隙水渗流、单元变形与泥砂沉积层固结沉降耦合关系形成黏性泥砂大变形自重固结数值模型。泥砂自重作为固结荷载,数值模型假定沉积泥砂各向同性且固结沉降应变、孔隙水渗流仅发生于竖直方向,为一维单向沉积固结过程;采用泥砂沉降柱试验确定泥砂非线性扩展幂次函数关系参数。模型应用中,划分竖向固结单元,由沉积泥砂固结本构关系确定各固结单元有效应力及超孔隙水应力,通过超孔隙水应力时间维度上的消散过程及各固结参数间的耦合关系计算泥砂固结沉降。数值模型计算结果表明,沉积黏性泥砂自重固结初期表现为有效应力调整过程,初始有效应力与孔隙比根据固结本构关系匹配调整为扩展幂次函数关系;沉积泥砂应变与应力固结度存在20%左右误差,泥砂固结沉降发展快于超孔隙水应力消散过程,证明沉积泥砂固结沉降变形的发展与超孔隙水应力消散并非同步耦合。计算模型应用于室内沉降柱试验模拟淤积黏性泥砂自重固结沉降预测中,模型输出与试验结果符合良好。     

Model for large strain consolidation under constant rate of strain

A numerical model, called CCRS1, is presented for one‐dimensional large strain consolidation under constant rate of strain loading conditions. The algorithm accounts for vertical strain, general constitutive relationships, relative velocity of fluid and solid phases, changing compressibility and hydraulic conductivity during consolidation, and an externally applied hydraulic gradient acting across the specimen. Soil compressibility is rate independent, and as such, the current model is most appropriate for less‐structured clays. Verification checks show excellent agreement with analytical and numerical solutions for small and large strain conditions. A series of numeric examples indicates that compressibility and hydraulic conductivity constitutive relationships can have an important effect on constant rate of strain consolidation response. Results also indicate that analytical solutions obtained using small strain theory can be in significant error for large strain conditions with changing coefficient of consolidation. Copyright © 2012 John Wiley & Sons, Ltd.     

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

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

软土大变形固结计算中参数识别的遗传算法

含水量很高的饱和软土在外荷载作用下沉降很大,小变形分析误差太大,必须通过大变形非线性固结计算来模拟。遗传算法是一种全局优化和搜索的仿生算法。近年来随着工程领域中复杂的大规模非线性系统的出现,遗传算法日益得到青睐,目前已经广泛应用到各个领域中。本文对遗传算法做了改进,主要体现在杂交算子的选取和轮盘赌模型的模拟退火拉伸,并将其用于饱和软土的大变形固结分析,解决其中的关键参数识别问题。研究表明该方法是行之有效的,值得进一步研究探索。     

Model for large strain consolidation with non-Darcian flow described by a flow exponent and threshold gradient

Many pieces of evidence have confirmed that the seepage in fine-grained soils can deviate from traditional Darcy's law that contributes to the consolidation theory. In this paper, a numerical model, referred to as consolidation with non-Darcian flow 2 (CNDF2), is developed for the 1D large strain consolidation of a saturated porous medium with the non-Darcian flow. The algorithm accounts for vertical strain, general constitutive relationships, the relative velocity of the fluid and solid phases, variable compressibility and permeability relation during consolidation, time-dependent loading, and unload/reload effects. Compared with the CS2 model proposed in previous studies, some parts of the CNDF2 model are modified in order to adapt to the non-Darcian flow, especially the equivalent series hydraulic conductivity and fluid flow. The verification examples of the CNDF2 model demonstrate excellent accuracy for both small strain and large strain consolidation. According to the applicable conditions of non-Darcian flow law, the CNDF2 model is most suitable for the fine-grained soil. The development of CNDF2 is first presented, followed by examples to analyze the influences of the non-Darcian flow on the consolidation behavior.     

CS2: A PIECEWISE-LINEAR MODEL FOR LARGE STRAIN CONSOLIDATION

This paper presents a piecewise-linear finite-difference model for one-dimensional large strain consolidation called CS2. CS2 is developed using a fixed Eulerian co-ordinate system and constitutive relationships which are defined by discrete data points. The model is dimensionless such that solutions are independent of the initial height of the compressible layer and the absolute magnitude of the hydraulic conductivity of the soil. The capability of CS2 is illustrated using four example problems involving small strain, large strain, self-weight, and non-linear constitutive relationships. In each case, the performance of the model is comparable to other available analytical and numerical solutions. Using CS2, correction factors are developed for the conventional Terzaghi theory which account for the effect of vertical strain on computed values by elapsed time and maximum excess pore pressure during consolidation. © 1997 John Wiley & Sons, Ltd.     

结构性土一维非线性大应变固结半解析解

以结构性较强的天然饱和软黏土为研究对象,考虑了沉积作用对其自重应力的影响,以及压缩性和渗透性的非线性变化,推导了任意加载条件下结构性土一维大应变固结控制方程,并采用半解析的方法对方程进行求解计算。再将其退化为无结构性的饱和软黏土固结解,与已有的大应变固结解进行了对比,验证了该解的正确性。最后将该半解析解计算结果与小应变固结理论解、不考虑结构性的固结理论解计算结果进行对比分析。结果表明:大应变固结理论的沉降计算值大于小应变固结理论的计算值,且二者的差值随着荷载的增加而增加;当考虑土体的结构性时,地表沉降计算值小于不考虑土体结构性的沉降计算值。     

双层饱和软土地基一维大应变固结研究

研究了双层软土地基一维非线性大应变固结问题,同时考虑了土体受沉积作用影响的自重分布以及渗透性在固结过程中的变化等因素。对此复杂问题,利用半解析法进行了求解。通过与已有的大应变固结解析解对比,检验了半解析法在大应变固结分析中的有效性。最后较为详细地分析了双层软土地基一维大应变固结性状,并与相应条件下的小应变固结作了比较。     

碎石桩复合地基大变形非线性固结模型

基于分段线性差分法,建立了一种碎石桩复合地基大变形非线性固结模型SC1。该模型可以进行等应力和等应变条件下桩土复合地基固结分析,同时能够考虑涂抹区、井阻、桩土自重、变荷载,以及固结过程中桩土参数的非线性变化和大变形问题,并编制了 Fortran 计算程序。分别采用基于等应变假设的桩土复合地基解析解和室内试验数据验证了模型的正确性,小变形条件下,该模型数值解与有无井阻及涂抹区的桩土复合地基解析解均能较好吻合;大变形非线性条件下,模型计算结果与室内试验数据基本一致。采用该模型预测实际工程,随固结时间延长,模型等应力与等应变的沉降量计算值相差逐渐增大,现场实测值介于等应力和等应变条件下模型计算结果之间。     

软黏土层一维有限应变固结的超静孔压消散研究

根据土力学固结理论计算分析软黏土层固结过程的超静孔隙水压力值,确定软黏土体固结过程的强度增长,对排水固结法处理软土地基至关重要。软黏土层固结过程中土体变形较大时,有限应变固结理论和小应变固结理论计算分析软黏土固结所得结果差异较大。利用非线性有限元法及程序,通过对软黏土层固结工程算例的计算结果分析,研究了有限应变固结理论和小应变固结理论计算分析软黏土层一维固结超静孔压值消散的差异;探讨了软黏土体一维固结过程中,几何非线性、土体渗透性变化和压缩性变化对超静孔隙水压力消散的影响。研究结果表明,当土体的变形较大时,有限应变固结理论计算出的超静孔压要比小应变固结理论得到的值消散的更快。考虑土体固结过程中渗透性的变化时,超静孔压消散变慢;可用软黏土渗透性变化指数ck 反映渗透性变化对超静孔压消散的影响,渗透性变化指数ck值越小、超静孔压消散越慢。固结过程中软黏土压缩性的大小及变化也影响超静孔压的消散,可用软黏土的压缩指数cc反映固结过程中压缩性的大小及变化对超静孔压消散的影响,软黏土的压缩指数cc越小,固结过程软黏土层中的超静孔压消散越快。     

饱和软土大应变自重固结非线性分析

饱和软土的自重固结涉及许多复杂的非线性问题,而现有分析方法一般基于线性假设,将固结控制方程简化,在沉降变形较大的情况下并不合理。基于Gibson大变形固结方程,推导在单面排水和双面排水条件下的边界条件,综合考虑饱和软土的自重特性与非线性压缩、渗透关系,采用稳定性较好的Crank-Nicolson型差分格式离散并求解非线性控制方程,并从孔隙比、沉降、孔隙水流速3个方面将数值解与CS2固结数值模型进行对比验证。结合工程实例,研究土层不同初始厚度和初始孔隙比对自重固结的影响。结果表明：初始厚度与初始孔隙比较大土层最终沉降也较大。与传统方法相比,该方法所得出的数值解精度较高,实现了高度非线性固结方程的快速求解,其成果更能广泛运用于实际工程的分析之中。     

非饱和土拟二维平面应变固结问题的解析计算方法

基于Fredlund非饱和土一维固结理论,建立了二维平面应变条件下的固结方程组,并得到了单层非饱和土平面应变条件下的解析解。基于相关理论,假设体变系数和渗透系数都为常量,同时考虑到瞬时加载条件下,沿着土体深度方向上产生均匀或者线性分布的初始超孔隙压力,建立了二阶二元偏微分方程组。求解时,引入函数方法来降低方程的阶数,然后通过分离变量法获得方程的通解。在此基础上,结合一个针对单面排水条件下二维平面应变问题案例,通过与数值解对比,验证了所提方法的正确性。并采用所提方法计算获得了二维平面下超孔隙水压力、气压力沿垂直和水平方向消散的等时线,通过计算对比,分析了不同线性分布情况下,初始超孔隙压力对固结消散过程的影响。研究结果表明:初始超孔隙压力的不同分布对超孔隙气压力消散的影响几乎可以忽略,而对超孔隙水压力消散的影响更大。     

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

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

Finite volume approach for finite strain consolidation

The paper presents the finite volume formulation and numerical solution of finite strain one‐dimensional consolidation equation. The equation used in this study utilises a nonlinear continuum representation of consolidation with varying compressibility and hydraulic conductivity and thus inherits the material and geometric nonlinearity. Time‐marching explicit scheme has been used to achieve transient solutions. The nonlinear terms have been evaluated with the known previous time step value of the independent variable, that is, void ratio. Three‐point quadratic interpolation function of Lagrangian family has been used to evaluate the face values at discrete control volumes. It has been shown that the numerical solution is stable and convergent for the general practical cases of consolidation. Performance of the numerical scheme has been evaluated by comparing the results with an analytical solution and with the piecewise piecewise‐linear finite difference numerical model. The approach seems to work well and offers excellent potential for simulating finite strain consolidation. Further, the parametric study has been performed on soft organic clays, and the influence of various parameters on the time ate consolidation characteristics of the soil is shown. Copyright © 2015 John Wiley & Sons, Ltd.     

Axisymmetric consolidation of a poroelastic soil layer with a compressible fluid constituent due to groundwater drawdown

Axisymmetric consolidation of a poroelastic soil layer with a compressible fluid constituent induced by groundwater drawdown was studied based on Biot’s axisymmetric consolidation theory. Laplace and Hankel transforms were employed to solve the governing equation. Explicit analytical solutions are obtained in the Laplace–Hankel transform domain when groundwater drawdown is induced by a constant pumping well. Based on the solutions, numerical computations were performed to study the influences of the compressibility of the fluid constituent on the consolidation behavior of the soil layer.     

基于衬砌长期渗漏水影响的隧道施工扰动诱发超孔隙水压消散及地层固结沉降解

基于Terzaghi-Rendulic固结理论,采用保角映射和分离变量法,计算得到隧道衬砌半渗透漏水边界条件下,盾构施工扰动引起的隧道周围土体超孔隙水压力消散解和地表土体固结沉降,通过工程实例进行验证,发现理论方法与实测数据趋势一致;此外,利用参数分析获得了土体固结沉降和超孔隙水压力的分布影响规律。结果表明:衬砌与土体的相对渗透比越大,固结沉降的初始速率越大,但不影响最终收敛值;土体弹性模量越小,最终固结沉降量越大;土体中超孔隙水压力随着时间增加,在隧道开挖后较短时间里以较大幅度逐渐消散,消散到约为初始值的1/10时减幅放缓;距离衬砌外壁越远,土体初始超孔隙水压力越小,其消散速度也越慢。     

双面排水非线性固结超静孔压的非对称性解析

为了分析一般意义的非线性固结过程中超静孔压的发展规律,基于欧拉描述的基本方程和构形转换关系,提出了非线性固结过程中超静孔压与孔隙比之间的严格函数表达式。在Morris提出的线性有限应变一维固结解析解答基础上,得到双面排水情况下正常固结和自重固结土体的超静孔压解析解。分析和讨论了超静孔压的非对称分布特征及其影响因素。结果表明,无量纲参数N在超静孔压分布形态上起决定作用。在N不接近于0的情况下,按经典Terzaghi固结理论的超静孔压对称性分布假设将不再适用于非线性固结问题。     

地铁荷载下不同固结度软黏土的孔压试验模型

地铁运营前下卧层土体存在不同程度固结,这会影响到运营期隧道周围土体中孔压的变化,进而产生不均匀沉降。通过GDS循环三轴试验系统对杭州饱和软黏土进行动力测试,研究了固结度、固结应力、循环应力比对孔压发展的影响。在现有研究成果的基础上建立了考虑初始固结程度的孔压模型。研究表明：循环荷载下饱和软黏土的孔压发展形态比应变发展形态更具规律性,固结程度愈高的土体孔压发展愈缓慢,并且随振动次数增加,在较低的孔压水平就可达到稳定。预测地铁长期荷载下土体的孔压和沉降时考虑孔压测试时的滞后现象会达到更好的效果,所得的孔压模型能够较好地模拟地铁运营荷载下孔压的发展。     

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.     

变渗透系数软土一维非线性固结沉降数值模拟

固结沉降计算是工程建筑中的重要内容之一。传统的固结沉降计算多是基于渗透数为常数的Terzaghi模型,这与实际的固结过程有较大的差别。采用Carman—Koze模型与一维非线性固结沉降模型进行耦合来模拟固结沉降过程。结果表明,固结导致的渗透系数的变化对固结有较大的影响,特别是在后期这种影响会越来越大。