静力开挖并应力平衡以后开始进行蠕变计算，当在静力开挖时设置set large 时，蠕变一直进行下去，不受控制，顶板穿透底板很大值。当静力开挖不设置 set large时，蠕变计算没多长时间就有单元格压坏导致无法进行蠕变，求高手指导。现贴出文件
new
plot block group
;建模 中砂岩
gen zone brick size 60 100 2 p0 0 0 0 p1 120 0 0 p2 0 200 0 p3 0 0 8.40 group 1
;粗砂岩
gen zone brick size 60 100 2 p0 0 0 8.40 p1 120 0 8.40 p2 0 200 8.40 p3 0 0 16.85 group 2
; 细沙岩
gen zone brick size 60 100 1 p0 0 0 16.85 p1 120 0 16.85 p2 0 200 16.85 p3 0 0 20 group 3
;8#煤
gen zone brick size 60 100 3 p0 0 0 20 p1 120 0 20 p2 0 200 20 p3 0 0 24.50 group 4
;砂泥岩
gen zone brick size 60 100 1 p0 0 0 24.50 p1 120 0 24.50 p2 0 200 24.50 p3 0 0 28.70 group 5
;粉砂岩
gen zone brick size 60 100 2 p0 0 0 28.70 p1 120 0 28.70 p2 0 200 28.70 p3 0 0 35.78 group 6
;粗砂岩
gen zone brick size 60 100 1 p0 0 0 35.78 p1 120 0 35.78 p2 0 200 35.78 p3 0 0 39.81 group 7
;煤
gen zone brick size 60 100 1 p0 0 0 39.81 p1 120 0 39.81 p2 0 200 39.81 p3 0 0 41.13 group 8
;中砂岩
gen zone brick size 60 100 2 p0 0 0 41.13 p1 120 0 41.13 p2 0 200 41.13 p3 0 0 53.18 group 9
;粗砂岩
gen zone brick size 60 100 3 p0 0 0 53.18 p1 120 0 53.18 p2 0 200 53.18 p3 0 0 64.08 group 10
;细砂岩
gen zone brick size 60 100 2 p0 0 0 64.08 p1 120 0 64.08 p2 0 200 64.08 p3 0 0 70.18 group 11
;砂泥岩
gen zone brick size 60 100 2 p0 0 0 70.18 p1 120 0 70.18 p2 0 200 70.18 p3 0 0 76.85 group 12
;5#煤
gen zone brick size 60 100 3 p0 0 0 76.85 p1 120 0 76.85 p2 0 200 76.85 p3 0 0 85.85 group 13
;砂泥岩
gen zone brick size 60 100 1 p0 0 0 85.85 p1 120 0 85.85 p2 0 200 85.85 p3 0 0 88.95 group 14
;中砂岩
gen zone brick size 60 100 2 p0 0 0 88.95 p1 120 0 88.95 p2 0 200 88.95 p3 0 0 100 group 15
;粗砂岩
gen zone brick size 60 100 2 p0 0 0 100 p1 120 0 100 p2 0 200 100 p3 0 0 110 group 16
;粉砂岩
gen zone brick size 60 100 2 p0 0 0 110 p1 120 0 110 p2 0 200 110 p3 0 0 120 group 17
;* material 1 中砂岩
model moh range group 1
property den 2670 s 7.3e9 co 3.6e6 bul 0.6e10 fric 20 ten 3.4e6 range group 1
;* material 2 粗砂岩
model moh range group 2
property den 2670 s 6e9 co 3.8e6 bul 3.7e9 fric 20 ten 24.2e6 range group 2
;* material 3 细砂岩
model moh range group 3
property den 2400 s 1.1e9 co 4e6 bul 1.7e9 fric 30 ten 2.2e6 range group 3
;* material 4 8#煤
model moh range group 4
property den 1460 s 0.3e9 co 3.2e6 bul 0.6e9 fric 20 ten 0.5e6 range group 4
;* material 5 砂泥岩
model moh range group 5
property den 2630 s 9e9 co 2e6 bul 1.0e9 fric 14 ten 3.8e6 range group 5
;* material 6 粉砂岩
model moh range group 6
property den 2670 s 9e9 co 4.8e6 bul 1.0e10 fric 17 ten 3.1e6 range group 6
;* material 7 粗砂岩
model moh range group 7
property den 2670 s 6e9 co 4.8e6 bul 4.7e9 fric 20 ten 2.2e6 range group 7
;* material 8 煤
model moh range group 8
property den 1460 s 0.3e9 co 3.2e6 bul 0.6e9 fric 20 ten 0.5e6 range group 8
;* material 9 中砂岩
model moh range group 9
property den 2670 s 7.3e9 co 4.6e6 bul 0.6e10 fric 20 ten 4.4e6 range group 9
;* material 10 粗砂岩
model moh range group 10
property den 2670 s 6e9 co 4.8e6 bul 4.7e9 fric 20 ten 2.2e6 range group 10
;* material 11 细砂岩
model moh range group 11
property den 2400 s 1.1e9 co 4e6 bul 1.7e9 fric 30 ten 2.2e6 range group 11
;* materia l2 砂泥岩
model moh range group 12
property den 2630 s 9e9 co 2e6 bul 1e9 fric 14 ten 2.8e6 range group 12
;* material 13 5#煤
model moh range group 13
property den 1460 s 0.3e9 co 3.2e6 bul 0.6e9 fric 20 ten 0.6e6 range group 13
;* materia l4 砂泥岩
model moh range group 14
property den 2630 s 9e9 co 2e6 bul 1e9 fric 14 ten 3.8e6 range group 14
;* material 15 中砂岩
model moh range group 15
property den 2670 s 7.3e9 co 4.6e6 bul 0.8e10 fric 20 ten 4.4e6 range group 15
;* material 16 粗砂岩
model moh range group 16
property den 2670 s 6e9 co 4.8e6 bul 4.7e9 fric 20 ten 2.2e6 range group 16
;* material 17 粉砂岩
model moh range group 17
property den 2670 s 9e9 co 4.8e6 bul 1e10 fric 17 ten 3.1e6 range group 17
save 留煤柱初始1新.sav
;边界条件
fix x range x -0.1 0.1
fix x range x 119.9 120.1
fix y range y -0.1 0.1
fix y range y 199.9 200.1
fix z range z -0.1 0.1
apply szz=-7.5e6 range z=120
initial szz=-7.5e6 grad 0 0 2.5e4
initial sxx=-3.75e6 grad 0 0 1.25e4
initial syy=-3.75e6 grad 0 0 1.25e4
set grav 0 0 -9.81
history unbal
pl add con szz
solve
save 留煤柱初始2.sav
ini state=0
ini xdisp=0 ydisp=0 zdisp=0
m null range x 20 100 y 28 38 z 20 24.50
m null range x 58 62 y 38 44 z 20 24.50
m null range x 20 100 y 44 54 z 20 24.50
m null range x 58 62 y 54 60 z 20 24.50
m null range x 20 100 y 60 70 z 20 24.50
m null range x 58 62 y 70 76 z 20 24.50
m null range x 20 100 y 76 86 z 20 24.50
m null range x 58 62 y 86 92 z 20 24.50
m null range x 20 100 y 92 102 z 20 24.50
m null range x 58 62 y 102 108 z 20 24.50
m null range x 20 100 y 108 118 z 20 24.50
m null range x 58 62 y 118 124 z 20 24.50
m null range x 20 100 y 124 134 z 20 24.50
m null range x 58 62 y 134 140 z 20 24.50
m null range x 20 100 y 140 150 z 20 24.50
m null range x 58 62 y 150 156 z 20 24.50
m null range x 20 100 y 156 166 z 20 24.50
m null range x 58 62 y 166 172 z 20 24.50
m null range x 20 100 y 172 182 z 20 24.50
m null range x 58 62 y 182 188 z 20 24.50
m null range x 20 100 y 188 198 z 20 24.50
set large
solve
save kaiwa6m.sav
group 333 range x 20 100 y 28 38 z 20 24.50
group 333 range x 58 62 y 38 44 z 20 24.50
group 333 range x 20 100 y 44 54 z 20 24.50
group 333 range x 58 62 y 54 60 z 20 24.50
group 333 range x 20 100 y 60 70 z 20 24.50
group 333 range x 58 62 y 70 76 z 20 24.50
group 333 range x 20 100 y 76 86 z 20 24.50
group 333 range x 58 62 y 86 92 z 20 24.50
group 333 range x 20 100 y 92 102 z 20 24.50
group 333 range x 58 62 y 102 108 z 20 24.50
group 333 range x 20 100 y 108 118 z 20 24.50
group 333 range x 58 62 y 118 124 z 20 24.50
group 333 range x 20 100 y 124 134 z 20 24.50
group 333 range x 58 62 y 134 140 z 20 24.50
group 333 range x 20 100 y 140 150 z 20 24.50
group 333 range x 58 62 y 150 156 z 20 24.50
group 333 range x 20 100 y 156 166 z 20 24.50
group 333 range x 58 62 y 166 172 z 20 24.50
group 333 range x 20 100 y 172 182 z 20 24.50
group 333 range x 58 62 y 182 188 z 20 24.50
group 333 range x 20 100 y 188 198 z 20 24.50
config creep
model cvisc range group 4
prop den 1460 s 0.3e9 co 3.2e6 bulk 0.6e9 fric 20 ten 0.5e6 dilation 10 mvis 1.8e12 kvis 1.8e12 range group 4
hist crtime
hist gp zdisp 50 97 24.50
hist gp xdisp 50 92 23
hist gp zdisp 50 97 20
hist gp zdisp 50 89 24.50
pl hist 3 v 2
set creep mindt=1 maxdt=6e3
set creep dt auto on
solve age 2592000
save 1yue.sav
solve age 7776000
save 3yue.sav
solve age 18144000
save 7yue.sav
solve age 63072000
save 2nian.sav
