求教！UDEC模拟煤层开挖 没法达到初始平衡

lizhi17

求教！！我建了一个模型 但是没法达到初始平衡，一直在循环，这个跟什么有关系的？除了边界条件，还跟什么有关？望高手不吝赐教~附件里有所有的命令。以下也是
round 0.05
blo 0,0 0,200 300,200 300,0
;底板
;01砂岩（69m）
cr 0,0 300,0
cr 0,69 300,69
;02煤层
cr 0,71 300,71
;03泥岩
cr 0,75 300,75
;04煤层
cr 0,76 300,76
;05砂岩（直接底）
cr 0,80 300,80
;06开挖煤层
cr 0,82 300,82
;顶板
;07泥岩（直接顶）
cr 0,86 300,86
;08砂质泥岩
cr 0,89 300,89
;09煤层
cr 0,90 300,90
;10砂岩
cr 0,93 300,93
;11泥砂岩互层
cr 0,104 300,104
;12砂岩
cr 0,200 300,200
;开挖分块
cr 75,80 75,82
cr 225,80 225,82
;节理划分
;01砂岩(69m)
jreg id 1 0,0 0,69 300,69 300,0
jset 0,0 300,0 14,0 14,0 0,14 range jreg 1
jset 90,0 14,0 14,0 20,0 20,0 range jreg 1
jset 90,0 14,0 14,0 20,0 10,14 range jreg 1
;02煤层(2m)
jreg id 2 0,69 0,71 300,71 300,69 range jreg 2
jset 90,0 2,0 2,0 15,0 7.5,69 range jreg 2
;03泥岩(4)
jreg id 3 0,71 0,75 300,75 300,71
jset 90,0 4,0 4,0 15,0 15,71 range jreg 3
;04煤层(1m)
jreg id 4 0,75 0,76 300,76 300,75
jset 90,0 1,0 1,0 10,0 10,75 range jreg 4
;05砂岩(4m)
jreg id 5 0,76 0,80 300,80 300,76
jset 0,0 300,0 2,0 2,0 0,78 range jreg 5
jset 90,0 2,0 2,0 10,0 5,76 range jreg 5
jset 90,0 2,0 2,0 10,0 10,78 range jreg 5
;06开挖煤层(2m)
jreg id 6 0,80 0,82 300,82 300,80
jset 0,0 300,0 300,0 0.5,0 0,80.5 range jreg 6
jset 90,0 0.5,0 0.5,0 1,0 0.5,80 range jreg 6
jset 90,0 0.5,0 0.5,0 1,0 1,80.5 range jreg 6
;07泥岩(4m)
jreg id 7 0,82 0,86 300,86 300,82
jset 0,0 300,0 300,0 1,0 0,83 range jreg 7
jset 90,0 1,0 1,0 1,0 2.5,82 range jreg 7
jset 90,0 1,0 1,0 1,0 5,83 range jreg 7
;08砂质泥岩(3m)
jreg id 8 0,86 0,89 300,89 300,86
jset 0,0 300,0 0,0 1.5,0 0,87.5 range jreg 8
jset 90,0 1.5,0 1.5,0 10,0 10,86 range jreg 8
jset 90,0 1.5,0 1.5,0 10,0 5,87.5 range jreg 8
;09煤层(1m)
jreg id 9 0,89 0,90 300,90 300,89
jset 90,0 1,0 1,0 10,0 10,89 range jreg 9
;10砂岩(3m)
jreg id 10 0,90 0,93 300,93 300,90
jset 90,0 3,0 3,0 15,0 7.5,90 range jreg 10
;11砂泥岩互层(11m)
jreg id 11 0,93 0,104 300,104 300,93
jset 0,0 300,0 0,0 4,0 0,97 range jreg 11
jset 90,0 4,0 4,0 15,0 15,93 range jreg 11
jset 90,0 4,0 4,0 15,0 7.5,97 range jreg 11
;12砂岩(96m)
jreg id 12 0,104 0,200 300,200 300,104
jset 0,0 300,0 0,0 14,0 0,118 range jreg 12
jset 90,0 14,0 14,0 20,0 20,104 range jreg 12
jset 90,0 14,0 14,0 20,0 10,118 range jreg 12
change con=3
;设置材料号
change mat 1 range jreg 1
change mat 2 range jreg 2
change mat 3 range jreg 3
change mat 2 range jreg 4
change mat 1 range jreg 5
change mat 2 range jreg 6
change mat 3 range jreg 7
change mat 4 range jreg 8
change mat 2 range jreg 9
change mat 1 range jreg 10
change mat 5 range jreg 11
change mat 1 range jreg 12
;设置材料参数
prop mat 1 dens=2780 bulk=11.5e9 shear=9.5e9 fric=36 coh=1.3e7 tens=1e6
prop mat 2 dens=1400 bulk=2.78e9 shear=2.08e9 fric=35 coh=0.75e6 tens=3e5
prop mat 3 dens=2560 bulk=4.29e9 shear=3.92e9 fric=32 coh=1e6 tens=0.8e6
prop mat 4 dens=2560 bulk=6.29e9 shear=6.92e9 fric=35 coh=1.2e7 tens=1e6
prop mat 5 dens=2660 bulk=6.5e9 shear=6.75e9 fric=35 coh=0.5e6 tens=0.5e6
;设置节理材料号
change jmat 1 range jreg 1
change jmat 2 range jreg 2
change jmat 3 range jreg 3
change jmat 4 range jreg 4
change jmat 5 range jreg 5
change jmat 6 range jreg 6
change jmat 7 range jreg 7
change jmat 8 range jreg 8
change jmat 9 range jreg 9
change jmat 10 range jreg 10
change jmat 11 range jreg 11
change jmat 12 range jreg 12
;设置节理参数
prop jmat 1 jkn=17.24e9 jks=7.33e9 jfric=36
prop jmat 2 jkn=27.73e9 jks=11.55e9 jfric=35
prop jmat 3 jkn=23.76e9 jks=10.33e9 jfric=32
prop jmat 4 jkn=55.46e9 jks=23.09e9 jfric=35
prop jmat 5 jkn=120.68e9 jks=51.31e9 jfric=36
prop jmat 6 jkn=110.93e9 jks=46.18e9 jfric=35      
prop jmat 7 jkn=95.04e9 jks=41.32e9 jfric=32
prop jmat 8 jkn=103.29e9 jks=47.06e9 jfric=35
prop jmat 9 jkn=55.46e9 jks=23.09e9 jfric=35
prop jmat 10 jkn=80.45e9 jks=34.2e9 jfric=36  
prop jmat 11 jkn=38.69e9 jks=14.71e9 jfric=35
prop jmat 12 jkn=17.24e9 jks=7.33e9 jfric=36
;施加初始应力条件  r=25000  按开挖煤层上表面埋深400米计算
bound stress 0,0,-7.05e6 range -0.1,300.1 199.9,200.1
insitu stress -9.64e6,0,-7.05e6 szz=-9.64e6 & ygard -2e4,0,0 zgrad 0,-2e4
;施加初始位移条件
bound xvel=0.0 range -0.1,0.1 -0.1,200.1
bound xvel=0.0 range 299.9,300.1 -0.1,200.1
bound yvel=0.0 range -0.1,300.1 -0.1,0.1
set grav 0.0 -9.81
hist unbal
solve rat 1e-5

mary2882671

第一：你没有生成单元（zone），所以你的块体是刚体，所加边界条件失效。
第二：你的insitu stress命令用的不太对，如果你想迅速收敛的话，最好一层一层岩层的使用，这样初始不平衡力才会比较小，收敛速度才会快一些。

cumtjh

gen edge n

tianjl619

楼上是个牛人。学习下。