各位同行，以下是我施加在滑坡体上的渗流边界。但是发现一个问题就是，渗流边界是施加上了。但是坡体里面没有水的存在，希望哪位能帮我看下。
命令流如下：
new
restore 1.sav;导入模型
fix x range x 331557.8 331558.0
fix x range x 334585.5 334585.7
fix y range y 3493813.9 3493814.1
fix y range y 3498467.8 3498468.0
fix x y z range z 499.9 500.1
model ela
pro bulk 1.5e12 shear 1.2e12 dens 2570 range group 3
pro bulk 3.9e10 shear 2e10 dens 2120 range group 2
pro bulk 1.6e10 shear 7.6e9 dens 1830 range group 1
set grav 0,0,-10
solve
save chushidiyinli.sav
ini xvel=0
ini yvel=0
ini zvel=0
ini xdisp=0
ini ydisp=0
ini zdisp=0
config fluid ;配置流体计算参数
set fluid on ;设置流体计算为启用
set mech on ;设置力学计算为启用
;========================
;流体模型参数
;=========================
model fl_iso ;设置材料为各项同性的等效孔隙介质
prop perm 8.68e-20 poro 0.4 ran group 3 ;设置粉砂质泥岩岩层的岩体渗透率、孔隙率
prop perm 3.12e-6 poro 0.45 ran group 2 ;设置松散堆积体的渗透率、孔隙率
prop perm 1.33e-6 poro 0.35 ran group 1
set fl biot off ;设置材料的比奥固结属性
ini fdensity 1e3 ;初始化流体密度，水的密度1000kg/立方米
ini sat 0 ran group 3 ;初始化流体饱和度，设为0.3
ini sat 0.27 ran group 2
ini sat 0.32 ran group 1
ini fmod 2e9 ftens -1e10 ;初始化流体弹性模量和抗拉强度，水的抗拉强度设为负值
def derive
;define bulk(K)(单位:N,kg,m,pa)
b_mod1=E_mod1/(3*(1.0-2*p_ratio1)) ;3
b_mod2=E_mod2/(3*(1.0-2*p_ratio2)) ;2
b_mod3=E_mod3/(3*(1.0-2*p_ratio3)) ;1
;define shear(G)
s_mod1=E_mod1/(2.0*(1+p_ratio1)) ;3
s_mod2=E_mod2/(2.0*(1+p_ratio2)) ;2
s_mod3=E_mod3/(2.0*(1+p_ratio3)) ;1
end
set E_mod1 =2.7e+11 p_ratio1 =0.21 ;3
set E_mod2 =5.2e+9 p_ratio2 =0.28 ;2
set E_mod3 =2e+9 p_ratio3 =0.30 ;1
derive
model mohr range group 3
model mohr range group 2
model mohr range group 1
prop bulk b_mod1 shear s_mod1 range group 3
prop bulk b_mod2 shear s_mod2 range group 2
prop bulk b_mod3 shear s_mod3 range group 1
prop fric 39 coh 1e+8 tens 1.5e+8 range group 3
prop fric 25 coh 1.2e+4 range group 2
prop fric 25 coh 1.0e+4 range group 1
ini dens 2570 range group 3
ini dens 2120 range group 2
ini dens 1830 range group 1
model fl_null range group 3
apply pwell 0.001 range group 2
apply pwell 0.03 range group 1
set grav 0,0,-10
solve
save jiangyu2.sav
