土的工程特性及其本构理论2010

土的强度、本构理论及固结方程大连理工大学岩土所唐小微第1章概论1.1

内容土的工程特性包括土的应力应变关系，强度理论。土的本构理论（constitutiverelationship

）是反映材料的力学性状的数学表达式，表示形式一般为应力-应变-强度-时间的关系。本构关系也称为:本构定律（constitutivelaw）本构方程（constitutiveequation)数学模型（mathematicalmodel

）

1.2

土的构成土是岩石经风化，搬运，沉积，固节和地壳运动而成的碎散矿无颗粒的集合体。土包含三相:固相液相气相

1.3

土力学的产生与发展诞生，《土力学和地基基础》K.Terzaghi，1925。前期基础:莫尔-库仑强度理论有效应力原理饱和黏土的一维固结理论

主要问题变形问题弹性理论稳定问题极限平衡分析刚塑性和理想塑性理论土的本构关系研究的历史与发展20世纪60年代，高重建筑物及深厚基础问题,及计算机技术发展为土的本构关系研究建立了必要性和可能性；80年代达到高潮；目前的发展方向：土的结构性、非饱和土、循环加载、动力本构模型等。数值分析连续介质有限元，有限差分，边界元，无单元不连续离散元DEM，不连续变形分析DDA，流行元MEM，颗粒流PFC第2章土的工程特性2.1土工基本实验--三轴仪与试验(1)1930

A.Casagrande

提出圆柱试样(2)1933Seffert

用三轴仪研究固结(3)1934Rendulic

利用测定土的强度参数(4)1959

黄文熙，汪闻韶研制动三轴试验(5)三轴仪的发展：动三轴、大尺寸三轴仪、高压三轴、非饱和土三轴仪、应力与应变路径控制三轴仪……(6)变形与强度试验(7)应力状态明确，简单；排水条件明确图2－1三轴仪简图图2－2常规三轴压缩试验的应力状态简图两个参数p,q

对于三轴应力状态：

=

这两个参数对于三轴应力状态使用比较简便图2－3三轴试验中的剪切破坏情况2.2

土的应力变形特性特性：非线性、弹塑性、压硬性、剪胀性、各向异性、结构性、流变性、应变硬（软）化、减载体缩。影响因素：应力水平、应力路径、应力历史(3S:stresslevel,stresspath,stresshistory)、其它。2.3.1研究历史4.广义密塞斯（Mises）和广义屈雷斯卡（Tresca）5.现代的强度理论：破坏是应力应变关系的最后状态：包括在本构关系模型之内6.与时间有关、拉伸、断裂及孔隙水压力：水力劈裂1.1776年，库仑（Coulomb)公式：2.1900年，莫尔(Mohr):3.土的抗剪强度

f是作用在其破坏面上的正应力

n的单值函数2.3土的强度2.3.2土的强度的特点1.土是碎散颗粒的集合，颗粒之间的相互联系是一般相对薄弱的。所以土的强度主要是由颗粒间的相互作用力决定，而不是由颗粒矿物的强度本身决定的。2.土的破坏主要是剪切破坏，其强度主要表现为抗剪（摩擦）强度。3.粘聚力：颗粒间的连接－粘聚力。4.三相组成，固体颗粒之间的液体、气体及液、固、气间的界面对于土的强度有很大影响：孔隙水压力、吸力（毛细力）。5.地质历史造成土强度强烈的多变性、结构性和各向异性。6.土强度的这些特点体现在它受内部和外部、微观和宏观众多因素的影响，成为一个十分复杂的课题。1.屈服与强度：刚塑性弹－完全塑性应变软化断裂弹塑性图2－23土的几种本构关系模型2.土的强度和土体破坏1)土达到屈服不一定达到破坏2)在土体中，局部土达到强度，不一定引起土体的破坏3)渐进破坏与崩塌、断裂塑性区部分土体达到强度（屈服），地基并不一定破坏。图2－24

土中的塑性区3.土破坏的判断1）破坏是应力体变过程的最后阶段，这时微小的应力增量将会引起很大的，或者不可控制的应变增量；2）土的破坏主要是剪切破坏；3）有时用应力比和应力差判断破坏是不一致的。

1－

3峰值强度残余强度图2－25土的几种破坏形式定义断裂一定应变值破坏是应力体变过程的最后阶段，这时微小的应力增量将会引起很大的，或者不可控制的应变增量。4

影响土强度的因素

e为土的孔隙比；

C代表土的组成，component；

H代表应力历史，history；

T表示温度，temperature；

和

分别表示应变和应变率；

