循环荷载作用下超固结软黏土变形特性试验研究_刘国清

江西理工大学学报ＪｏｕｒｎａｌｏｆＪｉａｎｇｘｉＵｎｉｖｅｒｓｉｔｙｏｆＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ第36卷第3期2015年6月Vol.36,No.3Jun.20150引言我国东南沿海地区广泛分布着深厚的软黏土，软土地基上修建的高速公路、铁路、跨海大桥等交通工程在循环交通荷载作用下，往往产生过大的工后沉降，严重影响行车舒适度甚至危及行车安全.在实际交通工程中，由于软黏土低强度、收稿日期：２０１５－０２－０６基金项目：国家自然科学基金资助项目(51408441)；浙江省自然科学基金项目（LQ14E080011）作者简介：刘国清（1986-），男，硕士研究生，主要从事土动力学方面的研究，E-mail ：877986008@.通信作者：曾芳金（1965-），男，博士，教授，主要从事地基处理等方面的研究，E-mail ：zfj-hf@.文章编号：2095-3046（2015）03-0036-06DOI:10.13265/ki.jxlgdxxb.2015.03.007循环荷载作用下超固结软黏土变形特性试验研究刘国清1，曾芳金1，郭林2，郑敏1（1.江西理工大学建筑与测绘工程学院，江西赣州341000；2.温州大学建筑工程学院，浙江温州325035）摘要：利用GDS 振动三轴仪对温州饱和软黏土在不同超固结比下进行了循环加载试验，研究了超固结比对温州饱和软黏土应变发展规律、应力-应变滞回曲线特征、累积应变预测方程等的影响.试验结果表明：在相同循环应力比下，超固结软黏土的应变随超固结比的增大而减小，且随超固结比的增大，使软黏土破坏所需的循环应力比呈现出减小的趋势.由于超固结土先期承受了较大的固结压力，随着超固结比的增大，软黏土的应力-应变关系逐渐由粘弹性转变为近似弹性特征.基于不同超固结比下的试验结果建立了累积应变预测方程，阐述了方程中各参数的物理意义及确定方法.关键词：软黏土；超固结比；累积应变；应力-应变曲线中图分类号：TU443文献标志码：AExperimental research on strain behavior of overconsolidated softclay under cyclic loadingLIU Guoqing 1,ZENG Fangjin 1,GUO Lin 2,ZHENG Min 1(1.School of Architectural and Surveying ＆Mapping Engineering,Jiangxi University of Science and Technology,Ganzhou 341000,China;2.Architecture and Civil Engineering College,Wenzhou University,Wenzhou 325035,China )Abstract:Cyclic triaxial tests were conducted on Wenzhou soft clay under different overconsolidation ratio (OCR)through GDS dynamic triaxial apparatus to study the effects of OCR on the strain development,stress-strain hysteresis loop and permanent strain prediction equation.Test results show that:under the same cyclic stress ratio (CSR),the aixial strain deceases with the increase of OCR,and the CSR needed to reach failure tends to decrease with the increase of OCR.With the increase of OCR,the stress -strain relation gradually translates from visco -elastic to nearly elastic because larger confining pressure is subjected for overconsolidated clay.A new permanent strain prediction equation is established based on the test results.The physical meaning and determination methods the equation parameters are discussed in the research.Key words:soft clay;overconsolidation ratio;permanent strain;stress-strain hysteresis loop刘国清，等：循环荷载作用下超固结软黏土变形特性试验研究低渗透性、高压缩性等特点，超载预压正成为一种非常重要的软土地基处理方式，经超载预压处理过的软土表现为超固结特性.为较好预测路基软黏土在交通荷载作用下的沉降，迫切需要开展循环荷载作用下超固结软黏土变形特性的研究.目前，国内外已有越来越多的学者进行了饱和软黏土在循环荷载作用下的变形特性研究.王常晶等[1]利用GDS双向振动仪进行不排水三轴试验研究了双向循环荷载下饱和软黏土变形特性.刘添俊等[2]对珠江三角洲的典型淤泥质饱和软黏土进行了室内循环三轴试验研究了循环荷载作用下饱和软黏土的应变速率.黄茂松等[3]在上海地区典型饱和软黏土不排水循环三轴试验的基础上分析了饱和软黏土的不排水循环累积变形特性.王军等[4]研究了循环荷载作用下超固结软黏土软化孔压模型.国内外还有许多学者[5-10]对软土地基在循环交通荷载下的应变特性进行了深入研究.然而这些研究对象大多数都是正常固结土，超固结软黏土在循环荷载作用下的动力特性研究较少.且现有研究没有阐述超固结软黏土在循环荷载作用下的应力-应变滞回特性[11]，也没有详细的建立考虑超固结比的累积应变计算公式.基于此，笔者利用GDS振动三轴仪对温州饱和软黏土在不同超固结比下进行了循环加载试验，研究了超固结比对温州饱和软黏土应变发展规律、应力-应变滞回曲线特征、累积应变预测方程等的影响.1试验设计1.1试验土样本次试验土样为典型温州饱和软黏土，为减少取土时对土样的扰动，利用特制的薄壁管（直径160mm，高度300mm）在同一深度获取土样，然后将薄壁管两端蜡封运回实验室，在恒温恒湿箱内保存备用.经测定，土样的基本试验指标如下：天然密度ρ=1.57～1.63g/cm3，含水量w=58%～63%，比重G s =2.68～2.71,塑性指数I p=38,液限W L=68，为获取土样的先期固结压力，利用环刀取3个土样分别在WG-2A固结仪中按照12.5→25→50→100→200→400→800kPa的逐级加载次序进行压缩试验，每级加载持续5d，最后分别作出3个土样的e-lg P压缩曲线，利用卡萨格兰德法求出其先期固结压力均值为30kPa（求得3个土样先期固结压力分别为30.1kPa、29.5kPa、30.4kPa）如图1所示.1.2试验仪器本试验所用的仪器是英国GDS公司生产的电机控制动三轴试验系统，主要由驱动装置、压力室、平衡锤、围压控制器、反压控制器等系统组成，控制软件为GDS公司开发的一套非常高端、灵活的GDSLAB控制和数据采集软件，可以实现高精度的循环加载试验.1.3试验方案设计将备用的土样切成直径50mm，高100mm 的圆柱形试样，然后装入GDS三轴压力室进行反压饱和，反压值取280kPa，利用B值检测试样饱和程度，当B值大于0.98时视为试样完全饱和.本次试验方案一共设计了16组对比试验，将试验编号为CO1~CO4、CO5~C10、C11~C16的土样分别在100kPa、150kPa、200kPa的围压下进行初次各向同性固结，固结完成（每小时的排水量小于60mm3）后，将每组试验围压都降到100kPa 再次进行各向同性固结，固结完成后便得到超固结比为1.0、1.5、2.0的试样，最后都在100kPa围压下对每种超固结比下的试样分别施加不同的动应力进行循环加载试验.循环加载试验采用应力控制，波形为半正弦波，频率选用0.1Hz，试样总应变达到20%时视为试样破坏，终止试验.试验方案如表1所示，表中p0、p′0、q cyc分别表示初次固结围压、再次固结及循环加载阶段围压、循环动应力，三者单位均为kPa.定义循环应力比CSR=q cyc/p′0.1.81.61.41.21.00.81101001000图1试验土样e-lg P曲线P c=36kPalg Pr mine第３６卷第３期37试验编号P 0P ′0OCR q cyc CSR C01100100 1.0300.3C02100100 1.0500.5C03100100 1.0600.6C04100100 1.0700.7C05150100 1.5300.3C06150100 1.5500.5C07150100 1.5700.7C08150100 1.5800.8C09150100 1.5850.85C10150100 1.5900.9C11200100 2.0300.3C12200100 2.0500.5C13200100 2.0700.7C14200100 2.0900.9C15200100 2.0100 1.0C162001002.01101.1表1循环加载试验方案2试验结果及分析2.1典型试验曲线循环荷载作用下轴向应变随循环次数的发展曲线如图2所示，其中以N =1500为例作如下定义，总应变是指第1500次循环荷载作用下的应变峰值，用εa,t 表示，总应变在卸载过程中可以恢复的部分称为回弹应变，表示为εa,r ，另一部分不可恢复且随循环次数的增加不断累积的应变称为累积应变，即图中εa,p .为了更清晰地表达总应变、回弹应变、累积应变三者之间的关系，作出了单向循环加载下应力-应变关系曲线如图3所示，在此，我们可以更直观地看出总应变、回弹应变、累积应变之间的关系即εa,t =εa,r +εa,p .2.2不同超固结比下软黏土应变发展规律为了研究循环荷载作用下不同超固结比软黏土应变发展规律，分别作出了OCR 为1.0,1.5,2.0时不同循环应力比下应变随循环次数发展曲线如图4所示.图3单向循环加载下应力-应变关系曲线B CAεa,tεa,rεa,pq c y cε5001000150020002500201612840CSR =0.7CSR =0.6CSR =0.5CSR =0.3εa ,t /%（a ）OCR =1.0CSR =0.9CSR =0.85CSR =0.8CSR =0.7CSR =0.5CSR =0.305001000150020002500（b ）OCR =1.5201612840εa ,t /%图4应变与循环次数关系曲线CSR =1.1CSR =1.0CSR =0.7CSR =0.9CSR =0.3CSR =0.505001000150020002500201612840εa ,t /%（c ）OCR =2.0图2循环荷载下轴向应变随循环次数发展曲线5001000150020002500ε/%N =1500εa ，pεa ，tεa ，r江西理工大学学报2015年6月N /次N /次N /次N /次38表2CSR=0.5和CSR=0.7下的试验数据对比试验编号OCR CSR N/次εa,t C02 1.00.51000 5.5% C06 1.50.51000 1.0% C12 2.00.510000.5% C04 1.00.764020% C07 1.50.76407.0% C13 2.00.7640 2.8%0.00.30.60.9 1.2 1.5 1.880604020N=10N=100εa,t（a）OCR=1.0qcycN=10N=100N=10000.00.20.40.60.8 1.0 1.2εa,t80604020qcyc（b）OCR=1.5图5CSR=0.7时不同超固结比下应力应变关系曲线N=10N=100N=10000.00.20.40.60.8εa,t80604020qcyc（c）OCR=2.0从图4可以看出，当循环应力比较小时，循环荷载作用下，软黏土的应变随循环次数的增加而增大，且这种增大的速率随循环次数的增加逐渐减缓.当循环应力比增大到一定程度，土样则会在较小的循环次数下变形急剧增长，迅速破坏.随着超固结比的增加，试样达到破坏所需的循环应力比逐渐增大.图4（a）中超固结比为1.0的土样在CSR= 0.7时即发生破坏，而图4（b）中超固结比为1.5的土样在CSR=0.85的发生破坏，图4（c）中超固结比为2.0的土样在CSR=1.0的情况下才发生破坏.为了清晰地比较同种循环应力比相同循环次数下，总应变随超固结比的变化情况，根据试验结果，整理出CSR=0.5和CSR=0.7下的6组具有代表意义的试验数据如表2所示.从表2中我们可以发现相同循环应力比下，试样的超固结比越大，一定循环次数下产生的变形越小.如循环应力比为0.5时，经过1000次循环后，OCR=1.0的试样产生总应变为5.5%，而OCR=1.5和2.0时，对应的总应变分别为1.0%和0.5%，减小了82%和91%，递减速度十分明显.循环应力比同为0.7的情况下，超固结比为1.0的试样应变增长的速率明显高于超固结比为1.5、2.0的试样，且达到破坏（总应变达到20%）.而此时超固结比为1.5的土样总应变为7.0%，超固结比为2.0的土样总应变不到3%（仅为2.8%）.2.3不同超固结比下应力-应变滞回曲线特征如前文典型试验曲线图3（单向循环加载下应力-应变关系曲线）中所示，单次循环荷载作用下，软黏土应力应变关系形成一个不封闭的滞回圈，随着循环次数的增大，应力应变滞回曲线逐渐向轴向应变增大的方向移动.为了忽略累积应变的影响，更加直接的比较不同循环次数下的应力应变滞回特性，以CSR=0.7为例，把不同圈数下应力应变曲线起点统一为0点即得到忽略累积应变的应力应变曲线如图5所示.从图5（a）中我们可以明显地看到，随着循环次数的增加，应力应变滞回曲线逐渐向远离纵坐标的方向倾斜，这说明，随着循环次数的增加，超固结土样发生软化，单次循环荷载作用下产生的应变逐渐增大，回弹模量[12]（循环动应力与回弹应变的比值）也相应增加.比较图5（a）、（b）、（c）我们发现OCR=1.0时的各圈应力应变曲线是明显分散开的，OCR=2.0各圈的应力应变曲线几乎重合在一起，OCR=1.5的各圈应力应变曲线分散程度介于两者之间，可见，相同循环应力比下，随着超固结比的增大土样弹性性质表现的愈明显，因为超固结比大的土样曾受过较大的先期固结压力，即本试验中超固结比为2.0刘国清，等：循环荷载作用下超固结软黏土变形特性试验研究第３６卷第３期39OCR =1.0OCR =2.0OCR =1.505001000150020002500201612840εa ,t /%（a ）CSR =0.5图6不同超固结比下累积应变与循环次数关系曲线OCR =1.0OCR =2.0OCR =1.50500100015002000250020151050εa ,t /%（b ）CSR =0.7图7累积应变与循环次数双对数曲线OCR =1.0OCR =2.0OCR =1.50101001000100001001010.1l g a ,plg εa,p （b ）CSR =0.7OCR =1.0OCR =2.0OCR =1.50101001000100001010.1（a ）CSR =0.5l g a ,p图8参数b 与OCR 的关系曲线CSR =0.5试验数据CSR =0.7试验数据拟合直线0.81.01.21.4 1.61.82.02.20.320.300.280.260.240.220.20bOCR的土样在200kPa 的大围压下各向同性固结后强度得到较大提高，土体表现为弹性，随循环次数的增加，单次荷载作用下应变的幅值较稳定,而超固结比相对较小的土样则在相同的循环应力比下随循环次数的增加软化明显，表现出粘塑性,随循环次数的增加，单次荷载作用下应变幅值呈现出增大趋势.2.4不同超固结比下累积应变预测方程绘制出CSR =0.5和CSR =0.7时不同超固结比下累积应变随循环次数的发展曲线如图6，可以看出，相同循环应力比下，超固结比越大，累积应变越小，超固结比对累积应变的发展具有明显的影响.从图6还可以看出不同加载条件下的累积应变发展呈现出典型的指数关系，关于循环荷载作用下超固结软黏土累积应变计算问题目前多采用的是Monismith [13-14]（1975）提出的指学模型：εp ＝C OCR ηm N b（1）对（1）式两边同时取对数有：lg εp =lg C OCR ηm+b lg N（2）式（2）中C OCR 为与OCR 有关的变量，m 为常数，η为循环应力比CSR ，b 为与OCR 有关的变量.将图6（a ）（b ）分别用双对数曲线对应表示为图7（a ）（b ）如下.图7中累积应变εp 与循环次数N 的双对数曲线呈现出较好的线性关系，其中图7（b ）中OCR =1.0的情况下直线短于OCR =1.5和OCR =2.0情况下的原因是试样在N =640时已剪破.这种线性关系恰好可以表示为（2）式的函数.对图中CSR =0.5和CSR =0.7下的两组数据进行线性拟合如图8，可以求出：b =-0.072OCR +0.359（3）此时（2）式可以表示为：lg εp =lg C OCR +m lg η+（-0.072OCR +0.359）lg N（4）即：江西理工大学学报2015年6月N /次N /次N /次N /次40lg εp =lg C OCR +m lg （q cyc p ′0）+（-0.072OCR +0.359）lg N （5）作出OCR =2.0时不同循环应力比下累积应变随循环次数变化曲线如图9所示.图9中直线的截距即为等式（4）中右侧前两项之和，设此截距为H ，则：H ＝lg C OCR +m lg （q cycp ′0）（6）将H 、q cyc /p ′0等值代入（6）式即可求出OCR =2时m =0.195，用同样的方法可以求出OCR =1.5时m =0.212，OCR =1.0时m =0.237，最后取3种超固结比下m 的平均值0.215.对图7中的CSR =0.5和CSR =0.7下的两组数据进行分析，求得不同超固结比下的C OCR ，如图10所示.可见试验的三组数据呈乘幂式C OCR ＝0.4154OCR -1.513（7）故可以得到累积应变计算公式如下：εp =C OCR +ηm N b ＝C OCR （q cyc p ′0）0.215N b （8）其中：C OCR ＝0.4154OCR -1.513，b ＝-0.072OCR ＋0.359综上εp =0.4154OCR -1.513（q cyc p ′0）0.215N -0.072OCR ＋0.359（9）为验证（9）式的可行性，作出CSR 为0.3时，累积应变随循环次数变化的试验结果与计算结果对比图如图11所示，可见利用本文建立的超固结土累积应变计算公式计算出的结果与试验数据吻合较好，建立的考虑超固结比的软黏土累积应变计算模型是可行的.3结论本文利用GDS 振动三轴仪对温州饱和软黏土在不同超固结比下进行了循环加载试验，研究了超固结比对温州饱和软黏土应变发展规律、应力-应变滞回曲线特征、累积应变预测方程等的影响，可以得到以下结论:1）在相同循环应力比下，超固结软黏土的应变随超固结比的增大而减小，且随超固结比的增大，使软黏土破坏所需的循环应力比呈现出增大的趋势.2）相同循环应力比下，由于超固结土先期承受了较大的固结压力，随着超固结比的增大，软黏土的应力-应变关系逐渐由黏弹性转变为近似弹性特征.3）基于不同超固结比下的试验结果建立了累积应变预测方程，阐述了方程中各参数的物理意义及确定方法.图9OCR =2.0时不同循环应力比下累积应变与循环次数关系CSR =0.7CSR =0.51010.10.0111010010001000lg Nb图10参数C OCR 与OCR 的关系曲线CSR =0.5试验数据CSR =0.7试验数据拟合曲线0.50.40.30.20.10.81.01.21.4 1.61.82.02.2OCRC O C R图11CSR =0.3时试验结果与计算结果的对比曲线3210500100015002000OCR =1.5OCR =1.0OCR =2.0计算值εa ,p /%N（下转第73页）刘国清，等：循环荷载作用下超固结软黏土变形特性试验研究第３６卷第３期41参考文献：[1]王常晶，陈云敏.双向循环荷载下饱和软黏土变形特性的试验研究[J].土木工程学报，2010，43(增刊1):573-576.[2]刘添俊，莫海鸿.循环荷载作用下饱和软黏土的应变速率[J].华南理工大学学报，2008，36(10):38-42.[3]黄茂松，李进军.饱和软黏土的不排水循环累计变形特性[J].岩土工程学报，2006，28(7):892-895.[4]王军，蔡袁强.循环荷载作用下超固结软黏土软化孔压模型的研究[J].岩土力学，2008，29(12):3217-3222.[5]边学成，曾二贤，陈云敏.列车交通荷载作用下软土路基的沉降[J].岩土力学，2008，29(11):2991-2995.[6]陈颖平，黄博，陈云敏．循环荷载作用下软黏土不排水累积变形特性[J]．岩土工程学报，2008，30(5):764-768.[7]Miura N,Fujikawa K,Sakai A,et al.Field measurement ofsettlement in Saga airport highway subjected to traffic load [J].Tsuchi-to-kiso,1995,43(449):49-51.[8]Parr G B.Some aspects of the behaviour of London clay underrepeated loading[D].UK:University of Nottingham,1972.[9]蔡袁强，刘新峰，郭林，等.飞机荷载作用下超载预压软土地基的沉降[J].浙江大学学报，2013，47(7):1157-1163.[10]金解放，钟海兵，吴越，等.静载荷与循环冲击作用下岩石损伤变量定义方法的选择[J].有色金属科学与工程，2013，4(4):85-90.[11]郭林，蔡袁强，谷川，等.循环荷载下软黏土回弹和累计变形特性[J].浙江大学学报2013，47(12):2111-2117.[12]郭林，蔡袁强，王军，等．循环荷载下温州结构性软黏土的应变特性研究[J].岩土工程学报，2012，34(12):2249-2254.[13]Monismith C L,Ogawa N,Freeme C R.Permanent deformationcharacteristics of subgrade soils due to repeated loading [J].Transportation Research Record,1975,12(16):1-17.[14]Li D,Selig E T.Cumulative plastic deformation for fine-grainedsubgrade soils [J].Journal Soils of Geotechnical Engineering,1996,122(12):1006-1013.参考文献：[1]Malki H A,Karayiannis N B,Balasubramanian M.Short ‐termelectric power load forecasting using feedforward neural networks[J].Expert Systems,2004,21(3):157-167.[2]马莉，张德丰，许勇.基于动态模糊神经网络的生物工程算法研究[J].计算机工程与科学，2010，32(3):137-140.[3]Chen S,Cowan C F N,Grant P M.Orthogonal least squareslearning algorithm for radial basis function networks [J].Neural Networks,IEEE Transactions on,1991,2(2):302-309.[4]杨兴，朱大奇，桑庆兵.专家系统研究现状与展望[J].计算机应用研究，2007，24(5)5-9.[5]张毅，张思博，李铮.改进的遗传灰色RBF 模型的短期电力负荷预测[J].电测与仪表，2014，51(5):1-4.[6]蒋喆.支持向量机在电力负荷预测中的应用研究[J].计算机仿真，2010，27(8):282-285.[7]伍世虔,徐军.动态模糊神经网络-设计与应用[M].北京:清华大学出版社，2008.[8]王继选，刘小贞，韩中合.基于多约束条件的电力负荷预测模型的研究[J].华东电力，2013，41(3):643-645.[9]张悍东，张莉，汤敏.云模型优化LSSVR 的短期电力负荷预测[J].自动化仪表，2013，34(11):1-5.[10]罗勇，郑金，宁美凤.基于相似日搜索的PSO-WNN 组合模型在短期电力负荷预测中的应用[J].信息与控制，2013，42(3):373-375.[11]杨胡萍，毕志鹏.粒子群优化的灰色模型在中长期负荷预测中的应用[J].电测与仪表，2011,48(2):40-43.[12]张国强，张伯明.基于组合预测的风电场风速及风电机功率预测[J].电力系统自动化，2009，33(18):92-95.[13]冯志国，赖浩.大型住宅小区供配电设计[J].江西理工大学学报，2013，34(5):54-59．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．（上接第41页）任金霞，等：基于动态模糊神经网络的短期电力负荷预测第３６卷第３期73。