非饱和土的本构模型研究

引用平均土骨架应力的概念,研究推导出非饱和土的刚度参数随吸力变化而变化的关系式,进而推导得到用平均土骨架应力表述的非饱和土LC屈服面函数以及硬化规律。从土力学原理推导,得到土样由于在净应力和吸力作用下产生体积变形引起土样饱和度变化的关系式。由平均土骨架应力推广,得到三轴应力状态的椭圆屈服函数,这一非饱和土本构模型的优点在于考虑了应力作用后土样饱和度的变化,通过对已有试验数据的初步验证,表明提出的非饱和土本构模型的合理性和适用性。     

非饱和土的清华弹塑性模型

为研究土在非饱和增湿情况下的应力、应变和强度关系,提出将含水率引入清华模型的硬化参数中建立非饱和土的清华弹塑性模型。通过在干土中掺加冰屑的方法进行增湿试验,表明该模型可以预测不同含水率的非饱和土的应力-变形和强度关系,从天然风干状态增湿到其他含水率的应力-应变全过程的计算结果也与试验结果相符合。将含水率直接引入弹塑性模型而不去研究复杂的基质吸力,可能是非饱和土本构模型工程应用的一条简便途径。     

非饱和土热力学参数确定的探讨

研究非饱和土热力学的力学性质,必须建立能量守恒方程,表征能量的两个基本物理力学参数（体积比热系数和热传导系数）的确定是求解能量守恒方程的重要因素。对于不同颗粒组成的土体,对其干燥颗粒有效热力学参数进行了分析,对固体、液体和气体三相共存的热力学参数进行了研究,并对计算结果进行了对比分析,在此基础上,提出了较适合计算热传导系数的方法,并将此方法和体积平均法、自洽法及Hashin-Strikman方法进行了比较。经与试验结果对比表明：提出的方法能较好地描述非饱和土的热力学参数特性。     

控制吸力循环荷载下非饱和黏性土的弹塑性双面模型

建立非饱和土在动态荷载作用下的本构模型有助于解决相关的工程问题。在塑性增量流动理论框架内,将BBM(Barcelona basic model)模型与笔者提出的塑性硬化准则相结合,考虑控制吸力的循环荷载作用,建立了一个描述非饱和黏性土的弹塑性双面模型。该模型包含一个边界面和一个几何相似的加载面,二者结合记忆中心的概念,共同描述应力路径周期性改变加载方向时非饱和土的非线性、循环塑性和变形累积特性。采用动三轴剪切试验确定了一种非饱和黏质粉土的动态力学特性,分析了吸力、净围压和动应力幅值对循环塑性的影响。模型预测与非饱和土在静、动态加载下的试验结果的对比显示,二者吻合较好,表明所建立的弹塑性双面模型能够模拟非饱和土在控制吸力循环荷载作用下的力学行为。     

非饱和土剑桥模型的基本框架

分析了用有效应力原理和用两个应力变量建立非饱和土力学理论的区别与联系，借助于有效应力原理将饱和土的剑桥模型推广到非饱和土，并与现有的几种塑性模型作了对比分析。     

非饱和土蠕变力学特性试验及经验模型研究

为了预测基质吸力对非饱和粉质黏土蠕变行为的影响,开展其非饱和蠕变经验模型的研究。基于粉质黏土的非饱和蠕变试验,尝试改进Mesri模型和Log-Modified模型。其中Mesri模型和Log-Modified模型的应力-应变关系均采用双曲线函数描述,通过初始切线模量建立其与基质吸力的联系,从而得到应力-基质吸力-应变关系。而Mesri模型和Log-Modified模型的应变-时间关系分别通过幂函数、三参数幂函数描述,分别整合两个模型,得到改进后的Mesri模型和Log-Modified模型。改进后的两个模型为应力-基质吸力-应变-时间模型,利用ε/D-ε关系曲线和幂次关系拟合得到应力-基质吸力-应变关系参数,依据lnε-lnt关系曲线确定改进Mesri模型的应变-时间关系参数,采用BFGS算法和通用全局优化法求解改进Log-Modified模型的应变-时间关系参数。将关系参数分别代入两个改进后的模型,得到相应的非饱和蠕变经验模型。对非饱和蠕变曲线进行辨识,对比两种改进模型的预测精度。结果表明:改进后的Log-Modified模型的预测效果较差,改进后的Mesri模型能较为准确地描述非饱和粉质黏土的蠕变特性。     

