考虑渗流的有效应力分析

3DEC流体-力耦合分析(Coupled hydro-mechanical analysis)】，两个关键命令行如下：

model config fluid model mechanical off

默认情况下，当打开流体模型后，假定初始饱和度是1，既全饱和。然后随着计算的进行，单元会变得不饱和。手册上说将所有的饱和度设置为0也能得出相同的结果，但目前没有进行测试。

zone fluid cmodel assign isotropic

zonefluidpropperm 1e-10poros 0.3

zone water density 1000

• 设置初始孔隙压力(Pore Pressure)  
  
• 降雨对边坡稳定的影响---FLAC的两相流模拟  
  
• FLAC2D 9.0 孔隙压力(Pore Pressure)  
  
• 3DEC孔隙压力计算(Pore Pressure)  
  

(4) 为了使模型尽快达到稳态流状态，设置水的体积模量为一个较低的数值10000。

zone gridpoint initialize fluid-modulus 10000 zone gridpoint initialize fluid-tension 0

(5) 使用model solve 命令并为流体指定求解条件来确定稳定流动状态。当地下水流动比率低于设定值0.001时，就实现了稳态流动。

model solve fluid ratio-flow 1e-3

(6) 一旦产生了孔隙压力，即可打开力学功能，关闭流体计算，并将水的体积模量设定为零，以防止由于力学变形产生孔隙压力：

model mech active onmodel fluid active offzone gridpoint initialize fluid-modulus 0

图中的黄色虚线是计算出来的浸润线，箭头代表流量变化，不同的充填色代表网格点孔隙压力。