土的剑桥模型发展综述.docx

土的剑桥模型发展综述.docx

《土的剑桥模型发展综述.docx》由会员分享，可在线阅读，更多相关《土的剑桥模型发展综述.docx（7页珍藏版）》请在冰豆网上搜索。

土的剑桥模型发展综述

土的剑桥模型发展综述

土体本构理论是岩土工程学科的重要基础理论。

随着对土体力学特性的不断深入,塑性理论逐渐被应用于土体本构关系的研究中来。

Roscoe[1]于1963年提出著名的剑桥粘土模型,是应用塑性理论的代表,被看做现代土力学的开端。

在本构理论研究发展过程中,各种建模思想不断涌现,出现了各种不同形式的土体本构模型,但弹塑性模型中得到公认的还只有剑桥模型。

现在国际岩土本构的一大发展趋势是又回到剑桥模型,在剑桥模型基础上进行改进和修正,本文简要介绍了剑桥模型,并对剑桥模型的发展作了较为系统的评述。

1.关于剑桥模型及修正剑桥模型

1958-1963年间,英国剑桥大学的Roscoe等[1]根据正常固结粘土和弱超固结粘土的三轴试验,提出的剑桥粘土的本构模型,标志着人们在土体力学特性认识上的第一次飞跃。

他们将“帽子”屈服准则、正交流动准则和加工硬化规律系统地应用于Cam模型之中,并提出了临界状态线、状态边界面、弹性墙等一系列物理概念,构成了第一个比较完整的土塑性模型。

Roscoc和Burland[2]又进一步修正了剑桥模型,认为剑桥模型的屈服面轨迹应为椭圆,给出了现在众所周知的修正剑桥模型。

可以这样说,剑桥模型开创了土力学的临界状态理论。

试验证明,对于正常固结粘土和弱固结的饱和重塑粘土,孔隙比e与外力p,q之间存在有唯一的关系,且不随应力路径而发生变化。

该模型试图描述室内试验所观察到的现象,即从某一初始状态开始加载直到最终维持塑性常体积变形的临界状态,其基本组成如下:

（1）在（e,p）平面中,存在一条曲线,在正常固结粘性土中的所有应力遵循此路径,这被称为正常固结线（NCL）。

这条线提供了体积硬化规则,可以被广义化为一般应力条件。

（2）在（e,p,q）空间中存在一条线,所有的残余状态都遵循此路径,而与实验类别和初始条件无关。

这条线与（e,p）平面中的正常固结线平行,在此线上,剪切变形发生而没有体积变形发生。

（3）从固结排水和不排水实验中所得到的应力路径位于唯一的状态面,通称为Roscoe面。

事实上,在不排水路径中,土随着塑性体积应变的发展而硬化。

其中,体积应变的弹性和塑性应变增量之和保持常数。

Roscoe面价值在于给出了屈服面类型的一个选择依据。

模型基于对临界状态线、相关联塑性理论中屈服面与固结定律的假定。

该模型假定:

①屈服只与应力球量p和应力偏量q两个应力分量有关,与第三应力不变量无关;②采用塑性体应变硬化规律,以ev为硬化参数;③假定塑性变形符合相关联的流动法则,即g（s）=f（s）;④假定变形消耗的功,即塑性功为

式中:

为塑性偏应变增量。

由

（1）式得到的最初的剑桥模型屈服面形状为子弹头形,屈服函数为

后来提出了修正的假定式（3）来代替

（1）式,即

在此假定的基础上,由（3）式得到的修正剑桥模型屈服函数为椭圆,可以表示为

其中,硬化函数pc为

剑桥模型是当前在土力学领域内应用最广的模型之一,其主要特点有:

基本概念明确;较好地适宜于正常固结粘土和弱超固结粘土;仅有3个参数,都可以通过常规三轴试验求出,在岩土工程实际工作中便于推广;考虑了岩土材料静水压力屈服特性、剪缩性和压硬性。

2.岩土临界状态下的变形特征

所谓临界状态,即在常应力状态和常体积下岩土能够发生稳定的变形或流动;其孔隙比与应力状态之间存在着唯一的关系。

临界状态首先是由Casagrade提出,然后由Roscoe等发展并应用于粘土的本构模型中。

