求助：采空区顶板怎么会压入顶板 如何调整

chongling

采出1.5m煤层之后，采空区顶板下沉将近2m，如何控制避免这种情况的出现 不知道哪儿高手能指教一番？！谢谢！！

sch

没施加初始地应力，并检查力学参数

li_2001

没施加初始地应力，并检查力学参数

cumtsp

顶板设接触

chongling

4# cumtsp
请教，顶板如何设接触，能具体说一下么？谢谢

chongling

2# sch

初始地应力的施加应该没问题吧 我把地应力施加的代码贴出来，实在地下370-420m的模型
ini szz=-10.6e6 grad 0,0,2.5e4
ini sxx=-7.95e6 grad 0,0,1.875e4
ini syy=-7.95e6 grad 0,0,1.875e4
apply szz=-9.25e6 range z  39.9  40.1
apply szz=-10.6e6 range z -15.1 -14.9
apply sxx=-7.95e6  grad 0,0,1.875e4  range x  -30.1 -29.9
apply sxx=-7.95e6  grad 0,0,1.875e4 range x  59.9 60.1
apply syy=-7.95e6  grad 0,0,1.875e4  range y  -0.1  0.1
apply syy=-7.95e6   grad 0,0,1.875e4 range y  99.9  100.1
请教 有问题么
另外，力学参数于什么问题才会出现主题说的这种情况啊  
谢谢指教！

cloud328

可能参数太低。全部代码贴出来。

chongling

我是刚开始学习 请各位看看里面有什么错误 谢谢
new
plot set title text ""
set logfile mn1.log
set log on
plot add axes red

;巷道及其上下部分
gen zone brick size 4 20 2 p0 -1.6 0 0 p1 1.6 0 0 p2 -1.6 100 0 p3 -1.6 0 1.3 ;巷道
gen zone brick size 4 20 4 p0 -1.6 0 1.3 p1 1.6 0 1.3 p2 -1.6 100 1.3 p3 -1.6 0 6 ;巷道顶部
gen zone brick size 4 20 5 p0 -1.6 0 6 p1 1.6 0 6 p2 -1.6 100 6 p3 -1.6 0 11 ;group hl老顶
gen zone brick size 4 20 7 p0 -1.6 0 11 p1 1.6 0 11 p2 -1.6 100 11 p3 -1.6 0 40 ;group hs上覆岩层
gen zone brick size 4 20 2 p0 -1.6 0 -1 p1 1.6 0 -1 p2 -1.6 100 -1 p3 -1.6 0 0 ;group hzd直接底
gen zone brick size 4 20 5 p0 -1.6 0 -5.8 p1 1.6 0 -5.8 p2 -1.6 100 -5.8 p3 -1.6 0 -1 ;group hld老底
gen zone brick size 4 20 5 p0 -1.6 0 -15 p1 1.6 0 -15 p2 -1.6 100 -15 p3 -1.6 0 -5.8 ;group hx下部岩层

