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.     

Analysis of subsidence due to a point sink in an anisotropic porous elastic half space

Analytical solutions of the long-term consolidation deformation and excess pore water pressure due to fluid withdrawal from a saturated anisotropic porous elastic half space are presented. In the analysis, both the permeability and the elastic properties are considered cross-anisotropic. Cases of pervious and impervious half-space boundary are studied. The consolidation settlements as effected by the anisotropy and boundary conditions for pore water are illustrated and discussed.     

Frictionless contact formulation for dynamic analysis of nonlinear saturated porous media based on the mortar method

A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd.     

二维饱和土体动态孔隙率及相关动力响应特性研究

饱和多孔介质的动力响应研究在众多工程领域具有重要意义。充分考虑孔隙率的变化规律与影响因素,有利于合理揭示饱和多孔介质的相关力学行为。为此,将动态孔隙率模型与用于表征饱和多孔介质动力特性的u-U-p型方程结合,构建相应的非线性力学模型,利用Comsol Multiphysis PDE求取相应的数值解,以此研究不同透水条件下,受谐波载荷激励的二维饱和土体的孔隙率、变形量及孔隙水压力的变化规律。结果表明：孔隙率的变化与土骨架的体应变及孔隙水压力直接相关,土体压缩过程中,孔隙率相应减小,土骨架与孔隙流体的相互作用增强,土体运动时所受阻力增大,其无量纲竖向位移小于孔隙率被视为常数时的情况,在此条件下,由于土体的变形量减小,其孔隙水压力也相对减小。故充分考虑动态孔隙率,有利于更加精确地研究等饱土体和多孔介质的相关力学行为。此外,土体上表面透水条件下,孔隙流体可以从土体表面自由排出,土骨架承受的载荷更大,与不透水条件相比,土体孔隙率、竖向位移、孔隙水压力等变化更为显著。     

Large deformation elastic electro-osmosis consolidation of clays

This paper presents the theoretical background of an elastic electro-osmosis consolidation model for saturated soils experiencing large strains, which considers volumetric strains induced by changes in both the hydraulic and electric driven pore water flows. Three fully coupled governing equations, considering the soil mechanical behaviour, pore water transport and electrical field, and their numerical implementation within an updated Lagrangian finite element formulation, are presented. The proposed model is first verified against a classical one-dimensional analytical solution for electro-osmosis consolidation to demonstrate its accuracy and efficiency. Then, various numerical examples are investigated to study the deformation characteristics and time dependent evolution of excess pore pressure. Finally, the importance of considering large strains in a consistent and proper way is demonstrated, and differences compared to models based on small strain theory are highlighted.     

Coupled analysis of dynamically penetrating anchors

The development of a numerical procedure for the finite element analysis of anchors dynamically penetrating into saturated soils is outlined, highlighting its unique features and capabilities. The mechanical behaviour of saturated porous media is predicted using mixture theory. An algorithm is developed for frictional contact in terms of effective normal stress. The contact formulation is based on a mortar segment-to-segment scheme, which considers the interpolation functions of the contact elements to be of order N, thus overcoming a numerical deficiency of the so-called node-to-segment (NTS) contact algorithm. The nonlinear behaviour of the solid constituent is captured by the Modified Cam Clay soil model. The soil constitutive model is also adapted so as to incorporate the dependence of clay strength on strain rate. An appropriate energy-absorbing boundary is used to eliminate possible wave reflections from the artificial mesh boundaries. To illustrate the use of the proposed computational scheme, simulations of dynamically penetrating anchors are conducted. Results are presented and discussed for the installation phase followed by ‘setup’, i.e., pore pressure dissipation and soil consolidation. The results, in particular, reveal the effects of strain rate on the generation of excess pore pressure, bearing resistance and frictional forces. The setup analyses also illustrate the pattern in which pore pressures are dissipated within the soil domain after installation. Hole closure behind a dynamic projectile is also illustrated by an example.     

变荷载作用下轴对称饱和半空间均质地基Biot固结分析

