土体本构模型

1、土土 体体 本本 构构 模模 型型本构关系：材料的应力本构关系：材料的应力应变应变(时间时间)关系关系本构模型：反映材料的应力应变本构模型：反映材料的应力应变( (时间时间) )关系的数关系的数学模型，即数学表达式。当然，这种数学表达式可能学模型，即数学表达式。当然，这种数学表达式可能很复杂，而且包括一系列的数学表达式。很复杂，而且包括一系列的数学表达式。E1. 1.应力和应变应力和应变 Txzyzxyzyxppp(一）应力和应变分量的几种表示方法一）应力和应变分量的几种表示方法1.一般分量一般分量 土体中一点的应力状态，可以用处于该点的正六面体单元的表面上的6个（9个）应力分量来表示，即3个

2、正应力分量 ，3个剪应力分量 写成矩阵形式为zyx,zxyzxy, Tzxyzxyzyx应力偏量（1）矩阵或向量表示法偏应力1. 1.应力和应变应力和应变(一）应力和应变分量的几种表示方法一）应力和应变分量的几种表示方法1.一般分量图51中表示了单元体上的这6个应力分量。相应地，也有6个应变分量，以矩阵表示为Tzxyzxyzyx（1）矩阵或向量表示法1.1.应力和应变应力和应变图5-1（2）张量表示法如果某些量依赖于坐标轴的选取，并且，当座标变换时，它们的变换具有某种指定的形式，则这些量总称为张量。一点的应力分量就总称为应力张量。ijijijzzyzxyzyyxxzxyx在进行公式推导时，一般

3、尽量用一种表示方法：矩阵或张量，不宜混用1.1.应力和应变应力和应变（2）张量表示法ijSpppzzyzxyzyyxxzxyx注意：在进行公式推导时，一般尽量用一种表示方法：矩阵或张量，不宜混用偏应力张量 在弹性力学中，法向应力和应变以拉为正，压为负；而土体一般不能受拉，土力学中讨论的地基应力、土压力等，都是以压为正，拉为负。因此，土力学中，应力应变分量的正负规定就与弹性力学相反，即正面上的负向应力为正，负面上的正向应力为正。不仅正应力如此，剪应力也如此，以保持一致，并能套用弹性力学公式。 注意：注意： 1.1.应力和应变应力和应变2. 主应力应变分量主应力应变分量 在正六面体单元中可以找到3

4、个互相垂直的面，其上剪应力为0，只作用有正应力。这样的面叫正应力面，所作用的正应力叫主应力。3个面上的主应力按大小排列，分别为大主应力 、中主应力 和小主应力 。1231.1.应力和应变应力和应变应力状态表示方法之一：主应力方向余弦主应力与坐标轴的选择无关，应变也可用三个主应变分量表示，矩阵形式 Ti321 将坐标系的三个轴顺着三个主应力方向放，分别以1，2，3表示，如图5-2所示。再对这个坐标系的8个挂限分别作等倾面。8个挂限的等倾面围成了一个正八面体。这些等倾面叫八面体面。根据力的平衡关系可以推得正八面体面上的正应力和剪应力分别为 3. 八面体应力和应变八面体应力和应变图5-21.1.应力

5、和应变应力和应变32131OCT21323222131OCT21323222131OCT 剪应力OCT作用在八面体面上还有个方向问题。这决定于中主应力 接近大主应力 还是小主应力 。 与应力相应，还有八面体面上的应变，正应变和剪应变分别为 21332131OCT21323222132OCT1.1.应力和应变应力和应变 在土体本构模型理论中，常常也用球应力、偏应力以及或作为应力分量。球应力也称为平均正应力以p表示4. 球应力、偏应力及相应应变球应力、偏应力及相应应变32131p偏应力又叫广义剪应力，以q表示21323222121q1.1.应力和应变应力和应变注意：这里的偏应力和Sij的区别，建议

6、这里不用偏应力)(31zyx注意：此式也可用6个应力分量表示p 和q也可以用八面体应力来表示，如下p OCTOCTq23 反映了复杂应力状态下受剪的程度，因此常用来表示剪应力。当 时，如轴对称的三轴仪试样受力情况， 32311.1.应力和应变应力和应变可以推知相应的应变分量体积应变：321v偏应变： 其中 表示了复杂受力状态下的剪切变形。对于轴对称三轴试样的变形，有 s32312v332131vs1.1.应力和应变应力和应变231232221)()()(32s 球应力和偏应力，以及相应的应变分量，实际上与八面体应力和应变是等效的，仅仅是系数不同。但在分析能量时，要简单得多。可以推得： 体积变形

7、能 ：vvpW形变能 ：ssqW1.1.应力和应变应力和应变 对于一组确定的p和q，可以有许多种主应力分量的组合，解是不确定的。因此，要有第三个分量。第三个分量常取应力罗德（Lode）参数 313123121322 式中 ， ， 为三个应力摩尔圆的直径，见图5-3 3221311.1.应力和应变应力和应变213132b还有一个参数b也反映了中主应力接近大主应力的程度。若 ，b=1；若 ,b=01232图5-3相应地，也有应变罗德参数3131221.1.应力和应变应力和应变5. 应力不变量应力不变量第一应力不变量： 3211I第二应力不变量：1332212I第三应力不变量：3213I 此外，还有

8、下面两个偏应力不变量，它们须与第一应力不变量相结合形成三个独立的应力分量：1.1.应力和应变应力和应变第二偏应力不变量 ：第三偏应力不变量 ：213232221612J2131323212713222J不随坐标轴的选取而改变表示一点应力状态的方法？（二）应力空间和应力路径1应力和应变空间 为了表示应力状态，表示各应力分量的数值，常常以应力分量为坐标轴形成一个空间，叫做应力空间。该空间内的一点的几个坐标值就是相应的应力分量。 如果应力分量取三个主应力 ， 和 ，以三个主应力分量为坐标轴构成一个直角坐标系，叫主应力空间。这个空间内一点有三个坐标值，就代表了实际土体中一点的某种应力状态。图5-4中的

9、M点代表了应力状态 ， 和 。123M1M2M31.1.应力和应变应力和应变图5-41.1.应力和应变应力和应变弹塑性力学：Pi平面为过原点与空间主对角线垂直的平面平面 在主应力空间内，法线与空间主对角线重合的等倾面，被叫做 面。所谓空间主对角线，就是与3个坐标轴的夹角都相等的线。主应力空间中，在该线上有 321 八面体面是几何空间（长度坐标系）内的面，面是在应力空间内的面。两者坐标系不同，物理概念不同。再者，八面体面在几何空间内的八个挂限都有，而 面只存在于应力空间内的第一挂限和与其相对的挂限，其它挂限内的等倾面并不是面。空间主对角线也只存在于这两个挂限。1.1.应力和应变应力和应变 利用面

10、可以较好地反映应力状态。图5-4中点M的坐标代表主应力分量。通过M点作面。它到原点的距离为pOO331321在面上，M到空间主对角线的距离OM21323222131q32 它们分别与应力分量p和q有关。而点M在 面内的方位可反映第三个分量。将图5-4中的三个主应力坐标轴，以及代表应力状态的点M 投影到 面上，如图5-5所示。 1.1.应力和应变应力和应变 在该面上放一个二维的直角坐标系，令Y轴与2轴重合，X轴在1 的那一侧。定义到X轴的转角叫应力罗德角。它就是与第三应力分量有关的参数。可以证明，它与罗德参数间的关系为：3tan1.1.应力和应变应力和应变应力空间还可以用其他形式的应力分量为坐标

11、。如果以，和六个应力分量为坐标，则应力空间是六维空间，无法用图形表示，仅可以作抽象的理解。 p-q 平面1.1.应力和应变应力和应变 如果忽略第三应力不变量或应力罗德角对变形的影响，可以只用、两个分量来构成二维的应力空间，叫pq平面，如图5-6所示。在后面的本构模型理论中，常常会用到这种平面。图5-61.1.应力和应变应力和应变图5-61.1.应力和应变应力和应变表示应力状态或应力路径也有优点P204qp )(2131p)(2131q二维问题中， 与应力空间相应，以应变分量为坐标轴形成一个空间，叫做应变空间。该空间内的一点的几个坐标值就是应变分量。图5-8所示为主应变空间。它的三个坐标轴分别为

12、 ， 和 。123图5-81.1.应力和应变应力和应变2.应力路径 在应力空间内，代表应力状态的点移动的轨迹，叫应力路径。它表示应力变化的过程，或者加荷的方式。1.1.应力和应变应力和应变图91.1.应力和应变应力和应变设土体中一点初始应力状态如图5-9应力空间内点所示，受力后变化到。从到，可以有各种方式，如、和按比例增加；初期增加得多，和增加得少，而后期反过来。对于某种加荷方式，代表应力状态的点将从沿某种轨迹移动到。加荷过程中，不同的加荷方式可以用不同的应力路径来表示。 更常用的是用p-q平面的应力路径1.1.应力和应变应力和应变与其相应，当然也有应变路径。普通三轴应力状态下pq（三）应力应

13、变关系矩阵kD D广义虎克定律增量形式)()()(yxzzzxyyzyxxEEEEEE（三）应力应变关系矩阵广义虎克定律 DDDD（三）应力应变关系矩阵 复杂应力状态下的应力应变关系是多元化的，要表示出多元素与多元素之间的关系，就要用张量或矩阵。常用到的增量形式的应力应变关系的矩阵为D 式中D叫刚度矩阵，如果应力和应变分量取一般形式，各有6个分量，则矩阵D为66，共36个元素。如果用主应力和主应变分量，则矩阵D为33，共9个元素。二维问题的应力分量为 ，应变分量为 ，因此其矩阵D也是3 3 的，将上式展开可写成： xyyx,xyyx,1.1.应力和应变应力和应变xyyxxyyxDDDDDDDD

14、D333231232221131211 对于任一元素D i j，其意义为，要产生单位应变增量 而其它应变增量为0时，在应施加的应力增量 中的分量 即为Dij。显然，D i j 的值愈大，材料愈难变形，表示材料刚度愈大。ji1.1.应力和应变应力和应变xyyxxDDD131211 应力应变关系也可写成相反的形式，即: C 式中C叫柔度矩阵。对于二维问题，将其展开，可写成xyyxxyyxCCCCCCCCC333231232221131211 在复杂受力条件下，建立土的应力应变关系，实际上就是要给出矩阵C或D。 C或D互为逆矩阵1.1.应力和应变应力和应变Cij越大，材料越软xyyxxCCC1312

15、112.2.土体三维变形的试验土体三维变形的试验三轴仪应力变形试验 三轴仪的构造示意如图5-10所示。图5-10仪器构造： 中间为圆柱形土样。其下为透水石，透水石放在三轴仪底座上；试样顶部也放有透水石再上面是金属的试样帽。试验时，土样的上下两端与透水石接触处，分别放置滤纸。试样外侧包有薄橡皮膜，膜的下端扎紧于底座，上端扎紧于试样帽。所谓压力室就是能够施加水压力或气压力的密室，侧向为有机玻璃筒，上部为金属顶盖，下部固定于底座，其间设有密封圈防止漏水，顶盖的中央为一金属活塞杆传递竖向荷载。 2.2.土体三维变形试验土体三维变形试验实验原理： 试验时在压力室中充水并加压，这一压力叫围压。围压通过橡皮

16、膜从侧向传到试样上，也通过试样帽从竖向作用给土样，此时试样受各向相等的压力：小主应力。待固结稳定后再用加压设备竖向加荷。土样上增加的竖向应力叫偏应力q（轴向附件应力），此时竖向应力为大主应力，1+q。由于土样是圆柱形的，故中主应力2 。在加竖向荷载时，可以用测微表量测试样的竖向变形量，由此可推得轴向应变a /L0，式中L0为初始试样高度。0L2.2.土体三维变形试验土体三维变形试验 三轴仪中的试样是圆柱形的，其受力和变形是轴对称的，它有两个方向的应力 和 ，同时测得两种应变 和 ，由它们可推出侧向应变 13avr2avr2.2.土体三维变形试验土体三维变形试验 三轴仪试样的应力变形状态是轴对称

17、的，而实际工程问题中土体应力应变状态往往并非轴对称的，因此需要有相应的试验设备来研究更加复杂的应力状态。 2.2.土体三维变形试验土体三维变形试验2.平面应变试验岩土工程许多问题可以简化为平面应变问题平面问题 注意：（1）三个方向的尺寸；小主应力方向尺寸应较小，试样才能达到破坏； （2）试验仪器中难以处理的问题；不同方向的加荷板打架！2.2.土体三维变形试验土体三维变形试验 2 1 3 3.真三轴试验 真三轴仪的试样为立方体，从三个方向分别施加三个主应力分量。由于加荷方式的不同，产生了不同型式的真三轴仪： （）三个方向全为刚性板加荷。2.2.土体三维变形试验土体三维变形试验 （） 方向为刚性板

18、加荷，另两方向为气压或液压柔性加荷。 （） 方向柔性加荷，而 和 方向为刚性板。1312 2 1 3 图5-12是河海大学的真三轴仪示意，属于第（）种类型。 图5-12a.整体结构b.加荷与变形示意2.2.土体三维变形试验土体三维变形试验实验原理 2方向的传力块B是由多层金属板与橡皮相间复合而成。在竖向该传力块可与试样同步压缩，而在2向靠金属板传力保持刚性。竖向荷载由试样和传力块共同承担，但荷载2.2.土体三维变形试验土体三维变形试验传感器只量测试样上的荷载，从而可算得 1。传压块B上下有滚轮，可适应试样在2向的变形。这样2方向的加荷板不要予留空隙，可使2均匀作用于试样，且试样自始至终规整。小

19、主应力则用气压施加。4.空心扭剪仪土样P1P2a.空心圆柱试样b.扭剪仪整体结构图5-132.2.土体三维变形试验土体三维变形试验外室为什么是空心而不用实心？ 仪器所用的试样为空心的圆柱体，如图13(a)所示，仪器的整体结构如图13(b)所示。试样被包在内外橡皮膜之中。该仪器可以对试样施加种荷载，径向内压力 、径向外压力 、竖向压力z和环向扭剪应力 ，根据内外径向应力可以推算出环向应力 。1 r2rz2.2.土体三维变形试验土体三维变形试验注意：研究变形与研究强度的土工仪器各有哪些？3.3.土体三向变形的主要规律土体三向变形的主要规律 利用前面所讲的一些土体应力变形试验的仪器进行试验研究可以揭

20、示土体变形的许多规律。这是建立本构模型的依据。 1. 非线性和非弹性 图5-15(a)是金属和混凝土等坚硬材料的轴向拉压曲线，图5-15(b)为土的三轴试验得出的轴向应力 与轴向应变 之间的关系曲线。与金属材料不同的是，初始的直线阶段很短，对于松砂和正常固结粘土，几乎没有直线阶段，加荷一开始就呈非线性。13a 这种非线性变化的产生，是因为除弹性变形以外还出现了不可恢复的塑性变形。土体是松散介质，受力后颗粒之间的位置调整，在荷载卸除后，不能恢复，形成较大的塑性变形。 图15a.金属b.土体3.3.土体三向变形的主要规律土体三向变形的主要规律 如果加荷到某一应力后再卸荷，曲线将如图16所示。为加荷

21、段，为卸荷段。卸荷后能恢复的应变即弹性应变。不可恢复的那部分应变为塑性应变。 图163.3.土体三向变形的主要规律土体三向变形的主要规律 经过一个加荷退荷循环后，再加荷，将如图16中的段所示，它并不与线重合，而存在一个环，叫回滞环。回滞环的存在表示退荷再加荷过程中能量消耗了，要给以能量的补充。再加荷还会产生新的不可恢复的变形，不过同一荷载多次重复后塑性变形逐渐减小。非线性和非弹性是土体变形的突出特点。弹性、塑性、粘性（流变性）3.3.土体三向变形的主要规律土体三向变形的主要规律弹性、塑性、粘性能够恢复的变形；不能恢复的变形；状态随时间而变化的性质3.3.土体三向变形的主要规律土体三向变形的主要

22、规律剪胀性和塑性体积应变 土体受力后会有明显的塑性体积变形。图17为土样在三轴仪中逐步施加各向相等的压力后，再卸除所得到的与体积应变之间的关系曲线 图173.3.土体三向变形的主要规律土体三向变形的主要规律vpevpv 可见存在不可恢复的塑性体积应变 ，而且它往往比弹性体积应变更大。这一点与金属不同，金属塑性变形是由于晶格之间的错动滑移而造成的，只有形状改变，不产生体积变化。金属是很密实的材料，晶格间没有可压缩的孔隙，因此被认为是没有塑性体积变形的。土体的塑性变形也与颗粒的错位滑移有关。这种错动滑移不仅在受剪时发生，受压时也存在。在各向相等的压力作用下，从宏观上来说，是不受剪切的；但在微观上，

23、颗粒间有错动。 3.3.土体三向变形的主要规律土体三向变形的主要规律vp 剪切也会引起塑性体积变形 在三轴仪中对土样施加偏压力 的同时，减小围压 ，并令 ，使球应力保持不变，所得出的应力应变曲线将如图19所示。尽管体积应力不变，但图中仍有体积应变，此时测得的体积应变完全是剪切造成的。在图19（）中，体积应变 随偏应力 增大而增大。剪切引起的体积收缩，叫剪缩。软土和松砂常表现为剪缩。 3132313 v313.3.土体三向变形的主要规律土体三向变形的主要规律图5-19ab3.3.土体三向变形的主要规律土体三向变形的主要规律a 在图19（）中，开始阶段为剪缩，以后曲线向上弯曲，为负的体积应变，即体

24、积膨胀，这种现象叫做剪胀。紧密砂土，超固结粘土，常表现为剪胀。3.3.土体三向变形的主要规律土体三向变形的主要规律a.松砂b.密砂图5-203.3.土体三向变形的主要规律土体三向变形的主要规律 砂土受剪所产生体积变形可用图20来说明。假定土体沿水平向受剪切。对于松砂，受剪后某些颗粒填入原来的孔隙，体积减小；对于密砂，原来的孔隙体积较小，受剪时一些颗粒必须上抬才能绕过前面的颗粒产生错动滑移，于是体积膨胀。塑性剪应变 土体受剪发生剪应变。剪应变的一部分与骨架的轻度偏斜相对应，荷载卸除后能恢复，是弹性剪应变。另一部分则与颗粒之间的相对错动滑移相联系，荷载卸除后不能恢复，为塑性剪应变。不仅剪应力能引起

25、不仅剪应力能引起剪应变，体积应力也会引起剪应变剪应变，体积应力也会引起剪应变。三轴仪中的土样在应力 和 下变形稳定后，保持 不变而降低 ,见图21（a），则会发现，随着 减小，轴向应变不断增大，直至最后达到破坏。31331333.3.土体三向变形的主要规律土体三向变形的主要规律a.单元体应力变化 c.应力路径 3.3.土体三向变形的主要规律土体三向变形的主要规律331AspBBAqp破坏线AB破坏线b.摩尔圆变化 d.剪应变 等q试验 在这一应力变化过程中，应力摩尔圆直径不变，位置不断向左移动，如图21（b），摩尔圆从移动到。当围压降到一定值，摩尔圆与库仑破裂线相切，土样剪坏，这时剪应变已发展

26、到很大数值。由此可见，球应力的变化确实引起了不可恢复的剪应变。这种应力变化可以用图21（c）中坐标系中的线段来表示。还可点绘出剪应变 s随球应力减小而增加的关系曲线，如图21（d）中的段。3.3.土体三向变形的主要规律土体三向变形的主要规律硬化和软化 三轴试验测得的轴向应力与轴向应变的关系曲线有两种形态。图22（）所示曲线有一直上升的趋势直至破坏，这种形状的应力应变关系叫硬化型,软土和松砂表现为这种形态图22(a)3.3.土体三向变形的主要规律土体三向变形的主要规律 图22（）所示曲线前面部分是上升的，应力达到某一峰值后转为下降曲线，即应力在降低，而应变却在增加，这种形态叫做软化型。紧密砂和超

27、压密粘土表现为这种形态。图22(b)3.3.土体三向变形的主要规律土体三向变形的主要规律 将转换为，点绘曲线，其形式与图22也相似，存在硬化和软化两种形式。对于其他剪切试验（如直剪、单剪），得出的关系曲线也有硬化和软化的区别。 3.3.土体三向变形的主要规律土体三向变形的主要规律应力路径对变形的影响 岩土材料存在较大的塑性变形。沿不同的应力路径加荷，各阶段的塑性变形增量不同，累积起来，就有不同的应变总量。这就是应力路径对变形的影响。图23ab3.3.土体三向变形的主要规律土体三向变形的主要规律虚线表示排水试验的有效应力路径。实线表示先做不排水试验，其有效应力路径为，达到接近破坏的点后，排水固结

28、，保持不变而增加，应力路径为。两种应力路径初始和终了应力状态相同，两种路径所对的轴向应变大不同。实线，因点接近破坏线，必然产生较大的轴向应变，最终必然较大如图23（）中的线所示。虚线远离破坏线，其轴向应变必然较小，如图23（）中的1线所示。3.3.土体三向变形的主要规律土体三向变形的主要规律应力历史对变形的影响应力历史指历史上的应力路径。由于塑性变形不可恢复，历史上发生的变形将保存和积累起来。它无疑会影响今后的变形。 图24中，、两点具有相同的应力，然而点处于初始加荷曲线上，点处于再加荷曲线上，两点对应不同的，它们所处应力应变关系曲线的斜率也不同。如果施加同样的荷载增量，则对应状态的土体应变增

29、量大，而对应状态的土体应变增量小。因、两点有着不同的应力历史，加荷后就有不同的变形。超固结土比正常固结土变形小，也是这个缘故。3.3.土体三向变形的主要规律土体三向变形的主要规律图243.3.土体三向变形的主要规律土体三向变形的主要规律又如，压缩试验中，Pc的影响图24中，、两点具有相同的应力，然而点处于初始加荷曲线上，点处于再加荷曲线上，两点对应不同的，它们所处应力应变关系曲线的斜率也不同。如果施加同样的荷载增量，则对应状态的土体应变增量大，而对应状态的土体应变增量小。因、两点有着不同的应力历史，加荷后就有不同的变形。超固结土比正常固结土变形小，也是这个缘故。各向异性 地基土一般是水平向成层

30、。由于土沉积过程中水平和竖直方向条件不同，土的结构存在着差异，使土体在许多方面表现为各向异性。应力应变关系也不例外，这种叫原生各向异性。此外，各向应力状态不同，还能引起新的各向异性。重塑土本来不存在土体结构上的两向差异，但只要各向应力不等，在应力应变关系上就会表现为各向异性，称为应力引起的各向异性。-次生各向异性3.3.土体三向变形的主要规律土体三向变形的主要规律ab图25 应力引起的各向异性3.3.土体三向变形的主要规律土体三向变形的主要规律ABfxy13A3A1A平面应变条件下，从大主应力和小主应力方向分别加载刚度矩阵或柔度矩阵为非对称：C12C21; C11=C223.3.土体三向变形的

31、主要规律土体三向变形的主要规律13A3A1AxyyxxyyxCCCCCCCCC3332312322211312118、固结应力对变形的影响 高压三轴试验资料表明，土体在高围压（固结应力）下的变形性状与低围压情况下有所不同。主要有如下三个方面： （1）强度包线不呈直线，而是呈向下微弯的曲线，如图28（a）所示。这表示有效强度指标 随着固结压力的增加而降低了。为了反映这种变化，可以用折线来代替曲线，也就是在不同的压力范围用不同的强度指标。如图28（b）所示，压力低于 A，用 1，压力高于 A，用 2。 3.3.土体三向变形的主要规律土体三向变形的主要规律ap30lg8、固结应力对变形的影响3.3.

32、土体三向变形的主要规律土体三向变形的主要规律ap30lgf0 f如何确定？ （2）有些土，如紧密砂，在低压力下受剪时体积会发生膨胀。而在高压力下，所有土都表现为剪缩，如图29中所示的v-a曲线。 （3）软化现象一般也是在低压力下表现出来的。在高压下，通常（1-3）- a曲线是硬化型的，如图29所示。3.3.土体三向变形的主要规律土体三向变形的主要规律 图29 不同围压下的应力应变曲线3.3.土体三向变形的主要规律土体三向变形的主要规律3.3.土体三向变形的主要规律土体三向变形的主要规律9、中主应力对变形的影响中主应力对强度影响平面应变条件下测定的摩擦角比轴对称条件下大35度；中主应力对应力应变

33、曲线有影响b=0; b=0.5; b=1.0 4.4.弹性非线性模型弹性非线性模型 实际工程中的初始应力状态、应力增量和应力路径是千变万化的，试验无法模拟这种复杂的变化。因此，必须通过假定、推理、验证，建立某种符合实际变形规律的数学计算方法数学模型，将少量的特定条件下试验得出的结果推广到一般，运用于工程。这种数学模型就是本构模型。4.4.弹性非线性模型弹性非线性模型线性弹性模型（广义虎克定律）3213213211111111EEEEEEEEEEzyzyzxzxyzyzyxzzzxyyzyxxGGGEEE111111 D C4.4.弹性非线性模型弹性非线性模型线性弹性模型（广义虎克定律） D 2

34、10000002100000021000000100010001)21)(1 (vvvvvvED GGGGKGKGKGKGKGKGKGKGKD000000000000000000343232000323432000323234用E、v，v0.5,不能反映剪胀用K、G，其取值基本没有限制,能反映剪胀4.4.弹性非线性模型弹性非线性模型线性弹性模型（广义虎克定律） D 210000002100000021000000100010001)21)(1 (vvvvvvED用E、v，v0.5?4.4.弹性非线性模型弹性非线性模型线性弹性模型（广义虎克定律） )1 (2000000)1 (2000000)1

35、 (20000001000100011vvvvvvvvvEC C4.4.弹性非线性模型弹性非线性模型12EG213 EB4.4.弹性非线性模型弹性非线性模型 假定土体为线弹性材料会有较大误差，因此提出非线性弹性模型。式中只含有两个参数：弹性模量E和泊松比，它随应力状态变化, 通过试验得出弹性参数随应力而变化的规律，从而建立相应公式。 )1 (2000000)1 (2000000)1 (20000001000100011vvvvvvvvvEC C4.4.弹性非线性模型弹性非线性模型弹性模量E和泊松比，体积变形模量K和剪切模量G它随应力状态变化, 通过试验得出弹性参数随应力而变化的规律，从而建立相

36、应公式。 210000002100000021000000100010001)21)(1 (vvvvvvED GGGGKGKGKGKGKGKGKGKGKD000000000000000000343232000323432000323234（三）应力应变关系矩阵2个独立参数（一）弹性参数的确定1、弹性模型E 由应力应变全量的广义虎克定律广义虎克定律，当2 =3 = 0时， ，因此， 。E1111E4.4.弹性非线性模型弹性非线性模型23111EzyzyzxzxyzyzyxzzzxyyzyxxGGGEEE111111（一）弹性参数的确定1、弹性模型E 对于粘性土，做无侧限压缩试验，此时，2 =3=

37、0 ，加竖向应力，测相应的应变 ，见图32（a）。点绘出轴向应力与轴向应变 的关系曲线，如图32（b）所示；曲线上的点A所对应的大主应力1 =，大主应变 = ，故弹性模量E11aa1aasE它是曲线在点A 处的割线的斜率，故称割线模量，以Es表示。4.4.弹性非线性模型弹性非线性模型23111E4.4.弹性非线性模型弹性非线性模型 b：割线模量a图32 c：泊松比 对于增量广义虎克定律，用Et和 t分别表示模量和泊松比，则 无侧限压缩试验时， , 自然有 。因此， 曲线的切线斜率，就等于增量关系中的弹性模量，称切线模量， tttEddEdd3211032dd11ddEtaatddE4.4.弹性

38、非线性模型弹性非线性模型zyzyyxzzzxzxzxyyyzyzzyxxGEGEGE1111112、泊松比写出广义虎克定律的另一个式子EE2133做无侧限压缩试验， E13，而11E故 1313 量测侧向膨胀应变 ,它就是 。点绘 关系曲线，见图32（c），则线上一点割线的斜率就等于全量的泊松比，叫割线泊松比，以 s表示， . 。自然，曲线的切线斜率具有切线泊松比的意义，以 t 表示 r3ar)(ars4.4.弹性非线性模型弹性非线性模型artdd它用于反映增量的应力应变关系。图32 (c)4.4.弹性非线性模型弹性非线性模型复杂应力状态下E、v的定义？3、体积模量B 在三轴仪中对土样施加各向

39、相等的压力3，也就是球应力p，逐步增大，测相应的体积应变 可点绘出 关系曲线，如图5-34 所示，其割线斜率为割线体积模量vpv图5-34vspB切线的斜率为切线体积模量vtddpB4.4.弹性非线性模型弹性非线性模型tBsBpv4、剪切模量G 在三轴仪试验中对土样逐级施加偏应力 测相应的轴向应变 ，可求得偏应变为31qv3vas八面体剪应力 ：qoct32八面体剪应变：soct2 利用八面体面上的剪应力剪应变间的关系可确定剪切模量4.4.弹性非线性模型弹性非线性模型soctoctqG3soctoctqG3 点绘 关系曲线，可得割线剪切模量Gs 和切线剪切模量G t 如图35所示。 sq3图3

40、54.4.弹性非线性模型弹性非线性模型3qsGtGs( (二二) )、双曲线模型、双曲线模型( (邓肯邓肯- -张模型张模型) )一.切线弹性模量 康纳等人发现，在加荷时， 关系曲线可以用双曲线来拟合，如图5-36（）所示。对于某一小主应力 来说，可以用下式表示31a3aaba31式中，和为试验常数。上式也可写成aaba314.4.弹性非线性模型弹性非线性模型 若以 为纵坐标， 为横坐标，构成新的坐标系，则双曲线转换成直线，见图36（），其斜率为，截距为。31aaab图5-364.4.弹性非线性模型弹性非线性模型 邓肯和张利用上述关系推导出了弹性模量公式。在 不变条件下，由式 得增量的弹性模量

41、aatE313111将式 代入上式得2atbaaE由式 得baa3114.4.弹性非线性模型弹性非线性模型aaba313atddE1代入式 得23111baEt4.4.弹性非线性模型弹性非线性模型 现在来研究和的意义以及它们随应力的变化。由式 可见，当 时aaa4. 4.弹性非线性模型弹性非线性模型a而 是曲线 的初始切线斜率，其意义为初始切线模量，用Ei来表示，见图5-36（）。因此这表示是初始切线模量的倒数。试验表明， 随3 变化。如果在双对数纸上点绘 和 的关系，则近似地为一直线，如图5-37所示。iEa14. 4.弹性非线性模型弹性非线性模型图5-37这里 为大气压力。引入 是为了使纵

42、横坐标化为无因次量。直线的截距为 ，斜率为 。于是有 ， ,由此得 ，4. 4.弹性非线性模型弹性非线性模型aa由式 还可见，当 时这里用 表示当 时 的值，也就是 的渐近值。实际上， 不可能趋向无穷大，在达到一定值后试样就破坏了，这时的偏应力为 ，它总是小于 a 叫破坏比。将式 和式 代入式 ，并利用上式，得4. 4.弹性非线性模型弹性非线性模型ffRb3123111baEt令 ffuRb313114. 4.弹性非线性模型弹性非线性模型fs3131叫应力水平。它表示当前应力圆直径与相同小主应力 条件下破坏应力圆直径之比，反映了强度发挥程度。式（5-43）也可写成令4. 4.弹性非线性模型弹性

43、非线性模型 破坏偏应力与固结压力有关，由下图中的几何关系不难推出 将上式和式（5-40）代入式（5-43），得 Et随应力水平增加而降低，随固结压力增加而增加。包含个参数，、为强度指标，另外三个数K、和Rf4. 4.弹性非线性模型弹性非线性模型个参数确定：、为强度指标，另外三个数K、和 Rf的确定方法在推导中已作说明Rf对不同的3会有不同的值，取平均值 naaipKpE3二、切线泊松比 图5-35中所示为三轴试验测得的 关系曲线。利用式（5-21）可由体积应变推出侧向应变 。普通三轴仪竖向加荷，侧向为膨胀应变，故 为负值。点绘 关系曲线，如图5-41（）所示。库哈威和邓肯也用双曲线来拟合，与式

44、（5-34）相似，将有avrr）（rarraDf式中和为两个参数。将上式转换成rarDf4.4.弹性非线性模型弹性非线性模型ab图41 它表示，若以 为纵坐标， 为横坐标，则试验关系将为一直线，如图41（）所示 arr4.4.弹性非线性模型弹性非线性模型由上式可得aarDf1 如果将图5-41（a）纵横坐标交换，就成为 曲线。其斜率为增量泊松比ar ）（ 将式（5-49）代入上式，并利用式（5-35a）把所含 用应力来代替可得a4.4.弹性非线性模型弹性非线性模型arart（5-49）baa31121Aftsin2cos2sin11331331cRpKpDAfnaa其中4.4.弹性非线性模型弹

45、性非线性模型4. 4.弹性非线性模型弹性非线性模型 由式（5-47），当（- ）时， 。可见是 渐近值的倒数。当（- ）时 式中，为初始切线泊松比，即各向等压状态下的泊松比。对于不同的，有不同的 值，在半对数纸上点绘 与 关系曲线，近似为一直线，如图5-42所示。其截距为，斜率为。于是有初始切线泊松比为4. 4.弹性非线性模型弹性非线性模型图5-42最后得切线泊松比公式为 4. 4.弹性非线性模型弹性非线性模型 上式算得的 有时可能大于.，在试验中测得的值也确有可能超过.，这是由于土体存在剪胀性。然而，有限元计算中，若大于或等于.，劲度矩阵就出现异常。因此，实际计算中，当.时，令. 参数G、F

46、、D。三、切线体积模量 邓肯等人还提出了一种确定切线体积模量 的方法，并在有限元计算中使用 和 两参数。这习惯上被称做模型，而计算中使用 和 者被称作模型。对于E-B模型，二维问题的刚度矩阵可表示为tBtEtBtEt EEBEBEBEBEBBD0003303393 在三轴试验中施加偏应力 ，则平均正应力的变化为 。因此313131p4.4.弹性非线性模型弹性非线性模型vtB3131 邓肯等人假定，与应力水平无关，或者说与偏应力 无关，它仅仅随固结压力而变。对于同一， 为常量。由上式可见，这相当于假定 与 成比例关系。根据这种假定，对同一，如果点绘 /3 关系曲线，应为一直线，如图5-43所示。

47、31tBv3131v4.4.弹性非线性模型弹性非线性模型3131p7 . 07 . 0313svstB图5-434.4.弹性非线性模型弹性非线性模型邓肯等人取与应力水平.相应的点与原点连线的斜率作为平均斜率，即4. 4.弹性非线性模型弹性非线性模型 对于不同的， 也不同。在双对数纸上点绘 与 关系曲线，可得一直线，如图5-44所示，其截距为 ，斜率为，则四、回弹模量 对卸荷和再加荷的情况，试验表明应力应变关系曲线与加荷是不一样的，应该由回弹试验测定弹性模量。 图394.4.弹性非线性模型弹性非线性模型31aEurEt加荷卸荷再加荷 在图5-39中，为加荷状态的应力应变关系曲线，其斜率为t；而卸

48、荷与再加荷的曲线之间略有差异，中间本应存在一个回滞环，这里近似假定它们一致，且为一直线，如AB所示，其斜率为 。它具有卸荷再加荷情况下弹性模量的物理意义，叫回弹模量。urE 在曲线的不同位置卸荷，回弹模量略有不同。邓肯等人认为可以忽略这种差异，假定 不随 变化。但对于不同的围压 ，可测出显著不同的 , 即 须随而 变。 urE313urEurE34.4.弹性非线性模型弹性非线性模型 在双对数纸上点绘 与 关系曲线，可得一直线，如图5-40所示，其截距为 ，斜率为。这样回弹模量可由下式计算aurpElgap3lgurKnaaururppKE3图5-404.4.弹性非线性模型弹性非线性模型回弹标准

49、：4.4.弹性非线性模型弹性非线性模型 03131且时，用Eur；否则，用Et 加荷函数：433131)(flffl0.75(fl)max：Eur0.75(fl)maxfl(fl)max： 75. 011)(maxllturtffEEEE（四）参数的确定及参数的变化范围v模型：F、C、Rf、K、n、G、F、D、模型：用， Rf、 K，K和共个参数 作三轴排水试验测轴向应变 a和体积应变 v。点绘应力应变关系曲线。按以下方法确定参数： （）在坐标系中绘不同3时的极限摩尔圆。如图5-38所示。 4.4.弹性非线性模型弹性非线性模型如图5-38 （）对每一个3，点绘 - 关系曲线，近似取为直线，如图

50、5-36（）所示。 31aa4.4.弹性非线性模型弹性非线性模型图5-36（）ffuRb31311（）对各用式（5-42）计算。 4.4.弹性非线性模型弹性非线性模型 （3）点绘 与 关系曲线，如图5-35所示。 aipElgap3lgnaaipKpE3若值小，则曲线在高应力水平下仍然挺拔，坡度较大。 （）由试验关系，利用式（5-21）转换为关系。 4.4.弹性非线性模型弹性非线性模型GFiap3log （）如果要确定体积模量有关参数，可从 - - 关系曲线上，找出应力水平.时的 值近似地令v31avvB331 然后点绘 与 关系曲线，如图5-42所示从拟合的直线可定K和。atpBlgap3l

51、g4.4.弹性非线性模型弹性非线性模型图5-42参数数值范围说明（有效内摩擦角） 软土较低，砂土中等，堆石较高（有效粘聚力） .K502500软土较低，砂土中等，堆石较高n01.0Rf0.50.95G0.20.6F0.010.2D120.0Kb（0.33.0）Km01.0Kur（1.23.0）K（五）对双曲线模型的讨论1. 模型和模型的差异两种模型采用了相同的弹性模量，无论加荷还是卸荷。它们的差异仅在于泊松比 和体积模量。 在有限元计算中所得结果差异比较大；E-B模型实际上用了两种泊松比 t 和 ur 。这一点与v模型无论加荷还是卸荷都取一种泊松比是不同的。EB模型的泊松比较小；回弹状态时ur

52、 m4.4.弹性非线性模型弹性非线性模型 tmnabftpKKSR3231121BEt31212.双曲线模型反映的变形特征（3）模型将回弹模量与加荷模量区别开来。 （4）由于模型用于增量法有限元计算，故能在一定程度上反映应力路径对变形的影响。图5-46中有两条应力路径:起点A，终点B；ADB变形要大于ACB。 4.4.弹性非线性模型弹性非线性模型（）通过，或随应力的变化反映了非线性变形特性 ；实际计算结构时，各单元的E用Ev和E-B模型计算是不同的（）模型较恰当地反映了剪切变形随应力水平的增加而增加，随围压的增加而减小；体积变形随围压的增加而减小。 4.4.弹性非线性模型弹性非线性模型图5-4

53、6（4）应力路径对变形的影响：图5-46中有两条应力路径:起点A，终点B；ADB变形要大于ACB。 3.双曲线模型的存在问题（）不能反映剪胀剪缩性（）不能反映软化特性（）不能反映各向异性，模型的建立和参数的确定完全依据大主应力方向加荷的试验结果。 4.4.弹性非线性模型弹性非线性模型aarDf1rav2aaavDf12)1 (1aavDf（4） 模型对卸荷状态取回弹模量，以便区别加荷模量。但没有提出卸荷时的泊松比和卸荷时的体积模量。 4.双曲线模型实用性4.4.弹性非线性模型弹性非线性模型首先因为它能反映土体变形的主要特点，尽管还存在许多问题。 第三，在广泛应用中逐渐积累了使用它的经验。 其次

54、，该模型有一个重要特点就是简单。 1. KG模型4.4.弹性非线性模型弹性非线性模型( (二二) )、其它非线性模型、其它非线性模型qpGGpKKGGitKit式中， 是非线性KG模型的五个参数。一般地，0，0，0。可由各向等压固结试验的曲线确定，由等p的三轴固结排水剪切试验确定。5.5.弹塑性模型弹塑性模型 弹塑性模型把在荷载作用下所发生的变形分成两部分：一是弹性变形，即可恢复的变形；另一是塑性变形，即不可恢复的变形。 弹性变形可以用广义虎克定律来求解，塑性变形则要用塑性力学的方法来建立应力应变关系。 peddd破坏准则和屈服准则 硬化规律 流动法则 三大假定：5.5.弹塑性模型弹塑性模型

55、为了建立塑性变形的关系式，须要给出三个条件，即三方面的假定。 对这三个假定采用的具体形式不同就形成了不同的弹塑性模型。 一、屈服准则和破坏准则一、屈服准则和破坏准则 塑性变形并不是任何荷载下都发生，而是加到一定的荷载时才发生。究竟达到什么应力状态才发生，就要有一个判别标准，这个标准叫做屈服准则 。 荷载逐步增加，变形逐步发展，达到某种应力状态时，土体破坏。这种应力状态是一种极限应力状态。判断是否达到这一状态的标准，叫破坏准则。 屈服是发生塑性变形的下限应力状态，破坏则是其上限应力状态。因此讨论塑性变形不能不明确其上下限。 5.5.弹塑性模型弹塑性模型一、屈服准则和破坏准则一、屈服准则和破坏准则

56、屈服-准则 破坏-准则5.5.弹塑性模型弹塑性模型 ij ij （一）破坏准则 土体的破坏决定于应力状态，故破坏准则可写成fijkf 其中， 是应力分量的某种函数，叫破坏函数；k f是试验确定的常数。 ijf 破坏函数在主应力空间内代表一曲面，叫破坏面。若表示应力状态的点落在破坏面以内，材料不破坏；若落在破坏面上，材料破坏。应力状态永远不能超出破坏面。 5.5.弹塑性模型弹塑性模型fijkf 土体的破坏有两种型式：剪切破坏和拉裂破坏。但由于土的抗拉强度小，往往假定为0，不考虑拉裂破坏。下面所讲的破坏准则主要指剪切破坏准则 。建立破坏准则，其实就是将抗剪强度以某种应力分量的形式表达出来。 5.5

57、.弹塑性模型弹塑性模型 破坏准则主要有以下几种：1、屈雷斯卡 准则 这一准则是假定最大剪应力达到某一数值时破坏，即fk231如果主应力的大小不确定，可写成0222222322113231231ffffffkkkkkk5.5.弹塑性模型弹塑性模型 它在主应力空间内是以空间主对角线为中心轴的正六角柱面，见图5-48（a）。上式也可用应力不变量来表示，即06496362746242222332fffkJkJkJJ5.5.弹塑性模型弹塑性模型321图5-48常数fk平面形状？屈雷斯卡 准则:破坏与球应力无关 132Mohr-Coulunb Mises Trasca 屈雷斯卡 准则包含了破坏与球应力无关

58、的假定；对岩土体，破坏与球应力有关，故5.5.弹塑性模型弹塑性模型图5-481312Ikf广义屈雷斯卡 准则平面形状？饱和土不排水剪强度为uc231 其摩尔图包线如图5-48（b）所示，Cu为不排水粘聚力，与式（5-54）相比可见，kfCu，屈雷斯卡准则很好地体现了饱和土不排水条件下的强度特征。5.5.弹塑性模型弹塑性模型c2、米塞斯 准则 该准则是假定偏应力q达到一定值时破坏，即fkq21323222121 它在主应力空间为圆柱面如图5-49（a）所示。5.5.弹塑性模型弹塑性模型312平面形状？ 132Mohr-Coulunb Mises Trasca 子午面上破坏曲线的形状？b5.5.弹

59、塑性模型弹塑性模型312321广义米塞斯准则考虑球应力对强度的影响，将kf用I1 （I1=1+2+3）的函数代替，则成为广义米塞斯准则，在主应力空间为圆锥面，如图5-49（b）所示。)(12IkJfb图5-495.5.弹塑性模型弹塑性模型321破坏线Mpprq广义米塞斯准则在qp二维应力空间为一倾斜直线如图5-49（c）所示。剑桥模型中破坏准则为q=Mpq=M（p+pr）子午面上破坏曲线的形状？5.5.弹塑性模型弹塑性模型321DruckerPrager准则fkJI213、摩尔库仑 准则 对土体，摩尔库仑强度理论已将某一面上的抗剪强度转换为达到破坏时单元体主应力之间的关系，它是现成的破坏准则。

60、其表达式为cossin223131c2452245231ctgtg或5.5.弹塑性模型弹塑性模型在保持 相同的情况下，按摩尔库仓准则，三轴压缩和拉伸具有相同强度（不受 影响）。而实际试验结果，在 相同情况下，拉伸试验所得的强度常高于压缩试验测得的强度。 它表示破坏与破坏与2 2无关无关，在主应力空间破坏面是2轴平行的面，且投影到1轴与3轴构成的平面内，是一直线。如果1、2、3的大小不确定，则为6个面，它们在主应力空间构成不等角的六角锥面，如图5-50所示图5-505.5.弹塑性模型弹塑性模型 PI平面内是一等边不等角的六边形。图5-505.5.弹塑性模型弹塑性模型 132Mohr-Coulun