基于FLAC3D数值模拟确定合理封孔深度

基于FLAC3D的分段空场嗣后充填法采场结构参数优化发布时间：2022-05-12T05:37:26.971Z 来源：《科学与技术》2022年第3期作者：罗代祺[导读] 为实现矿山安全高效生产，确定合理采场结构参数成为矿山生产首要任务罗代祺云南铜业股份有限公司，云南昆明650000摘要：为实现矿山安全高效生产，确定合理采场结构参数成为矿山生产首要任务。

将经验法与数值模拟法相结合，基于Mathews稳定图法初步稳定性分析以及FLAC3D数值模拟6种方案，选取顶板位移、最大主应力作为评价因子，确定最优采场结构参数。

研究结果表明，落凼矿区最佳分段空场嗣后充填法采场结构参数为采矿跨度15m，采场长度32m，间柱宽度5m；暴露面积越小，顶板位移和最大主应力越小；暴露面积相等情况下，采场跨度越小，顶板位移和最大主应力越小。

关键字：分段空场嗣后充填法；采场；结构参数；数值模拟0引言在矿山生产中，采场结构参数与采场稳定性息息相关，直接影响采矿工作人员生命安全，也关乎矿山生产能力和企业效益。

数值模拟法确定最优采购结构参数具备多方案、可重复、可视化等优点，徐帅[1]基于FLAC3D进行数值计算，建立SOM模型对深埋厚大矿体采场结构参数进行优化研究。

郝益民[2]运用FLAC3D对4种方案进行数值模拟研究，确定了有利于阶段空场嗣后充填采场的结构参数。

陈晖[3]采用FLAC3D开展“三软”条件下金矿采场结构参数优化。

江飞飞[4]基于三维离散元程序3DEC对不同的一、二步骤采场宽度组合方案进行了模拟，确定了合理的采场结构参数。

龙科明[5]运用ANSYS对不同采场结构参数进行了数值模拟分析，确定了采场最优宽高组合。

因此，结合经验类比法，依靠数值模拟来确定采场结构参数，能够为拉拉铜矿采矿稳定性分析提供可靠依据。

1工程概况拉拉铜矿落凼矿区深部矿段矿体按倾角分为三类，＜15°的水平至缓倾斜矿体、15°~35°的缓倾斜至倾斜矿体以及＞35°的倾斜至急倾斜矿体。

基于FLAC（3D）的深部巷道围岩稳定性数值模拟研究贾晓亮【摘要】针对由于巷道剧烈底鼓、巷道断面变形造成巷道运输及维护困难等问题,采用数值模拟计算软件FLAC~（3D）,数值计算得到宏发煤矿＋1 650 m集中运输大巷巷道开挖后围岩垂直应力、水平应力分布及塑性区分布,分析并得到了巷道底鼓机理：距离巷道底板临空面更近的倾斜C_（17）煤层的软弱顶板岩层随掘进扰动而首先破坏,造成左侧局部底鼓,继而巷道底板围岩应力因底鼓而持续动态重新分布,并逐渐向右侧及深部发展,致使底鼓逐渐恶化并表现出左右非对称特征。

综合考虑施工、经济、效果、瓦斯等因素,提出采用＂构筑反底拱＋底板注浆＋底板锚杆＂联合支护的底鼓控制措施。

【期刊名称】《能源与环保》【年(卷),期】2017(039)006【总页数】5页(P18-22)【关键词】FLAC3D 底鼓数值模拟围岩巷道支护【作者】贾晓亮【作者单位】中煤炭科工集团重庆研究院有限公司,重庆400037【正文语种】中文【中图分类】TD862.1深部巷道围岩变形主要表现为围岩变形量大、变形速度快、巷道持续变形、流变、深部岩石扩容等特征[1-3]。

针对该情况，国内学者进行相关研究，张合超、房健、苏军等[4-6]采用FLAC3D模拟了不同支护参数条件下巷道支护效果，研究了地下开采对围岩的影响范围。

徐文彬、许梦国、余伟健等[7-9]分析了交叉型巷道群围岩由于不同的开挖顺序而产生不同的应力状态、位移变形以及塑性屈服变化特征，对比了不同布置形式的底部结构在开挖影响下的稳定性。

宁建国、刘泉声、何富连等[10-12]采用FLAC3D数值模拟软件分析空间交叉巷道的应力分布特点以及沿下方巷道轴向方向形成的应力增高和降低区。

牛学超、闫长斌、肖明等[13-15]提出了大型地下洞室施工开挖过程动态模拟的三维有限元数值分析方法，并通过工程实例分析爆破震动作用下地下硐室群的稳定性。

宏发煤矿地处云贵高原，地质构造极为复杂，巷道所处的煤层直接顶为粉砂质泥岩，顶板节理发育，表现出深部巷道围岩变形的特征，且其巷道原支护结构方式和参数不合理等原因使得巷道发生严重变形和破坏，其中深部高应力和围岩较低承载能力是该矿巷道变形破坏的主要原因。

基于FLAC3D的高速公路隧道两种开挖方式稳定性分析采用数值计算对高速公路隧道施工过程进行分析，得到了施工过程中围岩的变形和支护结构的受力状态。

结果表明，随着开挖的推进，中隔墙上部以及洞侧壁处围压变形和应力较为集中，易产生相关地质灾害。

所以，在相关工程施工过程中应该对以上部位进行高质量的支护防灾。

对于底部在中墙位置隆起较为严重，所以应该加强此处的加固措施，条件允许是可开挖较深施做浅基础。

标签：隧道；围岩；数值计算1 概述近年来，随着经济建设及基础设施的快速发展，我国对于高速公路的需求越为强烈。

为更好的选线施工，减少施工期间地质灾害的威胁，方便施工以及满足交通需求，必须开山造路或者凿洞通路，这也就面临隧道开挖相关复杂地质问题。

目前，山岭隧道越来越多，比如川兰铁路的山岭隧道占到了线路总长的一半，而山岭隧道所遇到的围岩复杂多变，构造运动影响强烈，隧道断面的跨度也越来越大，所以不同的断面形式在一定方面适应复杂地区隧道的开挖。

