FLAC3D滑坡模拟.doc

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