FLAC3D滑坡模拟

FLAC3D滑坡模拟一、源程序;CreateMaterialZonesgenzonebricksize555&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)genzonebricksize555p0(0,0,5)edge5.0groupMaterial;CreateBinZonesgenzonebricksize155&p0(4,1,0)p1add(3,0,0)p2add(0,3,0)&p3add(2,0,5)p4add(3,6,0)p5add(2,5,5)&p6add(3,0,5)p7add(3,6,5)genzonebricksize155&p0(6,1,5)p1add(1,0,0)p2add(0,5,0)&p3add(0,0,5)p4add(1,6,0)p5add(0,5,5)&p6add(1,0,5)p7add(1,6,5)genzonebricksize515&p0(1,4,0)p1add(3,0,0)p2add(0,3,0)&p3add(0,2,5)p4add(6,3,0)p5add(0,3,5)&p6add(5,2,5)p7add(6,3,5)genzonebricksize515&p0(1,6,5)p1add(5,0,0)p2add(0,1,0)&p3add(0,0,5)p4add(6,1,0)p5add(0,1,5)&p6add(5,0,5)p7add(6,1,5)groupBinrangegroupMaterialnot;Createnamedrangesynonymsrangename=BingroupBinrangename=MaterialgroupMaterial;AssignmodelstogroupsmodelmohrrangeMaterialmodelelasrangeBin;Createinterfaceelementsint1faceranplaneori(4,0,0)nor(-5,0,2)dist0.01z(0,5)y(1,6)int2faceranplaneori(0,4,0)nor(0,-5,2)dist0.01z(0,5)x(1,6)int1faceranx5.96.1y16z510int2faceranx16y5.96.1z510int1maxedge0.55int2maxedge0.55;Movebintowardmaterialinixadd-1.0rangeBininiyadd-1.0rangeBin;Assignpropertiespropshear1e8bulk2e8fric30rangeMaterialpropshear1e8bulk2e8rangeBininiden2000int1propks2e9kn2e9fric15int2propks2e9kn2e9fric15;Assignboundaryconditionsfixxrangex-0.10.1anyx5.96.1any;在x=0和x=6的两个面上加约束fixyrangey-0.10.1anyy5.96.1anyfixzrangez-0.10.1Bin;Monitorhistorieshistunbalhistgpzdisp(6,6,10)histgpzdisp(0,0,10)histgpzdisp(0,0,0);Settingssetlargesetgrav00-10;Cyclingstep4000savebin.sav二、分析结果1、原结构网格图2、z方向位移图3、不平衡力变化图4、histgpzdisp(6,6,10)zdisp(6,6,10)=-3.0928e-0035、histgpzdisp(0,0,10)zdisp(0,0,10)=-1.1589e+0006、histgpzdisp(0,0,0)zdisp(0,0,0)=-2.4947e+000