考虑Hansbo渗流的饱和黏土一维弹黏塑性固结分析

为进一步深入探讨饱和黏土的固结机制,引入考虑时间效应的统一硬化(UH)本构模型描述饱和黏土变形的弹黏塑性,以及Hansbo渗流方程描述固结过程中的非Darcy渗流,修正了饱和黏土一维固结方程,并利用有限差分法进行了求解。通过与文献中一维流变固结试验结果和理论计算结果的对比,验证了所提计算方法和UH模型的适用性。然后探讨了Hansbo渗流参数及UH模型参数对流变固结过程的影响。计算结果表明:土体的黏滞效应是引起流变固结初期地基不排水面附近出现孔压升高现象的主要原因,而渗流的非Darcy特性增强了这一效应,同时土体的这两个特性还延缓了固结中后期饱和黏土地基中孔压的整体消散和地基的沉降。另外,随着回弹指数和超固结参数的增大,或渗透指数的减小,上述流变固结特性更加明显。但是Hansbo渗流参数和渗透指数对地基最终沉降量均无影响。     

基于非牛顿指数渗流和分数阶Merchant模型的一维流变固结分析

为进一步深入分析饱和黏性土的一维流变固结机制,引入Koeller定义的弹壶元件修正Merchant模型,并以此描述土骨架的黏弹性变形行为;同时引入非牛顿指数渗流模型描述流变固结中的非达西渗流,推导出一个新的饱和黏性土一维流变固结方程,并用隐式有限差分法进行数值求解。通过与有关文献中退化模型结果的对比,验证了有限差分算法的有效性。在此基础上,分析非达西渗流模型参数以及分数阶Merchant流变模型参数对流变固结过程的影响。计算结果表明:分数导数阶数和黏滞系数对地基沉降的影响比对孔压消散的影响更为明显,且分数导数阶数越小或黏滞系数越大,地基沉降需要的时间越长。同时,比起达西渗流,非牛顿指数渗流会延缓地基中的孔压消散,使得沉降发展变慢。另外,考虑流变效应时地基沉降要慢于孔压的消散。     

基于分数阶Kelvin模型的饱和黏土一维流变固结分析

引入基于Caputo分数导数的弹壶元件修正Kelvin模型,以描述饱和黏土的一维流变本构关系。沿用Terzaghi饱和土一维固结理论的假设推导流变固结方程,引入Laplace变换和基于Fourier级数展开的Laplace数值逆变换解法进行数值求解。通过与基于整数阶导数模型解析解的对比,证明数值解法的有效性。通过对文献中一维流变固结试验结果的模拟,验证修正Kelvin模型的适用性。然后分析弹壶元件中分数导数阶数和黏滞系数对地基流变固结进程的影响。计算结果表明,在固结开始相当长的一段时间内,孔隙水压的整体消散速度要快于Terzaghi一维固结理论,但在固结后期则会慢于后者;而且在整个固结过程中,地基沉降速率都要慢于后者。总体来看,地基沉降滞后于孔压消散,并且分数导数阶数越小或黏滞系数越大,这种现象就越明显,而且沉降稳定需要的时间越长。     

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

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

饱和软土成层地基一维非线性固结解析解

假定土体在固结过程中压缩性和渗透性的变化成正比,基于 - 及 - 关系,推导出饱和软土成层地基一维非线性固结解析解,分别给出了按沉降定义和按有效应力定义的每层土平均固结度及整个土层总固结度的计算公式。采用Fortran语言编制了相应的计算程序,将计算得到的结果与已有双层地基一维非线性固结解析解计算结果进行比较,验证该解析解的正确性。利用该程序分析成层地基一维非线性固结性状,分别讨论了初始竖向渗透系数、初始体积压缩系数、荷载值及土层厚度对地基固结性状的影响。分析结果表明：在成层地基一维非线性固结过程中,初始竖向渗透系数对超静孔压的影响较为复杂,对上下层地基固结速率影响不同;初始体积压缩系数增大,超静孔压增大,固结速率变小;所加荷载值越大,超静孔压消散越慢,固结发展越慢;超静孔压消散速率不仅取决于土层厚度,同时取决于各层土渗透性的相对大小。     

结构性软黏土时效特性的一维弹黏塑性模拟

