第九章岩石力学的数值模拟导论.ppt

（4）进一步提出可采用形象法来寻找参考线、参考面，并给出了具体的寻找方法和步骤。利用形象法亦可迅速完成参考面的识别，该方法形象、直观，这是一种高效率的搜寻方法。 （5）建立了块体静动力学约束条件――块体间无侵入发生及不承受拉伸力的控制标准和位置锁定原则。 （6）实际工程中非连续面的大小是有限平面，它们在空间分布上具有很大的随机性。针对这一情况，作者分析了目前处理节理网络技术方面的不足，提出利用拓扑学单纯同调原理解决块体的形成问题。根据空间各有限平面，可以为有限元法、离散元法 块体理论以及其他要求具有一般形状的数值方法提供模型数据。该方法不仅可以生成规则的块体，而且还可以生成形状十分复杂的多边形、多面体，甚至是带孔洞的块体。 （7）根据总刚矩阵 的特点，提出了一种高效三角分解方法――图解法。这种借助于图解理论的非零元素存贮方法使得计算过程直观而高效。通过图形和连接矩阵的配合使用，使整个过程直观明了，具有存贮量小、计算量少、舍入误差小等优点。利用图解技术，事先对块体优化编号，确定存贮量，并使三角分解过程中产生的新元素个数最少，从而达到减少数据存贮和计算，提高解题效率的目的。 （8）根据非连续变形分析理论，编制和开发了相应的程序。并结合三峡船闸高边坡稳定性分析实例进行了研究和分析，其计算结果与实际观测结果大体接近。 （9）从理论上解决了非连续变形分析理论与有限元、边界元及传统极限平衡方法的耦合，并给出耦合原理及耦合的公式。 （10）针对非连续变形分析法存在的问题，提出解决方案。时步的大小事关非连续变形分析理论计算结果的正确性，通过众多的计算表明仅当时步小于0.005s时，计算结果是稳定可靠的。 （11）针对非连续变形分析理论自身的特点及高边坡滑坡位移复杂性，从理论上深刻分析、研究了边坡滑坡时产生混沌的机理及过程，并提出防止边坡进入混沌状态的实际措施。 （12）将不连续变形分析与AutoCAD动态技术结合，使块体系统非连续变形分析法的计算实施过程形象、直观，完成了块体系统的变形或破坏全过程的模拟。 2、不连续变形及其特点 非连续变形分析以各个块体的位移作为求知量，通过块体间的接触和几何约束形成一个块体系统，块体单元受非连续面的控制，在块体运动的过程中单元之间可以接触，也可以分离。在块体运动过程中，满足块体间不侵入和不承受拉伸力的条件，总体平衡方程是由系统的最小势能原理求得的。求解方程组即可得到块体当前时步的位移场、应力场、应变场及块体间的作用力，反复形成和求解总体平衡方程式，即可得到多外时步后块体的变形，应力及位移情况，也可求得块体系统最终达到平衡时应力场、位移场等情况以及运动过程中各块体的相对位置及接触关系。 当然，非连续变形分析方法也和其它任何一种数值方法一样，有其一定的适用范围和局限性，这主要表现在： （1）岩体的种类繁多，其性质极为复杂。在通常情况下，非连续岩体的力学表现为非均质、各向异性、非线性等复杂的特性，随着应用工程环境的不同，使岩体数值计算有许多非确定因素。这对非连续变形分析在对岩体参数（ ， ）取值时带来麻烦，岩体参数的取值直接影响到计算结果的正确性。目前，非连续变形分析的数值计算中，一般假定岩体是弹性体，对塑性、粘性体等尚不适用。 （2）计算时步的大小肯定会影响非连续变形分析的最终结果，时步过大将导致计算结果与实际测量值相差较大，从更深理论上研究时步的合理取值对非连续变形分析的计算具有较大的指导意义。计算时步的大小与实际问题联系起来，合理地考虑时间因素也是今后需要研究的问题。 （3）非连续变形分析方法比较适合模拟硬岩，对于软岩或软硬相间的情况在判定块体接触时会遇到困难。也就是说非连续变形分析比较适合于块体自身小变形的情况。 （4）在非连续变形分析实施的过程中，有大量的数据需要处理，需要大量的计算机时间来存贮、转换和调用，耗用了大量的计算机CPU时间，应优化和改进计算方法，缩短计算时间。 3、DDA今后的发展方向 DDA是近10年来迅速发展起来的一种较新的岩土数值计算方法，已于1996年在美国召开过一次国际性研讨大会。之后，在日本、英国相继举行了研讨会。若能将非连续变形分析和已发展成熟的有限元、边界元以及极限平衡等方法相结合，发挥各自的优势，则会使这一方法更趋于完善。对非连续变形分析方法的研究，应重点改进前面提及的不足之处。除上述工作外，正在进行和将要进行的工作还应包括如下的内容： （1）继续完善三维非连续变形分析理论以及计算程序的实现仍然是今后一个时期内需要解决的问题，优化计算机程序和设计方法是提高解题效率和推广非连续变形分析法应用的前提。 （2）从理论和计算机程序设计上解决颗粒形状单元在不同颗粒级配条件下的非连续变形分析。对常见的几种概率分布形式下的颗粒