主要模拟的库水位升降条件下库岸边坡的稳定状态，考虑水位的骤升，骤降。不耦合时先进行渗流场模拟，每次考虑水位骤升5M，具体做法是：在初始应力场基础上，关闭力学模式打开渗流模式，施加孔压边界，指定迭代步数得到上升不同水位的渗流场；到最高水位时迭代至稳定状态，然后逐步解除部分边界条件释放孔压，指定相应迭代步数得到下降不同水位的渗流场；得到的孔压跟流速矢量都还比较符合。然后每个工况关闭渗流打开力学模式，在弹性模型下计算至瞬态稳定（即考虑孔压对边坡的影响），然后在摩尔-库仑模型下用内置强度折减法求安全系数，上升过程的模拟到15M水位出现了局部破坏，塑性区不贯通。下降的工况算不出安全系数，迭代很长时间上限值一直未出来，不知道什么原因？FLAC-3D里面用PL CONT SZZ命令输出的是总应力吧？渗流场计算结束时的Z向应力图是正确的，但关闭渗流模式计算孔压对边坡的力学响应时应力图好象不对，因为被水淹没的坡面应力是零。难道是关闭渗流模式后，边界上的空压消失？在这个环节是否要给坡面施加水压力？ yV"ZRrjO'Z
