删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于固结理论和蠕变模型的软土固结蠕变理论

本站小编 Free考研考试/2020-03-23

孙明乾1, 王清1, 牛岑岑1, 孙铁2
1.吉林大学 建设工程学院,吉林 长春 130000;
2.建设综合勘察研究设计院有限公司,北京 100007
收稿日期: 2015-03-13
基金项目: 国家自然科学基金资助项目(41172236);国家自然科学基金青年基金资助项目(41402243);国土资源部地面沉降监测与防治重点实验室开放基金资助项目(KLLSMP201601).
作者简介: 孙明乾(1987-),男,辽宁丹东人,吉林大学博士研究生;
王清(1959-),女,江苏连云港人,吉林大学教授,博士生导师。

摘要: 通过分级加荷的方式对天津滨海新区经过真空预压处理后的软土进行一维压缩蠕变试验.根据试验所得到的成果,在Terzaghi一维固结理论的基础上,用修正的Singh-Mitchell经验蠕变模型建立了一个新的一维固结微分方程.该方程以一维压缩蠕变试验中所获得的土体力学指标作为计算参数,考虑了软土在主固结阶段结束后所产生的蠕变变形.利用该方程对不同深度和不同应力水平下的试样进行固结和蠕变计算,计算结果与试验结果很接近,表明该方程能够充分反映软土的应力、应变和时间的关系,适用于天津滨海新区的吹填软土地基.
关键词:软土地基蠕变Terzaghi一维固结理论经验蠕变模型Singh-Mitchell模型
Consolidation and Creep Theory Based on Terzaghi Consolidation Theory and Empirical Creep Model
SUN Ming-qian1, WANG Qing1, NIU Cen-cen1, SUN Tie2
1.Construction Engineering College, Jilin University, Changchun 130000, China;
2.China Institute of Geotechnical Investigation and Surveying Co., Ltd., Beijing 100007, China.
Corresponding author: WANG Qing, E-mail:wangqing@jlu.edu.cn
Abstract: Based on the one-dimensional compression creeping and step-loading test on the soft soil subjected to the vacuum preloading treatment in Tianjin Binhai New Area, the modified Singh-Mitchell empirical creep model was used to improve the Terzaghi consolidation theory, and a new one-dimensional consolidation equation was obtained. The new equation used the soft soil mechanics indexes, obtained from the one-dimensional compression and creep tests, as the calculation parameters, and considered the creep deformation generated after primary consolidation. The results obtained from the new equation on different depths of specimens with different stress levels are very close to the test results, which indicates that the new equation can reflect the stress-strain-time relationship of soft soil and is applicable for the soft soil foundation in Tianjin Binhai New Area.
Key Words: soft soil foundationcreepTerzaghi consolidation theoryempirical creep modelSingh-Mitchell model
向海要地一直是人类获取土地资源的重要手段,自从20世纪30年代荷兰围海造陆工程以来,世界上许多沿海地区都在大规模开展围海造陆工程,如日本、阿联酋等国家.为了满足经济快速增长的需要,缓解土地资源日益减少的问题,我国的围海造陆工程在近几十年也取得了非常大的发展.以天津市为例,天津滨海新区自上世纪末,进行了大规模的围海造陆项目,在天津滨海地区形成了大量的软土地基.但是由于围海造陆形成的软土具有明显的蠕变特性,使得天津滨海软土地基在实际工程中的工后沉降量较大,在较长时间内无法达到稳定,甚至产生不均匀沉降,对地基上的建筑物稳定性产生较大的影响.例如根据地面沉降监测,区内南疆港地区在吹填造陆十年后,每年仍有30 mm左右的地面沉降量[1].因此研究天津滨海新区吹填软土地基的固结与蠕变特性,掌握其应力和变形规律,将为天津滨海软土地基的工后沉降问题解决提供理论依据,为围海造陆工程方法的进一步发展解决实际问题,对软土固结与蠕变的机理研究进行进一步的探索.
土体固结理论是由Terzaghi最早提出的[2],这个固结理论有一定的前提条件,并不能完全反映实际工程中土体的应力应变关系.比如Terzaghi固结理论所针对的问题主要为土体的主固结阶段;而软土在次固结阶段中,土体中孔隙水压力已经消散,在施加应力不变的情况下仍然会产生应变,这一现象也可以称为蠕变变形.此后国内外许多学者为了进一步完善土体固结和蠕变理论,均进行了大量的实际工程和室内试验研究,提出了许多土体的固结和蠕变理论模型.Christe在土体应力应变关系分析中应用了多元件组合模型,用来模拟软土的蠕变变形,并对Terzaghi一维固结理论进行了改进[3].Singh和Mitchell基于大量的室内压缩蠕变试验数据,总结出土体蠕变过程中的应力应变变化规律,得到了著名的Singh-Mitchell经验蠕变模型[4].王清等对天津滨海新区围海造陆实际工程的施工工艺、软土的固结蠕变性质,以及软土在吹填沉积过程中微观结构变化等方面均作出了许多有意义的研究[5].
因此,为了研究软土在不同应力条件下的固结和蠕变规律,本文通过分析天津滨海新区吹填软土的一维压缩蠕变试验成果,利用Singh-Mitchell经验蠕变模型,在Terzaghi一维固结理论的基础上构造了一个反映软土次固结阶段应变的一维固结微分方程,对软土的固结和蠕变性质提供一种新的力学模型解释.
1 试验概述试验土样取自天津滨海新区中心渔港经过真空预压加固处理的吹填造陆地区.分别在深度为5,10和20 m处取土样,进行一维压缩蠕变试验.吹填软土地基主要由软黏土组成,基本参数见表 1.
表 1(Table 1)
表 1 基本物理力学参数Table 1 Basic properties of experimental soil
土性 埋深h/m 密度ρ/(g·cm-3) 含水率w/% 孔隙比e1 压缩系数a/MPa-1 液性指数
粉土 5 1.937 28.3 0.785 0.444 1.045
粉质黏土 10 1.975 24.2 0.742 0.305 0.923
粉质黏土 20 1.982 26.8 0.725 0.312 0.735


