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.     

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.     

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

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

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

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

Model for large strain consolidation with compressible pore fluid

A numerical model, called CCPF1 (C onsolidation with C ompressible P ore F luid 1 ), is presented for one‐dimensional large strain consolidation of a saturated porous medium with compressible pore fluid. The algorithm includes all the capabilities of a previous large strain consolidation code, CS2, written for incompressible pore fluid. In addition, fluid density and fluid viscosity are functions of fluid pressure in CCPF1. Generalization of the numerical approach to accommodate these functions requires several modifications to the CS2 method, including phase relationships, intrinsic permeability, pore pressure, fluid potential, and mass flux. Inertial forces are neglected and isothermal conditions are assumed. The development of CCPF1 is first presented, followed by an example that illustrates the effects of pore fluid compressibility on the mechanics of consolidation of saturated porous media. Copyright © 2004 John Wiley & Sons, Ltd.     

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

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

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.     

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

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

Electro-osmotic consolidation of soil with variable compressibility,hydraulic conductivity and electro-osmosis conductivity

In present study, the non-linear variations of soil compressibility, hydraulic and electro-osmosis conductivities were analyzed through laboratory experiments, and incorporated in a one-dimensional model. The analytical solutions for excess pore water pressure and degree of consolidation were derived, and numerical simulations were performed to verify its effectiveness. The results indicated that both the non-linear variations of hydraulic and electro-osmosis conductivities showed remarkable impacts on the excess pore water pressure and degree of consolidation, especially for soils with relative high compressibility. A further comparison with previous analytical solutions indicated that more accurate predictions could be obtained with the proposed analytical solutions.     

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

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

Analytical solutions of linear finite‐ and small‐strain one‐dimensional consolidation

Analytical solutions are presented for linear finite‐strain one‐dimensional consolidation of initially unconsolidated soil layers with surcharge loading for both one‐ and two‐way drainage. These solutions complement earlier solutions for initially unconsolidated soil layers without surcharge and initially normally consolidated soil layers with surcharge. Small‐strain solutions for the consolidation of initially unconsolidated soil layers with surcharge loading are also presented, and the relationship between the earlier solutions for initially unconsolidated soil without surcharge and the corresponding small‐strain solutions, which was not addressed in the earlier work, is clarified. The new solutions for initially unconsolidated soil with surcharge loading can be applied to the analysis of low stress consolidation tests and to the partial validation of numerical solutions of non‐linear finite‐strain consolidation. They also clarify a formerly perplexing aspect of finite‐strain solution charts first noted in numerical solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

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

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

指数形式渗流定律下软土非线性大变形固结分析

考虑土中指数形式渗流定律以及土体的非线性固结特性,以超静孔隙水压力为变量在拉格朗日坐标系内建立了软土一维大变形固结问题的控制方程及其求解条件,并运用有限差分法获取其数值解答。在指数形式渗流定律退化为达西定律下,通过将差分解与已有的半解析解进行对比,验证了数值计算的可靠性。最后对指数形式渗流定律下软土一维非线性大变形固结性状进行计算分析,结果表明： 1时,软土的非线性大变形固结速率会随外载增大而减慢; 1时,软土的非线性大变形固结速率会随着外荷载的增加而加快;软土非线性大变形固结速率要比非线性小变形固结速率快,且差别会随荷载增大而加剧;此外,大变形固结理论的最终沉降值要小于小变形固结理论,且差别会随着荷载的增大而加剧。     

基于摩尔-库仑模型的非线性本构模型的开发及其在应变局部化中的应用

以FLAC中的Mohr-Coulomb（M-C）本构模型为基础,在C++环境下实现了考虑拉伸截断的非线性本构模型的二次开发。针对应力及位移的分布规律,将理想弹塑性的Hoek-Brown（H-B）本构模型、M-C本构模型及自定义本构模型的数值解及H-B本构模型的理论解进行了比较,验证了自定义本构模型的正确性。采用弹-脆-塑性的自定义本构模型,以先加载,后挖洞的方式模拟了圆形巷道围岩的应变局部化过程,在计算模型的边界条件对称及均质本构参数的条件下,模拟出了4个对称的V型坑,模拟结果与静水压力条件下岩爆的现场观测及试验结果吻合,而采用接近脆性的H-B应变软化模型所得结果的对称性较差     

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

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

Jaumann率型一维大变形本构关系的双曲线模型

鉴于Jaumann率型本构关系难以得到应力应变全量之间的严格解析形式,采用有限元法研究了一维大变形过程中Kirchhoff应力和Green应变之间的变化规律。首先,基于虚功原理推导了完全拉格朗日描述的大变形有限元方程,然后提出了模拟大变形本构关系的双曲线模型和割线模量处理方法。算例分析结果表明,双曲线模型能够较好地拟合大变形本构关 系。大变形计算本构关系应当考虑初始孔隙比的影响。     

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

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

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

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

Theoretical investigation of the effects of consolidation on contaminant transport through clay barriers

Consolidation of clayey contaminant barriers such as landfill liners has been postulated as a cause of early breakthrough of contaminants. In this paper we theoretically investigate this proposition. For this purpose a sophisticated one‐dimensional, large‐deformation model of coupled mechanical consolidation and solute transport is employed. This new model is a generalization of existing coupled consolidation and solute transport models described in the literature. It takes into account both non‐linearities in geometry as well as constitutive relations. The latter relate the compressibility, hydraulic conductivity and coefficient of effective diffusivity to the deformation of the soil. The model is applied to a case study of a clay liner and geomembrane system. Results obtained from numerical solution of the model equations are compared with those from various simplified models, including a ‘diffusion only’ (i.e. a rigid soil) model traditionally used in contaminant barrier design. For barriers incorporating low compressibility soils (as for many well compacted clays), there is little difference between contaminant transit (i.e. breakthrough) times predicted by the two models. However, for contaminant barriers incorporating more compressible soils, consolidation is shown to significantly accelerate transport. These results indicate the potential importance of accounting for the effects of soil consolidation and highlight the limitations of existing models when modelling solute transport through composite barriers utilizing soft soils. Based on these limited results, we suggest a possible way of taking into account soil consolidation using simplified models. Copyright © 2008 John Wiley & Sons, Ltd.