FLAC3D滑坡模拟.docVIP

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 2、z方向位移图 3、不平衡力变化图 4、hist gp zdisp (6,6,10)