表 1 基本物理力学参数 Table 1 Basic properties of experimental soil

一维压缩蠕变试验采用WG型高压固结仪.本次试验采用分级加荷的方式,分级加载的加荷等级为25,50,100,200和300 kPa.在实际条件下软土的蠕变往往要持续很长的时间,但是由于进行软土压缩试验的室内条件有限,而且影响因素较多,因此进行较长时间的压缩蠕变试验误差较大,精度有限.在实际试验操作中,通常采用在较短时间内施加较大荷载的办法,缩短蠕变时间,减小误差,并规定一定的应变标准作为加载依据.一般取每级加荷时间间隔以变形稳定为标准,稳定标准取0.005 mm/d, 即1天内的压缩量不超过0.005 mm时,可施加下一级荷载[6].
2 反映软土固结和蠕变性状的应力-应变-时间关系式2.1 考虑蠕变的一维固结微分方程土体蠕变和固结理论是土力学研究中的重点问题,国内外学者都对这一问题做了很多研究.最早奠定土体固结理论研究基础的是Terzaghi在1925年《建立在土的物理学基础的土力学》一书中提出的Terzaghi一维固结理论.Terzaghi认为土体的固结过程即为土中孔隙水压力在应力作用下不断消散的过程.当孔隙水压力减小为零时,土体的应变也达到了最大值.然而对于主要由软黏土组成的软土地基来说,孔隙水压力完全消散、有效应力基本上不变的情况下,随时间发展沉降量继续变大,即发生蠕变现象.由于软土的固结过程与传统的Terzaghi一维固结理论并不相符,因此将土体中孔隙水压力消散前的固结阶段称为主固结阶段,孔隙水压力消散后继续产生应变的阶段称为次固结阶段.在次固结阶段软土继续产生应变的原因是土中的结合水以黏滞流动的状态逐渐排出,使得土颗粒之间的结合水膜厚度减小,土颗粒之间的连接发生改变,在宏观上表现为土体产生了应变,体积减小[7].为了改进传统Terzaghi一维固结理论无法阐述软土在次固结阶段产生应变的不足,本文尝试利用蠕变理论建立新的软土一维固结模型解决这一问题.
由土体常规压缩试验的e-p关系曲线可知
(1)

