土的剑桥模型发展综述

1、土的剑桥模型发展综述土体本构理论是岩土工程学科的重要基础理论。随着对土体力学特性的不断深入 ,塑性理论逐渐被应用于土体本构关系的研究中来。Roscoe1于1963 年提出著名的剑桥粘土模型,是应用塑性理论的代表 ,被看做现代土力学的开端。在本构理论研究发展过程中 ,各种建模思想不断涌现 ,出现了各种不同形式的土体本构模型 ,但弹塑性模型中得到公认的还只有剑桥模型。 现在国际岩土本构的一大发展趋势是又回到剑桥模型 ,在剑桥模型基础上进行改进和修正 ,本文简要介绍了剑桥模型 ,并对剑桥模型的发展作了较为系统的评述。1. 关于剑桥模型及修正剑桥模型1958 - 1963 年间 , 英国剑桥大学的 R

2、oscoe等1 根据正常固结粘土和弱超固结粘土的三轴试验 , 提出的剑桥粘土的本构模型 ,标志着人们在土体力学特性认识上的第一次飞跃。他们将 “帽子 ”屈服准则、正交流动准则和加工硬化规律系统地应用于 Cam 模型之中 ,并提出了临界状态线、状态边界面、弹性墙等一系列物理概念 ,构成了第一个比较完整的土塑性模型。Roscoc 和Burland2 又进一步修正了剑桥模型 ,认为剑桥模型的屈服面轨迹应为椭圆 ,给出了现在众所周知的修正剑桥模型。 可以这样说,剑桥模型开创了土力学的临界状态理论。试验证明 ,对于正常固结粘土和弱固结的饱和重塑粘土,孔隙比 e 与外力 p , q 之间存在有唯一的关系

3、,且不随应力路径而发生变化。该模型试图描述室内试验所观察到的现象,即从某一初始状态开始加载直到最终维持塑性常体积变形的临界状态,其基本组成如下 :(1) 在( e , p) 平面中 ,存在一条曲线 ,在正常固结粘性土中的所有应力遵循此路径 ,这被称为正常固结线 (NCL) 。这条线提供了体积硬化规则 ,可以被广义化为一般应力条件。(2) 在( e , p , q) 空间中存在一条线 ,所有的残余状态都遵循此路径 ,而与实验类别和初始条件无关。这条线与 ( e , p) 平面中的正常固结线平行 ,在此线上 ,剪切变形发生而没有体积变形发生。(3) 从固结排水和不排水实验中所得到的应力路径位于唯一

4、的状态面,通称为 Roscoe 面。事实上 ,在不排水路径中 ,土随着塑性体积应变的发展而硬化。其中 ,体积应变的弹性和塑性应变增量之和保持常数。Roscoe 面价值在于给出了屈服面类型的一个选择依据。模型基于对临界状态线、 相关联塑性理论中屈服面与固结定律的假定。该模型假定 : 屈服只与应力球量 p 和应力偏量 q 两个应力分量有关,与第三应力不变量无关 ; 采用塑性体应变硬化规律 ,以ev 为硬化参数; 假定塑性变形符合相关联的流动法则 ,即g (s) = f ( s) ; 假定变形消耗的功 ,即塑性功为式中 :为塑性偏应变增量。由(1) 式得到的最初的剑桥模型屈服面形状为子弹头形,屈服函

5、数为后来提出了修正的假定式(3) 来代替 (1) 式,即在此假定的基础上 ,由(3) 式得到的修正剑桥模型屈服函数为椭圆,可以表示为其中 ,硬化函数 p c 为剑桥模型是当前在土力学领域内应用最广的模型之一,其主要特点有:基本概念明确 ;较好地适宜于正常固结粘土和弱超固结粘土 ;仅有 3 个参数 ,都可以通过常规三轴试验求出 ,在岩土工程实际工作中便于推广 ;考虑了岩土材料静水压力屈服特性、剪缩性和压硬性。2. 岩土临界状态下的变形特征所谓临界状态 ,即在常应力状态和常体积下岩土能够发生稳定的变形或流动 ;其孔隙比与应力状态之间存在着唯一的关系。临界状态首先是由Casagrade 提出 ,然后

