三维模态变形体离散元方法及应用研究.pdfVIP

水电站设计 第26卷第3期 D H P S 2010年 9月 三维模态变形体离散元方法及应用研究 张 冲 ，金 峰 (1．中国水电顾问集团成都勘测设计研究院，四川成都 610072；2．清华大学 水利系，北京 100084) 摘 要：首先建立 了三维模态变形体离散元 (3MDEM)方法，阐述了其基本原理及特点，并应用一系列算例进行了验证。这些算例 涵盖了静力到动力、小变形到大位移、连续到非连续等各种不同情况，以综合验证模态变形体离散元实际的计算能力。结果表明， 3MDEM在小变形 以及连续介质力学领域可以给出和有限元相媲美的应力和位移结果，而在大变形、大位移强非线性领域可 以给 出与3DEC相媲美的计算结果。最后本文对 1981年梅花拱坝的溃坝过程进行仿真研究。通过分析得知，梅花拱坝 由于周边缝施 工设计不当，沿周边缝抗剪强度过小，导致蓄水后坝肩推力促使坝体沿周边缝发生整体上滑，拱坝失去支撑最终坝体断裂并溃决。 关键词：离散元(DEC)；非连续体 ；数值方法；大变形；仿真分析 中图法分类号：O241 文献标识码 ：A 文章编号：1003—9805(2010)03—00O1—09 边界作为不连续界面通过接触关系检索算法确定接 1 前 言 触关系 j，并根据线性或非线性接触本构模型规定 界面上的力和位移关系；块体的空间运动分解为刚 离散元法是2O世纪70年代初 由Cundall…提 性运动和变形的叠加，块体 内部进一步采用多种变 出的一种分析不连续岩体运动与变形的有效数值分 形模态将块体的变形分解为数种变形模态的叠加。 析方法。它建立在牛顿第二定律基础之上，可以分 计算时则使用显式时步步进 的动态松弛解法求解平 析刚体、变形体的力学行为，块体之间的力通过块体 动方程、转动方程和变形方程 ，从而得出块体的运动 间接触作用相互传递，块体之间可以接触或分离不 和变形。接着，再使用线弹性的或者其他本构模型 需要满足变形协调条件 ，适合于求解弹塑性、大变形 求解块体的应力。 和动力问题。在静力和动力问题求解过程中离散元 这种数值方法是一种可以实现从连续介质到非 采用时步步进 的动态松弛法，可以模拟结构的变形 连续介质、从小变形到大位移全过程仿真的分析方 和破坏过程 ，实现过程的仿真。之后 Cundall对离 法。同时该方法拥有较高的计算效率，不需要在块体 散单元方法不断改进和完善，发展了变形体的离散 内部进一步差分四面体单元，以六面体块体为例，当 元法 。J，并与 Itasca公司合作开发了二维离散元 采用21种变形模态离散块体的变形时(考虑N阶 大型软 件包 UDEC_4和三 维离 散元 大 型软 件 的变形)，再加上块体的平动和空间转动，每个块体仅 3DEC_5J， 为离散元方法的工程应用奠定了基础。 仅有 27个 自由度，相当于3DEC在六面体块体内部 目前采用的变形体离散元都是在块体内部进一 增加一个节点所带来的自由度，因此计算效率较高。 步进行四面体差分网格划分，在形成四面体之后，将 本文在给出方法本身之后，选取了若干有代表 质量、力等物理量分配到四面体的节点上。计算时 性的实例进行了模型验证，并结合梅花拱坝破坏实 求解节点运动方程，从而得出节点的位移和单元的 例，充分阐述了该方法在水 电工程中的应用。 应变。对于大型工程 ，这种方法导致节点数急剧增 加，从而计算代价太大。 2 模态变形体离散元方法 本文基于DEM 的基本理念，结合空间多体柔 性动力学基本方法 J，建立了一套完备的