数值分析方法在岩土重要性

① 数值分析方法能解决岩土工程中哪些问题

岩土的多相多场耦合。岩土的失稳与破坏。边坡的非侵入式有限元可靠度计算。土质堤防的渗流破坏。严寒区的路基的融沉计算。岩土与结构物的相互作用。

② 数值计算法

6.1.2.1 边坡数值计算的安全系数确定

数值分析方法考虑岩土体应力应变关系，克服了极限平衡方法的缺点，为边坡稳定分析提供了较深入的概念。

目前，数值计算的失稳判据主要有两类:一是以数值计算不收敛作为失稳的标志;二是以广义塑性应变或者等效塑性应变从坡脚到坡顶贯通作为边坡破坏的标志。而用数值分析结果获取边坡安全系数也主要有两种方法:强度折减法、数值计算与极限平衡的耦合分析法。

(1)强度折减法:首先选取初始折减系数，将岩土体强度参数进行折减，将折减后的参数输入，进行数值计算，若程序收敛，则岩土体仍处于稳定状态，然后需要再增加折减系数，直到程序恰好不收敛，此时的折减系数即为稳定或安全系数。_[52]

(2)数值计算与极限平衡的耦合分析法:首先采用数值分析法，计算边坡内的应力应变以及位移分布;然后将计算的应力分布结果，通过应力张量变换，求出指定滑动面上的应力分布;最后通过极限平衡方法求出与该滑动面对应的稳定性安全系数。_[52]

6.1.2.2 边坡数值计算方法存在的问题剖析

应该指出，尽管近年来数值模拟方法和理论方面取得了显着的进展，但仍不能很好的适应岩土工程的复杂情况，其主要原因有两方面:(1)数学模型的不确定性。由于岩体力学性质千变万化(弹性、塑性、流变、应变硬化及应变软化等)，且具有复杂的结构特性(岩体结构、岩体介质结构及地质结构等)，不但至今对岩体的失稳或破坏还缺少可靠的判据或准则，而且工程开挖方法、开挖步序对围岩的力学状态(应力和应变)及稳定条件具有重大的影响，在某些情况下还起到决定性的作用，这使得目前对于数学模型的建立，尤其是本构模型的给定还带有相当程度的盲目性。(2)参数的不确定性。岩体的物理力学性质、初始地应力等参数多变，仅通过有限的现场调查和室内试验来获得参数输入信息，数据往往具有很大的离散性，很难全面反映岩体真实情况。

“数学模型给不准”和“输入参数给不准”的困难已成为岩体力学数值分析应用的“瓶颈”问题。事实上，无论数值分析技术多么发达，它们总只是某种手段，关键还是对岩体基本特性的认识。

③ 数值计算方法概述

在采矿工程中，数值模拟方法不仅能模拟岩体复杂的力学和结构特征，还能很方便地解决现场监测过程中需要大量人力、物力而无法完成的、现有力学理论不能求解的复杂形体问题，并对矿山岩体稳定性进行预测与预报。

关于岩土工程的数值分析方法，很多学者都作过系统综述^{[53，68，72]}，笔者只拟简单介绍。岩土工程数值分析方法，主要分为三大类，如图7-1所示。

图7-1 边坡工程数值分析方法

(1)连续介质数值分析方法

连续介质数值分析方法的理论基础是弹(塑)性力学。因此，在该类数值分析方法公式的推导过程中，需要满足基本方程和边界条件。只是在求解手段上，采用了不同于弹性力学的各种近似解法。这类数值分析方法包括有限差分法、有限单元法和边界单元法等，它适用于连续介质体的地下工程围岩与结构的应力分析和位移求解。

(2)非连续介质数值分析方法

非连续介质数值分析方法的理论基础是牛顿运动定律，它并不满足结构的位移连续条件，但是可以求出结构在平衡状态下的位移或者在不可能处于平衡状态时的破坏模式。此外，尽管结构不受位移连续的约束，但应满足给定的单元和交界面的本构定律。这类数值分析方法主要有离散单元法和不连续变形分析(DDA)。这些数值分析方法可用于分析节理岩体可能发生的不连续变形，如洞室围岩附近岩块的分离与滑落等。

(3)混合介质数值分析方法