非饱和土毛细滞回与变形耦合弹塑性本构模型

在修正剑桥模型的基础上,提出了一个非饱和土毛细滞回与骨架变形耦合的弹塑性本构模型。该模型考虑了基质吸力与饱和度对屈服应力的影响,可以同时描述非饱和土的弹塑性变形特性与毛细循环滞回效应。根据塑性体变的产生使非饱和土进气值增大的特点,建立了变形对土-水特征曲线影响的数学描述。该模型有效地考虑了饱和度对前期屈服应力的作用,准确地反映了土体在不同土-水状态条件下（脱湿和吸湿过程）强度特性的变化,而且还可以有效地描述水力循环历史对土体变形的影响。通过与试验数据对比,证明了该模型能够模拟非饱和土的主要力学特性。     

非极限状态非饱和土主动土压力试验及理论分析

目前土压力研究大都以极限状态下的土体为研究对象,且假定土体处于饱和或干燥状态,未考虑墙体位移与土体非饱和特性对土压力的影响,在实际工程的应用中有局限性。鉴于此,开展主动平动模式下墙后不同含水量砂土的刚性挡墙土压力室内模型试验,并采用渗压计和土压力盒分别量测不同深度处土中的基质吸力和土压力,以及利用DIC图像关联技术观察不同挡墙位移时的土体位移情况。试验结果表明,当墙后土体处于非极限状态时,土体破坏面始终通过墙踵,且其形态接近于平面;其次,在此基础上,结合非极限状态下墙土摩擦强度发挥特性和非饱和土强度准则,提出位移相关的非饱和土强度模型,并建立非极限状态下非饱和土主动土压力计算模型,以及与室内试验结果对比验证了该模型的合理性;最后,针对提出的计算方法探讨了墙体位移和土体基质吸力对主动土压力的影响。参数分析结果表明,非饱和土主动土压力随挡墙位移量的增大而逐渐减小,而随着基质吸力的增大呈现先减小再增加的趋势,且存在一极小土压力值,朗肯土压力值和Fredlund扩展朗肯土压力值分别为该模型在饱和与非饱和情况下位移达到极限状态时的特殊值。     

非饱和黏性土抗剪强度的温度效应试验研究

为了系统地了解温度对黏性土工程性质的影响,本文采用南京地区三种不同矿物成分的黏性土,制成不同含水量和干密度的试样,在5～45℃条件下,开展了抗剪强度试验,获得了三种土的非饱和重塑试样的抗剪强度与温度的关系.试验结果表明:黏性土的黏聚力随温度升高呈线性变化;亲水矿物含量较高的黏性土抗剪强度对温度变化较敏感,黏聚力随着温度...     

非饱和土的水力和力学特性及其弹塑性描述

简单回顾了非饱和土本构模型研究的发展历程,总结了近几年非饱和土弹塑性本构模型最新研究成果,重点介绍了能统一模拟非饱和土水力性状和力学性状耦合的弹塑性本构模型。通过对建立模型过程中的几个核心问题讨论,较详细地说明该类模型的结构、性能以及相关问题。非饱和土水力性状的滞回性用假定存在饱和度弹性区间的弹塑性过程来模拟;该类耦合模型不仅考虑了吸力对非饱和土水力性状和力学性状的影响,还考虑了饱和度对应力-应变关系和强度的影响以及土体变形对土-水特征曲线的影响。用同一套模型参数,耦合模型可统一预测在吸力控制或含水率控制下沿各种应力路径下非饱和土的水力-力学特性,并简单介绍了膨胀性非饱和土的弹塑性本构模型以及耦合模型在有限元数值计算中的应用。     

