刘晓宇-CDEM流固耦合及块体破裂模型.ppt

CDEM的流固耦合及块体破裂模型 地质体的非连续性 地质体的流固耦合特性 地质体的渐进破坏 裂隙岩体渗流应力耦合模型 裂隙岩体渗流应力耦合模型 裂隙岩体渗流应力耦合模型 块体的破裂模型 弹塑性材料本构模型 弹塑性材料本构模型 弹塑性材料本构模型 强度折减方法计算边坡安全系数 问题由来 问题由来 不连续分叉条件 相关资料 张丽, 刘晓宇, 李世海. 裂隙岩体稳定/非稳定渗流数值模拟, 岩石力学与工程学报, 2009(9), 28, Supp.2 3409-3416. X.Y. Liu, Y. Liu and S.H. Li. Modelling of Hydro-mechanical Coupling for Fractured Rocks. Proceeding of International Symposium/RCUK-China Summer School on Discrete Element Methods and Numerical Modelling of Discontinuum Mechanics, pp 288-302, 2008/9/24. 韩永臣, 刘晓宇, 李世海. 模拟岩石材料脆性破裂过程的三维离散元模型, 力学与实践. 均质土坡，坡高H=20m，坡角?=30?，杨氏模量E=100MPa，泊松比?=0.3，容重?=20kN/m3，内聚力C=30kPa，内摩擦角?=15?。边坡左右两侧边界为法向约束，底边为双向固定约束。 Time Step: 120000 Time Step: 200000 Time Step: 240000 Time Step: 280000 Time Step: 320000 Time Step: 360000 水平位移：20000步 水平位移：30000步 水平位移：50000步 水平位移：80000步 水平位移：90000步 一、不连续地质体建模技术 二、裂隙岩体渗流应力耦合模型 三、裂纹扩展模拟技术 单元界面破裂 块体内部破裂 主要内容 Incomplete fracture geometry information in field surveys 裂纹扩展模拟技术——界面破裂 Ks：剪切刚度 Kn：法向刚度 当相邻两块体之间完全连续时，邻居节点采用连接型弹簧，遵循线弹性本构关系 当连接型弹簧满足张拉-压剪复合破坏准则时，连接型弹簧将发生张拉或压剪模式的弹脆性破坏，转变为接触型弹簧 其中，张拉模式的弹脆性破坏采用最大张拉力准则模拟 压剪模式的弹脆性破坏采用莫尔-库仑准则模拟 当相邻两块体之间存在结构面时，邻居节点采用接触型弹簧，遵循非线性本构关系 块体3 块体1 块体2 Incomplete fracture geometry information in field surveys 界面破裂算例 顺层岩体 反倾岩体 含软弱带岩体 软弱土层 一、不连续地质体建模技术 二、裂隙岩体渗流应力耦合模型 三、裂纹扩展模拟技术 单元界面破裂 块体内部破裂 弹脆性材料（岩体） 弹塑性材料（土体） 主要内容 Incomplete fracture geometry information in field surveys 裂纹扩展模拟技术——块体内部破裂 自适应网格模拟技术 扩展有限元技术 C 3 s 拉破坏面 3 s 3 s 拉伸破坏 通过相对位移选择滑动面： 3 s 3 s 3 s 1 s 1 s 1 s 1 s 剪切破坏 常规单元切割 切割点与角点重合 裂纹扩展切割 边界单元切割 块体一条边 与切割面关系 不相交 被剖分为两部分 该边至少有一个角点在切割面上 块体一个面 与切割面关系 不相交 该边两角点P1、P2将继承父节点与其它块体节点的混合节点关系 该面面中心点将继承父面中心点与其它块体节点的混合节点关系 被剖分为两部分 边被完全切割 边被不完全切割 边界边 边被完全切割 边被不完全切割 边界边 面被完全切割 面被不完全切割 边界面 Key point ? Block cutting technique 对于给定的切割平面方程，给出切割面Γ、切割面上方Ω+和切割面下方Ω-的level set函数 块体边与切割面的交点坐标： 利用level-set函数遍历块体所有面，检查各面与切割面之间的关系 常规单元切割 插值关系 混合节点 混合节点 接触节点 1。寻找混合节点； 2。寻找接触节点； 3。寻找插值节点； 边界单元切割 接触节点 混合节点 混合节点 接触节点 接触节点 自由节点 1。寻找混合节点； 2。寻找接触节点