按照命令流的理解,该代码的编写者对模型的处理基本流程如下:
(1)先用弹性模型 model elas 和流体模型 model fl_iso 获得初始平衡应力场,即只打开力学进程,关闭流体进程。在设定流体边界条件时,通过 ini pp 命令给定了初始水压,并用 fix pp 命令固定了模型x方向和z方向的孔压。
(2)然后使用摩尔库伦模型 model mo 和流体模型 model fl_iso ,将开挖部分设定为不透水模型 model fl_null 并设定孔压 pp 为0;打开力学进程和流体进程,采用主从进程法进行耦合计算。
现在需要反思的是:(1)在初始平衡场的时候,仅用弹性模型求解而不用摩尔库伦模型直接求解,是否合适。这是否与模型设定的初始应力场为静水应力场有关。(2)开挖分析中,对于开挖部分设定了不透水模型 model fl_null ,为何还需要对其孔压通过 ini pp 0 进行设定,这样的设定与 fix pp 的差别又在哪里。(3)实际效果没有相关图片验证,比较难以理解边界条件的合理性。