三维块体离散元可变形计算方法研究田振农.pdf

第 27 卷 增 1 岩石力学与工程学报 Vol.27 Supp.1 2008 年 6 月 Chinese Journal of Rock Mechanics and Engineering June ，2008 三维块体离散元可变形计算方法研究 1 1 1 1，2 田振农 ，李世海 ，刘晓宇 ，柳丙善 (1. 中国科学院 力学研究所，北京 100080 ；2. 中国石油大学 储运工程系，山东 东营 257062) 摘要：在三维刚性块体离散元的基础上，研究可变形块体离散元方法，并编制相应的可变形块体离散元程序。该 方法把块体单元看作弹性体，并对其划分有限元网格，对每个可变形单元的计算是假设在力边界条件下求解可变 形单元的节点位移；可变形块体之间力和位移的传递是通过构筑的节理单元实现的，节理单元反应岩土工程中实 际节理的变形性质和强度特征，结构面破坏前，体现了连续介质计算方法的特点，结构面如果发生破坏，预设节 理处的节理单元直接转变化为法向和切向弹簧，体现了离散元计算方法的特点。另外，节理单元有实际的厚度， 可以根据实际节理的性质确定计算参数，这一点也体现出了对刚性块体离散元方法的改进。与刚性块体离散元方 法比较，该方法可以反映岩体变形性质，如泊松效应、三维岩体中的真实波速，它还可以模拟结构面的破坏演化 过程；与有限元方法比较，它可以更好地描述岩体中节理的特征，允许单元有较大的位移；相对于以前的可变性 离散元方法，它既克服了有限差分格式带来计算上的复杂，又建立在三维模型的基础上，为有关岩土工程的数值 分析提供了一个可行的数值方法。 关键词：数值分析；离散元法；岩体；可变形单元 分类号：O 241 文献标识码 A 文章编号 1000 –6915(2008)增 1 –2832 –09 RESEARCH ON DEFORMABLE CALCULATION METHOD BASED ON THREE-DIMENSIONAL BLOCK DISCRETE ELEMENT TIAN Zhennong 1 1 1 1，2 ，LI Shihai ，LIU Xiaoyu ，LIU Bingshan (1. Institute of Mechanics ，Chinese Academy of Sciences，Beijing 100080，China；2. Department of Storage and Transportation Engineering，China University of Petroleum，Dongying，Shandong 257062，China) Abstract ：Deformable block discrete element model is studied on the basis of three-dimensional rigid block discrete element method，and corresponding program is developed for this method. Each block element is considered as an elastic body，and finite element grids are generated. Each node displacement in the deformable block is calculated