混合介质数值分析方法是连续和不连续分析方法的耦合。在地下结构的某些区域(如洞室附近)，围岩体由于开挖影响而发生块体的分离而不连续，在另外区域(如远离洞室)，则岩体一般仍相互联系而处于连续状态。因此，考虑两种不同力学介质的耦合分析很必要。目前常见的耦合方法有有限元与离散元的耦合、边界元与离散元的耦合等。混合介质吸取连续介质和非连续介质两种数值分析方法中的优点，在可能发生不连续变形的岩体，采用非连续介质方法模拟，而远离洞室的岩体一般仍处于连续状态，可采用连续介质模型分析。

本章分别采用有限元强度折减法、有限元和离散元相结合的CDEM法、FLAC差分法，开展安家岭露天矿露天井工联合开采的数值模拟分析，研究露天开采和井工开采的相互作用及影响规律。

④ 岩土数值软件的目前的现状和发展是怎样的

岩土数值软件的目前的现状和发展还可以。

岩土工程师主要研究岩土构成物质的工程特性。其首先研究从工地采集的岩土样本以及岩土样本中的数据，然后计算出工地上的建筑所需的格构。地基、桩、挡土墙、水坝、隧道等的设计都需要岩土工程师为其提供建议。

岩土工程师的发展可以分为两大块：一块在咨询公司，这是未来岩土工程师主要的服务企业。其主要负责勘察、测试、设计、检验、监测等与数据、论证、决策有关的工作。另一块在工程公司，其主要负责岩土工程的实施。

注意：

数值岩土力学是解决岩土工程问题的较有效的手段, 已被学术界和工程界广泛接受, 作为一种力学状态的分析工具, 它越来越多的应用于岩土体的稳定分析、岩土工程设计和岩土工程基本问题分析中。

随着岩土工程领域的不断扩展与延伸, 岩土力学数值分析方法得到了迅速发展, 出现了各种各样的数值分析方法。将岩土力学常用的数值分析方法分为连续变形分析方法与非连续变形分析方法两大类。

⑤ 除有限单元法外，岩土工程常用到哪些数值方法，并对比其优缺点