(2)
式中:Δe为固结产生的孔隙比的变化;Δε=εsεp为固结产生的应变增量,εs为土体在固结蠕变结束后最终的应变量;εp为土体在主固结阶段结束后的应变量;e1为土体固结蠕变开始前的孔隙比;ep为土体在主固结阶段结束后的孔隙比;es为土体固结蠕变结束后的孔隙比.
Terzaghi一维固结理论的提出是在一定前提条件下的.在土体蠕变固结过程中,在土体中取一单位体积的微元体作为研究对象.在固定的时间内,土体微元体体积的变化即孔隙体积的变化.由于土为三相介质,且Terzaghi一维固结理论以饱和土体为研究对象,土颗粒的变形可以忽略不计,因此孔隙体积的变化等于排出水量的体积.
dt时间内孔隙体积V的变化为
(3)
式中:为主固结沉降阶段微元体内孔隙体积的变化;为主固结沉降阶段结束后微元体内孔隙体积的变化;表示微元体体积,在土体固结过程中不发生变化.
根据Terzaghi一维固结理论的前提条件,在蠕变固结过程中土体孔隙比的变化量和有效应力的变化值为正相关关系,且二者比值为压缩系数.这一前提条件在土体蠕变的整个过程中均适用.根据这一前提条件可以将主固结沉降阶段微元体内孔隙体积的变化进行转化:
(4)
式中:σ为总应力;σ′为有效应力;u为孔隙水压力.
由式(2)可得
(5)
式中为土体在主固结沉降阶段结束后的应变速率, 则dt时间内微元体孔隙体积V的变化可以转化为
(6)
根据达西定律,dt时间内微元体内水量Q的变化为
(7)
式中:k为土体的渗透系数;γw为水的容重.
CV即为土的竖向固结系数.由dQ=dV,以及式(6)和式(7)得
(8)
2.2 根据修正Singh-Mitchell经验蠕变方程改进的一维固结微分方程Singh和Mitchell为了研究土的蠕变性质,进行了一系列的蠕变试验,得到了许多软土蠕变过程中应力和应变的变化规律,并提出了著名的Singh-Mitchell经验蠕变方程[8-9].Singh-Mitchell经验蠕变方程表明,软土在较低应力的作用下产生蠕变,土体应变速率和蠕变经历时间在对数坐标中呈现明显的线性相关关系,随着蠕变时间的增加,土体应变速率逐渐减小.Singh-Mitchell经验蠕变方程的基本表达式如下:
(9)
式中:D为应力水平,定义为D=(σ1-σ3)/σ1f-σ3f),f表示破坏时刻;t1为可取单位时间;Aαm为蠕变试验材料常数.
本研究对不同深度的、经过真空预压处理后的吹填软土试样进行一维压缩蠕变试验.根据一维压缩蠕变试验结果,得出了应变速率与时间的关系双对数曲线,如图 1所示.
图 1(Fig. 1)
图 1 应变速率与时间关系曲线Fig.1 Curves of strain rate vs time

根据本研究一维压缩试验结果所得到的关系曲线与Singh-Mitchell的研究结果十分符合,线性关系十分明显,经过拟合可以得到
(10)
式中:图 1中关系曲线的截距,表示软土在蠕变过程中单位时间的应变速率;m图 1中关系曲线斜率的负数.
图 1中可以看出(即曲线的截距)随所加荷载p的增大而增大,由此设想p存在某种联系.可以绘出与应力p的关系,如图 2所示.
图 2(Fig. 2)
图 2 曲线的截距与施加荷载p的关系曲线Fig.2 Relationship between and p

图 2中可以看出,应变速率与时间双对数曲线的截距与施加荷载p呈现很好的线性关系,可以得出其拟合方程:
(11)
式中βγ分别为图 2中应变速率与时间对数关系曲线中的斜率与截距,将式(11)与应变速率与时间对数关系曲线联立可求得
(12)

(13)
式(13)中的计算参数mβγ如之前所述,可以根据一维压缩蠕变试验的结果求出.m值采用式(14)进行计算,减小试验误差.
(14)
式中是由拟合方程式(11)得出的计算值.
则式(13)可以表示为
(15)
由式(8)和式(15)可得
(16)
3 固结微分方程的计算3.1 固结微分方程的参数值根据一维压缩蠕变试验的结果可以得出固结微分方程式(16)所需的参数αβγ.初始孔隙比e1及压缩系数a取值如表 1所示,固结系数Cv取值为天津滨海新区工程经验值0.02 m2/d.吹填软土试样在不同深度、不同荷载作用下的蠕变参数αβγ表 2所示.
表 2(Table 2)
表 2 固结微分方程参数值Table 2 Parameter values of consolidation differential equation
h/m α×103 β γ
25 50 100 200 300
5 -967 -972 -977 -988 -994 0.602 -1.788
10 -931 -953 -959 -970 -988 0.753 -1.673
20 -970 -974 -978 -985 -990 0.577 -0.776
??注:粗体数字20, 50, 100, 200, 300为施加载荷,单位是kPa.