;充填体及其上下部分
gen zone brick size 6 20 2 p0 1.6 0 0 p1 4.6 0 0 p2 1.6 100 0 p3 1.6 0 1.3 ;group ctt充填体
gen zone brick size 6 20 4 p0 1.6 0 1.3 p1 4.6 0 1.3 p2 1.6 100 1.3 p3 1.6 0 6 ;group ctz直接顶部
gen zone brick size 6 20 5 p0 1.6 0 6 p1 4.6 0 6 p2 1.6 100 6 p3 1.6 0 11 ;group ctl充填老顶
gen zone brick size 6 20 7 p0 1.6 0 11 p1 4.6 0 11 p2 1.6 100 11 p3 1.6 0 40 ;group cts上覆岩层
gen zone brick size 6 20 2 p0 1.6 0 -1 p1 4.6 0 -1 p2 1.6 100 -1 p3 1.6 0 0 ;group ctzd直接底
gen zone brick size 6 20 5 p0 1.6 0 -5.8 p1 4.6 0 -5.8 p2 1.6 100 -5.8 p3 1.6 0 -1 ;group ctld老底
gen zone brick size 6 20 5 p0 1.6 0 -15 p1 4.6 0 -15 p2 1.6 100 -15 p3 1.6 0 -5.8 ;group ctx下部岩层
;工作面及其上下部分
gen zone brick size 20 20 2 p0 4.6 0 0 p1 60 0 0 p2 4.6 100 0 p3 4.6 0 1.3 ;group m 煤体
gen zone brick size 20 20 4 p0 4.6 0 1.3 p1 60 0 1.3 p2 4.6 100 1.3 p3 4.6 0 6 ;group mz直接顶
gen zone brick size 20 20 5 p0 4.6 0 6 p1 60 0 6 p2 4.6 100 6 p3 4.6 0 11 ;group ml老顶
gen zone brick size 20 20 7 p0 4.6 0 11 p1 60 0 11 p2 4.6 100 11 p3 4.6 0 40 ;group ms上覆
gen zone brick size 20 20 2 p0 4.6 0 -1 p1 60 0 -1 p2 4.6 100 -1 p3 4.6 0 0 ;group mzd直接底
gen zone brick size 20 20 5 p0 4.6 0 -5.8 p1 60 0 -5.8 p2 4.6 100 -5.8 p3 4.6 0 -1 ;group mld老底
gen zone brick size 20 20 5 p0 4.6 0 -15 p1 60 0 -15 p2 4.6 100 -15 p3 4.6 0 -5.8 ;group mx下部岩层
;下区段工作面部分
gen zone brick size 15 20 2 p0 -30 0 0 p1 -1.6 0 0 p2 -30 100 0 p3 -30 0 1.3 ;group m 煤体
gen zone brick size 15 20 4 p0 -30 0 1.3 p1 -1.6 0 1.3 p2 -30 100 1.3 p3 -30 0 6 ;group mz直接顶
gen zone brick size 15 20 5 p0 -30 0 6 p1 -1.6 0 6 p2 -30 100 6 p3 -30 0 11 ;group ml老顶
gen zone brick size 15 10 7 p0 -30 0 11 p1 -1.6 0 11 p2 -30 100 11 p3 -30 0 40 ;group ms上覆
gen zone brick size 15 20 2 p0 -30 0 -1 p1 -1.6 0 -1 p2 -30 100 -1 p3 -30 0 0 ;group mzd直接底
gen zone brick size 15 20 5 p0 -30 0 -5.8 p1 -1.6 0 -5.8 p2 -30 100 -5.8 p3 -30 0 -1 ;group mld老底
gen zone brick size 15 10 5 p0 -30 0 -15 p1 -1.6 0 -15 p2 -30 100 -15 p3 -30 0 -5.8 ;group mx下部岩层
attach face
plot block state
save moxingjihe1.sav
mod mohr
;煤层
prop bulk=0.4e9     shear=0.133e9      fric=18   coh=2.4e6   tens=1.0e6   range z 0 1.3
ini density 1300 range z 0 1.3
;直接顶 泥岩
prop bulk=0.8e9    shear=0.4e9  fric=25    coh=6e6   tens=5e6    range z 1.3 6
ini density 2500 range z 1.3 6
;老顶 中砂岩
prop bulk=0.55e9    shear=0.32e9  fric=32    coh=7.8e6   tens=9.2e6   range z 6 11
ini density 2550 range z 6 11
;上覆岩层 砂岩
prop bulk=0.776e9    shear=0.558e9  fric=32    coh=8.0e6   tens=9.0e6   range z 11 40
ini density 2650 range z 11 40
;直接底 粉砂岩
prop bulk=0.55e9    shear=0.32e9  fric=32    coh=7.8e6   tens=9.2e6   range z -1 0
ini density 2550 range z -1 0
;老底 细砂岩左部
prop bulk=0.536e9     shear=0.369e9      fric=32   coh=8.7e6   tens=13e6  range   z -5.8 -1
ini density 2500 range z -5.8 -1
;下部岩层 砂岩
prop bulk=0.776e9    shear=0.558e9  fric=32    coh=8.0e6   tens=9.0e6   range z -15 -5.8
ini density 2650 range z -15 -5.8
set grav 0 0 -10
fix x range x -60.1 -59.9
fix x range x 59.9  60.1
fix z range z -15.1 -14.9
fix y range y -0.1 0.1
fix y range y 99.9 100.1
;施加应力边界条件
ini szz=-10.6e6 grad 0,0,2.5e4
ini sxx=-7.95e6 grad 0,0,1.875e4
ini syy=-7.95e6 grad 0,0,1.875e4
apply szz=-9.25e6 range z  39.9  40.1
apply szz=-10.6e6 range z -15.1 -14.9
apply sxx=-7.95e6  grad 0,0,1.875e4  range x  -30.1 -29.9
apply sxx=-7.95e6  grad 0,0,1.875e4 range x  59.9 60.1
apply syy=-7.95e6  grad 0,0,1.875e4  range y  -0.1  0.1
apply syy=-7.95e6   grad 0,0,1.875e4 range y  99.9  100.1
plot block state
solve
save 初始平衡.sav
;位移初始化
ini xdis=0
ini ydis=0
ini zdis=0
;挖巷道
model null range x -1.6 1.6 y 0 100 z -1 1.3
;定义二维数组函数addd，便于确定顶锚杆的施加位置，顶4根锚杆，端部锚杆斜打
def addd
array xz(2,4) zz(2,4)
xz(1,1)=-1.2
xz(1,2)=-0.4
xz(1,3)=0.4
xz(1,4)=1.2

