应用FLAC3D研究地下开采对围岩影响范围.docVIP

应用FLAC3D研究地下开采对围岩的影响范围 摘要：利用FLAC3D软件研究地下开采对围岩的影响范围。结果表明，地下工程的开挖对围岩造成的影响。在空区的上部、下部及四周．影响范围均大致为3倍采空区尺寸，远大于该尺寸即可近似视为无限边界。 关键词：采矿工程；顶板；FLAC3D；影响范匿；地下开采 由于地下工程的结构和岩石破坏过程十分复杂，传统的定量理论分析方法假设较多，相似材料模拟成本高且误差大，影响研究结果的真实性。而数值计算由于其可以综合考虑多种因素，且计算结果直观，在地下工程稳定性研究中得到了日益广泛的应用，已经成为岩石力学的一种重要研究手段，也是分析水平承载层(如顶板)稳定性的有力工具。运用FLAC3D软件对地下开采建模及计算过程分析中产生的若干问题进行研究。 1 FLAC3D程序简介 FLAC3D是Fast Lagrangian Analysisof Continua in 3 Dimensions的简写，是美国明尼苏达ITASCA软件公司编制开发的三维显式有限差分程序。FLAC3D程序自美国ITASCA咨询公司推出后，已成为目前岩土力学计算中的重要数值方法之一。FLAC3D是二维应用程序FLAC2D在三维空间的扩展，用于模拟三维土体、岩体或其他材料的力学特性，尤其是达到屈服极限时的塑性流变特性，广泛应用于边坡稳定性评价、支护设计及评价、地下洞室、 施工设计(开挖、填筑等)、河谷演化进程再现、拱坝稳定性分析、隧道工程、矿山工程等多个领域。FLAC3D采用拉格朗日差分公式来处理有限变形问题，将计算域划分为若干单元，单元网格可以随着材料的变形而变形，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析等领域有独到的优点，是一种功能强大的工程计算软件。FLAC3D程序于20世纪90年代被引入我国，主要应用于岩土力学分析，例如边坡稳定性、开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道围岩稳定性研究等，目前FLAC／FLAC3D软件已广泛应用于岩土工程、采矿工程、水利工程、地质工程等行业【1—-3】。 FLAC3D的基本原理和算法与离散元相似，在求解过程中采用了离散元的动态松弛法，不需要求解大型联立方程组(刚度矩阵)。FLAC3D采用最大不平衡力来描述计算的收敛过程，如果单元的最大不平稳力随着时步增加而逐渐趋近极小值，则计算是稳定的，否则计算不稳定。 在FLAC3D所采用的显式时间差分求解中，所有的矢量参数(力、速度、位移)都储存在网格节点上，所有的标量及张量(应力及材料特性)储存在单元的中心位置。通过运动方程由应力及外力可以求出节点的速度和位移，再由空间导数得出单元应变率，借助于材料的本构关系，由单元应变率可以获得新的单元应力。计算循环如图l所示。 FLAC3D具有强大的前后处理功能，其自动三维网格生成器定义了多种基本的单元形状，通过组合基本单元和利用内嵌的FISH语言自定义单元形态、材料分布、特殊的边界条件等，可以生成复杂的三维网格，如交叉隧道等。在计算过程中，用户可以随时用彩色或灰度图或数据文件输出结果，图形上可以表示出网格、结构以及有关变量的等值线、矢量图等，也可以给出任意截面上的变量图或等值线图，计算域可旋转以从不同的角度观测结果。 FLAC3D计算结果的可靠程度与岩体力学参数、原岩应力、岩体结构类型、计算域的结构效应及其开挖、加载等约束条件有关。 FLAC3D软件的不足之处制约其进一步推广应用【3】，主要表现在：(1)对网格尺寸十分敏感，同一模型采用不同尺寸的网格单元可能导致求解时间相差数倍；(2)对于与时间密切相关的物理过程(如固结、长期动力影响、材料流变等)，求解时间很长；(3)前处理功能较弱，尤其建立较复杂的地质模型时，费时费力，计算周期长、难度大。 2 建模过程影响因素分析 顶板岩层稳定性的影响因素很多，但以其强度、构造发育程度、分层厚度、埋藏深度、泊松比及弹性模量为主要因素，且前三者可视为决定性因素。随着开采深度的增加，深度对顶板稳定性的影响会愈加明显。岩层中的裂隙发育程度与其强度密切相关，因此，也与顶板稳定性密切相关。国内顶板分类方案中普遍考虑了岩层的岩石结构、岩石的物理力学性质、结构面的形态和发育程度、覆盖层厚度、含水情况等，此外还有时效特性和采矿工程因素。这 些因素以各自的作用机理对岩层顶板稳定性发挥作用。在采掘过程中应加强对开采地质条件的预测，以便采取有效措施控制顶板。 高峰【4】将顶板冒落的主控因素归纳为9个方面：地层与岩性、地质与构造、地应力、地下水、工程环境、采空区结构特征、开采方法、支护与空区处理、时间因素。其多方案数值模拟计算显示，顶板下开采矿体厚度、泊松比与顶板安全厚度成一次函数关系，岩石弹性模量、采矿推进距离、矿柱宽度、充填体厚度、内聚力与顶板安