--
Fｌac3D 中文手册
ＦLＡC3Ｄ的计算模式中是否需要做孔压分析取决于是否采用conｆig ｆlｕid命令。
１ 无渗流模式(不使用ｃoｎｆig ｆlｕid)
即使不使用命令config ｆｌuｉd，仍然可以在节点上施加孔压。这种模式下,孔压将保持为常量。如果采用塑性本构模型的话，材料的破坏将由有效应力状态来控制。
节点上的孔压分布可由initiａl pp命令或waｔer ｔaｂle命令来设定。
如果采用waｔer table命令，由程序自动计算水位线以下的静水孔压分布。此时,必须施加流体密度(ｗater density)和重力(seｔ gｒaviｔy)。流体密度值和水位位置可以用命令pｒinｔ wateｒ显示。如果水位线是由ｆace关键字来定义的,则可用命令ｐlot watｅr命令显示水位。
这两种情况,单元的孔压都由节点孔压值平均求出,并在本构模型计算中用作有效应力。这种计算模式下，体积力中不反映流体的出现:用户必须根据水位线以上或以下相应地指定干密度和湿密度。使用命令priｎt gp pp和priiｎt zｏne pｐ可分别得到节点或单元孔压。ｐlot coｎtour ｐp命令可绘出节点孔压云图。 ２ 渗流模式(使用ｃoｎｆig fｌuｉｄ)
如果使用命令conｆiｇ flｕid，则可进行瞬时渗流分析，孔压改变和潜水面的改变都可能出现。在conｆig ｆluid模式下，有效应力计算
--
--
(静态孔压分布）和非排水计算均被执行。除此之外,还可进行全耦合分析,这种情况下,孔压改变将使固体产生变形,同时体积应变反过来影响孔压的变化。
如果采用渗流模式，单元孔压仍由节点孔压平均求出。但这种模式，用户只能指定干密度（不论是水位以上还是以下)，因为FＬＡＣ3D将流体的影响考虑到了体积力的计算中。
采用渗流模式时，渗流模型必须施加到单元上,使用命令ｍｏｄｅl fl_isotｒopｉc模拟各向同性渗流，mｏｄｅl ｆl＿aｎisotropic模拟各向异性渗流,ｍodel ｆl_nulｌ模拟非渗透物质。注意，力学模型为空的单元并不代表渗流模型为空。
流体性质(参数)可施加到单元或节点上。各向同性渗透率、孔隙率、比奥系数和非排水热系数等单元流体性质由命令property施加。 对于各向同性渗流,渗透率通过perｍ关键字赋予。对各向异性渗流，渗透率的3个主值采用关键字ｋ1，k2,k3赋予,主方向由关键字fdｉp,fdd，ｆrｏt确定。渗透率的主方向服从右手系统。fdip和ｆｄd分别为k1和ｋ2确定的平面的倾向和倾角。frｏｔ为ｋ1轴和倾角矢量的旋转角。如果不特别指定,比奥系数默认为1，孔隙率默认为0.５。节点的渗流性质由命令initiaｌ指定。这些性质包括流体重度、流体体积模量、比奥模量、流体抗拉强度和饱和度。每种性质在空间上都可以变化。流体重度也可以用ｗatｅr命令给出。
在渗流模式里，有必要知道可压缩性被定义在以下两种参数中:（1）比奥系数和比奥模量；(2）流体体积模量和孔隙率。第一种参数表征
--
--
的是固体颗粒的可压缩性(对不可压缩颗粒，比奥系数设为１)。对第二种参数,固体颗粒被认为是不可压缩的。
单元属性可由命令pｒinｔ zｏne ｐropertｙ显示,节点属性由prinｔ gp命令显示。流体重度,如果随着水位位置被确定，则可由ｐriｎt ｗater命令显示。渗流性质可由命令ｐloｔ bcｏntour pｒopeｒｔy显示。对于各向异性渗流,渗透率的各球形分量可通过使用单元的属性关键字kxｘ,ｋyy，ｋzz,kxy,kｘz,kyz来显示(注意，这些球形分量不可被直接初始化）。
初始节点孔压分布的施加对于渗流模式和非渗流模式都是一样的(如,要么用iｎitial pp命令或用waｔeｒ table命令）。在指定节点可用命令ｆｉｘ ｐp或fｒee pp对孔压固定或释放。流体涌入或渗漏或可由命令ａpply施加。
渗流计算由命令set fluid和ｓoｌve控制。如,set flｕｉｄ on或ofｆ命令开启或关闭渗流计算模式。具体使用开启或关闭模式取决于渗流分析的耦合程度。
渗流分析结果以下面这些命令给出。命令ｐｒiｎt gp pp和print zｏne pp分别给出节点和单元孔压。节点和单元孔压历史可由命令ｈistroy gp pp和命令hiｓtｏry zone pp进行监测。对于瞬时计算,孔压与时间的关系可由命令histｏrｙ flｔime监测。命令plot contouｒ pp绘出节点孔压云图。命令ｐlｏt conｔoｕr satｕｒaｔioｎ绘出饱和度云图。命令pｌot ｆluid绘出流量矢量图。渗流模式的所有信息由命令pｒｉnt fluiｄ命令给出。FISH还提供了一些渗流变
--
Flac3D 中文流体计算
