Flac3D模型最大不平衡力一直不收敛.急急急

cumt_baicai

我模拟的是采矿的问题，由于flac3d的局限性，先由Autocad绘图，然后导入ansys建立模型，由于模型存在采空区和弱弱夹层，需要设置interface接触面，故根据“导来导去法”将模型分成六部分分别建立，最后导入flac3d进行模拟，在模型初始平衡后，进行采空区开挖，但是开挖后，采空区顶板一直不下沉，实体煤侧就像刚体一样没有变形，无论把参数调的多小，一直不收敛。下面是模型图及flac命令流，请懂得人帮忙分析一下。;-------------------接触面
inter 1 face range x 32.3846 110 y 0 6 z 25.9 26.1
inter 2 face range x 73.9231 110 y 0 6 z 34.9 35.1
model mohr
;-------------砂岩
prop bulk=9.0e9 shear=4.3e9 coh=2.7e6 fric=35 ten=3.8e6 range group 2
prop bulk=9.0e9 shear=4.3e9 coh=2.7e6 fric=35 ten=3.8e6 range group 8
prop bulk=9.0e9 shear=4.3e9 coh=2.7e6 fric=35 ten=3.8e6 range group 12
;-------------粉砂岩
prop bulk=6.4e9 shear=3.2e9 coh=2.1e6 fric=32 ten=2.9e6 range group 7
prop bulk=6.4e9 shear=3.2e9 coh=2.1e6 fric=32 ten=2.9e6 range group 10
;-------------砂质泥岩
prop bulk=5.1e9 shear=2.3e9 coh=1.5e6 fric=34 ten=2.5e6 range group 3
prop bulk=5.1e9 shear=2.3e9 coh=1.5e6 fric=34 ten=2.5e6 range group 9
;-----------泥岩
prop bulk=3.9e9 shear=1.8e9 coh=1.3e6 fric=39 ten=1.9e6 range group 1
prop bulk=3.9e9 shear=1.8e9 coh=1.3e6 fric=39 ten=1.9e6 range group 4
prop bulk=3.9e9 shear=1.8e9 coh=1.3e6 fric=39 ten=1.9e6 range group 11
;----------煤
prop bulk=1.2e9 shear=0.5e9 coh=0.3e6 fric=32 ten=0.7e6 range group 5
prop bulk=1.2e9 shear=0.5e9 coh=0.3e6 fric=32 ten=0.7e6 range group 6
;----------夹层
prop bulk=1.3e9 shear=0.5e9 coh=0.3e6 fric=20 ten=0.6e6 range group 13
prop bulk=1.3e9 shear=0.5e9 coh=0.3e6 fric=20 ten=0.6e6 range group 14
prop bulk=1.3e9 shear=0.5e9 coh=0.3e6 fric=20 ten=0.6e6 range group 15
;-----------密度赋值
;-------------砂岩
ini dens 2700 range group 2
ini dens 2700 range group 8
ini dens 2700 range group 12
;-------------粉砂岩
ini dens 2500 range group 7
ini dens 2500 range group 10
;-------------砂质泥岩
ini dens 2400 range group 3
ini dens 2400 range group 9
;-----------泥岩
ini dens 2200 range group 1
ini dens 2200 range group 4
ini dens 2200 range group 11
;----------煤
ini dens 1400 range group 5
ini dens 1400 range group 6
;----------夹层
ini dens 1800 range group 13
ini dens 1800 range group 14
ini dens 1800 range group 15
;----------------界面参数
int 1 prop kn 300e9 ks 300e9 tens 100e9 煤层上下接触面参数
int 2 prop kn 300e9 ks 300e9 tens 100e9
attach face
;-----------固定边界
fix x range x -0.1 0.1
fix x range x 109.9  110.1
fix y range y -0.1  0.1
fix y range y 5.9  6.1
fix z range z -0.1 0.1
;--------------初始化应力
ini szz  -20.875e6  grad  0 0 2.5e4
ini sxx  -25.05e6   grad  0 0 2.5e4
ini syy  -25.05e6   grad  0 0 2.5e4
apply szz -19.125e6 range z 69.9 70.1
set grav 0 0 -10
;--------------
hist unbal
plot hist 1
set large
solve
save canshu.sav
;------------
restore canshu.sav
initial xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
model null range x 44 110 y 0 6 z 26 35
solve
save caikongqu.sav
;------------
restore caikongqu.sav
model null range group 6
model null range group 15
solve
save zuizhong.sav