为了研究土体结构破坏对软黏土一维变形的时效特性的影响,在Bjerrum的等时间线体系基础上,提出了等黏塑性应变率线等黏塑性应变率线概念,建立了非结构性软黏土的一维弹黏塑性模型;为了描述土体结构渐进破坏特征,定义了结构性参数--结构应变,在非结构性模型的基础上推导了结构性软黏土一维弹黏塑性模型;讨论了通过试验法直接确定模型参数的方法,并利用新建模型对温州天然软黏土的一维常规压缩试验、天然Ariake 黏土的分级快速固结试验、结构性Berthierville黏土的一维等应变率压缩试验及长期蠕变试验进行模拟。模拟与试验结果的对比表明,该模型能较好地描述结构性软黏土一维压缩变形的时效特征。     

主应力轴旋转下K0固结饱和粉质黏土孔压及变形特性试验研究

主应力旋转将加速孔压和塑性应变的累积,影响土体的力学特性。为研究主应力轴循环旋转对 固结饱和粉质黏土孔压和变形的影响,利用GDS-HCA试验系统开展了不同应力路径的 固结饱和粉质黏土不排水循环三轴试验和循环扭剪试验,研究了主应力轴循环旋转对 固结饱和粉质黏土累积塑性应变、孔压等动力特性的影响。结果表明：动剪应力幅值、心形应力路径所包围的面积、循环偏应力幅值等均随循环剪应力比和循环动应力比的增大而增大。试样未发生破坏时, 固结饱和粉质黏土的孔压比随着振动荷载作用次数的增加而增大;当试样发生破坏时,孔隙压力迅速消散,孔压比急剧下降。循环扭剪试验时试件产生的轴向累积塑性应变始终大于循环三轴试验产生的累积塑性应变,说明主应力轴循环旋转会加速试件轴向应变的累积。试件的累积塑性应变随循环剪应力比和循环动应力比的增大而增大。循环动应力比越大时,主应力轴旋转造成的轴向累积塑性应变差异越显著。在Monismith幂次函数模型的基础上,建立了主应力轴循环旋转条件时K0固结饱和粉质黏土累积塑性应变方程,并对模型的可靠性进行了分析和验证。     

人工岛沉降的离心机模型试验与数值分析

采用离心机模型试验与数值方法,对某一人工岛深厚软黏土地基的沉降特性进行了研究。在模拟现场施工条件下,通过离心机模型试验,研究了60 m厚软黏土地基的沉降特性。基于离心机模型试验,建立了二维有限元模型,应用弹黏塑性本构模型对软黏土地基的沉降规律进行了分析。通过模型试验与数值计算结果的对比分析,结果表明,模型试验与数值计算得到的工后沉降曲线基本吻合,初期沉降速率较大,后期沉降速率相对稳定。在浅部土层,计算得到的超孔压消散快于模型试验结果,而深部土层两者超孔压消散规律基本一致。二者规律的一致性验证了所用数值方法的可靠性。     

考虑非达西渗流的双层软土地基大变形非线性固结分析

现有的双层软土地基大变形固结理论均假定土中渗流遵循达西定律,但软黏土中渗流在低水力坡降下会出现偏离达西定律的现象已逐渐为人们认识。综合考虑双层软黏土中的非达西渗流、软土的非线性压缩渗透特性及实际中的变荷载作用,在拉格朗日坐标系中建立以超静孔压为变量的软土大变形非线性固结模型,并给出其有限差分解。在此基础上,通过与已有的非达西定律下单层地基大变形固结数值解计算结果相对比,以验证其差分解的可靠性。最后着重分析上、下层非达西渗流参数m_1、i_(11)及m_2、i_(12)对固结性状的影响,并分析了在大、小变形不同几何假定下对双层地基超静孔压消散及固结沉降影响的异同。结果表明:上层土参数m_1、i_(11)对固结性状的影响要比下层土参数m_2、i_(12)显著;m_2、i_(12)的增大会引起上层土超静孔压消散速率的加快,但双层地基的固结速率会减慢;大变形假定下双层软土地基的固结速率要比小变形假定下快,但两种几何假定下双层地基的最终沉降量是相等的。     

考虑非Darcy渗流时循环荷载下饱和黏土一维固结分析

考虑到土体处于超固结状态下的压缩性一般要比正常固结状态下的低,引入描述非Darcy渗流的Hansbo方程,修正了Terzaghi饱和黏土一维固结方程,并将其推广到低频循环荷载作用的情况。采用有限体积法对该方程进行了求解,并讨论了非Darcy渗流参数、循环荷载周期以及超固结状态下土体压缩性对固结进程的影响。计算结果表明,在矩形波载作用下,按孔压定义和按变形定义的固结度都随时间增长而震荡增加,且按孔压定义的固结度的震荡幅度明显大于按变形定义的固结度。另外,非Darcy渗流减缓了地基的沉降速率,且循环荷载周期越短,或超固结状态下压缩性越高,地基的沉降速率就越慢。     