随着开挖期间的多期次开挖扰动使得岩体强度明显减弱，岩体破碎严重，其自身稳定性明显降低，同时也造成后期施工难度加大，软弱夹层地带更为明显，易产生类型及机制较为复杂的地质灾害。

本文在基于FLAC数值软件的计算下，主要对Ⅵ级围岩下的大跨度隧道在不同的开挖方式下其对应围岩稳定性进行对比分析研究。

2 计算方案2.1 计算模型某隧道长约270m，围岩为Ⅴ、Ⅵ级，开挖跨径最大为37.0m，埋深最大为75m。

研究区地貌主要为侵蚀山地和剥蚀山地，地层岩性主要为人工填土、碎石、粉砂泥质岩和泥岩组成，隧道入口所出露岩性为残破积亚粘土夹碎石。

整个隧道橫穿山岭地带，且隧道入口比较陡峭，坡度角约为50°，隧道出口较入口比较缓，坡角约30°。

以隧道入口所出露岩性为研究对象，运用FLAC3D软件建立三维地质及结构模型（见图1）。

由于数值模型的建立必然会受到边界效应的影响，所以在所建模型之时，将模型埋深取为48.5m，隧道左右内壁距模型边缘均为50m，沿隧道垂直走向方向长度为120m，隧道内壁顶部距底部边界为40m，沿隧道走向方向长度为10m。

本科生毕业论文（设计）题目：基于flac3D深基坑开挖模拟与支护设计指导教师: 职称:评阅人: 职称:摘要随着城市化过程中不断涌现的高层建筑和超高层建筑以及城市地下空间的开发，深基坑工程越来越多，深基坑工程项目的规模和复杂性日益增大，给深基坑工程的设计和施工带来了更大的挑战。

在这样的背景下，深基坑支护结构设计和变形量预测已成为岩土工程领域的重要研究课题之一。

本文以武汉市万达广场深基坑工程作为研究对象，利用勘查资料和深基坑支护结构设计要求，比选合理的基坑支护方案并进行相应的计算设计。

同时，本文针对深基坑工程变形量验算等难以解决的问题引用了flac3D数值模拟方法，对基坑开挖、支护结构施工进行全方位的模拟监测，将计算设计结果和模拟计算结果进行对比验算，得出比较合理的支护结构设计方案和变形量控制方案。

根据基坑实际情况和勘查资料，本文选择的围护方案为以大直径混凝土排桩、双排桩、角撑与对顶撑相结合的内支撑为主的多种联合支护方案，结合坡顶大面积卸土减载、坑内被动区加固的措施。

计算部分主要设计计算大直径混凝土排桩（钻孔灌注桩）桩长、内力和配筋，而对卸土减载、内支撑结构、坑内被动区加固和降水设计只进行了简要的说明；flac3D模拟部分主要从建立模型、设置大直径混凝土排桩、放坡开挖、放坡坡面土钉施工、预应力锚索（代替内支撑）施工和基坑主体开挖为顺序进行建模计算，最后进行变形量监测、分析，输出桩单元、锚单元的内力分布情况并给出相应的结论与建议。

本文以常规计算和数值模拟相结合的方式进行参考对比，常规计算和数值模拟分析结果非常接近，给出了有效合理的安全系数。

关键词：深基坑支护设计flac3D模拟数值模拟AbstractWith the urbanization process ,high-rise buildings and supertall buildings are continuously emerging .As a result ,underground space development project and deep excavation project become more and more. At the same time, the scale and complexity of deep excavation increasing bigger. they make the design and construction of deep excavation to face greater challenges. So structural design and deformation prediction of deep excavation has become an important research issue in the field of geotechnical engineering. In this paper, the deep excavation of Wanda Plaza, Wuhan is studied. And using survey data and structural design of deep excavation requirements to select reasonable foundation pit ,then to conduct the corresponding design. The meantime, as checking the deformation of deep excavation is a difficult problems ,it uses flac3D numerical simulation method to monitor the progress of deep pit’s excavation, construction .Then comparing the design results of the calculation and simulation results to obtained reasonable support structure design and control program of deformation.According to the actual situation and exploration data, the envelope of large diameter piles concrete piles, angle brace and top brace on the combination of a variety of internal support-based programs are selected, combined with slope Top large dump load shedding and the reinforcement measures of pit passive zone.1) The calculation part of the paper mainly introduce the design and calculation of large diameter concrete piles or bored pile, and the rest just briefly introduce the dumping load shedding, internal support structure, the pit design of passive zone strengthening and precipitation.2) With flac3D, successively study the model building, setting large diameter concrete piles, sloping excavation, soil nailing construction, pre-stressed cable (instead of internal support) construction and excavation for the foundation pit .Finally, conduct the deformation monitoring , output pile element, the internal force distribution analysis in anchorage unit .And then, provide the corresponding conclusions and recommendations.In this paper, conventional calculations and numerical simulation methods are used. And their results were very close. So it can give an effective and reasonable safety factor through the combination of these methods.Key words: deep excavation design flac3D numerical simulation目录第一章绪论 (1)第一节选题思路 (1)第二节设计流程 (1)第二章工程概况及场地工程地质条件 (3)第一节工程概况 (3)第二节场地工程地质条件 (4)第三章A-OPQRSA段基坑支护结构设计 (10)第一节设计依据 (10)第二节设计参数 (10)第三节A-OPQRSA段基坑支护方案选择 (11)第四节A-OPQRSA段基坑减载放坡设计 (13)第五节A-OPQRSA段基坑支护桩设计 (13)第六节A-OPQRSA段基坑地下水控制方案设计 (24)第四章基于flac3D基坑开挖模拟分析 (27)第一节关于flac3D的概述 (27)第二节基坑维护方案 (27)第三节计算模型及参数 (28)第四节初始应力计算 (29)第五节支护桩施工 (31)第六节模拟分层开挖和设定锚杆 (32)第七节设置采样记录变量 (34)第八节计算结果分析 (35)第五章结论与问题 (44)第一节结论 (44)第二节设计过程中存在问题 (45)致谢 (47)参考文献 (48)附录 (49)第一章绪论第一节选题思路深基坑工程设计是当今岩土工程界关注的热点话题，深基坑工程的难题在于对变形量的预测，基坑允许的变形、垂直位移的计算是比建筑物自身允许的沉降和沉降计算更为复杂的课题，但又是基坑工程尤其是在软土地区和工程地质、水文地质复杂地区无法回避的问题。

