Numerical simulation of elasto-plastic electro-osmosis consolidation at large strain

In this paper, a numerical solution for the electro-osmosis consolidation of clay in multi-dimensional domains at large strains is presented, with the coupling of the soil mechanical behaviour, pore water transport and electrical fields being considered. In particular, the Modified Cam Clay model is employed to describe the elasto-plastic behaviour of clay, and some empirical expressions are used to consider the nonlinear variation of the hydraulic and electrical conductivities of the soil mass during the consolidation processes. The implementation of the theoretical model in a finite element code allows for analysis of the evolution of the transient response of the clay subjected to electro-osmosis treatment. The proposed model is verified via comparison with data from a large strain electro-osmosis laboratory test, to demonstrate its accuracy and effectiveness. Various numerical examples are also investigated to study the deformation characteristics and time-dependent evolution of the excess pore pressure. Finally, a well-documented field application of electro-osmosis is simulated to provide further verification. The results show that the numerical solution is effective in predicting the nonlinear behaviour of clay during electro-osmosis consolidation.     

高岭土电渗固结特性及数值模拟研究

电渗固结是促进低渗透性软土排水固结的有效方法。为了揭示不同电势梯度影响高岭土电渗固结的基本规律,在自制电渗试验装置上对高岭土进行电渗试验。试验过程中测量电流、排水量、沉降量以及有效电压随时间的变化,并进行单位排水能耗分析。基于电渗固结多场耦合控制方程,实现土体电渗固结全耦合分析的有限元数值方法,计算结果与解析解吻合良好,验证了程序的有效性。为预测不同电势梯度下土体沉降量随时间的变化关系,分别对0.5,1.0,1.5 V/cm 3种电势梯度电渗固结试验进行数值模拟分析,获得模型表面沉降量分布、阳极超静孔隙水压力时空发展规律、阳极位置固结度等曲线,计算结果和试验结果吻合良好,可为实际电渗试验提供理论指导。     

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.     

Influence of gas generation in electro‐osmosis consolidation

This paper presents a numerical model for the elasto‐plastic electro‐osmosis consolidation of unsaturated clays experiencing large strains, by considering electro‐osmosis and hydro‐mechanical flows in a deformable multiphase porous medium. The coupled governing equations involving the pore water flow, pore gas flow, electric flow and mechanical deformation in unsaturated clays are derived within the framework of averaging theory and solved numerically using finite elements. The displacements of the solid phase, the pressure of the water phase, the pressure of the gas phase and the electric potential are taken as the primary unknowns in the proposed model. The nonlinear variation of transport parameters during electro‐osmosis consolidation are incorporated into the model using empirical expressions that strongly depend on the degree of water saturation, whereas the Barcelona Basic Model is employed to simulate the elasto‐plastic mechanical behaviour of unsaturated clays. The accuracy of the proposed model is evaluated by validating it against two well‐known numerical examples, involving electro‐osmosis and unsaturated soil behaviour respectively. Two further examples are then investigated to study the capability of the computational algorithm in modelling multiphase flow in electro‐osmosis consolidation. Finally, the effects of gas generation at the anode, the deformation characteristics, the degree of saturation and the time dependent evolution of the excess pore pressure are discussed. Copyright © 2016 John Wiley & Sons, Ltd.     

Electro-Osmotic Consolidation Studies on Peat and Clayey Silt Using Electric Vertical Drain

Peat and soft soil deposits are abundant in Malaysia. Studies using electrokinetic geosynthetics (EKG) such as prefabricated electric vertical drains (EVD or ePVD) have shown the potential for the application of these materials in electrokinetic strengthening of soft clays and increasing solids content in mining tailings and sludges. There are no studies reported on the effectiveness of electro-osmosis using EVD in peat. Laboratory experiments were carried out to investigate the effectiveness of electro-osmotic stabilization of two peat samples and a slightly organic clayey silt sample, from the North Sarawak region. Commercially available EVD was used to induce electro-osmosis and drainage of pore water. The influence of voltage gradient, organic content, initial water content, anode–cathode configuration and presence of roots on electro-osmotic consolidation was investigated. Electro-osmotic consolidation using EVDs was effective in the improvement of peat, over a large range of initial water contents. An increase in the voltage gradient generally improved consolidation. The maximum voltage gradient for optimum results appeared to be in the region 120 V/m. The 2anodes-1cathode configuration performed better than the 1anode-1cathode configuration. The presence of roots reduced the electro-osmosis effect.     

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

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

