new
config fluid
set fluid off
gen zon brick p0 0 0 -50 size 80 1 50
gen zon brick p0 50 0 0 size 30 1 20
group jiyan rang x 0 45 a z -25 -50
group louzuo range x 50 80 z -25 20
group manggou range x 45 50 z -10 0
group huitian range x 45 50 z -10 -25
m e
prop bu 3e7 sh 1e7
ini pp 0 grad 0 0 -10e3 ran z 0 -40
ini dens 2000
model fl_iso
prop por 0.5 perm 1e-10
ini fden 1000 ften -1e10
fix pp 0 range x 44.9 45.1 z -10 0
fix pp 0 range x 44.9 50.1 z -9.9 -10.1
model fl_null rang gro louzuo
model null rang gro manggou
fix z ran z -50
fix x ran x 0
fix x ran x 80
fix y
fix x ran x 44.9 45.1 z -10 0
fix x ran x 49.9 50.1 z -10 20
fix z ran x 44.9 50.1 z -9.9 -10.1
set grav 10
his id 1 gp pp 65 0 -25
his id 2 gp pp 65 0 -26
his id 3 gp pp 65 0 -27
his id 4 gp pp 65 0 -28
his id 5 gp pp 50 0 -20
his id 6 gp pp 50 0 -18
his id 7 gp pp 50 0 -16
his id 8 gp pp 68 0 -26
his id 9 gp pp 70 0 -26
his id 10 gp pp 72 0 -26
his id 11 gp pp 74 0 -26
his id 12 gp pp 76 0 -26
his id 13 gp pp 78 0 -26
his id 14 gp pp 79 0 -26
solve
save sy1.sav
new
rest sy1.sav
set fluid on mech off
group jiyan rang z -25 -50
ini fmod 2e3 ften 0.0 rang gro jiyan
ini fmod 2e3 ften 0.0 rang gro huitian
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
prop por 0.5 perm 5e-9 ran gro huitian
app pp 0 grad 0 0 -10e3 ran x -0.1 0.1 z 0 -50
;app vwell 1e-9 range id 1186 a 1187 a 1188 a 1189
app pp 0 ran x 44.9 45.1 z 0 -10
app pp 0 ran x 45 50 z -9.9 -10.1
solve
sav sy2_1.sav
程序如上，设置左侧边界为自由渗透边界，右侧为不透水，
凹槽的地方为孔隙水压是0，现遇到的问题是计算了4万步，不见收敛。
请各位个看看哪里出错了。谢谢