zz(1,1)=1.3
zz(1,2)=1.3
zz(1,3)=1.3
zz(1,4)=1.3


xz(2,1)=-1.9
xz(2,2)=-0.4
xz(2,3)=0.4
xz(2,4)=1.9

zz(2,1)=2.84
zz(2,2)=3
zz(2,3)=3
zz(2,4)=2.84

end

def add1
xz1=xz(1,index)
zz1=zz(1,index)

xz2=xz(2,index)
zz2=zz(2,index)
end

addd
;定义函数inscab，采用循环模式模拟计算模型内锚杆的排数
def inscab

y1=0.2
loop m(1,124)
loop index(1,4)
add1
command
;设定锚杆的位置
sel cable id 2 begin xz1 y1 zz1 end xz2 y1 zz2 nseg 4
;设定锚杆的物理力学参数
sel cable prop xca 3.8e-4 e 2e11 yt 1.56e5 gr_k 1.28e10 gr_c 2.25e5 range id 2
end_command
end_loop
;排距0.8m
y1=y1+0.8
end_loop

end
;定义addd1数组函数用来确定帮锚杆位置，每帮锚杆3根，帮上上下两根锚杆是倾斜打的
def addd1
array xs(2,3) zs(2,3) xy(2,3) zy(2,3)

xs(1,1)=-1.6
xs(1,2)=-1.6
xs(1,3)=-1.6

zs(1,1)=-0.5
zs(1,2)=0.3
zs(1,3)=1.1

xs(2,1)=-3.3
xs(2,2)=-3.3
xs(2,3)=-3.14

zs(2,1)=-0.5
zs(2,2)=0.3
zs(2,3)=1.8

xy(1,1)=1.6
xy(1,2)=1.6
xy(1,3)=1.6

zy(1,1)=-0.5
zy(1,2)=0.3
zy(1,3)=1.1

xy(2,1)=3.3
xy(2,2)=3.3
xy(2,3)=3.14

zy(2,1)=-0.5
zy(2,2)=0.3
zy(2,3)=1.8

end

def add2

xs1=xs(1,index)
zs1=zs(1,index)

xs2=xs(2,index)
zs2=zs(2,index)

xy1=xy(1,index)
xy2=xy(2,index)

zy1=zy(1,index)
zy2=zy(2,index)

end
addd1
def inscab1
y0=0.2
loop m(1,124)
loop index(1,3)
add2
command
;施加帮锚杆单元和设定物理力学参数
sel cable id 1 begin xs1 y0 zs1 end xs2 y0 zs2 nseg 3
sel cable id 1 begin xy1 y0 zy1 end xy2 y0 zy2 nseg 3  
sel cable prop xca 2.54e-4 e 2e11 yt 1.04e5 gr_k 1.28e10 gr_c 2.25e5 range id 1
end_command
end_loop

y0=y0+0.8
end_loop
end

inscab
inscab1
;2.5m处巷道变形情况依次为顶板 底板
hist gp zdis 0 2.5 1.3      
hist gp zdis 0 2.5 -1
;50m 处巷道变形情况依次为顶板 底板
hist gp zdis 0 50 1.3
hist gp zdis 0 50 -1
solve
save 巷道成型.sav
;位移初始化
ini xdis=0
ini ydis=0
ini zdis=0
;挖巷道
model null range x -1.6 1.6 y 0 100 z -1 1.3

;定义二维数组函数addd，便于确定顶锚杆的施加位置，顶4根锚杆，端部锚杆斜打
def addd
array xz(2,4) zz(2,4)
xz(1,1)=-1.2
xz(1,2)=-0.4
xz(1,3)=0.4
xz(1,4)=1.2

