最新FLAC3D数值模拟上机报告
FLAC_3D的锚杆拉拔数值模拟试验_江文武
图 3 网格剖分图
szz
Z sxx X
锚杆
X Y
sxx
沿锚杆轴 向施加固 定的速度 v
szz 7.5m
10 m 5 m 限制 Y 方向的位移
图 4 锚杆拉拔数值模型示意图
为了模拟锚杆拉拔过程中的影响因素，即影 响锚杆锚固的效应的因素：1） 模拟了在同样的外 部条件下，唯有浆体的摩擦角（ φg = 00 ，100 ，200 ， 300 ，400 ） 不同的条件作用下沿着锚杆轴向、径向 锚杆的应力与应变的分布规律以及锚杆的锚固 力、浆体界面上的剪应力分布特征；2） 模拟了在 同样的外部条件下，唯有浆体有效围压（ σm = 0， 2，4，6，8 MPa） 不同的条件作用下沿着锚杆轴向、 径向锚杆的应力与应变的分布规律以及锚杆的锚 固力、浆体界面上的剪应力分布特征；同时还模拟 了锚杆在拉拔过程中，锚索与岩体间的界面发生 剪切屈服、产生滑动直至拉拔破坏具体过程.
·130·
哈尔滨工业大学学报
第 41 卷
变形和强度起着重要的作用［1 ～ 4］. 加锚岩体的数 值模拟方法大都还是基于有限元法，但一般都过 低估计锚固效果. 然而 FLAC － 3D 即三维快速拉 格朗日分析方法的出现，又为锚杆在岩体锚固机 理提供了新的机遇. 本文就锚固体的摩擦角、有效 围压等对锚杆锚固性能的影响作了分析，对锚杆 拉拔过程中锚杆锚固失效的特点进行了探索，并 将现场试验与数值模拟计算进行了对比和分析.
3 数值模拟试验结果
通过多种方案的数值模拟试验可知图 5（ a） 是现场试验得到一系列的力与位移之间的曲线， 从图 5（ a） 中得知锚杆直径为 15. 2 mm 的锚杆锚 固力 = 17 t / m. 图 5（ b） 是根据现场的地质条件建 模后计算得到的锚杆所受力与位移之间的曲线， 图 5（b）中显示当锚杆自由端施加的力小于某一 值时，力与位移基本成正比关系，当力达到一定值 即锚 固 力 时，力 保 持 不 变，而 位 移 呈 无 限 增 大 趋 势，说明锚杆已经整体失稳，锚固作用失效，图 6
预应力锚索加固边坡的FLAC3D数值模拟分析
文献标识码 ： Ｃ
文章编号 ： １ ０ ０ ８— ３ ３ ８ ３ （ ２ ０ １ ３ ） ０ ３— ０ １ ０ ７— ０ ２
预应力锚索支护技 术是边坡 加固的先进技术之一 ， 在 国 内外工程 中得 到广泛 应用 ， 但其 理 论研 究 还相 对滞 后。因 此， 结合工程实际 ， 利用 Ｆ Ｌ Ａ Ｃ ３ Ｄ软件 ， 对边坡 的锚 固作 用效 果及边坡稳定性状态进行数值模拟分析 ， 对边坡加固工 程具 有重 要 的 现 实 意义 。
为０ ． ２— ０ ． ６ ＭＰ ａ 。
生产 厂 家 提 供 的 出 厂 证 明 取 得 ， 具 体参 数 如 下： Ａ＝ １ ４ ０ ｍ ｍ ， Ｔ ＝ ２ ． ６ ｅ Ｎ ， Ｅ＝ １ ． ９ ５ ｅ “ Ｐ ａ ， 灌 浆体 的参数 经现场 抗拔试验 获得 ， 具体参数为 ｇ ｒ ＿ ｃ ｏ ｈ ＝ １ ０ ｅ ， ｇ ｒ ＿ ｋ＝ ２ ｅ 。 ２ ． ２ ． ３ 边界条件 以及初 始条件的设置 根据 以上建立 的模 型和 实际情 况限制模 型底部任 何方 向的位移 和右侧水平方 向的位移 ， 模型上部 与边坡部位为 自
＝２ ０ ｋ Ｎ ／ ｍ ， Ｃ＝６ ０ ｅ Ｐ ａ， ＝２ ０ 。 ， Ｋ ＝３ ． ５ ７ ｅ ７ Ｐ ａ
，
为保证贵惠高速公路区间各 路段 高边坡 的稳定 ， 坡 比采 用１ ： ０ ． ５～ １ ： １ ， 坡高 １ ０ ｍ一级 ， 采用框架式锚杆 和框架预应
Ｃ＝
２ ． ０ ８ ｅ Ｐ ａ， 锚索钢绞线 的横截 面积 、 抗 拉强度 、 弹 性模 量从
力锚 索联 合支 护 等 处 治 措施 ， 锚杆长 ６ ． ０ ｍ， 锚 索长 ６— ２ ６ ｍ， 间距 ５ ｍ ｘ ４ ｍ。预应力 锚索采用无 粘结 钢绞线 Ａ Ｓ Ｔ － Ｍ Ａ ４ １ ６—８ ７ ａ 标准 ２ ７ ０级 ７ 中１ ５ ． ２ ４ ａ ｒ ｍ。锚 固段长度 ８ ｍ， 钻 孔孔径 ‘ ｐ １ ３ ０ ｍｍ， 锚索孑 Ｌ 内 自孔底 一次性 压满水 泥浆 ， 注浆 压力为 ０ ． ３ ５— ０ ． ６ Ｍ Ｐ ａ 。锚 索 自由段采用 防护油及 塑料 管 隔离 ， 每束锚索设计施加张拉力 ８ ５ ０ ｋ Ｎ 。锚杆材料采用 ２ ５ 水泥砂浆 锚 杆 ， 施 工时 下 倾 与 水 平 夹 角为 ３ ０ 。 ， 允 许 误 差 ±１ 。 ， 锚杆注浆 的水泥浆 强度必 须保证 ９３ ０ ＭＰ ａ ， 注浆压 力
FLAC^3D中锚杆支护的数值模拟研究综述
收稿日期：２０２００３?０９基金项目：贵州理工学院高层次人才科研启动经费支持项目（０２０３００１０１８０２９）；国家安监总局２０１８年安全生产重大事故防治关键技术科技项目（ｇｕｉｚｈｏｕ－０００５－２０１８ＡＱ）；国家留学基金资助项目（２０１９０８５２００１７）作者简介：杜学领（１９８６－），男，满族，河北承德人，博士，副教授，研究方向为煤岩体动力灾害机理及防治。
ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１００５－２７９８．２０２０．０９．００１ＦＬＡＣ３Ｄ中锚杆支护的数值模拟研究综述杜学领（贵州理工学院矿业工程学院，贵州贵阳 ５５０００３）摘 要：锚杆是目前煤矿、岩土等支护工程中常用的材料，在ＦＬＡＣ３Ｄ中可进行锚杆支护的多角度数值模拟研究，文章对ＦＬＡＣ３Ｄ锚杆支护的数值模拟研究进行总结和展望。
分析表明：目前ＦＬＡＣ３Ｄ中建立锚杆模型以Ｃａｂｌｅ、Ｐｉｌｅ两结构单元较多，可采用Ｂｅａｍ、Ｌｉｎｅｒ等构建金属钢、托盘、锚喷层等支护要素。
ＦＬＡＣ３Ｄ一般作为验证性手段，对实验、实践内容进行证明。
重点总结了锚杆及支护构件、预应力锚杆、煤矿巷道支护、岩土工程支护、特殊锚杆的实现等研究进展。
当前的研究应用中，存在复杂工程问题的动态还原能力相对较弱、时空与时步不对应、研究细节缺失及二次开发的非公开性、对理论研究的支撑作用相对薄弱等问题。
未来，ＦＬＡＣ３Ｄ依然作为重要的模拟手段应用于新型支护理论和支护技术的验证，并可在跨平台建模及跨平台研究、精细化建模与多因素耦合研究等方面取得新的突破。
关键词：ＦＬＡＣ３Ｄ；锚杆支护；数值模拟；巷道支护；文献综述中图分类号：ＴＤ３５３ 文献标识码：Ａ 文章编号：１００５２７９８（２０２０）０９?０００１?１５ＲｅｖｉｅｗｏｆＮｕｍｅｒｉｃａｌＳｉｍｕｌａｔｉｏｎｏｆＲｏｃｋＢｏｌｔＳｕｐｐｏｒｔｉｎＦＬＡＣ３ＤＤＵＸｕｅｌｉｎｇ（ＩｎｓｔｉｔｕｔｅｏｆＭｉｎｉｎｇＥｎｇｉｎｅｅｒｉｎｇ，ＧｕｉｚｈｏｕＩｎｓｔｉｔｕｔｅｏｆＴｅｃｈｎｏｌｏｇｙ，Ｇｕｉｙａｎｇ ５５０００３，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｔｈｅｂｏｌｔｉｓａｃｏｍｍｏｎｌｙｕｓｅｄｍａｔｅｒｉａｌｉｎｃｏａｌｍｉｎｅｓ，ｇｅｏｔｅｃｈｎｉｃａｌａｎｄｏｔｈｅｒｓｕｐｐｏｒｔｉｎｇｐｒｏｊｅｃｔｓ．ＴｈｅｍｕｌｔｉａｎｇｌｅｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｒｏｃｋｂｏｌｔｓｕｐｐｏｒｔｃａｎｂｅｃａｒｒｉｅｄｏｕｔｉｎＦＬＡＣ３Ｄ．Ｔｈｉｓｐａｐｅｒｓｕｍｍａｒｉｚｅｓａｎｄｆｏｒｅｃａｓｔｓｔｈｅｂｏｌｔｓｕｐｐｏｒｔｎｕｍｅｒｉｃａｌｓｉｍｕｌａ ｔｉｏｎｉｎＦＬＡＣ３Ｄ．ＴｈｅｒｅｓｅａｒｃｈｓｈｏｗｓｔｈａｔｃａｂｌｅａｎｄｐｉｌｅｓｔｒｕｃｔｕｒｅｅｌｅｍｅｎｔａｒｅｕｓｅｄｆｒｅｑｕｅｎｔｌｙｉｎｃｕｒｒｅｎｔｒｏｃｋｂｏｌｔｍｏｄｅｌｉｎＦＬＡＣ３Ｄ，ａｎｄｂｅａｍ，ｌｉｎｅｒｃａｎｂｅｕｓｅｄｔｏｂｕｉｌｄｅｌｅｍｅｎｔｓｓｕｃｈａｓｍｅｔａｌｓｔｅｅｌ，ｔｒａｙａｎｄａｎｃｈｏｒｓｐｒａｙｌａｙｅｒ．ＦＬＡＣ３Ｄｉｓｇｅｎｅｒａｌｌｙｕｓｅｄａｓａｖｅｒｉｆｉｃａｔｉｏｎｔｏｏｌｔｏｐｒｏｖｅｔｈｅｔｒｕｔｈｏｆｅｘｐｅｒｉｍｅｎｔａｌａｎｄｐｒａｃｔｉｃａｌｃｏｎｔｅｎｔ．Ｔｈｅｒｅｓｅａｒｃｈｐｒｏｇｒｅｓｓｏｆｂｏｌｔｓａｎｄｓｕｐｐｏｒｔｉｎｇｅｌｅｍｅｎｔｓ，ｐｒｅ?ｔｅｎｓｉｏｎｅｄａｎｃｈｏｒｓ，ｃｏａｌｍｉｎｅｒｏａｄｗａｙｓｕｐｐｏｒｔ，ｇｅｏｔｅｃｈｎｉｃａｌｅｎｇｉｎｅｅｒｉｎｇｓｕｐｐｏｒｔａｎｄｍｏｄｅｌｉｎｇｏｆｓｐｅｃｉａｌｂｏｌｔｓａｒｅｓｕｍｍａｒｉｚｅｄ．Ｉｎｔｈｅｃｕｒｒｅｎｔｒｅ ｓｅａｒｃｈａｐｐｌｉｃａｔｉｏｎｓ，ｔｈｅｒｅａｒｅｐｒｏｂｌｅｍｓｓｕｃｈａｓｒｅｌａｔｉｖｅｌｙｗｅａｋｄｙｎａｍｉｃｒｅｓｔｏｒａｔｉｏｎｏｆｃｏｍｐｌｅｘｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ，ｎｏｎ?ｃｏｒｒｅｓｐｏｎ ｄｅｎｃｅｂｅｔｗｅｅｎｒｅａｌｔｉｍｅａｎｄｓｔｅｐｓｉｎＦＬＡＣ３Ｄ，ｌａｃｋｏｆｒｅｓｅａｒｃｈｄｅｔａｉｌｓａｎｄｎｏｎ?ｄｉｓｃｌｏｓｕｒｅｏｆｓｅｃｏｎｄａｒｙａｐｐｌｉｃａｔｉｏｎ，ａｎｄｒｅｌａｔｉｖｅｌｙｗｅａｋｓｕｐｐｏｒｔｆｏｒｔｈｅｏｒｅｔｉｃａｌｒｅｓｅａｒｃｈ．Ｉｎｔｈｅｆｕｔｕｒｅ，ＦＬＡＣ３Ｄｓｔｉｌｌｃａｎｂｅｕｓｅｄａｓａｎｉｍｐｏｒｔａｎｔｓｉｍｕｌａｔｉｏｎｍｅｔｈｏｄｆｏｒｔｈｅｖｅｒｉｆｉｃａｔｉｏｎｏｆｎｅｗｓｕｐｐｏｒｔｔｈｅｏｒｙａｎｄｔｅｃｈｎｏｌｏｇｙ，ａｎｄｃａｎｍａｋｅｎｅｗｂｒｅａｋｔｈｒｏｕｇｈｓｉｎｃｒｏｓｓ?ｐｌａｔｆｏｒｍｍｏｄｅｌｉｎｇａｎｄｃｒｏｓｓ?ｐｌａｔｆｏｒｍｒｅｓｅａｒｃｈ，ａｎｄｒｅｆｉｎｅｄｍｏｄｅｌｉｎｇａｎｄｍｕｌｔｉｆａｃｔｏｒｃｏｕｐｌｉｎｇｒｅｓｅａｒｃｈ．Ｋｅｙｗｏｒｄｓ：ＦＬＡＣ３Ｄ；ｒｏｃｋｂｏｌｔｓｕｐｐｏｒｔ；ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ；ｒｏａｄｗａｙｓｕｐｐｏｒｔ；ｌｉｔｅｒａｔｕｒｅｒｅｖｉｅｗ 锚杆作为一种岩土锚固中的重要材料，被广泛应用于边坡、坝体、基坑、隧道、巷道等工程场合，锚杆的使用距今已有１００多年的历史［１］。
土木工程数值模拟(FLAC3D)课件第2-7章
第二章 网格划分
第二章 网格划分
Generate <关键字> zone 产生三维空间的单元体 surface 产生三维空间的面 point 在三维空间定义参考点以帮助单元体和面的生成 merge 使Gen zone产生的相邻网格合并连接在一起
2020/7/10
土木工程数值模拟（FLAC）
主要语句
条件语句 IF 条件表达式 [THEN] … [ELSE] … ENDIF
FISH中条件运算符没有“并”、“或”、“否”这样的符号
表达“1<aa<2”的条 件
if aa > 1.0 if aa < 2.0
执行语句
endif endif
主要语句
循环语句 LOOP var (exp1, exp2)
内部矩形巷道贴满单元体单元格 数6、12、8，体外环绕放射状网 格单元7
上机内容：直墙半圆拱
2020/7/10
土木工程数值模拟（FLAC）
10
第二章 网格划分
建立任何网格都要从两个方面考虑：一是重要区域精确解 所需要的单元体密度；二是网格边界定位对结果的影响。应 力、应变变化大的区域往往单元体密度大。
内部矩形巷道边长分别是3m 6m 4m, 单元格数size也是3、6、4
2020/7/10
土木工程数值模拟（FLAC）
9
第二章 网格划分
利用参数fill来生成需填充的网格
gen zone radbrick p0=(24,-20,0) & p1=(34,-20,0) & p2=(24,-10,0) & p3=(24,-20,10) & dimension 3 6 4 & size 6 12 8 7 & fill group inner
基于FLAC3D数值模拟的煤层群上行开采分析
基于 FLAC3D 数值模拟的煤层群上行开采分析摘要：针对近距离煤层群上行安全高效开采上部遗漏煤炭资源的问题，采用理论分析和FLAC3D数值模拟实验相结合的方法，分析了近距离煤层群上行开采的采动影响范围，模拟了煤层群下部7号煤层1377工作面和9号煤层1397工作面开采后对上部5号煤层1357工作面产生的采动影响。
通过现场工业性试验表明：在1357工作面回采时，在回风巷端头开始80m范围内倾角变化量较大，可达2°至3°，80m范围以外煤层倾角基本无变化量；回风顺槽的下沉量从停采线向切眼方向下沉量从0.2m逐渐增大到4.2m，而运输顺槽的下沉量基本未发生变化，实现上部遗漏煤炭资源的安全开采，获取良好的经济效益和社会效益。
关键词：近距离；煤层群；上行开采；采动影响；数值模拟实验中图分类号：TD 822文献标识码：A1工程概况钱家营煤矿位于河北省，井田煤系属于下统二叠系和上统石炭系，中奥陶统马家沟组石灰岩构成基底地层，厚度500m左右的煤系地层，煤系共有十几层含煤层，煤层总厚度达19.79m，含煤系数3.96%。
井田内含6层可采煤层，即煤5、煤7、煤8、煤9、煤11、煤12-1。
煤层平均厚度分别为1.42m、3.24m、1.25m、2.2m、0.78m、2.29m，煤层平均间距分别为32.56m、6.26m、10.69m、15.24m、62.09m。
目前该矿井正开采5号煤层的1357工作面，工作面范围内煤层厚度在0.4～2.4m之间，平均厚度1.4m，开切眼附近和工作面中部为薄煤层开采条件，煤层厚度小于1.3m，5号煤层整体结构简单，赋存较稳定；煤层倾角在4～10°之间，平均7°，煤层走向在N50°～70°之间。
5号煤层1357工作面倾斜下方距离32.56m为7煤层1377工作面，已回采完毕；下方距离64.75m为9煤层1397工作面，正在回采；下方12-1煤层暂无工程施工。
基于FLAC3D的深基坑支护三维数值模拟分析
Ｖｏ Ｉ ． ２ ９ Ｎｏ ． １，Ｆ ｅ ｂ ．２ ０ １ ７
ｄ ｏ ｉ ：１ ０ ． ３ ９ ６ ９ ／ ｊ ． ｉ ｓ ｓ ｎ ． １ ６ ７ １ — ８ ７ ９ ８ ． ２ ０ １ ７ ． ０ １ ． ０ ０ ７
基于 Ｆ Ｌ Ａ Ｃ ３ Ｄ 的 深 基 坑 支 护 三 维 数 值 模 拟 分 析
ｆ ｏ ｕ ｎ ｄ ａ ｔ ｉ ｏ ｎ ｐ ｉ ｔ ｓ ｕ ｐ ｐｏ ｒ ｔ ｉ ｎ ｇ ｂ ａ ｓ ｅ ｄ ｏ ｎ ＦＬＡＣ３ Ｄ
Ｃ Ｈ ＥＮ Ｇ Ｚｅ ｈａ ｉ 。Ｙ Ｕ Ｚｈｅ ｎｓ ｈｕ ａ ｉ （ Ｓ ｃ ｈ ｏ ｏ ｌ ｏ ｆ Ｃｉ ｖ ｉ ｌ Ｅｎ ｇ ｉ ｎ ｅ ｅ ｒ ｉ ｎ ｇ ａ ｎ ｄ Ａｒ ｃ ｈ ｉ ｔ ｅ ｃ ｔ ｕ ｒ ｅ ，Ｚ ｈ ｅ ｊ ｉ ａ ｎ ｇ Ｕｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ ｏ ｆ Ｓ ｃ ｉ ｅ ｎ ｃ ｅ ａ ｎ ｄ Ｔｅ ｃ ｈ ｎ ｏ ｌ ｏ ｇ ｙ，Ｈａ ｎ ｇ ｚ ｈ ｏ ｕ ３ １ ０ ０ ２ ３，Ｚ ｈ ｅ ｊ ｉ ａ ｎ ｇ，Ｃｈ ｉ ｎ ａ）
中图分类号 ： ＴＵ ４ ７ ６ ． ４ 文 献 标 志 码 ：Ａ 文 章 编 号 ：１ ６ ７ １ — ８ ７ ９ ８ （ ２ ０ １ ７ ） ０ １ — ０ ０ ３ ７ — ０ ６
Ｔｈｒ ｅ ｅ — ｄｉ ｍｅ ｎ ｓ ｉ ｏ ｎａ ｌ ｎｕ ｍｅ ｒ ｉ ｃ ａ ｌ ｓ ｉ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ ａ ｎ ａ ｌ ｙ ｓ ｉ ｓ ｏ ｆ ｄ ｅ ｅ ｐ
程 泽海 ， 于 振 帅
（ 浙 江科 技 学 院 土 木 与 建 筑 工 程 学 院 ， 杭州 ３ １ ０ ０ ２ ３ ）
摘 要 ： 为 研 究 基 坑 不 同 支 护 方 式 对 围 护 结 构 变 形 及 稳 定 性 的影 响 ， 利用 Ｆ Ｉ Ａ Ｃ ３ Ｄ 三 维 快 速 拉 格 朗 日差 分 方 法 对某地铁深基坑分步开挖与支护进行数值模 拟 ， 并 对 两 种 支 护 方 案 进 行 对 比 分 析 。研 究 结 果 表 明 ： 地 下 连 续 墙 最 大 水平 位 移 出 现 在 墙 顶 ， 且 位 于 地 下 连 续 墙 长 度 方 向 的 中部 ； 在分步开挖 时 ， 第 一 步 开 挖 时 地 下 连 续 墙 的 位 移
FLAC3D在深基坑工程开挖中的数值模拟分析
随着基坑开挖深度 的增 加而加大 ， 基坑 壁 向坑 内的水平位 移变 化趋势仍 然是 中间部分最 大 ， 边角处 最 小， 而且基坑壁的长边 由于开挖 的范 围相对较大 ， 其变形量 相对 于短边也增 大 ， 这充分体 现 了基坑 开挖
过程 中的时空效应 ， 数值模拟计算结果可 以为工程设计提供指导 和参考 。 关键词 ： Ｆ Ｌ Ａ Ｃ ３ Ｄ ； 深 基坑 ； 位移 ； 数值模拟 中图分类 号 ： Ｔ Ｕ ４ ７ ０ ． ３ 文献标 识码 ： Ａ 文章编 号 ： １ ６ ７ ２ ＿ｌ １ ４ ４ （ ２ ０ １ ３ ） ０ ４ —０ ｏ １ ７ —０ ４
第 １ １ 卷第 ４期
２ ０ １ ３年 ８月
水 利与 建筑工 程学 报
ｏ ｕ ｒ ｎ ａ ｌ ｏ ｆ Ｗａ ｔ ｅ ｒ Ｒｅ ｓ ｏ ｕ ｒ ｃ ｅ ｓ ａ ｎ ｄ Ａ ｒ ｃ ｈ ｉ ｔ ｅ ｃ ｔ ｕ ｒ ａ ｌ Ｅ ｎ ｎ ｅ ｅ 血ｌ ｇ
Ｖｏ １ ． １ １ Ｎｏ． ４
Ａｐ ｐ ｌ ｉ ｃ ａ ｔ ｉ ｏ ｎ ｏ ｆ ＦＬＡＣ３ Ｄ ｉ ｎ Ｎｕ ｍｅ ｒ ｉ ｃ ａ ｌ Ｓ ｉ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ Ａｎ ａ ｌ ｙ ｓ ｉ ｓ ｆ ｏ ｒ Ｄｅ ｅ ｐ Ｆｏ ｕ ｎｄ ａ ｔ ｉ ｏ ｎ Ｐｉ ｔ Ｅｘ ａｖ ｃ ａ ｔ ｉ ｏ ｎ
ｔ ｒ ａ ｌ ｂ ｓｅ ａ ｍｅ ｎ ｔ ，ｔ ｈｅ ｕ ｐ ｌ ｉ ｆ ｔ ｇ ｒ ｏ ｗ ｓ ｔ ｏ ｈｅ ｔ ｌ ａ ｒ ｇ ｅ ｓ ｔ ｖ ｌ ａ ｕ ｅ ，ａ ｎ ｄ ｎ ｅ ｒ ａ ｈｅ ｔ ｆ ｏ ｕ ｎ ｄ ａ ｔ ｉ ｏ ｎ ｐ ｉ ｔ ｗ ｌｌ ａ ，ｉ ｔ ｉ ｓ ｓ ｍａ ｌｌ ｅ ｒ ．Ａｔ ｔ ｈ ｅ ｓ ａ ｍｅ ｔ ｉ ｍｅ ， ｔ ｈ ｅ Leabharlann Ａ ｕ ｇ．， ２０１ ３
基于FLAC3D的特大型露天边坡稳定性数值模拟分析
基于FLAC3D的特大型露天边坡稳定性数值模拟分析提纲：第一章：引言- 研究背景和意义- 国内外研究现状和进展- 研究目的和内容第二章：理论基础和数值模拟方法- 岩土力学基础理论- 边坡稳定性分析方法- FLAC3D软件介绍及使用方法第三章：数值模拟分析- 选取模拟模型及边界条件- 调试模型参数和边界条件- 分析模型的动态响应及应力变形分布第四章：模拟结果分析及讨论- 不同荷载及边坡角度条件下边坡的稳定性分析结果- 分析影响稳定性的因素及其重要性- 建议边坡的设计和加固方式第五章：结论与展望- 结论总结- 存在问题及展望未来研究方向- 对边坡设计和加固的意义和应用前景的评价注意：此提纲为中文版，如需翻译成英文可使用在线翻译工具进行翻译。
第一章：引言随着城市化的加速和工业经济的不断快速发展，大型的露天开采工程在当今社会中已成为常态。
然而，这些巨型露天工程也面临着一系列的问题，其中最重要的问题之一是边坡稳定性问题。
由于不同地形条件和巨大的冲击力，这种问题极其棘手，需要进行彻底和全面的研究。
边坡稳定性数值模拟分析是一种非常重要的研究方法，可以帮助工程师理解边坡的工程行为和各种负荷受力情况。
在此过程中，FLAC3D软件已经得到广泛的应用，它可以通过数值计算法来模拟实际的边坡开挖和加固过程，有效预测边坡的稳定性情况。
在分析边坡稳定性的过程中， FLCA3D模拟技术已经成为一种有效和可靠的工具。
本文旨在通过FLAC3D软件，对大型露天边坡的稳定性进行数值模拟分析，并通过实验结果来探讨边坡稳定性的各种因素和影响，以此作为改进边坡设计和加固方案的依据。
本论文内容分为四个章节，除此外还有引言和结论部分。
在本章中，我们将首先介绍大型露天开采工程中边坡稳定性问题的背景和意义。
其次，我们将对国内外关于边坡稳定性问题的研究进行回顾和评价。
最后，我们将阐明本研究的目的和内容。
首先，随着城市化的加速和工业化的迅速发展，巨型露天开采工程已经成为当今社会的常态。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
FLAC 3D数值模拟上机报告计算模型分别如图1、2、3所示，边坡倾角分别为30°、45°、60°，岩土体参数为： 密度ρ＝2500 kg/m 3, 弹性模量E ＝1×108 Pa ，泊松比μ＝0.3，抗拉强度σt ＝0.8×106 Pa ，内聚力C ＝4.2×104 Pa ，摩擦角φ＝17°，膨胀角Δ＝20°试用FLAC 3D 软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC 3D 软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。
附 换算公式：1 kN/m 3= 100 kg/m 3剪切弹性模量：881100.38510()2(1)2(10.3)E G Pa μ⨯===⨯+⨯+ 体积弹性模量：881100.83310()3(12)3(120.3)E K Pa μ⨯===⨯-⨯-⨯ 一 坡度为30°的情况4025.36604010030°图1 倾角为30°的边坡(单位：m)算例分析： 命令流： new;========================================================== 建立网格模型gen zone brick p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10gen zone brick p0 40 0 40 p1 100 0 40 p2 40 2 40 p3 74.64 0 60 p4 100 2 40 & p5 74.64 2 60 p6 100 0 60 p7 100 2 60 size 30 1 10;=========================================================;设置边界条件fix x y z range z -0.1 0.1fix x range x 99.9 100.1fix x range x -0.1 0.1fix y;======================;初始地应力的生成model elasprop density 2000 bulk 3e9 shear 1e9set gravity 0 0 -10solveini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0;===================================;安全系数求解model mohrprop density 2500.0 bulk 8.3E7 shear 3.8E7 coh 42000.0 tens 0.8E6 friction 17 dilation 20 solve fos file slope3dfos1.sav associated安全系数：最终计算边坡稳定性系数为1.453图1 网格剖分图图2 速度矢量图图3 速度等值线图图4 位移等值线图图5 剪应变增量云图二 坡度为45°的情况1004060404045°图2 倾角为45°的边坡(单位：m)算例分析：命令流： new;========================================================== 建立网格模型gen zone brick p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10gen zone brick p0 40 0 40 p1 100 0 40 p2 40 2 40 p3 60 0 60 p4 100 2 40 &p5 60 2 60 p6 100 0 60 p7 100 2 60 size 30 1 10;========================================================= ;设置边界条件fix x y z range z -0.1 0.1 fix x range x 99.9 100.1 fix x range x -0.1 0.1 fix y;====================== ;初始地应力的生成 model elasprop density 2000 bulk 3e9 shear 1e9 set gravity 0 0 -10 solveini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0;=================================== ;安全系数求解 model mohrprop density 2500.0 bulk 8.3E7 shear 3.8E7 coh 42000.0 tens 0.8E6 friction 17 dilation 20 solve fos file slope3dfos1.sav associated 安全系数：最终边坡的稳定性系数为1.14图1 网格剖分图 图2 速度矢量图图3 速度等值线图 图4 位移等值线图图5 剪应变增量云图三 坡度为60°的情况100406048.454060°图3 倾角为60°的边坡(单位：m)算例分析：命令流： new;========================================================== 建立网格模型gen zone brick p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10gen zone brick p0 40 0 40 p1 100 0 40 p2 40 2 40 p3 51.55 0 60 p4 100 2 40 & p5 51.55 2 60 p6 100 0 60 p7 100 2 60 size 30 1 10;========================================================= ;设置边界条件fix x y z range z -0.1 0.1 fix x range x 99.9 100.1 fix x range x -0.1 0.1 fix y;====================== ;初始地应力的生成 model elasprop density 2000 bulk 3e9 shear 1e9 set gravity 0 0 -10 solveini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0;=================================== ;安全系数求解 model mohrprop density 2500.0 bulk 8.3E7 shear 3.8E7 coh 42000.0 tens 0.8E6 friction 17 dilation 20 solve fos file slope3dfos1.sav associated 安全系数：最终边坡的稳定性系数为0.928图1 网格剖分图图2 速度矢量图图3 速度等值线图图4 位移等值线图图5 剪应变增量云图四 坡度为60°的边坡开挖情况开挖后坡面原始坡面345°45°100406048.454060°算例分析：命令流： new;========================================================== 建立网格模型gen zone brick p0 0 0 0 p1 100 0 0 p2 0 2 0 p3 0 0 40 size 50 1 10gen zone brick p0 40 0 40 p1 100 0 40 p2 40 2 40 p3 50 0 50 p4 100 2 40 p5 50 2 50 p6 100 &0 50 p7 100 2 50 size 30 1 10gen zone brick p0 53 0 50 p1 100 0 50 p2 53 2 50 p3 63 0 60 p4 100 2 50 p5 63 2 60 p6 100& 0 60 p7 100 2 60 size 15 1 10gen zone brick p0 45.77 0 50 p1 53 0 50 p2 45.77 2 50 p3 51.55 0 60 p4 53 2 50 p5 51.55 2 &60 p6 63 0 60 p7 63 2 60 size 15 1 10 group exc1gen zone wedge p0 40 0 40 p1 50 0 50 p2 40 2 40 p3 45.77 0 50 p4 50 2 50 p5 45.77 2 50 &size 30 1 10 group exc2group section1 range y 0 2 group exc1 group section2 range y 0 2 group exc2 attach face;========================================================= ;定义本构模型 mod elaspro density 2500 bulk 0.83e8 she 0.38e8;================================================================ ;设置边界条件fix x y z range z -.1 .1 fix x range x 99.9 100.1 fix x range x -0.1 0.1fix y;======================;设置重力加速度set gravity 0 0 -10.0;===================================;设定初始条件ini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0;===================================;开挖mod mohrpro density 2500 bulk 0.83e8 she 0.38e8 fric 17 coh 4.2e4 ten 0.8e6 dila 20mod null range group section1mod null range group section2solve fos file slope3dfos1.sav associated安全系数：最终边坡的稳定性系数为1.36图1网格剖分图图2 速度矢量图图3 速度等值线图 图4 位移等值线图图5剪应变增量云图。