基于Biot固结理论,对随时间变化的变荷载作用下轴对称饱和半空间均质地基固结问题进行了分析。引入状态变量,利用Laplace-Haknel联合变换法,求解了状态方程,得到变荷载作用下饱和土骨架位移、应力、孔隙水压力及渗流量的一般积分形式解,通过算例详细分析了变荷载作用下二维地基固结问题。结果表明,对缓加荷载作用下地基,只有加载较快时,在地基较深处才有明显的Mandel-Cryer效应出现。在循环荷载作用下,多维固结情况下会出现负孔压现象。     

Hybrid time integration and coupled solution methods for nonlinear finite element analysis of partially saturated deformable porous media at small strain

The goal of the paper is to determine the most efficient, yet accurate and stable, finite element nonlinear solution method for analysis of partially saturated deformable porous media at small strain. This involves a comparison between fully implicit, semi‐implicit, and explicit time integration schemes, with monolithically coupled and staggered‐coupled nonlinear solution methods and the hybrid combination thereof. The pore air pressure p_a is assumed atmospheric, that is, p_a=0 at reference pressure. The solid skeleton is assumed to be pressure‐sensitive nonlinear isotropic elastic. Coupled partially saturated ‘consolidation’ in the presence of surface infiltration and traction is simulated for a simple one‐dimensional uniaxial strain example and a more complicated plane strain slope example with gravity loading. Three mixed plane strain quadrilateral elements are considered: (i) Q4P4; (ii) stabilized Q4P4S; and (iii) Q9P4; “Q” refers to the number of solid skeleton displacement nodes, and “P” refers to the number of pore fluid pressure nodes. The verification of the implementation against an analytical solution for partially saturated pore water flow (no solid skeleton deformation) and comparison between the three time integration schemes (fully implicit, semi‐implicit, and explicit) are presented. It is observed that one of the staggered‐coupled semi‐implicit schemes (SIS(b)), combined with the fully implicit monolithically coupled scheme to resolve sharp transients, is the most efficient computationally. Copyright © 2015 John Wiley & Sons, Ltd.     

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

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

开孔管桩复合地基排水固结解析解

现有复合地基的竖向增强体往往无法兼具排水能力与竖向刚度。基于该工程现状,提出一种桩身开孔、能够达到排水效果的新型预制管桩,并针对其开展地基固结分析。假设孔洞流体满足管流条件以及管桩内往孔压连续,并考虑褥垫层的协调作用,推导出排水桩等应变固结解析解。基于该解析解,首先对砂井、碎石桩、不排水桩以及排水桩在相同外径尺寸情况下的排水固结能力进行了对比,结果显示开孔排水桩具有最好的排水固结能力。接着分析了桩身模量、褥垫层模量以及井径比的影响。最后分析了桩身性质（桩身开孔和桩体弹性模量）的影响,计算结果表明,排水能力和模量是影响固结性能的两个主要因素,且排水能力最为重要。但在排水能力相同的情况下,桩身模量越大,则固结越快。     

Analysis of pumping a compressible pore fluid from a saturated elastic half space

A solution method is presented for the consolidation of a saturated, porous elastic half space due to the pumping of a pore fluid at a constant rate from a point sink embedded beneath the surface. It is assumed that the saturated medium is homogeneous and isotropic with respect to both its elastic and flow properties. The soil skeleton is modelled as an isotropic, linear elastic solid obeying Hooke's Law while the pore fluid is assumed to be compressible with its flow governed by Darcy's Law. The solution has been evaluated for a soil with a value of Poisson's ration of 0.25 and for a number of different cases of pore fluid compressibility. It is demonstrated that this compressibility can have a significant influence on the rate of consolidation around the sink. The solutions presented may have application in practical problems such as the extraction of groundwater and other fluids from compressible geological media.     

The axisymmetric consolidation of a semi‐infinite transversely isotropic saturated soil

