三维离散元法中的接触问题.pdfVIP

第六届全国岩士力学数值分析与解析方法讨论会，中国，广州 1998年8月4～6日 三维离散元法中的接触问题 焦玉勇 葛修润 (中崮科学院武汉岩土力学研究所 武汉 d301／71) 提要块体问的接触翔断是兰维离散元{去中最为关键柏环节．也屉制约计算速度和精度的最主要的因幂。本文详细探讨 r三摧离教元=i击中的接触问题，对榴部块俸的检索、接触力及其作用点的计算等做丁某些改进，在此基础上编制r三维离 敬元接程序，并对算倒进行r试葬。 荚量词三维离散元；接魅问题；枢邻块体；接触力 ￡i，i 1 引 言 离散元法是美国学者Cundall于j 971年提出的一种用于模拟节理岩体力学行为的数值方法。出于该 方法体现了岩体变形的非镯性和非连续性等物理本质，克服了常见的连续性方法无法模拟大位移、大变形 的缺陷，因而在近30年来得到了广泛的应用。离散元法最初是基于动态松弛的，即用动态方法统．求解静 态、准静态和动态问题，一开始被用来分析岩质边坡的稳定性问题。迄今为止，甲面问题的离散元法0经相 当成熟，其中商品化软件UDEC是岩土力学领域公认的经典离散元法程序，』￡成功的应用范例不胜枚举。 与此相对照，空闻闷题的离散元法却发展较为缓慢，一方面是众所周知的计算机水平的限制．另一方面是 团开发出了商品亿程序3DEC，该程序功能强大，在国外已有几个应用的例子，仍闪其售价太高，嘲内至今 还没有引进。我们认为，作为一种有实用价值的数值方法，不管其理沧深度她何．赴国内开展该方面的研究 还是有意义的。在这里，本文不准备涉及具体原理，只是对接触判断这一i维离散元次中的关键问题加肌 探讨。 2相邻块体的搜索 l 3 4 5 在一个块体系统中检查一个块体同哪些块体相邻 是不困难的，但要使检索次数正比于块体数Ⅳ丽不是 6 j／／＼s 9 i0 Ⅳ2却不是一件容易的事，因此必须采用适当的形体几 ／ 何描述表达块体和研究区域，建立行之有效的搜索葬 祛，才能在大范围几何空间内成功进行相邻块体的检 11 ＼12 。3 14 15 验。就有限元丽言，单元的几何形状由大量的结点坐标 ＼ 决定，单元的表面及单元间的空间关系由～个结点的 √ 有序排列矩阵——单元联结肄碲定，以贽单元间组合 J6 】7 18 19 20 起来决定刚度矩阵。在离散元法中，单元的数据存贮不 2l 22 23 24 25 仅能体现单元的几何位置，而且应该便于迅速的查找 和修改，因此采用动态数据链表对块体进行几何描述 是比较适宜的，也便于实现块体的“物理地址”和“内存 140 兰维离散元法中的接触问题 焦玉勇葛修润 地址”的统一。 有了块体的几何数据，将它们在所研究的区域中定位并由此进行检索就比较容易了，在这里我铘采用 一种称为“基于实体搜索”的方法，具体做法是：张比模拟区域稍大的范围内框出一个立方体，设其边长为 f，以此立方体为边界将整个空问划分为若干个小长方体(设有m个)。为直观起见我们用平面图形说明， 25)，空间有m个“盒子”组成，称为“盒子空闻”。根据块体的结点坐 如图1所示，m一25(三维空间中，n一1 标，计算块体的最大轮廓，为长方体形．轮廓面平行于各坐标面。将每个块体的轮廓面乎行的向外扩展给定 值CTRL，