zz(1,1)=1.3
zz(1,2)=1.3
zz(1,3)=1.3
zz(1,4)=1.3


xz(2,1)=-1.9
xz(2,2)=-0.4
xz(2,3)=0.4
xz(2,4)=1.9

zz(2,1)=2.84
zz(2,2)=3
zz(2,3)=3
zz(2,4)=2.84

end

def add1
xz1=xz(1,index)
zz1=zz(1,index)

xz2=xz(2,index)
zz2=zz(2,index)
end

addd
;定义函数inscab，采用循环模式模拟计算模型内锚杆的排数
def inscab

y1=0.2
loop m(1,124)
loop index(1,4)
add1
command
;设定锚杆的位置
sel cable id 2 begin xz1 y1 zz1 end xz2 y1 zz2 nseg 4
;设定锚杆的物理力学参数
sel cable prop xca 3.8e-4 e 2e11 yt 1.56e5 gr_k 1.28e10 gr_c 2.25e5 range id 2
end_command
end_loop
;排距0.8m
y1=y1+0.8
end_loop

end
;定义addd1数组函数用来确定帮锚杆位置，每帮锚杆3根，帮上上下两根锚杆是倾斜打的
def addd1
array xs(2,3) zs(2,3) xy(2,3) zy(2,3)

xs(1,1)=-1.6
xs(1,2)=-1.6
xs(1,3)=-1.6

zs(1,1)=-0.5
zs(1,2)=0.3
zs(1,3)=1.1

xs(2,1)=-3.3
xs(2,2)=-3.3
xs(2,3)=-3.14

zs(2,1)=-0.5
zs(2,2)=0.3
zs(2,3)=1.8

xy(1,1)=1.6
xy(1,2)=1.6
xy(1,3)=1.6

zy(1,1)=-0.5
zy(1,2)=0.3
zy(1,3)=1.1

xy(2,1)=3.3
xy(2,2)=3.3
xy(2,3)=3.14

zy(2,1)=-0.5
zy(2,2)=0.3
zy(2,3)=1.8

end

def add2

xs1=xs(1,index)
zs1=zs(1,index)

xs2=xs(2,index)
zs2=zs(2,index)

xy1=xy(1,index)
xy2=xy(2,index)

zy1=zy(1,index)
zy2=zy(2,index)

end
addd1
def inscab1
y0=0.2
loop m(1,124)
loop index(1,3)
add2
command
;施加帮锚杆单元和设定物理力学参数
sel cable id 1 begin xs1 y0 zs1 end xs2 y0 zs2 nseg 3
sel cable id 1 begin xy1 y0 zy1 end xy2 y0 zy2 nseg 3  
sel cable prop xca 2.54e-4 e 2e11 yt 1.04e5 gr_k 1.28e10 gr_c 2.25e5 range id 1
end_command
end_loop

y0=y0+0.8
end_loop
end

inscab
inscab1
;2.5m处巷道变形情况依次为顶板 底板
hist gp zdis 0 2.5 1.3      
hist gp zdis 0 2.5 -1
;50m 处巷道变形情况依次为顶板 底板
hist gp zdis 0 50 1.3
hist gp zdis 0 50 -1
solve
save 巷道成型.sav

hist reset
set logfile kaiwatanmo.log
set log on
;位移初始化
ini xdis=0
ini ydis=0
ini zdis=0
;提前开帮0-5m
ini sxx=0 range x 1.6 3.1 y 0 5 z 0 1.3
ini syy=0 range x 1.6 3.1 y 0 5 z 0 1.3
ini szz=0 range x 1.6 3.1 y 0 5 z 0 1.3
model null range x 1.6 3.1 y 0 5 z 0 1.3



;充填
model mohr range x 1.6 3.1 y 0 5 z 0 1.3

prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6  range x 1.6 3.1 y 0 5 z 0 1.3
ini density 3000 range x 1.6 3.1 y 0 5 z 0 1.3

