请教关于煤层开采后顶板嵌入底板，如何解决？

friendqgp

模拟工作面开挖过程中，超前支承压力对巷道围岩稳定性的影响。
模型尺寸216（长）*100（宽）*63（高），在x 28 32，x 184 188，y 0 100，z 16 19范围内开掘两条巷道，为模拟工作面的上下顺槽，平衡后工作面每5m一开挖，工作面为放顶煤，采高9m。
工作面每开挖5m循环1000步，在循环过程中，顶板向下运动，最终嵌入底板而导致计算终止。
分别尝试了在煤层与底板之间和煤层与顶底板之间设置接触面，最后一次设置的接触面渗透深度为0，均无法阻止顶板向底板嵌入而导致的计算终止，请教各位高手，该如何解决这个问题，确保计算不被终止呢？

这是设置了两个接触面（煤层与顶底板）的效果

这是设置了一个接触面（煤层与底板）的效果
该怎么解决呢？
求版主及高手给个提示，不然做不下去了啊。
附上模型命令：
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设置一个监测点，分别记录巷道顶板中部垂直位移和两帮水平位移，写入文件

friendqgp

顶顶，期待版主或高手来解答

xldcumt

你超前支承怎么做的，把支承的刚度设置大点。我不是高手，我也做这方面的工作，可以相互交流。qq:296116649

zuoruo69

期待各种神秘斑竹

黄花水

（1）参数设置是否合适
（2）考虑充填

takex

你好，请问你顶板嵌入底板的问题解决了么，求指教