S表示土的结构，Structure；

c和

为粘聚力及内摩擦角。其中各种因素并不独立，可能相互重叠。（1）材料的强度是指材料破坏时的应力状态。（2）定义破坏的方法（数学表达式）是破坏准则。破坏准则常常是应力状态的组合。（3）强度理论是揭示土破坏的机理的理论，它也以一定的应力状态的组合来表示。因而强度理论与破坏准则的表达式是一致的。2.3.3

土的强度理论一般表达式对于各向同性材料或者应力张量的主应力和应力不变量＝用主应力表示主应力：在三个没有剪应力的方向上的正应力4.

球应力张量与偏应力张量＝＋偏应力张量sij偏应力张量的不变量土力学中常用的三个应力（不）变量平均主应力p广义剪应力q应力洛德角

2.2.2应变1.与应力的情况相似2.体应变3.广义剪应变4.应变洛德角＝土力学中常用的应力应变关系表示线弹性：非线弹性：塑性：

土的经典强度理论1.特雷斯卡（Tresca）准则及其广义准则2.米泽斯（VonMises）准则及其广义准则3.莫尔－库仑（Mohr-Coulomb)强度准则1.特雷斯卡（Tresca）准则与广义特雷斯卡（extendedTresca）准则对岩土，考虑平均主应力的影响，广义形式六棱柱的表面：图2－26特雷斯卡与米泽斯准则图2－27广义的形式锥面2.米泽斯（VonMises）和广义米泽斯（extendedVonMises）准则广义米泽斯——Drucker-Prager准则考虑平均主应力对土抗剪强度的影响图2－28米泽斯和广义米泽斯准则

1

3

2圆柱面与圆锥面图2－29

平面上的各强度准则3.莫尔－库仑强度准则莫尔（Mohr），1900单值函数在一定的应力范围，线性关系－库仑公式三轴平面图2－30莫尔－库仑强度准则三维空间

平面近代的强度理论1.

莱特－邓肯(Lade-Duncan)强度准则2.

松冈元-中井照夫（Matsuoka-Nakai）破坏准则3.

双剪应力强度理论4.

隐式的破坏准则本构关系－应力应变与强度关系土的强度，或者破坏是其应力应变过程的最后阶段，即在微小的应力增量下，会产生很大（或者不可控制）的应变增量。因而破坏是应力应变关系的最后阶段。第3章土的本构理论3.1

土的弹性模型3.1.1概述3.1.2线弹性：（广义）胡克定律3.1.3非线弹性：增量胡克定律3.1.1概述1.线弹性模型一般不适用于土，但有时还是可以近似使用的：地基中应力计算；分层总和法（分段线性）。2.非线弹性模型使用最多，有很好的实用性：一般参数不多；物理意义明确；使用的试验比较简单；使用增量广义胡克定律的形式

(Duncan-ChangModel)。3.1.2线弹性模型（广义）胡克定律广义胡克定律（各向同性）土力学中常用的K、G形式张量表示矩阵表示3.1.3非线弹性－增量的广义虎克定律增量广义虎克定律（各向同性）E=Et

t切线模量1.Duncan-Chang双曲线模型

1）基本原理Kondner在1963年所做的三轴试验中，应力应变可用双曲线模拟。

0图3－1三轴试验的应力应变曲线在常规三轴压缩试验中：＝＝所以1/a

代表了双曲线的初始斜率（模量）

1

1

ult=1/b1/b

代表了曲线应力的极终值图3－2参数a和b的物理意义

(

)ult=1/b1Ei=1/a0

a

1b02）参数a,b确定利用常规三轴压缩试验数据图解确定参数a,b图3－3二者的线性关系与试验结果EtEur

1

1-

3图3－4

初始模量Et与卸载（再加载）模量Eu（Et、Eur随着围压