This paper presents an exact analytical solution to fully coupled axisymmetric consolidation of a semi‐infinite, transversely isotropic saturated soil subjected to a uniform circular loading at the ground surface. The analysis is under the framework of Biot's general theory of consolidation. First, the governing equations of consolidation are transformed into a set of equivalent partial differential equations with the introduction of two auxiliary variables. These partial differential equations are then solved using Hankel–Laplace integral transforms. Once solutions in the transformed domain have been obtained, the actual solutions in the physical domain for displacements and stress components of the solid matrix, pore‐water pressure and fluid discharge can be finally obtained by direct numerical inversion. The accuracy of the numerical solutions developed is confirmed by comparison with an existing exact solution for an isotropic and saturated soil that is a special case of the more general problem addressed. Numerical analyses are also presented to investigate the influence of the degree of material anisotropy on the consolidation settlement. Copyright © 2005 John Wiley & Sons, Ltd.     

A finite element formulation of gradient-based plasticity for porous media with C1 interpolation of internal variables

In this paper a new finite element formulation for numerical analysis of diffused and localized failure behavior of saturated and partially saturated gradient poroplastic materials is proposed. The new finite element includes interpolation functions of first order (C₁) for the internal variables field while classical C₀ interpolation functions for the kinematic fields and pore pressure. This finite element formulation is compatible with a thermodynamically consistent gradient poroplastic theory previously proposed by the authors. In this material theory the internal variables are the only ones of non-local character. To verify the numerical efficiency of the proposed finite element formulation, the non-local gradient poroplastic constitutive theory is combined with the modified Cam Clay model for partially saturated continua. Thereby, the volumetric strain of the solid skeleton and the plastic porosity are the internal variables of the constitutive theory. The numerical results in this paper demonstrate the capabilities of the proposed finite element formulation to capture diffuse and localized failure modes of boundary value problems of porous media, depending on the acting confining pressure and on the material saturation degree.     

Two formulations for dynamic response of a cylindrical cavity in cross‐anisotropic porous media

Two formulations for calculating dynamic response of a cylindrical cavity in cross‐anisotropic porous media based on complex functions theory are presented. The basis of the method is the solution of Biot's consolidation equations in the complex plane. Employing two groups of potential functions for solid skeleton and pore fluid (each group includes three functions), the u–w formulation of Biot's equations are solved. Difference of these two solutions refers to use of two various potential functions. Equations for calculating stress, displacement and pore pressure fields of the medium are mentioned based on each two formulations. Copyright © 2009 John Wiley & Sons, Ltd.     

条形荷载下高饱和土地基的固结与沉降

以Biot（1941）提出的理想土体的三维固结理论为基础,在水流连续条件中增加了考虑孔隙流体本身的压缩变形,并将这一因素以高饱和度土体的饱和度的变化为表征,写入文中建立的平面应变下的Biot二维固结方程。通过实施Laplace变换和Fourier正余弦变换获得基本方程组,求解得到初始函数表达式;再利用数值反变换得到任意时刻、任意一点的数值解,并侧重分析了高饱和土饱和度的变化对孔隙流体的压缩性能的影响,以及由此而产生的固结全过程的变化,表明了考虑孔隙流体的压缩对分析高饱和度土体的固结行为的意义。该解可以进一步推广到其它类型荷载作用在土体不同位置,以及不同边界条件下多层地基的二维固结分析。     

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

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

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

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

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr–Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils.     

饱和土地下源u-P形式解答动力响应计算

一直以来,由Biot孔隙弹性动力方程得到的饱和土地下源Green函数都是u-w形式（u为固相介质位移,w为流相相对于固相的平均位移）。应用两相介质纵波解耦理论,得到了饱和土半空间地下点源荷载的u-P形式（P为孔压）Green函数频域解答;克服了u-w形式Green函数在边界元（BEM）积分时的增根影响。再由Hankel反演,结合Somigliana表象积分,完成BEM计算。并以计算结果分析了地下集中力作用时,饱和土位移、孔压、排水量等动力特性,这对地铁等交通工程、地震工程、土-结构动力相互作用（SSI）的响应计算都具有较重要应用价值。