竖井地基热排水固结模型试验及有限元模拟

设计竖井地基热排水固结模型试验系统,选取宁波软黏土在水热循环温度70℃下开展热排水固结模型试验,分析了无堆载加热、分级堆载、恒载降温阶段各测点的温度、孔隙水压力及地表沉降变化规律,并以试验为原型,利用COMSOL软件建立热排水固结耦合模型,进行了有限元模拟分析。结果表明:模拟结果与试验结果可以相互验证;无堆载加热阶段,土体孔隙水压力先增大后消散,地表沉降先隆起后下沉,地基发生固结压缩;分级堆载阶段,耦合模型地基固结速率较退化模型快,地基平均固结度到达90%时,耦合模型所用时间较退化模型显著减少;恒载降温阶段,降温能使残余孔压快速消散,土体有效应力增加,土体沉降继续增大,地基进一步固结压缩。     

冲击荷载下饱和粘土孔压特性初探

冲击荷载作用下，孔压的增长与消散规律是研究动力固结法加固饱和粘土地基的基础与理论支持，就此问题进行了冲击荷载作用下饱和粘土孔压增长与消散规律的一维模型试验，试验中考虑了不同击数N和不同击能WH的影响，试验的结论新颖，对其进行了初步探讨。     

弹黏塑性模型及其在预压处理软基沉降计算中的应用

以准确计算预压荷载作用下软土路基的工后沉降为目的,分析了不同预压法加固软土地基时土体的应力变化。提出了用弱化应力路径变化过程的弹黏塑性模型表示软黏土的应力-应变关系。结合经典比奥固结理论,推导了预压荷载作用下土体沉降计算的有限单元法。通过具体算例,验证了所提方法。结果表明：（1）不同预压方法的加固机制和应力路径均不相同。堆载预压法加固软基是一不等向的正压固结过程,真空预压法则是在负压作用下的等向固结过程;（2）弹黏塑性模型是不受瞬时施加应力的大小及具体应力路径变化过程影响的软土流变模型,运用于预压荷载作用下的软基时具有不考虑应力路径和不划分主次固结的优越性;（3）在比奥固结理论中,用弹黏塑性模型反映软土的物理方程时,可计算软基考虑了弹黏塑变形性质的固结沉降和工后沉降;（4）弹黏塑性模型的软土体黏性参数? /V和塑性参数? /V对竖向位移都有比较明显的影响。     

考虑温度影响的改进西原流变模型及一维固结解

随着我国热能源开发,热力管道工程大规模应用,由此带来的岩土体应力、变形、渗流和温度耦合作用问题已经成为岩土工程领域的研究热点。在长期荷载作用下,岩土体的变形和强度随时间而变化,尤其是软土,具有明显的流变特征。传统的各流变模型并未考虑升温方式和温度对其各参数的影响。为了能更好地反映流变固结特性,考虑软黏土的黏弹塑性,在西原模型基础上,引入温度膨胀系数,建立了热力耦合作用下的改进西原模型,给出了其应力-应变关系,并利用Laplace变换和逆变换求解了一维热固结方程,得到了其解析解。计算结果表明:固结压力一定时,温度升高加快土体固结;温度一定时,土体中孔压随固结压力增大而减小;改进西原模型的弹性模量越大,孔压消散速率越大,土体固结越快;土颗粒受温度影响而产生的热膨胀对于孔压消散的作用非常小;黏滞系数越小,孔压消散更快,土体固结越快;试验结果与考虑土体的黏弹塑性的热力耦合改进西原模型计算结果吻合,改进西原模型可以较好地描述土体的热力耦合特性。     

基于一维和三维固结挤土桩沉降时效计算分析

饱和软黏土地基中的挤土桩,施工引起的桩侧土体超孔隙水压力消散和竖向荷载引起的桩底土体固结是导致其沉降的主要因素。基于三维固结理论,研究了桩侧地基土再固结沉降的变化规律,并从桩-土相互作用的原理出发,结合边界条件,推导出了在桩侧地基土再固结沉降的作用下挤土桩沉降的计算公式。运用太沙基一维固结理论并采用分层总和法计算竖向荷载作用下桩底土体的沉降。提出了饱和软土中挤土桩沉降时效的计算方法。并对工程实例进行了计算分析,其结果符合工程实际沉降规律。     

