非连续数值方法综述

非连续数值方法综述

杨凡

（河海大学水利水电学院，江苏南京210098）

摘要：非连续问题是岩土及水利工程中不可避免的一类难题，由于其对工程的影响巨大，

近几百年来特别近一个世纪以来一直是工程界研究的一个热门话题。从最早的非连续问题

解析解法—刚体极限平衡法出发，引申出近几十年来有关非连续问题研究的热点—非连续

问题的数值解法，然后对这些非连续的数值方法的基本原理和实际应用发展情况进行一一

综述。

关键词：非连续；数值方法；

岩石和土都是经历过变形的地质体，受其成因、组成、结构、年代等诸多因素的影响，岩土材料具有高度的非连续性、非均匀性和各向异性的特征，在力学性质上表现出强烈的非线性。岩土工程是一门综合应用岩石力学、土力学、工程地质学等基本知识解决实际工程中有关岩体与土体变形及稳定问题的学科[1]。岩土工程中的非连续变形问题主要是由岩石及土体中不连续面的存在引起的，岩土工程问题中的不连续面大致可分为两类，一类是指存在于岩体中的节理、软弱夹层以及土体中的剪切破坏面，另一类则是岩土结构如各类基础、挡土结构、地下结构等与岩土体之间的接触面。显然，不连续面对岩土体或结构的受力、变形有着重要的影响，因此为使计算结果真实地反映出岩土体及结构的受力和变形情况，在计算时不能忽视不连续面的存在[2]。

对于具有不连续面的结构，在承受荷载的过程中，不连续面的状态是在不断变化的，这将影响到两侧岩土体的应力和变形，从而影响到整个体系的应力场，而应力场的改变又影响到不连续面的状态。因此，解决岩土力学问题的关键在于对非连续变形的模拟，分析研究结构中各种不连续面的构造特点和力学性能，研究其受力状态的变化规律及其对结构整体性能的影响是工程设计中的关键研究课题之一，具有很大的学术意义和实用价值[3]。

几百年来，人们对非连续变形问题作了大量的研究工作。最早有关非连续问题的研究主要集中在寻求解析解的层面上。1773年，法国科学家库伦在大量实验基础上总结了著名的库伦土压理论，刚性楔体和静力平衡的应用也为后续研究奠定了一个基调。在此基础上，瑞典圆弧法的提出表明刚体极限平衡方法的正式诞生，在接下来的几十年间该方法得到了很好的发展并出现了Bishop、Janbu、Sarma等改进的方法。在发展过程中，极限平衡法在水利工程上得到了很好的运用。但其缺点亦非常明显，即不考虑土体的变形，得到的结果偏安全。在许多的工程中，其运用受到局限。

近年来，伴随着计算机技术的高速发展，有关非连续问题的数值分析逐渐发展成为热门。非连续分析的数值方法主要包括：界面单元有限元法(FEM with interface elements)、刚性有限元法（RFEM）、离散单元法(DEM：Distinct Element Method)、非连续变形分析(DDA：Discontinuous Deformation Analysis)、数值流形方法（NMM：Numerical Manifold Method）、无单元法（(Element—Free Method）、耦合方法（Couple Methods）以及渐进破坏模型(Progressive Fracture Models)等[4]。以下对各个方法研究进展情况进行综述。

1 界面单元有限元法(FEM with interface elements)

有限单元法在连续介质上面的成功应用，使其广泛地被人们接受，但在解决前处理问题、应力与应变解答不连续问题计算等方面还存在着一些局限。为使其能够处理简单的非连续性问题，人们提出了各种能够反映非连续性质的简单力学元件或特殊界面单元用于模拟单一非连续界面的力学行为[4]。这类方法主要包括有：联结单元(Linkage Element)（Ngo和Scordelis，1967)，无厚度接触单元法或节理单元法(Joint Element)(Goodman，1968)、薄层单元法(Thin-Layer Element)(Desai，1982)、接触—摩擦单元(Contact—Friction Interface Element)(雷晓燕和swoboda，1994)，其中联结单元法与节理单元法是在性质差异悬殊的两相邻单元间界面上设置法向与切向弹黄，一旦法向弹黄受拉则视界面发生分离，若切向相互作用力超过界面的极限强度则认为界面产生相对滑移。而薄层单元法假定界面是由特殊介质所组成的厚度很薄的实体单元，这种单元的本构性质取决于周边介质力学性态，由试验直接测取。通常认为界面的法向性质与周围岩土材料相应特性相同，Desai与Nagaraj(1988)对此曾提出采用非线性方程来考虑法向应力—应变关系。关于切向响应特性，Desai等研制了一套循环直剪试验仪来测定界面剪力与相对滑移距离之间的依赖关系，并建议用多项式函数或五参数Ramberg —Osgood模型加以描述.所谓接触—摩擦单元与接触单元的区别在于它采用节点接触应力代替节点接触力，对于接触面的几何方面与静力学方面约束条件用附加条件纳入到系统方程中，从而能够通过计算自动满足接触界面的滑移、分离以及粘结真实状态。

这类方法将岩土介质视为准连续介质，仍以连续分析为主，但可以对个别具有控制作用的宏观非连续面的变形与破坏等力学效应给予重点分析，在工程中得到了一定应用，但由于下述自身的缺陷而受到了限制，在处理复杂的非连续变形问题时显得无能为力[4]。例如，(l)只能对原生的非连续界面进行计算，对次生的非连续界面无法处理；(2)界面单元数目不能设置得太多；(3)界面弹簧刚度的选取较为困难。

2 刚性有限元法（RFEM）

刚性有限元法最早由日本Kwaai教授于1976年提出，该方法将离散后的块体视为刚体，块体之间用界面上的法向弹簧和切向弹簧相联结，以块体形心处的刚性位移为基本变量，用分片的刚体位移模式逼近实际整体位移场，以块体间的联结弹簧反映结构内部的弹性，并用界面应力表征结构内部的应力[5]。刚性有限元中的模型主要包括刚体—弹簧模型(RBSM)、刚性有限元、分块刚体位移—界面应力元和块体—弹簧模型等。此类模型的界面特性均服从Coulomb摩擦定律，对体系的静力学约束条件考虑得比较充分，在连续状态的应力分析方面可给出较高的计算精度，对于临界状态能够估算出极限荷载，并可有效的用于少量块体界面间的摩擦接触分析。但这类方法过分强调岩土体结构面的作用，对结构体的变形没有给予足够重视。虽然这种方法可以用于原生界面的破坏分析，但不能模拟实际岩土体的破坏发展过程，也无法模拟破坏发展导致的次生界面的非连续变形行为和块体失稳后的运动过程。

3 离散单元法(DEM：Distinct Element Method)

离散单元法(DEM：Distinct Element Method)是由美国Cundall等(1971)提出，它特别适用于含