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

水平位移： 20000 步 水平位移： 30000 步 水平位移： 50000 步 水平位移： 80000 步 水平位移： 90000 步 一、不连续地质体建模技术 二、裂隙岩体渗流应力耦合模型 三、裂纹扩展模拟技术 ? 单元界面破裂 ? 块体内部破裂 主要内容 Incomplete fracture geometry information in field surveys 裂纹扩展模拟技术 —— 界面破裂 Ks ： 剪切刚度 Kn ： 法向刚度 ? 当相邻两块体之间完全连续时，邻居节 点采用连接型弹簧，遵循线弹性本构关系 ? 当连接型弹簧满足张拉 - 压剪复合破坏准 则时，连接型弹簧将发生张拉或压剪模式 的弹脆性破坏，转变为接触型弹簧 其中，张拉模式的弹脆性破坏采用最 大张拉力准则模拟 c t n A F ? ? 压剪模式的弹脆性破坏采用莫尔 - 库 仑准则模拟 c n s A c F F 0 tan ? ? ? ? 当相邻两块体之间存在结构面时，邻居 节点采用接触型弹簧，遵循非线性本构关 系 块体 3 块体 1 块体 2 Incomplete fracture geometry information in field surveys 界面破裂算例 顺层岩体 反倾岩体 含软弱带岩体 软弱土层 一、不连续地质体建模技术 二、裂隙岩体渗流应力耦合模型 三、裂纹扩展模拟技术 ? 单元界面破裂 ? 块体内部破裂 ? 弹脆性材料（岩体） ? 弹塑性材料（土体） 主要内容 ? ? ? ? Incomplete fracture geometry information in field surveys 裂纹扩展模拟技术 —— 块体内部破裂 自适应网格模拟技术 扩展有限元技术 块体的破裂模型 C 3 ? 拉破坏面 3 ? 3 ? ? 拉伸破坏 通过相对位移选择滑动面： 3 ? 3 ? 3 ? 1 ? 1 ? 1 ? 1 ? t ? ? ? 1 0 tan 0 ? ? ? ? ? ? n c ? 剪切破坏 常规单元切割 切割点与角点重合 裂纹扩展切割 边界单元切割 块体一条边 与切割面关系 不相交 被剖分为两部分 该边至少有一个角 点在切割面上 块体一个面 与切割面关系 不相交 该边两角点 P1 、 P2 将继承父节点与 其它块体节点的混合节点关系 该面面中心点将继承父面中心点与其 它块体节点的混合节点关系 被剖分为两部分 边被完全切割 边被不完全切割 边界边 边被完全切割 边被不完全切割 边界边 面被完全切割 面被不完全切割 边界面 Key point ? Block cutting technique ? ? ? ? ? ? ? ? ? ? ? ? ? 0 ) ( 0 ) ( 0 ) ( in on in x x x ? ? ? 对于给定的切割平面方程，给出切割面 Γ 、切割面上方 Ω + 和 切割面下方 Ω - 的 level set 函数 ? ? ? ? ? 0 ) ( min ) ( min : 0 ) ( min ) ( min : ? ? ? ? ? ? ? ? i I i i I i i I i i I i el el el el ? ? ? ? 不切割块体 切割块体 ) ( 1 2 2 1 1 1 x x x x ? ? ? ? ? ? ? 块体边与切割面的交点坐标： 利用 level-set 函数遍历块体所有面，检查各面与切割面之间 的关系 常规单元切割 插值关系 混合节点 混合节点 接触节点 1 。寻找混合节点； 2 。寻找接触节点； 3 。寻找插值节点； 边界单元切割 接触节点 混合节点 混合节点 接触 节点 接触节点 自由节点 1 。寻找混合节点； 2 。寻找接触节点； 3 。寻找插值节点； 裂纹扩展切割 混合节点 接触节点 或 混合节点 1 。寻找混合节点； 2 。寻找接触节点； 3 。寻找插值节点； 混合节点 块体一条边被切割面剖分为两部分 1 。建立新节点 S1 、 S2 与角点 P1 、 P2 之间的插值关系； 2 。将位于切割点的所有新节点定义为克隆关系； Y 该边两角点 P1 、 P2 将继承父节点与其它块体节点的克隆关系 该边是否为边界边 Y 将新产生的两个节点 S1 、 S2 定义为接触关系； N 该边是否为完全被切割 N Y 1 。消除原有新节点的插值关系； 2 。 将每一对新节点（ S1,S2 ）都定义为接触关系； 3 。将位于接触面上方的所有新节点定义为克隆关系； 4 。将位于接触面下方的所有新节点定义为克隆关系； 块体一个面被切割面剖分为两部分 将面心节点看作自由节点 Y 建立面心节点与相邻面内节点之间