考虑起始坡降双层地基的大应变非线性固结

天然软土成层分布特性及土中渗流存在起始水力坡降的现象已被人们熟知。但变荷载下能同时考虑黏土中起始水力坡降、软土非线性压缩渗透特性及大应变特性的双层地基固结理论还鲜见报道。在拉格朗日坐标系中建立以超静孔压为变量的双层软土地基大应变非线性固结模型并给出其有限差分解。通过与考虑起始水力坡降的单层地基大应变非线性固结数值计算结果对比,验证了差分解的可靠性。着重分析了上、下土层起始坡降无量纲参数R1、R2对双层地基固结性状的影响,分析在大应变与小应变假定下双层地基超静孔压消散及固结沉降变形的异同。结果表明：上层土无量纲参数R1对双层地基固结性状的影响程度较下层土无量纲参数R2显著;大应变假定下双层地基渗流前锋的下移速度要快于小应变假定下的移动速度;大应变假定下考虑起始水力坡降的双层软土地基超静孔压消散速率要比小应变假定下快,且大应变假定下考虑起始水力坡降的双层地基最终沉降量要比小应变假定下大。     

循环荷载作用下考虑非达西渗流的软黏土一维流变固结分析

基于Hansob渗流理论和Merchant流变模型,将微分型流变方程耦合到一维固结理论中,建立了循环荷载作用下双变量表述的软黏土一维流变固结耦合方程组。采用FlexPDE软件对方程组进行了求解,通过与已有文献结果对比,验证了算法的可靠性,在该基础上研究了循环荷载作用下饱和软黏土的一维流变固结特性,分析了Hansob渗流参数、Merchant流变模型参数以及循环荷载类型对饱和软黏土固结特性的影响。研究结果表明,按超孔压定义的平均固结度 和按变形定义的平均固结度 均呈现循环变化特征,且 > ;当循环周次较大时, 进行稳定循环状态,由于流变效应, 峰值仍呈逐渐增大趋势;随着Hansob渗流参数 和 的逐渐增大,地基固结速率逐渐降低,同一循环周次内 、 的峰值也随之降低;随着流变参数F的逐渐增大,同一循环周次内 的峰值逐渐降低, 的峰值显著提高;随着 的逐渐增大,同一循环周次内 、 的峰值逐渐降低。循环荷载作用下地基的平均固结度始终处于循环状态,无法达到固结完成状态,其最大平均固结度也远小于线性荷载的情况。     

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

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

一个描述软黏土时效特性的一维弹黏塑性模型

为了描述软黏土一维应力-应变关系的时效特性,基于Bjerrum的等时间线体系,提出等黏塑性应变率线概念,推导了黏塑性应变率与黏塑性应变增量的关系,建立了软黏土的一维弹黏塑性模型;从理论上分析了新建模型与3种典型的一维弹黏塑性时效本构模型的内在联系,表明新建模型与其他3种模型在本质上是等效的,且形式更简洁,物理意义更明确;利用新建模型对软黏土的固结-蠕变耦合效应、应变率效应、应力松弛效应等时效特性进行了理论分析,并得到了相应的解析解;结合宁波软黏土的一维固结试验,阐述了模型参数的确定方法,并用新建模型对宁波软黏土的固结-蠕变试验、温州软黏土的一维多级等应变率试验、香港海相软黏土的一维应力松弛试验进行模拟,验证了新建模型的有效性。研究结果表明,新建模型能很好地模拟软黏土的一维时效特性。     

循环荷载作用下黏弹性地基一维固结性状研究

针对单层黏弹性地基模型,对循环荷载作傅立叶级数展开,由固结方程求得循环荷载作用下不透水 (透水) 边界饱和软黏土一维固结解析解.通过不同渗透系数、黏滞系数和弹性模量以及黏弹性与弹性不同地基模型比较等对有效应力的影响,进行了固结性状分析.结果表明、在循环荷载作用下地基中各点的有效应力并不随荷载的变化而同步变化,而是按一定规律滞后发展;有效应力最终达到一个稳定状态,每一个加载卸载循环下有效应力幅值在一个振荡周期的平均值变化趋近于0;软黏土的流变特性对固结的影响主要发生在中、后期.