我的做法是:(1)初始应力场:先在无水的情况下计算初始应力平衡。(2)计算瞬态渗流场:然后关闭力学模式打开渗流模式,施加坡面迎水面孔压边界条件,迭代指定时间得到瞬态渗流场,通过PL CON PP 及PL ADD FLOW RED可以看到库水倒灌坡内地下水后的渗流场。(3)计算渗流的力学响应:关闭渗流场,并设流体模量为0(以免力学计算引起孔压的变化)弹性模型下计算至平衡状态。(4)用FLAC-3D内置强度折减法求安全系数。计算中不考虑土体的湿化效应,即不考虑抗剪强度参数随含水量的变化(因为没试验很难量化,FLAC-3D里容易实现),也不考虑非饱和区基质吸力对稳定性的贡献。不知道我的计算过程对不对?水位上升时安全系数是一直递增的,感觉不怎么合理,应像我上面分析的原因一样先递减再递增。下降过程的安全系数是先递减然后递增,还算合理!望版主能指点!谢谢!