弹塑性问题的有限单元法PPT讲稿.ppt

弹塑性问题的有限单元法第1页，共97页，编辑于2022年，星期日线弹性材料（线弹性材料（Elastic Material）n在外力作用下，材料的应力应变呈线性关系，并且在某一应力在外力作用下，材料的应力应变呈线性关系，并且在某一应力水平卸载后，材料恢复到原来的状态，即变形为零。水平卸载后，材料恢复到原来的状态，即变形为零。Elastic implies that the element will return to its original size and shape once all applied loads are removed,if not stressed past its yield point(y).第2页，共97页，编辑于2022年，星期日 岩土体材料岩土体材料（Rock and SoilRock and Soil）n岩土体不是完全的线弹性体，应力应变很少呈线性关系，且当外力岩土体不是完全的线弹性体，应力应变很少呈线性关系，且当外力消除以后，往往有不可恢复的塑性变形。因此，本章考虑材料的非消除以后，往往有不可恢复的塑性变形。因此，本章考虑材料的非线性，主要讨论岩土体的弹塑性模型。线性，主要讨论岩土体的弹塑性模型。If stressed past its yield point(plastic range),the material will take a permanent set and will not fully return to its original size and shape,when all loads are removed.第3页，共97页，编辑于2022年，星期日第一节第一节 应力与应变分析应力与应变分析一、一点的应力状态与应力张量一、一点的应力状态与应力张量二、应力张量的分解二、应力张量的分解三、应力不变量三、应力不变量四、主应力空间、四、主应力空间、平面与罗德（平面与罗德（LodeLode）角）角五、应变张量五、应变张量主要内容：主要内容：第4页，共97页，编辑于2022年，星期日一、一点的应力状态与应力张量一、一点的应力状态与应力张量n应力张量可以采用分量记法，即用矩应力张量可以采用分量记法，即用矩阵或张量下标表示阵或张量下标表示n 任一点的任一点的6个应力分量：个应力分量：=x，y，z，xy，xz，yzT（3-1）式中 i,j=x,y,z第5页，共97页，编辑于2022年，星期日张量张量引入张量：引入张量：0 0 0 0阶张量：阶张量：阶张量：阶张量：3 3 3 30 0 0 01 11 1阶张量：阶张量：阶张量：阶张量：3 31 1 1 13 32 2 2 2阶张量：阶张量：3 32 2 2 29 93 3 3 3阶张量：阶张量：3 33 3 3 32727应力和应变是二阶张量应力和应变是二阶张量标量标量标量标量矢量矢量矢量矢量第6页，共97页，编辑于2022年，星期日用二阶张量在用二阶张量在用二阶张量在用二阶张量在x,y,z 坐标系表示坐标系表示或写成：或写成：第7页，共97页，编辑于2022年，星期日矩阵记法矩阵记法（3-1）张量下标记法张量下标记法式中式中 i,j=x,y,zi,j=x,y,z 应当注意应当注意：张量分量的矩阵记法与二维数组的区别。：张量分量的矩阵记法与二维数组的区别。应力张量可以用其分量表示成应力张量可以用其分量表示成33的对称方阵，而写成矩的对称方阵，而写成矩阵形式阵形式33的二维数组却不一定是张量，更不一定是应力张量，的二维数组却不一定是张量，更不一定是应力张量，只不过是按一定顺序排列的数组而已。只不过是按一定顺序排列的数组而已。第8页，共97页，编辑于2022年，星期日二、应力张量的分解二、应力张量的分解 张量可以合成和分解张量可以合成和分解n应力张量应力张量 分解为分解为 球张量球张量 和和 偏张量偏张量=+=+第9页，共97页，编辑于2022年，星期日=+（3-2）n物理意义物理意义n应力张量的球张量分量应力张量的球张量分量作用在该点的平均应力作用在该点的平均应力(m)或静水压力或静水压力(p)n应力张量的偏斜分量应力张量的偏斜分量作用在该点的偏应力作用在该点的偏应力(si)和剪和剪应力应力(sij)(ij)。第10页，共97页，编辑于2022年，星期日张量分解在塑性理论中的意义张量分解在塑性理论中的意义n在弹性力学中，应力球张量只产生弹性应变在弹性力学中，应力球张量只产生弹性应变(应变球张量应变球张量)，应力偏，应力偏张量只产生弹性剪应变张量只产生弹性剪应变(应变偏张量应变偏张量)，本构关系非常简单。，本构关系非常简单。n在金属塑性理论中假设体应变为弹性的，故体应变只有弹性分量而与塑性在金属塑性理论中假设体应变为弹性的，故体应变只有弹性分量而与塑性无关，剪应变有塑性分量。无关，剪应变有塑性分量。n将应力分解为球张量与偏张量，不仅使它们与体应变和剪应变之将应力分解为球张量与偏张量，不仅使它们与体应变和剪应变之间的关系相互对应，而且可以简化本构关系的分析。间的关系相互对应，而且可以简化本构关系的分析。n即使在应力球张量与偏应变、应力偏张量与体应变发生耦合作用的即使在应力球张量与偏应变、应力偏张量与体应变发生耦合作用的岩土塑性本构关系理论中，将应力分解为球张量与偏张量，也便于岩土塑性本构关系理论中，将应力分解为球张量与偏张量，也便于分析它们对塑性体应变与剪应变的各自影响，从而建立相应的本构分析它们对塑性体应变与剪应变的各自影响，从而建立相应的本构关系。关系。第11页，共97页，编辑于2022年，星期日三、应力不变量三、应力不变量n应力张量的分量表示法与坐标轴的选取有关，当应力张量的分量表示法与坐标轴的选取有关，当进行一点应力状态分析或建立弹塑性本构关系时，进行一点应力状态分析或建立弹塑性本构关系时，如果能够找到与坐标系选取无关的应力不变量来如果能够找到与坐标系选取无关的应力不变量来表示，则将会简捷得多。表示，则将会简捷得多。n下面就寻找这样的不变量。下面就寻找这样的不变量。第12页，共97页，编辑于2022年，星期日（一）应力张量不变量（一）应力张量不变量n在弹性力学中已经证明，通过一点可以找到相互在弹性力学中已经证明，通过一点可以找到相互垂直的三个主平面，在这些面上的剪应力为零。垂直的三个主平面，在这些面上的剪应力为零。n主平面上作用的正应力就称为主平面上作用的正应力就称为 主应力主应力。主平面。主平面的方向称为的方向称为 主方向主方向。对于一定的应力状态而言，。对于一定的应力状态而言，主应力和主方向是不变的。主应力和主方向是不变的。第13页，共97页，编辑于2022年，星期日现在直角坐标系中取一四面体，设斜面现在直角坐标系中取一四面体，设斜面ABCABC是主平面，则作用于该面的法向是主平面，则作用于该面的法向应力即为主应力，如图应力即为主应力，如图3-23-2所示。所示。（3-3）又设又设l、m、n代表主应力的方向余弦，代表主应力的方向余弦，l=cos(x，N)，m=cos(y，N)，n=cos(z，N)则主应力则主应力在三个坐标轴上的投影为在三个坐标轴上的投影为px=l,py=m,pz=n平衡方程ABOxyz N NC图图3-2 3-2 四面体应力四面体应力第14页，共97页，编辑于2022年，星期日n要使要使l l、m m、n n有非零解，则必有有非零解，则必有（3-5）第15页，共97页，编辑于2022年，星期日n式式(3-5)(3-5)的三个根即为主应力的三个根即为主应力1 1、2 2、3 3。n该点的主应力值不会因坐标选择而改变，因而该点的主应力值不会因坐标选择而改变，因而I I1 1、I I2 2、I I3 3的值也是不会改变的，它们分别称为应力张量的第一、的值也是不会改变的，它们分别称为应力张量的第一、第二、第三不变量第二、第三不变量。n如果坐标轴方向与主应力方向一致，则（如果坐标轴方向与主应力方向一致，则（无剪应力无剪应力）I I1 1=1 1+2 2+3 3 I I2 2=-(=-(1 12 2+2 23 3+3 31 1)I I3 3=1 12 23 3 (3-6)(3-6)（3-5）第16页，共97页，编辑于2022年，星期日（二）应力偏张量不变量（二）应力偏张量不变量n应力偏张量应力偏张量S Sijij也有三个不变量也有三个不变量J J1 1=s=sx x+s+sy y+s+sz z=0=0J J2 2=-(s=-(sx xs sy y+s+sy ys sz z+s+sz zs sx x)+s)+s2 2xyxy+s+s2 2yzyz+s+s2 2zxzx =1/6 (=1/6 (x x-y y)2 2+(+(y y-z z)2 2+(+(z z-x x)2 2+6(+6(2 2xyxy+2 2yzyz+2 2zxzx =1/2 (s =1/2 (s2 2x x+s+s2 2y y+s+s2 2z z)+s)+s2 2xyxy+s+s2 2yzyz+s+s2 2zxzxJ J3 3=s=sx xs sy ys sz z+2s+2sxyxys syzyzs szxzx-s-sx xs s2 2yzyz-s-sy ys s2 2zxzx-s-sz zs s2 2xyxy (3-7)(3-7)n式中式中 J J1 1、J J2 2、J J3 3分别称为应力偏张量的第一、第二、第三分别称为应力偏张量的第一、第二、第三不变量。不变量。第17页，共97页，编辑于2022年，星期日n当取坐标轴与主轴方向一致时，则当取坐标轴与主轴方向一致时，则 J J1 1=0=0 J J2 2=-(s=-(s1 1s s2 2+s+s2 2s s3 3+s+s3 3s s1 1)=1/6 (=1/6 (1 1-2 2)2 2+(+(2 2-3 3)2 2+(+(3 3-1 1)2 2 J J3 3=s=s1 1s s2 2s s3 3 第18页，共97页，编辑于2022年，星期日n应力张量不变量应力张量不变量I1、I2、I3n应力偏张量不变量应力偏张量不变量J2、J3、(J1=0)n在研究岩土弹塑性问题时我们更关注其中的在研究岩土弹塑性问题时我们更关注其中的I1、J2、J3n这是因为这是因为I1只与平均应力只与平均应力m或静水应力或静水应力p有关，有关，而而J2反映剪应力的大小，反映剪应力的大小，J3表示剪应力的方向。表示剪应力的方向。第19页，共97页，编辑于2022年，星期日四、主应力空间、罗德角四、主应力空间、罗德角n假设：岩土体为各向同性体，因此主应力的作用方向就假设：岩土体为各向同性体，因此主应力的作用方向就无关紧要，通常只要研究作用在一点的主应力大小就可无关紧要，通常只要研究作用在一点的主应力大小就可以了。以了。n三个主应力正好可以用三维空间来直观地描述。三个主应力正好可以用三维空间来直观地描述。n以三个主应力为轴而组成的以三个主应力为轴而组成的笛卡尔空间坐标系笛卡尔空间坐标系就称为就称为主主应力空间应力空间，如图所示。，如图所示。第20页，共97页，编辑于2022年，星期日OQQ以主应力表示的物体中一点的应力状态在主应力空间中对应一个点 Q Q(1 1，2 2，3 3)。n原点O与Q的连线OQOQ称为该点的应力矢量，它代表着岩土体中相应点的应力大小与方向。n在主应力空间中，与三个坐标轴成相等倾角的线称为在主应力空间中，与三个坐标轴成相等倾角的线称为线线(等压等压线）。线）。线的方程可以表示为线的方程可以表示为n 1 1=2 2=3 3 (3-8)(3-8)第21页，共97页，编辑于2022年，星期日OQQ偏平面偏平面()()的方程为的方程为式中 偏平面与原点的距离(3-9)而平面的方程为为了确定偏剪应力的方向引入罗德角的概念。MO偏剪应力与偏剪应力与OMOM线的夹角就定义为线的夹角就定义为罗德角，规定顺时针（罗德角，规定顺时针（-），逆时），逆时针（针（+）。这样）。这样就代表偏剪应就代表偏剪应力在偏平面上的作用方向。力在偏平面上的作用方向。第22页，共97页，编辑于2022年，星期日n 与等压线相正交的平面称为偏平面，通过坐标原点与等压线相与等压线相正交的平面称为偏平面，通过坐标原点与等压线相正交的平面称为正交的平面称为平面。可见平面。可见平面是一个特殊的偏平面。由偏平面的定平面是一个特殊的偏平面。由偏平面的定义可知，在一个偏平面内平均应力为常量，故偏平面的方程为：义可知，在一个偏平面内平均应力为常量，故偏平面的方程为：（3-9）第23页，共97页，编辑于2022年，星期日OQQn将应力空间中代表一点应力状将应力空间中代表一点应力状态的应力分量态的应力分量OQ向偏平面和向偏平面和等压线上投影，即可得到作用等压线上投影，即可得到作用在偏平面上的正应力分量在偏平面上的正应力分量和和偏剪应力分量偏剪应力分量r 所以应力空所以应力空间中任一点的应力状态可以间中任一点的应力状态可以用偏平面上的应力来表示。用偏平面上的应力来表示。第24页，共97页，编辑于2022年，星期日可以证明可以证明(3-11)n同样可以证明正应力和应力不变量同样可以证明正应力和应力不变量m、J2及及的关系的关系为：为：第25页，共97页，编辑于2022年，星期日五、应变张量五、应变张量n对于连续变形体来说，应力与应变是对偶的，对于连续变形体来说，应力与应变是对偶的，即有应力就必然产生应变，有应变就必然产生即有应力就必然产生应变，有应变就必然产生应力。应力。n应力和应变都属于二阶对称张量。应力和应变都属于二阶对称张量。n从张量的角度分析，它们有许多相似之处。例如从张量的角度分析，它们有许多相似之处。例如它们都可以分解为球张量与偏张量，都具有不变它们都可以分解为球张量与偏张量，都具有不变量等。量等。第26页，共97页，编辑于2022年，星期日一点的应变状态可以用几个应变分量描述，一般可以用应变张量表示：第27页，共97页，编辑于2022年，星期日3000.025 模拟软、硬顶板下开采煤层底板的变形和受力情况模拟软、硬顶板下开采煤层底板的变形和受力情况两侧各留80m，开采120m，分6步开挖完毕实例分析（用实例分析（用ADINA软件）软件）第28页，共97页，编辑于2022年，星期日岩性岩性E（MPa）（kg/m3）（0）C（MPa）t（MPa）硬硬顶顶356000.2026504112.706.21软顶软顶27000.232100408.131.80软软煤煤15000.381300201.200.64软软底底29000.242600408.621.93第29页，共97页，编辑于2022年，星期日第30页，共97页，编辑于2022年，星期日（a）不同不同顶顶板岩性板岩性对对底板底板竖竖直直应应力的影响力的影响 第31页，共97页，编辑于2022年，星期日（c）不同）不同顶顶板岩性板岩性对对底板最大剪底板最大剪应应力的影响力的影响 第32页，共97页，编辑于2022年，星期日第二节第二节 岩土介质弹塑性本构关系岩土介质弹塑性本构关系一、岩土介质本构关系基本类型一、岩土介质本构关系基本类型二、增量塑性理论简介二、增量塑性理论简介三、常用的弹塑性模型三、常用的弹塑性模型主要内容：主要内容：第33页，共97页，编辑于2022年，星期日一、岩土介质本构关系基本类型一、岩土介质本构关系基本类型n本构关系的基本概念：本构关系的基本概念：n岩土介质在外力（广义）作用下，产生应力与应变。岩土介质在外力（广义）作用下，产生应力与应变。n普遍的应力普遍的应力应变关系即称为应变关系即称为本构关系本构关系。n材料的本构关系可以表示为材料的本构关系可以表示为 ij=f(ij,t,T,应力历史，等应力历史，等)(3-16)式中式中 t加载历时；加载历时；T温度温度第34页，共97页，编辑于2022年，星期日 弹性力学中的广义胡克定律就是最简单的材料本构关系，它不弹性力学中的广义胡克定律就是最简单的材料本构关系，它不计时间、温度、应力历史和应力路径，应力和应变之间存在唯一的计时间、温度、应力历史和应力路径，应力和应变之间存在唯一的对应关系。当材料应力超出弹性范围而进入塑性阶段时，应力和应对应关系。当材料应力超出弹性范围而进入塑性阶段时，应力和应变之间就没有唯一的对应关系，而受应力历史或应力路径的影响，变之间就没有唯一的对应关系，而受应力历史或应力路径的影响，这时材料的应力应变关系就称为塑性本构关系。这时材料的应力应变关系就称为塑性本构关系。塑性本构关系要比弹性本构关系复杂得多，如果再考虑材塑性本构关系要比弹性本构关系复杂得多，如果再考虑材料应力应变关系随时间和温度变化，则其本构关系将更复杂。料应力应变关系随时间和温度变化，则其本构关系将更复杂。第35页，共97页，编辑于2022年，星期日ijij=f(=f(ijij,t,T,t,T,应力历史，等应力历史，等)弹塑性本构关系弹塑性本构关系弹性本构关系弹性本构关系第36页，共97页，编辑于2022年，星期日岩土介质应力应变关系岩土介质应力应变关系典型试验曲线典型试验曲线(1)线弹性(2)非线性弹性(3)线性弹塑性(4)非线性弹塑性tttt(5)粘弹性硬化硬化软化软化理想塑性理想塑性应变软化应变软化应变硬化应变硬化蠕变蠕变松弛松弛要理解掌握要理解掌握第37页，共97页，编辑于2022年，星期日二、增量塑性理论简介二、增量塑性理论简介（一）屈服条件（一）屈服条件（二）加载条件（二）加载条件（三）加载和卸载准则（三）加载和卸载准则（四）流动法则（四）流动法则（五）一般的弹塑性本构关系（五）一般的弹塑性本构关系第38页，共97页，编辑于2022年，星期日(一一)屈服条件屈服条件n屈服条件屈服条件是表示在复杂应力状态下，材料进入初始屈服是表示在复杂应力状态下，材料进入初始屈服时应力分量之间所必须满足的条件。时应力分量之间所必须满足的条件。n如果以应力如果以应力作为坐标轴，屈服条件用作为坐标轴，屈服条件用f()=0表示，表示，则应力空间中则应力空间中f=0的一个曲面称为的一个曲面称为屈服曲面屈服曲面。n当应力点当应力点位于曲面之内位于曲面之内(f0)，材料处于弹性状态；当，材料处于弹性状态；当应力点位于曲面上应力点位于曲面上(f=0)，材料开始屈服进入塑性状态。，材料开始屈服进入塑性状态。第39页，共97页，编辑于2022年，星期日n这样可以分清平均应力这样可以分清平均应力(或静水压力或静水压力)与偏应力与偏应力(或剪应或剪应力力)对体应变与偏应变的贡献。对体应变与偏应变的贡献。n屈服条件也可以表示为主应变或应变不变量的函数，由屈服条件也可以表示为主应变或应变不变量的函数，由于两种形式通过本构关系可以互换，因此以下将采用应于两种形式通过本构关系可以互换，因此以下将采用应力屈服函数的形式。力屈服函数的形式。n假定材料是各向同性的，屈服条件将与坐标轴的选取无关，假定材料是各向同性的，屈服条件将与坐标轴的选取无关，因此可以表示成只是应力不变量的函数，例如：因此可以表示成只是应力不变量的函数，例如：第40页，共97页，编辑于2022年，星期日岩土介质常用的屈服条件：n1.莫尔库仑（Mohr-Coloumb）屈服条件n2.德鲁克普拉格（Drucker-Prager）屈服条件第41页，共97页，编辑于2022年，星期日1.莫尔莫尔库仑库仑(M-C)屈服条件屈服条件n莫尔库仑屈服条件可以表示为n=c+ntan (3-18)n式中、n破坏面上的剪应力与法向正应力；c、材料的粘聚力和内摩擦角。第42页，共97页，编辑于2022年，星期日=c+ntg (3-18)第43页，共97页，编辑于2022年，星期日图 屈服面与屈服轨迹(a)空间屈服面；(b)偏平面屈服迹线莫尔库仑屈服面六棱锥面，棱边上导数的方向不定，在计算中带来不便。为此，德鲁克和普拉格对其进行了改进，提出用一个内切于莫尔提出用一个内切于莫尔库仑条件六棱锥面的圆锥面来作为屈服面库仑条件六棱锥面的圆锥面来作为屈服面第44页，共97页，编辑于2022年，星期日图 屈服面与屈服轨迹(a)空间屈服面；(b)偏平面屈服迹线第45页，共97页，编辑于2022年，星期日2.德鲁克德鲁克普拉格普拉格(D-P)屈服条件屈服条件德鲁克普拉格屈服函数可以表示为:第46页，共97页，编辑于2022年，星期日广义广义Mises屈服条件屈服条件第47页，共97页，编辑于2022年，星期日 综合分析岩土数值分析中常用的莫尔综合分析岩土数值分析中常用的莫尔-库仑屈服条件和德鲁克库仑屈服条件和德鲁克-普拉格屈服条件的相关函数表达式、数值模拟中需要的主要物理力学参数普拉格屈服条件的相关函数表达式、数值模拟中需要的主要物理力学参数和在主应力空间中的物理意义（结合相关屈服条件图进行综合分析）。和在主应力空间中的物理意义（结合相关屈服条件图进行综合分析）。3.M-C与与D-P屈服条件的比较屈服条件的比较n 屈服条件是表示在复杂应力状态下，材料进入初始屈服时应力分量屈服条件是表示在复杂应力状态下，材料进入初始屈服时应力分量之间所必须满足的条件。如果以应力之间所必须满足的条件。如果以应力作为坐标轴，屈服条件用作为坐标轴，屈服条件用f()=0表示，则应力空间中表示，则应力空间中f=0的一个曲面称为屈服曲面。当应力点的一个曲面称为屈服曲面。当应力点位于曲面之内位于曲面之内(f0)，材料处于弹性状态；当应力点位于曲面上，材料处于弹性状态；当应力点位于曲面上(f=0)，材料开始屈服进入塑性状态。，材料开始屈服进入塑性状态。第48页，共97页，编辑于2022年，星期日莫莫尔尔库仑库仑屈服条件可以表示屈服条件可以表示为为：（1）莫尔）莫尔库仑库仑(M-C)屈服条件屈服条件第49页，共97页，编辑于2022年，星期日第50页，共97页，编辑于2022年，星期日（2）德鲁克）德鲁克普拉格普拉格(D-P)屈服条件屈服条件第51页，共97页，编辑于2022年，星期日第52页，共97页，编辑于2022年，星期日第53页，共97页，编辑于2022年，星期日第54页，共97页，编辑于2022年，星期日FLAC3D中的本构模型第55页，共97页，编辑于2022年，星期日（二）加载条件（二）加载条件n在复杂应力状态下，材料进入塑性状态后卸载，当再加载时，在复杂应力状态下，材料进入塑性状态后卸载，当再加载时，屈服函数会随着以前发生过的塑性变形的历史而有所改变。屈服函数会随着以前发生过的塑性变形的历史而有所改变。n当应力分量满足某一关系时，材料将重新进入塑性状态当应力分量满足某一关系时，材料将重新进入塑性状态而产生新的塑性变形，这种现象叫做而产生新的塑性变形，这种现象叫做加载屈服加载屈服。n材料在初始屈服以后再进入塑性状态时应力分量间所必材料在初始屈服以后再进入塑性状态时应力分量间所必须满足的函数叫做须满足的函数叫做加载条件加载条件，用，用=0来表示，在应力空间来表示，在应力空间中相应的曲面称为中相应的曲面称为加载面加载面。对于理想塑性材料，加载面。对于理想塑性材料，加载面等于屈服面。等于屈服面。第56页，共97页，编辑于2022年，星期日加载函数加载函数n虽然仍可以用应力或应力不变量作变量，但还依赖于塑性应变的过程，一般可以表示为式中：H硬化参量。n常用的加载面的二种模型：(1)等向硬化模型 (2)随动硬化模型第57页，共97页，编辑于2022年，星期日(1)等向硬化模型，认为加载面在应力等向硬化模型，认为加载面在应力空间中的形状和中心位置保持不变，空间中的形状和中心位置保持不变，随着强化程度的增加，初始屈服面随着强化程度的增加，初始屈服面作形状相似的扩大。如图所示。作形状相似的扩大。如图所示。(2)随动硬化模型，认为在塑性变随动硬化模型，认为在塑性变形过程中，加载面的大小、形状形过程中，加载面的大小、形状与初始屈服面相同，只是在应力与初始屈服面相同，只是在应力空间中沿塑性变形方向作刚体移空间中沿塑性变形方向作刚体移动，如图所示。动，如图所示。第58页，共97页，编辑于2022年，星期日（三）加载和卸载准则（三）加载和卸载准则n材料进入屈服后，加载状态和卸载状态的本构关系是不材料进入屈服后，加载状态和卸载状态的本构关系是不同的。因此，在弹塑性分析过程中，需要随时了解材料同的。因此，在弹塑性分析过程中，需要随时了解材料屈服后的受载状态，以便确定相应的本构关系。屈服后的受载状态，以便确定相应的本构关系。n在单向受力状态下，可简单地根据应力的增减来判断；而复在单向受力状态下，可简单地根据应力的增减来判断；而复杂应力状态下，则应根据加载函数的增减来判断。杂应力状态下，则应根据加载函数的增减来判断。第59页，共97页，编辑于2022年，星期日理想弹塑性材料的加载与卸载准则理想弹塑性材料的加载与卸载准则n理想弹塑性材料的屈服面位置和形理想弹塑性材料的屈服面位置和形状是不变的。当应力保持在屈服面状是不变的。当应力保持在屈服面上时称为上时称为加载加载，这时塑性变形增大。，这时塑性变形增大。n当应力点从屈服面改变到屈服当应力点从屈服面改变到屈服面内部时，称为面内部时，称为卸载卸载。n nd d卸载卸载d d加载加载屈服面屈服面第60页，共97页，编辑于2022年，星期日应变强化材料应变强化材料d d加载加载d d中性变载中性变载d d卸载卸载n n加载面加载面第61页，共97页，编辑于2022年，星期日（四）流动法则（四）流动法则(增量理论增量理论)n塑性状态下的本构关系与初始弹性范围内本构关系塑性状态下的本构关系与初始弹性范围内本构关系最主要的区别在于应力应变间没有一一对应关系。最主要的区别在于应力应变间没有一一对应关系。n但在某一给定状态下，有一个应力增量，相应地但在某一给定状态下，有一个应力增量，相应地必有唯一的应变增量。必有唯一的应变增量。n因此，在一般塑性变形条件下，只能建立应力应变增因此，在一般塑性变形条件下，只能建立应力应变增量之间的关系。量之间的关系。n这种用增量形式表示的材料本构关系，称为这种用增量形式表示的材料本构关系，称为增量增量理论理论或或流动法则流动法则。第62页，共97页，编辑于2022年，星期日n假设材料在不同应力状态下含有不同的塑性应变能即塑性假设材料在不同应力状态下含有不同的塑性应变能即塑性位势位势Q；在主应力空间中，把同量塑性位势的点连起来得；在主应力空间中，把同量塑性位势的点连起来得到一个塑性势面；到一个塑性势面；n在塑性势面上任一点的塑性应变增量与塑性位势函数的在塑性势面上任一点的塑性应变增量与塑性位势函数的梯度方向一致，即梯度方向一致，即:n上式称为塑性位势理论或流动法则。上式称为塑性位势理论或流动法则。第63页，共97页，编辑于2022年，星期日n其中其中Q为塑性位势函数，是应力分量和硬化参量的函数。为塑性位势函数，是应力分量和硬化参量的函数。n如果如果Q取为加载函数取为加载函数或屈服函数或屈服函数f，称为相适应的流动，称为相适应的流动法则；反之，称为不相适应的流动法则。法则；反之，称为不相适应的流动法则。n目前在岩土工程与工程地质非线性分析中，大多采目前在岩土工程与工程地质非线性分析中，大多采用相适应的流动法则。用相适应的流动法则。第64页，共97页，编辑于2022年，星期日硬化规律硬化规律n是指塑性变形（流动）以后，屈服面的位置发生了变化，如何修是指塑性变形（流动）以后，屈服面的位置发生了变化，如何修正屈服面。对于土壤和岩石材料，采用正屈服面。对于土壤和岩石材料，采用DruckerPrager模式，经模式，经常加一个帽，表示静水压力下的屈服。常加一个帽，表示静水压力下的屈服。（五）一般的弹塑性本构关系（五）一般的弹塑性本构关系第65页，共97页，编辑于2022年，星期日 其中弹性应变增量可由线弹性本构关系确定，塑性应变增量则由流动其中弹性应变增量可由线弹性本构关系确定，塑性应变增量则由流动法则确定。即：法则确定。即：(3-33)（五）一般的弹塑性本构关系（五）一般的弹塑性本构关系第66页，共97页，编辑于2022年，星期日 若材料屈服时，应力分量满足下列一般条件。即：若材料屈服时，应力分量满足下列一般条件。即：将上式写成微分形式得：将上式写成微分形式得：(3-37)第67页，共97页，编辑于2022年，星期日把式（把式（3-33）变换后代人式（）变换后代人式（3-37）：）：(3-41)参数参数A称为硬化函数，是硬化参量称为硬化函数，是硬化参量Ha的函数。的函数。第68页，共97页，编辑于2022年，星期日三、常用弹塑性模型举例三、常用弹塑性模型举例n理想塑性的弹塑性模型。理想塑性的弹塑性模型。n假设：岩土体为具有理想塑性状态的材料，假设：岩土体为具有理想塑性状态的材料，A=0；满足相适应的流动法则，；满足相适应的流动法则，Q=f。n临界状态模型。临界状态模型。第69页，共97页，编辑于2022年，星期日（一）基于理想塑性的弹塑性模型（一）基于理想塑性的弹塑性模型第70页，共97页，编辑于2022年，星期日（二）临界状态模型（二）临界状态模型(帽盖模型帽盖模型)n临界状态模型中最简单的是德鲁克等人临界状态模型中最简单的是德鲁克等人1957年首先提年首先提出的仅反映岩土硬化特性的帽盖模型出的仅反映岩土硬化特性的帽盖模型n固定屈服面是表示材料屈服的极限状态，达到该状固定屈服面是表示材料屈服的极限状态，达到该状态时，材料已临近破坏。它可用莫尔态时，材料已临近破坏。它可用莫尔库仑或德鲁库仑或德鲁克克普拉格屈服函数来表示，写成一般形式，即普拉格屈服函数来表示，写成一般形式，即帽盖屈服面函数的一般形式为帽盖屈服面函数的一般形式为第71页，共97页，编辑于2022年，星期日硬化规律硬化规律n是指塑性变形（流动）以后，屈服面的位置发生了变化，如何修是指塑性变形（流动）以后，屈服面的位置发生了变化，如何修正屈服面。对于土壤和岩石材料，采用正屈服面。对于土壤和岩石材料，采用DruckerPrager模式，经模式，经常加一个帽，表示静水压力下的屈服。常加一个帽，表示静水压力下的屈服。第72页，共97页，编辑于2022年，星期日第三节第三节 非线性问题的有限元解法非线性问题的有限元解法一、迭代法（一、迭代法（讲解讲解）二、增量法二、增量法三、弹塑性问题的有限元解法三、弹塑性问题的有限元解法主要内容：主要内容：第73页，共97页，编辑于2022年，星期日一、迭代法一、迭代法（一）基本概念（一）基本概念（二）牛顿迭代法（二）牛顿迭代法（三）修正牛顿迭代法（三）修正牛顿迭代法第74页，共97页，编辑于2022年，星期日（一）基本概念（一）基本概念n简单迭代法又称逐次迭代法，基本思想是构造不动点方程，以求得近简单迭代法又称逐次迭代法，基本思想是构造不动点方程，以求得近似根。即由方程似根。即由方程 g()=0变换为变换为=(),然后建立迭代格式，然后建立迭代格式，n n当给定初值当给定初值 0 后后,由迭代格式可求得数列由迭代格式可求得数列 n。如果。如果 n收敛于收敛于，则它就是方程的根。因为：，则它就是方程的根。因为：n 第75页，共97页，编辑于2022年，星期日迭代法收敛的几何意义迭代法收敛的几何意义第76页，共97页，编辑于2022年，星期日迭代法发散的几何意义迭代法发散的几何意义第77页，共97页，编辑于2022年，星期日（二）牛顿迭代法（二）牛顿迭代法 用迭代法可逐步精确方程用迭代法可逐步精确方程 根的近似值，但根的近似值，但必须要找到必须要找到 的等价方程的等价方程 ,如果如果 选得不选得不合适合适,不仅影响收敛速度不仅影响收敛速度,而且有可能造成迭代公式发散。而且有可能造成迭代公式发散。能否找到一种迭代方法能否找到一种迭代方法,既结构简单既结构简单,收敛速度快收敛速度快,又不存又不存在发散的问题。这就是本节要介绍的牛顿迭代法在发散的问题。这就是本节要介绍的牛顿迭代法牛顿迭代法的基本思想牛顿迭代法的基本思想 牛顿迭代法一种重要和常用的迭代法牛顿迭代法一种重要和常用的迭代法,它的基本思想它的基本思想是将非线性函数是将非线性函数f(x)逐步线性化逐步线性化,从而将非线性方程从而将非线性方程f(x)=0近似地转化为线性方程求解。近似地转化为线性方程求解。第78页，共97页，编辑于2022年，星期日 对于方程对于方程 ，设其近似根为设其近似根为 ，函数函数f(x)可在可在 附近作泰勒展开附近作泰勒展开 忽略高次项忽略高次项,用其线性部分作为函数用其线性部分作为函数f(x)f(x)的近似，的近似，设设 的根的根 ,则有则有 ,即即 将右端取为将右端取为 ,即即 是比是比 更接近于更接近于 的近似值的近似值 这就是著名的牛顿迭代公式这就是著名的牛顿迭代公式第79页，共97页，编辑于2022年，星期日牛顿迭代法的几何解释牛顿迭代法的几何解释 方程方程f(x)=0)=0的根的根x*是曲线是曲线y=f(x)与与x轴交点的横坐标轴交点的横坐标,设设xk k是根是根x*的某个近似值的某个近似值,过曲线过曲线y=f(x)的横坐标为的横坐标为xk k的点的点Pk=(xk,f(xk)引切线交引切线交x轴于轴于xk+1,并将其作为并将其作为x*新的近似值新的近似值,重复上重复上述过程述过程,可见一次次可见一次次用切线方程来求解用切线方程来求解方程方程f(x)=0)=0的根的根,所所以亦称为牛顿切线法。以亦称为牛顿切线法。第80页，共97页，编辑于2022年，星期日（三）修正牛顿迭代法（三）修正牛顿迭代法 为了避免在每一个迭代步骤中都要重新计算刚度矩阵及其逆矩为了避免在每一个迭代步骤中都要重新计算刚度矩阵及其逆矩阵，可采用修正牛顿迭代法。即在迭代过程中恒将阵，可采用修正牛顿迭代法。即在迭代过程中恒将 取为第一取为第一次循环的次循环的 ，则迭代公式可变为：，则迭代公式可变为：第81页，共97页，编辑于2022年，星期日二、增量法二、增量法 增量法是将作用在受力体系上的荷载划分为若干荷载增量，这增量法是将作用在受力体系上的荷载划分为若干荷载增量，这些增量可以相等，也可以不等，例如分成些增量可以相等，也可以不等，例如分成n级施加，即：级施加，即：第82页，共97页，编辑于2022年，星期日三、弹塑性问题的有限元解法三、弹塑性问题的有限元解法n求解时，采取逐步增加荷载的办法，每次增加的载荷求解时，采取逐步增加荷载的办法，每次增加的载荷尽量小，用求解一系列线性问题来代替求解非线性问题尽量小，用求解一系列线性问题来代替求解非线性问题n求解的方法即增量法求解的方法即增量法n增量切线刚度法增量切线刚度法n增量初应力法增量初应力法n增量初应变法增量初应变法第83页，共97页，编辑于2022年，星期日n弹塑性应变增量和应力增量之间的关系可近似地弹塑性应变增量和应力增量之间的关系可近似地表示为：表示为：nDep 和和E、应力水平有关，而与应力的增量无关，、应力水平有关，而与应力的增量无关，所以（所以（3-68）可以近似看成线性的）可以近似看成线性的(一一)增量变刚度法增量变刚度法第84页，共97页，编辑于2022年，星期日 在起初受载时，物体内部产生的应力和应变还是弹性的，因此可在起初受载时，物体内部产生的应力和应变还是弹性的，因此可以用弹性理论计算。以用弹性理论计算。一旦单元进入屈服，就要采用增量加载方式，此时一旦单元进入屈服，就要采用增量加载方式，此时 位移列阵位移列阵 0 应变应变 0 应力应力 0然后作用载荷增量然后作用载荷增量 R1对于尚在弹性的单元，单元刚度矩阵对于尚在弹性的单元，单元刚度矩阵对于塑性区域的单元，单元刚度矩阵对于塑性区域的单元，单元刚度矩阵第85页，共97页，编辑于2022年，星期日Dep 取决于当时的取决于当时的 刚度矩阵组合成刚度矩阵组合成K0 也和当时的应力水平有关，则有也和当时的应力水平有关，则有可求得可求得第一次载荷增量后的位移、应变、应力可求得第一次载荷增量后的位移、应变、应力的新水平。的新水平。第86页，共97页，编辑于2022年，星期日继续加载，重复上述计算，直至全部载荷加完为止。平衡方程可以写