Dynamics of unsaturated soils using various finite element formulations

Unsaturated soils are three‐phase porous media consisting of a solid skeleton, pore liquid, and pore gas. The coupled mathematical equations representing the dynamics of unsaturated soils can be derived based on the theory of mixtures. Solution of these fully coupled governing equations for unsaturated soils requires tremendous computational resources because three individual phases and interactions between them have to be taken into account. The fully coupled equations governing the dynamics of unsaturated soils are first presented and then two finite element formulations of the governing equations are presented and implemented within a finite element framework. The finite element implementation of all the terms in the governing equations results in the complete formulation and is solved for the first time in this paper. A computationally efficient reduced formulation is obtained by neglecting the relative accelerations and velocities of liquid and gas in the governing equations to investigate the effects of fluid flow in the overall behavior. These two formulations are used to simulate the behavior of an unsaturated silty soil embankment subjected to base shaking and compared with the results from another commonly used partially reduced formulation that neglects the relative accelerations, but takes into account the relative velocities. The stress–strain response of the solid skeleton is modeled as both elastic and elastoplastic in all three analyses. In the elastic analyses no permanent deformations are predicted and the displacements of the partially reduced formulation are in between those of the reduced and complete formulations. The frequency of vibration of the complete formulation in the elastic analysis is closer to the predominant frequency of the base motion and smaller than the frequencies of vibration of the other two analyses. Proper consideration of damping due to fluid flows in the complete formulation is the likely reason for this difference. Permanent deformations are predicted by all three formulations for the elastoplastic analyses. The complete formulation, however, predicts reductions in pore fluid pressures following strong shaking resulting in somewhat smaller displacements than the reduced formulation. The results from complete and reduced formulations are otherwise comparable for elastoplastic analyses. For the elastoplastic analysis, the partially reduced formulation leads to stiffer response than the other two formulations. The likely reason for this stiffer response in the elastoplastic analysis is the interpolation scheme (linear displacement and linear pore fluid pressures) used in the finite element implementation of the partially reduced formulation. Copyright © 2008 John Wiley & Sons, Ltd.     

Use of coupled finite element analysis in unsaturated soil problems

This paper presents a variation of Biot's consolidation theory for analysing problems involving unsaturated soils, and implemented using the finite element method. The numerical method is applied to a few geotechnical problems as examples and the results obtained are compared to some published data. The illustrative examples show how the numerical method can be used to analyse seepage and consolidation problems associated with unsaturated soils and demonstrate the flexibility and applicability of the presented method. Copyright © 2000 John Wiley & Sons, Ltd.     

非饱和黄土的弹塑性软化本构模型

通过对陕西杨陵地区的非饱和黄土进行等应力比压缩试验和常规三轴试验,在对试验结果的应力-应变曲线特性进行分析的基础上,由试验结果直接确定屈服函数,根据Drucker公设,即相适应的流动法则,选择合适的硬化参数,建立了非饱和黄土的弹塑性本构模型,并将模型计算结果与试验结果进行比较。结果表明,该模型能够较好地模拟非饱和黄土与固结应力有关的应变软化现象,计算的变形规律与试验结果基本吻合。     

非饱和土水-力本构模型及其隐式积分算法

在已有工作基础上建立了水力-力学耦合的非饱和土本构模型,在硬化方程中考虑饱和度的影响,同时在土水特征曲线中考虑了塑性体变的影响,从而使模型可以反映非饱和土中的毛细现象与土中弹塑性变形现象的耦合行为。采用隐式积分方法,建立了非饱和土耦合模型的数值模型,并推导了得到了水力-力学耦合的非饱和土的一致切线模量。利用该算法编制了本构模型计算的子程序,使其能向外输出切线刚度矩阵,用于有限元计算。为了验证该算法和程序的正确性,用所编制程序对不同路径下的土体行为进行了预测。通过预测结果与试验结果相对比,表明程序预测结果与试验数据相吻合,模型可以较好地模拟土体的水力-力学耦合行为特性。     

