剑桥模型推导

《剑桥模型推导》由会员分享，可在线阅读，更多相关《剑桥模型推导（12页珍藏版）》请在装配图网上搜索。

1、比容的定义:正常固结线(NCL)方程：v=N-Inp临界状态线(CSL)方程：v=T-Inpf(3)回弹线(SL-swellline)方程：v=v-kInpf(4)K注意：在lnp-v平面上，回弹线SL尽管穿过了CSL线，但并不意味等压卸载过程中应力点曾达到CSL线上，因为此坐标系中CSL为空间CSL曲线的投影，而SL始终在lnp-v平面上,并不能达到空间的CSL线上的应力状态。CSLCSL无拉力墙0CSLPHvorslev面Roscoe面NCL图1土的物态全界面v无拉力墙Hvorslev面CSLRoscoe面归一化后土的物态全界面NCL和艸切碗三向图？4正常固结土的胸态边界面圈厶35弹性墙上

2、费性比护积严惟一性和剑桥模型的區服面在上图2-34中AR为卸载回弹线（其方程如式（4）,过其作的竖直曲面，此曲面位于物态全界面（Roscoe面、Hvorslev面及无拉力墙构成）以下的阴影部分，即为一弹性墙，此弹性墙交物态边界面Roscoe面于AF，在AR线上荷载变化时，无塑性体积变化，亦即在弹性墙上，塑性体应变昇呆持为常数。如果选择塑性体应变为硬化参数，那么等塑性体应变面就V是屈服面，等塑性体应变线AF就是屈服轨迹。AF在p-q平面上的投影AF为屈服面在p-q平面上的屈服轨迹。在图2-35中回弹曲线与比容轴截距代表其塑性比容vp，在同一弹性墙上,0(8)(9)(10)(11)(12)(13)

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

4、vsdWppdp,qdp(7)vs关于弹塑性变形能，Roscoe作了如下的假设：(1)假定一切剪切应变都是不可恢复的,亦即无弹性剪应变,只有不可恢复的塑性剪应变(总剪应变等于塑性剪应变)de0sdpdss(2)假定弹性体应变可从各向等压固结试验中所得的回弹曲线求取,即由式(4)可得dPe型pdvedpUev1+e1+epdWe=pdedpp1,e故:dpddedvvvv1+ep(3)假定全部耗散能(塑性变形能)等于由摩擦产生的能量耗散,即:dWp卩pdpMpdpss式中卩为内摩擦系数，其值等于p-q平面上临界状态线CSL的斜率M(14)于是pd&qdsvsdpMpdsp1es(19)dp八1e

5、PlsM-6血(三轴压缩)(15)3一sin或6sinM(三轴伸长)(16)3+sin所以dW，dWedWp，dpMpdsp1es(17)而单位体积的土在p,q应力作用下如产生应变ds和ds,变形能为vsdW，pds+qdsvs(18)则由式(17)和式(18)可得能量方程:将式(13)代入上式,则pdsp=Mp一qJdspvs(20)dspq或-，M-，M一耳dspp式(20)实际表示了流动法则，即表示了塑性应变增量在p-q平面上的方向，与这一方向正交的轨迹就是在这个平面上土的屈服轨迹(相适应的流动法则),如图2-34所示.设此屈服轨迹的方程为:f(p,q,h)=o(21)则df-dfdp+

6、dfdq+苗dH-0诃dqfdHlint-f-丁卞l/、.、rr、小rt“(22)因为在冋一丿出服面上硬化参数为吊数，所以dH，0,则df需,如pq,dqq,0(23)根据相适应流动法则dsp，d色(24)”dps(25)(26)(27)(28)将以上两式代入式(23),则得dpdp+dqdp=0vs将式(20)代入上式,则得将此微分方程变换可得到-Mqfdpf+Mpfdqfdp=(Mp)2pf积分得到+Inp=CMpf式中C为积分常数.利用p轴上起始各向等压固结试验点A,对应p=p,q=0代入上式,0则得C=lnp,将之代入式(28),则得得到湿粘土(正常固结和轻超固结土)的屈服轨迹0方程为

7、f=q-MInp=0(29)pp其在p-q平面上的形状如图2-34和图2-35(a)所示，像一个”帽子”，是子弹头形，以p为硬化0参数.由于NCL上每一个p都对应于一个Vp(或p),所以实际上这一模型是以塑性体应变p00vv为硬化参数.对于重超固结土，可得到类似的屈服面，只是对应的p不同.0空间无拉力墙的方程为Mhvq=3p0pexp(30)3h(九丿Hvoslev面的方程为rvq=hp+(Mh)exp()(31)九式中h为Hvoslev线的斜率.空间Roscoe面的方程为:q=(Nv-九Inp)(32)入一K由湿粘土对应p=p,v=v的不排水试验路径在p-q平面上的投影或归一化的Roscoe

8、面，00由式(2)得v=v=N一,Inpf(33)00将式(33)代入式(32)，则得对应不排试验路径在p-q平面上的方程为半In4-0(34)p，一p也为指弹头形，但显然此不排水路径与屈服轨迹并不重合，不排水路径在屈服轨迹以外.剑桥模型增量型应力-应变本构关系将式(32)微分，可得,(dq_ndp，)+dpMpp因由式(11)知dv=(1+e)dvdv=-所以,1-K/,1+eMp,-+e)Mp,一ndp+dqq”,一K又因n=qpdn=dq_-ndppp于是ds=丄v1+edn+，璽Mp将式(38)代入能量方程(19)，可得,-Kd&=s(1+e)Mpdp+Mdqm_n于是剑桥模型的弹塑性

9、矩阵可表示为:,-+e)Mpdpdq(35)(36)(37)(38)(39)(40)M一耳修正剑桥模型:1965年,英国剑桥大学的Burland采用了一种新的能量方程形式，得到了修正剑桥模型.他建议以下式代替式(14)即假定总的塑性变形能等于塑性体变能和由摩擦耗散能的算术平方根,以之代替式(19)右边第二项,则2、2s丿+M2-qp故可得:d,p_M2-mvd,p2qs此即修正剑桥模型的流动法则.将其代入式(26)，得到在p-q平面上的屈服轨迹方程为或pM2pm2+mc或(p，)2+()2M2或:p一p0/212+q0丿0Mp/2丿021(42)(43a)(43b)(43c)(43b)即为椭圆

10、方程.其顶点在qMpf线上，以pr(,p)为硬化参数，即0vH(p)pp(,p)H(,p)000vv因为pexp0其增量型应力-应变关系为d，二丄v1+e(K)勿內九型M2+mp2,d,dp、_K)2+,2)pfd8丄厂d8(1+s(M2+,2)2,(_k)耳2,M2-,2dpdq+1+eM221M2+,2p?于是修正剑桥模型的弹塑性矩阵可表示为：然而,有限元等数值计算中,常按如下一般弹塑性矩阵式de=妙0勝阍-+楼来表示,由等塑性硬化规律帥dH別帥dp別A=0dHd8pdpdpd8pdpv0vH=H(8p)二pv0有：A=ddHdQdHd8pdpv按相适应的流动法则(p)(1+e)(2pp)1+eA0exp(k丿ddHdddpd0”dHd8pdpdpd8pdpv0vMi+eaap(-K)i