求大神指导，flac求解大坝模型施加水压力后，求解不收敛

jiangxiaolu0598

建立了一个大坝模型。上游水位高（z方向力作用范围是275-380）高105m ,作用在上游面X的变化范围是0-344,下游水位高（z方向力作用范围是280-300）高20m ,作用在下游面X的变化范围是476-704,整个模型共有7种材料属性，共14005个单元。小弟试着编写如下程序 ：
impgrid Flacmodel_haitang.Flac3D
DEF para  
   MK1=10.21e9
   MK2=3.69e9
   MK3=1.2e9
   MK4=3.33e9
   MK5=7.01e9
   MK6=7.87e9
   MK7=11.91e9
  

G1=8.74e9
G2=1.67e9
G3=0.26e9
G4=1.54e9
G5=3.61e9
G6=4.72e9
G7=8.20e9
F1=56.31
F2=40.36
F3=22.29
F4=32.84
F5=40.79
F6=44.66
F7=50.20


C1=2e6
C2=0.75e6
C3=0.12e6
C4=0.45e6
C5=0.8e6
C6=1e6
C7=1.4e6
D1=2.4e3
D2=1.6e3
D3=2.3e3
D4=2.6e3
D5=2.6e3
D6=2.6e3
D7=2.6e3
END
para  
model mohr
fix x y z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 895.9 896.1
fix y range y -0.1 0.1
fix y range y -30.1 -29.9

pro den D1 bulk MK1 shear G2 coh        C1        fri        F1 range group 1
pro den D2 bulk MK2 shear G2 coh        C2        fri        F2 range group 2
pro den D3 bulk MK3 shear G3 coh        C3        fri        F3 range group 3
pro den D4 bulk MK4 shear G4 coh        C4        fri        F4 range group 4
pro den D5 bulk MK5 shear G5 coh        C5        fri        F5 range group 5
pro den D6 bulk MK6 shear G6 coh        C6        fri        F6 range group 6
pro den D7 bulk MK7 shear G7 coh        C7        fri        F7 range group 7

app nstress 38e3 grad 0 0 -10e3 ran z 275 380 x 0 344 ··上游施加水压力
app nstress 30e3 grad 0 0 -10e3 ran z 280 300 x 476 704 ··下游施加水压力

set grav 0 0 -10

solve
save tiaoshi.sav  最后用实验室的比较快的计算机算了一晚上，540000多步，发现算到某一处算不动了，mach.rato在8e-6左右来回徘徊。 请问各位大神，为什么现在算的收敛不了呢，是不是施加的水压力不对，还是少了对模型的一些处理。期待大神的莅临指导啊！！