受采动影响厚煤层沿顶巷道围岩稳定性模拟.pdfVIP

磊泰辜器镕餐蕞 中国煤炭学会第六届青年科技学术研讨会论文集 受采动影响厚煤层沿顶巷道 围岩稳定性模拟? 山东科技大学谭云亮王春秋 1 引言． 、 一 叫随着巷道围岩控制理论与技术的发展，人们越来越认识到对于厚煤层开采而肓，顺槽 巷道沿顶布置有利于顶板的稳定性，但对受采动影响时．巷道的底板及两帮的破坏与变形 的控制常常不被引起重视，为此。本文采用FEM与BEM耦台技术，对华丰矿2406上平 巷围岩的破坏及变形特征进行散值计算。) ／ 2 FEM与BEM耦合采法原理 如图1所示，把区域D(边界r)分成两部分。在m上应用FEM方法．在m上应用 BEM方法。设m与璐相连的内边界F．，区域珥由边界n+n包围，区域m由边界r】 +f、包围。在珥内及边界f、+r1上，采用FEM算法，则有刚度方程为： [K]fU)一f，)+fD)(1) 式中；D——与体力对应的节点力矢量。 在岛内及边界n+rl上，采用BEM算法，则BEM的 离散方程可以表示为 (2) [H]{u}一{G{{P)+{丑} 式中：丑——对应体力的载荷力量。 将方程(2)左乘[M][G]_。，并令暖’]=[M][G]-1[H]， 圈l FEM与BEM耦台 [村]=[G]“[B]=CD’]，[，]=[村][P]，则有： 葬法示意图 [K’]{U)=[F’]+[D’](3) 即将BE；M方程转挟戒等价的FEM方程形式。将式(1)、(3)联立即求得耦合方程。 这种耦合算法对于分折地下巷道稳定性显示出了突出的优点，它可以充分利用FEM 与BEM各自舯优势，即在巷道周围附近区域可以用FEM解算，以适应巷道复杂的材料 特性{而远处可以划分BEM来解算，它可以比较简单地处理地应力条件，从而改善整个 模型的计算效率。 3本构模型及算法 ·基童珥目，目家自髂科学基金项目和山东省基空项目赍助 14 谭云亮等： 受采动影响厚煤层沿顶巷逋围岩稳定性模拟 3．1本构模型 岩石材料是现有材料最为复杂的～种，对其本构关系的描述，国内外有许多种．每一 种本构模型都有它的优缺点。从工程应用的角度来讲，本构关系应既符合岩石材料特点， 又能使得模型中的物理力学参数易通过实验获得，简单应用。为此，特选取典型的M。hf． Coulomb准则作为判定巷道围岩是否发生破坏的准则．即： 卜扣升瓴cos一一兴蛐i旷c。婶=。 1 — ㈤ 5Q一鲁stna十!挚一。 式中：，，Q——屈服势函数和应变势函数； c——粘聚力； 旷一内磨擦角； a——剪胀角。 3．2算法过程 对于任何时步的不平衡载荷妒可由下式绐出： 妒。J。Brcrdv～R ㈣ 式中：a一应力张量{ R——当前时步的蓑荷总重。 计算选代过程为： 1)由弹性矩阵给出的位移增量和弹性应力增量为： gu3+j-[K膏2咖一1 如．=DB△鑫{ ㈨㈣ 2)由前一时步的位移增量经加速因子a修正后为： 甜．=8H醯。 ㈣