第43卷第6期　2009年6月上海交通大学学报J OU RNAL OF SHAN GHA I J IAO TON G UNIV ERSIT YVol.43No.6　J un.2009　收稿日期:2008207213基金项目:国家自然科学基金资助项目(50679041);上海市科学技术委员会资助项目(08201200903);江西省教育厅科学技术研究资助项目(G JJ 09367)作者简介:丁勇春(19792),男,江苏大丰市人,博士生,从事基坑与地下工程数值仿真方面的研究.王建华(联系人),男,教授,博士生导师,电话(Tel.):021*********,E 2mail :wjh417@. 文章编号:100622467(2009)0620976205基于FL AC3D 的基坑开挖与支护三维数值分析丁勇春1,　王建华1,　徐　斌1,2(1.上海交通大学土木工程系,上海200030;2.南昌工程学院土木工程系,南昌330029)摘　要:采用三维快速拉格朗日方法(FL AC3D )建立了考虑基坑分步开挖与支护全过程的三维动态计算模型,土体采用修正剑桥模型模拟,考虑了支护结构与土体的接触滑移作用,分析了基坑施工中围护墙变形、地表沉降、坑底隆起、坑外深层土体变形的基本特性.计算得到的地表沉降曲线与已有文献的经验沉降曲线基本一致,验证了计算结果的适用性.分析结果可为类似基坑工程的设计和施工提供有益参考.关键词:基坑开挖;数值分析;三维快速拉格朗日方法;修正剑桥模型;土与结构相互作用中图分类号:TU 473.2 文献标识码:AThree 2Dimensional Numerical Analysis of Braced ExcavationBased on F LAC 3DD I N G Yong 2chun 1,　W A N G J i an 2hua 1,　X U B i n1,2(1.Depart ment of Civil Engineering ,Shanghai Jiaotong University ,Shanghai 200030,China ;2.Depart ment of Civil Engineering ,Nanchang Instit ute of Technology ,Nanchang 330029,China )Abstract :Fast Lagrangian Analysis of Continua in 3Dimensions (FL AC3D )was employed to investigate t he deformatio n characteristics of a staged excavated and supported foundation pit.Modified Cam 2clay model was adopted to model t he soil behavior and t he interaction between soils and retaining st ruct ures was also taken into account in t he numerical model.The deflection of retaining walls ,t he settlement of gro und surface ,t he heave of excavation bottom ,and t he movement of t he deep st rata out side t he excavation were analyzed.The ground surface settlement curves of t he numerical model are basically consistent wit h t he empirical ones f rom t he existing literat ures ,so it s feasibility is proved.The numerical result s p rovide a usef ul reference for t he design and const ruction of similar deep excavation project s.Key words :braced excavation ;numerical analysis ;fast Lagrangian analysis of continua in 3dimensio ns (FL AC3D );modified Cam 2clay model ;soil 2st ruct ure interaction 基坑工程除了要保证基坑的整体性和支护结构的稳定安全外,还必须确定支护结构变形和基坑内外地层的变形,这样才能合理评估基坑施工对周围环境可能造成的不利影响,采取相应的工程措施,确保基坑施工的顺利进行,从而实现基于变形控制的基坑设计[1].目前基坑工程中常用的分析方法可以分为两类:行业和地方基坑设计规程所广为采用的竖向弹性地基梁法[2,3];基于连续介质力学的数值分析方法,如有限元法、有限差分法等.竖向弹性地基梁法只能对围护结构的内力和变形进行分析,不能计算支撑体系的内力和变形,更无法分析坑外地表沉降和坑底隆起变形,因此无法对基坑开挖引的环境影响进行分析和评价.数值分析方法可以弥补常规弹性地基梁法的不足.本文采用三维快速拉格朗日方法(FL AC3D )建立考虑支护结构与土体相互作用的基坑三维计算模型,土体采用修正剑桥模型,并考虑支护结构(围护墙和立柱桩)与土体间的接触问题,实现基坑分步开挖及支护的三维全过程动态分析.1　数值计算原理1.1　修正剑桥模型目前在岩土工程数值计算中常用的土体本构模型有两大类:一类是以Duncan 2Chang 模型为代表的非线性弹性模型;另一类是弹塑性模型,如Mohr 2Coulomb 模型、Drucker 2Prager 模型及修正剑桥模型等.Mohr 2Coulomb 模型和Drucker 2Prager 模型由于对加载和卸载采用同一模量,在得到合理围护结构侧向变形的同时,往往导致不合理的坑底隆起变形.修正剑桥模型能够反映土体加载与卸载模量的差异,考虑土体材料静水压力屈服特性和压硬性,在软土地基开挖分析中应用非常广泛.修正剑桥模型的屈服函数可表示为q 2p ′2+M 21-p ′0p ′=0(1)式中:q 为偏应力;p ′为平均有效应力;p ′0为先期固结压力;M 为p ′-q 平面内临界状态线(CSL )的坡度.修正剑桥模型屈服轨迹如图1所示.(a )p ′-q 平面内(b )主应力空间内图1　修正剑桥模型屈服轨迹Fig.1　Y ield locus of modified Cam 2clay model 修正剑桥模型的初始状态参数确定方法如下:(1)确定土体单元的竖向有效应力σ′v ;(2)根据现场静止侧压力系数K 0,确定土体单元侧向有效应力σ′h ,σ′h =K 0σ′v (2) (3)对每个土体单元,计算平均有效应力p ′和偏应力q :p ′=13(σ′1+σ′2+σ′3)q =12(σ′1-σ′2)2+(σ′2-σ′3)2+(σ′3-σ′1)2(3) (4)根据式(1)和超固结比OCR 计算先期固结压力p ′0,p ′0=p ′1+qM p ′2OCR (4) (5)确定正常固结线上p ′=1kPa 作用下单元的初始比体积υe ,0=υCSL +(λ-κ)ln 2(5)式中:λ为υ-ln p 平面内压缩线坡度;κ为υ-ln p 平面内回弹线坡度;υCSL 为临界状态线在p ′=1kPa 下的比体积.如根据土工试验能确定土样的不排水抗剪强度c u 和比体积υ0,则可根据下式确定参数υCSL[4],c u =M 2expυCSL -υ0λ(6) 如果已知土体颗粒比重G s 、土体液限w L 和土体塑性指数I p ,也可按下列经验公式确定参数[4]:υCSL =1+G s (w L +0.3I p )(7)1.2　土体与支护结构的接触采用数值方法进行基坑开挖与支护分析时涉及到是否需要考虑支护结构与土体间的相互接触问题.文献[5]中的研究结果表明,不考虑接触作用将导致计算结果与实际情况不符.主要原因可归结为:①无法反映支护结构与土体间由于接触面或土体的破坏引起的相互滑移;②结构与土体间存在拉应力,不能模拟结构与土体的脱离效应.FLAC3D 中的接触面单元具有单面特性,不同于G oodman 接触面单元.在每一计算时间步Δt 内,接触面节点与目标面之间的绝对法向侵入位移u n及相对剪切位移增量Δu s 均被计算,将其代入接触面本构方程就可以确定下一时刻(t +Δt )接触面上的法向力F (t +Δt )n和切向力F (t +Δt )s,F (t+Δt )n =k n u n A +σn A F (t+Δt )s=F (t )s+k s Δu (t+Δt/2)sA +τs A(8)式中:k n 为法向刚度;k s 为切向刚度;τs 为附加剪应力;A 为接触面节点的代表面积.接触单元服从库仑779　第6期丁勇春,等:基于FL AC3D 的基坑开挖与支护三维数值分析 剪切破坏准则和拉伸破坏准则.1.3　支护结构模型FLAC3D 提供了6种支护结构单元:梁单元、索单元、桩单元、壳单元、土工格栅单元和衬砌单元.对基坑而言,板桩式围护结构可采用衬砌单元模拟,桩列式和重力式挡土结构可采用实体模拟,水平内支撑如结构层楼板和混凝土或钢支撑可分别采用壳单元和梁单元模拟,坑外土钉和锚杆可采用索单元模拟,坑内立柱与立柱桩可采用桩单元模拟.2　计算模型及参数确定2.1　计算模型及边界条件分析模型取一方形基坑为研究对象,基坑平面尺寸为56m ×56m ,考虑模型的对称性后取1/4模型进行计算.基坑最大开挖深度H max =20m ,分5步开挖,每步开挖4m.1/4计算模型的三维尺寸为128m ×128m ×100m ,如图2所示.土体采用8节点6面体模拟,土体区段总数为13520,网格点总数为15309.模型外边界采用侧向约束,中心对称面采用对称边界,模型底部全约束.(a )整体模型(b )支护结构图2　计算模型Fig.2　Numerical model 计算中不考虑土体的分层以及基坑降水的影响,采用总应力法计算,相应的土体计算参数采用总应力指标.土体本构模型采用修正剑桥模型,不考虑地下水作用.土体重度γ=17.15kN/m 3,孔隙比e =1.2,侧压力系数K 0=0.5,κ=0.01,λ=0.14,M =1.2,泊松比μ=0.35.模型参数相当于上海软土地区第3层土体[6]的参数,按正常固结考虑(OCR =1).2.2　支护结构及接触面参数水平梁板支撑按刚度等效原则简化为壳单元,围护墙采用衬砌单元模拟,立柱和立柱桩采用桩单元模拟.结构单元总数为5175,结构单元节点数为2658.支护结构强度按C30混凝土考虑,考虑80%强度折减[7]后混凝土的弹性模量E c =24GPa ,μc =0.2,ρc =2500kg/m 3.坑边x =28m 和y =28m 处为两道连续墙,连续墙深度40m.坑内共设9根桩,桩直径0.8m ,间距8m ,桩长60m.共设5道水平支撑,支撑板厚0.12m.不同结构单元间采用共用节点法实现内力的传递.接触计算为物理非线性问题,有限元法在计算中往往存在收敛困难问题.FL AC3D 由于基于显式的全过程动力计算,因此,对不稳定问题不存在计算上的障碍.FLAC3D 中的衬砌单元与土体间的切向相互作用也具有单面特性,因而不能同时考虑围护墙与内外两侧土体的相互接触算法.在计算模型中采用以下近似处理办法:衬砌单元建立在墙外土体区域的外表面上以模拟围护墙与墙外土体的相互接触作用,围护墙与坑内土体的相互接触采用在坑内外土体间建立接触面单元,墙底处衬砌单元节点与坑内外土体网格点自由度耦合,假定墙底处结构单元节点与网格点变形协调.桩单元可直接实现桩与土界面的接触算法,通过加入桩端屈服弹簧可考虑桩的端承效应.参考文献[5],接触界面摩擦系数取为0.25,最大剪应力取为20k Pa.3　计算结果分析3.1　围护墙变形不同开挖深度(H )和不同平面位置(y )围护墙的侧向变形(δh )如图3所示.开挖至坑底后围护墙整体上抬约15.0mm ,墙体最大侧移-45.3mm,基坑开挖期间围护墙变形主要为侧向变形.(a )不同开挖深度 (b )不同平面位置图3　围护墙变形Fig.3　Deflection of retaining walls 由于第1道支撑的设置先于第1次4m 深土体的开挖,故第1次挖土时围护墙顶受到第1道水平支撑强大刚度的约束,墙顶位移几乎为零.后续阶段,墙体最大侧移随着开挖深度的增加而增大,同时最大侧移点位置亦随着开挖深度的增加而逐渐下879上　海　交　通　大　学　学　报第43卷　移,开挖至坑底后,墙身最大侧移点位于坑底开挖面附近.围护墙的最大侧移随着基坑开挖深度的增加呈非线性增长,但增长的幅度却有下降的趋势.围护墙的侧向变形具有明显的三维空间效应,基坑中心处(y =0)对称面墙体侧移最大,而基坑拐角位置(y =28m )墙体侧移最小,且为向坑内整体刚性平移.3.2　地表沉降不同开挖深度(H )和不同平面位置(y )坑外地表沉降(δv )如图4所示.图中,d 为计算点至坑边的距离.可见,地表未出现隆起,这是由于考虑了围护墙与土体间的接触特性,墙体与土体之间出现了滑移,墙体的上抬未引起地表的隆起.第1次挖土(H =4m )最大地表沉降点距坑边较近;随着挖土深度的不断增加,最大沉降点位置逐渐远离基坑,基坑开挖深度达12m 后,最大地表沉降点位置几乎不再变化.图4　地表沉降Fig.4　Ground surface settlement 地表沉降也具有空间效应,靠近基坑侧边中点(y =0)的地表沉降较大,最大地表沉降值为-36.3mm ;而靠近基坑拐角(y =28m )的地表沉降较小,最小地表沉降值仅为-16.4mm.地表沉降与围护墙的变形是相关联的,两者在空间上的变化规律保持一致. 图4(c )所示为中心对称面上不同开挖阶段相对地表沉降(δv /δv ,max )和相对距离(d/H )关系曲线与文献[8,9]中提出的经验沉降包络曲线的对比.可见,数值模拟结果与文献[9]中的更为吻合;d/H >2时本文与文献[8]中的相差较大.原因如下:文献[8]中假定地表沉降影响范围仅为2倍基坑开挖深度,而根据本文数值模拟结果和上海软土地区相关基坑的现场实测[10],上海软土地区地表沉降的影响范围一般为3～4倍基坑开挖深度.因此文献[8]中关于地表沉降影响范围的假定对上海软土地区的基坑而言显然偏小.由于计算模型的外截断边界采用侧向约束,边界面竖向可自由变形,引起了d/H >1.5范围计算所得相对地表沉降稍大于文献[9]中的包络曲线,这部分差异是由于计算模型的边界条件简化引起的.3.3　坑底隆起图5所示为不同开挖深度中心对称面(y =0)上坑底土体的隆起变形(δvb ).坑底土体的隆起主要由竖向开挖卸载效应引起,另外,墙体向基坑内的侧向变形也会进一步挤推坑内土体,造成坑底土体的回弹[11].坑底土体的隆起值随着开挖深度的增加呈非线性增长,但回弹增量有减小的趋势.图5　坑底隆起Fig.5　Heave of excavation bottom 开挖至坑底后坑底土体的竖向整体隆起量约160mm ,而围护墙的整体回弹仅15mm 左右,表明坑底土体与围护墙间发生了较大的相对滑移.如不考虑接触滑移作用,坑底回弹将引起过大的围护墙上抬,这显然不合理.另外,地表与围护墙相交处由于也考虑了接触作用,墙体的上抬才未引起相应位置地表的隆起,从而得到了与经验曲线基本吻合的地表沉降曲线.因此,在数值计算模型中有必要考虑支护结构与土体间的相互接触作用.3.4　坑外深层土体变形图6所示为开挖至坑底(H =20m )后中心对称面(y =0)上坑外深层土体竖向和侧向的位移场分布.可见,坑外浅层土体产生沉降,但最大沉降点位于距坑边一定距离处;坑外深层土体产生隆起,最大979　第6期丁勇春,等:基于FL AC3D 的基坑开挖与支护三维数值分析 隆起点位置位于土体与围护墙的接触处.邻近基坑围护墙处土体的侧向变形与墙体的变形相似,均为深层凸出型,最大侧向变形点深度位于基坑最终开挖面附近;随着至坑边距离的增加,土体最大侧向变形点位置逐步向地表面过渡,当距离超过一定范围,土体的最大侧移点位于地表面.(a )竖向变形(b )侧向变形图6　坑外深层土体变形Fig.6　Movement of strata outside excavation 对比图6中坑外土体的竖向变形与侧向变形可见,当基坑首层支撑刚度较大且先于浅层土体开挖而架设时,坑外浅层土体的变形以沉降为主,侧向变形较小.因此,对基坑周围建(构)筑物上部结构、浅埋地下市政管线及城市道路等而言,地表沉降和差异沉降是引起其损坏的主要原因;对基坑最终开挖面深度的坑外土体而言,侧向变形占主导,竖向变形较小,且存在隆起与下沉两种不同的形态.因此对深埋地下结构如高层建筑桩基础、城市高架桥梁桩基础、地铁区间隧道等而言,坑外深层土体的侧向变形是引起其损坏的主要原因.4　结　论(1)计算得到的坑外地表沉降曲线与已有文献的经验沉降包络曲线基本一致,验证了FL AC3D 在基坑开挖分析中的适用性和有效性.(2)基坑围护结构变形、地表沉降、坑底土体位移、坑外深层土体位移是相互关联的有机整体,基坑变形表现出明显的空间特性.(3)只要合理确定土体模型参数及土体与支护结构间的接触关系,三维数值分析在基坑工程设计施工方案优化和环境影响评估的应用中具有明显优越性.本文计算模型进行了一定简化,未考虑土体的分层,也未考虑降水引起的渗流和固结对基坑变形的影响,有待今后进一步研究.参考文献:[1]　刘建航,侯学渊.基坑工程手册[M ].北京:中国建筑工业出版社,1997.[2]　J G J 120299,建筑基坑支护技术规程[S].[3]　DBJ 08261297,基坑工程设计规程[S].[4]　Wood D M.Soil behaviour and critical state soil me 2chanics [M ].Cambridge :Cambridge UniversityPress ,1990.[5]　范　巍,王建华,陈锦剑.连续墙与土体接触特性对深基坑变形分析的影响[J ].上海交通大学学报,2006,40(12):211822121.FAN Wei ,WAN G Jian 2hua ,CH EN Jin 2jian.The evaluation of deformation induced by excavation con 2sidering the properties of diaphragm 2soil interface [J ].Journal of Shangh ai Jiaotong U niversity ,2006,40(12):211822121.[6]　D G J 0823722002,岩土工程勘察规范[S].[7]　谢百钩.粘土层开挖引致地盘移动之预测[D ].中国台湾:国立台湾科技大学,1999.[8]　Clough G W ,O ′Rourke T D.Construction inducedmovements of in situ walls [C]//Proceedings ,ASCE Conference on Design and Perform ance of E arth R etai 2ning Structures .New Y ork :ASCE ,1990:4392470.[9]　Hsieh P G ,Ou C Y.Shape of ground surface settle 2ment profiles caused by excavation [J ].C anadian G eotechnical Journal ,1998,35(6):100421017.[10]　丁勇春,王建华,徐中华,等.上海软土地区某深基坑施工监测分析[J ].西安建筑科技大学学报(自然科学版),2007,39(3):3332338.DIN G Y ong 2chun ,WAN G Jian 2hua ,XU Zhong 2hua ,et al .Monitoring analysis of a deep excavation inShanghai soft soil deposits [J ].Journal of Xi ’an U ni 2versity of Architecture &T echnology (N atural Science Edition),2007,39(3):3332338.[11]　刘国彬,黄院雄,侯学渊.基坑回弹的实用计算法[J ].土木工程学报,2000,33(4):61267.L IU Guo 2bin ,HUAN G Yuan 2xiong ,HOU Xue 2yuan.A practical method for calculating a heave of excavated foundation [J ].China Civil E ngineering Journal ,2000,33(4):61267.89上　海　交　通　大　学　学　报第43卷　。

