[发明专利]一种煤层开采覆岩沉陷三维半解析分析方法

说明书

技术领域

本发明属于煤矿开采覆岩运动地表变形预计领域，特别涉及一种适用于计算覆岩沉陷三维演变规律的煤层开采覆岩沉陷三维半解析分析方法。

背景技术

开采覆岩沉陷研究发展100多年来，各种理论分析方法复杂多样。其中，波兰学者Litwinszyn基于开采引起的岩层和地表移动规律与作为随机介质的非连续颗粒体介质模型所描述的规律在宏观上相似原理，于1954年提出的随机介质理论方法，由于其数理概念明确，实用性强，得到了广泛的应用，我国在50年代末由刘宝琛、廖国华等学者引入此方法，并逐渐发展成为目前在国内外应用最广泛，理论与实践上最为成熟的概率积分法。但由于该方法未能给出材料介质位移和应力响应方程，使得这种主要用于地表移动变形分析的方法对于开采覆岩沉陷问题的岩体内部应力、应变计算产生了困难。

半解析数值方法是指在数值分析中采用与引入部分解析解或解函数，而计算得到的仍是一系列离散化数值结果的分析方法。它通过应用解析手段来帮助减少数值计算的维数，即由人完成诸如一维解函数、基本解、通解等较为规则、简单情况下的理论分析，而同时又通过数值离散手段用计算机来处理降维后的低维离散化方程的求解，弥补人为所取解析函数的不足，以逼近真实解。基于上述原因，本发明所要解决的技术问题是，如何利用半解析数值分析方法来解决层状岩体矿山开采覆岩沉陷的三维分析。

发明内容

本发明的目的在于提出一种煤层开采覆岩沉陷三维半解析分析方法，该分析方法不仅能够克服纯解析的理论分析在数学上的困难及应用上的局限性，又能够大大降低基于全离散原理的纯数值方法的计算工作量。

为了实现上述目的，本发明采用如下技术方案：

一种煤层开采覆岩沉陷三维半解析分析方法，包括如下步骤：

s1、建立层状覆岩沉陷三维半解析模型，其中：

模型上部进行单元层分割，表土为水平层，视为横观各向同性材料，为水平条单元，采用二维解析的层单元，部分的局部坐标与整体坐标相重合；

模型中部进行三角块分割，呈三角状的倾斜块，离散为三棱柱单元；

模型下部进行单元层分割，包括开挖在内的呈倾斜状的各岩层，为倾斜条单元，采用二维解析的层单元，部分的局部坐标与整体坐标有一倾角，在局部坐标系下二者的计算公式有相同的形式；

s2、利用一些层、柱单元的标准解析解，并结合有限元的离散化分析，完成条单元、柱单元和层单元之间的转化，实现对层状覆岩沉陷三维半解析模型各部分的位移、单元刚度矩阵、单元载荷矩阵和应力的计算。

在上述步骤s2中，位移、单元刚度矩阵、单元载荷矩阵和应力的计算公式如下：

(1)使用模型层单元位移计算公式，如下：

{f}=uvw=Σm=1pΣn=1q[N]m|Nyz{δ}mn|---(1)]]>

其中，