FLAC3D滑坡模拟
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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