暴雨诱发滑坡的高效无网格大变形数值模拟-Core.PDF

第27 卷 第4 期 应 用 力 学 学 报 Vol.27 No.4 2010 年 12 月 CHINESE JOURNAL OF APPLIED MECHANICS Dec.2010 文章编号：1000- 4939(2010) 04-0755-06 暴雨诱发滑坡的高效无网格大变形数值模拟 1 1,2 王东东 李凌 (厦门大学 361005 厦门) 1 (广东省水利水电科学研究院 510610 广州) 2 摘要：渗流与土体的大变形耦合损伤破坏问题是暴雨诱发滑坡数值模拟的重要内容，目前对滑坡 的全过程大变形非线性损伤破坏的仿真分析仍是一个难点。本文基于稳定节点积分再生核质点伽 辽金无网格法和大变形增量型本构关系，提出了一套能模拟暴雨诱发非饱和土质边坡大变形损伤 破坏的高效分析方法。该方法应用拉格朗日稳定节点积分，保证了数值积分的稳定性和高效性； 同时采用非饱和土的物理力学性质来模拟降雨过程中边坡土体性质的演化，并以各项同性损伤函 数和 Drucker-Prager 屈服准则反映土体的弹塑性损伤状态。该方法的非局部近似特性可以有效地 模拟从剪切带起始到发生滑坡的整个大变形损伤破坏过程。通过典型算例验证了方法的有效性， 并研究了滑坡发生时间与降雨强度和初始饱和度的关系。 关键词：暴雨；大变形；滑坡；无网格法；稳定节点积分 中图分类号：O242.21 文献标识码： A 的十几年间有了快速的发展，其近似函数具有不依 1 引 言 赖于节点之间的有序连接的特性，因此能有效地缓 解大变形时网格畸变这一问题，对于传统有限元法 我国的滑坡灾害频繁发生，其中多数系暴雨引 求解较困难的问题，譬如大变形的网格扭曲、动态 发，每年都造成重大的经济损失和人员伤亡[1-2] 。为 裂纹扩展和应变局部化的剪切带等多尺度现象，有 了有效预防滑坡灾害，边坡分析设计中采用了实验 着独特的优势[13]-14] 。 法、解析法和数值仿真等多种研究方法。实验法和 本文所研究的无网格法是采用移动最小二乘 解析法受限于模型的尺寸和流固耦合的复杂性，无 Moving Least Square MLS （ ， ）或再生核 法具备较广泛的适用性；而数值方法可以模拟丰富 Reproducing Kernel RK （ ， ）近似的伽辽金无网格 的地质条件，工作量也比实验方法小得多，近年来 法，即EFG[15]或RKPM[16] 。这是目前无网格方法中 得到了广泛应用[3-7] 。文献[8-10]在液相-固相耦合方 应用较多的两种方法，具有稳定性好、精度高等优 面开展了系统的研究工作，文献[11-12]对多孔介质 点，当采用多项式基函数计算形函数时，两者是等 材料的应变局部化理论进行了深入研究。然而边坡 价的。尽管伽辽金无网格法有着诸多优点，却还没 稳定性仿真往往涉及大变形情况下的流固耦合模 有得到大规模的工程应用，一个主要原因是伽辽金 拟，受拓扑网格的限制，对于边坡损伤破坏时材料 无网格法大都基于高阶的高斯积分，计算效率较 的大变形非线性分离，有限元法模拟起来相当困 低，并且不能严格满足分片试验的要求，无法充分 难，精度相对较低，也不易进行自适应分析。另一 发挥无网格法的优势，因此也没有在岩土类材料的 方面