UDEC指令

;new

title
moni ;定义名称moni
round=0.025
set ovtol=0.4 ;常规设置

bl 0,0 0,100 300,100 300,0 ;定义范围，四点坐标，顺时针方向

cr 0,15 300,15 ;画直线，两点坐标


jreg id 1 0,0 0,15 300,15 300,0 delete ;定义节理，命名1，定义范围，四点坐标，delete为常规语言
jset 90,0 5,0 5,0 5,0 0,0 rang jreg 1 ;定义节理角度90°；节理的长度；节理的距离，即纵向间隔；横向距离，即横向隔5一个；起始点坐标为0,0
jset 90,0 5,0 5,0 5,0 2.5,5 rang jreg 1 ;定义节理的另外一项

gen quad 10 range 0 300 0 100 ;基本上都有
zone model mo range 0 300 0 100 ;定义摩尔库伦模型的范围，X的范围，Y的范围
change jcons=5 range 0 300 0 100 ;定义节理的某个属性


change mat=1 range reg 0,0 0,15 300,15 300,0 ;定义物质1的范围，即赋予这个范围为物质1,只赋名，没有定义属性
change mat=2 range reg 0,15 0,20 300,20 300,15
change mat=3 range reg 0,20 0,23 300,23 300,20
change mat=4 range reg 0,23 0,28 300,28 300,23

change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围


set jmatdf 2 ;定义节理属性，使符合摩尔库伦准则普适
save ch.sav ;保存文件
cal prop-25.txt ;调出文件


res ch.sav ;调出前面保存的文件接着计算

;参数设置

prop mat=1 dens=2720 ;物质的密度


zone k=12.12e9 g=10.26e9 fric=33 coh=5.77e6 ten=3.6e6 range mat=1
;定义物质1的体积模量、剪切模量、内摩擦角、粘聚力、抗拉强度

prop jmat=1 jkn=11e8 jks=12e7 jfric=12 jcoh=0 jten=0
;定义节理的属性法向刚度、切向刚度、内摩擦角、粘聚力、抗拉强度

;边界条件
boun stress 0 0 0 range 0,300 99.9 100.1 ;定义边界条件，三向应力为0，上边界，范围，纵向波动范围
boun stress 0 0 -8.6095397e6 range ;定义边界应力条件，竖直方向的应力为，负号表示方向向下
insitu str -7.5e6 0 -15e6 szz -7.5e6 ygrad 1.1018e4 0 2.203607e4 ;定义应力属性，sxx方向即水平应力，sxy方向即剪切应力，syy方向即竖直应力，szz方向的应力，ygrad应力梯度
set gravity 0,-10 ;重力加速度
bound xvel=0 range -0.1 0.1 0 100 ;X位移边界,位移波动范围，范围
bound xvel=0 range 299.9 300.1 0 100 ;X位移边界，同上
bound yvel=0 range -0.1 300.1 -0.1,0.1 ;Y位移边界

hist solve_ratio type 1 ;设置不平衡力的精度，普适
solve ;计算
save 25-1.sav ;保

存
cal 25-2.txt ;调出文件

res 25-1.sav ;调出前面保存的文件
del 80,110 15,20 ;永久删除这个范围，即开挖，不可再充填，X的范围，Y的范围
model null range ;变性删除这个范围，即之后可以充填物质，X的范围，Y的范围
model m range ;充填这个范围，
solve ;计算
save 25-301.sav ;再保存，周而复始


set pline 25,20 300,20 10 ;定义测线，起始点坐标，测线分段

reset hist
set log on ;定义hist,保存在一个位置，在程序中输入hist1即可调出

hist xdisp 10,20 ;测该点的X方向位移
hist sxx 10,20 ;测该点X方向的应力
hist syy 10,20 ;测该点Y方向的应力

reset disp
reset hist ;位移清零，历史记录清零

his unbal
solve ratio 1e-5 ;记录不平衡力，计算直到精度达到1e-5

ch cons=3 range 77.5,80 0,3
ch mat=7 range 77.5,80 0,3 ;应变软化模型，这个是做充填用的

cable (73.4,2.91) (71.15,4.99) 10 12 380e-6 11 ;锚杆，起始点坐标，后面四个数值可固定

prop mat 12 cb_dens 7800 cb_ycomp 430e6 cb_yield 160e3 cb_ymod 1.3e11

prop mat 11 cb_kbond 6.3e9 cb_sbond 6e5 ;定义锚杆的属性，密度、抗拉强度、屈服强度

pl smax fill, not sman