岩土工程常用的数值方法包括：有限差分法、边界元法、离散元法、颗粒元法、不连续变形分析法、流形元法、模糊数学方法、概率论与可靠度分析方法、灰色系统理论、人工智能与专家系统、神经网络方法、时间序列分析法。
有限单元法的优缺点：有限单元法的理论基础是虚功原理和基于最小势能的变分原理，它将研究域离散化，对位移场和应力场的连续性进行物理近似。有限单元法适用性广泛，从理论上讲对任何问题都适用，但计算速度相对较慢。即，物理概念清晰、灵活、通用、计算速度叫慢。
有限差分法：该方法适合求解非线性大变形问题，在岩土力学计算中有广泛的应用。有限差分法和有限单元法都产生一组待解方程组。尽管这些方程是通过不同方式推导出来的，但两者产生的方程是一致。另外，有限单元程序通常要将单元矩阵组合成大型整体刚度矩阵，而有限差分则无需如此，因为它相对高效地在每个计算步重新生成有限差分方程。在有限单元法中，常采用隐式、矩阵解算方法，而有限差分法则通常采用“显式”、时间递步法解算代数方程。
边界元法：该方法的理论基础是Betti功互等定理和Kelvin基本解，它只要离散求解域的边界，因而得到离散代数方程组中的未知量也只是边界上的量。边界元法化微分方程为边界积分方程，离散划分少，可以考虑远场应力，有降低维数的优点，可以用较少的内存解决较大的问题，便于提高计算速度。
离散元法：离散元法的理论基础是牛顿第二定律并结合不同的本构关系，适用对非连续体如岩体问题求解。该方法利用岩体的断裂面进行网格划分，每个单元就是被断裂面切割的岩块，视岩块的运动主要受控于岩体节理系统。它采用显式求解的方法，按照块体运动、弱面产生变形，变形是接触区的滑动和转动，由牛顿定律、运动学方程求解，无需形成大型矩阵而直接按时步迭代求解，在求解过程中允许块体间开裂、错动，并可以脱离母体而下落。离散元法对破碎岩石工程，动态和准动态问题能给出较好解答。
颗粒元法：颗粒元方法是通过离散单元方法来模拟圆形颗粒介质的运动及其相互作用，它采用数值方法将物体分为有代表性的多个颗粒单元，通过颗粒间的相互作用来表达整个宏观物体的应力响应，从而利用局部的模拟结果来计算颗粒群群体的运动与应力场特征。 不连续变形分析方法：该方法是并行于有限单元法的一种方法，其不同之处是可以计算不连续面的错位、滑移、开裂和旋转等大位移的静力和动力问题。此方法在岩石力学中的应用备受关注。
流形元法；该方法是运用现代数学“流形”的有限覆盖技术所建立起来的一种新的数值方法。有限覆盖是由物理覆盖和数学覆盖所组成的，它可以处理连续和非连续的问题，在统一解决有限单元法、不连续变形分析法和其他数值方法的耦合计算方面，有重要的应用前景。
无单元法：该方法是一种不划分单元的数值计算方法，它采用滑动最小二乘法所产生的光滑函数去近似场函数，而且又保留了有限单元法的一些特点。它只要求结点处的信息，而不需要也没有单元的信息。无单元法可以求解具有复杂边界条件的边值问题，如开裂问题，只要加密离散点就可以跟踪裂缝的传播。它在解决岩石力学非线性、非连续问题等方面具有重要价值和发展前景。
混合法：对于复杂工程问题，可采用混合法，即有限单元法、边界元法、离散元法等两两耦合来求解。
模糊数学方法：模糊理论用隶属函数代替确定论中的特征函数描述边界不清的过渡性问题，模糊模式识别和综合评判理论对多因素问题分析适用。 概率论与可靠度分析方法：运用概率论方法分析事件发生的概率，进行安全和可靠度评价。对岩土力学而言，包括岩石（土）的稳定性判断、强度预测预报、工程可靠度分析、顶板稳定性分析、地震研究、基础工程稳定性研究等。
灰色系统理论：以“灰色、灰关系、灰数”为特征，研究介于“黑色”和“白色”之间事件的特征，在社会科学及自然科学领域应用广泛。岩土力学中，用灰色系统理论进行岩体分类、滑坡发生时间预测、岩爆分析与预测、基础工程稳定性、工程结构分析，用灰色关联度分析岩土体稳定性因素主次关系等。
人工智能与专家系统：应用专家的知识进行知识处理、知识运用、搜索、不确定性推理分析复杂问题并给出合理的建议和决策。岩石力学中，可进行如岩土（石）分类、稳定性分析、支护设计、加固方案优化等研究。 神经网络方法：试图模拟人脑神经系统的组织方式来构成新型的信息处理系统，通过神经网络的学习、记忆和推理过程进行信息处理。岩石力学中，用于各种岩土力学参数分析、地应力处理、地压预测、岩土分类、稳定性评价与预测等。
时间序列分析法：通过对系统行为的涨落规律统计，用时间序列函数研究系统的动态力学行为。岩石力学中，用于矿压显现规律研究、岩石蠕变、岩石工程的位移、边坡和硐室稳定性等、基础工程中降水、开挖、沉降变形等与时间相关的问题。

⑥ 应力场数值模拟方法

近30年来，人们采用现场测试、实验室试验、理论分析与模型试验等多种方法，使岩土力学研究取得很大进展^{［162～166］}。如今随着计算机技术的快速发展，岩土力学的研究进入了一个新的阶段，其中数值计算方法已成为解决岩土力学问题的重要手段之一。

6.1.1 概述

许多工程分析问题，如固体力学中的位移场和应力场分布分析、电磁学中的电磁场分析、振动特性分析、传热学中的温度场分析以及流体力学中的流场分布等，都可以通过在给定边界条件下对其控制方程进行求解得到，但是利用解析方法只能求出一些方程性质比较简单且几何边界相当规则的极少数问题。对于大多数实际工程技术问题，由于物体的几何形状比较复杂或者问题的某些特性是非线性的，因而一般无解析解。为了解决此类问题，一般采用两种处理方法：一种是进行简化处理，将方程和边界条件简化为能够处理的问题，从而得到在简化情况下的解，但这种方法应用非常有限，且假设过多将会导致错误的解；另一种是在广泛接收现代数学和力学理论的基础上，借助于计算机和计算软件来获得工程上要求的数值解，这就是目前应用非常广泛的数值模拟方法。

