3D+DEM在节理岩体研究中应用初探.pdfVIP

中国CSRM鞋岩工程专业委员会第二届学术大会论文集 3D DEM在节理岩体研究中的应用初探 焦玉勇 刘泉声 郭启军许锡昌 (中国科学院武汉岩土力学研究所武汉430071) 摘要根据3DDEM的基本原理，蟪制7相应的c语言程序，并对‘程算 例进行了模拟，计鼻结果表明；本文的程序是有效的，3DDEM是模拟节理岩体力 学行为的行之有效的敷值方法。 关麓调3DDEM，C程序，节理岩体 1引言 离散单元法(DEM)是1971年美国学者Cundal[提出的一种模拟岩土体非连续大变形 的数值方法，由于抓住了裂隙岩体变形的非连续性这一本质，固而它一出现就为国际岩土工 程界所接受，并在近30年中发展成为一个重要的岩土工程数值计算分支领域。目前，二维离 散单元法的应用已经比较成熟，二维程序UDEC是岩土力学领域公认的行之有效的非连续 问题分析程序。三维离散单元法(3DDEM)的理论在80年代末期才趋于完善，商品化程序 3DEC由美国ITASCA咨询集团开发成功。由于多年来3DEC一直保持着较高的价位，其应 用实例并不多见。我国在三维离散单元法的研究才刚刚起步，整体上落后于国外10年的时 间。鉴于我国的国情和三维离散单元法的巨大潜力，在国内开展这一方面的研究是根有意义 的。 ， 2 3DDEM的基本原理 离散单元法以天然节理、裂隙切割岩体所形成的块体为单元：计算每个单元上的各种作 用力，如重力、相邻块体的接触力、地下水压力、锚杆作用 力等，求出合力和合力矩，列出动力学方程式，通过差分 格式解出～个微小时段的速度和位移；在空间域上对所 有的块体实施上述求解过程，并对时域积分，岩体工程的 大位移、大变形就可求得了。这一求解格式避开了方程组 求解过程中对矩阵的存贮和处理，计算效率大大提高。除 位移求算外，三维离散单元法中较为复杂的是块体间的 接触判断，详情请参考文献[1，6]，下面仅给出基本原理 田1相接触的两块体 的矢量形式推导。， 2．1块体问的相对位移 一 如图1所示，设块体i，J在c点接触，它们的形心位置矢量分别为： 艿， X“1一{z“’，Y“1，z“’)，x“’={z“’，Y‘”，z“’)，形心的平动位移和转动位移分别为：甜” ={“。’，v“’，甜“1)，瑚唧一{越‘”，△以‘”，触‘‘1)，△x‘”={“叫，口∽，w“’}，出‘，1一{触‘”，肋，c”， △以“’)。若接触点的位置矢量为jrc={∞，，z}，则接触点C按照块体i和j的形心求得的位 移分别为： △jr91=△Y“’+△日‘‘’×(zc一工∽)] 世扩=埘“’+瑚‘p×(k—jr∽)f ‘1’ 在c点块体i相对于块体』的位移为： 蚪c=ⅡP’一衄y’ (2) 将式(1)代入式(2)并整理得； △k=[B]{U) (3) 其中， ． ■oo o H“’一(，一y“’)一10 0 0 一(z--z“’)y--y“’] 01o—k—z“’) o X--X。’ o一1 o =一￡∞ o [B]=1 一杠一。“’)l LO 0 l 0 J y--y‘‘’一(x-x‘。’)0 0—