第三土的弹塑性模型演示文稿

第三土的弹塑性模型演示文稿目前一页\总数五十九页\编于十七点（优选）第三土的弹塑性模型目前二页\总数五十九页\编于十七点

模型~§3.4弹塑性模型简介目前三页\总数五十九页\编于十七点§3.4弹塑性模型简介

一、模型

1、切线弹性模量该模型是通过非线性数学表达式，建立应力状态改变过程中弹性常数、与偏应力（σ1-σ3）之间的双曲线函数，从而更客观地反映非线性材料的应力应变关系。其中邓肯-张（Duncan-Chang）模型被国内外较为广泛的应用。首先，考特纳(Kondner)建议采用双曲线表示（σd–ε1）曲线，对于常规三轴压缩试验，当σ3为常数时，有：其目的就是通过常规三轴压缩试验曲线，找出一个共同的数学表达式，并从这个表达式出发，推导出切线模量Et的计算公式，以供增量弹性分析之用,我们简称E——μ模型。目前四页\总数五十九页\编于十七点§3.4弹塑性模型简介

一、模型

或

式中：a，b为试验常数。由常规三轴试验，可以得到图示曲线11三轴压缩试验曲线渐近线目前五页\总数五十九页\编于十七点§3.4弹塑性模型简介

一、模型

图中为应力应变曲线的渐近线，即为其破坏的极限值，该模量为初始弹性模量，在（-）~曲线上任取一间隔段求导，即为该增量区段内的切线模量：-当极小趋近于零时，起始弹性模量Ei目前六页\总数五十九页\编于十七点§3.4弹塑性模型简介

一、模型

1b即弹性模量的倒数。如将图中的坐标系纵轴变为斜率：式中为破坏比。Rf的取值在0.75~1.0之间。该直线纵轴截距为a目前七页\总数五十九页\编于十七点§3.4弹塑性模型简介

一、模型

在三轴试验中，随取值不同，（-）~曲线将不同。或式中：为的函数，随σ3而变动,可采用下式计算：邓肯-张提出这些试验曲线可用下式表达:目前八页\总数五十九页\编于十七点式中：pa为大气压，K、n为试验常数，§3.4弹塑性模型简介

一、模型

对不同的土类，K值可能小于100，也可能大于3500；n值一般在0.2~1.0之间。

因为切线模量根据式中：S为应力水平，反映强度发挥的程度将此式回代得目前九页\总数五十九页\编于十七点§3.4弹塑性模型简介

一、模型

目前十页\总数五十九页\编于十七点

根据莫尔-库仑破坏准则，有§3.4弹塑性模型简介

一、模型

目前十一页\总数五十九页\编于十七点此时:可见，a，b均是的函数。由求解的公式可知，该模量与所研究土体的c、φ强度指标，与K、n、Rf有关，是固结压力σ3的函数。此时，弹性常数的线性表达式已不适用，而应采用受荷过程中的增量弹性模量表达式。§3.4弹塑性模型简介

一、模型

其中的K、n如何确定呢？目前十二页\总数五十九页\编于十七点关系曲线

~两侧取对数，即使幂函数转换为直线函数，将公式直线斜率为n，于是切线弹性模量可求。§3.4弹塑性模型简介

一、模型

K、n的确定方法如图所示此时图中直线截距为logK，n目前十三页\总数五十九页\编于十七点二、切线泊桑比泊桑比是在受荷过程中，土体侧向应变εr与轴向应变εa的比值。类同于弹性模量的方法，当侧向应变趋于无穷小时，其切线泊桑比μi为初始泊桑比，μt为切线泊桑比。库哈威（Kulhawy）与邓肯将三轴试验中εr

、εa数据用双曲线函数表达，此关系曲线如图所示。-关系曲线~——轴向应变。——侧向应变。§3.4弹塑性模型简介

一、模型

s3=常数目前十四页\总数五十九页\编于十七点初始泊桑比是固结压力σ3的函数，其表达式为式中：G、F——为试验常数。切线泊桑比公式为式中：A为试验常数。§3.4弹塑性模型简介

一、模型

如仿照前面的作法时,则为图中直线的截距与斜率。目前十五页\总数五十九页\编于十七点三、回弹模量动荷载作用下，土体由加荷到卸荷过程中，（-）~

曲线如图所示。OA为加荷状态的应力应变关系曲线，其斜率为；AB为卸荷曲线，其斜率为。邓肯等人假定：不随（）变化，仅随而变。

§3.4弹塑性模型简介

