基于平衡颗粒法生成地质体复杂节理面网格方法研究.pdfVIP

2012 2012 9 16 -18 颗粒材料计算力学会议 年 月 日 日 National Conference on Computational Mechanics of Granular Materials (CMGM-2012) 湖南 · 张家界 论文编号：CMGM2012-D06 ISBN ：978-7-5611-7240-7 基于平衡颗粒法生成地质体复杂节理面 * 网格的方法研究 段文杰1† 2 2 李世海 冯春 1 , 300401 河北工业大学 天津 2 , 100190 中国科学院力学研究所 北京 摘要：宏观地质体中错综复杂的节理拓扑关系给节理面网格的划分带来了困难，进而制约了对含复杂节理地 质体的力学效应研究。文章提出了一种利用颗粒生成节理面网格节点的方法——平衡颗粒法。该方法首先将 复杂节理系统按照节理交线进行分割，而后通过节理面内颗粒的随机排布产生密实、微嵌入的颗粒系统，进 而采用全量法计算颗粒间的接触力，同时引入指向节理面中心的牵引势函数，在粘性阻尼的控制下，通过显 式迭代，逐渐实现接触力产生的斥力与势函数产生的引力之间的平衡。算例表明，基于平衡颗粒法产生的节 理网格节点排布均匀，在节理交线处未出现奇异节点，表明了该方法的正确性。 , , 关键词：节理 网格划分 颗粒 1 前言 节理裂隙是地质体的特有属性，其发育状况、空间展布特征、节理间的拓扑关系、节理面 [1, 2] 。数值分析是研究节理 内的充填胶结情况，直接影响着地质体的宏观渗透特性及力学性能 裂隙系统复杂物理力学特性的有效手段之一，而节理单元的网格剖分是进行数值分析的第一 步。 有限元网格生成技术可以按所生成的单元类型将它们划分为生成结构化网格的方法和生 成非结构化网格的方法[3] 。生成结构化网格的方法在边界上网格质量较低，不适用于节理面网 格的生成。生成非结构化网格的方法主要有布点及三角化方法，拓扑分解法，几何分解法， 基于栅格法几种。目前已经有利用野外露头上采集到的裂隙的观测数据，通过计算机处理最 终形成三维裂隙网络的人工几何模型的方法[4]，也已经有多种在节理网络图上生成有限元网格 [5, 6] [7, 8] [5] AFM(Advancing Front 的方法 以及生成节理岩体网格的方法 。其中文献 使用的是 Method) , AFM 法，首先将具有节理网络的计算域转化成单连通域 然后在单连通域内按 法生 成三角形单元。这种方法在节理面附近的网格质量很差，而且目前还不能用于三维问题。文 [7] 献 采用的是流形切割的方法，将三维有限元网格覆盖到物理模型上，然后对覆盖网格进行 剖分，形成最后的网格。这种方法可以用于三维问题，但是切割出的网格质量无法控制。 Bagi 提出了用大小不同的颗粒填满一个空间的方法[9]，Jianfei Liu 等人提出了用球形填充 (Sphere packing)来构造有限元网格的方法[10] 。受这两种方法启发，本文提出了一种利用颗粒生