基于FLAC-3D数值模拟分析逆作法的深基坑变形杨益飞;关群;孙若晗;高菊;张力【摘要】In this paper, the three-dimensional model is established by using ANSYS software, and the model is imported into FLAC-3D for numerical simulation and calculation of the deformation process of deep foundation pit excavation by top-down construction method.The results show that the deformation of diaphragm wall is in S shape and the maximum horizontal deformation is located above the 1/2-1/3 depth of the excavation surface.The adjacent soil settlement is in inverted bow shape and the maximum settlement is 6-8 m away from the pit.The bottom soil upheaval gradually increases with the soil excavation.The vertical displacement of pile top is roughly first slightly moved downwards and then uplifted with the excavation.The study results can provide reference for similar engineering.%文章利用ANSYS软件建立三维模型,导入FLAC-3D中对某逆作法的深基坑开挖变形过程进行数值模拟与计算.结果表明:深基坑的地下连续墙变形呈S形,最大水平变形位于开挖面以上1/2~1/3深度处;周边土体沉降呈倒弓形,距基坑6~8 m处沉降最大;坑底回弹随土体开挖逐渐增大;桩顶垂直位移随开挖大致呈先向下稍移再上抬的趋势.研究结果可为类似深基坑工程提供参考.【期刊名称】《合肥工业大学学报（自然科学版）》【年(卷),期】2017(040)004【总页数】6页(P522-527)【关键词】逆作法;深基坑;变形;数值模拟【作者】杨益飞;关群;孙若晗;高菊;张力【作者单位】合肥工业大学土木与水利工程学院,安徽合肥230009;合肥工业大学土木与水利工程学院,安徽合肥230009;合肥工业大学土木与水利工程学院,安徽合肥230009;合肥工业大学土木与水利工程学院,安徽合肥230009;合肥工业大学土木与水利工程学院,安徽合肥230009【正文语种】中文【中图分类】TU470.3近年来,逆作法在基坑、隧道、地铁等工程中得到了快速的发展。