Multi-dimensional chemo-osmotic consolidation of clays

A model of chemo-osmotic consolidation of clays in multi-dimensional domains is presented, with volumetric strains induced by both changes in the chemistry and osmotically driven pore water flow considered. Three fully coupled governing equations considering force equilibrium, pore water transport and chemical transport are presented and solved using the transient finite element method. The proposed approach is verified via consideration of chemo-osmotic consolidation of a compacted clay landfill liner and then applied to investigation of a hypothetical case with a local leachate leak in the compacted clay liner. An assessment of the impact of the two-dimensional nature of the system is made. The consolidation process is found to be dominated by osmotic consolidation in the early stages and subsequently by chemical consolidation. It is found that the surface settlement and the leachate concentration in the compacted clay liner may be highly overestimated by a one-dimensional analysis. Moreover, the peak negative excess pore pressure predicted by the two-dimensional solution remains in a shallow region under the leak while in the one-dimensional solution it moves progressively downwards.     

Electro-Osmotic Consolidation Experiments on an Organic Soil

A series of electro-osmotic consolidation experiments were carried out on an organic soil from the central Sarawak region of Malaysia. Commercially available EVD was used to induce electro-osmosis and drainage of pore water. The influence of the method of drainage and pumping interval of drained water on electro-osmotic consolidation was investigated. The combination of electro-osmotic consolidation and chemical stabilization techniques was investigated by installing lime or cement columns near the anode end of the test bed. The binder factor in the lime and cement columns was varied. The paper presents the experimental set-up, test procedures, experimental program, results and conclusions of the study. Electro-osmosis, using EVDs, was effective in the electrokinetic treatment of the organic soil. The two types of drainage, namely drainage at bottom and drainage tubes, were effective. Electro-osmotic consolidation improved with decrease in pumping interval. The lime and cement columns, used in combination with electro-osmosis, increased the undrained strength further.     

絮凝−真空−电渗联合加固滩涂软土的模型试验研究

针对真空预压作用下排水板淤堵与排水条件受限等问题,提出絮凝−真空−电渗联合加固法。首先通过沉降柱试验确定合适的有机絮凝剂,然后采用该絮凝剂,分别在 48 h（开始介入真空预压,固结度为0 ）、60 h（排水速率明显下降,固结度为60%）及 84 h（排水速率近乎 0,固结度为 80%）时介入电渗,开展不同电渗介入时间的絮凝−真空−电渗联合加固试验。试验从排水量、十字板剪切强度、含水率与孔压等对比分析联合加固的有效性,确定其最佳电渗介入时间。试验结果表明：当固结度为 80% 时介入电渗,絮凝−真空−电渗联合加固法能够有效地抑制排水速率减小的趋势,增长有效排水时间。同时,土体的抗剪强度和承载力亦得到大幅提升,孔压消散更加均匀。此外,在阳离子聚丙烯酰胺絮凝剂的作用下,初始排水速率快,在一定程度上使土体的渗透性得到提升,有效地解决了排水板淤堵问题,说明絮凝−真空−电渗联合加固法具有较强的优越性。     

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

基于分段线性差分法,建立了一种非饱和土一维大变形固结模型。该模型可考虑土性参数非线性变化,可计算与分析大变形问题,并编制了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.     

Simulation of the progressive failure of an embankment on soft soil

A pragmatic strain-softening constitutive model, which is based on Modified Cam Clay, was applied to the simulation of the progressive failure of an embankment constructed on a deposit of sensitive (strain-softening) clay in Saga, Japan. A comparison of the predictions for this case indicates that if softening is ignored, only relatively small deflections and consolidation settlements are predicted, especially after construction. In contrast, for the case where softening is included in the analysis, progressive failure within the clay induces large shear deformations and finally failure of the embankment is predicted. This comparison suggests that softening-induced progressive failure should be considered in the design of embankments on such soils, and the residual strength of the deposit may have an important influence on the overall factor of safety of the construction. Detailed analyses of predicted excess pore water pressures, shear strains and shear stress levels in the ground indicate that considering the strain-softening process: (a) is associated with the buildup of excess pore water pressure; (b) promotes strain localization; and (c) results generally in a larger zone of soil involved in the failure.     

