1 引言
3DEC能够模拟非常复杂的流体流动以及流体与力的相互作用。总的来说,可以计算如下四大类问题:
(1) 孔隙压力(Pore Pressure)
(2) 基质流体流动(Matrix Fluid Flow)
(3) 节理流体流动(Joint Fluid Flow)
(4) 支撑剂(Proppant)
大多数实践的岩土工程稳定性分析不必考虑流体流动,只考虑(1)静态孔隙压力(Pore Pressures)即可,尽管(2)(3)(4)可以计算流体-力的单向或双向相互作用,但它们主要的应用领域是水文分析和石油工程,因此我们暂不考虑。这意味着当进行静态孔隙压力计算时,不需要开启如下的流体流动命令:
model configure fluid onmodel configure matrixflow
在这种模式下,不必设置孔隙率和渗透率等流体属性,也不必设置流体边界条件。单元内的孔隙压力用来计算有效应力以确定单元是否破坏,孔隙压力以相同的命令分配给节理接触,节理的有效法向应力用来确定接触是否破坏。网格点上的孔隙压力不发生变化,但在使用塑性本构模型时,可能会诱发由有效应力状态控制的破坏。本文简要回顾了3DEC孔隙压力计算使用的命令。
block gridpoint initialize pore-pressure
block gridpoint list pore-pressure
block insitu pore-pressure 1e6 gradient 0 0 9810
block zone list pore-pressure
model newmodel random 10000 block create brick 0,100 0,100 0,100block cut joint-set dip 60 dip-direction 90 spacing 15block zone generate edgelength 10model gravity 0 0 -9.81block water density 1000block water plane normal 0.17,0,0.98 origin 50,50,90
model solve ratio 1e-7
block gridpoint initialize displacement 0 0 0
block fluid property density 1000.0
block water plane origin 0 0 10 normal 0 0 1
model solve