FLAC3D滑坡模拟
FLAC 3D滑坡模拟
一、源程序
; Create Material Zones
gen zone brick size 5 5 5 &
p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,5) &
p4 (3,3,0) p5 (0,5,5) p6 (5,0,5) p7 (5,5,5)
gen zone brick size 5 5 5 p0 (0,0,5) edge 5.0
group Material
; Create Bin Zones
gen zone brick size 1 5 5 &
p0 (4,1,0) p1 add (3,0,0) p2 add (0,3,0) &
p3 add (2,0,5) p4 add (3,6,0) p5 add (2,5,5) &
p6 add (3,0,5) p7 add (3,6,5)
gen zone brick size 1 5 5 &
p0 (6,1,5) p1 add (1,0,0) p2 add (0,5,0) &
p3 add (0,0,5) p4 add (1,6,0) p5 add (0,5,5) &
p6 add (1,0,5) p7 add (1,6,5)
gen zone brick size 5 1 5 &
p0 (1,4,0) p1 add (3,0,0) p2 add (0,3,0) &
p3 add (0,2,5) p4 add (6,3,0) p5 add (0,3,5) &
p6 add (5,2,5) p7 add (6,3,5)
gen zone brick size 5 1 5 &
p0 (1,6,5) p1 add (5,0,0) p2 add (0,1,0) &
p3 add (0,0,5) p4 add (6,1,0) p5 add (0,1,5) &
p6 add (5,0,5) p7 add (6,1,5)
group Bin range group Material not
; Create named range synonyms
range name=Bin group Bin
range name=Material group Material
; Assign models to groups
model mohr range Material
model elas range Bin
; Create interface elements
int 1 face ran plane ori (4,0,0) nor (-5,0,2) dist 0.01 z (0,5) y (1,6) int 2 face ran plane ori (0,4,0) nor (0,-5,2) dist 0.01 z (0,5) x (1,6) int 1 face ran x 5.9 6.1 y 1 6 z 5 10
int 2 face ran x 1 6 y 5.9 6.1 z 5 10
int 1 maxedge 0.55
int 2 maxedge 0.55
; Move bin toward material
ini x add -1.0 range Bin
ini y add -1.0 range Bin
; Assign properties
prop shear 1e8 bulk 2e8 fric 30 range Material
prop shear 1e8 bulk 2e8 range Bin
ini den 2000
int 1 prop ks 2e9 kn 2e9 fric 15
int 2 prop ks 2e9 kn 2e9 fric 15
; Assign boundary conditions
fix x range x -0.1 0.1 any x 5.9 6.1 any ;在x=0和x=6的两个面上加约束fix y range y -0.1 0.1 any y 5.9 6.1 any
fix z range z -0.1 0.1 Bin
; Monitor histories
hist unbal
hist gp zdisp (6,6,10)
hist gp zdisp (0,0,10)
hist gp zdisp (0,0,0)
; Settings
set large
set grav 0 0 -10
; Cycling
step 4000
save bin.sav
二、分析结果
1、原结构网格图
2、z方向位移图
3、不平衡力变化图
4、hist gp zdisp (6,6,10)
zdisp (6,6,10)=-3.0928e-003
5、hist gp zdisp (0,0,10)
zdisp (0,0,10)=-1.1589e+000
6、hist gp zdisp (0,0,0)
zdisp (0,0,0)=-2.4947e+000。
呷爬滑坡稳定性的FLAC_3D数值模拟分析_张雪东
呷爬滑坡位于锦屏一级水电站三滩坝址上游
N20°E/SE∠40°，岩性为三叠系杂谷脑组深灰－
灰黑色变质粉砂质板岩夹变质细砂岩（T32-3z ）。 滑坡平面展布呈箕形，主滑方向为 N62°W，
其圈谷地貌保存完好。剪出口高程 1 655 m，与枯 水期河水位持平，后缘高程 2 120 m，前后缘高差 465 m。滑坡纵长约 900 m，平均宽约 500 m，面积 约为 0.45 km2。据钻孔揭露，滑体最厚达 60～ 100 m，平均厚度约 60 m，滑坡残体体积约 13×106 m3。
计算模型（图 1，2）范围, 顺河方向（ Z 轴方 向，指向上游为正）宽为 950 m，即顺河方向包括 整个滑坡的宽度；垂直河谷方向（X 轴方向，指向 滑坡内为正）长为 1 040 m，模型的最大高度（Y 轴 方向，垂直向上为正）1 108 m（海拔 2 308 m）。计 算时，仅考虑自重力的作用，不考虑构造应力场。 因此，除坡面（临空面）外，将模型四周边界、底 边界设置为单向约束边界[4]。
表 1 计算模型中的各种介质力学参数 Table 1 Mechanical parameters of materials in calculation model
模型中的组 （Group）
1 3 5 7 4 6 8 9 10 11 12
岩土体
T32-3z 滑带I 滑带II 滑带III 滑体I 滑体II 滑体III 冲积物 崩坡积物 小滑坡 泥石流堆积
该滑坡发育于雅砻江右岸反向坡，滑体下部基 岩为深灰－灰黑色粉砂质板岩，岩石较软弱，正常
岩体产状 N10～30°E/SE∠40°，地表岩体产状则 为 N15～25°E/NW∠20～30°，表明岩体已产生
强烈变形，岩体松驰。在谷坡岩体重力作用下，岩 体向临空面方向易产生倾倒变形，板岩岩质较软， 容易产生弯曲变形，岩体越向地表，弯曲变形越大， 随着倾倒变形的加剧，在最大弯曲部位，产生楔形 拉裂缝，发生弯曲－拉裂变形[1]。由于河谷下切及河 水对岸坡前缘的冲刷，弯曲－拉裂进一步加剧，当 岩体弯曲折断裂隙贯穿，形成滑带产生滑动。
基于Flac3D的某土质边坡稳定性分析
Value Engineering0引言滑坡作为一种常见的地质灾害在我国频繁发生并且分布广泛，不仅破坏基础设施，影响各工程的施工与建设，还阻碍国家的经济建设与发展进程，并对人民群众的生命财产安全造成严重的影响。
滑坡地质灾害如果处理不及时，将会造成一定的设施破坏、财产损失甚至是人员的伤亡[1-3]。
因此，研究分析滑坡地质灾害的成因和机理及滑坡的防治措施等一直是工程地质领域的热点问题。
瞬时暴雨或长期降雨等条件下诱发滑坡是土质边坡中最容易发生的类型[4]。
针对降雨条件下的边坡稳定性，众多学者主要从多个方面进行了相关的研究。
目前边坡稳定性分析评价的主要方法是极限平衡法和数值模拟法。
数值模拟分析方法包括有限单元法、有限差分法、离散单元法等。
李安润等[5]通过极限平衡法和有限元数值分析方法，对降雨条件下某堆积体边坡进行了稳定性分析，并提出了合理的防治措施。
Chang等[6]采用PFC数值模拟软件分析研究了某黄土滑坡在地表水入渗条件下边坡的失稳破坏过程。
回恒酉等[7]采用传统极限平衡理论的条分法和数值模拟的Flac3D方法进行对比分析，得出Flac3D 法分析条件更加完善，在理论上更加可靠，而条分法计算理论理想化，计算结果偏保守的结论。
李振江等[8]通过GeoStudio软件对暴雨工况下的南京某下蜀土滑坡进行了模拟分析，研究了暴雨条件下边坡的孔隙水压力和位移变化，并对应急治理措施进行了检验分析。
张树轩等[9]利用Flac3D模拟分析了甘肃天水红旗山黄土滑坡的稳定性，为潜在强震区地震滑坡的变形机理及防震减灾研究提供了可靠依据。
Flac3D数值模拟法是近年来比较流行的边坡稳定性计算分析方法，主要应用于土质滑坡，在岩质滑坡方面，相对应用较少。
本文以江苏西南部一土质边坡为研究对象，结合现场监测数据，采用Flac3D软件对边坡进行稳定性评价及变形破坏机理分析，为滑坡防治提供参考。
1工程概况1.1边坡基本特征江苏省西南部某一典型土质边坡现状如图1所示。
(完整word版)FLAC3D滑坡模拟
FLAC 3D滑坡模拟一、源程序; Create Material Zonesgen zone brick size 5 5 5 &p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,5) &p4 (3,3,0) p5 (0,5,5) p6 (5,0,5) p7 (5,5,5)gen zone brick size 5 5 5 p0 (0,0,5) edge 5.0group Material; Create Bin Zonesgen zone brick size 1 5 5 &p0 (4,1,0) p1 add (3,0,0) p2 add (0,3,0) &p3 add (2,0,5) p4 add (3,6,0) p5 add (2,5,5) &p6 add (3,0,5) p7 add (3,6,5)gen zone brick size 1 5 5 &p0 (6,1,5) p1 add (1,0,0) p2 add (0,5,0) &p3 add (0,0,5) p4 add (1,6,0) p5 add (0,5,5) &p6 add (1,0,5) p7 add (1,6,5)gen zone brick size 5 1 5 &p0 (1,4,0) p1 add (3,0,0) p2 add (0,3,0) &p3 add (0,2,5) p4 add (6,3,0) p5 add (0,3,5) &p6 add (5,2,5) p7 add (6,3,5)gen zone brick size 5 1 5 &p0 (1,6,5) p1 add (5,0,0) p2 add (0,1,0) &p3 add (0,0,5) p4 add (6,1,0) p5 add (0,1,5) &p6 add (5,0,5) p7 add (6,1,5)group Bin range group Material not; Create named range synonymsrange name=Bin group Binrange name=Material group Material; Assign models to groupsmodel mohr range Materialmodel elas range Bin; Create interface elementsint 1 face ran plane ori (4,0,0) nor (-5,0,2) dist 0.01 z (0,5) y (1,6) int 2 face ran plane ori (0,4,0) nor (0,-5,2) dist 0.01 z (0,5) x (1,6) int 1 face ran x 5.9 6.1 y 1 6 z 5 10int 2 face ran x 1 6 y 5.9 6.1 z 5 10int 1 maxedge 0.55int 2 maxedge 0.55; Move bin toward materialini x add -1.0 range Binini y add -1.0 range Bin; Assign propertiesprop shear 1e8 bulk 2e8 fric 30 range Materialprop shear 1e8 bulk 2e8 range Binini den 2000int 1 prop ks 2e9 kn 2e9 fric 15int 2 prop ks 2e9 kn 2e9 fric 15; Assign boundary conditionsfix x range x -0.1 0.1 any x 5.9 6.1 any ;在x=0和x=6的两个面上加约束fix y range y -0.1 0.1 any y 5.9 6.1 anyfix z range z -0.1 0.1 Bin; Monitor historieshist unbalhist gp zdisp (6,6,10)hist gp zdisp (0,0,10)hist gp zdisp (0,0,0); Settingsset largeset grav 0 0 -10; Cyclingstep 4000save bin.sav二、分析结果1、原结构网格图2、z方向位移图3、不平衡力变化图4、hist gp zdisp (6,6,10)zdisp (6,6,10)=-3.0928e-0035、hist gp zdisp (0,0,10)zdisp (0,0,10)=-1.1589e+0006、hist gp zdisp (0,0,0)zdisp (0,0,0)=-2.4947e+000。
FLAC3D在某古滑坡稳定性分析中的应用
文章编号 ： １ ６ ７ ３—９ ４ ６ ９ （ ２ ０ １ ３ ） ０ ４— ０ ０ ６ ６— ０ ４
ｄ ｏ ｉ ： ｄ ｏ ｉ ： １ ０ ． ３ ９ ６ ９ ／ ｊ ． ｊ ｓ ｓ ｎ ． １ ６ ７ ３— ９ ４ ６ ９ ． ２ ０ １ ３ ． ０ ４ 滑坡 稳 定性 分 析 中 的应 用
许 容 ， 张发 明 ， 张加 家 ， 眭敏磊
（ 河海大学 地球科学与Ｔ程学 院 ， 江苏 南京 ２ １ ０ ０ ９ ８ ）
摘要： 以 西南 某水 电站 左岸 古 滑坡 为研 究对 象， 利用 Ｆ Ｌ Ａ Ｃ ３ Ｄ通 过 强度 折 减 法 对 古 滑坡 堆积 体 变形和 稳 定性情 况进 行 了三 维数 值 仿真 。研 究 结 果表 明 ， 滑坡 体在 天 然状 态 下稳 定性 良好 ， 这
文献标 识码 ： Ａ
Ａｐ ｐ ｌ ｉ ｃ ａ ｔ ｉ ｏ ｎ ｏ ｆ Ｆ Ｌ ＡＣ３ Ｄ ｉ ｎ ｓ ｔ ａ ｂ ｉ ｌ ｉ ｔ ｙ ａ ｎ ａ ｌ ｙ ｓ ｉ ｓ ｏ ｆ ｔ ｈ ｅ ａ ｎ ｃ ｉ ｅ ｎ ｔ ｌ ａ ｎ ｄ ｓ ｌ ｉ ｄ ｅ ｂ ｏ ｄ ｙ
ａ ｓ ｔ ｈ ｅ ｒ ｅ ｓ ｅ ａ ｒ ｃ ｈ ｏ ｂ ｊ ｅ ｃ ｔ ． Ａ ３ Ｄ ｎ ｕ ｍｅ ｒ ｉ ｃ ａ ｌ ｓ ｉ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ ｍｏ ｄ ｅ ｌ ｉ ｓ ｄ ｅ ｖ ｅ ｌ ｏ ｐ ｅ ｄ ｔ ｏ ａ ｎ ａ ｌ ｙ ｚ ｅ ｔ ｈ ｅ ｄ ｅ ｆ ｏ ｒ ｍ ａ ｔ ｉ ｏ ｎ ｆ ｅ ａ —
ｔ ｕ ｒ ｅ ｓ ａ ｎ ｄ ｓ ｔ ａ ｂ ｉ ｌ ｉ ｔ ｙ ｏ ｆ ｔ ｈ ｅ ａ ｃ ｃ ｕ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ ｂ ｏ ｄ ｙ ｏ ｎ ｔ ｈｅ ｔ ｈｅ ｏ ｒ ｙ ｏ ｆ Ｔｈ ｅ ｓ ｔ ｒ ｅ ｎ ｇ ｔ ｈ ｒ ｅ ｄ ｕｃ ｔ ｉ ｏ ｎ ｍｅ ｔ ｈ ｏ ｄ ｂｙ Ｆ ＬＡＣ３ Ｄ．Ｔｈ ｅ ｒ ｅ ｓ ｅ ａ ｒ ｃ ｈ ｒ ｅ ｓ ｕｌ ｔ ｓ ｓ ｈ ｏ ｗ ｔ ｈ ａ ｔ ｔ ｈ ｅ ａ ｃ ｃ ｕｍｕ ｌ ａ ｔ ｉ ｏ ｎ ｂ ｏ ｄ ｙ ｉ ｓ ｓ ｔ ａ ｂ ｌ ｅ ｕ ｎ ｄ ｅ ｒ ｎ ａ ｔ ｕ ｒ ａ ｌ ｃ ｏ ｎ ｄ ｉ ｔ ｉ ｏ ｎ ｓ ． Ｔｈ ｅ ｒ ｅ ｓ ｕｌ ｔ ｓ ａ ｒ ｅ ａ ｌ ｓ ｏ ｃ ｏ ｎ ｓ ｉ ｓ ｔ ｅ ｎ ｔ ｗｉ ｔ ｈ ｔ ｈｅ ｓ ｐ ｏ ｔ ｇ ｅ ｏ ｌ ｏ ｇ ｉ ｃ ａ ｌ ｓ ｕ ｒ ｖ ｅ ｙ ｓ ｗｈ ｉ ｃ ｈ ｐ ｒ ｏ ｖ ｅ ｔ ｈｅ ａ ｐ ｐｒ ｏ ｐ ｒ ｉ ａ ｔ ｅ ｎ ｅ ｓ ｓ ｏ ｆ ｐ ａ － ｒ ａ ｍｅ ｔ ｅ ｒ ｓ ｃ ｈｏ ｏ ｓ ｉ ｎｇ ａ ｎｄ ｔ ｈ ｅ ｒ ｅ ｌ ｉ ａ ｂ ｉ ｌ ｉ ｔ ｙ ｏ ｆ ３ Ｄ ｃ ａ ｌ ｃ ｕ ｌ ａ ｔ ｉ ｏ ｎ ｍｏ ｄ ｅ １ ． Ｋｅ ｙ ｗｏｒ ｄｓ： ＦＬ ＡＣ３ Ｄ ；ｓ ｔ ｒ ｅ ｎ ｇ ｔ ｈ ｒ ｅ ｄ ｕ ｃ ｔ ｉ ｏ ｎ ｍｅ ｔ ｈ ｏ ｄ；ｎ ｕ ｍｅ ｒ ｉ ｃ ａ ｌ ｓ ｉ ｍｕ ｌ ａ ｔ ｉ ｏ ｎ；ｓ ｔ ａ ｂ ｉ ｌ ｉ ｔ ｙ
基于FLAC_3D_软件的露天矿边坡稳定性分析与研究
图 3 模型 20 竖直方向应力图
图 4 模型 20 水平方向应力图
图 2 模型 20 的最大不平衡力收敛曲线
［2］陈高峰，卢应发.边坡稳定性影响因素组成成分分析［J］. 金属矿山，2008.
［3］王浩.边坡稳定性影响因素及评价方法［J］.黑龙江交通 科技，2011.
［4］Matsui T，San K C.Finite element slope stability analysis by shear strength reduction technique［J］. Soils and Foun－ dations，1992，32（1） .
矿技术，2011. ［9］E vert Hoek John Bray. Rock Slope Engineering. London:
Traps.Institution of Mining and Metallurgy，1981. ［10］Dawson E M ，Roth W H ，Drescher A.Slope sTability
采矿工程
露天采矿技术 20天矿边坡稳定性分析与研究
白润才，王志鹏，王喜贤
（辽宁工程技术大学矿业学院，辽宁 阜新 123000）
摘 要：采用 FLAC3D 对露天矿边坡稳定性的数值模拟，阐明了边坡岩体应力场的演变规律，揭示了研究
区非工作帮边坡的滑坡机理，确定了滑坡模式。结果表明，滑坡模式为沿“圆弧滑面”滑动。基于 FLAC3D 的强
由位移图来看，非工作帮的边坡滑移主要在边坡 （下转第 5 页）
FLAC3D在滑坡稳定性分析评价中的应用
FLAC3D在滑坡稳定性分析评价中的应用赵忠海【摘要】Taking the Gaozhi Village landslide as an example, this paper discussed the application of FLAC3D in landslide stability analysis and evaluation. The author ifrst introduced the basic principle and main process of FLAC3D, and then simulated the Gaozhi Village landslide stability under natural state and heavy rain state by Mohr-Coulomb of FLAC3D. The simulation result intuitively showed the distribution characteristics of plastic deformation and horizontal displacement of the slope in the above two conditions, and clearly showed the potential sliding surface inside the slope. According to the result of numerical simulation, the article analyzed and evaluated the landslide stability, studied and discussed the deformation and failure forms of the landslide. By analyzing the plastic deformation and displacement of slope, the author considered that the landslide is basically stable under the natural state and is unstable and likely to slide and lose its stability in the local area under the continuous high intensity rainfall state. This is consistent with the investigation and analysis. This result shows it is feasible and more simple and more intuitive and convenient to analyze and evaluate the stability of landslide by FLAC3D. Compared with the traditional limit equilibrium method, FLAC3D has the incomparable advantage.%本文以高治村滑坡为例，论述了FLAC3D在滑坡稳定性分析评价中的应用。
基于FLAC(3D)的石灰溪滑坡稳定性数值模拟分析
基于FLAC(3D)的石灰溪滑坡稳定性数值模拟分析徐智【摘要】滑坡稳定性分析是评价滑坡稳定状态及有针对性选取治理措施的关键.目前有关滑坡稳定性的分析方法大多建立在基于刚体极限平衡理论的基础上,而工程实践表明,滑坡失稳很多时候是由于滑体变形过大所致.为此,本文以印江县峨岭镇石灰溪滑坡为研究对象,采用基于FLAC(3D)的数值模拟方法,从变形角度对该滑坡堆积体的稳定性及其发展趋势进行了分析.研究表明,滑坡中后部为启动段,中前部为锁固段,降雨状态下滑带塑性区贯通率达85%,滑坡整体处于加速变形阶段,需及时进行防治.【期刊名称】《科技风》【年(卷),期】2011(000)005【总页数】2页(P132-133)【关键词】FLAC(3D);滑坡稳定性;数值模拟;稳定性分析【作者】徐智【作者单位】贵州省地矿局一〇三地质大队,贵州铜仁,554300【正文语种】中文滑坡稳定性分析是在滑坡勘查，即地质调查、地表勘测、钻孔取芯分析、地下水检测、室内外试验基础上,对滑动面以上滑体运用数学、力学分析手段来评价滑体稳定状态。
目前有关滑坡稳定性的评价方法大多采用刚体极限平衡理论，难以动态反映滑坡堆积体在外界诱发因素作用下应力场的改变以及由其引起的变形情况，而大量的工程实践表明，滑坡失稳事件很多时候是由于滑体变形过大所导。
为此，有必要采用数值模拟方法对滑坡堆积体在外界诱发因素作用下的应力场和应变场进行深入研究，进而从变形的角度对滑坡堆积体的稳定性状况及发展趋势作出判断。
石灰溪滑坡位于贵州省印江县峨岭石灰溪村所在斜坡体上，平面形态呈上宽下窄的“锁口形”，主滑方向300°。
滑体纵向长度约540m，横向宽度约155m，面积约8.10×104m2，钻孔揭露滑体平均厚度约32m，滑体体积约260×104m3，属大型堆积体滑坡，如图1所示。
滑坡前缘剪出口标高约420m，后缘标高约630m，相对高差约210m。
FLAC-3D在黄土滑坡变形特征分析中的应用
FLAC-3D在黄土滑坡变形特征分析中的应用杜文举，景淑媛(四川建筑职业技术学院，四川德阳618000)摘要：暴雨是影响黄土滑坡稳定性的重要因素，本文以陕西-黄土滑坡为例，利用有限差分法FLAC-3D软件对在自然和暴雨条件下黄土滑坡的变形进行数值模拟分析。
根据数值模拟结果，得到滑坡体在自然状态下和暴雨状态下的变形特征和塑性分布特征，确定了滑坡体变形特征与滑坡稳定性间的联系，为黄土滑坡的稳定性评价提供必要的理论依据。
关键词：黄土滑坡；FLAC-3D;变形特征中图分类号：P642.22文献标识码：A文章编号：1006-0995(2020)04-0638-04DOI:10.3969/j.issn.l006-0995.2020.04.0231引言滑坡是一种较普遍的地质灾害，给人们生命财产带来了重大威胁和损失叫黄土滑坡是黄土地区特有常见地质灾害，全国约有1/3的滑坡灾害发生于黄土地区；以其多发性、危害性给生态环境和社会发展造成了严重的影响阴。
齐县余家坪乡新寨村一湿陷性黄土滑坡，牡体重力作用及WP秀延河侧蚀作用下形成忙滑坡西侧主滑方向314。
，东侧主滑方向354。
，整体上呈后缘缓，前缘陡，中部较缓，滑体前缘直达秀延河沟底，滑坡后缘高程2006-2008m,前缘高程在1992~1995m之间，相对高差12~13m,平均坡度15°-25°,该滑坡有下滑趋势，降雨将诱发滑坡的发生％对其滑体前缘的西气东输管道工程造成巨大影响和危害，严重威胁管道输气的安全旳。
滑体是由软塑状-流塑状黄土组成，上部黄土呈湿-很湿状态，中下部呈很湿-饱和状态，含水量明显高于上部，根据前、后缘的地表滑坡特征、地表张拉裂缝、钻孔所揭露滑动面位置高程(m)(黄海高程系统)滑坡堆积土断层角砾滑动面例三叠系砂泥岩马兰黄土图13-3,滑坡剖面图202520202015201020052000199519901985198019751970以及物探测试成果，判定滑动面剖面形状为弧形。
FLAC3D软件建立边坡的三维数值模型分析
FLAC3D软件建立边坡的三维数值模型分析结合FLAC3D软件的优点，以某公路工程的边坡为例，对FLAC3D软件在建立边坡三维数值模型中的应用进行了分析和探讨。
标签：FLAC3D软件边坡三维数值模型0前言对于公路工程而言，边坡的稳定性直接影响着工程施工的顺利进行，影响着整个工程的施工质量，在工程中的作用是十分巨大的。
影响边坡稳定性的因素是多种多样的，运用FLAC3D软件，结合相应的岩土勘察参数，可以建立边坡的三维数值模型，从而方便对边坡的应力场分布规律以及最大不平衡力的收敛情况进行分析，以实现对边坡的加固。
1 FLAC3D软件概述FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
FLAC3D采用了显式拉格朗日算法以及混合离散分区技术，可以非常准确地对材料的塑性破坏和流动进行模拟。
由于不需要形成刚度矩阵，可以在较小的内存空间中，求解大氛围的三维问题。
FLAC3D的优点包括以下几个方面：（1）混合离散法的应用，相比于有限元法中常用的离散集成法更加准确，更加合理；（2）采用动态运动方程实现对于静态系统的模拟，在模拟物理上的不稳定过程不存在数值上的障碍；（3）采用了“显式解”方案。
因此，显式解方案对非线性的应力-应变关系的求解所花费的时间，几乎与线性本构关系相同，而隐式求解方案将会花费较长的时间求解非线性问题。
而且，它没有必要存储刚度矩阵，这就意味着；采用中等容量的内存可以求解多单元结构；模拟大变形问题并不比小变形问题多消耗更多的计算时间，因为没有任何刚度矩阵要被修改。
2 FLAC3D软件建立边坡的三维数值模型2.1工程概况某城乡高速通道全长21km，从山林地区穿行，与该地区的一条河流并行，公路整体边坡呈岩质，部分路段为土质边坡。
受公路自身承载力等因素的影响，边坡工程出现了风化、滑坡等病害，严重影响了公路的行车安全，需要引起相关人员的重视和研究。
基于CAD与ANSYS的FLAC_3D_边坡模拟分析(1)
邮局订阅号：82-946360元/年技术创新软件时空《PLC 技术应用200例》您的论文得到两院院士关注基于CAD 与ANSYS 的FLAC 3D 边坡模拟分析The simulation analysis of slope by FLAC 3D based on CAD and ANSYS(辽宁工程技术大学)田树昆曹兰柱TIAN Shu-kun CAO Lan-zhu摘要:FLAC 3D 自带的FISH 语言虽然可以完成复杂模型的建立,但是需要耗费大量时间,同时采用FISH 语言生成复杂模型比较困难,因此提出了一种快速建模方法,即以AutoCAD 与ANSYS 为基础建立FLAC 3D 模型。
采用AutoCAD 及VBA 二次开发程序进行基础数据预处理,在ANSYS 中完成复杂模型建立、材料属性定义、网格划分,并通过ANSYS_TO_FLAC 3D 接口程序,实现FLAC 3D 复杂模型的快速建立。
通过露天矿边坡三维模型的建立与数值模拟分析,验证了该方法的快速性、有效性及可行性。
该方法是综合应用不同程序的优点解决复杂工程地质问题的典型范例。
关键词:数值分析;接口程序;模型转换中图分类号:TP391.9文献标识码:BAbstract:Although it can complete complex model by FISH which in FLAC 3D ,but need to spend a lot of time.At the same time gen -erate complex models using FISH difficult.Therefore,a rapid modeling method proposed.That is to establish model of FLAC 3D based on AutoCAD and ANSYS.The procedures which secondary development using AutoCAD and VBA be used to basic data pre-process -ing.The model established completed in ANSYS and distributes material property,meshing,thereby achieving the purpose of building the complicated model in FLAC 3D by using the interface program of ANSYS_TO_FLAC 3D .Then three -dimensional models of surface mining established and numerical simulation analysis and the case study shows that the fastness,effectiveness and feasibility,can be realized.The procedure gives an example to build complicated engineering geological model with the advantages of different programs.Key words:numerical analysis;interface program;model transformation文章编号:1008-0570(2010)11-3-0259-021引言FLAC 3D 是由美国Itasca 公司开发的有限差分软件,可以进行岩土和其他材料的三维结构受力分析,通过调整三维网格中的多面体单元来拟合实际结构,单元材料可采用线性或非线性本构模型。