表 2 固结微分方程参数值 Table 2 Parameter values of consolidation differential equation

3.2 有限差分解法采用有限差分法对根据Singh-Mitchell经验蠕变方程改进的一维固结微分方程进行求解.为了保证微分方程在求解过程中格式稳定,采用Crank-Nicolson差分格式.求解后可以得到软土在蠕变过程中的孔压,进而求得软土的蠕变量.对微分方程式(16)取向前差分格式:
(17)
再取向后差分格式,可得
(18)
为了求解固结微分方程式,将采用两种差分格式分解的方程式分别乘以系数1/2,并将两式联立,可以得到式(19).
(19)
式中:
4 计算结果根据3.1节由试验得出土样在不同荷载作用下的参数αβγ,利用新建的一维固结微分方程式对天津滨海新区经过真空预压处理后的吹填软土试样分别作了不同深度、不同应力水平下的蠕变变形计算,绘制应变和时间的关系曲线,并与试验曲线进行了比较,如图 3所示.从图中可以看出利用新建一维固结微分方程式得到的计算值与试验值很接近,说明该方程对于天津滨海新区经过真空预压处理后的吹填软土地基是适用的,能够充分反映软土固结和蠕变性状的应力、应变和时间的关系.
图 3(Fig. 3)
图 3 新建一维固结微分方程计算值与试验值对比曲线Fig.3 Comparison between test data and calculated results obtained from consolidation differential equation (a)—取自5 m深处土样;(b)—取自10 m深处土样; (c)—取自20 m深处土样.

5 结语本文以天津滨海新区吹填形成的软土地基作为研究对象,通过分级加荷的方式对不同深度的软土进行一维压缩蠕变试验.利用试验所得到的一系列软土力学性质参数,对Terzaghi一维固结理论进行了改进.改进参考了Singh-Mitchell经验蠕变模型,并考虑了软土地基次固结沉降阶段发生的蠕变变形.新建的一维固结微分方程得到的计算值与试验值拟合较好,验证了该方程的有效性,而且该方程具有参数少、适用性较强等特点.
参考文献
[1]杜东菊, 杨爱武, 赵建军, 等. 天津滨海软土[M].北京: 科学出版社, 2102.
( Du Dong-ju, Yang Ai-wu, Zhao Jian-jun, et al. Tianjin coastal soft soil[M].Beijing: Science Press, 2102.)(0)
[2] Terzaghi K.Principles of soil mechanics[M].[S.l.]:Engineering News Record, 1925.(0)
[3]Christe I F. A re-appraisal of Merchant's contribution to the theory of consolidation[J].Géotechnique, 1964, 14(4) : 309–320.DOI:10.1680/geot.1964.14.4.309(0)
[4]Singh A, Mitchell J K. General stress-strain-time function for soils[J].Journal of the Soil Mechanics & Foundations Division, 1968, 94 : 21–46.(0)
[5]王清, 肖树芳. 海积软土的工程地质研究现状[J].世界地质, 2000, 19(3) : 253–257.
( Wang Qing, Xiao Shu-fang. A review on the engineering geology of marine soft soil[J].World Geology, 2000, 19(3) : 253–257.)(0)
[6]Yuan Xiao-qing.The laboratory simulation tests of the new vacuum preloading methods and research on settlement prediction[D].Jilin:Jilin University, 2102.
( 苑晓青.新型真空预压法室内模拟试验与沉降的预测研究[D].吉林:吉林大学, 2012.)(0)
[7]Le T M, Fatahi B, Khabbaz H. Numerical optimization to obtain elastic viscoplastic model parameters for soft clay[J].International Journal of Plasticity, 2015, 65 : 1–21.DOI:10.1016/j.ijplas.2014.08.008(0)
[8]Zhu Y B, Yu H M. An improved Mesri creep model for unsaturated weak intercalated soils[J].Journal of Central South University, 2014, 21(12) : 4677–4681.DOI:10.1007/s11771-014-2476-4(0)
[9]Gao H M, Chen G X. Modified Singh-Mitchell creep model for EPS composite soil[J].Advanced Materials Research, 2011, 311/312/313 : 339–343.(0)

