边坡稳定性的弹塑性有限差分简介.docVIP

边坡稳定性的弹塑性有限差分简介 　　摘要：有限差分法可以在极少的假定条件下能够得到边坡土体的应力-应变性质，预测边坡结构的稳定或失稳趋势，还可以动态地模拟边坡在施工条件下的稳定性，同时还能够得到边坡土体介质的应力、应变、位移、塑性区和安全系数等有用结果，克服了极限平衡法和工程地质类比法的缺陷。本文采用降雨渗流条件下有限差分程序，通过分析塑性区和剪应变率云图自动获取最危险滑动面，最后通过强度折减法获得边坡稳定性的安全系数值。 　　关键词：边坡稳定性；弹塑性；有限差； 　　中图分类号：TU74 文献标识码：A 文章编号： 　　 　　1概述 　　迄今为止，在岩土边坡稳定性的理论研究与工程应用中，广大科技工作者和工程技术人员几乎都采用刚体极限平衡法和工程地质类比法进行边坡稳定性及防护与加固等方面的研究。这两种方法都基于太多的假定条件，在计算过程中没有考虑材料本身的形变特征，甚至也不能得到边坡岩土介质的应力、应变、位移和塑性区等关键量的分布，因此无法定量、准确地描述边坡的破坏与稳定状态，并且在求解边坡安全系数时要通常事先假定滑裂面的形式，如圆弧形，对数螺旋性，折线形等形式，再在该滑动模式下进行安全系数计算，这样计算出来的结果用于稳定性评价带有很大的主观性。然而，岩土弹塑性有限差分数值计算方法可能够虑边坡岩土介质的非均质、非连续、大变形等特点，可以很方便、准确地得到边坡岩土介质的应力、应变、位移和塑性区等关键量的分布特征，使得计算结果更加精确有效，与工程实际较吻合。本章基于塑性力学基础理论，运用有限差分强度折减法分析边坡安全系数的数值计算概念，并采用有限差分强度折减法结合边界条件和初始条件，计算边坡的安全系数，并可以得到令人满意的结果，从而指导工程实践。 　　大量实验和实例表明：岩土材料的弹塑性变形主要是由可恢复的弹性变形和不可恢复的塑性变形组成。一般用虎克定律（或广义虎克定律）来计算弹性变形部分，在求解过程中通常认为：在承载能力极限条件下，当岩土材料的应力超过比例极限后或在正常使用条件下当岩土材料等效应变超过规范允许变形值时，应力应变关系不再遵从弹性虎克定律所规定的线形关系，即应力应变为弹塑性变形，这时进行岩土结构计算，关于静力平衡方程（1）总刚度矩阵（位移向量的函数）是不能直接求解的，必须要借用非线性方程组的数值求解方法才能对该式进行解答。 　　（1） 　　现在简要地介绍下土体的塑性力学主要组成部分，即弹塑性应力应变关系，屈服准则和破坏准则，流动法则和加工硬化规律，并着重阐述在岩土边坡数值分析中如何应用塑性力学理论，建立岩土和结构材料的弹塑性应力应变关系，即本构关系。 　　2岩土塑性力学基础 　　弹性力学、塑性力学和材料力学都是连续介质力学的分支，广泛应用于岩土、混凝土等材料。岩土弹性力学主要研究岩土在弹性工作阶段的应力应变规律，岩土塑性力学则研究岩土在塑性工作阶段应力应变的规律。从实验得知：岩土材料在弹性阶段，其应力应变关系满足线性关系，可以用虎克定律（或广义虎克定律）求解。在塑性变形阶段，其应力应变关系满足非线性关系。岩土材料在弹性工作阶段时，卸除荷载后其变形可以全部恢复；然而，进入塑性变形后，其变形在卸载后是不能完全恢复的，其中不能恢复的残余变形有部分为塑性变形。另外，塑性阶段和弹性阶段的差别还在于加载和卸载性质的不同，而且在塑性变形阶段材料的应力应变关系还依赖于应力历史和应力路径。 　　根据岩土材料单轴抗压强度试验可以得到应力-应变曲线，其曲线形式不同（见图图1），基本上可以分为以下三大类： 　　（1）理想弹塑性特征； 　　（2）应变硬化特征； 　　（3）应变硬化-软化特征。 　　岩土材料可定为第（2）、（3）类。 　　2.1.屈服准则和破坏准则 　　在单向受压应力状态下，当应力值达到岩土材料屈服极限时，该类材料开始发生屈服变形，即产生塑性流动。因此，受压应力=，这是单向应力状态下的屈服准则，又称屈服条件(yield condition)。然而，在复杂应力条件下，岩土材料中某一点开始产生塑性变形，这时需满足的条件通常为屈服函数： 　　(2a) 　　式中，实验常数。 　　当材料处于各向同性状态时，其坐标方向与屈服无关，故有 　　(2b) 　　或 　　 (2c) 　　经典塑性力学通常认为在静水压力作用下不影响材料的屈服，此时屈服条件可表示为应力偏量不变量的函数 　　(2d) 　　岩土材料的屈服准则通常表为 　　 (2e) 　　 (2f) 　　屈服条件可以表示为一标量函数(scalar function) ，该函数在主应力空间的图形称为屈服面(fieldsurface)。从图中可以看出屈服面是由多个屈服应力点拟合起来所构成的一个空间曲面，在该空间曲面中，屈服面所包围的空间