基于FLAC3D的深基坑支护数值模拟应用基于FLAC3D的深基坑支护数值模拟应用摘要：本研究利用有限元数值模拟软件FLAC3D对深基坑支护进行了数值模拟。

通过对建筑施工现场实际参数的调查和相关文献的研究，我们建立了一个三维数值模型，并进行了不同支护方案的比对分析。

结果表明，采用加强钢支撑和地下连续墙的支护方案，能够有效地减少土壤变形、保证建筑物的安全性。

关键词：深基坑、支护、FLAC3D、数值模拟、加强钢支撑、地下连续墙1. 引言近年来，城市建设和基础设施建设快速发展，深基坑建设越来越普遍。

但是，深基坑工程施工过程中的安全问题一直备受关注。

其中，深基坑支护是工程中的一个重要环节。

为了确保施工期间的安全性，提高深基坑工程的质量和效率，数值模拟成为了深基坑支护研究的重要方法。

本文利用FLAC3D有限元软件，对深基坑支护进行了数值模拟，探究了不同支护方案对支护效果的影响。

2. 建立数值模型本文选取了某施工现场所需建设的深基坑作为研究对象，通过现场实际参数的调查和相关文献的研究，建立了一个三维数值模型，包括土体、加强钢支撑和地下连续墙等要素。