A fully coupled simple model for unsaturated soils

Different phenomena influence the strength and volumetric behavior of unsaturated soils. Among the most important are suction hardening, hydraulic hysteresis, and the influence of volumetric strain on the soil-water retention curves. Fully coupled hydro-mechanical models require including all three phenomena in their constitutive relationships. Among these phenomena, suction hardening is the most influencing as it determines the apparent preconsolidation stress, the position of the loading-collapse yield surface, and the shift of both the isotropic consolidation and the critical state lines. In this paper, a fully coupled hydro mechanical model is presented. It is based on the modified Cam-Clay model but includes a yield surface with anisotropic hardening that takes account of the shift of the critical state line with suction. For highly overconsolidated materials, the sub-loading surface concept has been included in order to increase the precision of the model for these materials. The shift of the retention curves produced by volumetric strains is simulated using a hydraulic model based on the grain and current pore size distribution of the soil.     

Finite element formulation and algorithms for unsaturated soils. Part I: Theory

This paper presents a complete finite‐element treatment for unsaturated soil problems. A new formulation of general constitutive equations for unsaturated soils is first presented. In the incremental stress–strain equations, the suction or the pore water pressure is treated as a strain variable instead of a stress variable. The global governing equations are derived in terms of displacement and pore water pressure. The discretized governing equations are then solved using an adaptive time‐stepping scheme which automatically adjusts the time‐step size so that the integration error in the displacements and pore pressures lies close to a specified tolerance. The non‐linearity caused by suction‐dependent plastic yielding, suction‐dependent degree of saturation, and saturation‐dependent permeability is treated in a similar way to the elastoplasticity. An explicit stress integration scheme is used to solve the constitutive stress–strain equations at the Gauss point level. The elastoplastic stiffness matrix in the Euler solution is evaluated using the suction as well as the stresses and hardening parameters at the start of the subincrement, while the elastoplastic matrix in the modified Euler solution is evaluated using the suction at the end of the subincrement. In addition, when applying subincrementation, the same rate is applied to all strain components including the suction. Copyright © 2003 John Wiley & Sons, Ltd.     

纳米硅水泥土弹塑性有限元分析

试验表明水泥土中掺入适量的纳米硅材料可以大幅度地提高水泥土的抗压强度。在试验的基础上,得出纳米硅水泥土的屈服条件后,采用相关联的流动法则和塑性功硬化条件,建立了纳米硅水泥土的弹塑性本构关系方程。在此基础上,以空间轴对称为研究对象,采用三角形截面的圆环单元,应用虚功原理建立了轴对称条件下纳米硅水泥土有限元方程,用增量刚度法求解,且用Matlab软件编制了有限元程序,在三轴试验条件下验证了该程序的合理性。以某机场水泥土搅拌桩地基处理方法作为工程背景,应用不同纳米硅掺量的水泥土作为地基处理的新材料,与传统水泥土进行了比较,计算结果表明,纳米硅作为水泥土外加剂用于地基处理可有效地减少地基沉降,提高地基承载力,在相同上覆压力下纳米硅水泥土处理后的地基沉降小于普通水泥土材料。     

一个高饱和度非饱和土的本构模型

高饱和度的非饱和土中由于气体处于封闭状态,其内部气压的变化必将对土体的行为产生影响。首先,对高饱和度非饱和土特性进行探讨和研究,随后,在已有非饱和土模型框架基础上,采用广义有效应力原理,建立一个适用于高饱和度条件下的非饱和土的弹塑性本构模型。模型中引入气相耗散的影响,在硬化方程中考虑封闭气体压力改变的影响。最后,利用已有的试验结果来对模型进行验证,并将模型预测结果与前人模型进行对比,表明模型预测可以很好地预测土体的行为,尤其是在高饱和度条件下其结果比其他模型更加接近实际情况。     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.