请教几个FLAC3D的问题！我论文做的考虑流固耦合的库区边坡极限分析。在非耦合时问题如下：
（1）渗流场计算：在初始应力平衡的基础上，关闭力学模式打开渗流模式，考虑水位的骤升、骤降，每次水位变化5M，指定一定的迭代步数（或时间）得到瞬态渗流场。思路及设置应该正确吧？我通过计算得到了不同上升水位及不同时间的瞬态渗流场，孔压及流速矢量分布合理。
（2）计算瞬态渗流场对应力场的响应。具体做法是：关闭渗流模式，打开力学模式，设置弹性模型并设置流体模量为0（以免力学作用引起孔压变化），淹没段迎水坡面施加对应水压力，计算至平衡状态。这时通过PL CONT SZZ观察Z向应力图，发现与（1）中的Z向应力差别很大，最大值变小，总应力与受算总应力有差别。不知道是否正确？
（3）在（2）基础上进行安全系数的求解。在摩尔-库仑模型下，指定土体抗剪强度属性用内置强度折减法计算安全系数，很稀疏的网格收敛仍很慢。不知道什么原因？但计算得到上升过程的安全系数是递增的，且大于无水时安全系数。我认为应该是不正确的，因为水位刚上升阶段，坡脚受到浮托力的作用，会使得阻滑力下降很多，下滑力影响不大。且没有考虑土体的湿化（即C跟PHI不变），渗透力很小，这时计算得到的安全系数应该变小。上升到一定水位后，浮力作用大使得下滑力减小很大，安全系数再递增。而下降过程的安全系数应该是递减的。计算结果我觉得不对，不知道什么原因？
望高手指点，感激不尽！！谢谢！
