弹塑性理论基本知识[共32页]

1、岩体数值计算原理与方法,张强勇,目 录 1.绪论及张量基本知识 2.弹塑性理论基本知识 3.有限元计算方法 4.节理岩体弹塑性损伤断裂有限元计算方法 5.节理岩体损伤锚固有限元计算方法 6.大型工程弹塑性损伤及加锚数值计算分析 7.大型三维弹塑性和弹塑性损伤有限元计算程序使用介绍,2.弹塑性理论基本知识 2.1 弹性力学基本知识 2.2 岩土塑性力学基本知识,2.1 弹性力学基本理论知识,1.应力方向规定,正面上的应力，以与坐标轴的正向一致为正，反之为负。 负面上的应力以与坐标轴的负向一致为正，反之为负。 中 j表示应力方向，i表示应力所在截面的外法线方向,2. 应力变换公式 取任意一个斜面，

2、设其外法线上的单位矢量为ni，它的三个分量是方向余弦，即： 且有,设该斜面上的应力矢量为 ，斜面的面积为S，S在X1OX2内投影为S3，在X1OX3内投影为S2，在X2OX3内投影为S1。由四面体的平衡条件 得到: 已知： 则,截面上的应力矢量的各分量,将Pj投影到任意新的m坐标轴上去，可得： 上式表明： 已知材料中任意一点的应力状态 ，则通过该点的任意截面 （单位外法线矢量为ni）上的应力矢量Pj在任意方向 上的投影 为： 由于满足 形式，表明应力是二阶张量,已知一点的应力，利用应力变换公式可求得通过该点的任意截面上的正应力和剪应力。 如下图，已知一点的平面应力状态 可求得任意斜截面上的正应

3、力 ，剪应力,利用应力转换公式 ，可求得,3.主应力和主应力张量不变量 主应力：主平面（剪应力为零的平面）上的正应力。 设任一主平面的外法线方向为ni（i=1，2，3） 由 因n1，n2，n3不全为零，则,由 可求得主平面上的主应力矢量在三个坐标轴上的投影 则,因为主应力和坐标系的选择无关（即用主平面上的主应力描述一点的应力状态不随坐标系而变化），因此 在坐标变换时也保持不变，故称 分别为应力张量的第一、第二、第三不变量,4.偏应力张量及其不变量 由于任何张量 总可以分解为球张量和偏张量两部分，即 球张量 ：不随坐标系而变化的张量。 偏张量 ：张量与球张量之差。 应力张量 也可分解为应力偏张量

4、 和应力球张量 在主应力状态下,与应力张量不变量类似可得到偏应力张量的三个不变量 : 式中: 为偏应力第一不变量； 为偏应力第二不变量； 为偏应力第三不变量,5. 应力空间 一般应力空间：应力张量有6个独立分量，可视为6维空间中的一个点，这6维空间中的每一个坐标代表一种应力分量（或说代表一点的应力状态），这个应力空间称为一般应力空间。 应力路径：研究固体中一个点的应力变化可用应力空间中的一个点的轨迹来形象加以描述，这个轨迹表示了应力变化的历史，称为应力路径。 主应力空间：确定一点的主应力需要6个独立变量，三个主应力值和三个决定主轴方向的参数，如果研究的材料是各向同性的，主应力的方向可以不考虑，

5、则一点的应力状态只由三个主应力值决定，如果用三个互相垂直的轴各表示一个主应力值，则得到一个三维的应力空间，称为主应力空间,主应力空间中任意一点的主应力矢量为 , 把 分解为球张量和偏张量两部份： 由 知表征应力球张量的单位矢量n 与三个主应力坐标轴的夹角相等， ， 我们称与三个主应力坐标轴成相同夹角 的直线为静水压力轴，在静水压力轴上每一点的 球张量的几何解释是静水压力轴上的一个分矢量,6.八面体应力和等效应力 八面体： 如果一个平面的法线和三个坐标轴的夹角都相等（即 ），则这 个平面称为等斜面，在主应力空间中，这样的平面共有八个，它们组成的几何体称为八面体。 八面体上的正应力 与第一应力不变

6、量I1有关，八面体上的剪应力 和第二偏应力不变量J2有关,等效应力： 许多文献为了把三维应力状态和单轴试验中的拉（压）应力联系起来，定义了等效应力 在单向拉伸时 ，则 为了把三维应力和纯剪试验中的剪应力联系起来，定义了等效剪应力,2.2 岩土塑性力学基本知识,1. 塑性材料分类 理想弹塑性材料 应变硬化材料 应变软化材料 塑性材料的特征： 加载、卸载遵循不同的应力应变路径，是不可逆过程。 应力与应变之间不是一一对应的关系，应力应变关系与加载历史有关，即： 式中 为内变量，反映加载历史和材料微观结构变化的量，其不能直接测量，因此称为内变量,2. 理想弹塑性材料 屈服：材料应力达到某一定值的时候变

7、形可以不断增加而应力不再提高，这种现象称为屈服。 屈服面：在应力空间中，每一种应力状态（ ）下都能找到一个确定屈服的临界 ，这些值所形成封闭曲面称为屈服面 屈服准则：反映材料产生屈服的条件 称为屈服条件（屈服准则,对于金属材料认为： 其屈服条件只与主应力 大小有关，与其方向无关。 其屈服条件只和应力偏张量有关，而与 平均应力 无关。 这就是说材料的屈服面在空间中与静水压力轴平行,3.Tresca准则和Mises屈服准则 Tresca屈服准则 : 认为材料的最大剪应力达到一定极限值k，材料就会产生屈服，其屈服条件只与 、 有关，和中间主应力 无关。 Tresca屈服面在应力空间为一六棱柱面，在平

8、面上的投影为一正六边形,Mises屈服准则 : 认为偏应力张量的第二不变量达到某一极限值C，材料产生屈服。 J2=C J2在平面上的投影为一个圆，因此Mises屈服 面为一平行静水压力轴的圆柱面,4.流动法则和加卸载准则 流动法则 : 式中G为塑性势函数 ；为比例常数，q为内变量（反映加载历史和塑性变形的变量） 加、卸载准则,5.塑性一致性条件 塑性一致性条件： 屈服面 在应力空间随塑性内变量q的变化 而变化，但材料在加载过程中出现的新的应力点（ ）将始终在后继加载面（屈服面） 上，即， 这种新的应力始终保持 在加载面上的条件称为塑性一致性条件。即,由于塑性变形的影响，材料加卸载弹性模量随塑性

9、变形的变化而变化，这种特性称为弹塑性耦合,6. 弹塑性耦合,DruckerPrager屈服准则 : D-P准则考虑了静水压力（围压效应）的影响，随静水压力增加的增加，屈服面半径不断扩大 D-P准则在应力空间中为圆锥面，是Mises准则的推广,7.岩土材料中的屈服准则,Mohr-Coulomb（M-C）屈服准则,为材料内摩擦角,为Lode角,M-C准则在应力空间为一六棱锥面，是Tresca准则的推广,注意：D-P准则和M-C准则只适用岩土受压情形，即，在受拉区和实验有较大出入,Zienkiewicz-pande（辛潘屈服准则,D-P准则适用于软弱夹层 Z-P准则适于模拟岩石、砼材料,8. 岩土材料的弹塑性本构模型 弹塑性力学响应特征: 式中：为内变量（反映加载历史和塑性变形的量）； 为比例常数；H为反映材料软、硬化特