newrestored:\Box1\meshgroup1;读取转换后模型generatezonebricksize161610&;生成“底”p0=(10,10,-20)&p1=(1610,10,-20)&p2=(10,1610,-20)&p3=(10,10,0)group2attachfacerangez=-0.01,0.01;融合;模型及参数============================modelmohrpropertybulk=2.08e8shear=9.6e7friction=25propertycohesion=1e10tension=1e10;边界条件=======================fixxrangex9.9910.01fixxrangex1609.991610.01fixyrangey9.9910.01fixyrangey1609.991610.01fixzrangez-19.9920.01setgrav0,0,-10inidensity=1860solvesaved:\Box1\z1.sav;开挖===========================modelnullrangex=410,710y=1010,1410z=-20,0modelnullrangex=710,1310y=410,1310z=-20,0savez2.savinitialxdisplacement=0ydisplacement=0zdisplacement=0propertycohesion=1e5tension=1e4;定义剩余煤层modelmohrrangez=-20,0propertybulk=3.846e8shear=7.258e8friction=32propertycohesion=5e5tension=2e5inidensity=2550setlargestep=500;计算步数savez3.sav;显示开挖位移变化=====================plotcreateDispContplotaddcontourzdisp;添加等值线z位移量plotaddaxesplotshow