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 将流体的影响考虑到了体积力的计算中。

采用渗流模式时，渗流模型必须施加到单元上，使用命令modelfl_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假设孔压只在完全饱和材料中存在）。

孔压不能被固定在低于拉力极限的值，如果出现这种情况，FLAC3D会将其设定到拉力极限值。

如前所述，边界条件不是任意的。

FLAC3D在进行计算前会“检查”并“修正”这些条件。

可使用fix pp命令将孔压固定在某个值，也可在外边界或内边界上使用命令apply pp。

如果边界条件被用于一个非表面节点，则必须加关键字interior。

apply命令具有可以用“历史”命令进行监测的优点。

渗流边界条件可以通过apply命令用在单个或部分节点、单元面或单元上。

命令apply pwell为边界节点指定了一个流入或流出井。

如果加上interior关键字，则该条件用于内部节点。

命令apply discharge 和apply leakage为边界单元的表面分别指定了涌出和渗漏边界条件。

命令apply vwell为指定区域内的单元提供一个流速。

这些边界条件除了apply leakage外，均可使用history监测命令。

具有固定孔压节点就好像是流入源或流出源。

没有直接的命令显示这些节点的流入或流出量。

但可通过FISH变量gp_flow来记录。

孔压的初始分布，孔隙率，饱和度和流体属性可通过命令initial或property施加。

如果还加了重力，则孔压初始分布应与重力梯度，水的重度和节点饱和度和孔隙率相容。

如果这些初始分布不相容，则计算开始时所有单元中将出现流体流动。

因此，应在模拟开始时设一定的计算步来检验初始条件是否相容。

如果模型中含有接触面，有效应力将沿着这些接触面进行初始化（即：在节点应力初始化时，认为接触面应力包含孔压）。

water lable命令将包含沿着接触面的孔压，这是因为定义在单元节点上的孔压也在接触面节点上。

如果接触面的上下两面连在一起，在没有阻力时，将发生穿越接触面的流体流动。

但程序不对沿着接触面的流体流动（裂隙流）进行计算。

单渗流与渗流耦合问题FLAC3D既能进行单渗流分析，也能进行固流耦合分析。

耦合分析可由FLAC3D内置力学模型完成。

但要注意，渗流模型中的空单元并不是力学空单元。

必须用命令model fl_null给单元赋予流体空属性。

对于耦合过程，FLAC3D提供了几种计算模式。

其中之一是假设孔压一旦被赋予便不再改变。

该方法并不要求任何额外空间存储计算过程。

除此之外涉及到渗流的计算模式都要求使用命令config fluid。

命令model fl_iso使所有单元中都能发生渗流。

不同的耦合计算模式在下面讨论。

一般情况下，在能跟所模拟问题的物理过程相似的情况下，应使用尽可能简单的模式。

计算模式的选择根据以下几个方面确定。

时间比例对所需模拟的渗流或耦合问题用FLAC3D估计与涉及的不同进程相关的时间比例是非常有用的。

对有关研究问题的时间度量和扩散性的认识有助于估计最大网格宽度、最小区域尺寸、时步大小和计算可行性。

如果不同进程的时间比例相差太大，则很可能采用一种简单的（非耦合）方法。

时间比例可用特征时间给出。

以下这些由量纲分析得出的定义，都是基于解析的连续源理论表达式。

它们可用于得出FLAC3D分析的大致时间比例。

力学过程特征时间、流体扩散过程特征时间流体扩散率FLAC3D中使用了取决于控制过程的储水系数的几种形式：流体存储系数、地下潜水相存储系数、弹性存储系数以上定义，有几点特性值得注意：（1）因为FLAC3D中显式的时步对应于最小区域中信息从一个节点传到下一节点所需要的时间，时步的大小可用计算特征时间公式中特征长度的最小区域来估计。

重要的是注意FLAC3D中在用流体扩散率（即使是在耦合模拟中）计算显式流体时步。

因此，时步的大小可用特征长度的最小区域尺寸来估计。

（2）在饱和流体问题中，简化的体积模量不但导致时步的增加，同样导致到达稳定状态时间的增加，所以总步数增加，该总步数可用模型和最小区域的特征长度来估计。

（3）在部分饱和流体流动问题中，可通过调整流体体积模量加速收敛以趋于稳定状态，但要注意不可将体积模量减小太多以至产生数值不稳定。

数值稳定条件能由流体储量在一个特征长度区域的高度上必须保持低于地下潜水储量的要求推导而出。

（4）为避免扩散问题中的边界效应，模型的特征长度必须大于某个尺度。

同样，最小模拟时间由某个关系式控制。

（5）在耦合流体问题中，实际扩散率由流体刚度与岩土介质的刚度比来控制。

完全耦合模拟方法的选择用FLAC3D进行完全耦合的准静态固流耦合分析通常要耗费大量时间，且有时候并不必要。

很多情况下，可使用不同程度的非耦合方法简化分析并加快计算速度。

下面的例子给出了对应于流固耦合的不同水平的FLAC3D模拟方法。

选择计算方法时有3个主要的因素需要考虑：（1）模拟时间比例和扩散过程的特征时间；（2）耦合过程中强制扰动特性；（3）流固刚度比。

时间比例首先通过从扰动的开始阶段计算时间来考虑时间比例因素。

定义分析所需要的时间（模拟时间），对应于耦合扩散过程的特征时间。

短期行为（不排水）如果对应于耦合扩散特征时间，分析所需时间非常短，在模拟结果中流体流动的影响几乎可以忽略不计，则可采用不排水模拟（config fluied,set fluid off）。