离散单元法程序研制和应用.pdfVIP

8 西部探矿工程 2010年第2期 离散单元法程序研制和应用 柏 署 ，张 磊 ，杨 燕。 (1．湖南交通规划勘察设计院，湖南长沙410008；2．广州市中交四航局第二工程有限公司，广东 广州 510300； 3．湖南城建职业技术学院长沙校区，湖南长沙410015) 摘 要：基于离散单元动态松弛法，详细阐述了在 C++语言下编制程序中数据结构、块体生成和接 触判断三大难点，对应的提出了新的算法来解决问题，通过算法分析和流程图对算法进行实现。比较 了有限元和离散单元法在程序编写中的差异，并给出了算例说明离散单元法在岩土工程中的优势。 关键词：离散单元法；动态松弛；接触判断；数据结构；有限元 中图分类号：TP39 文献标识码：A 文章编号：10O4—5716(2O10)O2一OOO8一O5 在岩石力学的数值研究中，现阶段常用的方法有以 方法更适合岩土工程这一特殊领域。目前，离散元理论 有限单元法为代表的连续模型方法和以离散单元法为 日趋成熟，在国外也有相关软件[2]，但国内离散元应用 代表的离散模型方法。由于岩、土体的成因和构造复 软件系统还没有商品化，研究者基本都是应用国外的软 杂，岩性多样，而且其 内部存在大量的地质界面，如层 件，偶尔有学者编制程序，但都是用来为某一种问题服 面、片理、节理、软弱夹层，以及断层破碎带等。正是因 务，具有很强的针对性。这种现象严重制约了离散单元 为它们的存在，导致应力或位移不连续。使得连续模型 法的应用和推广，对离散元方法的计算软件研制，不仅 方法在应用中有很大的局限性，虽然Goodman等人[1] 对工程应用的需求，同时对该方法本身的发展与成熟， 在反映位移连续的有限元模型中引入节理单元，创立了 都有非常重要的意义。 节理单元法。在一定程度上可以描述节理、断层和断层 1 离散单元法简介 带、破碎带以及其他软弱结构面的不连续性质，然而节 离散元单元方法是 (DiscreteElementMethod)是 理裂隙过多时，由于过多的节理单元会导致数值计算的 Cundall，P．A 于 1971年所提出的一种不连续数值计 不稳定，而且对大变形也难 以处理。相比之下，离散元 算方法。传统的离散单元法有动态松弛法[3]和静态松 表 1 基坑变形观测成果表 参考文献： 1993． [1] 余志成，施文华．深基坑支护设计与施工[M]．北京：中国 建筑工业出版社，1997． 收稿 日期：2009-05-08 [2] 高大钊，孙钧，陈忠汉，黄书秩，程丽萍．深基坑工程[M]． 作者简介：胡胜华(1972-)，男 (汉族)，工程师，现从事岩土工程勘察、施工 北京 ：机械工业出版社 ，1999． 及地质灾害勘查与设计等工作。 [3] 林宗元．岩土工程治理手册[M-1．辽宁科学技术出版社， 2010年第 2期 西部探矿工程 9 弛法 剖，后来在算法上得到改进的同步松弛方法 引。 但 目前，所有的算法当中，最流行和应用最成熟的仍是 动态松弛法，所以通常意义上称的离散单元法也就指动 态松弛法。其实质是对临界阻尼振动方程进行逐步积 分。这种带有阻尼项的动态平衡方程，利用有限差分法 按时步迭代求解，是一种显式解法，所以不需要求解大 型矩阵，不但节省了大量计算工作，而且允许块体累积 发生大的平移和转动，克服了连续体模型的小变形假 设，可以用来求解一些非线性问题[6]。 力一位移方程为： — K × (1) — K ×U