Been等将临界状态应用于砂土本构模型中[3],后来,Noveno等[4]建议将临界状态概念应用于所有的岩土材料（岩石、砂土、粘土和混凝土）中。

岩土材料力学特性的试验研究有大量的文献报道,对这些文献的研究分析表明:

岩土的应力O应变曲线可大致归结为2种典型的形式[5],如图1所示,曲线①代表起始状态疏松的砂或正常固结粘土,曲线②代表起始状态紧密的砂或超固结粘土。

在图1中,有2个坐标系,即qOe1坐标系和evOe1坐标系,q为剪应力,e1为轴向应变,ev为体积应变。

对初始状态疏松的砂土或正常固结粘土,如应力O应变曲线①所示。

在剪切过程中,随着轴向变形逐渐增加,剪应力逐渐增加,体积逐渐减小,即处于压缩状态,到达B点,随着轴向变形的继续增加,剪应力和体积应变逐渐稳定并接近一条水平线,即体积应变率为0,应力率也为0。

对初始状态紧密的砂土或超固结粘土,如应力O应变曲线②所示。

在剪切过程开始段,剪应力较小、应变不大,整个试件体积略微压缩,表明这个阶段颗粒主要被挤向更紧密的排列,与侧限压缩试验的情况差不多;剪应力再增大之后,试件便进入屈服阶段,应变逐渐加快增大,试件体积开始膨胀（或剪胀）,即体积变形先是呈压缩状态后呈膨胀状态;剪应力到达某一峰值A点后,由于颗粒与颗粒间的咬合作用逐渐丧失,剪应力无法继续上升,反而开始下降,剪应力和体积变化逐渐稳定并接近于一条水平线,即B点,此时轴向应变e1和体积应变ev很大,试件进入破坏阶段。

试验曲线②中的A点,常被称为峰值点,对应的抗剪能力称为峰值强度,曲线中的B点,常被称为临界状态点,对应的抗剪能力称为残余强度。

从图1可以看出:

临界状态是岩土的体积保持不变或塑性体积应变率为0的状态,通常对应的应变为10%或更大,对于应变硬化材料,它与破坏强度相对应;对于应变软化材料,它与残余强度相对应。

3.剑桥模型发展

3.1基于经典塑性理论框架的修正

帽盖模型是Dimaggi和Sandier（1976）[6]在剑桥模型的基础上提出的。

该模型不仅能描述塑性屈服前的非线性、剪胀性等特性,还能描述屈服后的各种破坏性状与塑性硬化性状。

国内,魏汝龙（1981）[7]根据不排水三轴压缩试验资料得到的正常固结粘土模型,比修正剑桥模型具有更大的适应性,修正剑桥模型仅是它的特例。

杨林德、张向霞（2005）[8]针对CamOclay模型计算中存在没有充分考虑剪切变形的缺点,首先采用理论分析的手段,指出产生此现象的原因,同时提出由剪应力和体应力引起的应变分量分别采用不同的分配比例系数的思想。

在此基础上对CamOclay模型进行了改进。

3.2非关联流动

Banerjee和Stipho（1978）[9]基于剑桥模型框架,采用塑性增量理论,分析各向同性正常固结及微超固结粘土的不排水应力O应变反映,建立了关联和非关联流动的弹塑性本构模型。

Newson（1998）[10]提出一个软粘土非关联流动的临界状态本构模型,其屈服面与修正剑桥模型相同,塑性势函数为一经验函数,屈服面与塑性势面在偏平面上从低应力时的圆形到破坏时的松岗元屈服面。

3.3次塑性理论与临界状态土力学

Masin（2005）[11]将次塑性理论与临界状态土力学相结合,建立了一个既适用于正常固结也适用于超固结粘土的本构模型。

Masin（2005）[12]还利用次塑性模型预测修正剑桥模型的状态边界面的存在。

3.4有限应变的剑桥模型

Yatomi等（1989）[13]提出了考虑有限应变的剑桥模型,模型中加入了非共轴因素,可模拟土的局部剪切带。

Callari和Auricchio（1998）[14],Borja和Ronaldo等（1996,1998）[15,16]也分别建立了有限应变剑桥本构模型,Ortiz和Pandolfi（2004）[17]基于剑桥模型框架,建立非粘性土本构模型,具有指数型压缩曲线,屈服面及塑性势为椭圆,采用体积应变硬化,模型预测从小应变扩展到有限应变。