一、模型

加荷与卸荷曲线

AB加荷卸荷目前十六页\总数五十九页\编于十七点回弹模量可由下式计算：

关系曲线§3.4弹塑性模型简介

一、模型

在双对数纸上点绘关系曲线，得一直线，如图所示，其截距为，斜率为n。n目前十七页\总数五十九页\编于十七点邓肯-张模型采用非线性弹性常数，反映了土体变形的主要规律，其双曲线函数由常规三轴试验确定、参数少，适用于粘性土、砂土及σ3接近常数的岩土工程问题。不宜用于密砂及严重超固结土，它不能考虑应力路径对该模量的影响。2003年以来，国内河海大学等单位，对常规三轴仪进行改进，完成了轴向卸载、侧向加载及侧向卸载不同应力路径下切线模量公式的推导工作，对双曲线模型进一步给予完善，并通过试验进行了验证。§3.4弹塑性模型简介

一、模型

目前十八页\总数五十九页\编于十七点◆侧向卸载条件下的切线弹性模量式中：——轴向固结压力，——侧向应力式中，、分别为轴向、径向应力，为侧向固结压力。§3.4弹塑性模型简介

一、模型

目前十九页\总数五十九页\编于十七点◆侧向加载条件下的切线弹性模量Et式中：c、φ来自侧向加载试验。同样，若假设不同应力路径的强度指标一致，Ko固结侧向加载条件下的切线弹性模量公式如下：式中各符号意义同前其中Rf仍为§3.4弹塑性模型简介

一、模型

目前二十页\总数五十九页\编于十七点◆轴向卸载条件下的切线弹性模量Et式中：c、φ来自轴向卸载试验。若假设不同应力路径的强度指标一致，对于K0固结，上式可简化为：式中各符号意义同前其中Rf仍为§3.4弹塑性模型简介

一、模型

目前二十一页\总数五十九页\编于十七点§3.4弹塑性模型简介

一、模型

以上的研究工作，为邓肯——张模型在基坑支护工程等问题中的应用创造了良好条件。目前二十二页\总数五十九页\编于十七点K~G模型

目前二十三页\总数五十九页\编于十七点K~G模型

该模型又称为（Domaschuk-Vallianppam）模型，它以弹性体积模量K、剪切模量G为参数，以半对数曲线整理三向固结压缩试验，双曲线整理三向等压固结排水压缩试验结果，建立K、G参数的本构模型。土样在（各向等压）静水压力的条件下各向等压固结，得到孔隙比与压力的曲线，该曲线用半对数坐标系表示为体积模量K值的测定目前二十四页\总数五十九页\编于十七点λ——试验常数

已知土样体积应变，则对该式两侧求导体积模量§3.4弹塑性模型简介二、模型

式中：——初始孔隙比，目前二十五页\总数五十九页\编于十七点§3.4弹塑性模型简介二、模型

不同压力p下的Kt可求。三向等压固结的试验曲线由于土体应力历史的存在，因而定义前期固结压力pc与对应的体应变之比为初始体积模量Ki：

λ由下图中直线斜率（即压缩指数）确定目前二十六页\总数五十九页\编于十七点§3.4弹塑性模型简介二、模型

此时，切线体积模量为Kt式中：n为试验常数。在水利电力部土工试验规程（1979）中建议半对数坐标系内直线方程为：对左式微分，得：目前二十七页\总数五十九页\编于十七点§3.4弹塑性模型简介二、模型

式中：为p应力范围内的压缩指数。剪切模量G值的测定

八面体应力：与其相适应的应变为体积应变：目前二十八页\总数五十九页\编于十七点§3.4弹塑性模型简介二、模型

◆复杂受力状态下的剪切应变为：八面体剪应变的几何意义是:八面体剪应力的指向与八面体法线所交直角的改变量。工程中为便于描述塑性现象，经常将八面体剪应变乘以，得到广义应变的虚拟值。目前二十九页\总数五十九页\编于十七点然后，再按pc等于另一不同值的土样求取另一组曲线G值的测定采用p=常数的三向等压固结排水压缩试验，土样先在三向等压下固结，然后在固结比1.0、0.8、0.6、0.4、0.2五种不同情况下，进行三向等压固结排水试验，求取一组（5条）应力-应变曲线。§3.4弹塑性模型简介二、模型

目前三十页\总数五十九页\编于十七点所有的曲线都用考特纳(Kondner)建议的双曲线关系式：将其微分，得§3.4弹塑性模型简介二、模型

