Flac3D中文流体计算

Flac3D 中文手册

FLAC3D 的计算模式中是否需要做孔压分析取决于是否采用config fluid 命令。

1 无渗流模式（不使用config fluid ）即使不使用命令config

fluid ，仍然可以在节点上施加孔压。这种模式下，孔压将保持为常量。如果采用塑性本构模型的话，材料的破坏将由有效应力状态来控制。

节点上的孔压分布可由initial pp 命令或water table 命令来设定。

如果采用water table 命令，由程序自动计算水位线以下的静水孔压分布。此时，必须施加流体密度（water density ）和重力（set

gravity ）。

流体密度值和水位位置可以用命令print water 显示。如果水位线是由face 关键字来定义的，则可用命令plot water 命令显示水位。

这两种情况，单元的孔压都由节点孔压值平均求出，并在本构模型计算中用作有效应力。这种计算模式下，体积力中不反映流体的出现：用户必须根据水位线以上或以下相应地指定干密度和湿密度。使用命令print gp pp 和priint zone pp 可分别得到节点或单元孔压。plot contour pp 命令可绘出节点孔压云图。

2 渗流模式（使用config fluid ）如果使用命令config fluid ，则可进行瞬时渗流分析，孔压改变和潜水面的改变都可能出现。在config fluid 模式下，有效应力计算（静态孔压分布）和非排水计算均被执行。

除此之外，还可进行全耦合分析，这种情况下，孔压改变将使固体产生变形，同时体积应变反过来影响孔压的变化。

如果采用渗流模式，单元孔压仍由节点孔压平均求出。但这种模式，用户只能指定干密度（不论是水位以上还是以下），因为FLAC3D

将流体的影响考虑到了体积力的计算中。

采用渗流模式时，渗流模型必须施加到单元上，使用命令model

fl_isotropic 模拟各向同性渗流，model fl_anisotropic 模拟各向异性渗流，model fl_null模拟非渗透物质。注意，力学模型为空的单元并

不代表渗流模型为空。

流体性质（参数）可施加到单元或节点上。各向同性渗透率、孔隙率、比奥系数和非排水热系数等单元流体性质由命令property 施加。

对于各向同性渗流，渗透率通过perm 关键字赋予。对各向异性渗流，渗透率的3 个主值采用关键字k1,k2,k3 赋予，主方向由关键字

fdip,fdd,frot 确定。渗透率的主方向服从右手系统。fdip 和fdd 分别为

k1和k2确定的平面的倾向和倾角。frot为k1轴和倾角矢量的旋转角。

如果不特别指定，比奥系数默认为1，孔隙率默认为0.5。节点的渗

流性质由命令initial 指定。这些性质包括流体重度、流体体积模量、比奥模量、流体抗拉强度和饱和度。每种性质在空间上都可以变化。

流体重度也可以用water 命令给出。

在渗流模式里，有必要知道可压缩性被定义在以下两种参数中：（1 ）比奥系数和比奥模量；（2）流体体积模量和孔隙率。第一种参数表

征的是固体颗粒的可压缩性（对不可压缩颗粒，比奥系数设为1）。

对第二种参数，固体颗粒被认为是不可压缩的。

单元属性可由命令print zone property 显示，节点属性由print gp 命令显示。流体重度，如果随着水位位置被确定，则可由print water

命令显示。渗流性质可由命令plot bcontour property 显示。对于各向异性渗流，渗透率的各球形分量可通过使用单元的属性关键字

kxx,kyy,kzz,kxy,kxz,kyz 来显示（注意，这些球形分量不可被直接初始化）。

初始节点孔压分布的施加对于渗流模式和非渗流模式都是一样的

如，要么用initial pp 命令或用water table 命令）。在指定节点可用命令fix pp 或free pp 对孔压固定或释放。流体涌入或渗漏或可由命令apply 施加。

渗流计算由命令set fluid 和solve 控制。如，set fluid on 或off 命令开启或关闭渗流计算模式。具体使用开启或关闭模式取决于渗流分析

的耦合程度。

渗流分析结果以下面这些命令给出。命令print gp pp 和print zone

pp

分别给出节点和单元孔压。节点和单元孔压历史可由命令histroy gp pp 和命令history zone pp 进行监测。对于瞬时计算，孔压与时间的关系可由命令history fltime 监测。命令plot contour pp 绘出节点孔压云图。命令plot contour saturation 绘出饱和度云图。命令plot fluid

绘出流量矢量图。渗流模式的所有信息由命令print fluid 命令给出。

FISH 还提供了一些渗流变量。其中一个与节点有关的变量gp_flow ，只能通过FISH 函数使用。该变量描述了通过节点的净流入或流出量。

因为可以提供一个系统总的流入或流出量，这些流量的统计在孔压固定的边界是很有用的。

渗流边界条件，初始条件

FLAC3D 默认为不透水边界，即认为所有节点上的孔压随着从邻近单元流入或流出的量发生自由变化。可以使用命令fix pp 将节点上的孔压设为

“自由”，也可使用free pp 使节点上的孔压“固定”。如果孔压固定，流体可以在外边界上流入或流出节点。下面总结这两种边界条件的影响：1，孔压自由

这是默认的不透水边界条件。节点与外界之间不发生流量交换。系统

根据当前饱和度值和流体是否形成涡凹现象来计算压力和饱和度变化。

2，孔压固定

这是一种流体通过外界流入或流出的边界条件。如果设定孔压为0，

饱和度才可能变化。否则，饱和度被设为1(FLAC3D 假设孔压只在完全饱和材料中存在)。孔压不能被固定在低于拉力极限的值，如果出现这种情况，