6、由 Roscoe 等发展并应用于粘土的本构模型中。Been 等将临界状态应用于砂土本构模型中3 ,后来 ,Noveno 等4 建议将临界状态概念应用于所有的岩土材料(岩石、砂土、粘土和混凝土 ) 中。岩土材料力学特性的试验研究有大量的文献报道 ,对这些文献的研究分析表明 :岩土的应力 O应变曲线可大致归结为 2 种典型的形式 5 , 如图1 所示 ,曲线代表起始状态疏松的砂或正常固结粘土 ,曲线代表起始状态紧密的砂或超固结粘土。在图 1中,有2 个坐标系 ,即qO e1 坐标系和ev O e1 坐标系 , q为剪应力 , e1 为轴向应变 ,ev 为体积应变。对初始状态疏松的砂土或正常固结粘土

7、, 如应力 O应变曲线所示。在剪切过程中 ,随着轴向变形逐渐增加 ,剪应力逐渐增加 ,体积逐渐减小,即处于压缩状态 , 到达 B 点,随着轴向变形的继续增加 ,剪应力和体积应变逐渐稳定并接近一条水平线 ,即体积应变率为 0 ,应力率也为 0 。对初始状态紧密的砂土或超固结粘土 ,如应力 O应变曲线所示。在剪切过程开始段 ,剪应力较小、应变不大 ,整个试件体积略微压缩 ,表明这个阶段颗粒主要被挤向更紧密的排列 ,与侧限压缩试验的情况差不多 ;剪应力再增大之后 ,试件便进入屈服阶段 ,应变逐渐加快增大 ,试件体积开始膨胀 ( 或剪胀 ) ,即体积变形先是呈压缩状态后呈膨胀状态 ;剪应力到达某一峰值

8、 A 点后 ,由于颗粒与颗粒间的咬合作用逐渐丧失 ,剪应力无法继续上升 ,反而开始下降 ,剪应力和体积变化逐渐稳定并接近于一条水平线 , 即B 点, 此时轴向应变 e1 和体积应变 ev 很大 ,试件进入破坏阶段。试验曲线中的 A 点,常被称为峰值点 ,对应的抗剪能力称为峰值强度,曲线中的 B 点,常被称为临界状态点 ,对应的抗剪能力称为残余强度。从图 1 可以看出 :临界状态是岩土的体积保持不变或塑性体积应变率为 0的状态 , 通常对应的应变为 10 %或更大 ,对于应变硬化材料 ,它与破坏强度相对应 ;对于应变软化材料 ,它与残余强度相对应。3. 剑桥模型发展3. 1 基于经典塑性理论框架

9、的修正帽盖模型是 Dimaggi 和Sandier ( 1976)6 在剑桥模型的基础上提出的。该模型不仅能描述塑性屈服前的非线性、 剪胀性等特性 ,还能描述屈服后的各种破坏性状与塑性硬化性状。国内 , 魏汝龙 (1981)7 根据不排水三轴压缩试验资料得到的正常固结粘土模型 ,比修正剑桥模型具有更大的适应性 ,修正剑桥模型仅是它的特例。 杨林德、张向霞 (2005)8 针对 CamOclay 模型计算中存在没有充分考虑剪切变形的缺点 ,首先采用理论分析的手段 ,指出产生此现象的原因 , 同时提出由剪应力和体应力引起的应变分量分别采用不同的分配比例系数的思想。在此基础上对CamOclay 模型

10、进行了改进。3. 2 非关联流动Banerjee 和Stipho ( 1978 )9 基于剑桥模型框架 ,采用塑性增量理论,分析各向同性正常固结及微超固结粘土的不排水应力 O应变反映 ,建立了关联和非关联流动的弹塑性本构模型。 Newson (1998)10 提出一个软粘土非关联流动的临界状态本构模型 ,其屈服面与修正剑桥模型相同 , 塑性势函数为一经验函数 ,屈服面与塑性势面在偏平面上从低应力时的圆形到破坏时的松岗元屈服面。3. 3 次塑性理论与临界状态土力学Masin (2005)11 将次塑性理论与临界状态土力学相结合,建立了一个既适用于正常固结也适用于超固结粘土的本构模型。Masin(

11、2005)12 还利用次塑性模型预测修正剑桥模型的状态边界面的存在。3. 4 有限应变的剑桥模型Yatomi 等(1989)13 提出了考虑有限应变的剑桥模型 ,模型中加入了非共轴因素 ,可模拟土的局部剪切带。 Callari 和Auricchio ( 1998 )14 , Borja 和Ronaldo 等(1996 ,1998)15 ,16 也分别建立了有限应变剑桥本构模型 ,Ortiz 和Pandolfi (2004)17 基于剑桥模型框架 ,建立非粘性土本构模型 ,具有指数型压缩曲线 ,屈服面及塑性势为椭圆 ,采用体积应变硬化 , 模型预测从小应变扩展到有限应变。3. 5 时间相关的剑桥