目前三十一页\总数五十九页\编于十七点取：Gi=1/a为初始剪切模量；为应力偏张量最大值。引入剪应力的破坏比R§3.4弹塑性模型简介二、模型

该方程即是切线剪切模量关系式。目前三十二页\总数五十九页\编于十七点

试验证明，切线剪切模量与破坏比Rf，初始孔隙比e0，前期固结压力pc有关，其表达式为：式中：α、β为试验常数§3.4弹塑性模型简介二、模型

根据水利电力部土工试验规程，Gt简化求解时，土体中八面体偏应力q与相应的剪切应变仍采用双曲线方程表示为：目前三十三页\总数五十九页\编于十七点式中：与前面(Domaschuk)的推导相似，得出由三轴等向固结排水试验得出式中：§3.4弹塑性模型简介二、模型

目前三十四页\总数五十九页\编于十七点§3.4弹塑性模型简介二、模型

由于E~μ、K~G，模型是建立在广义虎克定律的基础上，由常规三轴试验资料得到了弹性常数与应力状态的非线性关系表达式，进行有限元计算，因此这两个模型是非线性弹性模型。目前三十五页\总数五十九页\编于十七点（剑桥模型）CambridgeModel目前三十六页\总数五十九页\编于十七点剑桥模型的中心是状态边界面，是英国剑桥大学的罗斯科（，1958~1963)等人为正常固结和弱超固结的重塑粘土（也就是所谓的“湿粘土”）创制的弹、塑性应力——应变关系的数学模型，是一种“帽子”模型。外荷作用下土体发生变形，随荷载增加土体变形增大；当土体由静止状态发生失稳或变形过大不能保证上部结构正常工作时，土体达到极限破坏，将不能连续承受荷载，这就是“临界状态”。临界状态§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）◆临界状态线和状态边界面目前三十七页\总数五十九页\编于十七点通过大量三轴排水与固结不排水实验证明，施加偏应力，在p~q平面内，使土体产生破坏的应力点均落在同一直线OD上如图，其方程为：在平面状态下,若取为横坐标，取为纵坐标，建立坐标系。pq0q～p平面内的临界状态线D§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）目前三十八页\总数五十九页\编于十七点在e~p平面内，破坏点均落在同一曲线CD上，如图pe0初始三向等压固结线CDABe～p平面内临界状态线e=ec-λ.lnp在e~p~q三维空间内，初始等向压缩曲线必在e~p面内，设其为曲线AB，随着偏应力的施加，土体中剪应力增大，孔隙比改变，曲线在三维空间坐标系中脱离原水平面向上方移动，达到破坏时，对应的空间曲线叫临界状态线。初始状态时对土样进行等向压缩，此时q=0,§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）目前三十九页\总数五十九页\编于十七点初始三向等压固结线临界状态线e～p平面上临界状态线投影e=ec-λ.lnpq～p平面上临界状态线投影q=M.pABDC因此，临界状态线是破坏点的空间轨迹，其在平面内投影为一直线，简称破坏线。qpe0屈服曲面§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）目前四十页\总数五十九页\编于十七点在三维p-q-e空间坐标系内，对应每一个p、q应力状态，均有相应的孔隙比e值，这样p、q、e三个坐标就定义了一个空间点。在该面以下为弹性区，该面以上是土体不可能达到的区域，沿着状态边界面将发生塑性变形。随着应力状态的变化，大量的空间点移动形成一个从AB到CD的空间曲面，这个空间曲面叫状态边界面。初始三向等压固结线临界状态线ABDCqpe0屈服曲面目前四十一页\总数五十九页\编于十七点当已知p、q、e时，即可确定空间点的位置，其位置在状态面以下为弹性状态，处于状态边界面上即为屈服，根据该曲线的数学表达式，可建立屈服准则。图中AB是状态边界面与eop平面的交线，初始三向等压固结线e在卸荷过程中（例基坑内土体开挖），图中EH为等向膨胀曲线，膨胀曲线e=ek-k.lnpE偏应力增大到极限时，应力状态到达临界状态线CD。AB是等向压缩曲线（q=0）初始三向等压固结线临界状态线ABDCqpe0屈服曲面H目前四十二页\总数五十九页\编于十七点初始三向等压固结线临界状态线q～p平面上临界状态线投影q=M.pABDCqpe0屈服曲面膨胀曲线e=ek-k.lnpEHGNF应力空间中状态边界面EG线上的N点，它在侧面p-q坐标系内投影为N′点,也是屈服轨迹E′G′上的N′点。在EG上与EH曲线间组成的平行于q轴的竖向平面内，均为弹性变形。（因为假定塑性势面与屈服面重合），在p-q平面上的屈服轨迹即为塑性势线。屈服轨迹在p-q平面上的投影塑性势线目前四十三页\总数五十九页\编于十七点1956年起

