岩体非连续变形数值分析方法的特点及比较.pdfVIP

2007年第4期 西部探矿工程 27 岩体非连续变形数值分析方法的特点及比较 张绪涛，陆晓敏 (河海大学土木工程学院，江苏南京210098) 摘要：岩体稳定性分析，有几种常用的非连续介质分析方法。概述了几种常用的数值分析方法，对 这几种分析方法作了比较，得出它们的区别与联系，并总结了理论存在的不足。 关键词：离散单元；DDA；拉格朗日元；块体单元 中图分类号：TU45文献标识码：B文章编号：1004—5716(2007)04—0027一03 岩体稳定性分析一直是岩土工程中一个重要问题。 统。在块体运动的过程中，严格满足块体阎不侵入。将 在工程中，岩体稳定性分析的结果正确与否直接关系着 边界条件和接触条件等一同施加到总体平衡方程。总 该工程的成败。目前岩体稳定性分析方法有很多种，特 体平衡方程由最小势能原理求得。求解方程组即可得 别是近年来，有学者提出基于不连续介质的力学分析模 到块体当前时步的位移场、应力场、应变场及块体间的 型，如离散单元、DDA、拉格朗日元、块体单元等[1]。主作用力‘6|。 要阐述这几种方法的区别与联系。 1 非连续变形数值分析方法 1．1 离散单元法 方法分析了重力坝抗滑稳定，并与有限元方法进行比 离散单元法是由Cundall于20世纪70年代提出较，文献[9]将该方法用于地表变形分析。 的。离散单元法是针对节理岩体提出的一种适用于模 1．3拉格朗日元法 ofContin— 拟岩土体大位移的数值方法。它是在块体准刚性假设 拉格朗日元法(FastI．agragianAnalysis 的前提下，以牛顿第二定律为理论基础建立起来的，以 微小运动状态的求解模拟岩体的大位移。其基本思想 集团于1986年提出的。该方法起初应用于流体力学， 是将所研究的物体看作被节理在内部分割的离散状块 研究每个流体质点随时问而变化的情况，即着眼于某一 体单元，块体间通过边界接触力相互联系，边界接触力 流体质点，研究它在任意一段时间内轨迹、所具有的速 用物理方程来描述。当所研究的块体在力系作用下或 度、压力等。后来该法被移植到固体力学中，将所研究 边界条件发生变化时，单元就在牛顿第二运动定律下发 的区域划分成网格，网格的节点相当于流体的质点，然 生平移或转动，允许调整各个块体单元的接触关系，接 后按时步研究网格节点的运动。这种方法适合于求解 触关系的调整产生新的接触力，在该力的作用下产生新 非线性的大变形问题，其按时步采用动力松弛的方法求 的的位移。如此循环：在不平衡力作用下，块体产生运 解，不需要形成刚度矩阵，不用求解大型联立方程，可用 动。不断调整位移、调整接触力，最终块体可能达到新 微机解决较大的工程问题[1州。 的平衡状态，也可能一直运动下去[1]。 拉格朗日元法对于岩土工程和采矿工程是一种设 离散单元法适用于节理岩体的应力分析，文献[2— 计和分析的强有力的方法。文献[113将该方法用于边 3]基于该方法模拟了边坡稳定性分析，并得出与工程地 坡稳定性分析；文献[12]是基于该方法对基坑进行优化 质特征相吻合的结论；文献[4]介绍了离散单元法在钢 设计。同样在文献[6]中也详细介绍了拉格朗日法的应 筋混凝土构件非线性分析中的应用；文献[5]详细介绍 用。 了离散单元法在工程领域中的应用。 1．4 刚体弹簧元的改进——块体单元法 1．2 DDA方法 块体单元法由任青文等改进的。该方法是在刚体 D