在做一个土工膜斜墙堆石坝的渗流例子，在分级加载的基础上做单渗流，命令流如下。
restore 加载分级编号17.sav
config fluid
set mech off
model fl_iso
set fl biot off
prop por 0.2 perm 1.02e-8 range group rockfill ；堆石体赋e和k
prop por 0.2 perm 5.1e-7 range group drain ；排水体赋e和k
prop por 0.2 perm 5.0e-9 range group cushion ；垫层赋e和k
prop por 0.5 perm 1.02e-8 range group soilddown ；地基1赋e和k
prop por 0.5 perm 1.02e-9 range group soildtop ；地基2赋e和k
prop por 0.2 perm 1.02e-16 range group panel ；土工膜等效后的土体赋e和k
ini sat 1.0 range z 0 125 ；水位以下包括地层赋予初始饱和度
ini fdens 1000 ften 1e-10 fmod 2e9
apply nstress -1225000 grad 0 0 9800 range x 0.1 200 y -75 427 z 35.588 125 ；加水荷载，为了避免加到模型的地层周围的面上，用XYZ控制范围。水位高程125m
fix pp 0 ran x 349.509 400 z 39.588 ；下游设置自由透水面
step 30000 ；30000步
这样计算的时候，发现不收敛，最大不平衡力一直处于一个常数，而且为8.9e5，这肯定是错误的。我问过陈育明老师，他说应该是孔压边界的问题，可是我不知道哪里的边界条件出错了。
