下载文档原格式
1.剑桥模型（Cam-clay Model ）剑桥模型是由英国剑桥大学Roscoe 等于1963年提出的，这个模型基于正常固结土和超固结土试样的排水和不排水三轴实验基础上，提出了土体临界状态的概念，并在实验基础上，再引入加工硬化原理和能量方程，提出剑桥模型。
这个模型从试验和理论上较好的阐明了土体弹塑性变形特征，尤其考虑了土的塑性体积变形，因而一般认为，剑桥模型的问世，标志着土本构理论发展的新阶段的开始。
（1） 剑桥模型。
剑桥模型基于传统塑性位势理论，采用单屈服面和关联流动法则屈服面形式也不是基于大量的实验而提出的假设，而是依据能量理论提出的。
依据能量方程，外力做功dW 一部分转化为弹性能e dW ，另一部分转化为耗散能（或称塑性能）p dW ,因而有dW =e dW +p dW（1-154） e dW =e eV qd d p γε+' （1-155）p pV p qd d p dW γε+'= （1-156）剑桥模型中，由各向等压固结实验中回弹曲线确定弹性体积变形p p d e k d eV ''+=1ε （1-157）式中，k 为膨胀指数，即 p In e '-回弹曲线的斜率。
同时，假设弹性剪切变形为零，即0=e d γ （1-158）则弹性能 p d ekp p d k dW e '+=''=1υ （1-159）剑桥模型中还建立如下的能量方程，即塑性能等于由于摩擦产生的能量耗散，则有 p p p V d p qd d p γνγε'=+'- （1-160）式中第一项改用负号，是因为p V d ε取以压为正。
代入式（1-161）⎪⎭⎪⎬⎫==ij pij p d s d d λεεθθσ （1-161） 并考虑式（1-158），则有γγγνd p M d p M d p dW p p p '='='= （1-162） 式中，M 为q p '-'平面上的破坏线的斜率，即ϕϕ'-'=sin 3sin 6M （1-163） 式中ϕ'为土体有效摩擦角。
土的剑桥模型发展综述----97dbe56c-6ead-11ec-b680-7cb59b590d7d土体本构理论是岩土工程学科的重要基础理论。
随着对土体力学特性的不断深入,塑性理论逐渐被应用于土体本构关系的研究中来。
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）平面上的正常固结线，在该平面上发生剪切变形而不发生体积变形。
土的弹塑性模型近年来，根据弹塑性理论建立的土的弹塑性模型发展很快，各国学者提出的弹塑性本构模型很多。
下面几节分别介绍剑桥模型，修正剑桥模型，Lade-Duncan 模型，以及清华模型的基本概念。
一．剑桥模型英国剑桥大学Roscoc 和他的同事（1958 ~ 1963 ）在正常固结粘土和超固结粘土试样的排水和不排水三轴试验的基础上，发展了Rendulic （1937）提出的饱和粘土有效应力和孔隙比成唯一关系的概念，提出完全状态边界面的思想。
他们假定土体是加工硬化材料，服从相关联流动规则，根据能量方程，建立剑桥模型。
剑桥模型从理论上阐明了土体弹塑性的变形特性，标志着土的本构理论发展新阶段的开始。
1.临界状态线和Roscoe 面各向等压固结过程中，孔隙比e 或比容()1e υυ=+与有效应力的关系可用下式表示： ln N p υλ'=- （1） 式中 N —— 当 1.0p '=时的比容。
因此exp N p υλ-⎛⎫'= ⎪⎝⎭（2）（a ）,p q ''平面（b ）,ln p υ'平面图1 临界状态线正常固结粘土排水和不排水三轴试验表明：它们有条共同的破坏轨迹，与排水条件无关。
破坏轨迹在,p q ''平面上是一条过原点的直线，在,ln p υ'平面上也是直线，目与正常固结线平行，分别如图（a ）和（b 〕 所示。
破坏轨迹线可用下式表示：cs csq Mp '= （3）ln cs csp υλ'=Γ- （4） 式中 CS ——表示临界状态；M——,p q''平面上临界状态线斜率；Γ—— 1.0p'=时土体的比容；csυ'平面上临界状态线斜率。
λ——,ln p一旦土体的应力路径到达这条线，土体就会发生塑性流动。
这时土体被认为处于临界状态，破坏轨迹被称为临界状态线。
临界状态线在,,''空间为一条空间曲线，如下图2所示。
基于ABAQUS二次开发状态相关砂土本构模型的研究司海宝;杨为民;黄伟【摘要】剪胀性是砂土最重要的特性之一,其变形主要取决砂土的当前状态,而基于状态相关砂土本构模型能较好反映砂土剪胀性.ABAQUS作为有限元计算平台,具有强大的非线性计算和前后处理功能,标准ABAQUS程序仅包含一些通用土体本构模型,缺少状态相关砂土本构模型,因而有必要针对土木工程具体问题进行二次开发.利用二次开发工具和UMAT数据接口程序,采用隐式积分算法将基于状态相关砂土本构模型嵌入到标准ABAQUS程序中,并将有限元数值计算结果与实验数据及理论计算值比较,进行单元敏感性分析.计算结果表明子程序运行稳定可靠,不仅充分利用ABAQUS方便、快捷的前后处理和强大的非线性求解平台,还可以完成更复杂,且切合实际土体应力应变关系的有限元数值计算,拓展ABAQUS在土木工程领域应用能力.【期刊名称】《安徽工业大学学报（自然科学版）》【年(卷),期】2013(030)003【总页数】7页(P301-307)【关键词】状态相关砂土模型;ABAQUS;隐式积分算法;二次开发【作者】司海宝;杨为民;黄伟【作者单位】安徽工业大学建筑工程学院,安徽马鞍山 243032;港口航道重点实验室南京水利科学研究院,江苏南京 210045;中国地质科学院地质力学研究所,北京100081;安徽工业大学建筑工程学院,安徽马鞍山 243032【正文语种】中文【中图分类】TU476剪胀性是砂土最重要的特性之一，砂土在剪切条件下除了发生剪切变形外，还会产生体积变形，砂土的体积变形可能是膨胀，也可能是压缩，主要取决于砂土的当前状态[1]。
Li和Dafalias在三轴压缩空间内建立弹塑性本构模型，进一步拓展了状态相关剪胀理论，较好模拟了砂土在各种状态下变形特征[2-3]。
ABAQUS是目前土木工程中应用较为广泛的非线性有限元分析软件，具有很强的非线性计算功能和前、后处理能力，同时其强大的二次开发功能使得其在工程领域具有良好的开放性。
10分钟认识剑桥模型王川第一节：认识“临界状态”首先，大家一定接受以下两张图（无数实验已经证明过）：图1 摩尔库伦强度理论图2 土的压实曲线（e为孔隙比，p’为有效应力）那么，如果把τ换成偏应力q（其中q=σ1-σ3），把σ换成平均主应力p（其中p=(σ1+2σ3)/3，p’表示其有效应力），就得到三轴实验中的p-q曲线：图3 p-q曲线土样的体积由固体颗粒和空隙组成，由于固体颗粒不可压缩，故土样体积的变化完全取决于空隙的变化，即土样体积v和孔隙率e描述的物理意义等价。
那么，将图2中e替换为v，就得到v-logp曲线：图4 v-logp曲线与图1和图2一样，图3和图4同样经历了无数实验的验证，属于“事实”。
基于图3和图4的定量分析以及实验观察，可以得出一个结论，这个结论就是临界状态（critical state）：无论土样的初始状态和经历的应力路径如何，在剪切的最终阶段，只有剪应变在持续增加，而土样所受的有效应力和体积趋于不变。
临界状态由图3和图4同时确定，因此图3和图4中的曲线也叫临界状态线CSL （Critical State Line）。
将临界状态现象翻译成数学语言：（1）体积不变对应于，为p’引起的体积的改变；（2）剪应变在变对应于，为q引起的剪应变；（3）有效应力不变等价于q与p’的比值为常量。
若令在一般情况下，有（被叫做应力比），则可以定义临界状态下的应力比：（被叫做临界状态应力比）。
从图3中能看出，M为常量，即“有效应力不变”。
◆第二节：剑桥模型假设（1）所有的剪应变都不可恢复，即（为弹性剪应变），（为塑性剪应变）。
（2）假定塑性变性能增量可表示为：（这一假设看不懂没关系，继续往后看）。
（3）相关联流动法则：（与塑性力学中关联流动一致）。
◆第三节：剑桥模型推导从能量角度推导屈服函数：应变能的增量等于主应力p’和偏应力q所做的功，即（式1）因为：（此处用了假设1）所以：（式2）（此处用了假设2）由式1和式2得：（式3）根据（假设1），整理式3得：d为剪胀系数，表示塑性应变的方向（因为d体现了与的相对大小，与塑性力学中流动法则表达得意义一致）；为剪胀方程。
