岩石蠕变-渗流耦合作用规律研究.pptVIP

目前研究不足 1、岩石蠕变-渗流耦合作用规律，尤其进行系统地物理试验和数值模拟研究方面很少； 2、岩石蠕变-渗透率全过程演化规律研究未曾清楚。 3、蠕变破裂判断准则及损伤分析和工程应用研究不足； 4、未基于软件平台进行编制蠕变破裂和蠕变-渗流耦合有限元模块进行工程应用； 5 层状边坡蠕变-渗流耦合数值模拟 海州露天矿国家矿山公园，东西长3.9Km，南北宽1.8Km，深达350m，面积近7.02 Km2。从1953—2009年之间共发生了90余次滑坡，平均每年发生1.6次滑坡，其中有36次因软弱夹层导致滑坡；降雨，地表水入侵诱发29次；夹层具有明显的流变特性，又在水的影响下，导致露天矿边坡滑动。为保证矿山公园的顺利建设和游客的生命安全，所以对海州矿露天层状边坡的稳定性分析和蠕变破裂-渗流耦合作用数值模拟研究尤为必要。 海州露天矿国家矿山公园 5 层状边坡蠕变-渗流耦合数值模拟 42.48 109.26 33.07 1.15 24.36 0.286 44.35 111.49 基 岩 4.28 11.94 28.00 0.5 20 0.323 4.00 11.82 黄 土 5.00 13.96 17.75 0.175 17.99 0.395 5.00 13.96 4#5#6# 36.39 94.59 33.07 1.15 24.36 0.286 37.13 95.50 8#上岩 35.02 90.08 33.07 1.15 24.36 0.286 35.74 91.92 粉砂岩 1.23 4.42 33.07 1.15 24.36 0.395 1.25 4.49 9#上岩 16.80 44.21 17.99 0.286 17.14 44.09 9#弱层 4.9 13.96 19.71 0.45 17.99 0.286 5 13.96 8#弱层 Ks /(kN·m-3) Kn/(kN·m-1) φ/(°) C/105Pa /(kN·m-3) /MPa /MPa 岩石 名称 （1）物理力学参数 稳定性分析 层状边坡蠕变破裂数值模拟 0.187 10.61 5 层状边坡蠕变-渗流耦合数值模拟 （2）有限元模型 剖面E20，有阜新电厂、海州矿博物馆和指挥部，在不影响计算结果的前提下，模型经过简化如图所示。宽1 290 m，高345 m，下边界标高为▽-150 m，坑底标高为▽-90 m。图形显示模型的横向与纵向比例为1:1，标高如图所示。 +15m +80m -90m 5 层状边坡蠕变-渗流耦合数值模拟 （3）稳定性分析结果 由 FOS=1.857下降到FOS=1.526。导致边坡滑动的原因主要为6#、8#和9#弱层，其中8#弱层为主要滑移面，边坡总体相对稳定，由于滑移面与地面贯通，若水入渗到各弱层使得其随时间的推移整个边坡会沿弱层发生蠕动，导致整体滑坡，所以建议对6#、8#和9#弱层对应地面处采取相应的隔水处理措施。所以海州矿边坡的稳定与否，与排水系统和防水措施正常运转关系密切。 5 层状边坡蠕变-渗流耦合数值模拟 受法向力作用下蠕变破裂过程模拟 （1）应用FEPG软件，编制相应的蠕变破裂程序，在进行有限元分析前，确定离散化网格时，充分考虑岩质边坡形成过程中的各种力学变化，使得在煤岩体流变数值分析过程中适时适量增量加载、逐步改变研究体系结构，计算模型尺寸如图。划分11 716个单元、6 054个结点。 （2）流变模型参数 （3）蠕变破裂模拟 5 层状边坡蠕变-渗流耦合数值模拟 时步15个月 时步30个月 时步45个月 每一个弱层蠕变破裂都是在其顶部与表土接触的地方出现裂缝。由于岩石抗压不抗拉，且弱面存在，在重力的作用下，弱面顶部首先发生蠕变拉破坏，产生裂纹，形成松散区；再有其他作用力时，会在弱面形成松散区的部位破坏；有最大主应力的地方，由于其强度比较高，虽然最大，但是并没有达到材料的抗拉强度。 5 层状边坡蠕变-渗流耦合数值模拟 层状煤岩体的长期变形过程中，岩层离层蠕变断裂过程可认为是在变形过程中岩层层面受拉破坏或受剪滑移并沿层面切向方向扩展和发育的复杂过程。因此，应用复合型最大拉应力理论或剪切破坏理论，即当岩层或节理变形过程中出现最大拉应力值大于岩石抗拉强度或节理壁面抗拉强度或节理面上的剪应力大于节理剪切强度时，认为岩石单元或节理单元沿该方向被拉坏或剪坏，此方向原来所承受的拉应力或剪应力被释放，转变为等效过量节点力，该方向刚度趋近于零。 离层蠕变断裂过程模拟 计算是分两个阶段进行,首先要进行试计算,得出各个结点的应力数值，根据前处理所获数据及计算模型形成初始地应力场。第二步 就是应用程序进行计算,判断应力