案例请教

ZHANZIQIANG

想模拟一个煤层开挖导致裂隙演化的模型，该模型分六层，从上至下依次是：上覆岩层、基本顶、直接顶、煤、直接底、基本底，命令如下。想问问各位前辈命令是不是有错误的地方？还有就是怎么出裂隙的图？
new
round=0.01
set ovtol=3
block 0,0 0,45.2 180,45.2 180,0

crack 0,3 180,3
crack 0,4.5 180,4.5
crack 0,10.2 180,10.2
crack 0,11.2 180,11.2
crack 0,21.2 180,21,2
crack 0,45.2 180,45.2

gen quad 1.6 1.2 range 0,180 0,11.2
gen quad 13 5.5 range 0,180 11.2 45.2

change cons=3 range 0,180 0,45.2
change jcons=2 range 0,180 0,45.2

change materia=1 range region 0,0 0,3 180,3 180,0
change materia=2 range region 0,3 0,4.5 180,4.5 180,3
change materia=3 range region 0,4.5 0,10.2 180,10.2 180,4.5
change materia=4 range region 0,10.2 0,11.2 180,11.2 180,10.2
change materia=5 range region 0,11.2 0,21.2 180,21.2 180,11.2
change materia=6 range region 0,21.2 0,45.2 180,45.2 180,21.2

change jmat=1 range region 0,0 0,3 180,3 180,0
change jmat=2 range region 0,3 0,4.5 180,4.5 180,3
change jmat=3 range region 0,4.5 0,10.2 180,10.2 180,4.5
change jmat=4 range region 0,10.2 0,11.2 180,11.2 180,10.2
change jmat=5 range region 0,11.2 0,21.2 180,21.2 180,11.2
change jmat=6 range region 0,21.2 0,45.2 180,45.2 180,21.2

property mat=1 dens=2500,k=2e10,g=1.1e10,fric=40,coh=8.5e6,ten=5.6e6
property mat=2 dens=2500,k=17e9,g=8.9e9,fric=34,coh=6e6,ten=2.4e6
property mat=3 dens=1390,k=13.05e9,g=6.3e9,fric=28,coh=4e6,ten=2.1e6
property mat=4 dens=2500,k=2e10,g=1.1e10,fric=33,coh=6.2e6,ten=2.3e6
property mat=5 dens=2700,k=2.4e10,g=1.3e10,fric=40,coh=9.5e6,ten=5.4e6
property mat=6 dens=2500,k=2e10,g=1.1e10,fric=35,coh=8.4e6,ten=4e6

property jmat=1 jkn=14e9,jks=8.7e9,jcoh=6.5e6,jfric=35,jten=0
property jmat=2 jkn=9.5e9,jks=6.4e9,jcoh=4.8e6,jfric=28,jten=0
property jmat=3 jkn=6e9,jks=4.7e9,jcoh=2.8e6,jfric=20,jten=0
property jmat=4 jkn=8.9e9,jks=5.9e9,jcoh=4.9e6,jfric=27,jten=0
property jmat=5 jkn=16e9,jks=9.5e9,jcoh=7.8e6,jfric=36,jten=0
property jmat=6 jkn=14e9,jks=8.7e9,jcoh=6.5e6,jfric=29,jten=0

set gravity 0,-10
boundary stress 0,0,-1.375e7 range 0 180 45.1,45.3
insitu stress -1.5144e7,0,-1.262e7 szz=-1.5144e7 ygrad 3e4 0 2.5e4 zgrad 0 3e4
solve

delete range region 40,4.5 40,10.2 50,10.2 50,4.5step 50000
save kw-10.sav

res kw-10.sav
delete range region 40,4.5 40,10.2 60,10.2 60,4.5
step 50000
save kw-20.sav

res kw-20.sav
delete range region 40,4.5 40,10.2 70,10.2 70,4.5
step 50000
save kw-30.sav

res kw-30.sav
delete range region 40,4.5 40,10.2 80,10.2 80,4.5
step 50000
save kw-40.sav

res kw-40.sav
delete range region 40,4.5 40,10.2 90,10.2 90,4.5
step 50000
save kw-50.sav