元宝湾矿3煤开采顶板FLAC数值模拟.doc

工程地质模型示意图模型尺寸：长420m，高度138m，单位宽度5m。开采煤层为31煤层，主要模拟采动顶板的冒裂带高度。表 地质模型各岩层物理力学参数参考数据岩石名称弹性模量（MPa）泊松比粘聚力（MPa）内摩擦角（）抗拉强度（MPa）密度（kg/m3）泥岩6000.282271.62100粗粒砂岩26000.206.8215.42700煤层3000.30.8220.651700中砂岩20000.206.4205.22600砂质泥岩7000.282.2271.92300本次模型采用Mohr-Coulomb本构关系。为此，模型材料选项中需要输入密度、体积模量K、剪切模量G、内聚力c、内摩擦角以及抗拉强度t。数值模拟过程中用到的参数，体积模量K、剪切模量G，可根据下面的公式进行换算： 式中： K为体积模量，G为剪切模量、E为弹性模量，v为泊松比。new;-model-gen zone brick p0 0 0 0 p1 420 0 0 p2 0 5 0 p3 0 0 34 size 60 1 17 gro 1gen zone brick p0 0 0 34 p1 420 0 34 p2 0 5 34 p3 0 0 80 size 60 1 46 gro 2gen zone brick p0 0 0 80 p1 420 0 80 p2 0 5 80 p3 0 0 138 size 60 1 29 gro 3gro 中粗粒砂岩 ran z 0 12 gro 1gro 砂质泥岩 ran z 12 19 gro 1gro 4煤 ran z 19 25 gro 1gro 砂质泥岩 ran z 25 28 gro 1gro K3粗砂岩 ran z 28 34 gro 1gro 砂质泥岩 ran z 34 50 gro 2gro 3煤 ran z 50 54 gro 2gro 砂质泥岩 ran z 54 80 gro 2gro K4粗砂岩 ran z 80 85 gro 3gro 砂质泥岩 ran z 85 106 gro 3gro 粗粒砂岩 ran z 106 114 gro 3gro 砂质泥岩 ran z 114 138 gro 3pl bl group;-岩石属性-def setup zs_dens=2800 m_dens=1400 sn_dens=2580 k3_dens=2780 k4_dens=2780 cs_dens=2780 zs_fric=29 m_fric=22 sn_fric=20 k3_fric=26 k4_fric=26 cs_fric=29 zs_coh=7e6 m_coh=1.2e6 sn_coh=2e6 k3_coh=3e6 k4_coh=3e6 cs_coh=7e6 zs_dil=0 m_dil=0 sn_dil=0 k3_dil=0 k4_dil=0 cs_dil=0 zs_tens=34e6 m_tens=1e6 sn_tens=2e6 k3_tens=34e6 k4_tens=34e6 cs_tens=24e6 zs_E=3e9 m_E=2e7 sn_E=2e8 k3_E=2e9 k4_E=2e9 cs_E=3e9 zs_v=0.22 m_v=0.32 sn_v=0.26 k3_v=0.24 k4_v=0.24 cs_v=0.23 zs_bulk=zs_E/(3*(1-2*zs_v) m_bulk=m_E/(3*(1-2*m_v) sn_bulk=sn_E/(3*(1-2*sn_v) k3_bulk=k3_E/(3*(1-2*k3_v) k4_bulk=k4_E/(3*(1-2*k4_v) cs_bulk=cs_E/(3*(1-2*cs_v) zs_shear=zs_E/(2*(1+zs_v) m_shear=m_E/(2*(1+m_v) sn_shear=sn_E/(2*(1+sn_v) k3_shear=k3_E/(2*(1+k3_v) k4_shear=k4_E/(2*(1+k4_v) cs_shear=cs_E/(2*(1+cs_v) endsetup m m prop bulk zs_bulk shear zs_shear coh zs_coh fri zs_fric tens zs_tens range group 中粗粒砂岩 prop bulk m_bulk shear m_shear coh m_coh fri m_fric tens m_tens range group 4煤 prop bulk m_bulk shear m_shear coh m_coh fri m_fric tens m_tens range group 3煤 prop bulk sn_bulk shear sn_shear coh sn_coh fri sn_fric tens sn_tens range group 砂质泥岩 prop bulk k3_bulk shear k3_shear coh k3_coh fri k3_fric tens k3_tens range group K3粗砂岩 prop bulk k4_bulk shear k4_shear coh k4_coh fri k4_fric tens k4_tens range group K4粗砂岩 prop bulk cs_bulk shear cs_shear coh cs_coh fri cs_fric tens cs_tens range group 粗粒砂岩;-边界条件-fix z ran z -0.1 0.1fix x ran x -0.1 0.1fix x ran x 419.9 420.1fix y ran y -0.1 0.1fix y ran y 4.9 5.1 apply szz -2.5e6 ran x 0 300 y 0 5 z 137.9 138.1set grav 0 0 -9.8hist unbalsolve ratio 1.0e-3save ini-model.savrest ini-model.savini xd 0 yd 0 z