实习(调研)报告实(调研)报告.doc

实习（调研）报告 一、国内外研究现状及发展趋势 人们对边坡稳定性的认识最早是从滑坡现象开始的，早期对边坡的研究多以长期观测的资料为基础，采用地质历史分析方法进行定性描述。20世纪60年代初，人们认识到边坡稳定性的研究必须将地质分析与力学机制分析结合起来，这期间，主要采用的是刚体极限平衡法[1-4]。进入80年代以来，由于计算技术的发展，使一系列数值分析方法发展起来。近年来一些新理论也被引入到边坡的研究中，它们从不同的角度，对边坡体的结构及稳定性影响因素等进行分析，从而极大地推动了边坡稳定性的研究。 在极限平衡方法理论体系的形成过程中，出现过一系列简化的方法，常见的有瑞典法、Bishop法，Janbu法，Spencer法，Morgenstern - Price法，Sarma法[5]等等。在这以后，又出现了一些改进方法，如改进的Morgenstern - Price法[6]等。另外，Enoki和Yagi （1990）考虑在滑楔间的界面上引入局部强度发挥度，运用滑块离散格式提出了广义极限平衡法。栾茂田（1995）提出了改进的滑楔模型。朱大勇（1997）提出了边坡全局临界滑动场理论[7]。近年来发展起来的基于极限平衡理论的有限元法[8]，是在求解复杂条件下岩土边坡应力和位移分布的基础上，采用各类搜索技术获得临界滑裂面及相应的安全系数，从而评价边坡的稳定性。仍有一些学者致力于极限平衡方法的研究工作，但从目前看还没有完善的已为公认的评价方法，其安全系数的各种表述时常不具有明确的物理意义。因此，鉴于极限平衡方法在计算中所作的种种假定，而且在最危险滑裂面的搜索以及最小安全系数的确定上均存在许多的困难，随着理论体系更为严格的数值分析在岩土工程中的大量应用，一定程度上减小了极限平衡方法的应用范围。 目前有限元在边坡稳定分析中常用的方法有圆弧搜索法和强度折减法等。有限元圆弧搜索法沿用了极限平衡法的思想，只是在假定滑裂面上的下滑力和抗滑力是由有限元计算较为精确地确定。而强度折减法是在有限元计算中，通过不断折减土体强度参数， 当计算不收敛时的折减系数即认为是边坡稳定安全系数。强度折减法一个主要的特征是不需要事先假定滑裂面，从目前应用情况看，是一种比较有发展前途的边坡稳定分析方法。 由于岩土是一种材料特性异常复杂的非连续体，因此也相应产生了一些能处理这种离散非连续体的数值分析方法。如 CundaII（1979）提出了分析颗粒状材料性质的DEM 方法（离散单元法）。近年来得到很大发展的FLAC法（快速拉格朗日法）就是在此基础上发展起来的，它是一种处理非连续、非均匀、各向异性和非线性的新型数值方法。在离散元的基础上，石根华（1988）提出了分析边坡稳定的DDA方法[9]（非连续变形分析方法），非常适合处理非连续的节理化岩体，但由于忽略了岩块的变形过程，因此无法模拟裂隙产生及失稳过程。另外石根华还应用流形中的有限覆盖技术建立了流形方法[10]（MM），解决了在数学上对连续与非连续材料域总体近似函数统一构造的问题，从而使得连续变形与非连续变形问题的统一分析成为可能，能够进行岩体裂隙张开、滑移等特征的分析。 近年来，无单元方法或无网格计算技术发展很快，以 BeIytschko（1994）等基于滑动最小二乘法提出的无单元 GaIerkin法[11]（EFGM）应用最为普遍。此外，田荣、栾茂田等人利用流形方法的有限元覆盖技术，结合无单元法，发展了有限覆盖-无单元法，有效解决了岩土介质渐进变形与破坏、开裂等一系列复杂难题。 随着计算机技术、多媒体和网络技术的发展，基于概率理论的摄动随机有限元法、 基于GIS与GPS的计算机仿真技术也正在边坡稳定评价中发挥越来越重要的作用。 人们在利用极限平衡理论的条分法进行边坡稳定性分析与研究时，首先使用枚举法搜索临界圆弧滑面。随着计算技术的发展，临界滑裂面的搜索开始转向数值分析方法， 也就是通常所说的最优化方法，这一方法主要又可分为以下： 模式搜索法是按照一定的模式，比较不同自变量间的目标函数，确定最优的搜索方向，最终找到最小值。Lefebvre（1971）[12]，Huang（1983）[13]都采用了模式搜索方法来确定临界滑动面，后来莫海鸿（1999）[14]在此基础上提出了改进的模式搜索法，还有江见鲸等[15] 基于简化 Bishop方法，采用消去法中的二分法求算临界滑面。另外， Giam 和Don-aId（1988）[16] 提出了一种由已知的有限元应力场确定临界滑移面的模式搜索法，它是根据有限元计算得到的应力场，由应力水平较高点出发，形成一个由坡顶到坡底的破坏面。这类方法一般是用来确定圆弧滑动面，而且该法的搜索点是在搜索过程中产生，容易陷于安全系数的局部极小值。 牛顿体系的分析方法借鉴了变分法中把安全系数看成是临界滑动面