开挖引起的煤岩变形及界面破裂Abaqus有限元数值模拟.docx

开挖引起的煤岩变形及界面破裂Abaqus有

限元数值模拟

侯雪沈新普

（沈阳工业大学计算力学研究所，沈阳市沈辽西路111号574信箱，邮编：110870）

摘要：本文应用ABAQUS大型有限元计算软件，模拟开采过程中，倾斜煤层处于不同倾角

时，煤岩变形及上覆岩石与煤层界面的损伤失效情况。计算中采用了Drucker-Prager条件

作为煤岩的屈服条件，TractionSeparateLaw准则作为上覆岩石与煤层之间界面单元的损

伤判定标准。数值模拟结果对煤层在开采过程中所导致的煤岩破裂与瓦斯突出等矿山灾害的

预防和防治具有重要的参考价值。

关键词：煤岩变形；界面单元；数值模拟；Abaqus

1引言

煤是我国经济发展建设中不可缺少的主要能源之一，随着煤矿开采规模的

不断扩大以及开采深度的增加，煤岩开采导致的煤岩破裂与瓦斯突出的预防和治

理是目前煤炭行业的一大难题。虽然国内外学者对综放过程煤岩移动及顶煤破坏

规律等问题进行了大量的理论分析与试验研究,取得了一定的成果,但这些研究

工作主要是针对近水平煤层进行的,通常地下煤层都具有一定的倾斜角，所以研

究倾斜煤层在开采过程中煤岩的应力位移变化情况，对预防和防治开采过程中引

起的煤岩破裂与瓦斯突出等矿山灾难具有重要的意义。

由于目前监测手段的局限性，对于煤岩在开采过程中的变形特征研究现场监

测和相似材料模拟试验均不能获取分析问题所需的足够信息。随着计算机技术的

发展，数值模拟已成为分析地下开采引起的岩层移动及地表变形过程的常用手

段。段东等学者用RFPA2D一GasFlow软件建立数值模拟模型，分析了水平煤岩在

开采过程中的应力位移变化[1]；信晓东利用数值模拟的方法研究了急倾斜煤层在

开采过程中煤岩的变形特征[2]；王树仁，王金安等应用FLAC3D计算程序，分析

了大倾斜角煤层在综合开采过程中，采场覆岩与煤层的运移及相应的应力场、破

坏场的变化特征[3]；熊祖强，王文等学者采用UDEC离散元数值模拟研究方法对

远程保护层开采的上覆没岩层运移规律进行了研究[4]。

本文在参照前人大量研究的基础上，为了较真实的模拟倾斜煤层在开采过程

中的应力位移变化情况，应用ABAQUS大型有限元计算软件，建立二维有限元模

型，模拟埋深1000m的倾斜煤层在倾角分别为10°，20°，30°，45°时，开

挖进行过程中，煤岩变形及上覆岩石与煤层界面的损伤失效情况。计算中煤岩采

用Drucker-Prager条件作为煤岩的屈服条件，TractionSeparateLaw准则作

为上覆岩石与煤层之间界面单元的损伤判定标准。数值模拟结果对深入研究煤岩

破裂与瓦斯突出等矿山灾害的预防和防治具有重要的参考价值。

2有限元模型

将倾斜煤层四个角度的模型简化为二维平面模型，以煤层倾角20°模型图

为例，如图1，整体模型宽100m,高200m，煤层与岩石之间设置一层几何厚度为

0.1m的界面单元，煤层实际埋深约1000m，模型上层顶板承受约1000m的岩石自

重压力。

图1煤层倾角20°时整体模型图

3模型材料属性

计算中煤岩采用Drucker-Prager屈服准则，并且在开挖过程中激活了

Drucker-Prager中的Creep效应，使Drucker-Prager塑性与蠕变耦合，蠕变曲

线如图2所示[5]，蠕变法则采用“时间硬化”幂函数定义为：

??cr?A?ntm(1)

eq

其中：ε?cr为等效蠕变应变率张量，

?n为等效蠕变应力，t是时间，模型参

eq

数A,n和m经过计算转换[6]，取值为A?2.5?10-20,n?3,m?-0.2。具体煤

层的力学参数如表1：

1000

800

600

400

200

0246810121416

t/h

图2文献中蠕变试验曲线

表1煤层的力学参数

材料属性弹性模量（GPa）泊松比密度（kg/m3）摩擦角剪胀角

煤层2.70.3715002020

上部岩石8.00.2725803030

下部岩石8.00.2725803030

界面单元采用TractionSeparateLaw曲线作为损伤判定标准如图2，最大

名义应变及线性损伤演化来判定单元的损伤发生和失效过程，如公式（2）和公

式（3）

?????

maxn,s,t

??

?o?o?o

??

nst

?

(2)

?,,分别为法向方向名义应变，第一剪切方向和第二剪切方向名义

n??式中

st

应变；?n,?,?o分别为：界面单元发生变形后的法向方向应变，第一剪切方向和

oo