等人为正常固结和弱超固结的重塑粘土（也就是所谓的伦敦“湿粘土”）创制的弹、塑性应力——应变关系的数学模型，◆

剑桥模型屈服函数的推导:按照能量理论：土体在p、q作用下发生了应变增量:dεv和dg

，相应的能量变化为：dW=p.dεv+q.dγ，这些能量等于储存在土体内可恢复的弹性能dWe和不可恢复的耗散能（塑性能）dWp之和目前四十四页\总数五十九页\编于十七点即：由于Roscoe假定§3.4弹塑性模型简介三、剑桥模型（CambridgeModel），所以从正常固结线上卸荷，使其回弹，在该区域内，单元体在p、q作用下发生弹性应变增量，该增量值可用e～lnp关系曲线的膨胀段来代替。各向等压固结试验膨胀段（回弹段）的曲线方程为：e=ek-k.lnp，<1>k为e~lnp关系曲线膨胀段的斜率。假设一切剪应变都是不可恢复的弹性剪应变为“零”目前四十五页\总数五十九页\编于十七点式中ek

、k均为试验常数，微分后得弹性应变孔隙比的变化量，并用角标“e”表示弹性变化，则若体积应变是用孔隙比定义的，<2>eo——为初始孔隙比，是常数。则：<3>将<2>式代入<3>后得：<4>目前四十六页\总数五十九页\编于十七点将<4>式代入<1>后得：剑桥模型又假设在同一屈服轨迹上即为常数，所以前面已讲过：q=M.p沿着屈服面移动时，塑性势均相等也就是说，剑桥模型认为塑性能等于塑性势面与屈服面重合目前四十七页\总数五十九页\编于十七点所以在屈服轨迹（塑性势）F（p，q，k)=0上，其屈服参数k为常数，则由正交定律得<5><6><7>代入<6>得：<7><8>§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）1928年，米色斯提出塑性应变增量的方向与塑性势面的法线方向一致，即与塑性势面正交目前四十八页\总数五十九页\编于十七点根据代入<8>式<9>根据<5>式得<10><10>代入<9>得§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）目前四十九页\总数五十九页\编于十七点因为所以即解该微分方程，得解如下：其中，c为积分常数设该屈服轨迹与坐标轴p的交点为A，则A点坐标为(po,0，0)<11><12>§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）目前五十页\总数五十九页\编于十七点代入后得所以<13>即也就是说，相应的屈服函数F为<14>若设该屈服函数与破坏线交于（px，qx，ex或vx），则§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）屈服函数!!目前五十一页\总数五十九页\编于十七点相应的屈服函数为因为所以<15><15>式又可写成<16>po——屈服轨迹与p轴交点处的p值（硬化参数）px——屈服轨迹与破坏线交点处的p值屈服轨迹的完整形状是一个对称的闭合曲线，象一个弹头形状，随着压力的增加而扩展。如前所述，正常固结粘土的状态边界面是原始固结曲线和临界状态线围成的部分，所以，根据状态边界面上任一点的屈服轨迹和正常固结线、膨胀曲线可以推导状态§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）目前五十二页\总数五十九页\编于十七点边界面方程；如屈服轨迹沿正常固结线与临界状态线移动所发生的曲面是状态边界面，如图qpe初始三向等压固结线临界状态线膨胀曲线e=ek-k.lnpe～p平面上临界状态线投影e=ec-λ.lnpq～p平面上临界状态线投影q=M.p屈服曲面B§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）目前五十三页\总数五十九页\编于十七点对于图中的B点(正常固结线、膨胀曲线和屈服曲面的公共点，且该点在e～p平面内），设其坐标为p=po，e=eo，正常固结线上q=0，则过B点的正常固结线与膨胀曲线为，正常固结线回弹曲线式中的ec、ek分别为正常固结曲线、膨胀曲线上p=1.0时的孔隙比，<17><18><17>-<18><19>将其代入到膨胀曲线的一般方程e-ek+k.lnp=0中得，<20>§3.4弹塑性模型简介三、剑桥模型（CambridgeModel）目前五十四页\总数五十九页\编于十七点与相应的屈服轨迹方程联立消去po得到状态边界面方程如下：即<21