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
