udce计算应力图，高手看看不光滑、连续？是参数设置还是

lydia

是参数设置还是网格网格划分的问题

lydia

round=0.1
set ovtol=3
bl 0,0 0,100 300,100 300,0
;基本底10m(2.5*2.5)
crack 0,10 300,10

jregion id 1 0,0 0,10 300,10 300,0 delete
jset 0,0 300,0 0,0 1,0 300,0 range jregion 1
jset 90,0 1,0 1,0 2,0 300,0 range jregion 1
jset 90,0 1,0 1,0 2,0 1,1 range jregion 1


;3米的煤层
crack 0,13 300,13


jregion id 2 0,10 0,13 300,13 300,10 delete
jset 0,0 5,0 0,0 1,0 300,10 range jregion 2
jset 90,0 1,0 1,0 1,0 300,10 range jregion 2
jset 90,0 1,0 1,0 1,0 0,11 range jregion 2




;2米煤层
crack 0,15 300,15

jregion id 3 0,13 0,15 300,15 300,13 delete
jset 0,0 5,0 0,0 1,0 0,13 range jregion 3
jset 90,0 1,0 1,0 1,0 0,13 range jregion 3
jset 90,0 1,0 1,0 1,0 0,14 range jregion 3




;直接顶 砂质泥岩 5m
crack 0,20 300,20
jregion id 26 0,15 0,20 300,20 300,15 delete
jset 0,0 300,0 0,0 1,0 0,15 range jregion 26
jset 90,0 1,0 1,0 2,0 300,15 range jregion 26
jset 90,0 1,0 1,0 2,0 2.5,16 range jregion 26

;老顶 细砂岩 6m
crack 0,26 300,26
jregion id 27 0,20 0,26 300,26 300,20 delete
jset 0,0 300,0 0,0 2,0 0,20 range jregion 27
jset 90,0 2,0 2,0 9,0 300,22 range jregion 27
jset 90,0 2,0 2,0 9,0 5,20 range jregion 27

;泥岩 2m
crack 0,28 300,28
jregion id 28 0,26 0,28 300,28 300,26 delete
jset 0,0 300,0 0,0 1,0 0,26 range jregion 28
jset 90,0 1,0 1,0 2,0 300,23 range jregion 28
jset 90,0 1,0 1,0 2,0 3,26 range jregion 28

;中砂岩10m
crack 0,38 300,38
jregion id 29 0,28 0,38 300,38 300,28 delete
jset 0,0 300,0 0,0 2.5,0 0,28 range jregion 29
jset 90,0 2.5,0 2.5,0 4,0 300,30.5 range jregion 29
jset 90,0 2.5,0 2.5,0 4,0 3,28 range jregion 29

;细砂岩 8m

crack 0,46 300,46
jregion id 30 0,38 0,46 300,46 300,38 delete
jset 0,0 300,0 0,0 2,0 0,38 range jregion 30
jset 90,0 2,0 2,0 9,0 300,40 range jregion 30
jset 90,0 2,0 2,0 9,0 5,42 range jregion 30

; 砂质泥岩 10m
crack 0,56 300,56
jregion id 31 0,46 0,56 300,56 300,46 delete
jset 0,0 300,0 0,0 1,0 0,46 range jregion 31
jset 90,0 1,0 1,0 2,0 300,47 range jregion 31
jset 90,0 1,0 1,0 2,0 2.5,46 range jregion 31


; 粉砂岩 25m

crack 0,81 300,81
jregion id 32 0,56 0,81 300,81 300,56 delete
jset 0,0 300,0 0,0 5,0 0,56 range jregion 32
jset 90,0 5,0 5,0 20,0 300,61 range jregion 32
jset 90,0 5,0 5,0 20,0 5,56 range jregion 32




; 泥岩 19m

crack 0,100 300,100

jregion id 33 0,81 0,100 300,100 300,81 delete
jset 0,0 300,0 0,0 1,0 0,81 range jregion 33
jset 90,0 1,0 1,0 2,0 300,82 range jregion 33
jset 90,0 1,0 1,0 2,0 3,81 range jregion 33


;中砂岩 20m
;crack 0,120 300,120
;jregion id 34 0,100 0,120 300,120 300,100 delete
;jset 0,0 300,0 0,0 2.5,0 0,100 range jregion 34
;jset 90,0 2.5,0 2.5,0 4,0 300,100 range jregion 34
;jset 90,0 2.5,0 2.5,0 4,0 3,102.5 range jregion 34



;黄土30m

;jregion id 35 0,120 0,150 300,150 300,120 delete
;jset 0,0 300,0 0,0 2.5,0 0,100 range jregion 35
;jset 90,0 2.5,0 2.5,0 4,0 300,100 range jregion 35
;jset 90,0 2.5,0 2.5,0 4,0 3,102.5 range jregion 35





save inir（hf）.save
res inir（hf）.save
gen quad 1.5  1.5 range  0,300 0,10
gen quad 1  1 range  0,300 10,21
gen quad 10  5 range 0,300 21,100

zone model mo range 0,300 0,100 ;( change cons=3 range 0,300 0,100)
change jcons=2 range  0,300 0,100

;设定不同岩层、煤层
change mat=1 range reg 0,0 0,10 300,10 300,0
change mat=2 range reg 0,10 0,15 300,15 300,10
change mat=3 range reg 0,15 0,20 300,20 300,15
change mat=4 range reg 0,20 0,26 300,26 300 20
change mat=5 range reg 0,26 0,28 300,28 300,26
change mat=6 range reg 0,28 0,38 300,38 300,28
change mat=7 range reg 0,38 0,46 300,46 300,38
change mat=8 range reg 0,46 0,56 300,56 300,46