目前在工程技术领域内常用的数值分析方法包括：有限单元法、边界元法、离散单元法以及有限差分法。最初常用的是有限差分法，它可以处理一些相当复杂的问题。但对于几何形状复杂的边界条件，其解的精度受到影响。20世纪60年代出现并得到广泛应用的有限单元法，使经典力学解析方法难以解决的工程力学问题都可以用有限元方法求解。它将连续的求解域离散为一组有限个单元的组合体，解析地模拟或逼近求解区域。由于单元能按各种不同的联结方式组合在一起，且单元本身又可有不同的几何形状，所以能适应几何形状复杂的求解域。但有限单元法需要的存贮容量常非常巨大，甚至大得无法计算。由于相邻界面上只能位移协调，对于奇异性问题（应力出现间断）的处理比较麻烦，这是有限单元法的不足之处。70年代末期，出现了另一种重要的数值方法为边界元法。边界元方法是把求解区域的边界剖分为若干个单元，将求解简化为求单元结点上的函数值，通过求解一组线性代数方程实现求解积分方程。上述两种数值方法的主要区别在于，边界元法是“边界”方法，而有限元法是“区域”方法，它们都是针对连续介质，只能获得某一荷载或边界条件下的稳定解。对于具有明显塑性应变软化特性和剪切膨胀特性的岩体，无法对其大变形过程中所表现出来的几何非线性和物理非线性进行模拟，这就使得人们去寻求适合模拟节理岩体运动变形特性的有效数值方法。

1971年Cundall，P.A^［167］提出了一种不连续介质数值分析模型——离散单元法。该方法优点在于适用于模拟节理系统或离散颗粒组合体在准静态或动态条件下的变形过程。离散单元法的基本原理不同于基于最小总势能变分原理的有限单元法，也不同于基于Betti互等定理的边界单元法，而是建立在牛顿第二运动定律基础上。最初的离散元法是基于刚性体的假设，由于没有考虑岩块自身的变形，在模拟高应力状态或软弱、破碎岩体时，不能反映岩块自身变形的特征，使计算结果与实际情况产生较大出入。Maini，T.，Cundall，P.A.^{［168～169］}等人针对刚体单元没有考虑岩块自身变形的缺点，利用差分方法提出了考虑岩石自身变形的改进的离散单元法，编制了通用的离散元程序UDEC（Universal Discrete Element Code），将离散元推广到模拟岩体破碎和变形情况，推动了离散元的进一步发展。我国学者也相继开展这方面的研究，王泳嘉教授^［170］等将离散单元法应用于采矿工程方面的研究。

6.1.2 FLAC数值模拟方法

（1）概述

数值模拟技术通过计算机程序在工程中得到广泛的应用。一直到20世纪80年代初期，国际上较大型的面向工程的通用程序有：ANSYS、NASTRAN、FLAC、UNDEC、ASKS以及ADINA等程序。它们功能越来越完善，不仅包含多种条件下的有限元分析程序，而且带有功能强大的前、后处理程序。

连续介质快速拉格朗日差分法（Fast Lagrangian Analysis of Continua，简写FLAC）是近年来逐步成熟完善起来的一种新型数值分析方法。把拉格朗日法移植到固体力学中，即将所研究的区域划分为网格，节点相当于流体质点，然后按照时步用拉格朗日方法来研究网格节点的运动，这就是固体力学变形研究中的拉格朗日数值研究方法。

FLAC与基本离散元法相似，但它克服了离散元法的缺陷，吸取了有限元法适用于各种材料模型及边界条件的非规则区域连续问题解的优点。FLAC所采用的动态松弛法求解，不需要形成耗机时量较大的整体刚度矩阵，占用计算机内存少，利于在微机的工程问题。同时，FLAC还应用了节点位移连续的条件，可以对连续介质进行大变形分析。

（2）数学模型

显式有限差分法的基本方程主要包括：平衡方程、几何方程、物理方程和边界条件。在FLAC^3D2.0中采用的拉格朗日描述方程，一般规定介质中一点由向量分量x_i，u_i，v_i，dv_i/dt（i=1，2，3）来表征，其分别代表位置、位移、速度和加速度分量。

其基本原理和基本公式简单叙述如下：

空间导数的有限差分近似

三维FLAC方法中采用了混合离散方法，区域被划分为常应变六面体单元的集合体；而在计算过程中，又将每个六面体分为常应变四面体，变量均在四面体上进行计算，六面体单元的应力、应变取值为其四面体的体积加权平均。

