FLAC3D滑坡模拟.doc

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