change mat=9 range reg 0,56 0,81 300,81 300,56
change mat=10 range reg 0,81 0,100 300,100 300,81




change jmat=1 range reg 0,0 0,10 300,10 300,0
change jmat=2 range reg 0,10 0,15 300,15 300,10
change jmat=3 range reg 0,15 0,20 300,20 300,15
change jmat=4 range reg 0,20 0,26 300,26 300 20
change jmat=5 range reg 0,26 0,28 300,28 300,26
change jmat=6 range reg 0,28 0,38 300,38 300,28
change jmat=7 range reg 0,38 0,46 300,46 300,38
change jmat=8 range reg 0,46 0,56 300,56 300,46

change jmat=9 range reg 0,56 0,81 300,81 300,56
change jmat=10 range reg 0,81 0,100 300,100 300,81



prop mat=1 dens=2430;砂质泥岩
prop mat=2 dens=1520;2#煤
prop mat=3 dens=2430;砂质泥岩
prop mat=4 dens=2380;细砂岩
prop mat=5 dens=2360;中砂岩
prop mat=6 dens=2300;泥岩
prop mat=7 dens=2380;细砂岩
prop mat=8 dens=2430;砂质泥岩

prop mat=9 dens=2380;粉砂岩
prop mat=10 dens=2300;泥岩



zone k=8.7e9,    g=10e9,    fric=40,  coh=2e7,   ten=2.4e6      range mat=1;砂质泥岩
zone k=2.7e9,    g=7e9,    fric=35,  coh=1e7,   ten=1.2e6      range mat=2;煤层
zone k=20e9,    g=11e9,    fric=40,  coh=2e7,   ten=2.4e6     range mat=3;砂质泥岩
zone k=40e9,    g=20e9,    fric=40,  coh=3.9e7, ten=5.4e6     range mat=4;细砂岩
zone k=8.7e9,   g=10e9,    fric=40,  coh=1.2e7, ten=2.4e6      range mat=5;泥岩
zone k=26.8e9,  g=7e9,     fric=45,  coh=3.5e7, ten=4e6      range mat=6;中砂岩
zone k=40e9,     g=20e9,    fric=40,  coh=3.9e7, ten=5.4e6     range mat=7;细砂岩
zone k=8.7e9,    g=10e9,    fric=40,  coh=2e7,   ten=2.4e6     range mat=8;砂质泥岩
zone k=15.6e9,  g=10.8e9,    fric=40,  coh=3.9e7, ten=5.4e6     range mat=9;粉砂岩
zone k=20e9,    g=11e9,    fric=40,  coh=2e7,   ten=2.4e6     range mat=10;泥岩





prop jmat=1 jkn=14e9,  jks=8.7e9, jcoh=6.5e6 , jfric=35,jten=0
prop jmat=2  jkn=7.5e9,  jks=6.7e9, jcoh=2.8e6 , jfric=20,jten=0
change jcons=5
prop jmat=3 jkn=12e9,  jks=7.5e9, jcoh=6.5e6 , jfric=35,jten=0
set jcondf=5 jmatdf=3
prop jmat=4 jkn=14e9,  jks=8.7e9, jcoh=4.9e6, jfric=38,jten=0
set jcondf=5 jmatdf=4
prop jmat=5 jkn=14e9,  jks=8.7e9, jcoh=6.5e6 , jfric=36,jten=0
prop jmat=6  jkn=14e9,  jks=8.7e9, jcoh=6.5e6 , jfric=29,jten=0
prop jmat=7 jkn=8.9e9,  jks=5.9e9, jcoh=4.9e6, jfric=27,jten=0
prop jmat=8 jkn=14e9,  jks=8.7e9, jcoh=6.5e6 , jfric=35,jten=0

prop jmat=9 jkn=8.9e9,  jks=5.9e9, jcoh=4.9e6, jfric=27,jten=0
prop jmat=10 jkn=14e9,  jks=8.7e9, jcoh=6.5e6 , jfric=35,jten=0



;工作面埋深150米
set gravity 0,-10

bound stress 0,0,-0.1e6 ygrad -1.34e4 0 -1.34e4
insitu str  -0.1e6 0 -0.1e6  szz=-0.1e6 ygrad -1.34e4 0 -1.34e4 zgra 0 -1.34e4


bound xvel=0 range -0.1 0.1 0 100
bound xvel=0 range 299.9 300.1 0 100
bound yvel=0 range 0 300 -0.1 0.1

solve

save hfyuanyan.save

plot pen block
plot pen block ydisplace fill


res hfyuanyan.save
reset disp
reset vel

;开掘巷道5
set pline 0,15 300,15 600
set pline 0,10 300,10 600

;del range reg 30,13 30,15 35,15 35,13
;del range reg 35,13 35,15 40,15 40,13
;del range reg 40,13 40,15 45,15 45,13
;del range reg 45,13 45,15 50,15 50,13
del range reg 50,13 50,15 55,15 55,13
solve
;del range reg 55,13 55,15 60,15 60,13
del range reg 60,13 60,15 65,15 65,13
solve
;del range reg 65,13 65,15 50,15 70,13
del range reg 70,13 70,15 75,15 75,13
solve
;del range reg 75,13 75,15 80,15 80,13
del range reg 80,13 80,15 85,15 85,13
solve

solve

lydia

找到问题了，是初始应力的事情，看来还的好好研究研究命令了，不能急于求成。

被窝里的天

lydia 发表于 2017-5-11 19:08
找到问题了，是初始应力的事情，看来还的好好研究研究命令了，不能急于求成。 ...

亲，想请问你gen quad 1.5 1.5 range 0,300 0,10这个命令里面两个1.5分别代表什么，为什么这么设啊