3的增加而加大)(Eur>Et

）lg(

pa)lg(Ei/pa)(

)1(

)3(

)21K1nKur初始模量与卸载（再加载）模量（三个围压的试验在双对数坐标下的结果）图3－5模量参数的确定：破坏比Rf=

1

1-

3(

1-

3)ult(

1-

3)f图3－6破坏比Rf15%0=模型的切线模量Et3）E-B模型采用切线变形模量和体积模量表示Et:切线变形模量Bt:体积模量胡克定律的一般公式：E-B模型:

对于同一围压的试验假设体积模量B是常数

70%(

-

)f

v70%图3－7体积模量B的确定对于同一个围压

3，B为常数，对于不同的围压，它与

3成指数关系试验参数Kb,m3.2土的弹塑性模型的一般原理3.2.1塑性理论在土力学中的应用3.2.2屈服准则与屈服面3.2.3流动规则与硬化定律3.2.4弹塑性本构模型的模量矩阵的一般表达式3.2.1塑性理论在土力学中的应用1.1776年库仑公式与土压力理论—刚塑性；2.借鉴金属塑性理论，弹性－理想（完全）塑性；3.1960s，弹塑性理论模型：在增量意义上是弹－塑性的。弹塑性理论回顾：Drucker假说屈服准则流动—正交法则硬化规律弹性－完全塑性elastic-perfectlyplastic增量弹塑性-incrementalelastic-plastic图3－9几种塑性模型刚塑性rigid-plastic不同塑性模型的应用1)刚塑性理论-极限平衡法：刚体滑动法、各种条分法、滑移线法（不计变形，不计过程）。2)弹－塑性理论：在一定范围为弹性，超过某一屈服条件为塑性变形。数值计算中出现“塑性区”。3)（增量）弹塑性理论模型：一开始就是弹塑性变形同时发生，屈服面不断发展。3.2.2屈服准则与屈服面1.屈服准则2.屈服函数3.屈服面与屈服轨迹4.土的屈服面与屈服轨迹的一般形式5.土的屈服面与屈服轨迹的确定1.屈服准则(yieldcriterion)判断是否发生塑性变形的准则－判断加载与卸载的准则A

BABA

B为屈服点；A´非屈服点图3－10屈服－弹塑性应变的判断准则A、B在屈服面上，A

B

不在屈服面上B

B

A

2.屈服函数

(yieldfunction,yieldequation)屈服准则的数学表达式对于刚塑性和弹性－塑性模型：H为常数；对于弹塑性模型：H是塑性应变的函数。加卸载的判断（应变硬化情况）为中性变载，只产生弹性变形为卸载，只产生弹性变形为加载，同时发生弹性、塑性变形f<0在屈服面之内，弹性变形3.屈服面与屈服轨迹屈服面－屈服准则在应力空间中的几何表示：1)三维应力空间：屈服面2)二维应力空间：屈服轨迹4.土的屈服面与屈服轨迹的一般形式1）由于土是一种摩擦材料，人们认为只是在应力比变化时颗粒间才会相对滑动位移（Mohr-Coulomb,；广义Mises；广义Tresca:锥形屈服面)qp图3－11锥形屈服面与射线屈服轨迹pq2)又由于土在各向等压条件下也会发生颗粒相对运动，土变密实，所以出现各种“帽子”屈服面（Cam-clay,；清华模型）qpqppq图3－12帽子屈服面3）二者的联合形式qP-图3－13普遍形式的屈服面5.土的屈服面与屈服轨迹的确定1)假设屈服面与屈服函数2)通过试验试加载勾画屈服轨迹3)通过试验确定塑性应变增量的方向和Drucker假说确定塑性势面＝屈服轨迹3.2.3流动规则与硬化定律1.流动规则(flowrule)2.硬化定律(strain-hardeninglaw)1.流动规则（flowrule)：用以确定塑性应变增量向量的方向的规则（或者确定塑性应变增量的各个分量间的比例关系）－塑性应变增量向量正交于塑性势面。所以也称为正交规则。相适应（相关联）的流动规则（associatedflowrule)：根据Drucker假说，塑性势面必须与屈服面重合，即f=g。不相适应（不相关联）的流动（nonassociatedflowrule)：塑性势面不必与屈服面重合f

g。2.加工（应变）硬化定律(strain-hardeninglaw)：是确定在一定的应力增量作用下引起的塑性应变增量大小的规律。硬化参数H(

pij):

是土在发生了一定的塑性应变后，其排列与组构变化的尺度。A:塑性硬化模量3.2.4弹塑性本构模型的模量矩阵的一般表达式{d}不相适应f

g相适应f=g4.水土耦合固结方程4.1单向固结的普遍方程4.2太沙基（Terzaghi）单向固结理论4.3三向压缩比奥（Biot）固结理论4.4太沙基（Terzaghi）-伦杜立克（Rendulic）准三维固结理论（扩散方程）4.5两种三向固结论理的比较4.1单向固结的普遍方程1.竖向力平衡条件：以dz,单位面积孔隙水为隔离体1）渗透力的反作用力Fz2）孔隙水自重＋骨架浮力的反作用力：

w

dz3）水压力增量：图5－32土单元的竖向力1.平衡方程v＝出逸流速(虚拟）（dischargevelocity)（1）2.连续性方程连续性条件（2）3.土骨架的应力应变关系（3）（2）H1z（4

