剑桥模型推导讲课稿

1、精品文档比容的定义：v=vl=.m_A =1.e正常固结线(NCL )方程：v = N - £ ln p'(2)临界状态线(CSL)方程：v=，Jnp(3)回弹线(SLswell line )方程：v =十一父In p'(4)注意：在lnp'-v平面上，回弹线SL尽管穿过了 CSL线，但并不意味等压卸载过程中应力点 曾达到CSL线上，因为此坐标系中CSL为空间CSL曲线的投影，而SL始终在lnp'-v平面上, 并不能达到空间的CSL线上的应力状态。q*iv图1 土的物态全界面无拉力墙图334正常固结土的物态边界面R2-35弹性墙上塑性比体积/惟一性 和

2、剑桥模型的屈服面在上图2-34中AR为卸载回弹线（具方程如式（4）,过其作的竖直曲面，此曲面位于 物态全界面（Roscoe面、Hvorslev面及无拉力墙构成）以下的阴影部分，即为一弹性墙，此弹 性墙交物态边界面Roscoe面于AF,在AR线上荷载变化时，无塑性体积变化，亦即在弹性 墙上，塑性体应变 寸保持为常数。如果选择塑性体应变为硬化参数，那么等塑性体应变面就 是屈服面，等塑性体应变线AF就是屈服轨迹。AF在p'-q'平面上的投影A'F'为屈服面在p'-q' 平面上的屈服轨迹。在图2-35中回弹曲线与比容轴截距代表其塑性比容 v°p

3、,在同一弹性墙上，或同一屈服线上，弹性墙的塑性比容 vp=v0p =const,也就是说其塑性体应变 4为常数。剑桥模型基于传统塑性位势理论，采用单屈服面和相关联流动法则。屈服面形式 （方程）A'F' 不是基于试验而提出的，上面已根据物理意义在几何上表示出屈服面 A'F',但还无法用数学 表达式表示，剑桥模型是依据能量理论得出的其屈服面方程，实质上是一种假设。依据能量方程，外力（荷载）做功dW一部分转化为变形体的弹性变形能dWe（可储存在变 形体内，外力或荷载卸除时，可完全释放出来），另一部分转化为耗散能（或称塑性变形能， 外力或荷载卸除时，不能再释放出来）dW

4、p ,因而有dW =dW e dWp（5）两种变形能可表示如下：(6)dWe = pd qd e关于弹塑性变形能，Roscoe作了如下的假设:（1）假定一切剪切应变都是不可恢复的，亦即无弹性剪应变，只有不可恢复的塑性剪应变（总 剪应变等于塑性剪应变）(8)d sp = d s（2）假定弹性体应变可从各向等压固2试验中所得的回弹曲线求取 ，即由式（4）可得(10)(11)(12),e dp dv F pdvedpdWe = pd ；： =dp1 - e故:ped v 二d v d v =d vdp(13)假定全部耗散能（塑性变形能）等于由摩擦产生的能量耗散，即:(14)式中 N为内摩擦系数，其值

5、等于p'-q'平面上临界状态线CSL的斜率M精品文档所以6si n（三轴压缩）(15)6si n M =-3 sin（三轴伸长）(16)dW-dWe dWP(17)而单位体积的土在p',q'应力作用下如产生应变d8和d/,变形能为dW=pd v qd s(18)则由式（17）和式（18）可得能量方程:于是p d v q ds =1 edp M'#sd(19)加 时dpp dsv -二11+e p J将式（13）代入上式，则或 d vp - m - - = M -d sp(20)式（20）实际表示了流动法则，即表示了塑性应变增量在p'-q'

6、平面上白方向，与这一方向正交的轨迹就是在这个平面上土的屈服轨迹（相适应的流动法则），如图2-34所示.设此屈服轨迹的方程为:f P,q,H =0(21)f；f；fdf = dp -dq dH = 0pqH(22)因为在同一屈服面上硬化参数为常数，所以dH =0,则开开df dP dq =0pq(23)根据相适应流动法则(24)d sp = d -：q将以上两式代入式（23），则得dpd vp dqd sp = 0将式（20）代入上式，则得dq qM =0dp P将此微分方程变换可得到Mq dp Mp dq dp n2 二 0（Mp）2 p积分得到-q- ln p =CMp(25)(26)(27

7、)(28)式中 C为积分常数.利用p'轴上起始各向等压固结试验点 A,对应p'= p0, q' = 0代入上式,则得方程为C =ln p0,将之代入式（28）,则得得到湿粘土（正常固结和轻超固结土）的屈服轨迹(29)其在p'-q'平面上的形状如图2-34和图2-35（a）所示，像一个“帽子"，是子弹头形，以p0为硬化参数.由于NCL上每一个p0都对应于一个v0P （或个）,所以实际上这一模型是以塑性体应变 1P为硬化参数.对于重超固结土，可得到类似的屈服面，只是对应的p0不同.空间无拉力墙的方程为,M h V q' = 3p'

8、0 <pf<M-hexpV3-hi 九，(30)Hvoslev面的方程为q = h p ( M-)hex)(31)式中 h为Hvoslev线的斜率.空间Roscoe面的方程为：(32)Mp , q ( N - v l rp ) 九一£由湿粘土对应p'= pO, v = v0的不排水试验路径在p'-q'平面上的投影或归一化的 Roscoe面， 由式(2)得v = v0 = N - ln p0(33)将式(33)代入式(32)，则得对应不排试验路径在p'-q'平面上的方程为q_lnE=0(34)p /. ： p也为指弹头形，但显然此不排

9、水路径与屈服轨迹并不重合，不排水路径在屈服轨迹以外.剑桥模型增量型应力-应变本构关系将式(32)微分，可得dv = -dq - dp dp _Mpp(35)因由式(11)知 dv= T1 +e cv所以九 1父/九 zdp(dq -ndp )+-1 e Mp If - dp dq(36)又因二qF pdq(37)于是1 e IL Md dp(38)将式(38)代入能量方程(19),可得(39)九一K I d、= TTIMp dp于是剑桥模型的弹塑性矩阵可表示为d v _1 - ''ds 一 1 e Mpdp(40)修正剑桥模型:1965年，英国剑桥大学的Burland采用了一种

10、新的能量方程形式，得到了修正剑桥模型.dWp他建议以下式代替式(14)(41)精品文档，以之代替式（19）右边第即假定总的塑性变形能等于塑性体变能和由摩擦耗放能的算术平方根 二项，则cn 2c 2cpd V3 = pd vpMpd .qd ；：1精品文档d ,pd7.p 5 d<d ®p sM2 -q- pd%p "2故可得:M 2 - 2-2此即修正剑桥模型的流动法则.将其代入式（26），得到22dq M - n0dp 2在p'-q'平面上的屈服轨迹方程为p M2pC =M22p -p。/ 2po / 2IMp0 / 2 )=1(42)(43a)(4

11、3b)(43c)(43b)即为椭圆方程.其顶点在q'=Mp'线上，以p0（）为硬化参数,即H(Po) = P0 =p0( J)=H( ；J)因为1佳=exp 二eo其增量型应力-应变关系为M 2+M p'.精品文档精品文档2s= TT %FLr 2" d"dp方方"才"P7于是修正剑桥模型的弹塑性矩阵可表示为d v d s2 二）1 e M22 p九（M 2+%n2（Z-k ）-乜仃= De1-Q I f IT eDe ;Q >DeA cgDe Q- /ds=De 一 Dpd)_ |Dep id ；)来表示，由等塑性硬化规律有： A =;:H :Q按相适应的流动法则H :- P0A =二-r -7一H 二 vP PPo 二 J