多层急倾斜煤岩层中的集中煤仓三维有限元仿真分析.docxVIP

多层急倾斜煤岩层中的集中煤仓三维有限元仿真分析李永和( 上海大学)王更雨 汉宁明( 窑街矿务局)摘要 以窑街矿务局三矿多层急倾斜煤岩层中的煤仓与巷道交叉处复杂结构为研究对象 多层急倾斜煤岩层中的集中煤仓 三维有限元仿真分析 李永和 ( 上海大学) 王更雨 汉宁明 ( 窑街矿务局) 摘要 以窑街矿务局三矿多层急倾斜煤岩层中的煤仓与巷道交叉处复杂结构为研究对象 , 应 用三维非平行岩层的有限元模型和 CAD 仿真技术 , 对该煤仓各典型部位的内力和变形进行了系 统的研究 , 得到了许多可视化非常直观的计算结果 , 这对煤仓或交叉处采取进一步的维护和加固 具有重要的指导意义 . 关键词 非平行岩层 煤仓交叉处 三维有限元仿真 中图分类号 TD 32 在采矿工程中 , 多层非平行岩层中的煤仓及巷道与煤仓交叉处这类大型复杂结构 , 由于开挖引起的内 力和变形及其稳定性的研究一直是地下工程中的重要课题 . 目前对这类复杂问题多采用平面有限元法进行 处理1 ～4 , 但很难真实反映巷道与煤仓交叉处的平面假设 , 因为交叉处是一个完全的三维问题 . 另外 , 传统的有限元虽然能给出内力和位移各种表现结果 , 但是很难给出研究对象任一部位可视化非常直观的结 果3 ,4 . 针对这些缺陷 , 本文对窑街矿务局三矿多层急倾斜煤岩层中煤仓及其交叉处这种复杂的实际工程 进行了三维有限元 CAD 仿真分析 , 给出了非常受现场工程技术人员欢迎的 、可视化非常清楚的结果 , 进一步搞清这类复杂结构的内力分布及其指导结构的加固与维护有着重要的指导意义. 对 有限元计算模型 1 111 非平行煤岩层的本构方程 为了能够分析多层非平行煤岩层煤仓结构 , 首先建立了如图 1 所示的层状模型. 这一模型是把岩层中的任一层看作是由相互平行 的子层组组成 , 其厚度为 d , 子层间的岩石为各向同性 , 可用杨氏 模量 E 和泊松比ν来描述. 在这个模型中 , 由于 y′方向变形远小于 其它 2 个方向的变形 , 故可通过变形叠加求得柔度矩阵关系为 ε = Cσ , (1) 式中 , n 图 1 层状模型 Fig11 L ayered model C = Cr + ∑T C T T ; ( 2) j ij j i = 1 1 - ν2 - ν( 1 + ν) 0 - ν( 1 + ν) 1 - ν2 0 0 0 2 ( 1 + ν) Cr = ; ( 3) 收稿日期 : 1998 - 10 - 16 1 E 煤炭学报1999 年第 24 卷1380000L Kn0 sin2β co s2β00L KsCi j=1;( 4)1co s2βsin2β- 2si 煤 炭 学 报 1999 年第 24 卷 138 0 0 0 0 L Kn 0 sin2β co s2β 0 0 L Ks Ci j = 1 ; ( 4) 1 co s2β sin2β - 2sin βco s β sin βco s β - sin βco s β Tj = ; ( 5) co s2β - sin2β 2sin βco s β ε= [εx , εz , εx z ] T , σ = σx , σz , σx z ] T ; Cr 为单独的完整土层柔度矩阵 ; Ci j 为单独的子层系统柔度 矩阵 ; Tj 为坐标变换矩阵 ; Ks 和 Kn 分别为单独层 (或子层) 的切向和法向刚度. 112 层状岩层弹塑性本构模型 在荷载水平 相 对 比 较 高 时 , 任 一 岩 层 均 可 看 作 弹 塑 性 问 题 , 其 屈 服 条 件 和 屈 服 面 可 采 用 Drucker2 Prager 屈服准则 3 ～5 : J 2 - K = 0 , J 2 为 应 力 偏 量 的 第 二 不 变 量 , J 2 = 2 S kk ; ( 6) a = f = a I 1 + 式中 , I 1 为 应 力 张 量 的 第 一 不 变 量 , I1 = σkk ; S ij tan φ 9 + 12tan2φ; K = 3 C 9 + 12tan2φ. 其中 : K 可用内聚力 C 和摩擦角φ 来描述. 参数 a , 在这种情况下 , 合理的本构关系可用增量形式表示为 dσ = D p dε. 对于增量应力 - 应变矩阵 D p 可用下式表达 , 即 (7) D p = D ( I - a r A T - a A T - , ?, - a A T) , (8) r 1 1 n