一个UDEC用于采矿或地下开挖工程的示例

先介绍一下udec的运算流程和常用命令（udec是一款命令驱动的软件）

1. 流程：
   切割块体—gen–赋块体模型和节理模型—指定边界条件和初始条件—迭代至平衡–模拟开挖—迭代至平衡—检查模型响应

应力保存在单元格上，速度和位移保存在节点上（大的叫网格，小的叫单元格）

有用的例子：稍等，在下面。 4-9为常用命令，分号后代表注释。

4. change cons=0; null, 开挖模型，即开挖
   

change cons=1; elastic model弹性模型，

change cons=3; Mohr-coulomb model一种塑性模型

5. change jcons=1 点接触库仑滑动模型
   

change jcons=2 节理面接触库仑滑动模型

change jcons=5 节理面接触有残余强度的库仑滑动模型

6. change mat…划定块体材料模型作用范围
   

pro mat。。。赋一种材料模型

7. change jmat range。。。划定节理模型作用范围
   

pro jmat。。。赋一种节理模型

8. plot bl/sress/disp或者plot
   

bl sress disp—画出混合图

9. plot mat joint 画出节理图，
    plot block slip lmag open yell识别滑移和张开位置 
   

10. 例子
    new
    round 0.01
    set ovtol=3
    set delc on

    block 0,0 0,207 500,207 500,0

    crack (0,199) (500,199)
    jreg id 1 0,199 0,207 500,207 500,199 delete
    jset 0,0 500,0 0,0 1,0 0,199 range jreg 1
    jset 90,0 1,0 1,0 1,0 0,199 range jreg 1
    jset 90,0 1,0 1,0 1,0 0.5,200 range jreg 1

    crack (0,179) (500,179)
    jreg id 2 0,179 0,199 500,199 500,179 delete
    jset 0,0 500,0 0,0 2,0 0,179 range jreg 2
    jset 90,0 2,0 2,0 3,0 0,179 range jreg 2
    jset 90,0 2,0 2,0 3,0 1.5,181 range jreg 2

    crack (0,175) (500,175)
    jreg id 3 0,175 0,179 500,179 500,175 delete
    jset 0,0 500,0 0,0 2,0 0,175 range jreg 3
    jset 90,0 2,0 2,0 4,0 0,175 range jreg 3
    jset 90,0 2,0 2,0 4,0 2,177 range jreg 3
    ;…
    crack (0,20) (500,20);煤层
    jreg id 14 0,20 0,29 500,29 500,20 delete
    jset 0,0 500,0 0,0 2,0 0,20 ran jreg 14
    jset 90,0 2,0 0,0 2,0 0,20 range jreg 14

    crack (0,0) (500,0);底板
    jreg id 15 0,0 0,20 500,20 500,0 delete
    jset 0,0 500,0 0,0 5,0 0,0 range jreg 15
    jset 90,0 5,0 5,0 10,0 0,0 range jreg 15
    jset 90,0 5,0 5,0 10,0 5,5 range jreg 15

    gen quad 8 range 0,500 0,207

    change cons=1
    change jcons=2
    change mat=1 ran jreg 1
    change mat=2 ran jreg 2
    change mat=3 ran jreg 3
    change mat=4 ran jreg 4
    change mat=5 ran jreg 5
    change mat=6 ran jreg 6
    change mat=7 ran jreg 7
    change mat=8 ran jreg 8
    change mat=9 ran jreg 9
    change mat=10 ran jreg 10
    change mat=11 ran jreg 11
    change mat=12 ran jreg 12
    change mat=13 ran jreg 13
    change mat=14 ran jreg 14
    change mat=15 ran jreg 15

