模拟工作面开挖过程中，超前支承压力对巷道围岩稳定性的影响。
模型尺寸216（长）*100（宽）*63（高），在x 28 32，x 184 188，y 0 100，z 16 19范围内开掘两条巷道，为模拟工作面的上下顺槽，平衡后工作面每5m一开挖，工作面为放顶煤，采高9m。
工作面每开挖5m循环1000步，在循环过程中，顶板向下运动，最终嵌入底板而导致计算终止。
分别尝试了在煤层与底板之间和煤层与顶底板之间设置接触面，最后一次设置的接触面渗透深度为0，均无法阻止顶板向底板嵌入而导致的计算终止，请教各位高手，该如何解决这个问题，确保计算不被终止呢？
18 分钟前 上传
下载附件 (26.09 KB) 这是设置了两个接触面（煤层与顶底板）的效果
18 分钟前 上传
下载附件 (29.06 KB) 这是设置了一个接触面（煤层与底板）的效果
该怎么解决呢？
求版主及高手给个提示，不然做不下去了啊。
附上模型命令：
rest 模型5-无支护巷道开挖.sav
set large
hist reset
pause
define face_cave
aa = 0.0
bb = 5.0
loop m (1,10)
command
model null range x 32 184 y aa bb z 16 25
end_command
command
inter 1 delete
inter 1 face range plane norm (0,0,1) origin (30,0,16) dist 0.3
;inter 2 face range plane norm (0,0,-1) origin (30,0,25) dist 0.3
inter 1 ctol 0.0
;inter 2 maxedge 4
inter 1 prop kn 20e11 ks 20e11 coh 10e8 fri 30；接触面设置的参数，这个参数应该很大了吧？
;inter 2 prop kn 20e9 ks 20e8 coh 10e7 fri 30
step 1000
end_command
aa = aa + 5.0
bb = bb + 5.0
end_loop
end
define hist_zdis
a = 1
b = 0.0
loop n (1,19)
command
hist id a nstep 40 gp zdis 30 b 19
end_command
a = a + 1
b = b + 5.0
end_loop
end
define hist_left_xdis
c = 21
d = 0.0
loop n (1,19)
command
hist id c nstep 40 gp zdis 28 d 17.5
end_command
c = c + 1
d = d + 5.0
end_loop
end
define hist_right_xdis
e = 41
f = 0.0
loop n (1,19)
command
hist id e nstep 40 gp zdis 32 f 17.5
end_command
e = e + 1
f = f + 5.0
end_loop
end
define hist_roof_smin
g = 61
h = 0.0
loop n (1,19)
command
hist id g nstep 40 zone smin 30 h 20
end_command
g = g + 1
h = h + 5.0
end_loop
end
define hist_roof_smax
i = 81
j = 0.0
loop n (1,19)
command
hist id i nstep 40 zone smax 30 j 20
end_command
i = i + 1
j = j + 5.0
end_loop
end
define histwrite_zdis
k = 1
loop n (1,19)
command
hist write k file zdis_1接触面.txt
end_command
k = k + 1
end_loop
end
define histwrite_left_xdis
l = 21
loop n (1,19)
command
hist write l file left_xdis_1接触面.txt
end_command
l = l + 1
end_loop
end
define histwrite_right_xdis
p = 41
loop n (1,19)
command
hist write p file right_xdis_1接触面.txt
end_command
p = p + 1
end_loop
end
define histwrite_roof_smin
q = 61
loop n (1,19)
command
hist write q file roof_smin_1接触面.txt
end_command
q = q + 1
end_loop
end
define histwrite_roof_smax
r = 81
loop n (1,19)
command
hist write r file roof_smax_1接触面.txt
end_command
r = r + 1
end_loop
end
set grav 0 0 -10
set large
hist_zdis
hist_left_xdis
hist_right_xdis
hist_roof_smin
hist_roof_smax
face_cave
solve
histwrite_zdis
histwrite_left_xdis
histwrite_right_xdis
histwrite_roof_smin
histwrite_roof_smax
save 模型5-无支护巷道采场开挖-1接触面.sav
return
;采场每5m开挖一次，每5m设置一个监测点，分别记录巷道顶板中部垂直位移和两帮水平位移，写入文件