如图6.1所示，所研究区域任一四面体，节点编号为1～4，规定与节点n相对的面为第n面，设定其内任一点的速度分量为v_i，则由高斯散度定理得

煤岩动力灾害力电耦合

式中：V——四面体体积，m³；S——四面体外表面，m²；n_j——外表面单位法向向量分量。

图6.1 四面体

对于常应变单元，n_j在每个面上为常量，因此通过上式积分可得

煤岩动力灾害力电耦合

式中上标f表示f面的变量值，对于为线性分布的速率分量，速度分量的平均值为

煤岩动力灾害力电耦合

式中上标l表示节点l的变量值。将（6.3）式代入（6.2）式可得

煤岩动力灾害力电耦合

经过变换可得节点速率计算公式：

煤岩动力灾害力电耦合

1）平衡方程（运动方程）

显式有限差分法采用的平衡方程就是人们熟知的牛顿第二运动定律，即

煤岩动力灾害力电耦合

式中：F_i——节点合力在i方向分力，N；m_i——节点质量，kg；a_i——节点加速度在i方向分量，m/s²。

作用于各个节点的合力：外力（集中力、均布力、重力等）和内力（单元变形引起的应力在单元节点上的分量）。节点质量是根据节点相邻单元的面积（体积）和密度，按照面积（体积）加权求出。

FLAC^3D以节点为计算对象，将力和质量均集中在节点上，然后通过运动方程在时域内进行求解。节点运动方程可以表示为如下形式：

煤岩动力灾害力电耦合

式中：（t）———t时刻l节点在i方向的不平衡力分量，可以由虚功原理导出；m^l———l节点的集中质量，在分析静态问题时，采用虚拟质量；而在分析动态问题时，则采用实际的集中质量。

将（6.7）式左端用中心差分来近似，则可得

煤岩动力灾害力电耦合

2）变形协调方程——几何方程

作为连续介质力学，变形体之间必须满足变形协调方程（几何方程），否则变形体就会出现分离或嵌入。变形协调方程反映了位移与应变间的关系，对于某一时步的单元应变增量可由下式确定：

煤岩动力灾害力电耦合

求出应变增量后，即可由本构方程得到应力增量，各时步的应力增量叠加即可得到总应力，在大变形时，还需根据本时步单元的转角对本时步前的总应力进行旋转修正，然后即可由虚功原理求出下一时步的节点不平衡力，进入下一时步的计算。

3）物理方程——本构关系

物理方程反映应力与应变之间的关系，在程序中通常被称为材料模式或材料模型。在FLAC^3D2.0中提供了10种基本材料模型，它们是：①Null；②Elastic，isotropic；③Elastic，transversely isotropic；④Druck-Prager plasticity；⑤Mohr-Coulomb plasticity；⑥Ubiquitous joint plasticity；⑦Strain-hardening／softening Mohr-Coulomb plasticity；⑧bilinear strain-hardening/softening ubiquitous-joint plasticity；⑨Modified Cam-clay plasticity 和⑩elastic，orthotropic。

本文进行应力场数值模拟时采用的是Mohr-Coulomb应变硬化软化破坏准则，在FLAC^3D2.0中，Mohr-Coulomb 模型的破坏准则以主应力σ₁，σ₂，σ₃来描述，相应的应变为三个主应变ε₁，ε₂，ε₃。根据Hooke定律，应力、应变增量具有如下表达形式：

煤岩动力灾害力电耦合

式中α₁，α₂为材料常数，可以由体积模量K和剪切模量G确定：

煤岩动力灾害力电耦合

不失一般性，令σ₁≥σ₂≥σ₃，摩尔—库仑准则为

其中：

煤岩动力灾害力电耦合

式中C，φ分别为煤岩的粘聚力和内摩擦角。

FLAC^3D2.0的Mohr-Coulomb 破坏准则如图6.2所示。

图6.2 FLAC^3D的Mohr-Coulomb 破坏准则

本着作中就是选用上述的Strain-hardening／softening Mohr-Coulomb plasticity模型，对单轴压缩煤岩以及矿山地下煤岩独巷掘进时围岩的变形破坏过程进行模拟。

4）阻尼力

对于静态问题，FLAC^3D2.0在式（6.7）的不平衡力中加入了非黏性阻尼，以使系统的振动逐渐衰减直至达到平衡状态（即不平衡力接近零），此时节点运动方程变为：