step 50
;监测位移
hist n=5
;顶板下沉、底臌变形
hist gp zdisp 0 2.5 1.3 ;1 顶沉
hist gp zdisp 0 2.5 -1  ;2 底臌
;四角位移 左上角 右上角 左下角 右下角
hist gp zdis -1.6 2.5 1.3  ;3 左上角
hist gp xdis -1.6 2.5 1.3  ;4 左上角
hist gp zdis 1.6 2.5 1.3   ;5 右上角
hist gp xdis 1.6 2.5 1.3   ;6 右上角
hist gp zdis -1.6 2.5 -1   ;7 左下角
hist gp xdis -1.6 2.5 -1   ;8 左下角
hist gp zdis 1.6 2.5 -1    ;9 右下角
hist gp xdis 1.6 2.5 -1    ;10 右下角
;两帮移近量 左 右
hist gp xdisp -1.6 2.5 0   ;11 左帮
hist gp xdisp -1.6 2.5 0.6  ;12 左帮
hist gp xdisp 1.6 2.5 0    ;13 右帮
hist gp xdisp 1.6 2.5 0.6  ;14 右帮

;监测应力
;充填体应力
hist zone szz 1.7 2.5  1.3                   ;15      内侧顶压
hist zone szz  3.0 2.5 1.3                   ;16     外侧顶压
hist zone szz 1.7 2.5  0                   ;17      内侧底压
hist zone szz  3.0 2.5 0                   ;18     外侧底压
;煤体应力监测                              
hist zone szz -1.7 2.5  1.3                  ;19      近煤壁处顶压力
hist zone szz -1.7 2.5  0                   ;20      近煤壁处底压
hist zone szz -3.7 2.5  1.3                   ;21      2m内部压力
hist zone szz  -6.6 2.5 1.3                   ;22     5m内部压力
hist zone szz  -10.6 2.5 1.3                   ;23     9m内部压力
hist zone szz 0 2.5 1.4                        ;24 顶板压力
hist zone szz 0 2.5 -0.9                        ;25 底板压力
hist gp zdis 0 35 1.3                        ;26 35m处顶板位移
hist gp zdis 0 35 -1                          ;27 35m处底板位移
hist zone szz 1.7 35 1.3                       ;28 35m处工作面侧垂直压力
hist zone szz -1.7 35 1.3                     ;29 35m煤壁侧垂直压力
;开挖上工作面0-5m煤层
model null range  x 3.1 60 y 0 5 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=8e6   tens=6e6    range x 3.1 60 y 0 5 z 1.3 6
ini density 2500 range x 3.1 60 y 0 5 z 1.3 6
set large
solve ratio 4e-4
save 开挖充填5m.sav

;提前开帮5-10m
ini sxx=0 range x 1.6 3.1 y 5 10 z 0 1.3
ini syy=0 range x 1.6 3.1 y 5 10 z 0 1.3
ini szz=0 range x 1.6 3.1 y 5 10 z 0 1.3
model null range x 1.6 3.1 y 5 10 z 0 1.3

model mohr range x 1.6 3.1 y 5 10 z 0 1.3

prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6    range x 1.6 3.1 y 5 10 z 0 1.3
ini density 3000 range x 1.6 3.1 y 5 10 z 0 1.3


step 50
;开挖上工作面5-10煤层
model null range  x 3.1 60 y 5 10 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=6e6   tens=2.5e6    range x 3.1 60 y 5 10 z 1.3 6
ini density 2500 range x 3.1 60 y 5 10 z 1.3 6
solve ratio 4e-4
save 开挖充填10m.sav
;提前开帮10-15m
ini sxx=0 range x 1.6 3.1 y 10 15 z 0 1.3
ini syy=0 range x 1.6 3.1 y 10 15 z 0 1.3
ini szz=0 range x 1.6 3.1 y 10 15 z 0 1.3
model null range x 1.6 3.1 y 10 15 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 10 15 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6    range x 1.6 3.1 y 10 15 z 0 1.3
ini density 3000 range x 1.6 3.1 y 10 15 z 0 1.3
step 50
;开挖上工作面10-15煤层
model null range  x 3.1 60 y 10 15 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 10 15 z 1.3 6
ini density 2500 range x 3.1 60 y 10 15 z 1.3 6
solve ratio 4e-4
save 开挖充填15m.sav




;提前开帮15-20m
ini sxx=0 range x 1.6 3.1 y 15 20 z 0 1.3
ini syy=0 range x 1.6 3.1 y 15 20 z 0 1.3
ini szz=0 range x 1.6 3.1 y 15 20 z 0 1.3
model null range x 1.6 3.1 y 15 20 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 15 20 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6   range x 1.6 3.1 y 15 20 z 0 1.3
ini density 3000 range x 1.6 3.1 y 15 20 z 0 1.3