prop mat=1 d=2400 bulk=20e9 shear=12e9 coh=2e6 ten=1.0e6 fric=25
prop mat=2 d=2100 bulk=18.3e9 shear=8.5e9 coh=3.3e6 ten=1.8e6 fric=28
prop mat=3 d=2500 bulk=20e9 shear=12e9 coh=4.0e6 ten=3.0e6 fric=45
prop mat=4 d=2100 bulk=25e9 shear=11.5e9 coh=3.0e6 ten=2.0e6 fric=30
prop mat=5 d=2400 bulk=20e9 shear=12e9 coh=2.9e6 ten=3.0e6 fric=43
prop mat=6 d=2100 bulk=18.3e9 shear=8.5e9 coh=3.3e6 ten=1.8e6 fric=28
prop mat=7 d=2500 bulk=20e9 shear=12e9 coh=4.0e6 ten=4.0e6 fric=45
prop mat=8 d=2500 bulk=24.2e9 shear=12.5e9 coh=5.5e6 ten=5.0e6 fric=45
prop mat=9 d=1800 bulk=13.2e9 shear=5.7e9 coh=2.0e6 ten=1.8e6 fric=25
prop mat=10 d=2500 bulk=20e9 shear=12e9 coh=4.0e6 ten=3.0e6 fric=45
prop mat=11 d=1800 bulk=13.2e9 shear=5.7e9 coh=2.0e6 ten=1.8e6 fric=25
prop mat=12 d=2500 bulk=20e9 shear=12e9 coh=4.0e6 ten=3.0e6 fric=45
prop mat=13 d=1800 bulk=13.2e9 shear=5.7e9 coh=2.0e6 ten=1.0e6 fric=25
prop mat=14 d=1350 bulk=13.3e9 shear=4.4e9 coh=1.3e6 ten=1.0e6 fric=23
prop mat=15 d=2500 bulk=20e9 shear=12e9 coh=2.9e6 ten=5.0e6 fric=43

change jmat=1 ran jreg 1
change jmat=2 ran jreg 2
change jmat=3 ran jreg 3
change jmat=4 ran jreg 4
change jmat=5 ran jreg 5
change jmat=6 ran jreg 6
change jmat=7 ran jreg 7
change jmat=8 ran jreg 8
change jmat=9 ran jreg 9
change jmat=10 ran jreg 10
change jmat=11 ran jreg 11
change jmat=12 ran jreg 12
change jmat=13 ran jreg 13
change jmat=14 ran jreg 14
change jmat=15 ran jreg 15

prop jmat=1 jkn=6.1e9 jks=2.4e9 jco=1.7e6 jfric=20
prop jmat=2 jkn=6.1e9 jks=2.4e9 jco=1.7e6 jfric=20
prop jmat=3 jkn=6.7e9 jks=2.7e9 jco=2.0e6 jfric=32
prop jmat=4 jkn=8.3e9 jks=3.2e9 jco=1.5e6 jfric=21
prop jmat=5 jkn=6.7e9 jks=2.7e9 jco=1.5e6 jfric=30
prop jmat=6 jkn=6.1e9 jks=2.4e9 jco=1.7e6 jfric=20
prop jmat=7 jkn=6.7e9 jks=2.7e9 jco=2.0e6 jfric=32
prop jmat=8 jkn=8.1e9 jks=3.2e9 jco=2.8e6 jfric=32
prop jmat=9 jkn=4.4e9 jks=1.7e9 jco=1.0e6 jfric=18
prop jmat=10 jkn=6.7e9 jks=2.7e9 jco=2.0e6 jfric=32
prop jmat=11 jkn=4.4e9 jks=1.7e9 jco=1.0e6 jfric=18
prop jmat=12 jkn=6.7e9 jks=2.7e9 jco=2.0e6 jfric=32
prop jmat=13 jkn=4.4e9 jks=1.7e9 jco=1.0e6 jfric=18
prop jmat=14 jkn=4.4e9 jks=1.6e9 jco=0.6e6 jfric=16
prop jmat=15 jkn=6.7e9 jks=2.7e9 jco=1.5e6 jfric=30

bound xvel=0 range -0.01,0.01 -0.01,207.01
bound xvel=0 range 499.99,500.01 -0.01,207.01
bound yvel=0 range -0.01,500.01 -0.01,0.01
insitu stress -2.59e6,0,-5.175e6 xgrad 0 0 0 ygrad 1.25e4 0 2.5e4 szz -2.8e6 zgrad 0 1.25e4
set gravity 0,-10

solve

save chushiyingli.sav

dele range 100 110 20 29

supp 107.75 24.5 angle=90 seg=5 width=4.5 mat=3
prop m=3 sup_kn=-1
table 1 0 0 50e-2 5e6 100e-2 6e6 10 6e6
pl bl support red
damp auto

solve ratio 1.0E-5
sav 10.sav
end