我们选用FLAC3D软件，采用三维非线性、非弹性有限元法建立了深基坑支护数值模型。

3. 数值模拟分析本文通过数值模拟的方法，分别对三种支护方案进行了分析和比较。

根据实际工程情境和可行性，将基坑侧壁加强钢支撑和地下连续墙结合起来，分别分析了它们分别对基坑支护的影响。

3.1 仅加强钢支撑采用钢支撑作为支护方案，计算结果表明，在基坑侧壁进行局部加强支撑的情况下，土体变形量和基坑下沉量都可控制在较小的范围内。

但是，当钢支撑的纵向间距较大时，局部土体变形较大。

3.2 仅连续墙支护采用地下连续墙作为支护方案，计算结果表明，连续墙的设置是很有必要的。

连续墙的加固作用可以有效地控制土体侧向位移和基坑下沉量。

但是，如果连续墙质量不好，可能会导致工程安全事故发生。

3.3 加强钢支撑和连续墙结合支护我们采用加强钢支撑和地下连续墙结合的支护方案，计算结果表明，在相同的施工条件下，结合支护方案的基坑下沉量更小，变形量也更小。

基于FLAC^(3D)的深基坑开挖变形数值模拟研究
杨博维
【期刊名称】《水资源与水工程学报》
【年(卷),期】2014(25)2
【摘要】结合沈阳北站综合交通枢纽改扩建工程的地下新建工程Ⅱ区基坑工程实例,运用有限差分软件FLAC3D模拟了基坑在4种不同支护方案下开挖过程,对4种开挖过程所引起的基坑坡顶和坡底的位移、坑壁水平方向的位移进行了分析。

