在用flac3d模拟充填采场的时候,初始应力平衡一直不收敛,算了几万次最大不平衡应力还有几个兆帕。
希望大家能帮我看下命令流哪里出错了
以下是命令流
new;建立矿体模型
gen zone brick p0=(0,0,0) size 56 40 8 ;底柱
gen zone brick p0=(0,0,48) size 56 40 8 ;顶柱
group dingdizhu ;组合成顶底柱对象
gen zone brick p0=(0,0,8) size 8 40 40 ;矿柱
gen zone brick p0=(16,0,8) size 8 40 40 ;矿柱
gen zone brick p0=(32,0,8) size 8 40 40 ;矿柱
gen zone brick p0=(48,0,8) size 8 40 40 ;矿柱
group kuangzhu range group dingdizhu not ;组合成矿柱对象
;建立充填矿房
;建立人工假底
gen zone brick p0=(8,0,8) size 8 40 2 ;人工假底
gen zone brick p0=(40,0,8) size 8 40 2 ;人工假底
group rengongjiadi range group dingdizhu not group kuangzhu not;组合人工假底对象
;建立胶结充填体
gen zone brick p0=(8,0,10) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,14) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,18) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,22) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,26) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,30) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,34) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,38) size 8 40 3 ;胶结充填
gen zone brick p0=(8,0,42) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,10) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,14) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,18) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,22) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,26) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,30) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,34) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,38) size 8 40 3 ;胶结充填
gen zone brick p0=(40,0,42) size 8 40 6 ;胶结充填
group jiaojie range group dingdizhu not group kuangzhu not group rengongjiadi not ;组合胶结充填对象
;建立浇面充填
gen zone brick p0=(8,0,13) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,17) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,21) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,25) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,29) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,33) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,37) size 8 40 1 ;浇面充填
gen zone brick p0=(8,0,41) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,13) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,17) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,21) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,25) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,29) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,33) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,37) size 8 40 1 ;浇面充填
gen zone brick p0=(40,0,41) size 8 40 1 ;浇面充填
group jiaomian range group dingdizhu not group kuangzhu not group rengongjiadi not group jiaojie not;组合浇面充填对象
;建立初始应力平衡矿房
gen zone brick p0=(24,0,8) size 8 40 40 ;建立初始模拟矿房
;初始应力平衡
;模型参数
mo mohr
pro bulk 14.58e9 she 8.75e9 fric 43.9 coh 3.56e6 ten 2.18e6 dil 3 range group dingdizhu
pro bulk 14.58e9 she 8.75e9 fric 43.9 coh 3.56e6 ten 2.18e6 dil 3 range group kuangzhu
pro bulk 14.58e9 she 8.75e9 fric 43.9 coh 3.56e6 ten 2.18e6 dil 3 range x 24 32 z 8 48
pro bulk 1.83e9 she 0.9e9 fric 32.5 coh 0.55e6 ten 1e6 dil 3 range group rengongjiadi
pro bulk 1.23e9 she 0.57e9 fric 30 coh 0.45e6 ten 0.5e6 dil 3 range group jiaojie
pro bulk 1.67e9 she 0.77e9 fric 30 coh 0.31e6 ten 1.9e6 dil 3 range group jiaomian
;设置边界条件
fix x range x -.1 .1
fix x range x 59.9 56.1
fix y range y -.1 .1
fix y range y -40.1 -39.9
fix z range z -.1 .1
;设置密度
ini density 3800 range group dingdizhu
ini density 3800 range group kuangzhu
ini density 3800 range x 24 32 z 8 48
ini density 2400 range group rengongjiadi
ini density 2200 range group jiaojie
ini density 2200 range group jiaomian
hist unbal
;设置初始应力条件
apply szz=-22e6 range z=(-.1,.1)
apply szz=-21e6 range z=(55.9,56.1)
apply sxx=-18e6 range x=(-.1, .1)
apply sxx=-18e6 range x=(55.9,56.1)
apply syy=-24e6 range y=(-.1,.1)
apply syy=-24e6 range y=(-39.9,-40.1)
set mec fo 50
set grav 0 0 -10
solve
;绘制初始应力图
plot
add bcontour szz
add axes black
show