step 50
;开挖上工作面15-20煤层
model null range  x 3.1 60 y 15 20 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 15 20 z 1.3 6
ini density 2500 range x 3.1 60 y 15 20 z 1.3 6
solve ratio 4e-4
save 开挖充填20m.sav

;提前开帮20-25m
ini sxx=0 range x 1.6 3.1 y 20 25 z 0 1.3
ini syy=0 range x 1.6 3.1 y 20 25 z 0 1.3
ini szz=0 range x 1.6 3.1 y 20 25 z 0 1.3
model null range x 1.6 3.1 y 20 25 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 20 25 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6    range x 1.6 3.1 y 20 25 z 0 1.3
ini density 3000 range x 1.6 3.1 y 20 25 z 0 1.3
step 50
;开挖上工作面20-25煤层
model null range  x 3.1 60 y 20 25 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 20 25 z 1.3 6
ini density 2500 range x 3.1 60 y 20 25 z 1.3 6
solve ratio 4e-4
save 开挖充填25m.sav

;提前开帮25-30m
ini sxx=0 range x 1.6 3.1 y 25 30 z 0 1.3
ini syy=0 range x 1.6 3.1 y 25 30 z 0 1.3
ini szz=0 range x 1.6 3.1 y 25 30 z 0 1.3
model null range x 1.6 3.1 y 25 30 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 25 30 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6   range x 1.6 3.1 y 25 30 z 0 1.3
ini density 3000 range x 1.6 3.1 y 25 30 z 0 1.3

step 50
;开挖上工作面25-30煤层
model null range  x 3.1 60 y 25 30 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 25 30 z 1.3 6
ini density 2500 range x 3.1 60 y 25 30 z 1.3 6
solve ratio 4e-4
save 开挖充填30m.sav






;提前开帮30-35m
ini sxx=0 range x 1.6 3.1 y 30 35 z 0 1.3
ini syy=0 range x 1.6 3.1 y 30 35 z 0 1.3
ini szz=0 range x 1.6 3.1 y 30 35 z 0 1.3
model null range x 1.6 3.1 y 30 35 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 30 35 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6     range x 1.6 3.1 y 30 35 z 0 1.3
ini density 3000 range x 1.6 3.1 y 30 35 z 0 1.3

step 50
;开挖上工作面30-35煤层
model null range  x 3.1 60 y 30 35 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 30 35 z 1.3 6
ini density 2500 range x 3.1 60 y 30 35 z 1.3 6
solve ratio 4e-4
save 开挖充填35m.sav

;提前开帮35-40m
ini sxx=0 range x 1.6 3.1 y 35 40 z 0 1.3
ini syy=0 range x 1.6 3.1 y 35 40 z 0 1.3
ini szz=0 range x 1.6 3.1 y 35 40 z 0 1.3
model null range x 1.6 3.1 y 35 40 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 35 40 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6    range x 1.6 3.1 y 35 40 z 0 1.3
ini density 3000 range x 1.6 3.1 y 35 40 z 0 1.3

step 50
;开挖上工作面35-40煤层
model null range  x 3.1 60 y 35 40 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 35 40 z 1.3 6
ini density 2500 range x 3.1 60 y 35 40 z 1.3 6
solve ratio 4e-4
save 开挖充填40m.sav

;提前开帮40-45m
ini sxx=0 range x 1.6 3.1 y 40 45 z 0 1.3
ini syy=0 range x 1.6 3.1 y 40 45 z 0 1.3
ini szz=0 range x 1.6 3.1 y 40 45 z 0 1.3
model null range  x 1.6 3.1 y 40 45 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 40 45 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6    range x 1.6 3.1 y 40 45 z 0 1.3
ini density 3000 range x 1.6 3.1 y 40 45 z 0 1.3

step 50
;开挖上工作面40-45煤层
model null range  x 3.1 60 y 40 45 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 40 45 z 1.3 6
ini density 2500 range x 3.1 60 y 40 45 z 1.3 6
solve ratio 4e-4
save 开挖充填45m.sav


;提前开帮45-50m
ini sxx=0 range x 1.6 3.1 y 45 50 z 0 1.3
ini syy=0 range x 1.6 3.1 y 45 50 z 0 1.3
ini szz=0 range x 1.6 3.1 y 45 50 z 0 1.3
model null range x 1.6 3.1 y 45 50 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 45 50 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6    range x 1.6 3.1 y 45 50 z 0 1.3
ini density 3000 range x 1.6 3.1 y 45 50 z 0 1.3