12、模型Hsieh 等(1990)18 ,19 利用剑桥塑性理论将粘性土的应力 O应变关系加入与时间相关的因素 ,把总变形分为瞬时应变和粘性应变部分。 瞬时应变利用 2个屈服面来计算 ,一个屈服面为修正剑桥模型的椭圆形 ,另一个用 Von Mises 的圆柱形屈服面来描述 ,Arai 等(1988)20 ,Namikawa (2001)21 分别结合剑桥模型建立了依赖时间的正常固结粘土的塑性模型,Arai 的模型还考虑了粘土各向异性 K 0 固结对应力 O应变关系的影响。 Liao 等(2000)22 基于修正剑桥模型与 Perzyna 粘塑性方程 ,运用相关联与不相关联的塑性流动法则 ,建立了一

13、个反映时间效应和剪胀行为的本构模型。 Yin 和Graham (1999)23 基于修正的剑桥模型建立了一个弹粘塑本构模型 ,这个模型能模拟加速蠕变、 卸荷再加载、松驰等软粘土变形行为。3. 6 考虑各向异性及结构性日本学者太田和关口 (1979)24 提出了反映各向异性和应力轴旋转的关口 O太田模型 ,并被广泛应用于日本的工程实践中。其精华在于引进新的应力比 ,使剑桥模型成为其特例,能考虑 K 0 固结引起的应力各向异性,和主应力轴旋转产生的塑性变形。但不能考虑中主应力的影响 ,即不能考虑真正三维的应力状态。剑桥模型沿球应力轴 ( p 轴) 为等向塑性体变硬化 ,在日本广泛采用的关口 O太田

14、模型 24 沿初始固结线 ( K 0 线) 为不等向塑性体变硬化。三轴试验数据表明 :自 K 0 状态向伸长方向剪切时 ,前者方法计算的体积应变偏小 ,而后者方法计算的体积应变偏大。孙德安、姚仰平、殷宗泽(2000)25 提出一种介于上述两者之间、考虑初始应力各向异性(如K 0固结 ) 的不等向塑性体变硬化弹塑性模型。模型的剪切屈服准则使用SMP 准则 ,该模型能在三维应力下较好地反映土的强度和变形特性,模型的土性参数与剑桥模型一样。Voyiadjis 等( 2000)26 基于 Dafalias 的各向异性修正剑桥模型建立了一个新的考虑体土的各向异性和结构性的本构方程,利用了一个称为塑性旋转

15、张量的内变量,这个塑性旋转张量是引起各向异性的背应力的函数 , 模型增加了 2 个背应力参数。 Wheeler 等(2003)27 利用多阶段三轴排水试验的数据 ,改进修正剑桥模型 ,建立了饱和软粘土各向异性本构模型。Liu 和Carter (2002)28 将土的结构性加入修正剑桥模型,3 个新的参数用来描述土的结构性,如果所模拟的土无结构性,该模型等同于修正剑桥模型。 John 和Carter 等(1995)29 提出一个土的结构性剑桥模型,这个模型基于临界状态框架,能考虑天然粘土的当前应力状态、应力历史、当前孔隙比和当前土的结构,采用塑性体应变硬化。针对修正剑桥模型仅适用于正常固结及弱固

16、结粘土,而不适用于严重超固结土的情况 , Amerasin和Kraft (1983)30 ,Banerjee 等(1986)31 建立了严重超固结土的剑桥模型。 Mita 和Dasari 等(2004)32 基于严重超固结土的三轴压缩 O伸长试验及平面应变试验 ,建立了一个伏斯列夫 O 修正剑桥 ( HvorsleOMCC) 模型 , 模型具有伏斯列夫面 , 采用 MohrOCoulomb 屈服准则 ,并被扩展到一般三维应力空间。刘元雪 ,施建勇 (2003)33 从岩土类材料极限应力状态线所致的各向异性出发 ,提出了应力空间变换思想。 以修正剑桥模型屈服面的中心为映射中心 ,给出了重塑土的应