数值模拟结果表明:坑底隆起是由于垂直方向卸载而引起的,采用混凝土围护桩和预应力锚索联合支护可以有效地控制基坑的变形。

研究结果和工程实践证明本工程的支护方案是有效可行的。

【总页数】5页(P24-28)
【关键词】基坑;FLAC3D;变形分析;支护结构
【作者】杨博维
【作者单位】辽宁工程技术大学土木与交通学院
【正文语种】中文
【中图分类】TV314;TV551.4
【相关文献】
1.基于FLAC-3D数值模拟分析逆作法的深基坑变形 [J], 杨益飞;关群;孙若晗;高菊;张力
2.基于FLAC3D的深基坑开挖与支护数值模拟应用 [J], 徐凌;陈格际;刘帅
3.基于FLAC3D数值模拟的地铁异型深基坑开挖阳角效应研究 [J], 沈启炜
4.基于FLAC^3D的地铁车站深基坑开挖与支护数值模拟 [J], 孔汇川;田淑朋;白晓明;
5.基于FLAC 3D的大型地下洞室群分步开挖稳定性数值模拟研究 [J], 张亚勤;孔维秋;栗剑;岳好真;李勇
因版权原因，仅展示原文概要，查看原文内容请购买。

基于FLAC 3D数值模拟的煤矿掘进巷道优化支护设计摘要：煤矿掘进巷道的顶板支护是煤巷安全管理的重点，如何合理优化支护设计，做好顶板安全管理是煤巷管理人员的工作重心。

本文从笔者的工作实际出发，通过观察现场地质条件，制定优化支护设计技术措施，并优化支护设计方案，结合FLAC 3D数值模拟技术，设计出更为合理的支护方式。

关键词：FLAC 3D；数值模拟；支护设计一、优化支护设计技术措施依据某工作面地质力学评估与巷道围岩变形分析结果提出以下支护优化方向。

1、顶板采用预应力锚索全长锚固技术由于顶板发生离层的位置主要在顶板软弱夹层，加上锚索预紧力大且延伸率低，使锚索受力高，易发生破断现象。

预应力锚索全长锚固技术，是在原中空注浆锚索的基础上进行了改进升级，它采用专用的无机复合锚固材料和高压注浆设备进行施工。

锚索先在迎头进行树脂端锚[1]，不用封孔，张拉预紧施加预应力，滞后通过锚索的中空结构注入无机复合锚固材料，当孔口流出无机锚固剂时停止注浆实现锚索的全长锚固。

采用该项支护技术的优点是顶板锚杆和锚索全部实现了全长锚固，锚杆、锚索同步承载，协调一致。

全长锚固锚索与端锚锚固锚索相比，提高了锚索的抗剪切能力和系统的刚性，消除了端锚锚索在非锚固段应力集中，将载荷进行了分散，也消除了锚索断裂弹出的现象。

围岩变形控制效果好，可降低锚索的支护密度，提高施工速度。

2、巷帮底角下扎45°锚杆控制底鼓底角锚杆的可起到三方面的作用：（1）底角锚杆将限制底角处岩层向巷道内发生水平位移，随着底板的松动及变形，底板上的水平应力也将逐渐传递到底角锚杆上。

这种情况下，底板岩石受的应力较无底角锚杆时会小很多，从而底板岩石不易发生破坏。

（2）施工底角锚杆有利于底板水平应力向底板深部转移，起到一定的卸压作用。

施加45°的底角锚杆，底板岩层承担水平应力将降低，水平应力将向深部转移，也有利于底板的稳定。

（3）施加底角锚杆有利于巷帮的稳定。

FLAC3D数值模拟在高位钻孔参数优化中的应用
陈康;杨志超;王福海;王威;王庆牛
【期刊名称】《能源技术与管理》
【年(卷),期】2015(040)005
【摘要】利用FLAC3D软件对芦岭煤矿Ⅱ817工作面顶板跨落情况进行模拟分析,初步判定了采空区上覆岩层冒落带、裂隙带和弯曲下沉带的高度,综合判定得出裂隙带岩层范围为18～37 m,合理确立了高位钻孔终孔位置为距8煤顶板20～32 m 岩层范围内.在抽采效果考察中,平均每个钻孔的抽采参数瓦斯体积分数均在35％以上,钻孔抽采瓦斯体积分数最大为55.3％,采空区瓦斯治理效果显著.
【总页数】3页(P195-197)
【作者】陈康;杨志超;王福海;王威;王庆牛
【作者单位】安徽省煤炭科学研究院,安徽合肥232000;中北煤化工有限公司色连二矿,内蒙古鄂尔多斯市017000;安徽省煤炭科学研究院,安徽合肥232000;安徽省煤炭科学研究院,安徽合肥232000;安徽省煤炭科学研究院,安徽合肥232000【正文语种】中文
【中图分类】TD712
【相关文献】
1.瓦斯抽采高位钻孔参数优化技术研究与应用 [J], 刘伟
2.高位钻孔布置参数优化及其应用 [J], 成建林;张必同;胡新社;成宽亮;原双丰;张海军
3.综放工作面高位定向钻孔层位参数优化与应用实践 [J], 彭冬; 段会军; 赵永哲
4.顶板高位钻孔抽放瓦斯参数优化及应用 [J], 王雪瑞
5.水压预裂工作面瓦斯抽采高位钻孔参数优化及应用 [J], 祝琨;刘文岗;刘庆林;党业龙;余杰;刘辉;董浩;吴越;王宁
因版权原因，仅展示原文概要，查看原文内容请购买。