3.5时间相关的剑桥模型

Hsieh等（1990）[18,19]利用剑桥塑性理论将粘性土的应力O应变关系加入与时间相关的因素,把总变形分为瞬时应变和粘性应变部分。

瞬时应变利用2个屈服面来计算,一个屈服面为修正剑桥模型的椭圆形,另一个用VonMises的圆柱形屈服面来描述,Arai等（1988）[20],Namikawa（2001）[21]分别结合剑桥模型建立了依赖时间的正常固结粘土的塑性模型,Arai的模型还考虑了粘土各向异性K0固结对应力O应变关系的影响。

Liao等（2000）[22]基于修正剑桥模型与Perzyna粘塑性方程,运用相关联与不相关联的塑性流动法则,建立了一个反映时间效应和剪胀行为的本构模型。

Yin和Graham（1999）[23]基于修正的剑桥模型建立了一个弹粘塑本构模型,这个模型能模拟加速蠕变、卸荷再加载、松驰等软粘土变形行为。

3.6考虑各向异性及结构性

日本学者太田和关口（1979）[24]提出了反映各向异性和应力轴旋转的关口O太田模型,并被广泛应用于日本的工程实践中。

其精华在于引进新的应力比,使剑桥模型成为其特例,能考虑K0固结引起的应力各向异性,和主应力轴旋转产生的塑性变形。

但不能考虑中主应力的影响,即不能考虑真正三维的应力状态。

剑桥模型沿球应力轴（p轴）为等向塑性体变硬化,在日本广泛采用的关口O太田模型[24]沿初始固结线（K0线）为不等向塑性体变硬化。

三轴试验数据表明:

自K0状态向伸长方向剪切时,前者方法计算的体积应变偏小,而后者方法计算的体积应变偏大。

孙德安、姚仰平、殷宗泽（2000）[25]提出一种介于上述两者之间、考虑初始应力各向异性（如K0固结）的不等向塑性体变硬化弹塑性模型。

模型的剪切屈服准则使用SMP准则,该模型能在三维应力下较好地反映土的强度和变形特性,模型的土性参数与剑桥模型一样。

Voyiadjis等（2000）[26]基于Dafalias的各向异性修正剑桥模型建立了一个新的考虑体土的各向异性和结构性的本构方程,利用了一个称为塑性旋转张量的内变量,这个塑性旋转张量是引起各向异性的背应力的函数,模型增加了2个背应力参数。

Wheeler等（2003）[27]利用多阶段三轴排水试验的数据,改进修正剑桥模型,建立了饱和软粘土各向异性本构模型。

Liu和Carter（2002）[28]将土的结构性加入修正剑桥模型,3个新的参数用来描述土的结构性,如果所模拟的土无结构性,该模型等同于修正剑桥模型。

John和Carter等（1995）[29]提出一个土的结构性剑桥模型,这个模型基于临界状态框架,能考虑天然粘土的当前应力状态、应力历史、当前孔隙比和当前土的结构,采用塑性体应变硬化。

针对修正剑桥模型仅适用于正常固结及弱固结粘土,而不适用于严重超固结土的情况,Amerasin和Kraft（1983）[30],Banerjee等（1986）[31]建立了严重超固结土的剑桥模型。

Mita和Dasari等（2004）[32]基于严重超固结土的三轴压缩O伸长试验及平面应变试验,建立了一个伏斯列夫O修正剑桥（HvorsleOMCC）模型,模型具有伏斯列夫面,采用MohrOCoulomb屈服准则,并被扩展到一般三维应力空间。

刘元雪,施建勇（2003）[33]从岩土类材料极限应力状态线所致的各向异性出发,提出了应力空间变换思想。

以修正剑桥模型屈服面的中心为映射中心,给出了重塑土的应力空间变换、应力增量变换的公式,并考虑了应力洛德角的影响。

在变换应力空间中对修正剑桥模型进行了重新表述与改进,该模型可以较好地反映各向异性的影响,反映三轴伸长等应力路径的应力应变特性,也反映某些路径所致的软化现象。

3.7循环荷载下的临界状态模型

Carter和Booker等（1979）[34],Hirai和Hi2royoshi（1987）[35]分别建立了循环荷