17、力空间变换、应力增量变换的公式 , 并考虑了应力洛德角的影响。 在变换应力空间中对修正剑桥模型进行了重新表述与改进 ,该模型可以较好地反映各向异性的影响 ,反映三轴伸长等应力路径的应力应变特性 ,也反映某些路径所致的软化现象。3. 7 循环荷载下的临界状态模型Carter 和Booker 等( 1979 )34 , Hirai 和Hi2royoshi (1987)35 分别建立了循环荷载下的临界状态模型,后者采用非关联流动法则,屈服面及塑性势均为剑桥模型形式。Nova 和Roberto ( 1983)36 基于剑桥模型 ,建立了模拟循环荷载下粘土或砂土变形特性的弹塑性本构模型,可模拟土的迟滞、

18、模量衰化及液化等行为。 Rouainia 和Wood (2000)37 提出了一个天然粘土的率相关本构模型,该模型是剑桥模型的扩展,建立在边界面塑性的运动硬化框架内,可考虑塑性应变引起的结构损伤。Borja 和Ronaldo等(2001)38 将经典修正剑桥模型与具有椭圆加载函数的各向异性边界面塑性模型结合,该模型还耦合非线性超弹性模型。3. 8 扩展到一般三维应力空间Nakai 等(1986)39 基于 SMP 屈服准则 ,建立了一个在三维应力空间模拟粘土各种不同应力路径下行为的本构模型,模型不仅考虑了中主应力的影响 ,也反映了应力路径对土体应力与应变关系的影响。模型参数确定方法与剑桥模型相

19、同。 Matsuoka 等(1999 ,2005)40 ,41 采用 SMP (空间滑动面 ) 准则对剑桥模型进行了扩展 ,将剑桥模型从三轴试验的轴对称应力状态扩展至一般应力状态 ,其修正的方法是引入变换应力 ,重新定义变换后的广义正应力和剪应力 ,然后代入剑桥粘土模型中 ,实现从 Mises 破坏准则向 SMP 破坏准则的扩展。 姚仰平等 (2000)42 将Lade 屈服准则引入最初的剑桥模型 ,以实现模型的三维化 ,采用变换应力将 Lade 屈服准则变为扩展的 Mises 屈服函数类型。孙德安 ,姚仰平 (2002)43 又提出了一个适用于粒状材料的简单而实用的弹塑性模型。 该模型通过对

20、修正剑桥模型进行改进 ,能够反映三维应力状态下的剪胀剪缩性和变形、强度的平均应力依存性。 姚仰平等 (2006)44 还把 Mohr O Coulomb 准则与剑桥模型相结合 ,提出了一个变换应力张量 ,在主应力空间上 ,将Mohr O Coulomb 准则变为圆锥形 ,能描述一般应力空间土的变形特性。3. 9 扩展到砂土、非饱和土Borja 和Ronaldo45 ,将经典修正剑桥模型从饱和粘土扩展到非饱和土。 Robinet 等46 建立了适用于膨胀土 ,具有两个屈服面的修正剑桥模型形式的本构方程。 砂土的变形行为不仅依赖于其相对密度也依赖于所处的固结压力。同一固结压力下,密砂剪胀 ,松砂剪

21、缩 ;同一密度的砂 ,在低固结压力下可能剪胀,而在高固结压力下则可能剪缩。 在临界状态理论的框架内 ,为了统一反映密度和压力对砂土变形特性的耦合影响, Been和Jefferies47 建议了一个状态参数 ,把围压和密度的耦合作用加入到本构方程中 ,其中 , Yu48 提出一个对粘土和砂土都适用的统一临界状态模型,采用应力比与状态参数的关系来描述土的状态边界面 ,重新定义了屈服函数及塑性势函数 ,模型适用范围广 ,相对简单 ,与剑桥模型相比 ,仅增加2 个参数。 Gajo 和Wood49 基于剑桥临界状态概念 ,小应变弹性区假定,状态参数概念 ,剑桥模型的流动法则 ,边界面概念和运动硬化塑性 ,建立了一个反映较大范围孔隙比和应力水平的粒状材料本构模型。4. 剑桥模型的局限性纵观剑桥模型 40 多年的发展 ,总结其局限性主要有 :(1) 受制于经典塑性理论 ,采用 Drucker 公设和相关联的流动法则 ,在很多情况下与岩土工程实际状态不符 ;破坏面有尖角 ,该点的塑性应变方向不易确定。(2) 因为屈服面只是塑性体积应变