煤岩动力灾害力电耦合

式中阻尼力（t）由下式确定：

煤岩动力灾害力电耦合

上式中α为阻尼系数，其默认值为0.8；而：

煤岩动力灾害力电耦合

5）初始条件与边界条件

边界条件包括面积力、集中载荷等应力边界条件和位移边界条件。此外也可加载体力和初始应力。在编写程序代码时，一般所有的应力和节点速度初始化为零，然后指定初始化应力。集中载荷则加载在面节点上，位移边界条件则以运动方程形式施加到相应的边界节点上。

边界条件分为应力边界条件和位移边界条件，应力边界条件为：

煤岩动力灾害力电耦合

式中：F_i———作用于节点i上的力；——作用于边界上的应力；n_j———边界上的法线沿j方向的矢量大小；Δs———边界的长度。

若是位移边界条件，应将边界条件以运动方程的形式施加到相应的边界节点上。

FLAC^3D2.0^［171］与FLAC^2D3.3也是由美国Itasca Consulting Group Inc开发的三维显式有限差分法程序，它可以模拟岩土或其他材料的三维力学行为。FLAC^3D2.0的计算循环过程如图6.3所示。

图6.3 FLAC^3D2.0的计算循环

6.1.3 FLAC数值模拟方法在采矿工程中的应用^{［172～179］}

采矿过程中围岩活动规律及巷道围岩稳定性问题涉及岩体力学特性、围岩压力、支护围岩相互作用关系及巷道与工作面时空关系等一系列复杂力学问题。随着我国经济建设的高速发展，岩土工程稳定性分析问题日益突出，除采矿工程外，在水利、交通（铁道和公路）、高层建筑的地基等行业也都存在着大量的岩土力学数值计算分析问题。能否用计算机数值模拟分析采矿岩层控制问题和岩土工程问题已成为一个大学岩层控制技术和岩土力学学科水平高低的标志之一。

与ANSYS、ADINA相比，FLAC 和UDEC的最大特点是计算分析岩土工程中的物理不稳定问题，因而特别适用于岩土工程中几何和物理高度非线性问题的稳定性分析，如采场的采动影响规律，软岩巷道的大变形问题，采动后的地表沉陷，露天矿的边坡稳定，水坝的稳定性等问题。

从力学计算方法上讲其主要特点

1）可以直接计算非线性本构关系；

2）物理上的不稳定问题不会引起数值计算的不稳定；

3）开放式程序设计（FISH），用户可以根据需要自己设计程序；

4）既可以分析连续体问题（FLAC），也可以分析非连续体问题（UDEC）；

5）可以模拟分析很大的工程问题；

6）高度非线性问题不增加计算时间。

在采矿工程中，许多学者利用FLAC软件对采矿过程中围岩活动规律及巷道围岩稳定性问题涉及到岩体力学特性、围岩压力、支护围岩相互作用关系及巷道与工作面的时空关系等一系列复杂的力学问题进行了一系列的研究，取得了显着的效果。梅松华等以施工期监测结果为基础，在正交设计原理的基础上，选定反演参数与水平，采用二维显式差分法FLAC进行弹塑性位移反分析。朱建明等在分析FLAC有限差分程序的基础上，提出了变弹性模量方法模拟时间因素对巷道围岩稳定性影响的衰减曲线，为揭示巷道围岩变形机理和有效指导围岩支护提供了有效的分析方法。来兴平等探讨了岩石力学非线性计算软件FLAC^2D3.3在地下巷道离层破坏数值计算中的应用。康红普对回采巷道锚杆支护影响因素进行了FLAC分析，认为FLAC^2D3.3在分析几何非线性和大变形问题方面性能优越。

在煤岩动力灾害预测中，这些方法的优点

1）可以提前知道煤与瓦斯突出、冲击矿压等煤岩动力灾害防治的重点区域；

2）可以得到大范围内的空间信息；

3）可以提前预测预报煤岩动力灾害的危险性；

4）可以确定在采掘过程中，应力的分布状况和集中程度。

在煤岩动力灾害预测中，这些方法也具有以下缺点

1）对实际问题均进行了简化处理；

2）对于煤岩体的力学特性，如弹性模量、泊松比等力学参数，也进行了简化，没有考虑其局部非均质性和各向异性；

3）只能作为一种近似方法使用。