基于FLAC3D的深基坑支护数值模拟应用随着城市现代化建设脚步的加快,城市人口数量的不断增多,城市建筑开始向更高和更深层次发展。

越来越多的高层建筑,超高层建筑,地铁,地下停车场,地下超市等在城市中建造,伴随而来的是深基坑工程的出现。

基坑研究由来已久,能够因地制宜选择合适合理的基坑支护方案以及施工开挖方式方法等是检验每个工程技术人员技术能力的最高标准。

本文依托呼和浩特市轨道交通1号线地铁博物馆站深基坑工程。

运用有限差分软件FLAC3D对深基坑开挖各阶段进行数值模拟计算,结合现场实际监测数据。

研究分析了地铁深基坑支护结构、支撑轴力、周边地表沉降变形变化特征。

研究发现:该深基坑开挖引起的支护结构地下连续墙变形曲线为“弓形”变形形态,引起的周边地表沉降变形曲线为“漏斗形”变形形态。

数值模拟结果与实际监测数据不同,但总体趋势基本一致,表明深基坑开挖的数值模拟具有一定的预测效果。

结合现场实际监测工作,能够更好的指导施工。

该深基坑支护方案及深基坑开挖方式方法合理可行,对今后该地区类似工程有一定的参考价值。

基于Ｆｌａｃ３Ｄ数值模拟的深部采场结构参数优化及应用收稿日期：２０２１－１１－３０；修回日期：２０２２－０３－１０基金项目：中国博士后科学基金面上项目（２０２１Ｍ７０２０１５，２０２１ＭＤ７０３８７４）作者简介：周亚博（１９８６—），男，河南郑州人，工程师，硕士，从事金属矿山开采理论研究等工作；内蒙古自治区锡林郭勒盟东乌珠穆沁旗满都镇，锡林郭勒盟山金阿尔哈达矿业有限公司，０２６３９９；Ｅ ｍａｉｌ：ｚｈｏｕｙａｂｏ＠ｓｄ ｇｏｌｄ．ｃｏｍ周亚博（锡林郭勒盟山金阿尔哈达矿业有限公司）摘要：针对阿尔哈达铅锌矿上向水平分层充填采矿法采场结构参数的确定开展相关研究，采用Ｆｌａｃ３Ｄ数值模拟软件建立了留顶柱不留底柱、留底柱不留顶柱共２种数值模型，分步模拟计算各方案在回采过程中的应力和位移分布规律，分别从拉应力和底板上鼓量、顶板下沉量分析比较方案的优缺点。

研究结果表明：首采分层开挖后，２种方案采空区的顶底板均出现了拉应力集中，采空区两侧围岩出现压应力集中区域；留底柱不留顶柱的方案拉应力出现明显下降，仅为２．３０ＭＰａ左右，因此选择留底柱不留顶柱的方案作为阿尔哈达铅锌矿上向水平分层充填采矿法的采场结构布置形式。

关键词：数值模拟；深部开采；采场结构参数；顶柱；底柱；充填采矿法 中图分类号：ＴＤ８５３．３４文献标志码：Ａ开放科学（资源服务）标识码（ＯＳＩＤ）：文章编号：１００１－１２７７（２０２２）０５－００５１－０５ｄｏｉ：１０．１１７９２／ｈｊ２０２２０５１０引　言采场结构参数不仅影响采矿生产安全，更与采矿作业效率密切相关［１］。

进行采场稳定性分析及结构参数优化研究，推荐合理的采场结构参数，对矿山安全生产具有重要的指导意义［２－３］。

学者们对采场结构参数的优化进行了大量的研究，刘文胜等［４］基于正交试验法，采用相似模拟试验的手段对眼前山铁矿的采场结构参数进行了优化；李守强等［５］采用数值模拟的方法对喀拉通克铜镍矿Ⅰ号矿体采场结构进行了设计优化；李江等［６］采用数值模拟的方法对大尹格庄金矿９种采场结构参数进行了分析，选择出最优方案；马春德等［７］采用ＡＮＳＹＳ与Ｆｌａｃ３Ｄ耦合模拟的方法分析了湘西金矿沃溪矿区的采场结构参数，优选出深部采场最优跨度。

基于FLAC3D的千米埋深的煤层底板影响深度模拟
董东林;洪益清;钱增江
【期刊名称】《煤炭工程》
【年(卷),期】2010(000)002
【摘要】在获取力学参数的基础上,基于FlAC3D建立了底板突水渗流-损伤耦合数值模型.确定了12煤层底板破坏深度范围,分析了水压、开拓进度(100m步距)对破坏深度的影响.研究成果认为,范各庄12煤的底板影响深度随着开采的步距的加大发生有规律的变化,在开采20a时,底板影响深度达到23m左右,这与范各庄矿的12煤底板的影响深度为18m的实际非常接近,为该矿的后续开采提供科学依据.
【总页数】4页(P71-74)
【作者】董东林;洪益清;钱增江
【作者单位】中国矿业大学煤炭资源与安全开采国家重点实验室,北京,100083;开滦集团有限责任公司,河北,唐山,067018;开滦集团有限责任公司,河北,唐山,067018【正文语种】中文
【中图分类】TD325
【相关文献】
1.煤层底板破坏规律的数值模拟研究——基于FLAC3D技术 [J], 王志光;郝青山
2.煤层底板破坏规律的数值模拟研究——基于FLAC3D技术 [J], 王志光;郝青山;
3.FLAC3D数值模拟大采深煤层底板采动破坏深度 [J], 张勤
4.基于FLAC3D煤层底板采动变形规律的数值模拟研究 [J], 刘磊;王娟;柳汉丰;刘
远滨
5.基于FLAC3D煤层底板扰动规律的数值模拟研究 [J], 田磊;田豪杰;殷德耀因版权原因，仅展示原文概要，查看原文内容请购买。