set large 设与不设对蠕变的影响

nizaina

静力开挖并应力平衡以后开始进行蠕变计算，当在静力开挖时设置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

qq384681566

把参数调大大点试试

hnlgdxnyckwl

"顶板穿透底板很大值"表明需要设置interface