jset 90 0 5.0 0 0 0 5 0 50 12.0 ran jreg 33 ;jset 90 0 1.0 0 1.0 0 10 0 55 13.0 ran jreg 33
;第2层,2.0m泥岩,直接底
jset 0 0 250 0 0 0 2.0 0 0 10.55 ran jreg 2
;第1层,10.0m粉砂岩,老底 jset 0 0 250 0 0 0 5.0 0 0 0 ran jreg 1 gen edg 5
;保存建模文件 sav 1.sav cal fb.txt q new res 1.sav
change cons=1 change jcons=2
set jcondf=2 【命令的含义】
这个命令是对新产生的节理设置参数，一般可以不用。具体可以参见command reference 中的set jcondf n 这个命令，一般可以不用。 ;定义岩石材料标号 ;煤
change mat=1 ran jreg 3 change mat=2 ran jreg 6 ;粉砂岩
16
change mat=3 ran jreg 1 change mat=4 ran jreg 4 ;泥岩
change mat=5 ran jreg 2 change mat=6 ran jreg 7
;粉细砂岩互层
change mat=7 ran jreg 8
;中细砂岩互层 change mat=8 ran jreg 5
;遗留煤层
change mat=9 ran jreg 34
;定义节理材料标号 ;煤
change jmat=1 ran jreg 3 change jmat=2 ran jreg 6 ;粉砂岩
change jmat=3 ran jreg 1 change jmat=4 ran jreg 4 ;泥岩
change jmat=5 ran jreg 2 change jmat=6 ran jreg 7
;粉细砂互层
change jmat=7 ran jreg 8
;中细砂岩互层 change jmat=8 ran jreg 5
;遗留煤层
change jmat=9 ran jreg 34
;定义岩石材料力学特性 ;煤
prop mat=1 d=1.4e3 k=4e9 g=3.5e9 prop mat=2 d=1.4e3 k=4e9 g=3.5e9 ;粉砂岩
pro mat=3 d=2.60e3 k=15e9 g=15e9 pro mat=4 d=2.60e3 k=15e9 g=15e9 ;泥岩
pro mat=5 d=2.50e3 k=6e9 g=9e9 pro mat=6 d=2.50e3 k=6e9 g=9e9
;粉细砂岩互层
pro mat=7 d=2.65e3 k=12e9 g=12e9
;中细砂岩互层
pro mat=8 d=2.75e3 k=15e9 g=15e9
;遗留煤层
pro mat=9 d=1.50e3 k=5e9 g=4e9
;定义节理材料力学特性
;煤 prop
jmat=1
jkn=2e9,jks=0.5e9,jcoh=0.2e6,jfric=3,jtens=0.1e6 prop
jmat=2
jkn=3e9,jks=0.5e9,jcoh=0.2e6,jfric=6,jtens=0.1e6 ;粉砂岩 pro
jmat=3
jkn=6e9,jks=1.5e9,jcoh=0.2e6,jfric=15,jtens=0.3e6 pro
jmat=4
jkn=6e9,jks=1.5e9,jcoh=0.2e6,jfric=15,jtens=0.3e6 ;泥岩 pro
jmat=5
jkn=7e9,jks=0.8e9,jcoh=0.2e6,jfric=8,jtens=0.2e6 pro
jmat=6
jkn=7e9,jks=0.8e9,jcoh=0.2e6,jfric=8,jtens=0.2e6
;粉细砂岩互层 pro
jmat=7
jkn=9e9,jks=1.2e9,jcoh=0.3e6,jfric=12,jtens=0.4e6
;中细砂岩互层 pro
jmat=8
jkn=9e9,jks=1.5e9,jcoh=0.3e6,jfric=15,jtens=0.4e6
;遗留煤层 pro
jmat=9
jkn=3e9,jks=0.5e9,jcoh=0,jfric=3,jtens=0
set grav 0 -10
17
;边界条件
bound stress 0,0,-12e6 range 0 250 65.4 65.6 bound xvel=0 range -.1,.1 0,65.5 bound xvel=0 range 249.9,250.1 0,65.5 bound yvel=0 range 0,250 -.1,.1 应力应该施加多大？ ;原岩应力
insitu stress -6e6 0.0 -12e6 szz=-6e6【为什么施加Z向应力】
首先你要知道你所用的二维UDEC，所建立的模型的力学模型是平面应变模型而不是平面应力模型，平面应变模型在z轴方向是没有位移的，因此在z轴方向要施加应力。 ;初次平衡
hist solve_rat type 1 【什么样的命令】 这个命令是记录UDEC的运行模式，每次当你运行UDEC时后，光标处不是有一行不断变化的参数吗？这个命令就是指的这个模式，solve rat 是指在模型运行的时候模型达到平衡的时候的不平衡力的极限比率。具体参见command reference 中的solve ratio命令（P1-219）。UDEC运算默认的都是这个模式类型，你用这个就行了，大家都习惯用这个的。 hist vmax
;第一测线 hist xdis 40.5 49.75 hist xdis 80.5 49.75 hist xdis 100.5 49.75 hist xdis 120.5 49.75 hist xdis 160.5 49.75
18
hist xdis 180.5 49.75
hist ydis 40.5 49.75 hist ydis 80.5 49.75 hist ydis 100.5 49.75 hist ydis 120.5 49.75 hist ydis 160.5 49.75 hist ydis 180.5 49.75
hist xvel 40.5 49.75 hist xvel 80.5 49.75 hist xvel 100.5 49.75 hist xvel 120.5 49.75 hist xvel 160.5 49.75 hist xvel 180.5 49.75
hist yvel 40.5 49.75 hist yvel 80.5 49.75 hist yvel 100.5 49.75 hist yvel 120.5 49.75 hist yvel 160.5 49.75 hist yvel 180.5 49.75
hist sxx 40.5 49.75 hist sxx 80.5 49.75 hist sxx 100.5 49.75 hist sxx 120.5 49.75 hist sxx 160.5 49.75 hist sxx 180.5 49.75
hist syy 40.5 49.75 hist syy 80.5 49.75
hist syy 100.5 49.75 hist syy 120.5 49.75 hist syy 160.5 49.75 hist syy 180.5 49.75
hist nstr 40.5 49.75【是法向应力吗？意义】 在这个点40.5 49.75（x y）附近接触的法向应力。具体可以参见command reference 中的hist nstr x y 命令（P1-81） hist nstr 80.5 49.75 hist nstr 100.5 49.75 hist nstr 120.5 49.75 hist nstr 160.5 49.75 hist nstr 180.5 49.75
hist sstr 40.5 49.75 hist sstr 80.5 49.75 hist sstr 100.5 49.75 hist sstr 120.5 49.75 hist sstr 160.5 49.75 hist sstr 180.5 49.75
hist sdis 40.5 49.75 hist sdis 80.5 49.75 hist sdis 100.5 49.75 hist sdis 120.5 49.75 hist sdis 160.5 49.75 hist sdis 180.5 49.75
hist svel 40.5 49.75 hist svel 80.5 49.75 hist svel 100.5 49.75
19
hist svel 120.5 49.75 hist svel 160.5 49.75 hist svel 180.5 49.75
hist sxy 40.5 49.75 hist sxy 80.5 49.75 hist sxy 100.5 49.75 hist sxy 120.5 49.75 hist sxy 160.5 49.75 hist sxy 180.5 49.75
hist ndis 40.5 49.75 hist ndis 80.5 49.75 hist ndis 100.5 49.75 hist ndis 120.5 49.75 hist ndis 160.5 49.75
hist nvel 40.5 49.75 hist nvel 80.5 49.75 hist nvel 100.5 49.75 hist nvel 120.5 49.75 hist nvel 160.5 49.75 hist nvel 180.5 49.75
;第二测线 hist xdis 40.5 27 hist xdis 80.5 27 hist xdis 100.5 27 hist xdis 120.5 27 hist xdis 160.5 27 hist xdis 180.5 27
hist ydis 40.5 27 hist ydis 80.5 27 hist ydis 100.5 27 hist ydis 120.5 27 hist ydis 160.5 27 hist ydis 180.5 27
hist xvel 40.5 27 hist xvel 80.5 27 hist xvel 100.5 27 hist xvel 120.5 27 hist xvel 160.5 27 hist xvel 180.5 27
hist yvel 40.5 27 hist yvel 80.5 27 hist yvel 100.5 27 hist yvel 120.5 27 hist yvel 160.5 27 hist yvel 180.5 27
hist sxx 40.5 27 hist sxx 80.5 27 hist sxx 100.5 27 hist sxx 120.5 27 hist sxx 160.5 27 hist sxx 180.5 27
hist syy 40.5 27 hist syy 80.5 27 hist syy 100.5 27 hist syy 120.5 27
hist syy 160.5 27 hist syy 180.5 27
hist nstr 40.5 27 hist nstr 80.5 27 hist nstr 100.5 27 hist nstr 120.5 27 hist nstr 160.5 27 hist nstr 180.5 27
hist sstr 40.5 27 hist sstr 80.5 27 hist sstr 100.5 27 hist sstr 120.5 27 hist sstr 160.5 27 hist sstr 180.5 27
hist sdis 40.5 27 hist sdis 80.5 27 hist sdis 100.5 27 hist sdis 120.5 27 hist sdis 160.5 27 hist sdis 180.5 27
hist svel 40.5 27 hist svel 80.5 27 hist svel 100.5 27 hist svel 120.5 27 hist svel 160.5 27 hist svel 180.5 27
hist sxy 40.5 27
20