FLAC3D强度折减程序
基于FLAC_3D_的强度折减法和点安全系数法对比
图 1 空间应力状态及摩尔库伦强度包络线 强度储备法的点安全系数公式推导如下 :
Kp ( Fs1 ) = Kp ( Fs1 ) = | FS | | FS | = | BF| | AB | cos < ( 8) ( 9)
φ]tanφ + 2 c [ (σ 1 +σ 3 ) - (σ 1 - σ 3 ) sin (σ φ 1 - σ 3 ) cos
Flac3D软件的快速拉格朗日差分法及编制 Fish 强度折减和点
55
将式 ( 4) 和式 ( 5) 带入式 ( 3) , 可得 :
Kp =
2 2 ( l2σ 1 + mσ 2 + nσ 3) f + c 2 2 lσ 1 2 2 2 2 2 2 2 2 + mσ 2 + nσ 3 - (lσ 1 + mσ 2 + nσ 3)
c′= c/ Fs ( 1) ( 2)
式 ( 7) 中 , m = 0 , n = ±
σ f 1 + c 。 σ f (σ + 1 3) + 2 c
<′= arctan (tan </ Fs )
目前强度折减法中判断岩土体失稳破坏的标准通常有 : 迭代求解的不收敛性 、 剪应变增量和塑性区的连通状态 、 岩 土体控制点位移与折减系数的关系曲线等 。这些判别标准 都是根据边坡局部或整体状态的突变现象来认定边坡的稳 定性 。如位移突变可以表征整个边坡体的一种状态变化 ,由 位移缓慢增长急剧增大的临界点 ,即边坡最大位移与折减系 数关系曲线上的转折点 ( 由位移缓慢增长急剧增大的临界 点) 可较好地判定边坡达到极限平衡状态 。
2 . College of Civil and Hydraulic Engineering , Tsinghua University , Beijing 100084 , China)
岩坡稳定的三维强度折减法分析
第25卷 第12期岩石力学与工程学报 V ol.25 No.122006年12月 Chinese Journal of Rock Mechanics and Engineering Dec .，2006收稿日期：2005–11–28；修回日期：2006–01–13基金项目：国家重点基础研究发展规划(973)项目(2002CB412702)；国家自然科学基金资助项目(50279051)作者简介：陈 菲(1980–)，女，2002年毕业于武汉大学水利水电建筑工程专业，现为硕士研究生，主要从事边坡稳定方面的研究工作。
E-mail ：feierchen_sunny@岩坡稳定的三维强度折减法分析陈 菲1，2，邓建辉3(1. 中国科学院 武汉岩土力学研究所，湖北 武汉 430071；2. 国电大渡河流域水电开发有限公司，四川 成都 610016；3. 四川大学 水力学与山区河流开发保护国家重点实验室，四川 成都 610065)摘要：目前，强度折减法基本上用于边坡的二维稳定分析。
为研究其在岩石边坡三维稳定分析中的适用性，选择平面滑动和楔形体滑动2个经典算例，运用强度折减法求解其安全系数，并与E. Hoek 和J. W. Bray 给出的解析解进行对比。
计算结果表明，数值解与解析解非常接近，模拟的滑动方向与理论假定也基本一致。
算例中，结构面使用实体单元模拟，为此进一步探讨结构面单元厚度和网格密度对计算精度的影响；结构面均为平面，其厚度为块体高度的1/10～/200，厚度与单元数对安全系数的计算结果影响不大。
对特征点位移准则和塑性区贯通准则进行比较分析，认为结构面塑性区贯通收敛准则在理论上与极限平衡理论较为吻合，实际应用上也易于操作，作为收敛准则较为合理。
三滑面五面体滑动算例进一步验证该方法的适用性。
关键词：边坡工程；岩坡；三维稳定性分析；强度折减法中图分类号：P642.22；O 242 文献标识码：A 文章编号：1000–6915(2006)12–2546–06THREE-DIMENSIONAL STABILITY ANALYSIS OF ROCK SLOPE WITHSTRENGTH REDUCTION METHODCHEN Fei 1，2，DENG Jianhui 3(1. Institute of Rock and Soil Mechanics ，Chinese Academy of Sciences ，Wuhan ，Hubei 430071，China ；2. Dadu River Hydropower Development Co .，Ltd .，Chengdu ，Sichuan 610016，China ；3. State Key Laboratory of Hydraulics and Mountain River Engineering ，Sichuan University ，Chengdu ，Sichuan 610065，China )Abstract ：The strength reduction method(SRM) is mainly used for two-dimensional slope stability analysis. In order to verify its feasibility for the stability analysis of three-dimensional issues ，two benchmark problems ，i.e. planar slide and wedge slide ，are selected ；and the safety factors are calculated by SRM. The obtained safety factors agree well to the analytical solutions by E. Hoek and J. W. Brown ；and the modeled sliding direction is also almost the same as the hypothetical one. In the analysis ，solid elements are used to model the behavior of structural planes ，so the influences of their thickness and element intensity on modeling precision are further investigated. For blocks cut by planar structural planes ，the influence is limited when the thickness is 1/10–1/200 of the block height. Compared with the displacement criterion of characteristic points ，the entirely run-through criterion of plastic zone is more rational for its theoretical coherence to the theory of limit equilibrium and its easy use in application. The stability of a pentahedron block is further analyzed to verify the applicability of the method. Key words ：slope engineering ；rock slope ；three-dimensional stability analysis ；strength reduction method 1 引 言结构面是控制岩石边坡稳定的重要因素之一。
基于FLAC3D的三维边坡稳定性强度折减法计算效率改进算法及其应用
2、边坡的稳定性系数等指标，可以对边坡的稳定性进行评价； 3、破裂面的位置和形状，有助于确定治理措施和优化设计方案。
在结果分析过程中，需要以下几个方面：
1、破裂面的预测：FLAC3D可以准确预测破裂面的位置和形状，但需要考虑 地质条件的复杂性和计算精度的限制；
2、稳定性系数的敏感性：稳定性系数受到多种因素的影响，如土体性质、 荷载等，需要对这些因素进行敏感性分析；
根据现场调查和工程地质勘探，收集到了边坡的岩土力学参数、水文地质条 件、可能的荷载和边界约束条件等数据。在此基础上，采用有限元法建立了边坡 的计算模型，详细模拟了边坡内部的物理和力学过程。通过不断调整模型的材料 参数和边界条件，使得模型的计算结果与实际观测数据尽可能接近。
在模型验证可靠后，采用整体强度折减法对边坡的稳定性进行分析。根据计 算结果，该边坡的稳定系数为1.25，属于基本稳定状态。失稳破坏模式主要为边 坡顶部出现拉裂缝，同时局部岩体出现滑动。为了提高该边坡的稳定性，建议采 取适当的加固措施，如增加锚杆、喷射混凝土等。
引言
边坡稳定性分析是岩土工程领域中非常重要的一个方面，对于保障工程项目 安全具有至关重要的作用。在工程建设过程中，边坡失稳可能会导致严重的工程 事故，因此对边坡稳定性进行准确的分析和评估具有重要意义。随着计算机技术 的不断发展，数值模拟方法在边坡稳定性分析中得到了广泛应用。本次演示将探 讨基于FLAC3D（Fast Lagrangian Analysis of Continua in 3 Dimensions） 软件的边坡稳定性分析方法。
4、荷载施加：根据实际情况，施加一定的荷载，如自重、外力等。 5、数值模拟：利用FLAC3D进行数值模拟，得到边坡的位移、应力等响应。
6、稳定性评估：根据模拟结果，对边坡的稳定性进行评估，计算稳定系数 等指标。
边坡强度折减法
基本原理：强度折减法中边坡稳定的安全系数定义为：使边坡刚好达到临界破坏状态时，对岩、土体的抗剪强度进行折减的程度，即定义安全系数为岩土体的实际抗剪强度与临界破坏时的折减后剪切强度的比值。
强度折减法的要点是公式1、2来调整岩土体的强度指标C 和φ(式中，F C 为折减后的粘结力，F φ为折减后的摩擦角，trial F 为折减系数），然后对边坡稳定性进行数值分析，不断地增加折减系数。
反复计算，直至其达到临界破坏，此时得到的折减系数即为安全系数S F 。
公式如下： trial F F C C /= （1） t a n =F φ-1)/)((tan trial F φ（2)实现过程：目前尚无统一的边坡失稳判据，现行的边坡失稳判据主要有以下几种：1 以数值计算的收敛性作为失稳判据2 以特征部位位移的突变性作为失稳判据3 以塑性区的贯通性作为失稳判据在FLAC3D 中求解安全系数时，单次安全系数的计算过程主要采用的是第一种失稳判据。
假设数值计算模型所有非空区域都采用摩尔-库伦本构模型，便可使用命令Solve fos 来求解安全系数：首先，通过给粘结力设定一个大值来改变内部应力，以找到体系达到力平衡的典型时步r N ;接着，对于给定的安全系数s F ，执行r N 时步，如果体系不平衡力与典型内力比率R 小于10-3，则认为体系达到力平衡。
如果不平衡力比率R大于10-3，再执行r N时步，直至R小于10-3退出当前计算，开始新一轮折减计算过程。
除上述以力不平衡比率小于10-3作为终止条件外，FLAC3D还采用：1 前后典型时步运算结束时的不平衡力比率R差值小于10%2 强度折减后的计算过程已运行了6个典型时步r N作为计算终止条件计算过程中，只要满足上述三个标准中的任何一个，便退出当前计算。
这样做的目的只要是为了控制整个强度折减法循环计算过程中的求解时间。
可以从这几个方面判断：边坡沿滑动面产生滑动、软弱面处产生的沿X方向的位移是否最大、剪切应变增量云图、安全系数、剪切应变增量云图、变形矢量图及速度矢量图、水平位移、竖直位移、垂直应力、最大不平衡力、在坡顶边缘和坡脚处设置监测点（水平应力竖直应力位移）。
FLAC 3D强度折减理论在边坡稳定分析中的应用
FLAC 3D强度折减理论在边坡稳定分析中的应用季聪;佴磊;马宏;佘小光;刘录君【摘要】FLAC3D是岩土工程中广泛应用的软件.本文主要介绍强度折减理论的基本原理,利用基于强度折减法的FLAC3D软件,对某高速公路K377+ 230 ～ K377+ 520段滑坡体进行数值模拟.计算边坡在自重作用达到初始化平衡状态时,边坡的横向与竖直方向的应力与位移,分析出边坡的破坏机制为牵引式滑坡,通过强度折减理论计算出边坡的稳定系数为0.97,由剪切应变增量云图确定边坡的滑动面.%FLAC 3D is a widely applied software in geotechnical engineering. We mainly introduce the principle of strength reduction theory. Numerical simulation on the K377 + 230 ~ k377 + 520 section of a highway has been made by FLAC 3D sofeware based on strength reduction theory. To calculate the slope lateral and vertical stress and displacement when the slope gets to the initialization equilibrium condition in dead weight function, then to analyze the result that the slope failure mechanism was retrogressive landslide; after that to calculate the slope stability coefficient as 0. 97 through the strength reduction theory, finally to determine the sliding surface by the shear strain increment chart.【期刊名称】《世界地质》【年(卷),期】2013(032)001【总页数】7页(P158-164)【关键词】FLAC3D;强度折减理论;滑动面;稳定系数【作者】季聪;佴磊;马宏;佘小光;刘录君【作者单位】中水东北勘测设计研究有限责任公司,长春130021【正文语种】中文【中图分类】P642.220 引言边坡工程是岩土工程领域的一个分支，其稳定性研究一直是个热点课题。
基于FLAC3D数值模拟求解边坡安全系数
基于FLAC3D数值模拟求解边坡安全系数作者：刘丰来源：《山东工业技术》2014年第23期摘要：本文结合鹤大高速公路某段路堑边坡，利用FLAC3D软件模拟模拟岩土边坡，求解该边坡的安全系数，并将分析结果与简化毕肖普条分法分析结果进行比较验证，结果发现利用数值模拟分析可以有效的解决工程复杂等因素的影响，具有很好的适用性。
关键词：强度折减法；边坡安全系数；FLAC3D0 引言在公路的设计与建设工程中，岩土边坡的稳定性历来是岩土工程领域的一个热点研究课题，也是公路施工以及维护考虑的重点，公路边坡主要的荷载来源是其自身的重力以及公路的上部的各种荷载，人们通常采用安全系数来评价其稳定性状态。
安全系数是以极限平衡法为基础的一种评价指标，由于其原理简单、物理意义明确，而成为边坡稳定性分析中的重要、关键指标。
1 强度折减法强度折减法是在边坡刚好达到临界破坏状态的时候，对岩土体的抗剪强度进行折减并达到一定的程度[1]。
边坡的安全系数通常的定义是研究的岩土体实际抗剪强度与临界破坏时的折减后剪切强度的比值[3]。
强度折减法主要是通过公式（1）和（2）来调整岩土体的c和。
在数值模拟分析过程中，通过不断的增加折减系数，并进行反复的计算，知道其达到了临界的破坏，这时的折减系数就是安全系数Fs。
其中：为折减后的内聚力、为折减后的内摩擦角、为折减的系数。
2 数值模拟分析2.1 工程简介本文实例分析采用的是鹤大高速公路某段，线位带自东北向西南展布于长白山山脉南部山区腹地，地势总体东北高西南低。
研究区设计带属于湿带大陆性气候，四季分明，气温变幅较大，降水量随季节变化明显。
沿线地层构造属华北地层分区之辽东分区，地层以前震旦纪的混合岩、侵入岩和变质岩为主，山间河谷及其两侧上覆第四纪松散土体。
主要为上更新冲积层，上部岩性具有二元结构，上部为棕黄色、褐色亚砂土、亚粘土，下部为黄褐色、灰白色砂砾石、砾卵石和漂石层。
2.2 模拟网格模型模拟的边坡坡高为8m，基底土体厚度为3m，边坡 H/B = 1/1，坡角距离左边界2m，坡顶距离至右边界取10m。
基于极限平衡法和强度折减法的边坡稳定性对比分析
基于极限平衡法和强度折减法的边坡稳定性对比分析摘要：为分析某矿露天开采时最终边坡的稳定性，分别采用极限平衡法和强度折减法计算边坡的安全系数，采用理正软件和FLAC3D软件作为计算工具，建立边坡模型，分别运用Morgenstern-Price法和强度折减法对最终边坡的稳定性进行计算，依据《非煤露天矿边坡工程技术规范》(GB51016-2014)对最终边坡的稳定性进行评价。
分析结果表明：采用Morgenstern-Price法和强度折减法对边坡的稳定性分析结果基本一致，该矿山最终边坡稳定性较好。
关键词：边坡稳定性；Morgenstern-Price法；强度折减法0引言边坡稳定性一直是露天矿山面临的重大问题，时刻影响着矿山的安全生产，边坡稳定性分析中，先后发展了工程地质分析法、类比法、极限平衡法、数值分析方法和不确定性分析方法(可靠性方法、模糊数学方法、灰色理论方法、神经网络方法等)，随着计算机技术的发展，数值分析方法运用越来越广，目前国内外边坡稳定性分析法主要以极限平衡法和数值分析法为主。
极限平衡法主要有瑞典条分法、Bishop法、Janbu法、Spencer法、Morgenstern-Price法、Sarma法、平面直线法和不平衡推力传递法。
极限平衡法把边坡上的滑体视为刚体，利用滑体的静力平衡原理分析边坡在各种极限破坏模式下的受力状态，并以边坡滑体上的抗滑力和下滑力之间的比值定义为安全系数。
强度折减方法的基本原理是将岩土材料的黏聚力和内摩擦角等抗剪强度参数进行折减，用折减后的参数进行边坡的稳定性分析计算，不断降低强度参数直至边坡失稳破坏为止，破坏时的折减数值即为边坡的安全系数。
极限平衡法中，Morgenstern-Price法既能满足力平衡又满足力矩平衡条件，是国际公认的最严密的边坡稳定性分析方法[1]。
数值分析方法有如有限元法（ANSYS、Plaxis、ABAQUS）、离散元法（PFC、3DEC）、边界元法（BEM）和拉格朗日元法（FLAC），FLAC3D是基于连续介质快速拉格朗日差分法编制而成的数值模拟计算软件, 是目前岩土工程界应用最为广泛的数值模拟软件之一，该程序采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
基于FLAC3D的强度折减法程序改进研究
ｏ ｐ ｅ ｒ ａ ｔ ｉ ｏ ｎ（ ｓ ｔ ｒ ｅ ｎ ｇ ｔ ｈ ｒ ｅ ｄ ｕ ｃ ｔ ｉ ｏ ｎ ） ， ａ ｎ ｄ ｉ ｎ ｅ ａ ｃ ｈ ｉ ｔ ｅ ｒ ａ ｔ ｉ ｏ ｎ， ｉ ｎ ｉ ｔ ｉ ａ ｌ ｓ ｔ ｒ ｅ ｓ ｓ ｍｕ ｓ ｔ ｂ ｅ ｒ ｅ ｃ ａ ｌ ｃ ｕ ｌ ａ ｔ ｅ ｄ ， ｗ ｈ ｉ ｃ ｈ ｉ ｓ ｕ ｎ ｎ ｅ ｃ ｅ ｓ ｓ ａ ｒ ｙ ａ ｎ ｄ
Ａｂｓ ｔ ｒ ａ ｃ ｔ ： Ｔｒ ａ ｄ ｉ ｔ ｉ ｏ ｎ ａ ｌ ｓ ｅ ｌ ｆ － ｍａ ｄ ｅ ｓ ｔ ｒ ｅ ｎ ｇ ｔ ｈ ｒ ｅ ｄ ｕ ｃ ｔ ｉ ｏ ｎ ｐ ｒ ｏ ｇ ｒ ａ ｍ ｉ Ｓ ｂ ａ ｓ ｅ ｄ ｏ ｎ ｂｉ ｓ ｅ ｃ ｔ ｉ ｏ ｎ ｍｅ ｔ ｈ ｏ ｄ ａ ｎ ｄ ｓ ｔ ｒ ｅ ｎ ｇ ｔ ｈ ｒ ｅ ｄ ｕ ｃ ｔ ｉ ｏ ｎ ｔ ｅ ｃ ｈ ｎｉ ｑ ｕｅ ｔ ｏ ａ ｎ ａ ｌ ｙ ｚ ｅ ｓ ｌ ｏ ｐｅ ｓ ｔ ａ ｂ ｉ ｌ ｉ ｔ ｙ ．Ｈｏ ｗｅ ｖ ｅ ｒ ，ｔ ｈ ｉ ｓ ｍｅ ｔ ｈ ｏ ｄ ｈａ ｓ ｄ ｅ ｆ ｅ ｃ ｔ ｓ ： ｉ ｎ ｉ ｔ ｉ ａ ｌ ｓ ｔ ｒ ｅ ｓ ｓ ｃ ｈａ ｎ ｇ ｅ ｓ ａ ｆ ｔ ｅ ｒ ｔ ｈｅ ｆ ｉ ｒ ｓ ｔ ｉ ｔ ｅ ｒ ａ ｔ ｉ ｖ ｅ
强度折减法在FLAC 3D数值模拟中的应用
内，从而有效防止石灰凝结问题。另外，在开展搅拌工作时，所选 参考文献：
择的石灰内应确保 没 有 杂 质、异 物，其 中 所 含 有 的 ＣａＯ含 量 在 ［１］ 管义能．公路桥梁施工中软土地基施工技术及应用实践探
８０％以上。与此同时，施 工 人 员 要 妥 善 储 存 石 灰，储 存 期 限 应 在
究［Ｊ］．建材与装饰，２０１６（２５）：２５９２６０．
三个月之内，石灰所具备的液性指数应保持在 ７０％左右。在运用 ［２］ 姚 彦．道路桥梁工程中软土地基的施工处理方法研究
深层石灰搅拌技术时，施工人员同样要将前期施工准备工作落实
［Ｊ］．住宅与房地产，２０１６（３）：２１７．
Ｃｏｎｓｔｒｕｃｔｉｏｎｔｅｃｈｎｏｌｏｇｙａｎａｌｙｓｉｓｏｆｓｏｆｔｓｏｉｌｆｏｕｎｄａｔｉｏｎｉｎｂｕｉｌｄｉｎｇｐｒｏｊｅｃｔ
·８２·
第２０４４ １卷８年第１３００期 月
山 西 建 筑
ＳＨＡＮＸＩ ＡＲＣＨＩＴＥＣＴＵＲＥ
ＯＶｃｏｔｌ．．４ ４２Ｎ０１ｏ８．３０
文章编号：１００９６８２５（２０１８）３０００８２０３
强 度 折 减 法 在 ＦＬＡＣ３Ｄ数 值 模 拟 中 的 应 用
术，严格依据施工流程开展施工工作，对管道内部进行检查，并对 ４ 结语
管道内存在的 水 垢 进 行 清 理，在 确 保 管 道 顺 畅 的 基 础 上 开 展 施
在建筑行业迅速发展的当下，人们对建筑工程质量及性能予
工；另外，施工方要加大检查力度，对施工人员的操作情况进行检 以高度关注。影响房建工程质量的因素有很多，就软土地基而
章显斌１ 余以强２
（１．福州职业技术学院，福建 福州 ３５０１０８； ２．浙江省交通运输科学研究院，浙江 杭州 ３１１３０５）
FLAC3DV3.0边坡稳定剪切强度折减法的多线程并行算法
FLAC3DV3.0边坡稳定剪切强度折减法的多线程并行算法瞿秀祥;林杭;王敏;张道勇【摘要】The safety factor of slopes was determined by the critical reduction parameter when utilizing the strength reduction method, in the process, the cohesion and internal friction angle were reduced by the bisection method usually. To reduce the iterations and computing time, the multi-thread parallel computing technology was utilized in the strength reduction method. Through theoretical analysis and computing, the minimum iteration equations of the bisection method and multi-thread parallel computing method were derived with the range (Kmin, Kmax) of reduction parameter and the given error η determined. The FLAC3D command flow was embedded in Python scripts, subsequently, the multi-thread computing of FLAC3D software were performed by using the Python scripts. A specific computing example was conducted by 4 threads parallel computing and the bisection method, separately. From the study of this paper, it is concluded that the final results are the same, meanwhile, the iterations of 4 threads parallel method is 0.46 times less than that of the bisection method, and the time is 0.55. The computing example suggested that the multi-parallel computing technology can reduce iterations and computing time effectively.%强度折减法将边坡的安全系数定义为边坡刚好达到临界破坏状态时的折减系数,该过程一般采用二分法对边坡内聚力和内摩擦角进行迭代计算.为优化迭代计算次数以缩短计算时间,提出多线程并行运算强度折减法的计算策略,得到二分法和多线程并行运算在不同折减系数计算范围(Kmin,Kmax)、给定误差条件η下所需最小迭代次数的理论公式.将FLAC3D命令流内嵌于Python脚本中,通过Python脚本同时调用多个FLAC3D 应用程序,分别采用4线程并行运算和传统二分法对实例边坡进行强度折减法的计算.研究结果表明:在本文计算条件下,2种方法计算结果相同,采用4线程并行运算时迭代次数是二分法运算的0.46倍,其总的计算时间为二分法的0.55倍.采用多线程并行运算能有效减少迭代次数及总的计算时间.【期刊名称】《铁道科学与工程学报》【年(卷),期】2018(015)007【总页数】7页(P1693-1699)【关键词】Python脚本;FLAC3D;边坡;强度折减法;多线程并行运算【作者】瞿秀祥;林杭;王敏;张道勇【作者单位】中南大学资源与安全工程学院,湖南长沙 410083;中南大学资源与安全工程学院,湖南长沙 410083;中南大学资源与安全工程学院,湖南长沙410083;中机国际工程设计研究院有限公司,湖南长沙 410007【正文语种】中文【中图分类】TU457运用FLAC3D编制强度折减法计算程序在边坡稳定性分析中应用广泛，研究边坡强度折减法计算时间具有重要经济效益和实际意义。
基于 FLAC3D 强度折减法的土钉支护参数优化
基于 FLAC3D 强度折减法的土钉支护参数优化
赵洪月
【期刊名称】《中国新技术新产品》
【年(卷),期】2013(000)006
【摘要】基于强度折减法采用正交试验设计对高速公路边坡土钉支护参数进行优化，分析支护参数敏感性，结果表明：(1)土钉支护参数的敏感性从高到低依次为：垂直间距原长度原倾角原水平间距；(2)优化后的土钉支护参数为：土钉长度10m、倾角为15毅、垂直间距为3m 和土钉水平间距为1.2m.
【总页数】2页(P22-22,23)
【作者】赵洪月
【作者单位】高邮市交通建筑工程处，江苏高邮 225600
【正文语种】中文
【中图分类】TD3
【相关文献】
1.基于FLAC3D的复合土钉支护数值模拟分析 [J], 李晶晶;程祖锋;耿立立;张蕊蕊
2.基于强度折减法的基坑土钉支护稳定性分析 [J], 曾酉源;唐延贵
3.基于强度折减的土钉支护基坑稳定性评价 [J], 王娟;熊传祥;黄伟达
4.基于强度折减法的土钉支护结构稳定性分析 [J], 古海东;杨敏
5.基于Flac3D的土钉支护应力分析 [J], 丰土根;单亚;杨景森;杨才;胡康俊
因版权原因，仅展示原文概要，查看原文内容请购买。
用FLAC～3D软件分析高速公路边坡的稳定性
系数为11054。
(2)受力分析。
图3所示为平行线路的高边坡滑体最低部一条块的受力分析示意图与力多边形。
(3)可靠性分析。
根据所给定参数的标准偏差值(S 1D )对滑面的稳定性系数进行了破坏概率分析,由图4可以得出高边坡的破坏失稳概率为80%左右(F s ≤1120)。
4 结语 太平湖大桥铜陵端由拱座基坑开挖引起的高边坡是顺倾结构边坡,经过破坏模式分析和稳定性计算,评判为欠稳定边坡,破坏失稳概率为80%左右(F s ≤1120),故需对高边坡进行加固处理。
参考文献[1]徐邦栋1滑坡分析与防治1北京:中国铁道出版社,20011[2]蒋建平,等1优势结构面理论在岩土工程中的应用1工程地质学报,2000,8(4):438～4411[3]The sl ope /w document .GEO -S LOPE I nternati onal L td .2000收稿日期:2005-01-113西南交通大学峨眉校区科技基金资助项目:05D13021尹紫红,男,副教授。
用F LAC 3D软件分析高速公路边坡的稳定性3 尹紫红 周志林 姚令侃(西南交通大学峨眉校区 四川峨眉 614202)(西南交通大学土木学院)摘 要 F LAC 3D 是一个先进的岩土计算软件。
介绍F LAC 3D的基本原理及其使用方法,并对西(西昌)攀(攀枝花)高速公路K158+400～+580边坡开挖过程中应力、位移和稳定性进行分析。
通过计算结果,对边坡工程设计和防治提供了依据和指导意见。
关键词 F LAC 3D 位移变形 应力变形 稳定分析1 前言西攀高速公路K158+400～+580边坡位于四川盐边县新民村豆芽菁冲沟右侧,边坡高20m ,属中山区斜坡剥蚀地貌,台阶斜坡地形,倾向东南,坡度较平缓,坡角约30°。
场地揭露地层由新至老为:耕填土、低液限粘土、昔格达组粉砂岩夹薄层泥岩、泥岩夹粉砂岩。
边坡紧靠豆芽菁冲沟,枯水季节沟内无水,雨季汇水最深112m 。
用强度折减法和FLAC 3D计算边坡的安全系数
Ｃ
。
＝ｒ （ ａｃ
）
（） ２
其 中 ｃ 为实际 的抗剪强度参数 ， 、
为折减 系
收稿 日期（ ｅｅ ｅ ａ ）：０６－ ３－ ８ 改回日期（ ｃｅｔｄ ２０ ０ Ｒ ｃｉ ｄｄｔ ２ ０ ０ ０ ； ｖ ｅ Ａ ｃｐｅ ）：０６－ ８—２ 。 ｏ
作者简 介（ ｉ ｒ ｈ ） 李新坡 （９ ８ ， 汉 ） 山东萃县人 ， 士， Ｂｏ ａ ｙ ： ｇｐ １７ 一） 男（ ， 博 主要从 事建 筑、 环境岩土工程研究。 ［ ｉ ｉ ｏ １７ Ｌ Ｘｎ （９ ８一） ｍ ｌ Ｐ ． ， ｅ ｐ ，ａ ｈＤ Ｒ － ｅ ｓａ ｈ ｉ ｃ ｏ ： ｅｔ ｈ ｉ ｌ ｎｉ ｅ ｎ Ｉ ｒ ｉ ｃ ｒ ｅｒ ｒ ｔ ｎ Ｇ ｏ ｅｎｅ ｇ ｅｒ ｇ １ ｃ ｔ ｔ ｅ＆ｅｖ ｏｍ ｎ． ｍ ｉ Ｓ ｐｒ ＠ １６ ｃ ］ ｃ ｄ ｅｉ ｅ ａ ｅ ｎ ｉ Ｏ ｈ ｅ ｕ ａ ｎｉ ｎ ｅｔ Ｅ— ａ ：ｍ ｏ ｒ ｌ ｅ ２ ．∞
中图分类号 ：＂．５ １Ｉ ７ ７４ 文献标识码 ：Ａ
２ 应变与抗剪强度相关性的存在。 ． 随着计算机软硬件 和数值计算软件 的发展 ， 采 在边坡稳定分析 中， 对安全 系数 Ｆ的传统上 的 用数值模拟技术对边坡进行稳定性分析计算 ， 成为 当前边坡稳定性计算 的新方向。但是一般的数值分 定义为： 土体实际抗剪强度与使边坡保持稳定所要 ｊ ｕｅ ｎ Ｆ是 析 软件不 能直 接 计 算 出 “ 全 系数 ” 安 。通 常 的方 法 求的最小抗剪 强度 的比值 。Ｄ ｎａ 指 出， 是取某一参数 的不同的值进行试算 ， 取破坏时的参 使边坡达到临界破坏状态土体强度必须进行折减 的 数值与实际值的比值作为安全系数…。强度折减法 个系数。因此 ， 在用有限元或有 限差分程序计算 即是对土体或其他材料的强度进行折减降低使边坡 Ｆ ， 时 一个显然 的途径就是 降低土 的抗剪强度使边
浅谈用FLAC(3D)强度折减法计算边坡安全系数时特定因素对精度的影响
科
浅 谈 用 Ｆ ＡＣ Ｌ 如强 度 折 减 法 计 算 边 坡 安全 系数 时特 定 因素 对精 度 的影 响
黄 小 刚
（ 庆 交通 大 学 ， 庆 ４ ０ ７ ） 重 重 ００ ４
摘 要 ： 于 强度 折 减 法 真 正 兴起 才 １ 由 ０余 年 ， 对 于 大 家广 为 接 受 的 极 限平 衡 方 法 ， 少人 对 之 仍 持 怀 疑 态 度 。 本 文 将 主 要 讨 论 剪 胀 角 、 相 不 抗拉 强度 、 网格疏 密程度及边界 范围四方面因素对 Ｆ Ａ ｍ强度折减 法计 算结果的影响。 ＬＣ 关 键 词 ： Ｉ ，； 度 折 减 法 ； 度 影 响 ＦＡＣｏ强 精 剪 胀 角对 安 全 系数 的 影 响 １强 度折减 法 的基本原 理 ２ ４广义剪应变或广义塑性应变判断依据 ， 通 剪ｆ 。 ） 画出广义剪应变或广义塑性应 姊ｊ （ 目前数值分析方 法一般 只能得 出边坡的应 过计算边坡应变后， Ｌ ０４ ｌ ０ｌ Ｌ ∞ １１ ． ０ 安垒蘸鼓 力、 位移、 塑性区等， 无法得到边坡危险滑面以及相 变的等值线图。 当某一幅直的等值线从坡底贯通到 表 ２ 抗 拉 力对 安 全 系数 的 影 响 坡 顶时 ， 味着 边坡 达到极 限 状态 。 意 应的安全系数 。 ！ 模拟方法相结合 ， 在给定的评判指标下 ， 通过调整 的重 点 ，究竟 ！ 安全系数 ｌ 咕 １ ５ ｌ５ ０ ＿ １５ ０ ０ ． ｔ ５ Ｏ Ｌ 仿 １ ３ ． ｌ ０ 仍 ｌ 以哪个为标准 ， 采取哪个判据为好 ， 至今为止还没 折减系数确杰 嘲蝴 稳定 出 亍 分析 ， 求得边坡的最 小稳定安全系娄 。其基本原理是利用公式（ 、 形成统一， １ 来 ） 只能根据实际工程概况并综合运 调整岩土体的强度指标 ｃ ， ， 得到一组新的 ｃ ， 值， 然后对边坡 隐定性进行数值分析 ， 不断地 ３Ｆ ＡＣｏ Ｌ ３软件内嵌强度折减法 ＦＡ ｍ Ｌ Ｃ 软件内有—条 自 动查找安全系数的命 增加折减 系数， 以反复计算 ， 直至其达到ｌ 临界破坏 状态， 此时得到折减系数即为安全系数 Ｆ 令 一ｏｖ ｓ ｓｌ ｆ ，该命令就是对内聚力 ｃ 内摩擦角 ｅｏ 和 ｃ ＝ ｃ／Ｆ …， Ｆ （） １ 进 行折 减 ， ａ 一 （Ｉ ， ） ＝ｔ ｔＩ ｎ （ａ ） （） 安全系数 ， ２ 其实质就是强度折减法， 该命令通过内 式中Ｃ ，为折减后的粘接力 ， 为折减后 的 插逼近的办法来确定安全系数。 其计算过程如— ： ３ 计算特征时步数 先赋予内聚力一个 比 ． １ 磨擦角， 为折减系数。 ２边坡失稳破坏判断依据 较大的值， 使得边坡内部应力有较大的改变 ， 计算 随着计算机性能的提高以及各种软件 的开 出系统恢复平衡 昕需的时步数 ，就是特征时步数 发， 强度折减法受制于计算程序和计算精度 的局面 Ｎｒ 。 ３ 计算安全系数 Ｆ ． ２ ：对于一个给定的安全系 得到了根本性的改变，０ ２ 世纪 ９ ０年代成为岩土工 目 前碰到的困难是， 数 Ｆ运行—个 Ｎ 时步， ， ｒ 如果不平衡率小于 ｌ 则 ０ 现 用的 边坡 失稳判 剧 主 认为系统处于平衡状态；如果不平衡率大于 ｌ Ｏ， 要有 以下 几种 则执 下—个 ｒ ， 亍 步 如此循环 ， 直到不平衡率小于 ０ 结束循环。在不平衡率大于 ｌ 。的情况下， Ｏ 将 ２ 收敛性准则 ， ． １ 即有限元数值计算迭代过程 １ 。 是以收敛作为准则， 当计算收敛时边坡仍处于稳定 当前 Ｎ 时步范围内不平衡率的平均值与上一个 ｒ ｒ 状态， 计到不收敛时破坏 ， 此时的折减系数即为安 Ｎ 时步范围内不平衡率的平均值进行比较 ，若相 全 系 。 差小于 １ ％， ０ 认为系统处于不平衡状态（ 此时平均 该数值比相应的 ２ ２塑性区判据 ，强度折减法计算边坡稳定 不平衡率收敛于—个稳定的数值， 全 系数 与 剪 切 应 变增 量 云 图 时， 一般都是｝ 土体看成理想弹塑性材料 ， 巴 当土体 平衡状态下 的大 ，此时系统一定处于持续运动状 它们的磨擦角小得多， 通常在 ０ ２ ￣内变化 ， ｏ～ ０ 因 其位移便无限制的 态）则跳出当前循环， ， 利用新的不平衡状态下的 Ｆ 此 ， 胀角 剪 定 范 增大， 当整个边坡中的塑 陛区连贯时， 则此时就认 值重新开始循环计算 ； 如果相差大于 １％， ０ 继续执 围 的。 为 边坡 产生滑 动破 坏 。 行 Ｎ 时步循环计算 ， ｒ 直到遇到 以 下三种 清 况中任 ４ ２抗拉强度。 ２ 表 为采用不同的抗拉强度时 ， ２ 特征点位移判据 ， ３ 该判据认为当边坡达到 意—种 ： Ｅ ａ 述差另 于 １ ％ ｈ 、 Ｏ ； 已经执行了 ６ Ｎ 次 ｒ 虽然抗拉 极限状态时 ， 某些点的位移会趋于流动状态。但在 时步；不平衡率小于 １。 遇到前两种 睛 ｃ Ｏ。 况时都会 计算得到的安全系数。从表中可以看出， 艮 从 ｋ ａ １ｋ ａ但安全系数 选取特征点时应特别注意， 具体选用哪个点作为标 跳出当前循环 ， 利用新的不平衡状态下的 Ｆ值重新 强度的变化范围彳 大 ， ｌＰ ￣ ０ Ｐ ， 的变化去彳 ， ， 口艮Ｊ 这说明抗拉强度这一因素对用强度 、 准很重要 ， 在某些复杂情况下 ， 选取不同的点可能 开始循环计算 ， 遇到第三种 晴 况就结束循环。 循环结束后 的 Ｆ值就是求得 的边坡安全系 折减法计算边坡安全系数的影响是很小的。 会出现不同的安全系数， 目前大多是 以坡顶水平位 ４ ３建模网格疏密程度。图 ２为将基准方案中 移的突变作为判断依据。 数。 的网格加密一倍后 ， 计算得出的安全系数 、 剪切应 ４强度折减法精度影响因素探讨 变增 量云 图 （ ｉ 基准方 案 ） 图 为 。此 时安 全系数 为 尽管 强度 折减 法 只对 抗剪 强 度参 数 和进 行折 ．， ０ ０ ． 有一点差别 ； 剪切应变增 减， 但由于其基于数值模拟技术 ， 相对极 限平衡法 １ ４ 比基准方案的 １ ５ 方法而言， 要输 ＾更好的强度参数。由于在实际工 量云图的位置和分布范围却没什么变化。这说明，
FLAC3d程序使用手册pdf
FLAC3d程序使用手册pdf•FLAC3D程序概述•FLAC3D程序安装与配置•FLAC3D程序基本操作•FLAC3D程序高级功能应用目•FLAC3D程序输出结果处理与可视化•FLAC3D程序性能优化与调试技巧录01FLAC3D程序概述FLAC3D是Itasca公司开发的一款三维显式有限差分程序，用于模拟岩土和其他材料的力学行为。
该程序基于拉格朗日算法，适用于大变形和非线性问题。
FLAC3D具有强大的后处理功能，可以直观地展示模拟结果。
求解器该程序采用显式求解器，能够快速求解大规模问题。
FLAC3D 允许用户定义模型的初始应力场和位移场。
边界条件用户可以在模型中定义各种边界条件，如固定边界、自由边界等。
本构模型FLAC3D 提供了多种岩土本构模型，如摩尔-库仑模型、修正的剑桥模型等。
接触算法该程序可以模拟不同材料之间的接触行为，如土与结构物的接触。
FLAC3D 广泛应用于岩土工程领域，如边坡稳定性分析、基坑开挖模拟等。
岩土工程地质工程环境工程水利工程该程序可用于模拟地质工程问题，如隧道开挖、地下水资源评价等。
FLAC3D 可用于模拟环境工程问题，如垃圾填埋场稳定性分析、污染物运移模拟等。
该程序可用于水利工程领域，如大坝稳定性分析、水库渗漏模拟等。
FLAC3D 程序应用领域02FLAC3D程序安装与配置•获取安装程序：从官方网站或授权渠道下载FLAC3D程序的安装文件。
0203双击安装文件，启动安装向导。
阅读并同意软件许可协议。
安装步骤0101选择安装目录和组件。
02等待安装程序完成文件的复制和配置过程。
03验证安装：启动FLAC3D程序，检查是否能正常打开和运行。
确保计算机操作系统符合FLAC3D 程序的最低系统要求。
操作系统要求根据FLAC3D 程序的要求，安装必要的依赖库和运行时环境。
依赖库安装将FLAC3D 程序所在的目录添加到系统的环境变量中，以便在命令行中直接调用程序。
环境变量设置根据需要，修改FLAC3D 程序的配置文件，以调整程序的运行参数和性能。
FLAC3D强度折减法在排土场边坡稳定性分析中的应用_郝全明
露天采矿技术2013年第7期FLAC3D 是由美国Itasca 公司开发的一个三维有限差分程序，其原理是采用快速拉格朗日有限差分法。
目前，FLAC3D 软件已广泛应用于边坡稳定性评价、拱坝稳定分析、地下矿床开采、隧道工程、支护设计等多个领域[1]。
在计算边坡稳定性系数时，直接在FLAC3D 软件中输入solve 命令求解的方法，计算速度慢，所需要的运算时间长，不够直观。
本文则另辟他径，采用将FLAC3D 软件和强度折减理论相结合的方法来求解边坡稳定性系数。
1FLAC3D 软件简介FLAC3D 包含11种材料本构模型，分别是空单元模型、3种弹性模型（正交各向异性、各向同性、横向各向同性）、7种塑性模型（应变硬化/软化模型、摩尔-库仑模型、双线性应变硬化/软化多节理模型、Druker-Prager 模型、D-Y 模型、修正的剑桥模型、多节理模型）。
FLAC 程序采用快速拉格朗日方法,它基于显示差分法求解运动方程和动力方程,将计算区域划分为若干个单元。
FLAC3D 具有功能强大的自动三维网格生成器，内部定义了多种基本的单元形状，在生成比较复杂的三维网格时，可以利用FISH 语言来自定义单元形状，然后组合这些基本单元即可[2]。
2强度折减法的原理在采用强度折减理论分析边坡的稳定性时，不需要对滑动面的位置和形状做假定，强度折减法的原理是首先将边坡强度参数粘聚力c 和内摩擦角φ同时除以折减系数F ，得到一组新的c 和φ值，然后将新的c 和φ值作为一组新的材料参数输入计算，不断增加F 的值，直到边坡达到极限状态，发生失稳破坏，这时所对应的F 的值即为边坡的稳定性系数。
用公式表示如下：tan φ′=（tan φ）Fc ′=c /F式中：F 为折减系数；c 、φ分别为岩土边坡最初的粘聚力和内摩擦角；c ′、φ′分别为经过折减后的粘聚力和内摩擦角[3]。
3边坡发生失稳破坏的判据本文采用的强度折减法判断边坡是否失稳破坏的依据主要有：1）位移突变判据。
点安全系数与强度折减法对比
式（7）ቤተ መጻሕፍቲ ባይዱ，m＝0， n
图 1 空间应力状态及摩尔库伦强度包络线 Fig.1 Stress state and mohr columb strength envelope
参数 取值
C /kPa 3.00
φ
/ º 19.60
γ
/ kN·m 20.00
-3
E / kPa 1.00e4
ν
Bulk / kPa
监测点1 监测点2 监测点3 最危险滑动面
折减系数为0.994时的塑性区图
Fs=0.994
折减系数为0.994 时的剪应变增量云图
图 2 采用强度折减法分析的 ACADS 边坡考题得到的失稳临界状态及点安全系数 Fig. 2 The critical state and point safety by strength reduction method with the ACADS slope case
算安全系数，求得安全系数为 0.99。该命令执行方 便，但安全系数只能精确到小数点后两位，且计算 时间较长，不适用于复杂边坡的稳定计算中。本文 结合 Flac3D 软件编制了强度折减 Fish 命令流，折减 系数增量为 0.001， 计算效率较高， 适用于复杂边坡 稳定性计算中。当折减系数为 0.984 时，边坡顶部 和中部的监测点合位移增量与强度折减系数增量之 比出现了第一次递增（见图 2） ，此时滑面位置的塑
K P ( Fs 2)
K P ( Fs 2)
AC AB
AE cos AB
（10）
[2C 1 3 tan ]cos （11） 1 3
n l 2 1 m 2 2 n 2 3