step 50
;开挖上工作面45-50煤层
model null range  x 3.1 60 y 45 50 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 45 50 z 1.3 6
ini density 2500 range x 3.1 60 y 45 50 z 1.3 6
solve ratio 4e-4
save 开挖充填50m.sav



;提前开帮50-55m
ini sxx=0 range x 1.6 3.1 y 50 55 z 0 1.3
ini syy=0 range x 1.6 3.1 y 50 55 z 0 1.3
ini szz=0 range x 1.6 3.1 y 50 55 z 0 1.3
model null range x 1.6 3.1 y 50 55 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 50 55 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6     range x 1.6 3.1 y 50 55 z 0 1.3
ini density 3000 range x 1.6 3.1 y 50 55 z 0 1.3

step 50
;开挖上工作面50-55煤层
model null range  x 3.1 60 y 50 55 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 50 55 z 1.3 6
ini density 2500 range x 3.1 60 y 50 55 z 1.3 6
solve ratio 4e-4
save 开挖充填55m.sav

;提前开帮55-60m
ini sxx=0 range x 1.6 3.1 y 55 60 z 0 1.3
ini syy=0 range x 1.6 3.1 y 55 60 z 0 1.3
ini szz=0 range x 1.6 3.1 y 55 60 z 0 1.3
model null range x 1.6 3.1 y 55 60 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 55 60 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6    range x 1.6 3.1 y 55 60 z 0 1.3
ini density 3000 range x 1.6 3.1 y 55 60 z 0 1.3
;开挖上工作面55-60煤层
model null range  x 3.1 60 y 55 60 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 55 60 z 1.3 6
ini density 2500 range x 3.1 60 y 55 60 z 1.3 6

solve ratio 4e-4
save 开挖充填60m.sav

;提前开帮60-65m
ini sxx=0 range x 1.6 3.1 y 60 65 z 0 1.3
ini syy=0 range x 1.6 3.1 y 60 65 z 0 1.3
ini szz=0 range x 1.6 3.1 y 60 65 z 0 1.3
model null range x 1.6 3.1 y 60 65 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 60 65 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6     range x 1.6 3.1 y 60 65 z 0 1.3
ini density 3000 range x 1.6 3.1 y 60 65 z 0 1.3

step 50
;开挖上工作面60-65煤层
model null range  x 3.1 60 y 60 65 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 60 65 z 1.3 6
ini density 2500 range x 3.1 60 y 60 65 z 1.3 6
solve ratio 4e-4
save 开挖充填65m.sav

;提前开帮65-70m
ini sxx=0 range x 1.6 3.1 y 65 70 z 0 1.3
ini syy=0 range x 1.6 3.1 y 65 70 z 0 1.3
ini szz=0 range x 1.6 3.1 y 65 70 z 0 1.3
model null range x 1.6 3.1 y 65 70 z 0 1.3

;充填
model mohr range x 1.6 3.1 y 65 70 z 0 1.3
prop bulk=1e9  shear=0.8e9    fric=34  coh=10e6   tens=12e6    range x 1.6 3.1 y 65 70 z 0 1.3
ini density 3000 range x 1.6 3.1 y 65 70 z 0 1.3

step 50
;开挖上工作面65-70煤层
model null range  x 3.1 60 y 65 70 z 0 1.3
;采空区处理
prop bulk=0.4e9    shear=0.2e9  fric=25    coh=4e6   tens=3e6    range x 3.1 60 y 65 70 z 1.3 6
ini density 2500 range x 3.1 60 y 65 70 z 1.3 6
solve ratio 4e-4
save 开挖充填70m.sav

sch

你的错误出在：模型初始不要施加应力边界条件，应力边界条件可以在工况中施加；如果有构造应力情况下，可施加初始应力！

sch

还有：模型开始自重应力生成，应该是弹性模型条件下。

zhangruian

我觉得你的物理参数取的太小！！采空区的沉降量没有你计算的那么大！！安心调参吧！！把bu和sh值先调大看看，然后再看塑性区的分布理想吗？？再决定调内聚力，抗拉强度，内摩擦角等！！！

tg0215

抗拉强度太大了吧？？？弹性模量也太小了吧？

chongling

谢谢楼上诸位 参数是沿用了以前别人做过的模拟 因为是一个地方的项目 我也觉得偏小 在修改一下
再次感谢sch tg0215 zhangruian cloud328等

hnlglyx

我也遇到同样的问题