相关话题/模型 理论

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于运营系统的出租车出行需求短时预测模型
    林永杰1,邹难21.美国西北大学土木与环境工程系,美国埃文斯敦60208;2.山东大学控制科学与工程学院,山东济南250100收稿日期:2015-04-28基金项目:国家科技支撑计划项目(2014BAG03B04);中央高校基本科研业务费专项资金资助项目(2014JC036)..作者简介:林永杰(1 ...
    本站小编 Free考研考试 2020-03-23
  • 三叉晶界优先形核模型及其动态再结晶仿真
    关小军1,2,付杰1,21.山东大学材料液固结构演变与加工教育部重点实验室,山东济南250061;2.山东大学材料科学与工程学院,山东济南250061收稿日期:2015-06-16基金项目:教育部高等学校博士学科点专项科研基金资助项目(200804220021).作者简介:关小军(1952-),男, ...
    本站小编 Free考研考试 2020-03-23
  • 多时段露天矿生产计划整体优化模型
    孙效玉1,邓鹏宏2,赵明31.东北大学资源与土木工程学院,辽宁沈阳110819;2.鞍钢矿业公司决策支持部,辽宁鞍山114004;3.中煤科工集团沈阳设计研究院有限公司露天工程事业部,辽宁沈阳110015收稿日期:2015-11-16基金项目:十二五国家科技支撑计划项目(2015BAB15B01); ...
    本站小编 Free考研考试 2020-03-23
  • 采空区三维实体模型的交互接口技术
    秦亚光,罗周全,周吉明,汪伟中南大学资源与安全工程学院,湖南长沙410083收稿日期:2015-11-28基金项目:十二五国家科技支撑计划项目(2012BAK09B02-05);国家自然科学基金资助项目(51274250).作者简介:罗周全(1966-),男,湖南邵阳人,中南大学教授,博士生导师。摘 ...
    本站小编 Free考研考试 2020-03-23
  • 有限需求信息下基于最大熵原理的风险厌恶库存模型
    邱若臻,苑红涛,冯俏东北大学工商管理学院,辽宁沈阳110167收稿日期:2015-06-29基金项目:国家自然科学基金资助项目(71372186);教育部人文社会科学研究一般项目(11YJC630165);中央高校基本科研业务费专项资金资助项目(N150604005).摘要:针对风险厌恶的库存决策者 ...
    本站小编 Free考研考试 2020-03-23
  • 基于XML的不确定时空数据模型
    柏禄一,徐长明,刘旭龙,于长永东北大学信息科学与工程学院,辽宁沈阳110819收稿日期:2015-07-07基金项目:国家自然科学基金资助项目(61402087,61401080);河北省自然科学基金资助项目(F2015501049);中央高校基本科研业务费专项资金资助项目(N130323006,N ...
    本站小编 Free考研考试 2020-03-23
  • 基于大数据的C-Mn钢数据预处理及神经网络模型
    吴思炜,曹光明,周晓光,刘振宇东北大学轧制技术及连轧自动化国家重点实验室,辽宁沈阳110819收稿日期:2015-07-28基金项目:钢铁联合基金重点项目(U1460204);辽宁省自然科学基金资助项目(2015020180)。作者简介:吴思炜(1989-),男,辽宁阜新人,东北大学博士研究生;刘振 ...
    本站小编 Free考研考试 2020-03-23
  • 油页岩粉尘层着火的理论模型与实验研究
    杨红霞,李刚,苑春苗,于立富东北大学资源与土木工程学院,辽宁沈阳110819收稿日期:2015-07-28基金项目:国家自然科学基金资助项目(51474053,51374001);国家科技支撑计划项目(2015BAK40B01)。作者简介:杨红霞(1991-),女,河北承德人,东北大学博士研究生;李 ...
    本站小编 Free考研考试 2020-03-23
  • 高压水射流螺旋式切槽辅助松动爆破模型
    粟登峰1,康勇1,2,王晓川1,2,郑丹丹31.重庆大学煤矿灾害动力学与控制国家重点实验室,重庆400030;2.武汉大学水射流理论与新技术湖北省重点实验室,湖北武汉430072;3.四川国泰民安科技有限公司,四川成都610100收稿日期:2015-08-25基金项目:国家重点基础研究发展计划项目( ...
    本站小编 Free考研考试 2020-03-23
  • 马克思主义理论招生名额
    提问问题:马克思主义理论招生名额学院:马克思主义学院提问人:17***46时间:2019-09-2113:24提问内容:老师您好!最近看到学校招生名额已经公布,想请问老师,如果最终推免的招生人数比计划的名额多的话,会不会挤占统考的名额呢?回复内容:你好,没有特殊情况是不允许超过公布计划数录取的。 ...
    本站小编 上海交通大学 2019-11-25