电渗前后淤泥质土的物理力学性质对比研究

为深入研究电渗排水对土体工程性质的影响,采用铁丝为电极,在模型箱内对南京市某河道淤泥质土进行电渗排水 模拟试验。电渗前后土的物理力学性质对比分析表明：电渗后土体的液、塑限与塑性指数均降低;电渗前后土体在相同含 水量和干密度情况下进行直接剪切试验,发现电渗后土体的凝聚力和内摩擦角均降低,即电渗后土体在特定情况下出现了 强度衰减的现象。论文基于土中阳离子流失机理对土体在特定情况下的强度衰减现象进行了微观解释。另外,利用提出的 新的电渗机理从微观角度解释了电渗过程中土体开裂现象,阳极因快速失水,最早出现裂隙,但裂隙延展范围较小;阴极 含水量始终高于阳极,且失水速度较慢,裂隙延展范围较大。值得一提的是,在以往的研究中并未关注到电渗加固过程中 土中阳离子流失对电渗后土体性质的影响。     

复杂工况单向固结理论分析及可视化

以Terzaghi单向固结理论为基础,对均质土层、分层土层、土层荷载增加、土层厚度变化情况下的孔隙水压力进行了理论分析。建立了以孔隙水压力为未知量的差分形式的控制方程,并通过泰勒级数展开,推导建立了具有二阶精度的边界方程;同时采用MATLAB软件编制相应的通用程序,得到Terzaghi单向固结理论孔隙水压力的数值解。此外,基于按作用消散的平均孔隙水压力的定义通过已经得到的孔隙水压力的数值解计算得到各不同复杂工况下的平均固结度,并且利用MATLAB软件强大的数值计算能力和图形处理功能,实现可视化,对于实际工程中的平均固结度的求解及通过平均固结度来推断地面沉降量有一定的指导和参考价值。     

边坡电渗模型试验及能量分析法数值模拟

采用导电塑料制成的电动土工合成材料(EKG)进行了15.5 d的边坡电渗加固试验,测定了电渗之后土体等含水率分布曲线。试验结果表明,电渗加固效果是随时间从阳极到阴极逐渐扩展的,因此,越靠近阳极接入点的土体加固效果越好,越靠近阴极末端的土体加固效果越差。提出了电渗能量分析法,该方法仅要求黏性土体在排水固结开始的初始时刻是处于饱和状态,而在排水固结进行的过程中土体可以是饱和的,也可以是非饱和的。提出了基于能量分析法的电渗过程数值模拟方法,并对试验过程进行了数值模拟,模拟结果与实测结果能够较为吻合。     

Numerical analysis of an embankment on soft soils considering large displacements

The behaviour of an embankment built on a Portuguese soft soil is analysed considering the material and geometric non-linearity associated with a coupled soil–water formulation. The numerical predictions are compared with the field data in terms of settlements, horizontal displacements and excess pore water pressures. The repercussions of including the large displacements formulation are also studied. It is found that the analysis considering large displacements results in a decrease in settlements and an increase in the rate of excess pore pressure dissipation, both of which are related to the reduction of the thickness of a deformable layer.     

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

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

真空预压加固软土地基变形与固结计算研究

考虑真空度的衰减情况,对真空预压加固软基的变形和固结度计算方法进行了研究。先根据真空预压时砂井的真空压力状态建立了真空预压的空间轴对称变形模型,用位移法推导该模型在等应变条件下的变形及体积应变计算公式。在此基础之上结合Hansbo砂井地基固结理论和真空预压的边界条件,推导了忽略竖向渗流情况下的真空预压加固软土地基的固结解析解。比较文中计算方法与已有计算理论和现场试验实测资料结果表明,体积应变的计算对真空预压的孔压和固结度的计算有较大的影响,而直接引用传统堆载预压的体积应变计算公式计算会导致较大误差。在固结度计算中,采用Hansbo的近似方法能满足计算精度的要求,所得计算结果与实测结果吻合较好     

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

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

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

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