湛江软黏土三轴剪切蠕变特性研究

湛江软黏土三轴剪切蠕变特性研究

中国沿海地区广泛分布的软沉积物。随着道路建设的高速发展,需要对软黏土的力学性质进行深入研究,因此发展和完善湛江软黏土的变形及变形控制理论显得尤为重要,而解决这些问题的关键之一是建立恰当的力学本构关系模型。工程实例表明,土具有弹性、塑性和流变性,与其它类型土相比,软黏土的流变变形更为显著,建立合理准确的流变模型是土力学研究的重要任务。本文以湛江软黏土为研究对象,通过分析常规三轴剪切蠕变试验数据,建立了湛江软黏土流变本构关系的理论模型和经验模型。1土的破坏过程分析试验仪器为GDS应力路径三轴仪,采用分级加载,试验围压为0、0.1、0.2、0.4MPa。当围压为0.05MPa时,本文根据湛江软黏土的三轴剪切蠕变试验,可得出不同轴向偏应力时应变随时间的变化如图1所示。由图1可知,土样在加载瞬时均有一定量的弹性瞬时应变,其值随偏应力的增大而不变或略有下降;瞬时弹性应变后,蠕变呈衰减稳定,变形随时间增长而趋于稳定。当偏应力继续增大至破坏强度时,瞬时变形后历时较短,即进入加速蠕变阶段并发生破坏,且破坏前有较明显征兆,取破坏强度为0.16MPa。若认为蠕变曲线主要由线性黏弹性变形和线性黏塑性变形构成,对于围压为0.05MPa情况下的应力-应变等时曲线如图2所示,图2曲线可视为折线,折点约为0.028MPa,如果不考虑轴向偏应力为0.16MPa的急剧破坏曲线,则线性黏弹性变形和线性黏塑性变形占主要部分。由以上分析可知,土的屈服应力σs=0.028MPa。建立模型时,线性黏弹性变形用广义Kelvin体描述,线性黏塑性变形用Bingham体或者用胡克体与黏壶串联的Bingham体来描述。同理,围压为0.1MPa时,σs=0.029MPa;围压为0.2MPa时,σs=0.05MPa;围压为0.4MPa时,σs=0.25MPa。图2中的曲线为明显的折线形式,折点出现在第1级加载主应力差值之前,若不考虑该部分,可以发现图2近乎为直线,据此可以认为此围压下的蠕变变形曲线可以用线性黏弹性变形曲线来较好地逼近,即用广义Kelvin体描述。在其它围压下的情况基本相似。2汉江软粘土变形模型的建立2.1发展理论的结构关系模型线性黏弹性模型、西原模型、线性黏弹塑性模型,分别如图3～图5所示。(1)pcp的计算其蠕变变形可表示为:ε(t)=[1EH+1E1(1−exp(−E1tη1)+1E2(1−exp(−E2tη2)]σ(1)ε(t)=[1EΗ+1E1(1-exp(-E1tη1)+1E2(1-exp(-E2tη2)]σ(1)(2)eh+1e111e1表达水平当σ≥σs时,蠕变变形可表示为:ε(t)=[1EH+1E1(1−exp(−E1tη1)+σ−σsση2t]σ(2)ε(t)=[1EΗ+1E1(1-exp(-E1tη1)+σ-σsση2t]σ(2)(3)关于1emp1e的,2t,2t;2t2t;2t2t2t,2t2t,2t当σ≥σs时,蠕变变形可表示为:ε(t)=[1EH+1E1(1−exp(−E1tη1))+σ−σsση2t+σ−σsσE2]σ(3)ε(t)=[1EΗ+1E1(1-exp(-E1tη1))+σ-σsση2t+σ-σsσE2]σ(3)(4)e表达价2,1e3,3.2.3,3.2.3,3.3.3,3.2.3,3.2.3,3.2.3,3.2.3,3.2.3,3.2.3,3.2e1e3,3.2e3,3.2e3,3.2e1,3.23,3.23,3.23,3.23.3,3.2e3,3.2e3,3.2e3,3.2e3,3.2e3,3.2e3,3.2e3,3.3,3.3.3,3.2e3,3.3,3.3.3,3.2e3,3.3,3.3.3.3.3.3.3.3.3.3,3.2e3,3.2e3,3.2e3,3.2e3,3.2e3,3.2e3,3.3,3.2e3,3.3,3.2e3,3.3,3.3.3,3.2e3,3.3,3.3.3,3.2e3,3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3,3.3.3.3,3.3.3.3,3.3.3.3.3,3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3.3其蠕变变形可表示为:ε(t)=[1EH+1E1(1−exp(−E1tη1)+1E2(1−e−E2t/η2)+1E3(1−exp(−E3tη3)]σ(4)ε(t)=[1EΗ+1E1(1-exp(-E1tη1)+1E2(1-e-E2t/η2)+1E3(1-exp(-E3tη3)]σ(4)(5)e1e当σ≥σs时,模型蠕变变形可表示为:ε(t)=[1EH+1E1(1−exp(−E1tη1)+1E2(1−exp(−E2tη2)+σ−σsση3t]σ(5)ε(t)=[1EΗ+1E1(1-exp(-E1tη1)+1E2(1-exp(-E2tη2)+σ-σsση3t]σ(5)(6)e表征1e当σ≥σs时,模型蠕变变形可表示为:ε(t)=[1EH+1E1(1−exp(−E1tη1)+1E2(1−exp(−E2tη2)+σ−σsση3t+σ−σsσE3]σ(6)ε(t)=[1EΗ+1E1(1-exp(-E1tη1)+1E2(1-exp(-E2tη2)+σ-σsση3t+σ-σsσE3]σ(6)2.2各模型拟合精度分析模型建好后,对于各种理论模型由最小二乘原理按高斯-牛顿法进行曲线的非线性拟合,本文应用Matlab程序进行计算机拟合,其拟合参数见表1～表6所列。若用拟合的平均误差表示拟合精度,各模型的拟合平均误差见表7所列,从表7可以比较各模型的拟合精度。从表7可以看出,扩展的七元件线性黏弹性模型比五元件线性黏弹性模型的拟合精度要高,扩展的七元件西原模型比五元件西原模型拟合精度要高;扩展的八元件线性黏弹塑性模型比六元件线性黏弹塑性模型拟合精度要高。由此可看出,线性黏弹性模型的拟合精度最高,西原模型和线性黏弹塑性模型的拟合精度次之,这与蠕变变形中的线性黏塑性变形成份所占比例小、蠕变变形主要以线性黏弹性变形为主及七元件线性黏弹性模型曲线能更好地逼近蠕变曲线有关。由表7还可以看出,线性黏弹塑性模型的拟合精度高于西原模型,这与线性黏弹塑性模型对线性黏塑性变形的描述能力高于西原模型有关。由于线性黏弹性模型的拟合精度已经很高,故建议有条件地选用五元件或七元件线性黏弹性模型来描述0.05MPa围压下的湛江软黏土的蠕变变形。2.3广义kelren模型除了进行围压为0.05MPa的三轴剪切蠕变试验,还进行了围压为0.1、0.2、0.4MPa的试验。围压为0.1MPa时,从0.03～0.295MPa共分10级加载,围压为0.2MPa时,从0.065～0.4MPa分5级加载,围压为0.4MPa时,从0.27～0.75MPa分5级加载,广义Kelvin模型拟合的精度普遍较高,特别是在围压分别为0.05MPa和0.4MPa时,拟合精度略微降低出现在围压分别为0.1MPa和0.2MPa时,这是由于湛江软黏土的蠕变过程中,出现结构性损伤,从而导致产生非线性变形的因素增多,而当应力水平较低或出现高于结构屈服应力时,导致产生非线性变形的因素相对较少,而广义的Kelvin模型相对于西原模型和线性黏弹塑性模型,其描述非线性曲线的能力较弱。限于篇幅,本文不再一一列表详述。2.4模型模拟和相结合模型从以上分析可以看出,广义Kelvin模型描述湛江软黏土蠕变变形,反映了蠕变变形的最主要的特征,效果最好,西原模型和线性黏弹塑性模型能较全面地反映蠕变变形的变化规律,效果次之。若选择广义Kelvin模型描述湛江软黏土蠕变变形,可以通过运用最小二乘法原理得到0.1、0.2、0.4MPa围压下的拟合参数。3经验模型3.1经验发展的本质关系模型(1)应力水平拟合参数ε=aexp(bσD)tλ(7)ε=aexp(bσD)tλ(7)其中,σD为应力水平,20%≤σD≤80%;a、b、λ为拟合参数。Mesri模型对应的蠕变变形可表示为:ε=[aσD1−bσD]tλ(8)ε=[aσD1-bσD]tλ(8)其中,σD、a、b、λ同(7)式。(2)加载时的实际变形ε(t)=A0+Atλ+Bt(9)ε(t)=A0+Atλ+Bt(9)其中,A0为加载时的瞬时变形;λ为曲线走势的参数;B为亚稳定阶段蠕变速率。3.2参数调整同样运用最小二乘法原理,采用Matlab程序进行计算机拟合,其拟合参数见表8、表9所列。3.3软土地变形模型的拟合拟合的平均误差见表10所列。由表10可以看出,幂函数型经验模型比Mesri模型的拟合精度高、使用更为简单,能够拟合的曲线形式更多,若选择幂函数型经验模型描述湛江软黏土蠕变变形,同样通过运用最小二乘法原理,可以得到0.1、0.2、0.4MPa围压下的拟合参数。4结构参数转换工程算例取通过湛江软黏土地区的一条公路,路基宽度为30m,高度为6m,地基土水平方向分析范围为基底宽度的2倍,竖直方向分析范围为地表18m,距离路面中心线10m范围内承受0.150MPa的法向应力。求解思路是在Visualstudio2005上以FLAC3D自带的Burger蠕变模型为蓝本,编写五元件线性黏弹性模型,然后嵌入FLAC3D,用FISH语言编写的线性插值程序可求出每个单元体在不同应力条件下的参数值,应用文献的参数转换方法,将一维参数转化为三维参数后赋于每个单元体,模型中竖直方向的沉降量分布云图如图6所示。以路基底部中心点为研究对象,除了利用五元件线性黏弹性模型求沉降外,利用文中的幂函数经验模型分别求出中心点以下的18个单元体的变形量之和,具体做法是根据每个单元体的应力场和表8～表10,利用线性插值决定每个单元体幂函数经验模型参数的取值,进而得到每个单元体竖直方向的变形量,认为18个单元体竖直方向的变形量之和为路基底部中心点的沉降量,用该方法计算时间与沉降量的关系如图7所示。从图7可以看出,2种曲线相近而略有差别,这主要与2种模型具有不同的拟合效果及不同计算过程等因素有关。5理论模型与经验模型湛江软黏土具有明显的流变特性,其蠕变变形以线性黏弹