）(4

)H图5－33土层剖面式（3）代入式（1）得：得：(4

)H：土层厚度（4）4.单向固结的普遍方程(4)图5－34土层内的应力分布外荷载随时间变化土层厚度随时间变化渗透系数是深度的函数单向固结的普遍方程的意义4.2太沙基（Terzaghi)的单向固结理论太沙基（Terzaghi)的单向固结理论适合下列情况：（1）

土体是均质的，完全饱和的；（2）土粒与水均为不可压缩介质，土体变形完全是由孔 隙水排出和超静水压力消散所引起；（3）土的渗透性k与压缩性mv均为常量（应力与应变直 线关系）；（4）

外荷重瞬时加到土体上，在固结过程中保持不变；（5）土体中只引起单向的渗流与压缩；（6）

土中渗流服从达西定律。＝0（5）M=π(2m+1)/2式（5）的孔隙水压力的消散一点的固结度图5－35土层的固结初始孔压分布(1)均布

(2)直线

(3)半正弦

(4)正弦

图5－36初始孔压分布一层土的固结度的表示二者是完全等价的图5－58圆形基础下土层的三维固结曲线Terzaghi一维固结曲线一维（单向）与三维固结计算的区别4.3三向压缩比奥（Biot）固结理论zyx1.平衡方程

fi为体积力,以土体为隔离体图5－46单元体上的应力

以土骨架作隔离体的平衡方程(1)三个方向上的渗透力：ix

w,iy

w,iz

wu:为超静水压力时，

为浮容重

；u:为总水压力（包括静水压力）时，

为饱和容重

sat。2.位移协调条件：应变-位移条件：土骨架在x,y,z方向的位移(2)3.土骨架的应力应变关系-线弹性广义胡克定律(3)或者(3ˊ)或者(3

)平衡、变形协调及本构关系三方程叠加(4)方程及未知数个数未知数4个：us,vs

ws

：土骨架的位移u：孔隙水压力三个方程少一个条件(4)4.饱和土体的连续性方程流出水量＝体积压缩达西定律连续性方程(5)（1）

Cv3是三维固结系数；（2）

是时间t

的函数。比较：(5)单向固结微分方程(5)(4)5.二维与一维形式平面应变一维形式：单向渗流固结问题对于荷载一次施加，并且不变可见，此时比奥理论与太沙基单向固结理论一致6.比奥固结理论原理及其在数值计算中应用（1）未知变量：结点的us,vs,ws;u；（2）有效应力原理；（3）平衡方程；（4）连续性方程；（5）变形协调条件；（6）本构模型：线性，非线性，弹塑性；（7）时间：从t=0开始，每次增加

t；（8）应力应变的非线性：不同时刻参数随有效应力变化。4.4太沙基（Terzaghi)-伦杜立克（Rendulic）准三维固结理论(扩散方程)

根据一维固结论理，将固结方程进行重要的简化，解决二、三维固结问题。骨架体应变：假设：骨架体变率：1.变形条件骨架体应变率：连续性方程：微分方程：2.二维与一维的形式二维一维3.固结系数的比较太沙基一维固结理论二者的固结系数是一致的4.5两种固结论理的比较——原理与条件比奥固结理论太沙基一维固结理论1.两种固结微分方程比奥(Biot)太沙基（Terzaghi)-伦杜立克(Rendulic)2.理论假设的比较相同之处线弹性（?）小变形（小应变）达西定律连续性条件：饱和、不可压缩主要区别是否假设正应力之和在固结与变形过程中为常数；实际上为是否满足变形协调条件。平衡方程（＋有效应力原理）应力应变关系——线弹性模型（也可以是其他模型）应变－位移关系：变形协调条件连续性条件太沙基-伦杜立克胡克定律平衡方程连续性条件不满足2.

理论建立条件的比较比奥理论比奥理论可解得土体受力后的应力、应变和孔压的生成和消散过程，理论上是完整严密的。太沙基-伦杜立克扩散方程假设三个主应力（总应力）之和不变，不满足变形协调条件，（应力应变解不严密）。只能解出孔隙水压力u。两种固结理论的比较3．曼代尔-克雷尔效应(Mandel-CryerEffect)

在不变的荷重施加于土体上后的某时段内，土体内的孔隙水压力不是下降，而是继续上升，而且可能超过应有的压力值。

该现象由曼代尔（Mandel）和克雷尔(Cryer)发现，故称为曼代尔-克雷尔效应，或称应力传递效应。时间因数lgTv